Nghiên cứu đã sử dụng ảnh Sentinel 2 tại tháng 11/2016 và tháng 1/2021 trên địa bàn huyện Yên Châu, cùng với việc xây dựng tám nhân tố để kiểm tra mối quan hệ là độ cao, độ dốc, hướng phơi, khoảng cách tới đường giao thông, khoảng cách tới khu dân cư, nhiệt độ trung bình tháng, lượng mưa trung bình năm và loại đất. Mời các bạn cùng tham khảo!
Lâm học BIẾN ĐỘNG LƯỢNG THỰC VẬT CHE PHỦ VÀ MỐI QUAN HỆ VỚI CÁC NHÂN TỐ TỰ NHIÊN XÃ HỘI TẠI YÊN CHÂU, SƠN LA Bùi Mạnh Hưng1, Nguyễn Thị Bích Phượng1, Nguyễn Thị Thảo1 Trường Đại học Lâm nghiệp TÓM TẮT Biến động lượng thực vật che phủ mặt đất vấn đề quan tâm giải tỉnh Nghiên cứu sử dụng ảnh Sentinel tháng 11/2016 tháng 1/2021 địa bàn huyện Yên Châu, với việc xây dựng tám nhân tố để kiểm tra mối quan hệ độ cao, độ dốc, hướng phơi, khoảng cách tới đường giao thông, khoảng cách tới khu dân cư, nhiệt độ trung bình tháng, lượng mưa trung bình năm loại đất Kết tính toán số thực vật (NDVI) cho thấy giá trị trung bình số thực vật năm 2016 cao năm 2021 0,1 Nhìn chung, lượng thực vật có địa bàn huyện bị suy giảm Tổng diện tích có số thực vật nhỏ giảm từ năm 2016 đến 2021 khoảng 55,5 Đồng thời, tổng diện tích khu vực có số thực vật mức độ thấp (NDVI từ - 0,6) có xu hướng tăng mạnh, từ khoảng 42 nghìn lên đến 69 nghìn Ngược lại, diện tích khu vực có số thực vật từ 0,6 - 0,9 lại có xu hướng giảm Kết phân tích mối tương quan NDVI với nhân tố cho thấy mức độ biến động số thực vật từ 2016 đến 2021 có quan hệ chặt với yếu tố độ dốc hướng phơi, sau đến độ cao, nhiệt độ trung bình khoảng cách đến đường Các nhân tố lượng mưa bình quân tháng loại đất gần khơng có mối quan hệ Từ khóa: ảnh Sentinel 2, biến động, lượng thực vật, nhân tố tự nhiên, Yên Châu ĐẶT VẤN ĐỀ Ở Việt Nam, rừng biến động diện tích rừng vấn đề đáng quan tâm Bởi lẽ, thời gian gần đây, nhiều khu vực địa bàn nhiều tỉnh khác xảy vụ rừng, phá rừng nghiêm trọng để lại hậu thiệt hại lớn tài sản cho đơn vị quản lý cho nhà nước Mất rừng cịn ngun nhân ảnh hưởng tới mơi trường khơng khí, sói mịn đất chất lượng nước nhiều địa phương Đặc biệt, rừng gây lên lũ lụt, sạt lở ngày trầm trọng nhiều địa phương Đồng thời, rừng, suy thoái rừng góp phần vào gây biến đổi khí hậu hiệu ứng nóng lên tồn cầu Ảnh hưởng lớn đến đời sống người, đặc biệt cộng đồng sống gần rừng Trong thời kỳ 1945 – 1975 nước khoảng triệu rừng, bình quân 100.000 năm Quá trình rừng diễn nhanh giai đoạn 1975 – 1990: 2,8 triệu ha, bình quân 140.000 ha/ năm Tuy nhiên từ năm 1990 đến nay, công tác trồng rừng chương trình dự án nhà nước đẩy mạnh phần làm cho diện tích rừng tăng lên cách đơn Tuy nhiên, rừng diễn nhiều địa phương (Đỗ Trọng Hoàn Nguyễn Hải Vân, 2017; Hương Thảo, 2010) Biến động tài nguyên rừng có suy thối tang rừng chịu ảnh hưởng nhiều yếu tố tự nhiên xã hội khác Ảnh hưởng nhân tố kinh tế xã hội tới rừng, suy thoái rừng phân tích nhiều nghiên cứu nước Các yếu tố tự nhiên hạn hán, cháy rừng, nhiễm khơng khí, xói mịn đất có mối quan hệ mật thiết tới suy thối rừng Các yếu tố kinh tế xã hội canh tác nương rẫy, lấn chiếm rừng, thu gom củi, khai thác gỗ để sản xuất than góp phần lớn vào rừng suy thoái rừng, đặc biệt nước Đông Nam Á (Hồng Nhung, 2018; Hương Thảo, 2010) Một số nhân tố thu nhập bình quân, thu nhập từ rừng, nhu cầu sử dụng gỗ phát triển cơng nghiệp có ảnh hưởng lớn đến tỷ lệ rừng tỉnh Điện Biên (Lã Nguyên Khang Trần Quang Bảo, 2015) Tuy nhiên, tới thời điểm phân tích ảnh hưởng nhân tố tới suy thoái rừng, rừng tang rừng chủ yếu dựa vào kết điều tra mặt đất TẠP CHÍ KHOA HỌC VÀ CƠNG NGHỆ LÂM NGHIỆP SỐ - 2021 31 Lâm học vấn với dung lượng hữu hạn, điều phần làm giảm độ xác kết phân tích Sử dụng cơng nghệ viễn thám để phân tích ảnh vệ tính, ảnh mơ hình số độ cao từ phân tích mối quan hệ với nhân tố tự nhiên, xã hội nâng cao độ xác hơn, số lượng điểm mẫu lớn Huyện Yên Châu, tỉnh Sơn La huyện có tài nguyên rừng tương đối phong phú Trong năm gần đây, cơng tác quản lý bảo vệ rừng, phịng cháy chữa cháy rừng địa bàn huyện quan tâm triển khai từ huyện đến sở Song ảnh hưởng thời tiết, ý thức quản lý, ý thức bảo vệ rừng người dân, nhu cầu đời sống tăng cao, sản xuất nông nghiệp đẩy mạnh với chuyển đổi diện tích đât sử dụng nên tượng rừng, suy thoái rừng diễn phức tạp Bên cạnh đó, biện pháp bảo vệ rừng, trồng rừng, nhiều khu vực địa bàn huyện có cải thiện đáng kể cho tài nguyên rừng, đặc biệt mặt trữ lượng thực vật mặt đất Tuy nhiên, chưa có nghiên cứu cụ thể địa bàn toàn huyện biến động rừng mối quan hệ chúng với nhân tố tự nhiên chưa thực (Ngọc Thuấn, 2020; UBND Yên Châu, 2018) Từ lý trên, nghiên cứu biến động rừng phân tích mối quan hệ biến động với yếu tố tự nhiên huyện Yên Châu cần thiết cho công tác quản lý bảo vệ tài nguyên rừng địa bàn Nghiên cứu tập trung chủ yếu khai thác sử dụng phần mềm Qgis (một phần mềm 4M: miễn phí, mạnh, mở mới) (Bùi Mạnh Hưng Nguyễn Thanh Thủy Vân, 2019) công nghệ viễn thám để thực tốn phân tích khơng gian đối tượng Raster huyện Yên Châu nhằm mục đích: (1) xây dựng đồ trạng số thực vật năm 2016 2021; (2) xây dựng đồ biến động lượng thực vật che phủ cho giai đoạn từ 32 11/2016-1/2021 (3) phân tích mối quan hệ biến động rừng với nhân tố tự nhiên, xã hội địa bàn làm sở cho việc quản lý tài nguyên rừng bền vững khu vực nghiên cứu PHƯƠNG PHÁP NGHIÊN CỨU 2.1 Khu vực nghiên cứu Yên Châu huyện miền núi, biên giới tỉnh Sơn La, nằm trục Quốc lộ 6, cách Hà Nội 240 km theo hướng Tây Bắc, cách thị xã Sơn La 64 km phía Đơng, khu vực đệm Cao nguyên Nà Sản Mộc Châu, có toạ độ địa lý: 104o10’ – 104o40’ kinh độ Đông, 21o07’-21o14’ vĩ độ Bắc, phía Đơng giáp huyện Mộc Châu, phía Tây giáp huyện Mai Sơn, phía Bắc giáp huyện Bắc Yên, phía Nam tiếp giáp với nước CHDCND Lào, có 47 km đường biên giới với nước Lào (UBND Yên Châu, 2018) 2.2 Phương pháp thu thập số liệu Các số liệu điều kiện tự nhiên, dân sinh kinh tế tình hình bảo vệ, phát triển tài nguyên rừng địa bàn huyện kế thừa từ hạt Kiểm lâm huyện Các báo cáo tổng kết công tác hàng năm huyện kế thừa từ UBND nên giám thống kê tỉnh Sơn La Bản đồ kiểm kê rừng kết hợp với đồ cập nhật diễn biến tài nguyên rừng kế thừa từ Chi cục Kiểm lâm Sơn La Các lớp đồ chủ yếu sử dụng để phân tích đường giao thơng, khu dân cư, hệ thống sơng suối, vị trí nương rẫy địa bàn huyện Hệ thống ảnh vệ tinh mơ hình số độ cao sử dụng nghiên cứu bao gồm: Ảnh Sentinel tải từ trang web quan khảo sát địa chất Hoa Kỳ (USGS) mơ hình số độ cao (DEM) tải từ trang web viện Vật lý Địa cầu Đại học Alaska Fairbanks (ASF), lớp nhiệt độ, lượng mưa loại đất tải từ trang web Tổ chức Nông Lương Thế giới (FAO) cụ thể bảng TẠP CHÍ KHOA HỌC VÀ CƠNG NGHỆ LÂM NGHIỆP SỐ 3- 2021 Lâm học Bảng Dữ liệu ảnh Sentinel mơ hình số độ cao sử dụng nghiên cứu Độ phân giải TT Mã ảnh Loại liệu Nguồn (m) S2A_OPER_MSI_L1C_TL_SGS 20161105 Ảnh Sentinel 10 USGS L1C_T48QVJ_A028901_20210103 Ảnh Sentinel 10 USGS AP_15354_FBS_F0400_RT1 DEM 12,5 ASF AP_15354_FBS_F0410_RT1 DEM 12,5 ASF AP_15602_FBS_F0400_RT1 DEM 12,5 ASF 2.3 Phương pháp xử lý, tính tốn 2.3.1 Xử lý ảnh Sentinel tính số thực vật (NDVI) Ảnh Sentinel sau tải xử lý tính tốn số thực vật ảnh sau (Nguyễn Hải Hòa cộng sự, 2016; Nguyễn Trọng Cương cộng sự, 2016) - Bước 1: Xử lý ảnh Sentinel-2A bao gồm bước nhỏ sau: + Gộp kênh ảnh: Khi thu thập ảnh viễn thám từ vệ tinh ảnh nằm kênh phổ khác có màu đen trắng Vì vậy, để tính tốn số thực vật, nghiên cứu tiến hành tổ hợp kênh ảnh Cụ thể kênh ảnh 2, 3, gộp + Tăng cường chất lượng ảnh: Ảnh viễn thám sau tổ hợp tăng cường cho việc giải đoán ảnh tốt + Hiệu chỉnh hình học: Trước phân tích, giải đốn cần kiểm tra thơng tin hệ quy chiếu tham số địa lý ảnh Ảnh vệ tinh nắm chỉnh giảm thiểu sai số hình học qua cho độ xác cao + Cắt ảnh theo ranh giới khu vực nghiên cứu: Ảnh cắt theo khu vực nghiên cứu, cụ thể theo ranh giới huyện Yên Châu Các bước thực phần mềm Qgis 3.16 - Bước 2: Chỉ số thực vật NDVI Trong nghiên cứu, số thực vật NDVI sử dụng để xác định lượng thực vật địa bàn huyện hai thời điểm khác 2016 2021 Chỉ số thực vật hay số thực vật chuẩn hóa (NDVI - Normalized Difference Vegetation Index) đại lượng phản ánh xác số lượng thực vật điểm mặt đất Chỉ số thực vật tính tốn cho ảnh Sentinel theo cơng thức sau: (Kênh − Kênh 4) NDVI = (Kênh + Kênh 4) Trong đó: NDVI số thực vật; Kênh kênh cận hồng ngoại; Kênh kênh màu đỏ (Shivangi S Somvanshi Maya Kumari, 2020) Bảng phân loại giá trị NDVI theo loại trạng thái sử dụng bảng (Lê Thị Thu Hiền, 2013) Bảng Phân loại NDVI theo lớp phủ bề mặt Giá trị NDVI < 0,1 0,1-0,2 0,2-0,3 0,3-0,6 > 0,6 Lớp phủ bề mặt Khu vực cằn cỗi đá; cát; mặt nước; bê tông Đất đá cằn cỗi, bụi Cây bụi trảng cỏ; đất nông nghiệp để trống Trảng cỏ, trồng nông nghiệp, rừng thưa Rừng nhiệt đới 2.3.2 Xử lý mơ hình số độ cao (DEM) lớp đồ Mơ hình số độ cao khu vực nghiên cứu tải xử lý sau: Gộp DEM lại với (do huyện Yên Châu nằm rải DEM khác nhau) Sau đó, điểm trũng khơng có liệu kiểm tra lấp đầy giá trị Tiếp đó, DEM cắt theo ranh giới huyện Yên Châu Mơ hình số độ cao sau sử dụng để tính tốn độ cao cho điểm, tính độ dốc, hướng phơi cho điểm địa bàn huyện DEM sử dụng nghiên cứu có độ phân giải 12,5 m Đối với lớp đồ đường giao thông, khu dân cư… tạo buffer cho TẠP CHÍ KHOA HỌC VÀ CƠNG NGHỆ LÂM NGHIỆP SỐ - 2021 33 Lâm học đối tượng với khoảng cách khác Sau đó, lớp chuyển sang dạng Raster để phân tích hồi quy tuyến tính (Alex Mandel cộng sự, 2016) 2.3.3 Phân tích hồi quy tuyến tính biến động lượng thực vật nhân tố Hiệu số giá trị số thực vật NDVI hai năm 2016 năm 2021 tính sau: ∆NDVI = NDVI2021 - NDVI2016 Sau tính độ lệch ∆, ∆ < khu vực suy thoái lượng thực vật che phủ, tức hàm lượng thực vật bị giảm hai mốc thời gian Cịn ∆ ≥ khu vực tăng lượng thực vật, nghĩa lượng thực vật xanh tăng lên mốc thời gian Sau đó, lớp ∆ sử dụng để phân tích hồi quy tuyến tính biến động lượng thực vật với nhân tố tự nhiên xã hội Phương trình hồi quy tuyến tính đa biến sau: ∆ = a0 + a1.Xi Trong đó: ∆ độ lệch giá trị số thực vật hai mốc thời gian; Xi nhân tố tự nhiên xã hội độ cao, độc dốc, hướng phơi, khoảng cách tới đường, khu dân cư, nhiệt độ, lượng mưa, loại đất KẾT QUẢ VÀ THẢO LUẬN 3.1 Kết phân cấp lớp đồ 3.1.1 Kết số thực vật năm 2016 2021 Chỉ số thực vật tính tốn cho ảnh Sentinel năm 2016 năm 2021 Kết tính tốn phân cấp thể hình a) b) c) d) Hình Kết tính tốn phân cấp số NDVI: a) c) kết tính tốn phân cấp giá trị NDVI cho năm 2016; b) d) kết tính tốn phân cấp giá trị NDVI cho năm 2021 34 TẠP CHÍ KHOA HỌC VÀ CÔNG NGHỆ LÂM NGHIỆP SỐ 3- 2021 Lâm học Các giá trị đặc trưng thống kê tính tốn cho số thực thực vật hai năm thể bảng Bảng Thống kê giá trị NDVI năm Các đại lượng Năm 2016 Năm 2021 Giá trị nhỏ -0,4313 -0,2232 Giá trị lớn 0,8684 0,8125 Giá trị trung bình 0,5602 0,4620 Sai tiêu chuẩn 0,1820 0,1359 Như vậy, giá trị lớn năm 2016 cao năm 2021 0,05 Như vậy, lượng thực vật cao ghi nhận địa bàn có dấu hiệu bị suy giảm Giá trị trung bình số thực vật năm 2016 cao năm 2021 0,1 Điều cho chứng minh rằng, nhìn chung, lượng thực vật có địa bàn huyện bị suy giảm mức độ định diện rộng hơn, giá trị trung bình bị giảm Biến động dãy giá trị số thực vật năm 2016 lớn năm 2021 Điều cho thấy lượng thực vật rải phạm vi toàn huyện năm 2021 Bảng Diện tích tương ứng cấp NDVI năm Năm 2016 Năm 2021 Cấp Số lượng điểm ảnh Diện tích năm 2016 (ha) Số lượng điểm ảnh Diện tích năm 2021 (ha) Cấp (-0.8-0) 12248 122,48 6698 66,98 Cấp (0-0.3) 931954 9319,54 1127943 11279,43 Cấp (0.3-0.6) 3271971 32719,71 5812563 58125,63 Cấp (0.6-0.9) 4370372 43703,72 1639341 16393,41 Kết tính tốn diện tích theo cấp số thực vật địa bàn tồn huyện Tổng diện tích đất trống, sông hồ, nhà cửa đối tượng khơng có thực vật xanh phủ giảm từ năm 2016 đến 2021 khoảng 55,5 Đồng thời, tổng diện tích có thực vật phủ xanh mức độ thấp (từ - 0,6) có xu hướng tăng mạnh, từ khoảng 42 nghìn lên đến 69 nghìn Đây kết việc phủ xanh đất trống đồi núi trọc địa phương Là kết công tác thực dự án trồng rừng địa phương Đồng thời kết công tác bảo vệ phát triển rừng lực lượng kiểm lâm chuyên trách năm qua Tuy nhiên, điều đáng báo động khu vực có lượng thực vật cao (chỉ số thực vật từ 0,6-0,9), rừng thường có trữ lượng lớn (Nguyễn Quang Giáp, 2015; Mai Trọng Thịnh Nguyễn Hải Hòa, 2017) lại có xu hướng giảm Điều thể diện tích rừng có trữ lượng lớn bị giảm dần, rừng có dấu hiệu suy thối chất lượng Kết tương tự kết nghiên cứu nhiều tác giả nước (Nguyễn Hải Hòa cộng sự, 2016; Nguyễn Hải Hòa Nguyễn Văn Quốc, 2017) 3.1.2 Kết đồ nhân tố tự nhiên xã hội Từ việc tính tốn mơ hình số độ cao, lớp thơng tin điều kiện tự nhiên, xã hội địa bàn huyện, kết đồ nhân tố thể hình đây, bao gồm nhân tố là: độ cao, độ dốc, hướng phơi, khoảng cách tới đường giao thông, khoảng cách tới khu dân cư, nhiệt độ trung bình tháng, lượng mưa trung bình năm loại đất TẠP CHÍ KHOA HỌC VÀ CƠNG NGHỆ LÂM NGHIỆP SỐ - 2021 35 Lâm học a) b) c) d) e) f) g) h) Hình Bản đồ nhân tố: a) độ cao; b) độ dốc; c) hướng phơi (tính theo độ); d) khoảng cách tới đường giao thông; e) khoảng cách tới khu dân cư; f) nhiệt độ trung bình tháng; g) lượng mưa trung bình năm; h) loại đất 36 TẠP CHÍ KHOA HỌC VÀ CÔNG NGHỆ LÂM NGHIỆP SỐ 3- 2021 Lâm học Kết cho thấy độ cao toàn huyện dao động từ 150 m đến 1500 m Độ dốc dao động từ đến 75 độ Do huyện miền núi, phần lớn diện tích huyện có độ dốc tương đối lớn (khoảng từ 35 độ đến 75 độ) Địa bàn huyện có đầy đủ hướng phơi đơng, tây, nam, bắc Khoảng cách đến đường giao thông tối đa đến 2200 m đến khu vực đông dân cư 4000 m Nhiệt độ trung bình tháng dao động từ 19 đến 23 độ, lượng mưa trung bình năm chạy từ 91 mm đến 113 mm Trên địa bàn huyện bao gồm loại đất Ferralsal, Cambisol, Leptosol, Luvisol Vertisol Trong 80% diện tích đết Ferralsal 3.2 Kết biến động lượng thực vật che phủ Bản đồ kết phân cấp suy thoái gia tăng lượng thực vật che phủ địa bàn toàn huyện thể hình sau Với ngưỡng chia NDVI sau: suy thối mạnh: < -0,6, trung bình: -0,4 đến -0,2, yếu: -0,2 đến Tăng lượng thực vật mạnh: > 0,6, trung bình: 0,2 - 0,4, yếu - 0,2 Hình Bản đồ mức độ suy thoái gia tăng lượng thực vật che phủ Bản đồ cho thấy phần lớn diện tích địa bàn huyện diễn tượng suy thoái lượng thực vât che phủ (diện tích có tơng màu vàng cam đỏ) Có nghĩa lượng thực vật xanh mặt đất có xu hướng giảm diện rộng Cụ thể thể bảng Bảng Biến động diện tích suy thoái tăng lượng thực vật che phủ Biến động Suy thoái lượng thực vật Tăng lượng thực vật Mức độ Mạnh Trung bình Yếu Tổng Yếu Trung bình Mạnh Tổng Diện tích (m2) 9238906,25 169108750,00 499156875,00 677504531,25 152008593,75 27500625,00 1775312,50 181284531,25 Diện tích (ha) 923,89 16910,88 49915,69 67750,45 15200,86 2750,06 177,53 18128,45 TẠP CHÍ KHOA HỌC VÀ CƠNG NGHỆ LÂM NGHIỆP SỐ - 2021 37 Lâm học Bảng cho thấy rang vòng năm từ năm 2016 đến năm 2021 diện tích khu vực có lượng thực vật che phủ bị giảm lớn gấp 3,7 lần so với diện tích khu vực có lượng thực vật tăng Đây vấn đề lớn đặt cho nhà quản lý tài nguyên thực vật nói chung tài nguyên rừng nói riêng địa bàn huyện Trong khu vực suy thoái lượng thực vật suy thối mạnh chiếm 1,36%, mức độ trung bình chiếm 24,96% mức độ yếu chiếm 73,68% Tại khu vực lượng thực vật tăng mức độ mạnh chiếm 0,98%, cịn mức độ trung bình yếu chiếm tỷ lệ 15,17% 83,32% Như vậy, khu vực suy thoái tăng lượng thực vật tỷ lệ mức yếu chiếm tỷ lệ cao nhất, sau đến mức trung bình cuối mức mạnh 3.3 Mối quan hệ biến động lượng thực vật che phủ nhân tố tự nhiên xã hội Kết phân tích tương quan tuyến tính biến động giá trị số thực vật với yếu tố tự nhiên xã hội thể bảng Bảng Kết phân tích tương quan tuyến tính NDVI nhân tố Tham số Độ cao Độ dốc Hướng phơi Khoảng cách đến đường Khoảng cách đến dân cư Nhiệt độ trung bình tháng Lượng mưa trung bình tháng Loại đất a b R2 -0,0502 -0,0546 -0,0546 -0,0786 -0,0754 -0,4492 -0,1540 -0,0971 -0,0001 -0,0021 -0,0021 -0,0001 0,0000 0,0163 0,0006 0,0001 0,0146744 0,0254256 0,0254256 0,0100982 0,0092571 0,0139443 0,0002778 0,0000015 Kết cho thấy mức độ biến động số thực vật khoảng năm từ 2016 đến 2021 có quan hệ chặt với yếu tố độ dốc hướng phơi, sau đến độ cao, nhiệt độ trung bình khoảng cách đến đường Các nhân tố lượng mưa bình quân tháng loại đất gần khơng có mối quan hệ Điều có nghĩa mức độ suy thối tăng lượng thực vật che phủ chịu ảnh hưởng lớn nhân tố độ dốc, độ cao, khoảng cách đến đường khu dân cư Bởi lẽ nơi có độ dốc độ cao lớn người dân tác động đến rừng, công tác trồng rừng, công nghiệp, nông nghiệp khu vực tiến hành Phần lớn hoạt động người tác động vào thực vật mặt đất tiến hành khu vực phằng, độ dốc thấp Đồng thời, hoạt động canh tác ảnh hưởng tới rừng nông nghiệp thường tiến hành với cường độ số lượng lớn khu vực gần đường, gần khu dân 38 Số điểm ảnh phân tích 5493785 5478661 5478661 85790 85497 1096 1096 18023 Giá trị Sig F 0,000 0,000 0,000 0,000 0,000 0,001 0,04 0,3 cư Còn khu vực xa đường, xa khu dân cư thảm thực vật bị tác động hoạt động người Trong đó, yếu tố lượng mưa loại đất khơng có quan hệ tới biến động lượng thực vật che phủ Bởi lẽ, mức độ chênh lệch lượng mưa huyện không lớn phần lớn diện tích tồn huyện đất Ferralit Do vậy, nhân tố không ảnh hưởng đến biến động lớp phủ thực vật Mối quan hệ tương tự kết vấn Mường Chà Sơn La phân tích nhân tố ảnh hưởng tới cháy rừng (Bùi Mạnh Hưng Nguyễn Thanh Thủy Vân, 2019; Nguyễn Ngọc Thạch cộng sự, 2017) KẾT LUẬN Như vậy, từ việc phân tích số thực vật từ ảnh Sentinel năm 2016 2021 địa bàn huyện Yên Châu cho thấy giá trị trung bình số thực vật năm 2016 cao năm 2021 0,1 Điều cho chứng minh rằng, nhìn chung, lượng thực vật có địa bàn huyện bị suy giảm mức độ TẠP CHÍ KHOA HỌC VÀ CÔNG NGHỆ LÂM NGHIỆP SỐ 3- 2021 Lâm học định diện rộng Lượng thực vật rải phạm vi toàn huyện năm 2021 Kết tính tốn cho thấy tổng diện tích khu vực có số thực vật nhỏ giảm từ năm 2016 đến 2021 khoảng 55,5 Đồng thời, tổng diện tích có số thực vật mức độ thấp (từ - 0,6) có xu hướng tăng mạnh, từ khoảng 42 nghìn lên đến 69 nghìn Ngược lại, điều đáng báo động khu vực có lượng thực vật cao (chỉ số thực vật từ 0,6 - 0,9), rừng thường có trữ lượng lớn (Lê Thị Thu Hiền, 2013; Nguyễn Quang Giáp, 2015) lại có xu hướng giảm Điều thể diện tích có lượng thực vật che phủ lớn bị giảm dần Nghiên cứu xây dựng tám nhân tố tự nhiên xã hội để kiểm tra mối quan hệ với số thực vật độ cao, độ dốc, hướng phơi, khoảng cách tới đường giao thông, khoảng cách tới khu dân cư, nhiệt độ trung bình tháng, lượng mưa trung bình năm loại đất Kết phân tích mối tương quan NDVI với nhân tố cho thấy mức độ biến động số thực vật khoảng năm từ 2016 đến 2021 có quan hệ chặt với yếu tố độ dốc hướng phơi, sau đến độ cao, nhiệt độ trung bình khoảng cách đến đường Các nhân tố lượng mưa bình quân tháng loại đất gần khơng có mối quan hệ Điều có nghĩa mức độ suy thối tăng lượng thực vật che phủ chịu ảnh hưởng lớn nhân tố độ dốc, độ cao, khoảng cách đến đường khu dân cư Trong đó, yếu tố lượng mưa loại đất khơng có quan hệ tới biến động lượng thực vật che phủ Bởi lẽ, mức độ chênh lệch lượng mưa huyện khơng lớn phần lớn diện tích tồn huyện đất Ferralit TÀI LIỆU THAM KHẢO Đỗ Trọng Hoàn Nguyễn Hải Vân (2017) Lý thuyết diễn biến rừng số suy nghĩ phát triển lâm nghiệp Việt Nam, Tổ chức người thiên nhiên, Hà Nội, Việt Nam Xem ngày 9/2/2021 Có link sau: https://www.thiennhien.net/2017/03/24/ly-thuyet-vedien-bien-rung-va-mot-suy-nghi-ve-phat-trien-lamnghiep-viet-nam/ Hương Thảo (2010) Tài nguyên rừng nguyên nhân suy thoái rừng Việt Nam, Tổng cục Lâm nghiệp, Hà Nội, Việt Nam Xem ngày 9/2/2021 Có link sau: http://tongcuclamnghiep.gov.vn/LamNghiep/Index/tainguyen-rung-va-nguyen-nhan-suy-thoai-rung-tren-thegioi-195 Hồng Nhung (2018) Suy thoái rừng giới - Thực trạng giải pháp, Ban trung ương mặt trận tổ quốc Việt Nam, Hà Nội, Việt Nam Xem ngày 9/2/2021 Có link sau: http://tapchimattran.vn/thegioi/suy-thoai-rung-tren-the-gioi-thuc-trang-va-giaiphap-14560.html Hương Thảo (2010) Tài nguyên rừng nguyên nhân suy thoái rừng giới, Tổng cục Lâm nghiệp, Hà Nội, Việt Nam Xem ngày 9/2/2021 Có link sau: http://tongcuclamnghiep.gov.vn/LamNghiep/Index/tainguyen-rung-va-nguyen-nhan-suy-thoai-rung-tren-thegioi-194 Lã Nguyên Khang Trần Quang Bảo (2015) Các phương pháp phân tích thống kê đa biến số liệu nghiên cứu lâm nghiệp SAS Tạp chí Nơng nghiệp & PTNT Số 15/2015: 1-12 Ngọc Thuấn (2020) Tập trung làm tốt công tác bảo vệ, phát triển rừng, Báo Sơn La online, Sơn La, Việt Nam Xem ngày 9/2/2021 Có link sau: http://www.baosonla.org.vn/vi/bai-viet/tap-trung-lamtot-cong-tac-bao-ve-phat-trien-rung-34727 UBND Yên Châu (2018) Điều kiện tự nhiên, Cổng thông tin điện tử huyện Yên Châu, Sơn La, Việt Nam Xem ngày 9/2/2021 Có link sau: https://yenchau.sonla.gov.vn/1297/31378/58877/416970 /Dieu-kien-tu-nhien/Dieu-kien-tu-nhien.aspx Bùi Mạnh Hưng Nguyễn Thanh Thủy Vân (2019) Sử dụng QGis phân tích thứ bậc (ahp) để phân cấp nguy cháy rừng huyện Mường Chà, tỉnh Điện Biên Tạp chí Khoa học Cơng nghệ Lâm nghiệp Số 2/2019: 38-47 Nguyễn Hải Hòa, Nguyễn Thị Thu Hiền Lương Thị Thu Trang (2016) Sử dụng GIS ảnh Landsat đa thời gian xây dựng đồ biến động diện tích rừng xã vùng đệm Xuân Đài Kim Thượng, Vườn quốc gia Xuân Sơn Tạp chí Khoa học Lâm nghiệp 3: 4524-4527 10 Nguyễn Trọng Cương, Nguyễn Hải Hòa Trần Quang Bảo (2016) Thành lập đồ rừng ngập mặn ven biển tỉnh Thái Bình từ ảnh Sentinel-2 Tạp chí Khoa học Công nghệ Lâm nghiệp Số 6/2019: 57-66 11 Shivangi S Somvanshi Maya Kumari (2020) Comparative analysis of different vegetation indices with respect to atmospheric particulate pollution using sentinel data Applied Computing and Geosciences 7: 100032 12 Lê Thị Thu Hiền (2013) Áp dụng số thực vật (NDVI) ảnh Landsat đánh giá hoang mạc hóa tỉnh Bình Thuận Tạp chí Khoa học Trái Đất 35(4): 357-363 TẠP CHÍ KHOA HỌC VÀ CƠNG NGHỆ LÂM NGHIỆP SỐ - 2021 39 Lâm học 13 Alex Mandel, Victor Olaya Ferrero, Anita Graser Alexander Bruy (2016) QGIS cookbook Packt Publishing Ltd 14 Nguyễn Quang Giáp (2015) Thử nghiệm phương pháp xây dựng đồ trạng lớp phủ thực vật lưu vực từ ảnh vệ tinh Landsat Tạp chí Khoa học Cơng nghệ Lâm nghiệp 3: 40-50 15 Mai Trọng Thịnh Nguyễn Hải Hòa (2017) Sử dụng ảnh viễn thám đa thời gian đánh giá biến động diện tích rừng ngập mặn thị xã Quảng Yên, tỉnh Quảng Ninh Tạp chí Khoa học Công nghệ Lâm nghiệp 3: 101-112 16 Nguyễn Hải Hòa Nguyễn Văn Quốc (2017) Sử dụng ảnh viễn thám landsat gis xây dựng đồ biến động diện tích rừng vùng đệm Vườn quốc gia Xn Sơn Tạp chí Khoa học Cơng nghệ Lâm nghiệp 3: 46-56 17 Nguyễn Ngọc Thạch, Đặng Ngô Bảo Toàn Phạm Xuân Cảnh (2017) Ứng dụng viễn thám GIS thành lập đồ nguy cháy rừng phục vụ phòng chống, giảm thiểu thiệt hại cháy rừng tỉnh Sơn La, Việt Nam Tạp chí Khoa học ĐHQGHN: Các Khoa học Trái đất Môi trường 33(3): 53-66 VEGETATION AMOUNT CHANGES AND RELATIONS BETWEEN THE CHANGE AND NATURAL SOCIAL FACTORS IN YEN CHAU, SON LA Bui Manh Hung1, Nguyen Thi Bich Phuong1, Nguyen Thi Thao1 Vietnam National University of Forestry SUMMARY Vegetation amount changes on the ground is one of the most concerned issues to solve in many provinces This study used Sentinel images in 2016 and 2021 in Yen Chau district, together with a construction of eight factors to test the relationship of elevation, slope, aspect, distance to roads, distance to residential areas, average temperature, average annual rainfall and soil type Results of the vegetation index (NDVI) calculation showed that the average NDVI value in 2016 was higher than that of 2021 by 0.1 In general, the amount of vegetation available in the district was reduced Total area with NDVI less than decreased from 2016 to 2021 by 55.5 At the same time, total area with low NDVI values (0-0.6) also tended to increase sharply, from about 42,000 to 69,000 On the contrary, the areas with NDVI from 0.6-0.9 tended to decrease Analysis results of the correlation between NDVI and the above factors indicated that the variation in NDVI values from 2016 to 2021 was most strongly related to slope and aspect, the following is elevation, average temperature and the distance to roads Factors such as average monthly rainfall and soil type had almost no relationship with NDVI changes Keywords: natural social factors, Sentinel 2, vegetation amount changes, Yen Chau Ngày nhận Ngày phản biện Ngày định đăng 40 : 02/4/2021 : 07/5/2021 : 17/5/2021 TẠP CHÍ KHOA HỌC VÀ CƠNG NGHỆ LÂM NGHIỆP SỐ 3- 2021 ... cuối mức mạnh 3.3 Mối quan hệ biến động lượng thực vật che phủ nhân tố tự nhiên xã hội Kết phân tích tương quan tuyến tính biến động giá trị số thực vật với yếu tố tự nhiên xã hội thể bảng Bảng... động lượng thực vật với nhân tố tự nhiên xã hội Phương trình hồi quy tuyến tính đa biến sau: ∆ = a0 + a1.Xi Trong đó: ∆ độ lệch giá trị số thực vật hai mốc thời gian; Xi nhân tố tự nhiên xã hội. .. biệt mặt trữ lượng thực vật mặt đất Tuy nhiên, chưa có nghiên cứu cụ thể địa bàn toàn huyện biến động rừng mối quan hệ chúng với nhân tố tự nhiên chưa thực (Ngọc Thuấn, 2020; UBND Yên Châu, 2018)