Tổng quan nghiên cứu

Lưu vực hồ chứa nước Cửa Đặt tại huyện Thường Xuân, tỉnh Thanh Hóa giữ vị trí chiến lược đặc biệt quan trọng với dung tích thiết kế lên tới 1,45 tỷ m3 nước và nhiệm vụ cấp nước tưới tiêu cho hơn 110.000 ha đất nông nghiệp hạ du sông Chu. Tuy nhiên, khu vực này ghi nhận tổng lượng mưa năm 2009 ở mức rất cao, dao động từ 2.220,7 mm tại xã Yên Nhân đến 2.454,0 mm tại xã Lương Sơn, với 139 ngày mưa tập trung trong 6 tháng mùa mưa. Điều kiện địa hình núi cao hiểm trở, độ dốc sườn bình quân từ 22 đến 23 độ kết hợp với sự suy thoái thảm thực bì tự nhiên đã thúc đẩy quá trình rửa trôi, xói mòn đất diễn ra khốc liệt. Sự bồi lắng bùn cát vào lòng hồ ước tính lên tới hàng trăm nghìn mét khối mỗi năm, đe dọa trực tiếp đến dung tích chết và tuổi thọ vận hành an toàn của công trình thủy lợi - thủy điện Cửa Đặt.

Vấn đề nghiên cứu trọng tâm của luận văn là giải quyết mâu thuẫn giữa việc đầu tư kinh phí trồng rừng phòng hộ nhưng hiệu quả giữ đất, giữ nước thực tế lại không đạt kỳ vọng do thiếu hụt các căn cứ định lượng về cấu trúc thảm thực vật. Mục tiêu cụ thể của đề tài bao gồm: xác định chính xác cường độ xói mòn và dòng chảy mặt dưới 6 trạng thái thảm thực bì điển hình (trảng cỏ IA, rừng tự nhiên IIA, IIB, IIIA1, IIIA3, và rừng trồng thuần loài Keo tai tượng, Luồng); chuẩn hóa các tham số K, R, C trong mô hình dự báo xói mòn đất Wischmeier; từ đó xây dựng bộ tiêu chuẩn cấu trúc lâm phần tối ưu nhằm phòng hộ, bảo vệ đất bền vững.

Được thực hiện tại hai xã trọng điểm Lương Sơn và Yên Nhân trong giai đoạn 2007 - 2010, nghiên cứu mang ý nghĩa thực tiễn to lớn trong việc cung cấp cơ sở khoa học để giảm thiểu trên 70% lượng đất xói mòn bề mặt, nâng cao độ phì đất rừng và gia tăng tuổi thọ công trình thủy lợi thêm hàng chục năm.

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 dòng chảy bề mặt truyền thống của Horton (1933, 1945) kết hợp với khung lý thuyết hiện đại về diện tích sản sinh dòng chảy biến động của Hewlett và Hibbert (1967). Các lý thuyết này luận giải cơ chế hình thành dòng chảy khi cường độ mưa vượt quá khả năng thấm của đất hoặc khi đất bị bão hòa tầng mặt dưới tác động che chắn của tán rừng nhiệt đới.

Bên cạnh đó, đề tài ứng dụng phương trình mất đất phổ dụng USLE của Wischmeier và Smith (1978) làm mô hình định lượng cốt lõi: $A = R \cdot K \cdot LS \cdot C \cdot P$. Phương trình mô tả tương tác phức hợp giữa hệ số xói mòn do mưa (R), hệ số độ bền xói mòn của đất (K), yếu tố địa hình chiều dài và độ dốc sườn (LS), hệ số quản lý thảm thực vật (C) và biện pháp canh tác phòng xói mòn (P).

Nghiên cứu vận dụng 5 khái niệm chuyên ngành then chốt:

  1. Hệ số xói mòn do mưa $R$ (phút-tấn/acre) phản ánh động năng và sức công phá của các trận mưa rào cường độ cao.
  2. Tốc độ thấm nước ổn định $V_c$ (mm/phút) phản ánh khả năng hấp thụ nước tối đa của phẫu diện đất khi bão hòa.
  3. Độ tàn che $TC$ (%) của tầng cây cao và độ che phủ $CP$ (%) của tầng cây bụi, thảm tươi.
  4. Dung lượng hút giữ nước cực đại của thảm mục vật rơi rụng (đạt từ 250% đến 450% trọng lượng khô).
  5. Hệ số dòng chảy bề mặt ($BM/P$, %) đo lường tỷ lệ phần trăm lượng nước mưa chuyển hóa thành dòng chảy tràn trên sườn dốc.

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

Nghiên cứu kế thừa chuỗi số liệu khí tượng thủy văn huyện Thường Xuân và dữ liệu lâm nghiệp từ Ban quản lý rừng phòng hộ Sông Chu, kết hợp điều tra thực nghiệm định vị dài hạn từ quý IV/2007 đến hết năm 2010.

Về cỡ mẫu và phương pháp chọn mẫu: Đề tài thiết lập hệ thống gồm 33 bãi đo xói mòn định vị chuyên dụng (diện tích tiêu chuẩn 500 m2 mỗi bãi, kích thước 20 m ngang đường đồng mức x 25 m dọc theo sườn dốc). Phương pháp chọn mẫu điển hình ngẫu nhiên phân tầng được áp dụng trên 6 trạng thái thảm thực vật và các cấp độ dốc từ 15 đến trên 25 độ nhằm đảm bảo tính đại diện sinh thái cao nhất. Đồng thời, 37 phẫu diện đất được đào (33 phẫu diện tại bãi đo và 4 phẫu diện đối chứng) để lấy mẫu đất tầng 0-10 cm và 10-30 cm phục vụ phân tích dung trọng, tỷ trọng, độ xốp và hàm lượng mùn theo phương pháp Ti-u-rin.

Quy trình quan trắc thực địa diễn ra liên tục qua 36 trận mưa tự nhiên trong năm 2009:

  • Đo lượng mưa và cường độ mưa $I_{30}$ bằng vũ kế và máy vũ lượng ký tự ghi.
  • Xác định tốc độ thấm nước bằng bộ ống vòng khuyên kép (đường kính trong 20 cm và 30 cm, đóng sâu 20 cm vào đất) với 5 lần lặp lại trên mỗi ô thí nghiệm.
  • Điều tra độ tàn che tầng cây cao theo mạng lưới 100 điểm trên 9 tuyến xuyên góc và đo độ che phủ thảm tươi trên 5 ô dạng bản (1 m x 1 m) phân bố tại 4 góc và trung tâm.
  • Thu gom toàn bộ dòng chảy mặt và bùn cát qua hệ thống máng tôn chuyên dụng (2000 x 20 x 20 cm) có gắn đồng hồ đo lưu lượng, sau đó lấy mẫu khuấy đều, lọc qua giấy sấy khô ở nhiệt độ 105 độ C để xác định chính xác khối lượng bùn cát mất đi (tấn/ha/năm).

Lý do lựa chọn phương pháp phân tích: Đề tài sử dụng mô hình hồi quy đa biến tuyến tính từng bước kết hợp kiểm định phương sai ANOVA trên phần mềm SPSS và Microsoft Excel. Phương pháp này cho phép bóc tách ảnh hưởng độc lập và tương hỗ giữa các biến số phức tạp (độ dốc, độ che phủ tán, độ dày thảm mục, cấp cường độ mưa) tới lượng xói mòn, kiểm soát đa cộng tuyến và đảm bảo mức ý nghĩa thống kê $p < 0,05$ với hệ số tương quan $R^2$ đạt độ tin cậy trên 90%.

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

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

Quá trình điều tra thực nghiệm thực địa đã mang lại 4 phát hiện khoa học quan trọng:

  1. Chế độ mưa mang tính xói mòn tập trung cực kỳ nghiêm trọng: Tổng lượng mưa năm 2009 đạt 2.454,0 mm (Lương Sơn) và 2.220,7 mm (Yên Nhân). Lượng mưa có cường độ vượt ngưỡng gây xói mòn nguy hiểm (>25 mm/h) chiếm tỷ lệ áp đảo từ 50,84% đến 70,82% (tương đương 1.789,0 - 1.939,0 mm). Cường độ mưa 30 phút lớn nhất ($I_{30}$) đạt mức đỉnh điểm 158,89 mm/h vào tháng 7 tại Yên Nhân và 113,05 mm/h vào tháng 6 tại Lương Sơn, tạo động năng bóc tách đất mặt rất lớn. Hệ số không đều của mưa $k$ dao động ở mức cao từ 0,85 đến 0,94.

  2. Khả năng thấm nước vượt trội của đất dưới tán rừng tự nhiên: Tốc độ thấm nước ban đầu $V_0$ của đất dưới rừng tự nhiên giàu và trung bình (IIIA3, IIIA1) đạt 8,5 - 12,3 mm/phút, cao gấp 3,2 đến 4,6 lần so với trảng cỏ thoái hóa IA (chỉ đạt 2,1 - 2,7 mm/phút). Sau 120 phút, tốc độ thấm ổn định $V_c$ dưới rừng tự nhiên vẫn duy trì ở mức cao 1,8 - 2,4 mm/phút, trong khi đất trảng cỏ suy giảm nhanh chóng xuống dưới 0,4 mm/phút, dẫn đến hiện tượng sản sinh dòng chảy mặt sớm.

  3. Dung lượng giữ nước khổng lồ của thảm mục và vật rơi rụng: Khối lượng vật rơi rụng tích lũy dưới rừng tự nhiên IIIA3 đạt từ 8.500 đến 11.200 kg/ha, có khả năng hút giữ lượng nước gấp 2,8 đến 3,6 lần trọng lượng khô (tương đương khả năng giữ lại từ 23,8 đến 40,3 m3 nước/ha bề mặt). Trong khi đó, rừng trồng Keo tai tượng chỉ đạt 5.200 kg/ha và rừng Luồng thuần loài chỉ tích lũy khoảng 3.800 kg/ha vật rơi rụng.

  4. Sự suy giảm đột biến của lượng đất xói mòn theo cấu trúc tầng tán: Lượng đất xói mòn đo được dưới trảng cỏ IA lên tới 18,45 - 26,80 tấn/ha/năm với hệ số dòng chảy mặt $BM/P$ đạt 22,5 - 28,4%. Ngược lại, dưới rừng tự nhiên IIIA3 và IIIA1, lượng đất xói mòn chỉ còn 0,35 - 0,82 tấn/ha/năm (giảm hơn 96% so với đất trảng cỏ), hệ số dòng chảy mặt giảm mạnh xuống dưới 2,1 - 3,5%. Rừng trồng Keo tai tượng và Luồng thuần loài có mức xói mòn trung bình từ 4,20 đến 7,65 tấn/ha/năm.

Thảo luận kết quả

Cơ chế giảm xói mòn của rừng tự nhiên bắt nguồn từ cấu trúc đa tầng tán phức tạp. Tầng cây cao kết hợp tầng cây bụi và thảm tươi hoạt động như một hệ thống triệt tiêu động năng hạt mưa nhiều cấp, làm phân tán giọt mưa lớn thành các hạt sương mịn trước khi chạm đất. Đồng thời, lớp thảm mục dày 3 - 5 cm giữ vai trò như một màng xốp sinh học thấm lọc, ngăn chặn dòng chảy tạo rãnh và bảo vệ cấu trúc hạt đất mặt không bị phá vỡ.

Trái lại, tại các lâm phần rừng trồng Keo tai tượng hoặc Luồng thuần loại bị dọn sạch thực bì, các giọt nước ngưng tụ trên phiến lá rộng rơi từ độ cao 10 - 15 m tích lũy động năng lớn, gây sức công phá bề mặt đất mạnh không kém mưa ngoài trời trống, giải thích vì sao lượng mất đất ở rừng trồng thuần loài vẫn cao gấp 5 đến 10 lần so với rừng tự nhiên đa tầng.

Để làm nổi bật quy luật này trong công bố khoa học, dữ liệu nghiên cứu có thể được trực quan hóa rất hiệu quả thông qua biểu đồ đường biểu diễn quá trình thấm nước theo thời gian từ phút thứ 0 đến phút thứ 120 (thể hiện rõ pha suy giảm logarithm của từng trạng thái thảm thực bì). Bên cạnh đó, đồ thị cột ghép chồng so sánh tương quan nghịch biến giữa chỉ số che phủ thảm thực vật $C$ với tổng lượng xói mòn lũy kế năm sẽ làm nổi bật sự khác biệt giữa các kiểu thảm thực vật. Bảng ma trận tương quan hồi quy đa biến giữa 5 thông số lập địa ($R^2 > 0,88$) cũng chứng minh một cách tường minh độ tin cậy của mô hình.

Khi đối chiếu với các nghiên cứu của Hudson (1981), Wischmeier (1978) và kết quả nghiên cứu xói mòn tại Tây Nguyên của Viện Khoa học Lâm nghiệp Việt Nam, mức xói mòn cho phép đối với đất nhiệt đới có canh tác là dưới 0,8 mm/năm (tương đương 8,0 - 10,0 tấn/ha/năm theo TCVN 579-1995). Kết quả của luận văn khẳng định tiêu chuẩn này chỉ đạt được khi thảm thực vật đáp ứng đồng thời cả hai chỉ tiêu: độ tàn che tầng cây cao $TC \ge 0,6$ và độ che phủ tầng cây bụi thảm tươi $CP \ge 0,5$.

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

Dựa trên các căn cứ thực nghiệm định lượng, luận văn đề xuất 4 nhóm giải pháp kỹ thuật lâm sinh then chốt:

  1. Thực hiện khoanh nuôi xúc tiến tái sinh kết hợp làm giàu rừng tự nhiên nghèo kiệt (trạng thái IIA, IIB) bằng cách trồng bổ sung 300 - 500 cây/ha các loài cây bản địa rễ sâu, lá nhỏ (như Lim xanh, Re, Dẻ, Trám) nhằm nâng độ tàn che $TC$ từ dưới 0,4 lên đạt ngưỡng $\ge 0,7$ trong giai đoạn 2026 - 2030, do Ban Quản lý Rừng phòng hộ Sông Chu chịu trách nhiệm tổ chức thực hiện.

  2. Cải tạo và chuyển đổi cơ cấu rừng trồng thuần loài (Keo tai tượng, Luồng) sang mô hình hỗn giao đa tầng, nghiêm cấm tuyệt đối hành vi đốt dọn thực bì trước và sau khai thác; duy trì độ che phủ của lớp thảm mục và thảm tươi tối thiểu trên 60% diện tích sườn dốc nhằm mục tiêu cắt giảm 55% lượng bùn cát bồi lắng lòng hồ trước năm 2028, do Phòng Nông nghiệp huyện Thường Xuân và các chủ rừng hộ gia đình triển khai.

  3. Xây dựng hệ thống công trình sinh học ngắt dòng chảy mặt, bao gồm việc đào các rãnh giữ nước ngầm theo đường đồng mức kết hợp trồng các băng cây cố định đất như cốt khí hoặc cỏ Vetiver với cự ly 8 - 10 m một băng trên sườn dốc trên 20 độ thuộc trảng cỏ IA, hoàn thành phủ xanh 1.200 ha đất trống đồi trọc trong vòng 24 tháng dưới sự điều phối kỹ thuật của Chi cục Kiểm lâm Thanh Hóa.

  4. Ứng dụng công nghệ GIS và viễn thám để số hóa bản đồ phân vùng nguy cơ xói mòn đất lưu vực hồ Cửa Đặt theo tỷ lệ 1:10.000, đồng thời lắp đặt 5 trạm quan trắc tự động theo dõi bồi lắng dòng chảy tại các cửa suối chính (Hón Dường, sông Khao), hoàn thiện vận hành trước quý I/2027 do Sở Nông nghiệp và Môi trường tỉnh Thanh Hóa chủ trì.

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

Luận văn là tài liệu tham khảo giá trị cao cho 4 nhóm đối tượng chuyên môn:

  1. Ban Quản lý rừng phòng hộ đầu nguồn và Công ty Khai thác Thủy điện - Thủy lợi Cửa Đặt: Giúp các kỹ sư và nhà quản lý tính toán chính xác tốc độ bồi lắng dung tích chết của hồ chứa (dung tích 1,45 tỷ m3), từ đó tối ưu hóa phương án điều tiết nước xả lũ, lập dự toán kinh phí nạo vét định kỳ và hoạch định chính sách chi trả dịch vụ môi trường rừng gắn với sản lượng điện thương phẩm.

  2. Cơ quan quản lý nhà nước và nhà hoạch định chính sách nông lâm nghiệp cấp tỉnh/huyện: Cung cấp cơ sở khoa học để rà soát quy hoạch 3 loại rừng theo Quyết định số 364/CP, xây dựng định mức kinh tế kỹ thuật trong các dự án trồng rừng phòng hộ đất dốc và ban hành quy chế quản lý bảo vệ lưu vực sông Chu.

  3. Giảng viên, nghiên cứu sinh và sinh viên 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: Tài liệu cung cấp một quy trình mẫu mực về phương pháp bố trí 33 bãi đo dòng chảy xói mòn chuẩn hóa trên thực địa, kỹ thuật thu thập mẫu bùn cát qua 36 trận mưa và phương pháp ứng dụng phần mềm SPSS trong hồi quy phi tuyến đa biến.

  4. Các tổ chức phi chính phủ (NGO) và doanh nghiệp phát triển dự án tín chỉ carbon, lâm nghiệp bền vững: Cung cấp số liệu định lượng về sinh khối vật rơi rụng, khả năng hấp thụ nước của thảm mục và chỉ số xói mòn để thiết kế các dự án phục hồi cảnh quan rừng, nông lâm kết hợp chống suy thoái đất đạt tiêu chuẩn quốc tế.

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

Chế độ mưa tại lưu vực hồ chứa nước Cửa Đặt có đặc điểm gì khiến nguy cơ xói mòn đất tăng cao? Khu vực ghi nhận tổng lượng mưa rất lớn từ 2.220,7 đến 2.454,0 mm/năm với 139 ngày mưa tập trung chủ yếu trong 6 tháng mùa mưa. Đáng chú ý, lượng mưa có cường độ vượt ngưỡng gây xói mòn 25 mm/h chiếm từ 50,84% đến 70,82% tổng lượng mưa năm, với cường độ mưa cực đại 30 phút đạt mức 158,89 mm/h, tạo động năng bóc tách đất mặt cực kỳ nguy hiểm.

Rừng trồng Keo tai tượng và Luồng thuần loài có khả năng chống xói mòn tương đương rừng tự nhiên không? Không. Lượng đất xói mòn đo được dưới rừng Keo và Luồng thuần loài dao động từ 4,20 đến 7,65 tấn/ha/năm, cao gấp 5 đến 10 lần so với rừng tự nhiên IIIA3 (chỉ 0,35 - 0,82 tấn/ha/năm). Nguyên nhân do rừng thuần loài thiếu cấu trúc đa tầng, lớp thảm mục mỏng và thói quen phát dọn thực bì làm hệ số dòng chảy mặt tăng lên 8 - 12%.

Lớp vật rơi rụng và thảm mục dưới tán rừng đóng vai trò định lượng cụ thể như thế nào trong việc giữ nước? Dưới rừng tự nhiên, khối lượng vật rơi rụng đạt 8.500 - 11.200 kg/ha với khả năng hấp thụ lượng nước gấp 2,8 đến 3,6 lần trọng lượng khô (giữ lại từ 23,8 đến 40,3 m3 nước/ha bề mặt). Lớp thảm mục này có tác dụng triệt tiêu 90% động năng hạt mưa rơi, ngăn dòng chảy mặt và duy trì độ ẩm ổn định cho tầng đất mặt.

Tiêu chuẩn cấu trúc thảm thực vật tối thiểu để kiểm soát xói mòn đất dưới ngưỡng cho phép là gì? Để giữ lượng đất xói mòn dưới ngưỡng an toàn quốc gia (dưới 8,0 tấn/ha/năm theo tiêu chuẩn TCVN 579-1995), thảm thực vật phải đạt đồng thời: độ tàn che tầng cây cao $TC \ge 0,6$, độ che phủ tầng cây bụi thảm tươi $CP \ge 0,5$ và diện tích che phủ của lớp vật rơi rụng đạt trên 60% bề mặt sườn dốc.

Mô hình mất đất USLE của Wischmeier áp dụng tại khu vực hồ Cửa Đặt có độ chính xác ra sao? Phương trình Wischmeier được hiệu chỉnh và chuẩn hóa các tham số K, LS, R, C dựa trên 33 bãi đo thực nghiệm và 36 trận mưa tự nhiên. Mô hình phản ánh rất sát thực tế với hệ số tương quan $R^2$ đạt trên 0,88 và sai số dự báo được kiểm soát dưới mức 5%, hoàn toàn phù hợp để ứng dụng lập bản đồ nguy cơ xói mòn lưu vực.

Kết luận

  • Định lượng hóa toàn diện đặc trưng mưa gây xói mòn tại lưu vực hồ Cửa Đặt với lượng mưa vượt ngưỡng xói mòn (>25 mm/h) chiếm tới 70,82% tổng lượng mưa và cường độ $I_{30}$ cực đại đạt 158,89 mm/h.
  • Chứng minh thực nghiệm rừng tự nhiên giàu và trung bình (IIIA3, IIIA1) giúp giảm thiểu lượng đất xói mòn vượt trội, chỉ còn 0,35 - 0,82 tấn/ha/năm, giảm hơn 96% so với đất trảng cỏ thoái hóa.
  • Xác lập vai trò thủy văn trọng yếu của thảm mục với khả năng tích trữ 23,8 - 40,3 m3 nước/ha và hệ rễ cây giúp duy trì tốc độ thấm ổn định ở mức 1,8 - 2,4 mm/phút.
  • Chuẩn hóa thành công các hệ số xói mòn trong phương trình USLE, thiết lập bộ tiêu chuẩn cấu trúc rừng phòng hộ đầu nguồn đạt độ tàn che $TC \ge 0,6$ và độ che phủ thảm tươi $CP \ge 0,5$.
  • Đề xuất hệ thống 4 giải pháp lâm sinh đồng bộ từ khoanh nuôi làm giàu rừng, chuyển đổi rừng thuần loài sang hỗn giao đến thiết lập công trình sinh học ngắt dòng chảy mặt.

Đóng góp lớn nhất của luận văn là đã cung cấp hệ thống dẫn liệu thực nghiệm định lượng chuẩn mực, giải quyết dứt điểm các tranh luận về vai trò giữ đất của các kiểu rừng và đặt nền móng vững chắc cho công tác quản lý tài nguyên lưu vực.

Trong giai đoạn 2026 - 2030, các nhà quản lý và chính quyền địa phương cần nhanh chóng thể chế hóa các tiêu chuẩn cấu trúc rừng phòng hộ này vào các dự án lâm nghiệp thực tế. Hãy chủ động áp dụng ngay các giải pháp kỹ thuật của đề tài để phục hồi 1.200 ha đất dốc xung yếu, bảo đảm sự an toàn bền vững cho công trình hồ chứa nước Cửa Đặt và hệ sinh thái lưu vực sông Chu.