Tổng quan nghiên cứu

Thoái hóa đất do xói mòn đang là thách thức sinh thái nghiêm trọng hàng đầu tại các vùng đồi núi Việt Nam. Theo thống kê thổ nhưỡng quốc gia, cả nước hiện có khoảng 25 triệu hecta đất dốc với cường độ mất đất bình quân ước tính lên tới 10 tấn/ha/năm, trong đó 10% đến 20% tổng diện tích lãnh thổ chịu ảnh hưởng xói mòn từ mức độ trung bình đến rất mạnh. Tại nhiều khu vực thuộc vùng trung du và miền núi phía Bắc, tỷ lệ mất đất mặt ghi nhận thực tế từ 1% đến 2% mỗi năm, làm suy giảm nghiêm trọng độ phì nhiêu và đe dọa trực tiếp đến an ninh lương thực.

Huyện Văn Bàn, tỉnh Lào Cai sở hữu tổng diện tích tự nhiên 142.345,46 ha, mang đặc thù địa hình lòng chảo bị chia cắt mạnh với nhiều sườn dốc hiểm trở vượt quá 25 độ phân bố dọc theo lưu vực suối Chăn và suối Nhù. Trong bối cảnh biến đổi khí hậu diễn biến phức tạp, lượng mưa tại địa phương tập trung tới 70% vào mùa mưa dồn dập, kết hợp với việc suy giảm thảm phủ thực vật và canh tác nương rẫy chưa hợp lý đã gây ra hiện tượng rửa trôi bề mặt khốc liệt, lũ quét và sạt lở đất thường xuyên.

Luận văn thạc sĩ tập trung giải quyết bài toán định lượng nguy cơ mất đất thông qua việc xác định chính xác các hệ số xói mòn đơn tính, xây dựng bản đồ phân cấp nguy cơ xói mòn đất tỷ lệ 1:150.000 và đề xuất các mô hình sản xuất nông - lâm nghiệp bền vững cho huyện Văn Bàn. Kết quả nghiên cứu có ý nghĩa thực tiễn then chốt trong việc bảo vệ 89.436,32 ha đất quy hoạch lâm nghiệp (chiếm 62,83% diện tích tự nhiên), duy trì sản lượng lương thực có hạt đạt trên 58.637 tấn/năm, nâng cao giá trị canh tác đạt 65 triệu đồng/ha và hạn chế tối đa các tổn thất kinh tế do thiên tai từng gây thiệt hại hơn 18,3 tỷ đồng cho địa phương.

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 trên nền tảng lý thuyết động lực xói mòn thủy lực do Ellison khởi xướng năm 1944, phân tích hiện tượng xói mòn như một hàm số phụ thuộc vào 5 biến số: loại đất, độ dốc địa hình, mật độ che phủ của thảm thực vật, lượng mưa và cường độ mưa. Quá trình xói mòn được chia tách thành 3 pha cơ học liên hoàn: xung lực hạt mưa phá vỡ cấu trúc bề mặt đất, dòng chảy vận chuyển các hạt đất bở rời theo sườn dốc và quá trình sa lắng trầm tích ở vùng trũng hạ lưu.

Trọng tâm lý thuyết của đề tài là việc ứng dụng Phương trình mất đất phổ dụng USLE của Wischmeier và Smith (1978) kết hợp mô hình chỉ số xói mòn SEIM và mô hình vật lý WEPP để lượng hóa không gian. Luận văn chuẩn hóa 5 khái niệm chuyên ngành then chốt tương ứng với các biến số trong công thức tính lượng mất đất hàng năm A = R x K x L x S x C x P:

  1. Hệ số xói mòn do mưa (R): Thể hiện động năng và cường độ của các trận mưa tạo dòng chảy mặt.
  2. Hệ số tính dễ xói mòn của đất (K): Đại lượng nghịch đảo phản ánh độ bền cấu trúc và tính kháng xói mòn của từng loại thổ nhưỡng.
  3. Hệ số độ dốc và chiều dài sườn dốc (LS): Thể hiện tác động khuếch đại của địa hình đến vận tốc và lưu lượng dòng chảy.
  4. Hệ số quản lý lớp phủ thực bì (C): Chỉ số thể hiện năng lực cản hạt mưa và giữ đất của các tầng tán rừng và thảm cỏ.
  5. Hệ số biện pháp bảo vệ đất (P): Tỷ số phản ánh hiệu quả giảm dòng chảy của các kỹ thuật canh tác như làm ruộng bậc thang, tạo rãnh thoát nước hoặc trồng cây theo đường đồng mức.

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

Nguồn dữ liệu của luận văn được tổng hợp đa tầng từ chuỗi quan trắc khí tượng của 7 trạm đại biểu khu vực Tây Bắc (Sa Pa, Than Uyên, Mù Căng Chải, Minh Lương, Ngòi Nhù, Văn Bàn, Tả Van) với lượng mưa trung bình dao động từ 1.302 mm đến 2.844 mm/năm. Nguồn tài liệu bản đồ gồm bản đồ hiện trạng sử dụng đất năm 2017, bản đồ thổ nhưỡng tỷ lệ lớn và hồ sơ rà soát quy hoạch ba loại rừng quy mô 89.436,32 ha do Sở Nông nghiệp và Phát triển nông thôn tỉnh Lào Cai cung cấp.

Luận văn áp dụng phương pháp chọn mẫu khảo sát thực địa có chủ đích tại 14 khu vực sườn dốc xung yếu có độ dốc trên 25 độ và dải đất canh tác ven các lưu vực sông suối lớn. Cỡ mẫu phỏng vấn chuyên gia bao gồm 20 cán bộ kỹ thuật và chuyên viên quản lý thuộc Phòng Tài nguyên và Môi trường, Hạt Kiểm lâm và Phòng Nông nghiệp và Phát triển nông thôn huyện Văn Bàn nhằm chuẩn hóa các hệ số canh tác P và hệ số thực bì C theo điều kiện bản địa.

Phương pháp phân tích cốt lõi là công nghệ Hệ thống Thông tin Địa lý (GIS) trên nền phần mềm ArcGIS 10.4, sử dụng thuật toán phân tích không gian Spatial Analyst để chồng xếp các lớp dữ liệu Raster đơn tính (R, K, LS, C, P) với độ phân giải lưới ô vuông chuẩn hóa. Thuật toán nội suy Kriging và IDW được lựa chọn để mô hình hóa trường mưa liên tục từ dữ liệu điểm trạm. Việc lựa chọn phương pháp GIS kết hợp USLE giúp xử lý nhanh khối lượng dữ liệu khổng lồ trên quy mô 142.345,46 ha với độ tin cậy khoa học cao, tiết kiệm tối đa chi phí đo đạc thực địa định kỳ.

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 không gian và phân tích dữ liệu thực nghiệm đã chỉ ra 4 phát hiện khoa học mang tính quyết định về hiện trạng xói mòn tại huyện Văn Bàn:

Thứ nhất, hệ số tính dễ xói mòn của đất (K) trên toàn địa bàn dao động trong khoảng từ 0,10 đến 0,44. Trong đó, dải giá trị từ 0,23 đến 0,28 chiếm tỷ trọng lớn nhất với 58,33% tổng diện tích tự nhiên, tập trung chủ yếu ở nhóm đất đỏ vàng chiếm 58.151,0 ha (40,9% diện tích huyện) và đất mùn vàng đỏ chiếm 44.215,0 ha (31,1% diện tích huyện). Điều này chứng minh phần lớn quỹ đất vùng nghiên cứu có kết cấu cơ giới kém bền vững và dễ bị tổn thương khi mất lớp phủ.

Thứ hai, chế độ mưa tạo dòng chảy mặt diễn biến cực đoan. Với tổng lượng mưa bình quân 1.500 mm/năm nhưng có tới 70% lượng mưa trút xuống dồn dập trong 4 tháng (từ tháng 7 đến tháng 10). Hệ số xói mòn mưa R theo công thức thực nghiệm của Nguyễn Trọng Hà đạt mức biến thiên từ 650 đến trên 1.500 J/m2, hình thành các đợt lũ quét với lưu lượng dòng chảy tức thời đạt tới 1.000 m3/s tại các nhánh suối nhỏ.

Thứ ba, sự khác biệt về cấu trúc thảm thực bì tạo ra sự phân hóa sâu sắc của hệ số che phủ C. Rừng tự nhiên 3 tầng có độ tàn che 0,7 - 0,8 đạt hệ số bảo vệ tối ưu với C = 0,0070, trong khi đất trồng trọt nương rẫy và rừng thoái hóa có hệ số C tăng vọt lên 0,0186 đến 0,0194 (cao gấp gần 3 lần). Số liệu thực nghiệm cho thấy khi độ tàn che của rừng suy giảm từ mức 0,8 xuống 0,3 - 0,4 thì lượng xói mòn đất tăng thêm 42,2% và dòng chảy mặt tăng 30,4%.

Thứ tư, bản đồ phân cấp nguy cơ xói mòn tỷ lệ 1:150.000 xác định diện tích đất chưa sử dụng lên tới 31.814,33 ha (chiếm 22,35% diện tích tự nhiên) nằm trên các sườn núi dốc trên 25 độ là vùng chịu cấp độ xói mòn mạnh đến rất mạnh, có nguy cơ rửa trôi mất trắng lớp thổ nhưỡng mặt.

Thảo luận kết quả

Nguyên nhân cốt lõi dẫn đến tình trạng xói mòn cục bộ nghiêm trọng tại Văn Bàn xuất phát từ sự tương tác bất lợi giữa địa hình lòng chảo dốc đứng và các hoạt động kinh tế - xã hội. Việc chuyển đổi đất rừng sang đất canh tác cây hàng năm ngắn ngày trên sườn dốc cùng với 1.091,74 ha đất sử dụng cho hoạt động khoáng sản đã phá vỡ thảm thực vật tự nhiên, triệt tiêu khả năng điều tiết nước của lưu vực.

Khi so sánh với các nghiên cứu kinh điển về đất dốc tại vùng núi phía Bắc, kết quả tính toán hoàn toàn tương thích với quy luật: độ dốc địa hình tăng 2 lần làm cường độ xói mòn tăng từ 2 đến 4 lần; chiều dài sườn dốc tăng 2 lần khiến lượng đất mất tăng từ 2 đến 7,5 lần.

Dữ liệu nghiên cứu được biểu diễn trực quan qua bảng ma trận chồng lớp Raster và bản đồ phân vùng nguy cơ, cho thấy các điểm nóng xói mòn tập trung dày đặc tại các xã vùng cao như Nậm Xây, Nậm Xé, Nậm Tha. Sự suy giảm độ phì nhiêu đất không chỉ làm giảm năng suất cây trồng mà còn bồi lắng một lượng trầm tích khổng lồ xuống hệ thống lòng hồ của 12 công trình thủy điện đang vận hành trên địa bàn (tổng công suất 104,4 MW), làm giảm tuổi thọ công trình và đe dọa sinh kế hạ du.

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

Dựa trên kết quả đánh giá thực trạng, luận văn đề xuất 4 nhóm giải pháp kỹ thuật và quản lý đồng bộ nhằm giảm thiểu xói mòn đất:

  1. Thiết lập mô hình canh tác nông lâm kết hợp theo đường đồng mức: Chuyển đổi 6.609,83 ha đất trồng cây hàng năm kém hiệu quả trên đất dốc sang mô hình nông lâm kết hợp đa tầng tán (trồng xen chè lâu năm, thảo quả, cây ăn quả với băng cây họ đậu phân xanh). Mục tiêu hạ chỉ số canh tác P xuống mức 0,5 - 0,6 và giảm trên 50% lượng đất trôi rửa bề mặt. Đơn vị chủ trì: Phòng Nông nghiệp và Phát triển nông thôn huyện Văn Bàn, lộ trình thực hiện từ năm 2025 đến 2028.

  2. Phục hồi thảm phủ rừng và đóng cửa rừng tự nhiên: Tăng cường khoán bảo vệ 46.597 ha rừng và đẩy mạnh trồng rừng mới trên diện tích 31.814,33 ha đất chưa sử dụng, ưu tiên trồng rừng gỗ lớn hỗn giao bản địa có từ 2 đến 3 tầng tán để đạt hệ số C dưới 0,008. Mục tiêu nâng tỷ lệ che phủ rừng toàn huyện từ 63,6% lên 66,5% vào năm 2030. Chủ thể thực hiện: Hạt Kiểm lâm huyện Văn Bàn và các Ban Quản lý rừng phòng hộ.

  3. Xây dựng hệ thống công trình sinh học - cơ học chống xói mòn: Đầu tư xây dựng bờ kè đá xếp, bậc thềm ruộng bậc thang kiên cố và hệ thống mương gom nước phân tán năng lượng dòng chảy mặt tại các sườn dốc trên 15 độ thuộc lưu vực suối Chăn và suối Nhù. Target metric: Giảm thiểu 80% nguy cơ sạt lở đất ven suối và giảm thiểu thiệt hại thiên tai (hạn chế mức rủi ro 18,3 tỷ đồng/năm). Chủ thể phối hợp: Phòng Tài nguyên và Môi trường cùng UBND 23 xã, thị trấn, thực hiện trong giai đoạn 2025 - 2027.

  4. Hoàn thiện quy hoạch sử dụng đất gắn liền với kiểm soát khoáng sản: Rà soát nghiêm ngặt 1.091,74 ha đất hoạt động khoáng sản, bắt buộc các doanh nghiệp khai thác vàng sa khoáng, sắt, fenspat phải thực hiện đề án hoàn nguyên môi trường, trồng cây phủ xanh bãi thải ngay sau khi kết thúc khai thác nhằm triệt tiêu các mương xói nhân tạo. Timeline kiểm tra định kỳ hàng quý bắt đầu từ năm 2025.

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

Công trình nghiên cứu mang giá trị học thuật và ứng dụng thực tiễn cao, đặc biệt hữu ích cho 4 nhóm đối tượng:

  1. Cơ quan quản lý Nhà nước về tài nguyên và quy hoạch nông thôn: Lãnh đạo và chuyên viên Sở Tài nguyên và Môi trường, Sở Nông nghiệp và Phát triển nông thôn tỉnh Lào Cai, UBND huyện Văn Bàn có thể sử dụng trực tiếp bản đồ phân cấp nguy cơ xói mòn 1:150.000 để tích hợp vào quy hoạch sử dụng đất cấp huyện giai đoạn 2021 - 2030 và xây dựng phương án phòng chống thiên tai cấp cơ sở.

  2. Nhà khoa học, giảng viên và học viên cao học: Các nhà nghiên cứu chuyên ngành Khoa học Môi trường, Quản lý Đất đai, Địa lý học có thể tham khảo quy trình tích hợp mô hình USLE trên nền tảng ArcGIS 10.4 cùng bộ hệ số thực nghiệm (R, K, LS, C, P) đã được chuẩn hóa theo điều kiện miền núi nhiệt đới gió mùa Việt Nam.

  3. Doanh nghiệp năng lượng và khai thác khoáng sản: Chủ đầu tư 12 dự án nhà máy thủy điện (công suất 104,4 MW) và các doanh nghiệp mỏ trên địa bàn có thể khai thác số liệu dòng chảy, lượng bùn cát bồi lắng để tối ưu hóa thiết kế hồ chứa, lập báo cáo đánh giá tác động môi trường và triển khai các công trình hoàn nguyên đất dốc.

  4. Các tổ chức phát triển nông thôn và hợp tác xã nông nghiệp: Các kỹ sư khuyến nông và chủ trang trại vùng cao có thể ứng dụng trực tiếp các hướng dẫn kỹ thuật canh tác ruộng bậc thang, trồng cây che phủ theo đường đồng mức nhằm bảo tồn độ phì cho đất sản xuất nông - lâm nghiệp.

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

  1. Phương pháp tích hợp mô hình USLE và công nghệ GIS mang lại ưu điểm vượt trội gì trong đánh giá xói mòn? Phương pháp này cho phép lượng hóa nguy cơ mất đất trên quy mô không gian rộng lớn 142.345,46 ha trong thời gian ngắn với chi phí thấp. Việc số hóa các lớp dữ liệu R, K, LS, C, P dạng Raster giúp phân vùng xói mòn chi tiết đến từng ô lưới tọa độ, khắc phục hoàn toàn tính định tính của các phương pháp phân vùng truyền thống.

  2. Yếu tố tự nhiên nào chi phối mạnh nhất đến cường độ xói mòn đất tại huyện Văn Bàn? Lượng mưa và địa hình là hai yếu tố chi phối mạnh nhất. Với lượng mưa trung bình 1.500 mm/năm nhưng tập trung tới 70% vào mùa lũ, kết hợp địa hình lòng chảo có độ dốc phổ biến trên 25 độ đã tạo ra động năng dòng chảy bề mặt cực lớn, bóc tách nhanh chóng lớp đất mặt tại các sườn đồi dốc.

  3. Hệ số tính dễ xói mòn K của các nhóm đất chính tại Văn Bàn phân bố ra sao? Hệ số K tại Văn Bàn dao động từ 0,10 đến 0,44. Trong đó, khoảng giá trị từ 0,23 đến 0,28 chiếm tỷ trọng áp đảo với 58,33% diện tích tự nhiên, tập trung ở nhóm đất đỏ vàng Fs, Fa (chiếm 40,9%) và đất mùn vàng đỏ HF (chiếm 31,1%), phản ánh đặc tính đất dễ bị rửa trôi khi mất lớp thảm phủ.

  4. Mức độ suy giảm thảm thực bì ảnh hưởng định lượng như thế nào đến lượng đất bị xói mòn? Khi độ tàn che của rừng giảm từ 0,7 - 0,8 xuống mức 0,3 - 0,4, lượng đất bị xói mòn tăng lên 42,2% và lưu lượng dòng chảy mặt tăng 30,4%. Rừng tự nhiên đa tầng tán có hệ số C chỉ 0,0070, trong khi đất trọc hoặc nương rẫy có hệ số C lên tới 0,0194, khiến lượng đất mất tăng gấp gần 3 lần.

  5. Giải pháp canh tác nào được đánh giá là tối ưu nhất cho vùng đất dốc Văn Bàn? Giải pháp hiệu quả nhất là mô hình nông lâm kết hợp đa tầng tán theo đường đồng mức kết hợp xếp đá ngăn dòng. Kỹ thuật này giúp giảm hệ số canh tác P từ 0,9 xuống 0,5 - 0,6, hạn chế hơn 50% lượng đất trôi rửa mà vẫn đảm bảo thu nhập bình quân trên 65 triệu đồng/ha/năm cho người dân.

Kết luận

  • Luận văn đã hoàn thành xuất sắc việc đánh giá toàn diện hiện trạng xói mòn đất trên quy mô 142.345,46 ha của huyện Văn Bàn bằng công nghệ viễn thám và GIS hiện đại.
  • Chuẩn hóa thành công bộ 5 tham số thực nghiệm theo phương trình mất đất USLE phù hợp với điều kiện tự nhiên đặc thù của vùng núi Đông Bắc.
  • Xây dựng bản đồ phân cấp nguy cơ xói mòn đất tỷ lệ 1:150.000, nhận diện chính xác 31.814,33 ha đất chưa sử dụng trên sườn dốc cao là khu vực xung yếu cần bảo vệ khẩn cấp.
  • Đề xuất khung 4 giải pháp công trình và phi công trình có tính khả thi cao, hướng tới mục tiêu duy trì tỷ lệ che phủ rừng trên 65% và giảm thiểu 80% thiệt hại sạt lở đất.
  • Cung cấp cơ sở dữ liệu số hóa tin cậy phục vụ đắc lực cho công tác quy hoạch sử dụng đất và chiến lược phát triển nông - lâm nghiệp bền vững tại địa phương.

Về lộ trình tiếp theo, chính quyền địa phương và các cơ quan chuyên ngành cần sớm triển khai thử nghiệm các mô hình canh tác nông lâm kết hợp mẫu tại các xã Nậm Tha và Nậm Xây trong giai đoạn 2025 - 2026 trước khi nhân rộng toàn huyện. Kính mời các nhà quản lý, nhà nghiên cứu và bạn đọc quan tâm tham khảo toàn văn luận văn thạc sĩ để khai thác chi tiết các bản đồ chuyên đề và giải pháp kỹ thuật chuyên sâu.