Tổng quan nghiên cứu

Việt Nam sở hữu ba phần tư diện tích tự nhiên là vùng đồi núi, nằm trọn trong vành đai khí hậu nhiệt đới gió mùa và được đánh giá là một trong năm quốc gia chịu tác động nặng nề nhất của biến đổi khí hậu toàn cầu. Hiện tượng thời tiết cực đoan như lũ lụt và hạn hán xuất hiện với tần suất ngày càng dày đặc, gây tổn thất nghiêm trọng về người và tài sản. Bản chất của các hiểm họa này bắt nguồn từ sự mất cân bằng trong động thái thủy văn và biến động dòng chảy bất thường tại các lưu vực sông suối. Mối quan hệ tương tác giữa đặc điểm lưu vực bao gồm lớp phủ thực vật rừng, độ dốc, độ chênh cao, hình dạng lưu vực, diện tích và chế độ mưa có ý nghĩa quyết định đến cơ chế hình thành lũ.

Luận văn tập trung giải quyết vấn đề cốt lõi: Xác định quy luật biến động dòng chảy và định lượng mức độ ảnh hưởng của các đặc trưng lưu vực đến các chỉ số thủy văn quan trọng. Mục tiêu cụ thể là xây dựng hệ thống phương trình tương quan thực nghiệm giữa các yếu tố hình thái lưu vực với hệ số biến động dòng chảy, hệ số tăng lũ, hệ số giảm lũ, thời gian trễ lũ và tổng lưu lượng nước. Nghiên cứu được triển khai trên phạm vi không gian gồm 17 lưu vực vừa và nhỏ phân bố tại bốn vùng sinh thái trọng điểm của Việt Nam từ Bắc Bộ, Bắc Trung Bộ, Nam Trung Bộ đến Tây Nguyên, với diện tích dao động từ 3.343 ha đến 163.559 ha và độ che phủ rừng quy đổi từ 6,8% đến 71,74%. Kết quả nghiên cứu mang ý nghĩa thực tiễn to lớn khi xác lập cơ sở khoa học để hoạch định tỷ lệ diện tích rừng quy đổi tối ưu từ 20% đến 50%, nhằm kiểm soát hệ số biến động dòng chảy năm dưới mức 200%, giảm thiểu nguy cơ lũ quét và tối ưu hóa an ninh nguồn nước quốc gia.

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 tuần hoàn thủy văn rừng của John D. Hewlett và mô hình cân bằng nước lưu vực của Bruijnzeel cùng Thomas Dunne. Quá trình thủy văn lưu vực bắt đầu từ khi nước mưa tiếp xúc với thảm thực vật, diễn ra các quá trình chặn giữ trên tán lá, chảy men thân cây, thấm vào đất, hình thành dòng chảy ngầm, dòng chảy mặt và thoát hơi nước trở lại khí quyển.

Các khái niệm then chốt được lượng hóa trong nghiên cứu gồm:

  • Khả năng giữ nước của tán rừng: Tán cây có khả năng lưu giữ từ 11,4% đến 34,3% tổng lượng mưa đối với rừng nhiệt đới và lên đến 20% đến 40% đối với rừng lá kim, làm suy giảm đáng kể động năng hạt mưa.
  • Lớp thảm mục đáy rừng: Đóng vai trò như màng lọc cơ học và xốp hút nước, có khả năng tích trữ lượng nước gấp 1,9 đến 3,09 lần khối lượng khô của chính nó, giúp đất rừng duy trì tốc độ thấm ổn định trên 80 mm/giờ và tích giữ khoảng 641 đến 679 tấn nước trên mỗi hecta hàng năm.
  • Diện tích rừng quy đổi và độ che phủ rừng quy đổi: Chỉ số phản ánh trọng số chức năng thủy văn của từng trạng thái rừng (rừng giàu có hệ số 1,00; rừng trung bình 0,96; rừng nghèo 0,89; rừng phục hồi 0,81; rừng trồng 0,66).
  • Chỉ số hình dạng lưu vực: Đo lường mức độ kéo dài hay tập trung dòng chảy qua tương quan giữa chu vi và diện tích.

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

Nghiên cứu lựa chọn cỡ mẫu gồm 17 lưu vực vừa và nhỏ đại diện cho các vùng địa hình, khí hậu và thổ nhưỡng điển hình của Việt Nam. Phương pháp chọn mẫu có chủ đích dựa trên các tiêu chí khắt khe: diện tích từ 33,43 km2 đến 1.635,59 km2 (quy mô đủ để hình thành dòng chảy lũ nhưng không quá lớn làm triệt tiêu ảnh hưởng của thảm thực vật), có chuỗi trạm quan trắc khí tượng thủy văn đo đạc liên tục tối thiểu 5 năm và có đầy đủ bản đồ hiện trạng rừng.

Nguồn dữ liệu sử dụng bao gồm:

  • Bản đồ địa hình số hóa tỷ lệ 1:50.000 của Tổng cục Địa chính năm 2004 để trích xuất mô hình số độ cao DEM.
  • Bản đồ hiện trạng rừng tỷ lệ 1:50.000 của Viện Điều tra Quy hoạch Rừng năm 2007.
  • Chuỗi số liệu đo mưa và lưu lượng dòng chảy chi tiết từng giờ trong năm 2007 tại 17 trạm thủy văn tương ứng ở cửa ra lưu vực.

Về phương pháp phân tích, tác giả ứng dụng công nghệ GIS (phần mềm ArcGIS 9.2 với độ phân giải ô lưới raster 30x30m) kết hợp lập trình FoxPro để trích xuất các thông số lưu vực. Phân tích thống kê đơn biến và đa biến trên SPSS 13.0 được lựa chọn nhằm mô hình hóa chính xác mối quan hệ phi tuyến và tuyến tính giữa các biến số hình thái với các chỉ tiêu dòng chảy ở mức ý nghĩa thống kê cao với xác suất sai số p nhỏ hơn 0,05.

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

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

Thứ nhất, lượng mưa và dòng chảy tại các lưu vực phân hóa sâu sắc theo không gian và mùa vụ. Lượng mưa năm biến động mạnh từ 1.144 mm đến hơn 3.500 mm, trong đó các lưu vực Bắc Trung Bộ đạt mức cao nhất với trung bình 3.501,0 mm, Tây Nguyên đạt 2.289,3 mm và Bắc Bộ đạt 1.652,0 mm. Lượng mưa tập trung trong 6 đến 7 tháng mùa mưa, chiếm từ 78% đến 95% tổng lượng mưa cả năm. Tương ứng với đó, lưu lượng dòng chảy mùa lũ chiếm tới 80% tổng lưu lượng năm, tạo nên các đường thủy đồ dạng răng cưa đỉnh nhọn với thời gian tập trung lũ rất ngắn.

Thứ hai, diện tích lưu vực thể hiện tương quan tuyến tính đồng biến rất chặt chẽ với các chỉ số lưu lượng và quy mô lũ. Phương trình hồi quy giữa diện tích lưu vực (DT, tính bằng km2) với tổng lưu lượng dòng chảy năm (Qt, triệu m3) đạt hệ số xác định R2 = 0,7331 qua hàm Qt = 41,2 + 1,733 * DT. Tương tự, hệ số tăng lũ trung bình (Fin, m3/s) và hệ số giảm lũ (Fde, m3/s) cũng tương quan chặt với diện tích lưu vực, lần lượt đạt phương trình Fin = -1,24 + 0,019 * DT (R2 = 0,6283) và Fde = -1,15 + 0,016 * DT (R2 = 0,514).

Thứ ba, thời gian trễ lũ (Lt) giữa đỉnh mưa và đỉnh lũ có sự phân hóa rõ rệt theo khu vực địa lý. Ở các lưu vực miền Bắc và miền Trung, đỉnh lũ thường xuất hiện rất nhanh, chỉ chậm hơn đỉnh mưa từ 1 đến 2 ngày (Lt dao động từ 6 đến 24 giờ). Ngược lại, tại khu vực Tây Nguyên, đỉnh lũ xuất hiện muộn hơn đỉnh mưa từ 15 đến 20 ngày nhờ dung tích trữ nước ngầm rất lớn của tầng đất đỏ bazan dày và hệ thống thảm thực vật rừng tự nhiên.

Thứ tư, độ che phủ rừng quy đổi (FCqđ) tại 17 lưu vực dao động từ 6,8% đến 71,74%. Rừng tự nhiên tập trung chủ yếu ở đai cao từ 500m đến 1.000m (chiếm 46% diện tích rừng) và sườn dốc 15 đến 25 độ (chiếm 36%), đóng vai trò kiểm soát chính đối với hệ số biến động dòng chảy năm (FCV dao động từ 94% đến trên 200%).

Thảo luận kết quả

Cơ chế vật lý đằng sau các phát hiện trên cho thấy diện tích lưu vực càng lớn thì tổng dung lượng đón mưa càng cao, làm gia tăng động năng và lưu lượng đỉnh lũ. Tuy nhiên, sự hiện diện của thảm phủ rừng tự nhiên đóng vai trò như một hồ điều hòa sinh thái khổng lồ. Tán rừng giữ lại từ 10% đến 20% lượng mưa, thảm mục và hệ rễ làm tăng độ rỗng của đất, giúp chuyển hóa dòng chảy mặt thành dòng chảy ngầm ngấm sâu, từ đó làm giảm tốc độ tăng lũ và kéo dài thời gian tập trung nước về lòng dẫn.

Dữ liệu nghiên cứu được biểu diễn trực quan thông qua các bảng ma trận tương quan đa biến và đồ thị hồi quy phân tán (Scatter Plot) với hệ số R2 cao. Biểu đồ thủy đồ quan trắc theo giờ tại các trạm Sơn Diệm, Bình Tường và Thượng Nhật mô tả rõ nét dạng đỉnh nhọn trong các trận mưa lớn, chứng minh rằng tại các lưu vực có độ che phủ rừng thấp dưới 30% và địa hình dốc trên 15 độ, tốc độ tăng lũ diễn ra rất đột ngột. Kết quả này hoàn toàn nhất quán với các công trình nghiên cứu thủy văn rừng trên thế giới của Hewlett và Andreassian, đồng thời củng cố nhận định rằng rừng phòng hộ đầu nguồn có tác dụng làm chậm lũ và điều tiết dòng chảy kiệt hữu hiệu nhất trên quy mô các lưu vực vừa và nhỏ.

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

Thứ nhất, thiết lập và duy trì ngưỡng độ che phủ rừng quy đổi an toàn tối thiểu từ 35% đến 50% trên toàn bộ các lưu vực sông suối có độ dốc trung bình trên 15 độ và lượng mưa năm vượt 2.000 mm. Bộ Nông nghiệp và Phát triển nông thôn cùng Ủy ban nhân dân các tỉnh cần rà soát quy hoạch 3 loại rừng, đặt mục tiêu duy trì hệ số biến động dòng chảy FCV dưới mức 200% và giảm tốc độ tăng lũ từ 15% đến 25% trong giai đoạn 2025-2030.

Thứ hai, nâng cấp mạng lưới trạm quan trắc khí tượng thủy văn tự động và tích hợp mô hình không gian GIS độ phân giải cao 30x30m tại các lưu vực trọng điểm. Cục Khí tượng Thủy văn cần chủ trì triển khai hệ thống cảm biến đo mưa và mực nước theo thời gian thực tại 100% các lưu vực xung yếu miền Trung và Tây Nguyên, nhằm nâng cao độ chính xác của dự báo đỉnh lũ trước thời gian xuất hiện từ 12 đến 24 giờ, hoàn thành trước năm 2027.

Thứ ba, chuyển dịch cơ cấu rừng trồng sang phát triển rừng phòng hộ hỗn loài bản địa nhiều tầng tán tại các đai cao trên 500m. Các Ban quản lý rừng phòng hộ và doanh nghiệp lâm nghiệp cần ưu tiên trồng bổ sung cây bản địa lá rộng, hạn chế tối đa việc khai thác trắng rừng trồng trên sườn dốc để nâng dung tích giữ nước của đất lên ngưỡng 600 đến 680 tấn/ha/năm, thực hiện liên tục trong lộ trình 2025-2035.

Thứ tư, ban hành quy chế phối hợp liên ngành và liên tỉnh trong việc điều tiết xả lũ hồ chứa dựa trên các phương trình dự báo thủy văn lưu vực. Bộ Tài nguyên và Môi trường cùng Ủy ban Quản lý Lưu vực sông cần ứng dụng các thuật toán tương quan giữa diện tích lưu vực và lưu lượng dòng chảy cực đại để vận hành tối ưu hệ thống hồ chứa thủy điện, thủy lợi, phấn đấu giảm thiểu 30% thiệt hại do lũ quét gây ra từ năm 2025 trở đi.

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

Nhóm thứ nhất bao gồm các cơ quan quản lý nhà nước và nhà hoạch định chính sách thuộc ngành Nông nghiệp và Phát triển nông thôn, Tài nguyên và Môi trường. Tài liệu cung cấp cơ sở định lượng để xây dựng quy hoạch rừng phòng hộ đầu nguồn, phân bổ hạn mức đất lâm nghiệp và xây dựng khung chính sách chi trả dịch vụ môi trường rừng gắn với lưu vực thủy điện.

Nhóm thứ hai là các kỹ sư công trình thủy lợi, thủy điện và các chuyên gia tư vấn phòng chống thiên tai. Kết quả từ các hàm hồi quy về tổng lưu lượng dòng chảy và hệ số tăng lũ là thông số đầu vào tin cậy để tính toán lưu lượng thiết kế công trình, thiết kế dung tích phòng lũ cho hồ đập và xây dựng bản đồ phân vùng nguy cơ lũ quét.

Nhóm thứ ba gồm các nhà khoa học, giảng viên, học viên cao học và sinh viên các chuyên ngành Quản lý tài nguyên rừng, Thủy văn học, Biến đổi khí hậu và Khoa học Môi trường. Luận văn là tài liệu tham khảo chuẩn mực về phương pháp kết hợp công nghệ GIS, mô hình số hóa độ cao DEM và phân tích thống kê đa biến SPSS trong nghiên cứu thủy văn rừng.

Nhóm thứ tư là các tổ chức phi chính phủ (NGO) và các ban quản lý dự án phát triển bền vững. Nghiên cứu cung cấp luận cứ khoa học sắc bén để thiết kế các dự án phục hồi cảnh quan rừng, trồng rừng đầu nguồn thích ứng với biến đổi khí hậu tại 4 vùng sinh thái của Việt Nam.

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

Rừng có thể ngăn chặn hoàn toàn nguy cơ xuất hiện lũ lụt trên các lưu vực hay không? Rừng không thể triệt tiêu hoàn toàn lũ lụt trong các trận bão có lượng mưa cực đoan vượt quá dung tích trữ nước bão hòa của đất. Tuy nhiên, rừng tự nhiên nhiều tầng tán giúp giảm tốc độ tăng dòng chảy từ 2,5 đến 27 lần so với đất nông nghiệp và kéo dài thời gian trễ lũ, giúp giảm sức tàn phá của đỉnh lũ.

Diện tích rừng quy đổi mang ý nghĩa gì trong công tác quản lý lưu vực? Diện tích rừng quy đổi là chỉ số chuẩn hóa phản ánh chính xác năng lực điều tiết nước của từng loại rừng, trong đó rừng giàu và rừng trung bình có trọng số cao nhất (từ 0,96 đến 1,00), còn rừng trồng có trọng số thấp hơn (0,66). Chỉ số này giúp nhà quản lý đánh giá thực chất chất lượng phòng hộ của lưu vực thay vì chỉ nhìn vào diện tích che phủ thô.

Vì sao đỉnh lũ tại Tây Nguyên lại xuất hiện muộn hơn đỉnh mưa từ 15 đến 20 ngày? Hiện tượng này xuất phát từ đặc điểm địa chất thổ nhưỡng đặc thù của Tây Nguyên với tầng đất phong hóa bazan rất dày, độ xốp cao và hệ sinh thái rừng tự nhiên còn khá nguyên vẹn. Dung tích thấm giữ nước ngầm khổng lồ đã làm chậm quá trình dồn nước mặt, tạo ra độ trễ dòng chảy kiệt kéo dài.

Nhân tố nào của lưu vực chi phối mạnh nhất đến tổng lưu lượng dòng chảy năm? Diện tích lưu vực và tổng lượng mưa là hai nhân tố chi phối mạnh nhất. Phương trình thực nghiệm chứng minh diện tích lưu vực tương quan tuyến tính rất chặt chẽ với tổng lưu lượng dòng chảy (R2 = 0,7331), trong đó lưu vực càng lớn thì lượng nước sản sinh hàng năm càng cao.

Làm thế nào để ứng dụng các hàm toán học trong luận văn vào cảnh báo sớm thiên tai? Các cơ quan khí tượng có thể nhập trực tiếp số liệu dự báo mưa và diện tích lưu vực vào phương trình hồi quy Fin = -1,24 + 0,019 * DT để tính toán tốc độ dâng của lũ. Kết hợp với chỉ số độ dốc và hình dạng lưu vực, hệ thống sẽ tự động đưa ra cảnh báo sớm về thời gian xuất hiện đỉnh lũ trước 6 đến 24 giờ.

Kết luận

  • Nghiên cứu đã định lượng hóa thành công mối quan hệ giữa đặc trưng hình thái 17 lưu vực điển hình với các chỉ số biến động dòng chảy tại Việt Nam.
  • Thiết lập hệ thống phương trình hồi quy thực nghiệm tin cậy giữa diện tích lưu vực với tổng lưu lượng nước (R2 = 0,7331), hệ số tăng lũ (R2 = 0,6283) và hệ số giảm lũ (R2 = 0,514).
  • Làm sáng tỏ quy luật phân bố dòng chảy theo mùa, trong đó mùa lũ chiếm 80% lượng dòng chảy cả năm và thời gian trễ lũ dao động từ 6 giờ đến 20 ngày tùy vùng sinh thái.
  • Xác lập cơ sở khoa học cho việc duy trì độ che phủ rừng quy đổi từ 20% đến 50% nhằm đảm bảo an toàn sinh thái và phòng chống thiên tai lưu vực.
  • Đề xuất 4 nhóm giải pháp đồng bộ về quy hoạch lâm nghiệp, công nghệ quan trắc GIS và cơ chế quản lý liên vùng cho giai đoạn 2025-2030.

Luận văn là tài liệu khoa học nền tảng có giá trị lý luận và ứng dụng thực tiễn cao cho công tác quản lý tài nguyên rừng và an ninh nguồn nước. Các cơ quan quản lý, nhà nghiên cứu và kỹ sư thủy văn hãy khai thác triệt để các mô hình toán học này vào thực tiễn quy hoạch lưu vực ngay hôm nay.