Tổng quan nghiên cứu

Lưu vực sông 3S bao gồm ba dòng sông chính là Sê Kông, Sê San và Sêrêpôk, trải rộng trên tổng diện tích tự nhiên 78.650 km², phân bố qua ba quốc gia Việt Nam (chiếm 38%), Campuchia (chiếm 33%) và Lào (chiếm 29%). Đây là lưu vực xuyên biên giới quy mô lớn nhất trong toàn bộ hệ thống lưu vực sông Mê Kông, đóng vai trò như một huyết mạch sinh thái và thủy văn sống còn. Lưu vực sông 3S cung cấp khoảng 20% tổng lượng dòng chảy mặt và đóng góp từ 15% đến 40% tổng lượng phù sa bùn cát cho toàn bộ dòng chính sông Mê Kông. Nguồn tài nguyên nước và phù sa từ lưu vực này đóng vai trò quyết định trong việc duy trì an ninh sinh kế cho hơn 4 triệu cư dân bản địa, đồng thời trực tiếp nuôi dưỡng hệ sinh thái Biển Hồ Tonle Sap của Campuchia và vựa lúa Đồng bằng sông Cửu Long của Việt Nam với hơn 18 triệu dân.

Tuy nhiên, lưu vực sông 3S đang phải đối mặt với nhiều áp lực nghiêm trọng từ sự suy giảm độ che phủ rừng trung bình 0,53%/năm, sự phát triển của hơn 41 dự án thủy điện với tổng công suất tăng gấp 8 lần so với thập niên 1990, và việc mở rộng diện tích đất nông nghiệp lên tới hơn 378.000 ha. Biến đổi khí hậu toàn cầu làm gia tăng các hiện tượng thời tiết cực đoan, đe dọa trực tiếp đến an ninh nguồn nước và chất lượng môi trường nước. Mục tiêu trọng tâm của nghiên cứu là ứng dụng mô hình thủy văn bán phân bố SWAT nhằm đánh giá toàn diện ảnh hưởng của biến đổi khí hậu theo các kịch bản nồng độ đại diện đến chất lượng nước, tập trung vào ba chỉ tiêu chính gồm tổng chất rắn lơ lửng (TSS), tổng Nitơ (T-N) và tổng Photpho (T-P). Nghiên cứu cung cấp cơ sở khoa học định lượng giúp các nhà quản lý xây dựng chiến lược quản lý tổng hợp tài nguyên nước xuyên biên giới hiệu quả và bền vững.

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 tích hợp của các lý thuyết thủy văn lưu vực, động lực học dòng chảy và quá trình chuyển hóa chất ô nhiễm phi điểm trong môi trường đất - nước. Trọng tâm của khung lý thuyết là mô hình thủy văn SWAT (Soil and Water Assessment Tool), được phát triển bởi Cục Nghiên cứu Nông nghiệp thuộc Bộ Nông nghiệp Hoa Kỳ. SWAT vận hành dựa trên phương trình cân bằng nước tổng thể, mô phỏng chu trình thủy văn thông qua hai pha chính: pha đất điều khiển lượng nước, tải lượng trầm tích, dinh dưỡng chảy vào mạng lưới sông; và pha nước điều khiển quá trình truyền tải các chất qua mạng lưới kênh dẫn và hồ chứa.

Trong phân tích chất lượng nước, mô hình tích hợp phương trình mất đất phổ dụng cải tiến (MUSLE) để tính toán xói mòn và tải lượng tổng chất rắn lơ lửng (TSS), kết hợp với các phương trình động học chuyển hóa dinh dưỡng của QUAL2E để mô phỏng sự dịch chuyển và chu trình biến đổi của tổng Nitơ (T-N) và tổng Photpho (T-P). Về mặt khí hậu, nghiên cứu áp dụng khung lý thuyết đánh giá biến đổi khí hậu của Ban Liên chính phủ về Biến đổi Khí hậu (IPCC) trong Báo cáo Đánh giá lần thứ năm (AR5), cụ thể là kịch bản đường dẫn nồng độ đại diện RCP 4.5. Lưu vực được chi tiết hóa thành các đơn vị phản ứng thủy văn (HRU) dựa trên sự đồng nhất về loại đất, trong đó nhóm đất Acrisols chiếm 63,14% và nhóm đất Ferralsols chiếm 18,18% trên tổng số 12 nhóm đất chính của lưu vực, kết hợp với hiện trạng sử dụng đất và các cấp độ dốc địa hình.

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

Cơ sở dữ liệu của nghiên cứu bao gồm chuỗi số liệu quan trắc khí tượng thủy văn liên tục 28 năm (giai đoạn 1981–2008) thu thập từ hệ thống trạm quan trắc chuẩn của Ủy hội sông Mê Kông (MRC) và mạng lưới khí tượng quốc gia ba nước. Dữ liệu không gian đầu vào bao gồm mô hình số độ cao (DEM) độ phân giải chuẩn, bản đồ phân loại thảm phủ sử dụng đất và bản đồ thổ nhưỡng tỷ lệ 1:250.000 theo phân loại của FAO/UNESCO. Phương pháp chọn mẫu dữ liệu được thực hiện có hệ thống trên 8 trạm thủy văn và trạm quan trắc chất lượng nước tiêu biểu phân bố đồng đều từ thượng nguồn đến hạ lưu, bao gồm Voeunsai, Chantangoy, Bankamphun, Bản Đôn, Cầu 14, Siempang, Lumphat và Pleiku.

Lý do lựa chọn mô hình SWAT tích hợp phần mềm SWAT-CUP với thuật toán tối ưu hóa SUFI-2 là khả năng xử lý các bài toán thủy văn quy mô lớn với cấu trúc dữ liệu không gian phức tạp và điều kiện địa hình chia cắt mạnh. Quá trình mô hình hóa trải qua ba bước nghiêm ngặt: phân tích độ nhạy của 19 thông số thủy văn và chất lượng nước, hiệu chỉnh tự động và kiểm định độc lập mô hình theo dữ liệu lưu lượng ngày và dữ liệu chất lượng nước tháng. Phương pháp chi tiết hóa thống kê biến đổi khí hậu sử dụng phương pháp thay đổi hệ số Delta kết hợp với 5 mô hình hoàn lưu khí quyển toàn cầu (GCM) nhằm chuyển đổi các dự báo khí hậu toàn cầu thành dữ liệu nhiệt độ và lượng mưa cục bộ chi tiết cho từng tiểu lưu vực thuộc hệ thống sông Sê Kông, Sê San và Sêrêpôk.

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

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

Kết quả hiệu chỉnh và kiểm định mô hình SWAT cho thấy mức độ phù hợp và độ tin cậy vượt trội trong việc mô phỏng quá trình thủy văn và chất lượng nước trên lưu vực sông 3S. Trong mô phỏng dòng chảy mặt, giai đoạn hiệu chỉnh ghi nhận hệ số xác định R² đạt trên 0,70 và chỉ số hiệu quả Nash-Sutcliffe (NSE) đạt trên 0,63 trên toàn bộ 4 trạm thủy văn chính gồm Voeunsai (sông Sê San), Chantangoy (sông Sê Kông), Bankamphun và Bản Đôn (sông Sêrêpôk). Giai đoạn kiểm định dòng chảy tiếp tục khẳng định tính ổn định cao của mô hình với R² vượt mức 0,69 và NSE đạt trên 0,67.

Đối với các chỉ tiêu chất lượng nước, mô hình phản ánh sát thực tế biến động tải lượng ô nhiễm. Tại giai đoạn hiệu chỉnh, chỉ số R² đạt trên 0,78 và NSE đạt trên 0,74 tại hai trạm trọng điểm là Bản Đôn và Siempang cho cả ba thông số TSS, T-N và T-P. Trong giai đoạn kiểm định tại trạm Cầu 14 và Siempang, các giá trị R² và NSE duy trì ở mức cao, lần lượt đạt trên 0,68 và 0,62. Khi áp dụng kịch bản biến đổi khí hậu RCP 4.5, nghiên cứu phát hiện xu thế tải lượng TSS gia tăng rõ rệt vào mùa mưa trên toàn lưu vực do lượng mưa cực đoan tăng mạnh thúc đẩy quá trình rửa trôi đất mặt. Đồng thời, tải lượng T-N và T-P đều ghi nhận xu hướng tăng liên tục trong cả mùa mưa lẫn mùa khô tại cả ba phân lưu vực Sê Kông, Sê San và Sêrêpôk.

Thảo luận kết quả

Diễn biến tải lượng ô nhiễm có thể được trực quan hóa thông qua các biểu đồ đường thủy văn so sánh giữa giá trị quan trắc và giá trị mô phỏng theo trục thời gian, kết hợp cùng bảng ma trận biến thiên tải lượng ô nhiễm qua các mùa. Nguyên nhân cốt lõi khiến tải lượng TSS tăng đột biến trong mùa mưa bắt nguồn từ sự gia tăng cường độ mưa kết hợp với hiện tượng suy thoái rừng đầu nguồn (giảm 0,53%/năm), khiến lưu tốc dòng chảy mặt tăng và gây xói mòn nghiêm trọng các tầng đất đỏ Bazan màu mỡ.

Sự gia tăng của Nitơ tổng và Photpho tổng trong mùa khô được giải thích bởi hiện tượng suy giảm lưu lượng dòng chảy tự nhiên kết hợp với việc gia tăng sử dụng phân bón hóa học trong canh tác nông nghiệp (diện tích nông nghiệp dự kiến vượt 378.000 ha), dẫn đến việc cô đặc các hợp chất dinh dưỡng trong lòng dẫn. Khi so sánh với các nghiên cứu thủy văn tại lưu vực sông Bé có chỉ số R² đạt 0,77 và NSE đạt 0,75 hoặc nghiên cứu tại lưu vực Cedar Creek có NSE đạt từ 0,51 đến 0,66, kết quả mô phỏng trên lưu vực 3S cho thấy độ tương thích vượt bậc. Điều này khẳng định rằng sự kết hợp giữa biến đổi khí hậu và áp lực chuyển đổi mục đích sử dụng đất đang làm thay đổi căn bản cơ chế truyền tải dinh dưỡng và trầm tích, đe dọa trực tiếp chất lượng nguồn nước hạ lưu.

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

Dựa trên các phát hiện định lượng, nghiên cứu đề xuất 4 nhóm giải pháp chiến lược nhằm kiểm soát ô nhiễm và bảo vệ bền vững tài nguyên nước lưu vực sông 3S:

  • Thiết lập và phục hồi dải đệm thực vật ven sông kết hợp mô hình nông lâm kết hợp tại các tiểu lưu vực có độ dốc trên 15 độ, nhằm đạt mục tiêu giảm từ 25% đến 30% tải lượng bồi lắng TSS trước năm 2030, do Sở Nông nghiệp và Phát triển Nông thôn các tỉnh Tây Nguyên chủ trì phối hợp với Ban Quản lý Rừng phòng hộ thực hiện.
  • Áp dụng các quy trình thực hành nông nghiệp tốt (GAP) và quản lý dinh dưỡng cây trồng tích hợp, cắt giảm từ 15% đến 20% lượng phân bón hóa học chứa Nitơ và Photpho thất thoát vào nguồn nước mặt trước năm 2028, do Trung tâm Khuyến nông các địa phương hướng dẫn các hộ nông dân và hợp tác xã triển khai.
  • Hoàn thiện và ban hành quy trình vận hành liên hồ chứa thủy điện linh hoạt trong điều kiện biến đổi khí hậu, bắt buộc duy trì dòng chảy môi trường mùa kiệt đạt tối thiểu 80% giá trị tự nhiên trung bình nhiều năm trong giai đoạn 2026–2030, do Bộ Tài nguyên và Môi trường cùng các chủ đầu tư thủy điện phối hợp thực hiện.
  • Xây dựng mạng lưới trạm quan trắc chất lượng nước tự động và cơ chế chia sẻ thông tin thủy văn xuyên biên giới thời gian thực giữa Việt Nam, Lào và Campuchia trước năm 2027, dưới sự điều phối chính của Ủy hội sông Mê Kông quốc tế (MRC) và Ủy ban sông Mê Kông Việt Nam.

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

Nội dung và kết quả nghiên cứu của luận văn mang lại giá trị khoa học và ứng dụng thực tiễn cao cho 4 nhóm đối tượng trọng tâm:

  • Cơ quan quản lý nhà nước về tài nguyên và môi trường (Sở Tài nguyên và Môi trường, Cục Quản lý Tài nguyên nước, Ủy ban sông Mê Kông Việt Nam): Khai thác bộ cơ sở dữ liệu định lượng về tải lượng TSS, T-N, T-P để xây dựng quy hoạch phân bổ tài nguyên nước lưu vực sông và hoạch định chính sách thích ứng biến đổi khí hậu xuyên biên giới.
  • Các nhà khoa học, nghiên cứu sinh và giảng viên chuyên ngành Môi trường, Thủy văn, Biến đổi khí hậu và GIS: Sử dụng luận văn như tài liệu tham khảo chuẩn mực về quy trình thiết lập, hiệu chỉnh và kiểm định mô hình SWAT kết hợp công cụ SWAT-CUP (thuật toán SUFI-2) cho các lưu vực sông lớn xuyên biên giới.
  • Các đơn vị tư vấn kỹ thuật và doanh nghiệp vận hành công trình thủy điện, thủy lợi: Ứng dụng kết quả dự báo biến động dòng chảy và bồi lắng trầm tích để tối ưu hóa quy trình xả lũ, bảo dưỡng công trình và lập báo cáo đánh giá tác động môi trường định kỳ.
  • Các tổ chức phi chính phủ (NGO) và tổ chức phát triển quốc tế hoạt động trong lĩnh vực bảo tồn tài nguyên nước: Sử dụng bằng chứng khoa học từ nghiên cứu để thiết kế các dự án viện trợ sinh kế bền vững, bảo tồn hệ sinh thái rừng đầu nguồn và nâng cao năng lực ứng phó thiên tai cho cộng đồng địa phương.

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

Mô hình SWAT có ưu điểm gì vượt trội so với các mô hình thủy văn khác khi áp dụng cho lưu vực sông 3S?
Mô hình SWAT cho phép phân chia lưu vực rộng lớn 78.650 km² thành các tiểu lưu vực và đơn vị phản ứng thủy văn HRU chi tiết dựa trên địa hình, thổ nhưỡng và sử dụng đất. SWAT có khả năng mô phỏng liên tục dòng chảy và chất lượng nước trong điều kiện thiếu hụt số liệu đo đạc liên tục, đồng thời tích hợp công cụ SWAT-CUP giúp tự động hóa quá trình hiệu chỉnh và kiểm định thông số với độ chính xác cao (R² > 0,78).

Kịch bản biến đổi khí hậu RCP 4.5 tác động cụ thể ra sao đến tổng chất rắn lơ lửng (TSS) trên lưu vực?
Theo kịch bản RCP 4.5, lượng mưa cực đoan trong mùa mưa gia tăng đáng kể trên cả ba phân lưu vực Sê Kông, Sê San và Sêrêpôk. Cường độ mưa lớn làm tăng năng lượng bóc mòn đất mặt và lưu tốc dòng chảy tràn, dẫn đến sự gia tăng mạnh mẽ của tải lượng TSS trôi về hệ thống sông chính, làm gia tăng hiện tượng bồi lắng và đục hóa nguồn nước mặt.

Tại sao nồng độ Nitơ tổng (T-N) và Photpho tổng (T-P) lại có xu hướng tăng trong cả mùa khô?
Trong mùa khô, nhiệt độ tăng cao làm tăng bốc thoát hơi nước trong khi lượng mưa suy giảm nghiêm trọng, dẫn đến lưu lượng dòng chảy mặt suy giảm mạnh. Sự sụt giảm thể tích nước làm suy yếu khả năng tự làm sạch và pha loãng của dòng sông, kết hợp với các hoạt động xả thải và dư lượng phân bón nông nghiệp tích tụ, khiến nồng độ T-N và T-P bị cô đặc và gia tăng rõ rệt.

Thuật toán SUFI-2 trong công cụ SWAT-CUP đóng vai trò gì trong nghiên cứu này?
Thuật toán SUFI-2 (Sequential Uncertainty Fitting ver. 2) thực hiện phân tích độ nhạy của 19 thông số thủy văn và tự động tối ưu hóa các giá trị tham số trong không gian đa chiều. Công cụ này giúp xác định dải độ không đảm bảo của mô phỏng, đảm bảo kết quả hiệu chỉnh đạt chỉ số NSE > 0,74 và R² > 0,78 cho chất lượng nước mà không phụ thuộc vào phương pháp hiệu chỉnh thủ công tốn thời gian.

Kết quả nghiên cứu này đóng góp gì cho việc bảo vệ Đồng bằng sông Cửu Long?
Lưu vực sông 3S cung cấp 20% lượng nước và 15% đến 40% lượng phù sa cho hạ lưu sông Mê Kông. Kết quả dự báo của luận văn giúp các cơ quan quản lý Việt Nam chủ động dự báo lượng phù sa bồi đắp, nguy cơ sạt lở bờ sông và xâm nhập mặn tại Đồng bằng sông Cửu Long, từ đó xây dựng kế hoạch điều tiết nước và bảo vệ an ninh lương thực bền vững cho hơn 18 triệu dân.

Kết luận

  • Nghiên cứu đã ứng dụng thành công mô hình thủy văn bán phân bố SWAT kết hợp công cụ SWAT-CUP để mô phỏng dòng chảy và chất lượng nước (TSS, T-N, T-P) trên toàn bộ lưu vực sông 3S rộng 78.650 km².
  • Độ tin cậy của mô hình được kiểm chứng vượt trội qua các chỉ số thống kê tiêu chuẩn, đạt R² > 0,69 và NSE > 0,67 cho dòng chảy, cùng R² > 0,68 và NSE > 0,62 cho các chỉ tiêu chất lượng nước tại các trạm quan trắc trọng điểm.
  • Dưới tác động của kịch bản biến đổi khí hậu RCP 4.5, tải lượng TSS có xu thế gia tăng mạnh trong mùa mưa, trong khi tải lượng Nitơ tổng và Photpho tổng gia tăng trong cả mùa mưa và mùa khô trên cả ba lưu vực sông Sê Kông, Sê San và Sêrêpôk.
  • Luận văn đã cung cấp hệ thống 4 nhóm giải pháp kỹ thuật, phi công trình và quản lý xuyên biên giới với lộ trình cụ thể nhằm kiểm soát xói mòn và suy thoái chất lượng nước lưu vực.
  • Đóng góp khoa học then chốt của công trình là thiết lập quy trình công nghệ tích hợp mô hình hóa và GIS chuẩn xác, tạo tiền đề vững chắc cho công tác quản lý lưu vực sông xuyên biên giới trong bối cảnh biến đổi khí hậu.

Trong giai đoạn 2026–2030, các cơ quan quản lý và các tổ chức nghiên cứu cần nhanh chóng chuyển giao kết quả mô hình vào thực tiễn quản lý lưu vực, tích hợp hệ thống dữ liệu SWAT vào nền tảng giám sát thông minh của Ủy hội sông Mê Kông, nhằm chủ động bảo vệ nguồn tài nguyên nước và phát triển bền vững khu vực Tam giác phát triển Việt Nam - Lào - Campuchia.