Tổng quan nghiên cứu

Việt Nam nằm trong nhóm các quốc gia chịu tác động sâu sắc nhất của hiện tượng ấm lên toàn cầu và biến đổi khí hậu cực đoan. Theo ước tính từ các công bố quốc tế về suy thoái đất toàn cầu, khoảng 1.965 triệu ha đất trên thế giới đang bị thoái hóa ở các mức độ khác nhau, trong đó các đánh giá chuyên biệt tại Việt Nam chỉ ra tỷ lệ đất bị thoái hóa dao động từ 12,6% đến 26,8% tổng diện tích tự nhiên. Sự thay đổi bất thường của chế độ mưa theo không gian và thời gian phá vỡ quy luật sinh thái tự nhiên, trực tiếp dẫn đến sự suy giảm liên tục của thảm thực vật – dấu hiệu chỉ thị sớm của quá trình thoái hóa đất đai và hoang mạc hóa. Vấn đề suy giảm sức sản xuất của đất gây thiệt hại trực tiếp đến sinh kế của hơn 94 triệu dân và đe dọa an ninh lương thực quốc gia.

Nghiên cứu được thực hiện với mục tiêu trọng tâm là đánh giá định lượng biến động lớp phủ thực vật và chế độ mưa trong chuỗi thời gian 33 năm liên tục từ năm 1982 đến năm 2014 trên toàn bộ lãnh thổ đất liền Việt Nam. Đề tài xác định hệ số tương quan Pearson giữa lượng mưa trung bình năm và chỉ số thực vật chuẩn hóa, từ đó phân tích và khoanh vùng các khu vực có nguy cơ thoái hóa đất dài hạn do khí hậu hoặc tác động nhân sinh.

Phạm vi không gian nghiên cứu bao phủ toàn bộ diện tích đất liền của 7 vùng khí hậu đặc trưng, bao gồm Tây Bắc, Việt Bắc - Đông Bắc, Đồng bằng Bắc Bộ, Bắc Trung Bộ, Nam Trung Bộ, Tây Nguyên và Nam Bộ. Về mặt khoa học và thực tiễn, công trình cung cấp cơ sở dữ liệu số chuẩn hóa và bản đồ chuyên đề hỗ trợ trực tiếp cho các nhà hoạch định chính sách trong công tác quy hoạch sử dụng đất, giám sát 2.360 lưu vực sông chính và bảo tồn hơn 14.600 loài thực vật nhiệt đới trên cả nước.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Nghiên cứu ứng dụng khung lý thuyết toàn diện về suy thoái đất đai theo định nghĩa chuẩn của Công ước Liên Hợp Quốc về Chống Sa mạc hóa (UNCCD) năm 2004 và Tổ chức Nông nghiệp và Lương thực Liên Hợp Quốc (FAO) năm 2000. Theo đó, thoái hóa đất đai được xác định là sự suy giảm hoặc mất đi liên tục các dịch vụ hệ sinh thái đất, đặc biệt là năng suất sơ cấp thuần (NPP). Bên cạnh đó, đề tài tích hợp mô hình xói mòn đất phổ dụng cải biên (RUSLE), trong đó hệ số che phủ thực vật C được xác định trực tiếp thông qua chỉ số thực vật viễn thám theo công thức chuyên biệt để ước tính mức độ xói mòn đất hàng năm do dòng chảy bề mặt.

Bốn khái niệm cốt lõi được vận dụng xuyên suốt công trình bao gồm:

  • Chỉ số sai khác thực vật chuẩn hóa (NDVI): Đại lượng quang phổ phản ánh mật độ sinh khối và sức khỏe của thảm thực vật dựa trên tỷ lệ phản xạ giữa kênh cận hồng ngoại và kênh ánh sáng đỏ.
  • Năng suất sơ cấp thuần (NPP): Lượng carbon hữu cơ thảm thực vật cố định được từ khí quyển sau khi trừ đi lượng tiêu hao qua hô hấp, có tương quan tuyến tính chặt chẽ với NDVI chuỗi thời gian.
  • Hệ số tương quan tuyến tính Pearson (r): Thước đo định lượng mối liên hệ đồng biến hoặc nghịch biến giữa chuỗi số liệu lượng mưa và chỉ số sinh khối thực vật.
  • Phân vùng khí hậu tự nhiên: Phân chia 7 vùng khí hậu đặc thù dựa trên chế độ nhiệt ẩm, bức xạ mặt trời (105 đến 170 kcal/cm2/năm) và thời gian xuất hiện mùa mưa.

Phương pháp nghiên cứu

Nguồn dữ liệu thứ cấp của nghiên cứu bao gồm bộ dữ liệu ảnh vệ tinh VIPHEN_EVI2 giai đoạn 1982 – 2014 từ dự án của Đại học Arizona do NASA MEaSUREs tài trợ, độ phân giải không gian 5.600 m với 26 lớp dữ liệu phổ kế thừa từ cảm biến AVHRR trên các vệ tinh NOAA-07, NOAA-09, NOAA-11, NOAA-14 và cảm biến MODIS Terra. Dữ liệu lượng mưa vệ tinh CHIRPS độ phân giải 0,05 độ (khoảng 5,6 km) trong 33 năm được trích xuất và xử lý tự động qua nền tảng điện toán đám mây Google Earth Engine. Ranh giới hành chính các cấp được trích xuất từ cơ sở dữ liệu không gian GADM.

Cỡ mẫu không gian bao phủ toàn bộ lưới pixel trên 330.000 km2 lãnh thổ đất liền Việt Nam, với chuỗi quan sát thời gian gồm 33 mốc hàng năm (n = 33). Phương pháp chọn mẫu toàn diện theo mạng lưới ô vuông không gian đảm bảo tính đại diện và loại bỏ hoàn toàn sai số chọn mẫu cục bộ.

Phương pháp phân tích hồi quy tuyến tính bình phương nhỏ nhất (OLS) được lựa chọn để xác định hệ số góc A biểu thị xu thế biến động NDVI theo thời gian, kết hợp kiểm định giả thuyết t-Student ở mức ý nghĩa thống kê p < 0,10 với bậc tự do df = 31 (giá trị t tới hạn 1,309). Mối tương quan giữa lượng mưa và NDVI được phân tích qua hệ số Pearson r và chuẩn hóa phân phối Fisher Z trên công cụ Raster Calculator của phần mềm ArcGIS. Lý do lựa chọn tổ hợp phương pháp này là nhằm chuẩn hóa phân tích chuỗi thời gian dài hạn trên quy mô quốc gia, khắc phục triệt để nhược điểm tốn kém kinh phí và tính chủ quan của các phương pháp điều tra thổ nhưỡng truyền thống.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Quá trình xử lý và phân tích thống kê chuỗi dữ liệu 33 năm (1982 – 2014) mang lại 4 phát hiện quan trọng:

Thứ nhất, phân bố không gian chỉ số NDVI trung bình 33 năm thể hiện sự phân hóa sinh thái sâu sắc giữa các vùng địa hình. Dải rừng nguyên sinh và rừng đầu nguồn dọc biên giới Việt – Lào trải dài từ Hòa Bình, chạy qua Bắc Trung Bộ đến Quảng Nam và khu vực Tây Nguyên (Kon Tum, Lâm Đồng) duy trì giá trị NDVI rất cao, thường xuyên vượt ngưỡng 0,75. Ngược lại, hai vùng đồng bằng châu thổ trũng thấp là Đồng bằng Bắc Bộ (rộng 16.700 km2) và Đồng bằng Nam Bộ (rộng 40.000 km2) cùng dải đất hẹp duyên hải ven biển chỉ đạt giá trị NDVI trung bình dưới 0,40.

Thứ hai, xu thế biến động tuyến tính của thảm thực vật cho thấy sự gia tăng sinh khối rõ rệt tại một số vùng rừng phòng hộ miền núi phía Bắc sau các chương trình phục hồi rừng tự nhiên. Tuy nhiên, tại các khu vực đô thị hóa nhanh và vùng chuyển đổi đất nông nghiệp chuyên canh, chỉ số NDVI sụt giảm đáng kể với mức giảm sinh khối từ 15% đến 28% qua các giai đoạn quan sát.

Thứ ba, hệ số tương quan Pearson giữa lượng mưa trung bình năm (dao động từ 1.200 đến 2.000 mm) và chỉ số NDVI thể hiện mối quan hệ thuận chiều chặt chẽ (|r| > 0,45) tại các vùng khô hạn cục bộ như Nam Trung Bộ và các cao nguyên Tây Nguyên. Tại những nơi này, sự gia tăng lượng mưa kích thích trực tiếp sự phục hồi và phát triển của thảm thực vật.

Thứ tư, nghiên cứu đã tách biệt thành công các điểm nóng suy giảm NDVI do tác động tiêu cực của mưa lớn cực đoan gây xói mòn cơ học với các điểm nóng suy thoái do áp lực con người (khai thác gỗ, chuyển đổi mục đích sử dụng đất) với mức ý nghĩa thống kê p < 0,10.

Thảo luận kết quả

Tại các khu vực có số giờ nắng cao từ 2.000 đến 2.500 giờ/năm như Nam Trung Bộ và Tây Nguyên, độ ẩm là yếu tố sinh thái giới hạn hàng đầu đối với sự phát triển của hệ thực vật. Khi lượng mưa phân bổ thuận lợi, sinh khối thực vật phản hồi tích cực và gia tăng chỉ số quang phổ. Trái lại, tại các vùng địa hình dồi núi hiểm trở chiếm 3/4 diện tích lãnh thổ (đặc biệt là sườn đông dãy Trường Sơn và vùng núi Tây Bắc có độ cao trên 1.000 m), lượng mưa tập trung cao độ vào mùa lũ (chiếm 70% đến 80% tổng dòng chảy năm) lại là nguyên nhân chính thúc đẩy quá trình rửa trôi, xói mòn tầng đất mặt và trượt lở đất, làm phá hủy cấu trúc thảm thực vật tự nhiên.

Kết quả này hoàn toàn tương đồng với các phát hiện toàn cầu của NASA MEaSUREs và nghiên cứu thoái hóa đất tiểu vùng sông Mê Kông của Shrestha và Roy năm 2008. Tuy nhiên, công trình của tác giả Nguyễn Dân Trí mang lại độ chi tiết vượt trội so với các đánh giá ASSOD trước đây nhờ sử dụng ảnh vệ tinh phân giải liên tục, giúp định vị chính xác tọa độ các điểm nóng thoái hóa.

Về mặt trực quan hóa, toàn bộ hệ thống dữ liệu nghiên cứu được biểu diễn thông qua chuỗi 15 bản đồ chuyên đề ArcGIS (bản đồ lượng mưa CHIRPS, bản đồ NDVI 33 năm, bản đồ tương quan r Pearson, bản đồ suy giảm sinh khối p < 0,10) và các đồ thị hồi quy tuyến tính, cho phép các cơ quan chuyên môn theo dõi trực quan diễn biến biến động đất đai qua từng thời kỳ lịch sử.

Đề xuất và khuyến nghị

Dựa trên các kết quả phân tích thực nghiệm, 4 nhóm giải pháp chiến lược được đề xuất nhằm nâng cao hiệu quả quản lý tài nguyên đất và thích ứng với biến đổi khí hậu:

  • Thiết lập hệ thống giám sát và cảnh báo sớm suy thoái đất bằng công nghệ viễn thám tự động: Ứng dụng thuật toán phân tích chuỗi thời gian trên Google Earth Engine nhằm giảm 25% thời gian phát hiện điểm nóng suy thoái lớp phủ thực vật tại 7 vùng khí hậu. Mục tiêu hoàn thiện và chuyển giao công nghệ trong giai đoạn 2026 – 2028 dưới sự chủ trì của Cục Viễn thám Quốc gia phối hợp với Viện Thổ nhưỡng Nông hóa.
  • Tái cơ cấu cây trồng và phục hồi rừng đầu nguồn tại các khu vực nhạy cảm khí hậu: Triển khai các mô hình nông lâm kết hợp, trồng rừng đa tầng tán tại các vùng khô hạn Nam Trung Bộ và Tây Nguyên nhằm nâng tỷ lệ che phủ rừng thực chất lên trên 60% vào năm 2030, do Bộ Nông nghiệp và Phát triển Nông thôn chỉ đạo thực hiện.
  • Cập nhật định kỳ bản đồ phân vùng nguy cơ thoái hóa đất do mưa lũ cực đoan: Thực hiện điều tra định lượng theo chuẩn Thông tư 14/2012/TT-BTNMT với tần suất 5 năm một lần cho 2.360 hệ thống sông ngòi, giảm thiểu 20% thiệt hại xói mòn đất nông nghiệp, giao Sở Tài nguyên và Môi trường các địa phương chủ trì tổ chức.
  • Hoàn thiện cơ chế chính sách chi trả dịch vụ môi trường rừng gắn với lưu giữ carbon: Gia tăng nguồn thu bổ sung cho cộng đồng tham gia bảo vệ rừng đầu nguồn thêm 30% trong lộ trình 2026 – 2030, do Quỹ Bảo vệ và Phát triển Rừng Việt Nam điều phối và giám sát thực hiện.

Đối tượng nên tham khảo luận văn

Công trình luận văn là tài liệu khoa học chuyên sâu và mang tính ứng dụng thực tiễn cao, đặc biệt hữu ích cho 4 nhóm đối tượng sau:

  • Học viên cao học và nghiên cứu sinh chuyên ngành Biến đổi khí hậu, Địa lý tài nguyên, GIS và Viễn thám: Khai thác quy trình xử lý dữ liệu VIPHEN_EVI2 và CHIRPS chuỗi thời gian dài trên 30 năm làm khung tham chiếu phương pháp luận cho các đề tài nghiên cứu định lượng về hệ sinh thái và khí hậu.
  • Chuyên viên kỹ thuật và lãnh đạo tại Cục Biến đổi khí hậu, Cục Quản lý Đất đai (Bộ Tài nguyên và Môi trường): Sử dụng bản đồ phân vùng tương quan mưa và suy thoái sinh khối làm cơ sở dữ liệu đầu vào cho việc xây dựng Quy hoạch không gian biển và đất đai quốc gia tầm nhìn 2050.
  • Cán bộ quản lý nông nghiệp, kiểm lâm và quy hoạch phát triển nông thôn cấp tỉnh: Tham khảo số liệu biến động lớp phủ thực vật thực tế để thiết kế các giải pháp bố trí mùa vụ, giống cây trồng và phương án phòng chống sạt lở đất trong mùa mưa lũ.
  • Chuyên gia dự án tại các tổ chức phi chính phủ và cơ quan quốc tế (UNDP, FAO, UNCCD): Khai thác dữ liệu 33 năm để xây dựng các báo cáo quốc gia về thoái hóa đất, đồng thời lập đề xuất tài trợ cho các dự án phục hồi cảnh quan rừng và thích ứng biến đổi khí hậu tại dải ven biển miền Trung và Tây Nguyên.

Câu hỏi thường gặp

Tại sao nghiên cứu lại chọn chuỗi dữ liệu 33 năm từ năm 1982 đến năm 2014? Theo quy chuẩn của Tổ chức Khí tượng Thế giới (WMO), một chu kỳ chuẩn khí hậu cần tối thiểu 30 năm quan sát liên tục. Khoảng thời gian 33 năm (1982 – 2014) đáp ứng đầy đủ tiêu chí này, cho phép đánh giá chính xác quy luật biến đổi dài hạn của lượng mưa và lớp phủ thực vật, loại trừ các nhiễu động thời tiết ngắn hạn trong từng năm đơn lẻ.

Chỉ số sai khác thực vật NDVI phản ánh hiện tượng thoái hóa đất như thế nào? Theo khung lý thuyết của UNCCD và FAO, thoái hóa đất biểu hiện rõ nhất qua sự sụt giảm bền vững của năng suất sơ cấp thuần (NPP). Do NDVI có mối tương quan tuyến tính chặt chẽ với sinh khối xanh, việc suy giảm liên tục giá trị NDVI trên cùng một đơn vị không gian trong nhiều thập kỷ là bằng chứng xác thực về sự suy thoái độ phì và chức năng hệ sinh thái đất.

Vì sao lượng mưa tăng tại các vùng núi cao lại có thể làm giảm chỉ số thực vật NDVI? Tại các vùng đồi núi dốc trên 15 độ ở Tây Bắc hoặc miền Trung, lượng mưa hàng năm tập trung lớn với cường độ mưa cực đoan trong thời gian ngắn sẽ tạo dòng chảy mặt mạnh. Hiện tượng này cuốn trôi lớp đất mặt màu mỡ theo phương trình mất đất RUSLE, gây xói mòn sâu, trượt lở đất và làm bật gốc thảm thực vật, dẫn đến giảm mật độ sinh khối quang phổ.

Ưu điểm nổi bật của việc kết hợp dữ liệu CHIRPS và VIPHEN_EVI2 so với trạm quan trắc mặt đất là gì? Hệ thống trạm đo khí tượng truyền thống tại Việt Nam phân bố thưa thớt, không đồng đều ở vùng sâu vùng xa. Dữ liệu viễn thám CHIRPS và VIPHEN_EVI2 cung cấp thông tin không gian dạng lưới pixel liên tục với độ phân giải khoảng 5,6 km, đảm bảo tính đồng nhất trên toàn bộ 330.000 km2 lãnh thổ và giảm thiểu 90% chi phí điều tra thực địa.

Làm thế nào để phân biệt suy giảm thực vật do khí hậu với suy giảm do tác động của con người? Nghiên cứu sử dụng phân tích tương quan Pearson kết hợp kiểm định thống kê Fisher Z. Những khu vực có NDVI suy giảm mạnh nhưng hệ số tương quan với lượng mưa không có ý nghĩa thống kê (|r| < 0,45, p > 0,10) được tách biệt thành nhóm suy thoái do nguyên nhân nhân sinh (phá rừng, khai thác khoáng sản, đô thị hóa) thay vì do biến động thời tiết.

Kết luận

  • Luận văn đã xây dựng thành công bộ cơ sở dữ liệu không gian tích hợp chỉ số NDVI và lượng mưa CHIRPS trong 33 năm (1982 – 2014) trên toàn bộ 7 vùng khí hậu của Việt Nam.
  • Kết quả nghiên cứu chứng minh sự phân hóa rõ rệt của lớp phủ thực vật với giá trị NDVI vượt trên 0,75 tại các dải rừng biên giới và dưới 0,40 tại các vùng đồng bằng châu thổ.
  • Đề tài xác định được mối tương quan tuyến tính chặt chẽ (|r| > 0,45, p < 0,10) giữa lượng mưa và sinh khối tại các vùng khô hạn trọng điểm như Nam Trung Bộ và Tây Nguyên.
  • Phương pháp tiếp cận viễn thám đa thời gian đã tách biệt hiệu quả các điểm nóng thoái hóa đất do xói mòn mưa lũ với các khu vực chịu áp lực từ hoạt động phát triển kinh tế của con người.
  • Công trình đóng góp quy trình phương pháp luận liên ngành hiện đại, tạo cơ sở khoa học tin cậy để triển khai các cam kết quốc tế về chống sa mạc hóa và hỗ trợ thực thi Thông tư 14/2012/TT-BTNMT.

Trong giai đoạn tiếp theo từ năm 2026 đến năm 2030, nghiên cứu cần được tiếp tục mở rộng bằng cách tích hợp thêm các biến số nhiệt độ bề mặt và độ ẩm đất từ các thế hệ vệ tinh mới. Các nhà khoa học, cơ quan quản lý và địa phương hãy ứng dụng ngay nguồn tư liệu viễn thám chuẩn hóa này vào công tác quản trị đất đai nhằm bảo vệ bền vững tài nguyên thiên nhiên quốc gia.