Tổng quan nghiên cứu

Việt Nam có hơn 3/4 diện tích tự nhiên là đồi núi, trong đó trên 60% diện tích có độ dốc vượt quá 15°. Tình trạng khai thác đất dốc thiếu quy hoạch đã khiến khoảng 11 triệu ha đất trên cả nước bị xói mòn và thoái hóa nghiêm trọng. Tại khu vực Tây Bắc, lưu vực hồ thủy điện Nậm Chiến (thuộc huyện Mường La, tỉnh Sơn La) có đặc thù địa hình chia cắt mạnh với hơn 90% diện tích sở hữu độ dốc trên 8° và phân bố ở dải độ cao từ 900 m đến 1.400 m so với mực nước biển. Lượng mưa trung bình năm tại đây đạt khoảng 1.200 mm, trong đó 76% lượng mưa tập trung dồn dập vào các tháng 6, 7 và 8, tạo nên áp lực xói mòn và rửa trôi đất mặt vô cùng lớn.

Vấn đề nghiên cứu trọng tâm xuất phát từ tập quán canh tác nương rẫy đốt rừng của cộng đồng dân cư địa phương (chiếm hơn 90% là đồng bào dân tộc thiểu số), làm suy giảm nghiêm trọng độ che phủ rừng (vốn chỉ đạt khoảng 45,96% trên toàn huyện). Mục tiêu cụ thể của luận văn là:

  • Đánh giá định lượng đặc điểm cấu trúc và khả năng chống xói mòn, bảo vệ đất của 7 trạng thái thảm thực vật chính trên lưu vực;
  • Phân tích thực trạng xói mòn tiềm năng theo từng cấp độ dốc và độ cao;
  • Đề xuất hệ thống giải pháp quy hoạch sử dụng đất và kỹ thuật lâm sinh bền vững.

Nghiên cứu được triển khai thực nghiệm toàn diện trên lưu vực hồ thủy điện Nậm Chiến trong năm 2011. Kết quả nghiên cứu có ý nghĩa quyết định trong việc thiết lập ranh giới an toàn sinh thái, duy trì cường độ xói mòn đất dưới ngưỡng cho phép là 0,8 mm/năm (tương đương 8 tấn/ha/năm), góp phần giảm thiểu bồi lắng lòng hồ, bảo vệ tuổi thọ công trình thủy điện và ổn định sinh kế cho người dân địa phương.

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

Khung lý thuyết áp dụng

Luận văn vận dụng lý thuyết động năng hạt mưa của Ellison (1945) và Hudson (1981), xác định xói mòn đất diễn ra qua hai pha chính: pha hạt mưa bắn phá phá vỡ kết cấu đất và pha dòng chảy mặt cuốn trôi các hạt đất. Việc giảm thiểu động năng hạt mưa thông qua tán lá và thảm mục đóng vai trò quyết định trong bảo vệ đất.

Nghiên cứu tích hợp phương trình mất đất phổ dụng USLE/RUSLE (Wischmeier & Smith, 1978; Foster và cộng sự, 1982) kết hợp mô hình dự báo xói mòn đất nhiệt đới của Trường Đại học Lâm nghiệp (Vương Văn Quỳnh, 1994, 1999). Phương trình toán học biểu diễn mối quan hệ đa biến:

$$d = f(\alpha, TC, H, CP, TM, X, K)$$

Trong đó, 5 khái niệm chuyên ngành cốt lõi bao gồm:

  • Độ tàn che tầng cây cao ($TC$): Tỷ lệ che bóng của tán lá tầng cây ưu thế đối với mặt đất;
  • Cường độ xói mòn ($d$): Lượng đất bị mất đi tính bằng mm/năm hoặc tấn/ha/năm;
  • Độ che phủ thảm mục ($TM$): Tỷ lệ che phủ bề mặt đất của tầng cành rơi, lá rụng khô kiệt;
  • Dung trọng đất ($D$) và Tỷ trọng đất ($d_t$): Khối lượng thể tích tự nhiên và khối lượng riêng của hạt đất tính bằng g/cm³;
  • Độ xốp tầng đất mặt ($X%$): Tỷ lệ phần trăm thể tích các khe hở trong đất, được xác định qua công thức $X = \frac{d_t - D}{d_t} \times 100$.

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

Nguồn dữ liệu ngoại nghiệp được thu thập thông qua phương pháp điều tra thực nghiệm kết hợp phương pháp kế thừa số liệu khí tượng thủy văn giai đoạn 2005 - 2010.

Phương pháp chọn mẫu điển hình được áp dụng để thiết lập cỡ mẫu gồm 60 ô tiêu chuẩn (OTC) có diện tích 1.000 m² (kích thước 33 × 30 m, chiều dài bố trí chạy dọc theo đường đồng mức) tại các dải độ cao từ 900 m đến 1.388 m và độ dốc dao động từ 5° đến 42°. Cơ cấu 60 OTC phân bố đều trên 7 trạng thái: 12 ô rừng trồng (thông, trẩu, hỗn giao), 10 ô rừng nghèo, 15 ô rừng trung bình, 9 ô rừng phục hồi, 2 ô rừng tre nứa, 4 ô đất nông nghiệp và 8 ô đất trống có cây bụi.

Trong mỗi OTC:

  • Điều tra tầng cây cao bằng thước đo cao Blum-Leiss và thước vanh đo đường kính tại vị trí 1,3 m; đo độ tàn che bằng mạng lưới 80 điểm ngẫu nhiên hệ thống;
  • Thiết lập 5 ô tiêu chuẩn dạng bản (5 × 5 m) tại 4 góc và tâm để đo thảm tươi cây bụi;
  • Thiết lập 25 ô nhỏ kích thước 1 × 1 m để thu thập toàn bộ vật rơi rụng sấy khô kiệt xác định khối lượng thảm mục;
  • Đào 1 phẫu diện đất chuẩn kích thước $0,8 \times 1,5 \times 1,2\text{ m}$, lấy mẫu đất nguyên dạng bằng ống dung trọng 100 cm³ ở các tầng sâu cách nhau 20 cm (0-20 cm, 20-40 cm, 40-60 cm...).

Lý do lựa chọn phương pháp này là nhằm đảm bảo tính đại diện sinh thái cao nhất, thu thập chính xác các thông số vi cấu trúc thực bì và lý tính đất. Toàn bộ số liệu được phân tích thống kê hồi quy trên phần mềm Excel, tính toán chỉ số xói mòn mưa $K$ qua phần mềm Sinh khí hậu và số hóa bản đồ không gian xói mòn bằng phần mềm MapInfo 10.0 trong suốt quá trình hoàn thiện đề tài năm 2011.

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

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

Phân tích số liệu thực nghiệm từ 60 ô tiêu chuẩn đã chỉ ra những khác biệt mang tính quy luật giữa các trạng thái thảm phủ:

  • Khả năng hạn chế xói mòn vượt trội của rừng tự nhiên: Cường độ xói mòn trung bình dưới tán rừng trung bình và rừng phục hồi duy trì ở mức rất thấp, dao động từ 0,05 mm/năm đến 0,35 mm/năm (tương đương 0,5 - 3,5 tấn/ha/năm), hoàn toàn nằm dưới ngưỡng giới hạn cho phép 0,8 mm/năm. Ngược lại, tại trạng thái đất nông nghiệp canh tác nương rẫy và đất trống cây bụi, cường độ xói mòn tăng vọt lên mức 15,0 - 24,5 mm/năm (150 - 245 tấn/ha/năm), làm mất đi từ 0,9 cm đến 2,1 cm tầng đất mặt phì nhiêu mỗi năm.
  • Tác động của cấu trúc tầng tán và thảm mục: Lớp thảm khô dưới tán rừng tự nhiên có khả năng ngậm nước gấp 1,38 đến 3,21 lần trọng lượng khô của nó, giữ lại lượng nước tương đương 35.840 lít/ha sau một trận mưa 3,6 mm. Khi độ tàn che tầng cây cao ($TC$) giảm từ 0,8 xuống 0,3 và tỷ lệ che phủ thảm mục ($TM$) giảm dưới 30%, lượng đất bị xói mòn tăng mạnh thêm 42,2% và dòng chảy mặt tăng 30,4%.
  • Biến động tính chất vật lý đất: Đất dưới tán rừng trung bình có độ xốp tầng mặt đạt từ 58% đến 65%, dung trọng đạt mức thấp từ 0,95 g/cm³ đến 1,12 g/cm³, tạo điều kiện thấm nước tối đa. Trong khi đó, đất nương rẫy thoái hóa bị nén chặt với dung trọng tăng lên trên 1,38 g/cm³ và độ xốp suy giảm xuống dưới 42%, khiến hệ số thấm nước mặt giảm đi đáng kể.
  • Tác động khuếch đại của độ dốc: Trên các sườn dốc từ 25° đến 42° không có thảm thực vật bảo vệ, cường độ xói mòn đất tăng gấp 3,2 lần so với các sườn dốc dưới 15°, đẩy lưu vực vào tình trạng xói mòn tiềm năng đặc biệt nghiêm trọng.

Thảo luận kết quả

Nguyên nhân cốt lõi dẫn đến sự khác biệt về lượng mất đất giữa các trạng thái là do sự suy thoái tầng tán và lớp phủ thảm mục. Khi mất đi lớp che phủ thực vật, hàng triệu hạt mưa với đường kính từ 2 mm đến 5 mm rơi tự do từ độ cao lớn mang động năng cực đại, bắn phá làm vỡ vụn các kết cấu hạt đất mặt, tạo màng sét bít kín các khe hở mao dẫn, từ đó chuyển hóa phần lớn lượng mưa thành dòng chảy tràn cuốn trôi đất.

So sánh với các nghiên cứu của Võ Đại Hải (1996), Thái Phiên và Nguyễn Tử Siêm (1999) hay Trần An Phong (1996), kết quả nghiên cứu tại lưu vực Nậm Chiến hoàn toàn tương thích và làm rõ hơn quy luật thực bì nhiệt đới vùng cao. Nghiên cứu cũng xác nhận rằng việc trồng thuần loài các cây bản lá rộng như trẩu có lượng đất xói mòn cao hơn 263% so với trồng cây bản lá kim hoặc hỗn giao như thông ba lá và keo lá tràm, do giọt nước tích tụ trên phiến lá to tạo thành giọt lớn có động năng rơi mạnh hơn.

Trong báo cáo phân tích, các dữ liệu không gian được thể hiện trực quan thông qua Bản đồ phân vùng cường độ xói mòn đất tỷ lệ 1/25.000 trên phần mềm MapInfo 10.0 và Biểu đồ hồi quy phi tuyến biểu diễn mối quan hệ giữa ba trục tham số: độ dốc địa hình (từ 5° đến 42°), chỉ số cấu trúc bảo vệ đất $C_1$ và lượng mất đất thực tế. Dữ liệu bảng tổng hợp dung trọng ở các độ sâu (0-20 cm, 20-40 cm, 40-60 cm) giúp lượng hóa chính xác dung tích hấp thu nước của toàn lưu vực.

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

Dựa trên kết quả đánh giá khả năng bảo vệ đất, 4 nhóm giải pháp chiến lược và kỹ thuật được đề xuất nhằm kiểm soát xói mòn lưu vực hồ thủy điện Nậm Chiến:

  1. Quy hoạch không gian sử dụng đất theo phân cấp độ dốc:
    • Hành động: Chuyển đổi bắt buộc toàn bộ diện tích đất nương rẫy canh tác nông nghiệp có độ dốc trên 25° sang mục tiêu lâm nghiệp; khoanh nuôi xúc tiến tái sinh tự nhiên kết hợp làm giàu rừng.
    • Chỉ tiêu (Target Metric): Giảm diện tích đất nông nghiệp dốc nguy hiểm xuống dưới 10%, nâng độ che phủ rừng phòng hộ toàn lưu vực lên mức trên 60%.
    • Thời gian & Chủ thể: Thực hiện trong giai đoạn 2012 - 2016; chủ trì bởi Ủy ban nhân dân huyện Mường La phối hợp với Sở Nông nghiệp và PTNT tỉnh Sơn La.
  2. Áp dụng đồng bộ kỹ thuật canh tác đất dốc bền vững (SALT):
    • Hành động: Thiết lập các băng cây xanh chắn đất bằng cây phân xanh họ Đậu (Cốt khí, đậu Triều) hoặc cỏ Vetiver theo đường đồng mức; kết hợp tạo bờ bậc thang cục bộ trên các diện tích đất nương rẫy có độ dốc từ 15° đến 25°.
    • Chỉ tiêu (Target Metric): Giảm thiểu 75 - 85% lượng đất trôi mặt ruộng, giữ cường độ xói mòn nương rẫy dưới mức 2,5 tấn/ha/năm.
    • Thời gian & Chủ thể: Triển khai liên tục hàng năm từ vụ xuân hè 2012; Trạm Khuyến nông huyện Mường La và UBND xã Ngọc Chiến trực tiếp hướng dẫn kỹ thuật cho các hộ dân.
  3. Cải thiện cấu trúc rừng và nâng cao chất lượng rừng phòng hộ:
    • Hành động: Nghiêm cấm khai thác trắng; áp dụng biện pháp lâm sinh tỉa thưa gián tiếp, trồng bổ sung các loài cây bản địa hỗn giao nhiều tầng tán (thông, dẻ, re) thay thế mô hình độc canh cây lá rộng rụng lá mùa khô.
    • Chỉ tiêu (Target Metric): Đảm bảo duy trì độ tàn che tầng cây cao $TC \ge 0,7$, độ che phủ tầng cây bụi thảm tươi $CP \ge 60%$, thảm mục phủ kín $\ge 70%$ mặt đất.
    • Thời gian & Chủ thể: Giai đoạn 2012 - 2020; do Ban Quản lý Rừng phòng hộ Mường La và các chủ rừng hộ gia đình thực hiện.
  4. Thực thi chính sách chi trả dịch vụ môi trường rừng (DVMTR):
    • Hành động: Trích lập quỹ dịch vụ môi trường rừng từ nguồn thu bán điện thương phẩm của Nhà máy Thủy điện Nậm Chiến để chi trả tiền khoán bảo vệ rừng trực tiếp cho cộng đồng bản địa.
    • Chỉ tiêu (Target Metric): 100% diện tích rừng lưu vực được giao khoán bảo vệ có hỗ trợ kinh phí, chấm dứt hoàn toàn hiện tượng đốt nương làm rẫy trái phép.
    • Thời gian & Chủ thể: Ký kết cam kết và giải ngân thường niên bắt đầu từ năm 2012; Quỹ Bảo vệ và Phát triển rừng tỉnh Sơn La cùng Ban Điều hành Thủy điện Nậm Chiến phối hợp giám sát.

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

Công trình nghiên cứu mang giá trị thực tiễn và học thuật cao, là tài liệu tham khảo đắc lực cho 4 nhóm đối tượng chính:

  • Cơ quan quản lý nhà nước về Nông nghiệp và Phát triển nông thôn: Các cán bộ quy hoạch tại Sở Nông nghiệp & PTNT, Chi cục Kiểm lâm các tỉnh miền núi phía Bắc có thể sử dụng dữ liệu luận văn để xây dựng tiêu chí phân loại rừng phòng hộ xung yếu và lập đề án bảo vệ đất dốc cấp tỉnh.
  • Ban quản lý các dự án Thủy điện và Thủy lợi: Các kỹ sư công trình thủy văn có thể khai thác mô hình toán dự báo lượng bùn cát lắng đọng lòng hồ, từ đó tính toán chính xác dung tích chết, dung tích hữu ích và tuổi thọ vận hành kinh tế của các nhà máy thủy điện vừa và nhỏ.
  • Giảng viên, nghiên cứu sinh và học viên cao học chuyên ngành Lâm nghiệp: Tài liệu cung cấp quy trình phương pháp luận chuẩn mực về thiết lập 60 ô tiêu chuẩn đại diện, kỹ thuật lấy mẫu phẫu diện đất chuyên sâu ở độ sâu 1,2 m và phương pháp ứng dụng phương trình xói mòn cải tiến của Đại học Lâm nghiệp.
  • Các tổ chức phi chính phủ (NGO) và cán bộ khuyến nông địa phương: Cung cấp cơ sở khoa học để thiết kế các mô hình sinh kế nông lâm kết hợp, chuyển giao công nghệ canh tác chống xói mòn trên sườn dốc cho đồng bào các dân tộc thiểu số vùng Tây Bắc.

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

1. Cơ sở khoa học nào xác định giới hạn xói mòn cho phép là dưới 0,8 mm/năm? Ngưỡng 0,8 mm/năm (tương đương khoảng 8 tấn/ha/năm) được xác lập dựa trên tốc độ hình thành đất tự nhiên từ quá trình phong hóa đá mẹ trong điều kiện nhiệt đới ẩm có hoạt động canh tác theo tiêu chuẩn của Hudson (1981). Khi cường độ mất đất nhỏ hơn tốc độ sinh đất, tầng đất mặt mới duy trì được độ phì nhiêu và tính bền vững lâu dài.

2. Vì sao rừng tự nhiên trung bình lại chống xói mòn tốt hơn nhiều so với rừng trồng thuần loài? Rừng tự nhiên trung bình sở hữu cấu trúc không gian nhiều tầng liên tục (tầng vượt tán, tầng ưu thế, tầng cây bụi thảm tươi và lớp thảm khô dày ngậm nước tới 35.840 lít/ha). Cấu trúc này triệt tiêu hoàn toàn động năng hạt mưa nhiều cấp và giữ độ xốp đất mặt đạt trên 60%, trong khi rừng trồng thuần loài thường có kết cấu tán đơn điệu và thảm mục nghèo nàn.

3. Mức độ dốc ảnh hưởng như thế nào đến tốc độ xói mòn đất tại lưu vực Nậm Chiến? Tại lưu vực Nậm Chiến, với hơn 90% diện tích có độ dốc trên 8°, khi độ dốc sườn núi tăng từ 10° lên 30°, vận tốc và lưu lượng dòng chảy mặt tăng theo cấp số nhân, khiến lượng đất xói mòn bị cuốn trôi tăng gấp 2 đến 3,5 lần trên các diện tích đất không có lớp phủ thực bì dày bảo vệ.

4. Điểm vượt trội của phương trình xói mòn Đại học Lâm nghiệp so với phương trình USLE là gì? Phương trình do Trường Đại học Lâm nghiệp phát triển đã tích hợp trực tiếp các biến số cấu trúc vi khí hậu rừng nhiệt đới như độ tàn che ($TC$), tỷ lệ thảm mục ($TM$), thảm tươi ($CP$) và độ xốp đất mặt ($X$). Điều này giúp mô hình hóa chính xác điều kiện rừng tự nhiên Việt Nam mà không bị sai số lớn như phương trình USLE nguyên bản (vốn xây dựng chủ yếu cho đất nông nghiệp đồng bằng Mỹ).

5. Giải pháp lâm sinh nào là cấp thiết nhất cho vùng đầu nguồn thủy điện Nậm Chiến? Giải pháp cấp thiết nhất là khoanh nuôi bảo vệ nghiêm ngặt các diện tích rừng phục hồi hiện có, kết hợp trồng dặm bổ sung các loài cây bản địa hỗn giao nhiều tầng tán trên các diện tích nương rẫy có độ dốc trên 25° để nhanh chóng thiết lập lại độ tàn che tầng tán trên 0,7 và độ che phủ thảm mục trên 70%.

Kết luận

  • Luận văn đã hoàn thành đánh giá thực nghiệm toàn diện trên 60 ô tiêu chuẩn đại diện, lượng hóa chính xác khả năng bảo vệ đất của 7 trạng thái thảm phủ tại lưu vực hồ thủy điện Nậm Chiến.
  • Khẳng định vai trò sinh thái cốt lõi của thảm thực vật rừng tự nhiên và lớp thảm mục trong việc giảm thiểu hơn 90% lượng đất xói mòn so với đất canh tác nương rẫy.
  • Xác định nguyên nhân xói mòn nghiêm trọng tại địa phương bắt nguồn từ tập quán canh tác đốt nương trên sườn dốc lớn (>25°) và suy giảm tầng tán rừng.
  • Ứng dụng thành công hệ thống thông tin địa lý (MapInfo 10.0) kết hợp phần mềm Sinh khí hậu để xây dựng bản đồ phân vùng xói mòn tiềm năng phục vụ quản lý lưu vực.
  • Đề xuất hệ thống 4 nhóm giải pháp đồng bộ kết hợp hài hòa giữa quy hoạch không gian, kỹ thuật lâm sinh và cơ chế tài chính chi trả dịch vụ môi trường rừng.

Đóng góp lớn nhất của luận văn là cung cấp bộ cơ sở dữ liệu thực nghiệm chuẩn xác và luận cứ khoa học tin cậy cho bài toán quản lý tài nguyên đất dốc vùng thủy điện Tây Bắc. Lộ trình tiếp theo đòi hỏi các cơ quan chức năng huyện Mường La đưa các khuyến nghị quy hoạch vào kế hoạch phát triển kinh tế - xã hội địa phương giai đoạn 2012 - 2020.

Các cơ quan quản lý, doanh nghiệp thủy điện và cộng đồng địa phương cần tiếp tục chung tay phối hợp, mở rộng nghiên cứu định vị dài hạn và triển khai quyết liệt các mô hình lâm sinh bền vững nhằm bảo vệ nguồn tài nguyên đất và nguồn nước vô giá cho tương lai.