Tổng quan nghiên cứu

Lưu vực sông và hồ chứa thủy điện đóng vai trò then chốt trong việc bảo đảm an ninh năng lượng và cân bằng hệ sinh thái khu vực miền núi. Tuy nhiên, tình trạng suy thoái thảm phủ thực vật do mở rộng đất sản xuất nông nghiệp đang làm gia tăng mức độ xói mòn và đẩy nhanh tốc độ bồi lắng tại nhiều công trình thủy điện trọng điểm. Tại tỉnh Lâm Đồng, dù tỷ lệ che phủ rừng đạt khoảng 61,5% với hệ thống 57 dự án thủy điện vừa và nhỏ, chất lượng rừng đầu nguồn vẫn đối mặt với nhiều áp lực khai thác.

Vấn đề nghiên cứu trọng tâm tập trung vào việc đánh giá biến động lớp phủ thực vật, xác định mức độ mất đất và dự báo nguy cơ bồi lắng lòng hồ chứa. Mục tiêu cụ thể là xây dựng giải pháp công nghệ tích hợp hệ thống thông tin địa lý (GIS), kỹ thuật viễn thám và mô hình toán phương trình mất đất phổ dụng (USLE) nhằm cung cấp cơ sở khoa học định lượng phục vụ công tác quản lý lưu vực thượng nguồn hồ thủy điện Đa Nhim.

Phạm vi không gian nghiên cứu bao quát toàn bộ lưu vực thượng nguồn với tổng diện tích tự nhiên 73.860 ha, nằm trên địa bàn 2 huyện Lạc Dương (chiếm 77% diện tích với 56.813 ha), Đơn Dương (chiếm 13% diện tích với 9.447 ha) và thành phố Đà Lạt (chiếm 10% diện tích với 7.600 ha). Hồ thủy điện Đa Nhim sở hữu dung tích trữ nước 165 triệu m3 cùng công suất thiết kế 160MW, giữ vị trí then chốt ở bậc thang thượng nguồn sông Đa Dâng - Đa Nhim - Đồng Nai. Khung thời gian đối sánh dữ liệu tập trung vào chu kỳ 10 năm giữa hai mốc năm 2001 và năm 2010.

Ý nghĩa thực tiễn của công trình thể hiện qua việc lượng hóa nguy cơ bồi lắng lên đến 99.101 tấn phù sa mỗi năm và thất thoát 19.104 tấn carbon hữu cơ hàng năm. Dữ liệu số hóa này hỗ trợ trực tiếp cho các nhà hoạch định chính sách xây dựng phương án chi trả dịch vụ môi trường rừng, bảo vệ độ phì nhiêu của đất và kéo dài tuổi thọ vận hành công trình.

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 của 3 khối lý thuyết chuyên sâu: lý thuyết phản xạ phổ thực vật trong viễn thám, lý thuyết xói mòn đất do nước của phương trình mất đất phổ dụng (USLE) do Wischmeier và Smith phát triển, cùng nguyên lý lan truyền dòng chảy và bồi lắng trầm tích hồ chứa kế thừa từ mô hình SWAT.

Hệ thống khái niệm chính bao gồm:

  • Chỉ số thực vật khác biệt chuẩn hóa (NDVI): Phản ánh trạng thái sinh trưởng và sinh khối xanh dựa trên đặc tính hấp thụ ánh sáng đỏ tại kênh 3 (bước sóng 0,63 - 0,69 µm) và phản xạ mạnh ở dải cận hồng ngoại tại kênh 4 (bước sóng 0,76 - 0,90 µm) của hệ thống Landsat.
  • Mô hình mất đất phổ dụng USLE: Mô phỏng lượng đất mất bình quân hàng năm qua công thức tích hợp 5 biến số gồm hệ số xói mòn do mưa (R), hệ số độ bền xói mòn của đất (K), hệ số chiều dài và độ dốc sườn đồi (LS), hệ số quản lý thảm phủ thực vật (C), và hệ số biện pháp canh tác bảo vệ đất (P).
  • Cấp độ xói mòn đất: Phân loại tiềm năng xói mòn từ mức nhẹ đến rất mạnh theo quy chuẩn kỹ thuật quốc gia TCVN 5299:1995.
  • Hàm lượng chất rắn lơ lửng (SS): Nồng độ hạt phù sa mịn trôi nổi trong nước mặt, đại diện cho cường độ bào mòn tầng đất mặt.

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

Nguồn dữ liệu của đề tài bao gồm ảnh viễn thám vệ tinh Landsat đa thời đoạn chụp năm 2001 và năm 2010, hệ thống bản đồ số độ cao (DEM) trích xuất từ đường bình độ do Sở Tài nguyên và Môi trường tỉnh Lâm Đồng cung cấp, bản đồ thổ nhưỡng, bản đồ quy hoạch 3 loại rừng và chuỗi số liệu khí tượng thủy văn giai đoạn 2000 - 2009.

Cỡ mẫu phân tích không gian bao phủ toàn bộ 73.860 ha lưu vực, trong đó mẫu thống kê chi tiết biến động thảm phủ thực vật thực hiện trên diện tích 31.550 ha. Mẫu khảo sát thực địa được thu thập tại 4 tiểu lưu vực đại diện cho các phương thức sử dụng đất điển hình: rừng kín lá rộng thường xanh, rừng thông ba lá thuần loài, rừng thông xen kẽ đất nông nghiệp và vùng canh tác nông nghiệp tập trung.

Phương pháp chọn mẫu thực địa tuân thủ nguyên tắc phân tầng theo đặc trưng thổ nhưỡng và dạng lập địa địa hình dốc, kết hợp lấy mẫu nước mặt trực tiếp để phân tích chỉ tiêu SS trong phòng thí nghiệm. Lý do lựa chọn giải pháp tích hợp phần mềm ENVI 4.5 và ArcGIS 9.3 là khả năng tự động hóa việc chồng lớp ma trận dữ liệu raster, tính toán hệ số LS từ mô hình mạng lưới tam giác không đều (TIN) với độ chính xác cao, giúp giải quyết triệt để sự phức tạp của địa hình chia cắt mạnh trên cao nguyên Lang Biang. Tiến trình nghiên cứu được triển khai từ tháng 10 năm 2012 đến tháng 02 năm 2013.

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

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

Quá trình phân tích dữ liệu không gian và thực nghiệm hiện trường mang lại 4 phát hiện quan trọng:

Thứ nhất, phân tích biến động thảm phủ thực vật trên diện tích thống kê 31.550 ha cho thấy xu hướng biến động rõ nét trong giai đoạn 2001 - 2010. Tính toán cho toàn lưu vực dự báo trong chu kỳ 10 năm, diện tích rừng hỗn giao và rừng lá rộng bị suy giảm khoảng 88 ha, trong khi diện tích rừng thông giảm tới 4.506 ha do chịu sức ép chuyển đổi sang đất sản xuất nông nghiệp.

Thứ hai, mức độ xói mòn đất có sự phân hóa cực kỳ lớn giữa các loại hình che phủ. Năm 2010, mức xói mòn bình quân tính ở cận trên của phân cấp xói mòn trên diện tích rừng hỗn giao là 0,41 tấn/ha/năm, rừng lá rộng là 0,43 tấn/ha/năm, rừng thông đạt 0,51 tấn/ha/năm. Ngược lại, mức độ xói mòn trên đất sản xuất nông nghiệp năm 2001 lên tới 10,68 tấn/ha/năm, cao gấp hơn 20 lần so với đất có rừng thông che phủ.

Thứ ba, mô hình dự báo lượng bồi lắng xác định tổng lượng phù sa đổ về lòng hồ thủy điện Đa Nhim hàng năm đạt mức 99.101 tấn/năm trong kịch bản xói mòn cực đại. Đi kèm với lượng đất bị rửa trôi là sự thất thoát dinh dưỡng nghiêm trọng với 19.104 tấn carbon hữu cơ bị cuốn trôi mỗi năm khỏi bề mặt lưu vực.

Thứ tư, kết quả phân tích nồng độ chất rắn lơ lửng (SS) tại các nhánh sông suối ghi nhận hàm lượng phù sa từ lưu vực nông nghiệp thuần túy cao gấp 31 lần lưu vực rừng lá rộng tự nhiên, gấp 11 lần lưu vực rừng thông thuần loại và gấp 3 lần khu vực rừng thông xen kẽ đất sản xuất nông nghiệp.

Thảo luận kết quả

Nguyên nhân chính dẫn đến sự chênh lệch lớn về mức độ mất đất bắt nguồn từ sự suy giảm lớp phủ thực vật và tập quán canh tác nông nghiệp trên sườn dốc vượt quá 15 độ. Tán rừng tự nhiên có khả năng triệt tiêu khoảng 30% động năng va đập của hạt mưa, kết hợp lớp thảm mục và hệ thống rễ cây chằng chịt giữ ẩm, làm tăng độ thấm nước và duy trì độ bền liên kết của nhóm đất Humic Acrisol. Khi chuyển đổi sang đất nông nghiệp mà không áp dụng biện pháp che phủ, dòng chảy tràn mặt đất tăng vọt, cuốn trôi lớp mùn màu mỡ.

Kết quả định lượng phù sa của đề tài có sự tương đồng sâu sắc với công trình công bố năm 2011 của nhà khoa học Lee H. MacDonald khi nghiên cứu dòng chảy mặt và bồi lắng phục vụ cơ chế chi trả dịch vụ môi trường tại Lâm Đồng. Sự phù hợp này khẳng định độ tin cậy của phương trình USLE tích hợp viễn thám.

Để phục vụ quản lý trực quan, dữ liệu xói mòn được thể hiện qua các bản đồ chuyên đề phân vùng theo tiêu chuẩn quốc gia và hệ thống biểu đồ cột so sánh trực tiếp tỷ lệ phần trăm xói mòn giữa năm 2001 và 2010. Bảng ma trận chồng lớp chuyển đổi hiện trạng sử dụng đất giúp nhà quản lý xác định chính xác các điểm nóng sạt lở tại các tiểu khu lâm nghiệp thuộc huyện Lạc Dương.

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

Dựa trên kết quả nghiên cứu định lượng, 4 nhóm giải pháp trọng tâm được kiến nghị thực hiện:

Thứ nhất, thực hiện quy hoạch và phục hồi rừng đầu nguồn. Ban Quản lý rừng phòng hộ Đa Nhim cùng Vườn Quốc gia Bidoup - Núi Bà cần khẩn trương khoanh nuôi tái sinh 4.506 ha rừng thông đang bị suy thoái, thiết lập hành lang đai rừng phòng hộ ven hồ chứa nhằm mục tiêu cắt giảm ít nhất 40% lượng bùn cát bồi lắng về hồ trong giai đoạn 2014 - 2020.

Thứ hai, áp dụng đồng bộ các kỹ thuật canh tác nông nghiệp bền vững trên đất dốc. Sở Nông nghiệp và Phát triển nông thôn phối hợp chính quyền địa phương chuyển giao kỹ thuật canh tác theo đường đồng mức, xây dựng hệ thống mương cắt dòng phủ cỏ và thiết lập các dải băng phân xanh chắn đất tại các vùng dốc từ 15 đến 20 độ tại xã Đa Sar và Đa Chais. Giải pháp này hướng tới hạ mức xói mòn đất nông nghiệp từ 10,68 tấn/ha/năm xuống dưới 3,0 tấn/ha/năm trước năm 2018.

Thứ ba, hoàn thiện cơ chế xã hội hóa lâm nghiệp và chi trả dịch vụ môi trường rừng (PES). Ủy ban nhân dân tỉnh Lâm Đồng và Công ty Thủy điện Đa Nhim cần sử dụng nguồn thu từ 160MW công suất thương phẩm để tăng mức kinh phí khoán bảo vệ rừng cho hơn 90% hộ đồng bào dân tộc thiểu số tại chỗ, tạo sinh kế ổn định lâu dài và hạn chế áp lực phá rừng làm rẫy.

Thứ tư, xây dựng hệ thống quan trắc và cảnh báo xói mòn bằng công nghệ WebGIS. Sở Tài nguyên và Môi trường tỉnh Lâm Đồng định kỳ 6 tháng một lần cập nhật dữ liệu viễn thám đa thời gian để theo dõi diễn biến thảm phủ, vận hành cơ sở dữ liệu số nhằm dự báo sớm nguy cơ bồi lắng và biến động lòng hồ.

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

Luận văn mang lại giá trị ứng dụng thiết thực cho 4 nhóm đối tượng:

  • Cơ quan quản lý nhà nước về tài nguyên và lâm nghiệp: Sở Nông nghiệp và Phát triển nông thôn, Sở Tài nguyên và Môi trường, Ban Quản lý Vườn Quốc gia Bidoup - Núi Bà có thể khai thác bộ dữ liệu không gian 73.860 ha để rà soát quy hoạch 3 loại rừng và xây dựng định mức chi trả dịch vụ môi trường rừng chuẩn xác.
  • Doanh nghiệp quản lý và vận hành nhà máy thủy điện: Các kỹ sư vận hành nhà máy thủy điện Đa Nhim và cụm thủy điện lưu vực sông Đồng Nai sử dụng số liệu dự báo 99.101 tấn phù sa/năm để chủ động lập kế hoạch nạo vét, tối ưu hóa dung tích 165 triệu m3 và kéo dài tuổi thọ thiết bị cơ điện.
  • Cán bộ khuyến nông và chuyên gia nông nghiệp bền vững: Tài liệu cung cấp cơ sở kỹ thuật để lựa chọn loại cây trồng có hệ số che phủ C cao, thiết kế ruộng bậc thang và hướng dẫn người dân canh tác giảm thiểu bạc màu đất đồi.
  • Giảng viên, học viên và nghiên cứu sinh chuyên ngành Môi trường: Luận văn là tài liệu tham khảo mẫu mực về quy trình xử lý ảnh viễn thám trên ENVI kết hợp thuật toán tính hệ số LS trong ArcGIS và mô hình USLE phục vụ các đề tài đánh giá tai biến địa chất.

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

Tại sao giải pháp tích hợp GIS, viễn thám và mô hình USLE lại mang lại hiệu quả vượt trội trong đánh giá lưu vực Đa Nhim? Sự kết hợp này cho phép tính toán tự động các hệ số địa hình phức tạp trên diện tích rộng 73.860 ha từ mô hình số độ cao DEM và ảnh vệ tinh Landsat. Phương pháp này rút ngắn hơn 60% thời gian khảo sát thực địa, giảm chi phí nhân lực và cho phép cập nhật định kỳ biến động thảm phủ với độ tin cậy khoa học cao.

Thảm phủ thực vật có vai trò định lượng cụ thể ra sao trong việc hạn chế xói mòn đất? Số liệu nghiên cứu chứng minh độ che phủ thực vật giữ đất vượt trội: rừng hỗn giao chỉ mất 0,41 tấn/ha/năm và rừng thông mất 0,51 tấn/ha/năm, thấp hơn 20 lần so với đất nông nghiệp (10,68 tấn/ha/năm). Tán rừng làm giảm 30% động năng hạt mưa rơi, ngăn dòng chảy mặt bóc tách tầng đất.

Mô hình USLE sử dụng các nhóm thông số đầu vào nào để tính toán lượng đất mất? Phương trình USLE tích hợp 5 hệ số tự nhiên và nhân tác gồm: hệ số xói mòn do mưa R (tính từ lượng mưa bình quân năm 1.500 mm), hệ số xói mòn đất K theo từng nhóm thổ nhưỡng, hệ số độ dài - độ dốc sườn LS, hệ số thảm phủ C và hệ số bảo vệ đất P.

Vì sao hàm lượng phù sa tại vùng canh tác nông nghiệp lại cao gấp 31 lần vùng rừng lá rộng? Do vùng nông nghiệp thường xuyên bị cày xới, thiếu tầng cây che bóng và canh tác trên địa hình dốc đứng từ 15 đến 20 độ. Khi có mưa lớn kéo dài, đất tơi xốp nhanh chóng bị cuốn trôi trực tiếp vào các suối Đa Lang Biang và Krông Klet, đẩy nồng độ chất rắn lơ lửng tăng vọt.

Kết quả nghiên cứu hỗ trợ gì cho chính sách chi trả dịch vụ môi trường rừng tại Lâm Đồng? Công trình cung cấp chứng cứ khoa học chuẩn xác về việc 65.062 ha đất có rừng giúp ngăn chặn 99.101 tấn phù sa lắng đọng lòng hồ mỗi năm. Đây là căn cứ pháp lý để chủ đầu tư thủy điện trích nguồn thu từ 160MW công suất để chi trả kinh phí bảo vệ rừng cho người dân bản địa.

Kết luận

  • Hệ thống hóa hoàn chỉnh quy trình công nghệ kết hợp viễn thám Landsat, GIS và mô hình USLE để đánh giá biến động môi trường lưu vực.
  • Lượng hóa chi tiết tốc độ suy thoái với dự báo mất 4.506 ha rừng thông trong 10 năm và bồi lắng 99.101 tấn phù sa hàng năm về hồ chứa.
  • Chứng minh tính ưu việt của hệ sinh thái rừng tự nhiên với mức xói mòn chỉ khoảng 0,41 - 0,51 tấn/ha/năm so với 10,68 tấn/ha/năm của đất nông nghiệp.
  • Cung cấp cơ sở dữ liệu số hóa tin cậy trên 73.860 ha phục vụ công tác điều hành quy hoạch 3 loại rừng và chính sách chi trả dịch vụ môi trường rừng.
  • Đề xuất hệ thống giải pháp công trình và phi công trình đồng bộ với lộ trình thực hiện cụ thể trong giai đoạn 2014 - 2020.

Đóng góp khoa học then chốt của luận văn là thiết lập phương pháp luận định lượng xói mòn chi phí thấp, linh hoạt, phù hợp với đặc thù vùng đồi núi Tây Nguyên. Bước đi tiếp theo là mở rộng tích hợp cảm biến viễn thám độ phân giải cao và mở rộng mô hình cho toàn bộ hệ thống lưu vực sông Đồng Nai trong 2 năm tới. Các nhà quản lý lưu vực và đơn vị vận hành thủy điện cần sớm áp dụng quy trình này để đảm bảo khai thác tài nguyên nước an toàn và bền vững.