Tổng quan nghiên cứu

Lưu vực Sơn Diệm thuộc huyện Hương Sơn, tỉnh Hà Tĩnh giữ vai trò là lưu vực đầu nguồn đặc biệt xung yếu của hệ thống sông Ngàn Phố và sông Lam. Khu vực nghiên cứu sở hữu tổng diện tích tự nhiên là 78.938,5 ha với chu vi đường phân thủy đạt 154,6 km, chiều dài lưu vực 30 km và hệ số hình dạng lưu vực là 1,54. Địa hình khu vực mang tính phân cắt mạnh mẽ, trong đó đồi núi cao chiếm tới 75% tổng diện tích, đồi núi thấp chiếm 15% và 10% còn lại là các thung lũng ven sông suối. Độ cao địa hình dao động từ 11 m đến 1.911 m với độ dốc trung bình 16,3 độ, cá biệt có nhiều sườn dốc hiểm trở lên đến 61 độ.

Bên cạnh địa hình phức tạp, lưu vực chịu ảnh hưởng của chế độ khí hậu nhiệt đới gió mùa khắc nghiệt với lượng mưa bình quân năm dao động từ 2.100 mm đến 2.500 mm phân bổ trong 150 đến 160 ngày mưa, ghi nhận lượng mưa ngày cực đoan từng đạt mốc 657,2 mm. Tình trạng khai thác gỗ không bền vững và chuyển đổi mục đích sử dụng đất trong các giai đoạn trước đã làm suy giảm đáng kể diện tích rừng tự nhiên, dẫn đến nguy cơ xói mòn rửa trôi đất mặt nghiêm trọng, gây lũ quét, sạt lở bờ sông và bồi lắng các công trình thủy lợi hạ lưu.

Mục tiêu cốt lõi của nghiên cứu là lượng hóa cường độ xói mòn tiềm năng cho 8 trạng thái rừng và đất lâm nghiệp chính tại lưu vực Sơn Diệm, phân tích mối tương quan giữa các yếu tố cấu trúc thảm thực vật với khả năng giữ đất, từ đó thành lập bản đồ xói mòn tiềm năng theo tiêu chuẩn kỹ thuật số và xây dựng hệ thống giải pháp lâm sinh đồng bộ nhằm quản lý bền vững tài nguyên đất dố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 tổng hợp lý thuyết động lực học xói mòn đất nhiệt đới và các mô hình dự báo mất đất hiện đại:

  • Phương trình mất đất phổ dụng USLE/RUSLE: Kế thừa nguyên lý cơ bản của Wischmeier và Smith về sự tương tác giữa động năng mưa, tính dễ xói mòn của đất, chiều dài - độ dốc sườn đất, lớp phủ thực vật và biện pháp canh tác.
  • Mô hình xói mòn đất lâm nghiệp nhiệt đới của Vương Văn Quỳnh: Mô hình toán học thực nghiệm được chuẩn hóa cho điều kiện địa hình dốc và rừng Việt Nam theo công thức: $$d = \frac{2 \cdot K \cdot \sin^2\alpha}{\left[\left(\frac{TC}{H}\right) + CP + TM\right]^2 \cdot X}$$ Trong đó: $d$ là cường độ xói mòn đất (mm/năm); $\alpha$ là độ dốc bề mặt (độ); $TC$ là độ tàn che tầng cây cao; $H$ là chiều cao bình quân tầng cây cao (m); $CP$ là độ che phủ của cây bụi thảm tươi; $TM$ là tỷ lệ che phủ của lớp thảm khô, thảm mục; $X$ là độ xốp tổng số tầng đất mặt 0-5 cm; $K$ là chỉ số năng lực gây xói mòn của lượng mưa hàng tháng.
  • Tiêu chuẩn Quốc gia TCVN 5299:2009: Tiêu chuẩn phân cấp độ xói mòn đất phục vụ đánh giá mức độ suy thoái đất lâm nghiệp và nông nghiệp.
  • Hệ thống phân loại trạng thái rừng: Thực hiện theo Thông tư số 34/2009/TT-BNNPTNT của Bộ Nông nghiệp và Phát triển nông thôn, phân định rõ các trạng thái rừng tự nhiên lá rộng thường xanh (giàu, trung bình, nghèo, phục hồi), rừng hỗn giao gỗ tre nứa, rừng trồng các cấp tuổi và đất trống thoái hóa.

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

Phương pháp nghiên cứu kết hợp giữa điều tra thực địa theo ô tiêu chuẩn, phân tích lý hóa mẫu đất trong phòng thí nghiệm và mô hình hóa không gian bằng hệ thống thông tin địa lý GIS:

  • Cỡ mẫu và phương pháp chọn mẫu: Bố trí 80 ô tiêu chuẩn chính đại diện cho 8 trạng thái rừng theo phương pháp rút mẫu điển hình có hệ thống. Kích thước mỗi ô tiêu chuẩn là 1.000 m² (40m x 25m) đối với rừng tự nhiên và 500 m² (20m x 25m) đối với rừng trồng. Tiến hành khảo sát bổ sung tại 150 điểm thực địa nhằm kiểm chứng các cấp tuổi của rừng trồng keo, thông, cao su và hiện trạng đất trống.
  • Đo đếm cấu trúc tầng tán: Chiều cao vút ngọn được đo bằng máy đo cao kỹ thuật số Vertex IV. Độ tàn che ($TC$), độ che phủ thảm tươi ($CP$) và thảm mục ($TM$) được xác định bằng phương pháp ống ngắm thị sai ngẫu nhiên qua 80 điểm quan trắc phân bố đều trên mạng lưới ô tiêu chuẩn.
  • Lấy mẫu và phân tích tính chất đất: Trong mỗi ô tiêu chuẩn, thiết lập 5 ô dạng bản kích thước 1 m² (tại 4 góc và tâm ô) để lấy mẫu đất nguyên dạng tầng 0-5 cm bằng ống dung trọng thể tích chuẩn 100 cm³. Tổng hợp 400 mẫu đất đơn thành các mẫu hỗn hợp để sấy khô kiệt ở nhiệt độ 105°C trong 6 giờ, xác định dung trọng ($D$), tỷ trọng ($d$) qua bình picnomet và tính toán độ xốp tổng số ($X$).
  • Xử lý số liệu không gian và phân tích GIS: Xây dựng mô hình số độ cao DEM độ phân giải lưới 30m x 30m từ bản đồ địa hình tỷ lệ 1/25.000 và 1/50.000. Ứng dụng thuật toán Hydrology trong phần mềm ArcGIS kết hợp MapInfo để khoanh vẽ đường phân thủy lưu vực tự động dựa trên tọa độ điểm đầu ra tại ngã ba sông Con và sông Ngàn Phố. Chồng lớp raster đa nhân tố theo công thức giải tích để thành lập bản đồ xói mòn tiềm năng hoàn chỉnh.

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

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

  • Cấu trúc thảm thực vật khác biệt rõ rệt giữa các trạng thái: Rừng tự nhiên lá rộng thường xanh giàu và trung bình duy trì độ tàn che tầng cây cao trung bình từ 0,72 đến 0,85, chiều cao bình quân từ 16,5 m đến 22,4 m, độ che phủ thảm mục tích lũy đạt từ 65% đến 78% và độ xốp tầng đất mặt duy trì ở mức cao từ 56,4% đến 62,8%.
  • Hiệu quả khống chế xói mòn của rừng tự nhiên: Lượng xói mòn tiềm năng ở các trạng thái rừng tự nhiên giàu và trung bình chỉ dao động từ 0,25 đến 0,65 mm/năm (tương đương khoảng 3,1 đến 8,2 tấn/ha/năm), thấp hơn nhiều so với ngưỡng xói mòn cho phép của đất nhiệt đới là 0,8 mm/năm theo tiêu chuẩn sinh thái của Hudson.
  • Nguy cơ xói mòn cực đại tại đất trống và rừng non: Trạng thái đất trống có cây bụi và đất trống thoái hóa có lượng mất đất tiềm năng lên tới 11,8 đến 18,4 mm/năm (tương đương 145 đến 230 tấn/ha/năm), cao gấp 25 đến 35 lần so với diện tích có rừng tự nhiên che phủ. Rừng trồng thuần loài keo từ 1 đến 3 năm tuổi có cường độ xói mòn trung bình đạt 4,5 đến 7,2 mm/năm do tỷ lệ che phủ thảm mục chỉ đạt dưới 25% và cấu trúc tán đơn tầng.
  • Phân cấp xói mòn theo TCVN 5299:2009: Toàn bộ lưu vực Sơn Diệm được phân định thành 5 cấp xói mòn rõ rệt. Diện tích đất có nguy cơ xói mòn rất mạnh và đặc biệt mạnh (mất đất trên 50 tấn/ha/năm) chiếm khoảng 18,6% diện tích, tập trung chủ yếu tại các sườn dốc canh tác nông nghiệp nương rẫy, các khu vực đồi trọc thuộc xã Sơn Hồng, Sơn Lĩnh và dải đất giáp ranh đường Quốc lộ 8A.

Thảo luận kết quả

Dữ liệu nghiên cứu chứng minh mối quan hệ phi tuyến tính chặt chẽ giữa các chỉ tiêu sinh thái rừng và cường độ mất đất. Khi độ tàn che tầng cây cao ($TC$) kết hợp cùng độ che phủ thảm tươi ($CP$) và thảm mục ($TM$) vượt mốc tổng hòa 1,2 đơn vị quy ước, động năng rơi tự do của hạt mưa nhiệt đới bị triệt tiêu từ 80% đến 92%, ngăn chặn hoàn toàn hiện tượng bóc tách các hạt sét và cát mịn trên bề mặt.

Khi trình bày các kết quả qua biểu đồ phân tán đa biến và ma trận không gian raster, mối quan hệ giữa độ dốc sườn núi ($\alpha$) và lượng mất đất thể hiện rõ dạng hàm mũ. Trên các sườn dốc trên 25 độ, nếu lớp thảm mục mặt đất bị mất đi do đốt thực bì hoặc dọn sạch nương rẫy, cường độ dòng chảy mặt tăng vọt gấp 4,2 lần, kéo theo tốc độ cuốn trôi dinh dưỡng tầng mặt diễn ra tức thì. So sánh với các công trình nghiên cứu của Võ Đại Hải tại vùng Tây Bắc và Phạm Văn Điển tại lưu vực sông Đà, kết quả tại lưu vực Sơn Diệm hoàn toàn tương thích về quy luật suy giảm xói mòn theo chiều tăng của độ xốp đất và độ dày tầng thảm mục. Việc ứng dụng công nghệ GIS raster mapping đã khắc phục triệt để các sai số cục bộ của phương pháp chuyên gia truyền thống, cung cấp bức tranh định lượng trực quan phục vụ công tác quy hoạch lâm nghiệp đầu nguồn.

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

  • Thiết lập mô hình lâm sinh hỗn giao bản địa: Chuyển đổi phương thức trồng rừng keo thuần loài sang mô hình hỗn giao đa tầng giữa keo lai với các loài cây bản địa như lim xanh, mỡ, lát hoa và trám trắng; kiểm soát mật độ trồng ban đầu đạt từ 1.600 đến 2.000 cây/ha nhằm nâng độ tàn che lên trên 0,70 chỉ sau 3 năm, giảm thiểu 60% lượng xói mòn trong giai đoạn kiến thiết cơ bản.
  • Khoanh nuôi phục hồi và xúc tiến tái sinh: Thực hiện khoanh nuôi bảo vệ nghiêm ngặt trên 12.500 ha rừng nghèo kiệt và đất trống có cây tái sinh tại các xã thượng lưu như Sơn Kim 1, Sơn Kim 2 và Sơn Hồng; nghiêm cấm các hoạt động chăn thả gia súc tự do và chặt phá thảm tươi; mục tiêu đưa cường độ xói mòn về mức an toàn dưới 5 tấn/ha/năm trong giai đoạn 2025-2030 do Chi cục Kiểm lâm Hà Tĩnh và Ban quản lý Rừng phòng hộ Hương Sơn chỉ đạo thực hiện.
  • Áp dụng kỹ thuật canh tác bảo tồn đất dốc: Bắt buộc áp dụng biện pháp trồng cây theo đường đồng mức, xây dựng các băng chắn xanh bằng cỏ vetiver hoặc cây họ đậu trên toàn bộ các sườn dốc nông lâm kết hợp có độ dốc trên 15 độ; tuyệt đối cấm đốt dọn thực bì toàn diện khi khai thác rừng trồng thương ngạch nhằm giữ lại ít nhất 80% lượng thảm mục che phủ bề mặt.
  • Ứng dụng hệ thống giám sát cảnh báo xói mòn số: Xây dựng cơ sở dữ liệu GIS viễn thám tích hợp mô hình xói mòn tự động tại UBND huyện Hương Sơn; thực hiện cập nhật biến động thảm phủ thực vật định kỳ 12 tháng/lần qua ảnh vệ tinh độ phân giải cao để khoanh vùng chính xác các điểm nóng nguy cơ sạt lở trước mùa mưa bão hàng năm.

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

  • Học viên cao học, nghiên cứu sinh và nhà khoa học chuyên ngành Lâm học, Quản lý tài nguyên rừng, Thủy văn môi trường: Nguồn tài liệu học thuật giá trị về phương pháp ứng dụng phương trình xói mòn thực nghiệm Vương Văn Quỳnh kết hợp công cụ GIS raster để tính toán động lực dòng chảy đất dốc.
  • Lãnh đạo và cán bộ kỹ thuật tại Chi cục Kiểm lâm, Sở Nông nghiệp và PTNT tỉnh Hà Tĩnh: Cơ sở dữ liệu khoa học trực tiếp để xây dựng quy hoạch ba loại rừng, phân vùng xung yếu phòng hộ đầu nguồn sông Ngàn Phố và lập đề án chi trả dịch vụ môi trường rừng.
  • Ban quản lý Rừng phòng hộ, các Công ty Lâm nghiệp tại khu vực Bắc Trung Bộ: Tài liệu hướng dẫn thực hành trong việc lựa chọn cơ cấu cây trồng phòng hộ, xác định chu kỳ khai thác và thiết kế các đai rừng giữ đất chắn lũ hiệu quả.
  • Các tổ chức phi chính phủ và cơ quan quản lý dự án bảo tồn lưu vực sông: Cung cấp khung phân tích chuẩn mực về lượng hóa dịch vụ hệ sinh thái đất và nước phục vụ thẩm định các dự án phục hồi cảnh quan rừng bền vững.

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

  • Mô hình xói mòn đất của Vương Văn Quỳnh có ưu điểm gì nổi bật so với mô hình USLE quốc tế khi áp dụng tại Việt Nam? Mô hình của Vương Văn Quỳnh tích hợp trực tiếp các thông số cấu trúc lâm học nhiệt đới nhiều tầng tán như tỷ số tàn che trên chiều cao cây ($TC/H$), độ che phủ thảm tươi ($CP$) và độ xốp tầng đất mặt ($X$). Điều này giúp mô hình phản ánh sát thực tế cấu trúc rừng phức tạp tại Việt Nam thay vì chỉ dựa vào hệ số che phủ canh tác chung chung như phương trình USLE chuẩn.

  • Tại sao chiều cao tầng cây cao ($H$) lại tỷ lệ thuận với cường độ xói mòn trong công thức tính toán? Khi tầng cây có chiều cao vút ngọn càng lớn nhưng thiếu cấu trúc tầng tán trung gian và thảm mục phía dưới, các giọt nước mưa đọng trên tán lá sẽ tích tụ thành giọt có kích thước lớn hơn. Khi rơi tự do từ độ cao trên 15 m xuống đất trống, động năng va đập của giọt nước tăng mạnh, làm phá vỡ kết cấu viên xốp của đất mặt nhanh hơn.

  • Phương pháp khoanh vẽ lưu vực kết hợp giữa mô hình DEM và bản đồ địa hình đem lại lợi ích gì? Phương pháp tự động bằng DEM 30m x 30m giúp xác định chính xác mạng lưới tích lũy dòng chảy và diện tích lưu vực khách quan. Việc đối soát thủ công với bản đồ địa hình 1/25.000 giúp hiệu chỉnh các đường phân thủy đi qua địa hình yên ngựa hoặc đồi bát úp, hạn chế triệt để sai số biên lưu vực ở mức dưới 2%.

  • Tại sao rừng trồng keo thuần loài ở giai đoạn 1-3 năm tuổi lại có nguy cơ xói mòn cao? Trong 3 năm đầu, tán cây keo chưa giao nhau khiến độ tàn che chỉ đạt dưới 0,35. Đồng thời, tập quán dọn sạch thực bì và đốt nương của người dân làm mất đi lớp thảm khô che phủ mặt đất, khiến đất mặt bị phơi trần hoàn toàn dưới các trận mưa rào cường độ lớn trên 50 mm/giờ.

  • Làm thế nào để ứng dụng bản đồ xói mòn tiềm năng vào công tác giao đất giao rừng tại địa phương? Bản đồ phân cấp xói mòn cung cấp thông tin chi tiết đến từng lô rừng. Những diện tích nằm trong cấp xói mòn mạnh và rất mạnh (độ dốc trên 25 độ, đất mỏng) bắt buộc phải đưa vào diện rừng phòng hộ xung yếu để khoanh nuôi bảo vệ, tuyệt đối không quy hoạch làm rừng sản xuất khai thác trắng.

Kết luận

  • Lượng hóa toàn diện cường độ xói mòn tiềm năng trên toàn bộ 78.938,5 ha diện tích lưu vực Sơn Diệm bằng sự kết hợp chuẩn xác giữa mô hình thực nghiệm và công nghệ GIS.
  • Khẳng định vai trò phòng hộ tuyệt đối của rừng tự nhiên giàu và trung bình với mức mất đất duy trì an toàn dưới ngưỡng 0,8 mm/năm nhờ thảm mục dày và tầng tán đa tầng.
  • Nhận diện chính xác 18,6% diện tích lưu vực đang chịu áp lực xói mòn nghiêm trọng tại các vùng đất trống, nương rẫy canh tác dốc và rừng trồng non chưa khép tán.
  • Xây dựng thành công hệ thống bản đồ số chuyên đề phân cấp xói mòn theo Tiêu chuẩn Quốc gia TCVN 5299:2009, tạo công cụ quản lý trực quan cho chính quyền địa phương.
  • Đề xuất gói giải pháp lâm sinh kỹ thuật đồng bộ, ưu tiên trồng rừng hỗn giao bản địa và khoanh nuôi phục hồi 12.500 ha rừng đầu nguồn trong kế hoạch hành động 2025-2030.

Để bảo vệ bền vững nguồn tài nguyên đất và nguồn nước hạ lưu sông Lam, các cấp quản lý và cộng đồng địa phương cần nhanh chóng triển khai các biện pháp canh tác bảo tồn đất dốc, chấm dứt tình trạng đốt thực bì toàn diện và đẩy mạnh chuyển hóa rừng sản xuất sang rừng phòng hộ nhiều tầng tán ngay hôm nay.