Tổng quan nghiên cứu

Thoái hóa tài nguyên đất do xói mòn là thách thức môi trường đặc biệt nghiêm trọng tại các vùng đồi núi phía Bắc Việt Nam. Tại huyện Quản Bạ, tỉnh Hà Giang – khu vực vùng cao biên giới thuộc Công viên Địa chất Toàn cầu Cao nguyên đá Đồng Văn có tọa độ địa lý từ 22°57’ đến 23°10’ vĩ độ Bắc và 104°40’30” đến 105° kinh độ Đông – địa hình bị chia cắt hiểm trở với phần lớn diện tích có độ dốc trên 25 độ. Áp lực gia tăng dân số cùng tập quán canh tác ngô độc canh trên đất dốc đã khiến lớp đất mặt màu mỡ bị rửa trôi nhanh chóng. Hơn 70% đến 85% tổng lượng mưa trong năm tập trung dồn dập vào mùa mưa (từ tháng 5 đến tháng 10), kết hợp với độ che phủ rừng còn hạn chế ở nhiều sườn đồi, tạo nên động năng dòng chảy bề mặt cực lớn, thường xuyên gây lũ quét, sạt lở và làm suy giảm nghiêm trọng quỹ đất canh tác vốn chỉ chiếm tỷ lệ nhỏ toàn huyện.

Nhằm giải quyết bài toán cấp bách này, luận văn thạc sĩ khoa học chuyên ngành Khoa học môi trường của tác giả Nguyễn Thị Mai Hương (Trường Đại học Khoa học Tự nhiên – Đại học Quốc gia Hà Nội, dưới sự hướng dẫn của PGS. Lê Văn Thiện, hoàn thành năm 2015) đã thực hiện đề tài nghiên cứu định lượng nguy cơ xói mòn đất và đề xuất mô hình nông nghiệp thích ứng. Nghiên cứu hướng tới ba mục tiêu cụ thể: xác định các thông số xói mòn chuyên biệt cho vùng núi đá cao Quản Bạ; thành lập bản đồ nguy cơ xói mòn đất tỷ lệ 1:25.000 bằng công nghệ GIS; và xây dựng giải pháp canh tác nông nghiệp bền vững. Kết quả nghiên cứu mang ý nghĩa thực tiễn to lớn khi cung cấp cơ sở dữ liệu chuẩn xác để chính quyền địa phương kiểm soát lượng mất đất trung bình hàng năm, hướng tới mục tiêu duy trì xói mòn dưới ngưỡng cho phép 5 tấn/ha/năm, đồng thời nâng cao hiệu quả kinh tế trên 1 đơn vị diện tích canh tác cho đồng bào dân tộc thiểu số.

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 được xây dựng dựa trên nền tảng lý thuyết động lực học xói mòn do nước và Phương trình Mất đất Phổ dụng (USLE – Universal Soil Loss Equation) do Wischmeier và Smith phát triển năm 1978. Mô hình USLE là công cụ kinh nghiệm chuẩn hóa toàn cầu để lượng hóa lượng đất mất trung bình hàng năm theo công thức toán học: A = R × K × L × S × C × P.

Trong đó, năm nhóm nhân tố cốt lõi được lượng hóa đồng bộ:

  1. Hệ số xói mòn do mưa (R): Thước đo động năng và sức mạnh dòng chảy tràn bề mặt do mưa tạo ra.
  2. Hệ số mẫn cảm xói mòn của đất (K): Biểu thị tính dễ bị tổn thương và khả năng chống chịu rửa trôi của từng loại đất, phụ thuộc vào thành phần cấp hạt, cấu trúc đoàn lạp và hàm lượng mùn hữu cơ.
  3. Hệ số địa hình (LS): Đại lượng tích hợp chiều dài sườn dốc (L) và góc dốc (S), phản ánh mức độ gia tăng vận tốc dòng chảy mặt khi địa hình dốc đứng.
  4. Hệ số lớp phủ thực vật (C): Tỷ số đánh giá mức độ che chắn bề mặt của thảm thực bì tự nhiên hoặc cây trồng nông nghiệp nhằm hấp thu động năng hạt mưa.
  5. Hệ số biện pháp kỹ thuật chống xói mòn (P): Tỷ số phản ánh hiệu quả giảm dòng chảy của các công trình cơ học như ruộng bậc thang, rãnh phân lũ hoặc đường đồng mức.

Bên cạnh mô hình USLE, đề tài kế thừa lý thuyết phân cấp độ dốc trong quy hoạch nông nghiệp của Viện Quy hoạch và Thiết kế Nông nghiệp, kết hợp các quan điểm sinh thái học cảnh quan về mô hình canh tác nông lâm kết hợp bền vững trên đất dốc (SALT).

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

Để thực hiện mô hình hóa không gian, nghiên cứu thu thập và xử lý hệ thống dữ liệu đa nguồn: số liệu quan trắc khí tượng thủy văn nhiều năm tại các trạm đo mưa khu vực; bản đồ thổ nhưỡng tỷ lệ 1:25.000 do Viện Quy hoạch và Thiết kế Nông nghiệp thành lập; bản đồ hiện trạng sử dụng đất và hiện trạng rừng huyện Quản Bạ năm 2014; cùng dữ liệu mô hình số độ cao (Digital Elevation Model – DEM).

Quy trình phân tích được thực hiện trên phần mềm ArcGIS 10.0 thông qua các bước kỹ thuật:

  • Tính toán hệ số R: Áp dụng phương trình hồi quy thực nghiệm dựa trên tổng lượng mưa trung bình năm và phân bố mưa ngày.
  • Tính toán hệ số K: Xác định theo bảng chỉ số xói mòn thực nghiệm của các loại đất chính ở Việt Nam theo công thức ISSS và toán đồ chuẩn, phân loại theo thang Katrinski.
  • Tính toán hệ số LS: Xây dựng lưới hướng dòng chảy (Flow Direction) và tích lũy sườn dốc từ DEM, áp dụng thuật toán Toxopeus cho sườn dốc dưới 20 độ và phương trình Gaudasasmita cho sườn dốc trên 20 độ.
  • Tính toán hệ số C và P: Giải đoán lớp phủ thực vật và hiện trạng canh tác, gán trọng số C*P chi tiết cho từng loại hình rừng và cây trồng.
  • Chồng lớp không gian (Raster Overlay): Tích hợp 4 lớp bản đồ đơn tính để kết xuất bản đồ xói mòn đất A, phân thành 7 cấp độ nguy cơ từ rất yếu đến mãnh liệt.

Nghiên cứu áp dụng phương pháp điều tra mẫu bao phủ toàn bộ 13 đơn vị hành chính xã, thị trấn của huyện Quản Bạ kết hợp khảo sát thực địa tại các điểm nóng xói mòn, sườn đồi dốc cao và các khu vực ruộng bậc thang dọc thung lũng sông Miện. Việc lựa chọn công nghệ GIS kết hợp USLE đảm bảo tính chính xác cao, rút ngắn thời gian xử lý dữ liệu không gian phức tạp và tối ưu hóa chi phí nghiên cứu so với các phương pháp đo đạc ô tiêu chuẩn đơn lẻ.

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

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

Quá trình mô hình hóa và phân tích dữ liệu không gian mang lại những phát hiện quan trọng về thực trạng xói mòn đất tại huyện Quản Bạ:

  • Phân hóa rõ nét của hệ số xói mòn do mưa (R): Lượng mưa tập trung cao độ trong 6 tháng mùa mưa (chiếm 70% đến 85% tổng lượng mưa năm). Các trận mưa rào cường độ lớn trên 25 mm/h với vận tốc hạt mưa đường kính 5 mm đạt tới 9 m/s tạo áp lực bóc tách bề mặt rất mạnh, khiến hệ số R đạt giá trị cao nhất tại các xã vùng cao đón gió ẩm.
  • Hệ số mẫn cảm xói mòn đất (K) dao động mạnh: Đất mùn vàng trên núi đỏ có hệ số K thấp khoảng 0,15 đến 0,16 nhờ cấu trúc đoàn lạp bền vững và hàm lượng chất hữu cơ cao. Ngược lại, nhóm đất feralit đỏ vàng trên đá sét và đất xám feralit có hệ số K lên đến 0,20 đến 0,23, thể hiện mức độ dễ bị rửa trôi tầng mặt khi mất lớp phủ sinh học.
  • Tác động sâu sắc của hệ số canh tác và quản lý (CP): Rừng giàu tự nhiên có hệ số CP tối ưu ở mức 0,00047, trong khi đất trồng cây hàng năm trên nương dốc từ 15 đến 20 độ có hệ số CP lên tới 0,07800 (cao gấp hơn 165 lần so với rừng giàu). Đất trống có cây bụi và đất chăn thả gia súc có hệ số CP dao động từ 0,01500 đến 0,07500.
  • Sự tập trung diện tích nguy cơ xói mòn cao: Phân vùng nguy cơ xói mòn cho thấy diện tích mất đất cấp độ mạnh đến rất mạnh (>50 đến trên 100 tấn/ha/năm) phân bố chủ yếu ở các sườn núi dốc trên 25 độ bị mất rừng tự nhiên và chuyển sang trồng ngô, dong riềng không có công trình giữ nước. Ngược lại, vùng xói mòn rất yếu (<5 tấn/ha/năm) tập trung tại các cánh rừng đặc dụng và khu vực ruộng bậc thang ổn định.

Thảo luận kết quả

Nguyên nhân chính dẫn đến hiện tượng xói mòn nghiêm trọng tại Quản Bạ xuất phát từ sự kết hợp bất lợi giữa yếu tố tự nhiên và nhân sinh. Địa hình karst chia cắt mạnh với nhiều sườn dốc trên 25 độ làm gia tăng năng lượng dòng chảy mặt gấp 4 lần khi tốc độ dòng chảy tăng gấp đôi. Trong khi đó, thói quen canh tác nương rẫy truyền thống, phát dọn thực bì triệt để đầu mùa mưa và độc canh cây ngô đã để lộ mặt đất trần đúng vào thời điểm các trận mưa rào tháng 5 - tháng 7 diễn ra dồn dập.

So sánh với các nghiên cứu thực nghiệm tại vùng núi Tây Bắc và Tây Nguyên của các nhà khoa học trước đây, lượng mất đất trên sườn dốc 15 độ trồng cây ngắn ngày không có biện pháp che phủ có thể lên tới 200 - 300 tấn/ha/năm, làm trôi mất hàng chục tấn chất hữu cơ và dinh dưỡng N, P, K mỗi héc-ta. Kết quả nghiên cứu tại Quản Bạ hoàn toàn tương thích với các quy luật động lực học đất này. Toàn bộ các phát hiện được thể hiện trực quan qua bản đồ nguy cơ xói mòn đất 7 cấp màu chuẩn hóa (từ màu xanh lá mạ cho cấp rất yếu đến màu da cam và đỏ thẫm cho các cấp xói mòn nghiêm trọng), đi kèm bảng thống kê diện tích cụ thể của từng cấp độ, tạo công cụ trực quan hỗ trợ đắc lực cho công tác quy hoạch lãnh thổ.

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

Nhằm giảm thiểu thất thoát tài nguyên đất và thiết lập nền nông nghiệp bền vững tại Quản Bạ, luận văn đề xuất 4 nhóm giải pháp đồng bộ:

  1. Ứng dụng kỹ thuật canh tác nông lâm kết hợp trên đất dốc (SALT và SALT2): Triển khai mô hình SALT (Kỹ thuật canh tác đất dốc) và SALT2 (Nông lâm súc kết hợp) trên các diện tích đất dốc từ 15 đến 25 độ. Thiết lập các băng rào chắn cây họ đậu phân xanh (như muồng hoa vàng, cỏ vertiver, keo dậu) theo đường đồng mức với khoảng cách 4 đến 6 mét một băng. Biện pháp này giúp giảm tốc độ dòng chảy mặt trên 60%, giữ lại hơn 80% lượng bùn cát rửa trôi và bổ sung lượng lớn sinh khối mùn hữu cơ cho đất. (Chủ thể thực hiện: Phòng Nông nghiệp & PTNT huyện Quản Bạ phối hợp cùng các Hợp tác xã; Mục tiêu: Áp dụng trên 1.500 ha đất nương rẫy; Lộ trình: 2016 – 2020).

  2. Khoanh nuôi tái sinh và tăng cường bảo vệ rừng đầu nguồn: Dừng triệt để việc khai phá nương rẫy mới trên sườn dốc trên 25 độ theo đúng quy chuẩn lâm nghiệp. Chuyển hóa toàn bộ đất trống đồi trọc dốc cao sang đất lâm nghiệp thông qua khoanh nuôi tái sinh tự nhiên và trồng rừng hỗn giao bản địa. Động thái này giúp duy trì hệ số bảo vệ đất C dưới ngưỡng 0,00235 và nâng tỷ lệ che phủ rừng toàn huyện đạt trên 55%. (Chủ thể thực hiện: Hạt Kiểm lâm Quản Bạ và Ban Quản lý rừng phòng hộ; Lộ trình: Thực hiện thường xuyên giai đoạn 2016 – 2025).

  3. Xây dựng công trình cơ học hỗ trợ chống xói mòn: Quy hoạch mở rộng và kiên cố hóa hệ thống ruộng bậc thang tại các khu vực sườn dốc thoải có nguồn nước tự chảy; xây dựng bờ kè đá chân nương, rãnh tiêu nước phân tán đỉnh dốc nhằm triệt tiêu sức xói mòn dạng mương xói. Mục tiêu đưa lượng mất đất hàng năm tại các vùng chuyên canh xuống dưới mức an toàn 10 tấn/ha/năm. (Chủ thể thực hiện: UBND các xã và cộng đồng dân cư địa phương; Lộ trình: 2016 – 2022).

  4. Đổi mới cơ cấu cây trồng và tập huấn kỹ thuật canh tác che phủ: Chuyển đổi các diện tích ngô độc canh kém hiệu quả sang trồng cây lâu năm có tán che rộng kết hợp cây dược liệu bản địa (thảo quả, atisô, tam thất) dưới tán rừng. Tổ chức các lớp tập huấn kỹ thuật cày bừa theo đường đồng mức, giữ lại tàn dư thực vật che phủ gốc và ủ phân hữu cơ cho ít nhất 80% hộ nông dân tại các xã trọng điểm. (Chủ thể thực hiện: Trạm Khuyến nông huyện Quản Bạ; Lộ trình: 2016 – 2018).

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

Công trình nghiên cứu cung cấp giá trị học thuật và thực tiễn chuyên sâu cho 4 nhóm đối tượng chính:

  • Cơ quan quản lý Nhà nước và chính quyền địa phương: UBND huyện Quản Bạ, Sở Tài nguyên và Môi trường, Sở Nông nghiệp và PTNT tỉnh Hà Giang sử dụng hệ thống bản đồ số và cơ sở dữ liệu xói mòn để quy hoạch phân vùng sử dụng đất, phân bổ chỉ tiêu cây trồng và hoạch định chiến lược phòng chống thiên tai, giảm thiểu rủi ro lũ quét sạt lở.
  • Chuyên gia quy hoạch nông – lâm nghiệp và cán bộ địa chính xã: Ứng dụng các thông số kỹ thuật về độ dốc, hệ số LS, hệ số K và mô hình SALT để thiết kế phương án giao đất giao rừng, chuyển đổi cơ cấu kinh tế nông nghiệp thích ứng với điều kiện địa hình miền núi hiểm trở.
  • Giảng viên, nghiên cứu sinh và học viên cao học chuyên ngành Môi trường, Địa lý, Thổ nhưỡng: Sử dụng công trình như một tài liệu tham khảo mẫu mực về phương pháp tích hợp mô hình USLE trong môi trường GIS để định lượng suy thoái đất tại các khu vực có địa hình phức tạp, thiếu hụt mạng lưới trạm đo vi khí hậu.
  • Các tổ chức phi chính phủ (NGOs) và ban quản lý dự án phát triển nông thôn: Khai thác kết quả đề tài để xây dựng các dự án tài trợ phát triển sinh kế bền vững, bảo tồn tài nguyên đất và nguồn nước đầu nguồn gắn liền với mục tiêu xóa đói giảm nghèo cho đồng bào vùng cao.

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

  1. Mô hình USLE kết hợp GIS mang lại ưu thế vượt trội nào trong đánh giá xói mòn tại Quản Bạ? Sự kết hợp này cho phép định lượng hóa lượng đất mất theo không gian phân tán liên tục thay vì chỉ tính toán trên các điểm quan trắc cục bộ. Thông qua công nghệ xử lý ảnh và GIS tỷ lệ 1:25.000, nghiên cứu đã bóc tách chính xác 4 lớp nhân tố (R, K, LS, C*P) cho toàn bộ diện tích huyện, giúp trực quan hóa các điểm nóng xói mòn với độ tin cậy khoa học cao và chi phí tối ưu.

  2. Tại sao các khu vực có độ dốc trên 25 độ không được bố trí sản xuất cây nông nghiệp hàng năm? Theo quy chuẩn khoa học đất, khi độ dốc vượt 25 độ, động năng dòng chảy bề mặt tăng theo hàm mũ, khiến lượng đất mất hàng năm có thể vượt ngưỡng 100 đến 200 tấn/ha/năm nếu canh tác cây ngắn ngày. Việc xới xáo đất trên sườn dốc đứng phá hủy hoàn toàn tầng thổ nhưỡng trong 2 đến 3 vụ canh tác, biến đất thành núi trọc và gây sạt lở đất nghiêm trọng.

  3. Hệ số canh tác và quản lý (CP) đóng vai trò quyết định ra sao trong việc giảm xói mòn? Hệ số CP phản ánh mức độ can thiệp trực tiếp của con người. Kết quả nghiên cứu cho thấy rừng tự nhiên với độ che phủ dày có hệ số C*P cực thấp (0,00047), trong khi đất trồng hoa màu nương rẫy có hệ số lên đến 0,07800 (cao gấp hơn 165 lần). Điều này chứng minh việc duy trì thảm phủ và áp dụng biện pháp cơ học có thể triệt tiêu đến 95% nguy cơ xói mòn đất.

  4. Mô hình canh tác SALT mang lại lợi ích cụ thể gì cho đồng bào huyện Quản Bạ? Mô hình SALT sử dụng các băng cây họ đậu đa niên trồng dày theo đường đồng mức để tạo thành các bờ chắn tự nhiên. Mô hình này giúp giữ lại trên 80% lượng đất mặt giàu dinh dưỡng bị trôi, duy trì độ ẩm đất trong mùa khô và cung cấp hàng tấn phân xanh hữu cơ mỗi năm, giúp cây ngô và hoa màu đạt năng suất ổn định lâu dài.

  5. Lượng mưa mùa mưa tại Quản Bạ tác động thế nào đến quá trình bóc tách đất mặt? Hơn 70% đến 85% lượng mưa cả năm dồn dập vào 6 tháng mùa mưa với nhiều đợt mưa rào cường độ vượt 25 mm/h. Các hạt mưa kích thước lớn (đường kính lên tới 5 mm) khi rơi tự do đạt vận tốc chạm đất 9 m/s, tạo lực va đập mạnh làm vỡ vụn các cấu trúc đoàn lạp đất mặt, tạo điều kiện cho dòng chảy tràn cuốn trôi lớp đất màu mỡ.

Kết luận

  • Luận văn đã xác lập thành công bộ thông số xói mòn đất chuyên biệt cho huyện Quản Bạ theo mô hình USLE, lượng hóa chi tiết các hệ số R, K, LS và C*P trên nền tảng dữ liệu thực địa và GIS.
  • Xây dựng hoàn chỉnh bản đồ nguy cơ xói mòn đất tỷ lệ 1:25.000, phân định rõ ràng 7 cấp độ xói mòn và xác định chính xác các điểm nóng suy thoái đất cần ưu tiên xử lý.
  • Đề xuất hệ thống giải pháp công trình và phi công trình đồng bộ, trọng tâm là mô hình canh tác nông lâm kết hợp SALT, SALT2 và khoanh nuôi phục hồi rừng phòng hộ sườn dốc trên 25 độ.
  • Lộ trình triển khai khuyến nghị chia thành hai giai đoạn: giai đoạn 2016 – 2020 tập trung hoàn thiện chính sách giao đất giao rừng và thí điểm mô hình SALT; giai đoạn 2020 – 2025 nhân rộng trên toàn bộ 13 xã, thị trấn.
  • Chính quyền địa phương, các nhà khoa học và cộng đồng dân cư cần nhanh chóng phối hợp đưa cơ sở dữ liệu bản đồ này vào thực tế quy hoạch sản xuất nông nghiệp, nhằm bảo vệ tài nguyên đất đai quý báu và đảm bảo an ninh sinh kế bền vững cho vùng cao Quản Bạ.