Tổng quan nghiên cứu

Biến đổi khí hậu đang tạo ra những áp lực chưa từng có đối với hệ thống an ninh tài nguyên nước toàn cầu, đặc biệt tại Việt Nam nơi có hơn 70% dân số sinh sống tập trung tại các vùng châu thổ và lưu vực sông. Thượng lưu vực sông Cầu đóng vai trò là nguồn cấp nước sinh hoạt, nông nghiệp và công nghiệp trọng yếu cho vùng kinh tế Đông Bắc, với tổng diện tích toàn lưu vực đạt 6.030 km² và tổng lượng dòng chảy hàng năm đạt 4,50 km³/năm. Tuy nhiên, tình trạng gia tăng nhiệt độ cùng sự phân bố mưa cực đoan kết hợp với áp lực phát triển kinh tế xã hội đang khiến chất lượng nước suy giảm nghiêm trọng.

Luận văn tập trung giải quyết vấn đề đánh giá định lượng và dự báo xu thế biến đổi chất lượng nguồn nước mặt trước tác động kép của biến đổi khí hậu và phát triển kinh tế. Mục tiêu trọng tâm là ứng dụng hệ thống mô hình tích hợp GIBSI để mô phỏng chi tiết chế độ thủy văn, quá trình xói mòn và động học chất lượng nước trên diện tích 4.415,5 km² thuộc hai tỉnh Bắc Kạn và Thái Nguyên, từ đó xây dựng các kịch bản dự báo đến năm 2050 và đề xuất giải pháp kiểm soát ô nhiễm.

Phạm vi nghiên cứu bao quát mạng lưới 1.003 đoạn sông và 2 hồ chứa lớn, sử dụng chuỗi dữ liệu thủy văn khí tượng từ năm 1997 đến năm 2012 để hiệu chỉnh và kiểm định mô hình. Kết quả nghiên cứu mang ý nghĩa thực tiễn to lớn khi thiết lập bộ công cụ khoa học có độ tin cậy cao, giúp các nhà quản lý dự báo chính xác nguy cơ suy giảm oxy hòa tan từ 8% đến 16% và kiểm soát tải lượng bùn cát lơ lửng gia tăng trong các kịch bản phát thải khí nhà kính tương lai.

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 thủy văn phân bố kết hợp với động học chất lượng nước và xói mòn lưu vực. Khung lý thuyết chính tích hợp ba mô hình thành phần: mô hình cân bằng nước ba tầng đất BV3C kết hợp sóng động học Saint-Venant một chiều trong HYDROTEL; mô hình mất đất phổ quát hiệu chỉnh RUSLE của Wischmeier và Smith; và mô hình động học chất lượng nước QUAL2E bắt nguồn từ lý thuyết truyền tải cân bằng oxy hòa tan Streeter-Phelps.

Các khái niệm then chốt được áp dụng trong luận văn bao gồm: Đơn vị thủy văn tương đối đồng nhất (RHHU) giúp chia nhỏ lưu vực theo đặc tính địa hình và sử dụng đất; Nhu cầu oxy hóa học (COD) và sinh hóa (BOD5) phản ánh tải lượng ô nhiễm hữu cơ; Oxy hòa tan (DO) biểu thị sức khỏe hệ sinh thái thủy vực; Tổng chất rắn lơ lửng (TSS) đại diện cho mức độ bồi lắng và xói mòn; cùng các kịch bản biến đổi khí hậu B1 (phát thải thấp) và A2 (phát thải trung bình cao) theo tiêu chuẩn của Ban liên Chính phủ về Biến đổi Khí hậu.

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

Nghiên cứu sử dụng nguồn dữ liệu đa tầng bao gồm: bản đồ mô hình số độ cao DEM, dữ liệu viễn thám phân loại 11 nhóm cơ giới đất, diện tích 194.554 ha đất canh tác nông nghiệp, cùng chuỗi số liệu quan trắc khí tượng thủy văn tại 8 trạm khí tượng và các trạm thủy văn Gia Bảy, Chác, Thác Giềng. Lưu vực thượng du được cấu trúc thành 1.015 đơn vị thủy văn RHHU tương đối đồng nhất để phản ánh chính xác tương tác không gian.

Phương pháp mô phỏng số tích hợp trên hệ thống GIBSI được lựa chọn vì khả năng kết nối liền mạch giữa mô đun thủy văn HYDROTEL, mô đun xói mòn RUSLE, mô đun nông nghiệp SWAT và mô đun chất lượng nước QUAL2E. Mô hình được hiệu chỉnh trong giai đoạn 1997 - 2001 bằng phương pháp thử sai, đạt chỉ số hiệu quả Nash-Sutcliffe là 0,82 và sai số biên RCEQM đạt 1,122 m³/s. Quá trình kiểm định giai đoạn 2002 - 2008 và chuỗi đối chuẩn 2009 - 2012 cho chỉ số Nash-Sutcliffe đạt 0,75, khẳng định độ tin cậy vững chắc của bộ thông số trước khi chạy kịch bản dự báo năm 2050.

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

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

Kết quả đánh giá hiện trạng cho thấy sự phân hóa chất lượng nước sâu sắc theo không gian địa lý. Tại vùng thượng nguồn tỉnh Bắc Kạn, chất lượng nước còn trong lành với nồng độ TSS trung bình chỉ từ 27 đến 36 mg/l và DO duy trì ở mức cao từ 5,93 đến 8,89 mg/l. Ngược lại, tại khu vực đô thị thành phố Thái Nguyên, suối Cam Giá và suối Mỏ Bạch bị ô nhiễm hữu cơ nghiêm trọng khi hàm lượng amoni (NH4+) vượt từ 1,3 đến 5,45 lần quy chuẩn kỹ thuật quốc gia, đồng thời BOD5 vượt từ 1,11 đến 1,17 lần ngưỡng cho phép do tiếp nhận trực tiếp nước thải sinh hoạt và luyện kim chưa qua xử lý. Tại hạ lưu trạm Cầu Vát, giá trị TSS trung bình vọt lên mức 109 mg/l do ảnh hưởng cộng hưởng của hoạt động khai thác cát sỏi.

Dưới tác động của biến đổi khí hậu đến năm 2050 theo kịch bản B1 và A2, nhiệt độ không khí trung bình tăng từ 1,2°C đến 2,4°C kéo theo nhiệt độ nước sông tăng tương ứng. Quá trình này làm giảm độ hòa tan của các chất khí, dẫn đến tổng lượng DO suy giảm từ 8% đến 16% trong các tháng mùa khô. Trong khi đó, vào mùa mưa, lượng dòng chảy mặt tăng cao làm gia tăng tải lượng xói mòn đất từ 12% đến 25%, thúc đẩy nồng độ bùn cát lơ lửng và tổng nitơ rửa trôi từ đất nông nghiệp tăng thêm khoảng 15% đến 28% so với năm cơ sở.

Thảo luận kết quả

Diễn biến suy giảm chất lượng nước bắt nguồn từ cơ chế thủy nhiệt và thủy văn học phức tạp. Vào mùa kiệt, khi dòng chảy chỉ chiếm 15% đến 20% tổng lượng nước năm, nhiệt độ nước tăng làm tốc độ phân hủy chất hữu cơ của vi sinh vật hiếu khí tăng nhanh, tiêu hao lượng lớn oxy trong khi khả năng hòa tan oxy từ khí quyển lại suy giảm mạnh theo định luật Henry. Dữ liệu mô phỏng khi được thể hiện qua các biểu đồ tương quan lưu lượng - nồng độ và bản đồ phân vùng ô nhiễm GIS cho thấy rõ nguy cơ thiếu hụt oxy cục bộ tại các đoạn sông tiếp nhận nước thải đô thị.

So sánh với các nghiên cứu lưu vực sông trên thế giới như dự án khảo sát sông ngòi tại châu Âu hay các nghiên cứu tại lưu vực sông Nhuệ - Đáy, kết quả luận văn có sự tương đồng hoàn toàn về quy luật suy giảm tự làm sạch khi nhiệt độ tăng. Điều này chứng minh rằng hiện tượng thời tiết cực đoan kết hợp với nguồn thải nhân tạo sẽ tạo ra hiệu ứng suy thoái kép, đòi hỏi phải kiểm soát đồng thời cả nguồn thải điểm và nguồn thải diện.

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

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 Cầu trong bối cảnh biến đổi khí hậu, luận văn đưa ra bốn nhóm giải pháp hành động cụ thể:

Thứ nhất, xây dựng và đưa vào vận hành hệ thống thu gom, xử lý nước thải sinh hoạt tập trung tại thành phố Thái Nguyên và thị trấn Chợ Mới giai đoạn 2025 - 2030, hướng đến mục tiêu 100% nước thải đô thị đạt quy chuẩn trước khi xả thải, cắt giảm ít nhất 40% tải lượng BOD5 và 50% tải lượng amoni xả vào lưu vực sông. Cơ quan chủ trì thực hiện là Ủy ban nhân dân tỉnh Thái Nguyên và Sở Xây dựng.

Thứ hai, áp dụng các kỹ thuật canh tác nông nghiệp bền vững trên 194.554 ha đất canh tác ven sườn dốc, bao gồm trồng cây theo đường đồng mức và duy trì vành đai thực vật đệm ven sông nhằm giảm thiểu từ 25% đến 30% lượng đất mất do xói mòn và hạn chế rửa trôi phân bón hóa học trong giai đoạn 2025 - 2035. Cơ quan chịu trách nhiệm là Sở Nông nghiệp và Phát triển Nông thôn các tỉnh Bắc Kạn, Thái Nguyên.

Thứ ba, bảo vệ nghiêm ngặt 399.500 ha rừng tự nhiên hiện có và đẩy mạnh trồng rừng phòng hộ đầu nguồn tại các huyện Bạch Thông, Chợ Đồn, Đại Từ, nhằm nâng tỷ lệ che phủ rừng thêm 5% đến 8% trước năm 2030, giúp tăng dung tích trữ nước ngầm mùa khô và cắt giảm lưu lượng đỉnh lũ quét. Cơ quan thực hiện là Chi cục Kiểm lâm và Ban Quản lý rừng phòng hộ liên tỉnh.

Thứ tư, thiết lập mạng lưới quan trắc chất lượng nước tự động liên tỉnh kết nối dữ liệu thời gian thực với mô hình GIBSI phục vụ điều hành xả lũ, điều tiết hồ chứa và cảnh báo sớm ô nhiễm nguồn nước từ nay đến năm 2028. Cơ quan điều phối là Ủy ban Quản lý Lưu vực sông Cầu phối hợp cùng Bộ Tài nguyên và Môi trường.

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

Nội dung luận văn mang lại giá trị học thuật và ứng dụng thực tiễn cao cho bốn nhóm độc giả chuyên biệt:

Nhóm thứ nhất là các nhà quản lý môi trường, quy hoạch tài nguyên nước tại các Sở Tài nguyên và Môi trường tỉnh Bắc Kạn, Thái Nguyên và các tỉnh lân cận. Nghiên cứu cung cấp cơ sở dữ liệu số hóa và bộ công cụ định lượng hỗ trợ phê duyệt báo cáo đánh giá tác động môi trường, xác định hạn ngạch xả nước thải và xây dựng kế hoạch hành động ứng phó biến đổi khí hậu cấp địa phương.

Nhóm thứ hai là các nhà khoa học, nghiên cứu sinh và học viên cao học chuyên ngành Khoa học Môi trường, Quản lý Tài nguyên, Thủy văn học. Công trình là tài liệu tham khảo mẫu mực về quy trình thiết lập, hiệu chỉnh và kiểm định hệ thống mô hình phân bố phức hợp HYDROTEL - RUSLE - SWAT - QUAL2E trên nền tảng hệ thông tin địa lý.

Nhóm thứ ba là các kỹ sư tư vấn, chuyên gia tại các viện nghiên cứu và tổ chức phát triển môi trường. Tài liệu cung cấp các thông số thực nghiệm, hướng dẫn xây dựng giải pháp công trình và phi công trình để thiết kế các dự án phục hồi hệ sinh thái nước ngọt thích ứng thời tiết cực đoan.

Nhóm thứ tư là giảng viên và sinh viên tại các trường đại học khối ngành tài nguyên, môi trường và địa lý học. Luận văn đóng vai trò là học liệu chuyên sâu phục vụ giảng dạy môn học Mô hình hóa môi trường và Quản lý tổng hợp tài nguyên nước lưu vực sông.

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

Mô hình GIBSI có điểm gì ưu việt hơn so với các mô hình thủy văn truyền thống? GIBSI là hệ thống phần mềm tích hợp hoàn chỉnh liên kết trực tiếp bốn mô hình HYDROTEL, RUSLE, SWAT và QUAL2E trên nền tảng hệ thống thông tin địa lý GIS. Điểm vượt trội là khả năng mô phỏng đồng thời dòng chảy phân bố, tương tác cân bằng nước 3 tầng đất, xói mòn sườn dốc và động học chất lượng nước theo bước thời gian ngày cho từng tiểu lưu vực riêng biệt.

Tại sao việc phân chia lưu vực thành các đơn vị thủy văn RHHU lại giữ vai trò cốt lõi? Đơn vị thủy văn tương đối đồng nhất RHHU cho phép mô hình hóa chi tiết 1.015 tiểu vùng có cùng tính chất thổ nhưỡng, độ dốc và thảm phủ thực vật. Cách tiếp cận này giúp tính toán chính xác lượng mưa ngấm, dòng chảy mặt và lượng phân bón rửa trôi từ 194.554 ha đất nông nghiệp, hạn chế tối đa sai số gộp của các mô hình tham số tập trung.

Biến đổi khí hậu làm suy giảm nồng độ oxy hòa tan DO trong nước sông qua cơ chế nào? Theo định luật Henry, độ tan của oxy vào nước tỷ lệ nghịch với nhiệt độ. Khi nhiệt độ nước sông Cầu tăng từ 1,2°C đến 2,4°C vào năm 2050, khả năng hấp thụ oxy tự nhiên suy giảm mạnh, đồng thời nhiệt độ cao kích thích hoạt tính của vi sinh vật hiếu khí làm tăng tốc độ tiêu thụ oxy, dẫn đến DO giảm từ 8% đến 16%.

Nguyên nhân chính gây ô nhiễm hữu cơ và amoni vượt chuẩn tại thành phố Thái Nguyên là gì? Tình trạng ô nhiễm tại suối Cam Giá và Mỏ Bạch xuất phát từ việc tiếp nhận lưu lượng lớn nước thải sinh hoạt đô thị và nước thải công nghiệp luyện kim chưa qua xử lý chuẩn mực. Khi đổ vào dòng chính sông Cầu vào mùa khô với dòng chảy suy giảm, nồng độ NH4+ bị tích tụ và vượt từ 1,3 đến 5,45 lần quy chuẩn cột B2.

Quy trình ứng dụng mô hình trong luận văn có thể chuyển giao cho các lưu vực khác không? Phương pháp luận và quy trình tích hợp mô hình GIBSI hoàn toàn có thể nhân rộng cho các lưu vực sông khác tại Việt Nam. Khi áp dụng cho lưu vực mới, chỉ cần thu thập bản đồ số độ cao DEM, dữ liệu sử dụng đất, loại đất cùng chuỗi khí tượng thủy văn địa phương để tiến hành hiệu chỉnh lại 7 thông số thủy lực cơ bản.

Kết luận

  • Ứng dụng thành công hệ sinh thái mô hình phân bố GIBSI để mô phỏng chính xác chế độ thủy văn và động học chất lượng nước trên 4.415,5 km² lưu vực thượng du sông Cầu với chỉ số Nash-Sutcliffe đạt mức tốt 0,82.
  • Định vị rõ nét các điểm nóng ô nhiễm hữu cơ, suy giảm DO và dư lượng amoni cao tập trung tại khu vực đô thị, công nghiệp thành phố Thái Nguyên và vùng nông nghiệp dốc.
  • Định lượng hóa xu thế tác động của biến đổi khí hậu đến năm 2050, cảnh báo nguy cơ suy giảm từ 8% đến 16% oxy hòa tan mùa kiệt và gia tăng từ 12% đến 25% bùn cát xói mòn mùa mưa.
  • Đóng góp bộ cơ sở dữ liệu số hóa 1.015 tiểu lưu vực RHHU và hệ số thực nghiệm tin cậy phục vụ công tác quy hoạch phân vùng xả thải môi trường nước.
  • Đề xuất khung giải pháp toàn diện kết hợp giữa xử lý nước thải đô thị, canh tác chống xói mòn, phục hồi rừng phòng hộ và thiết lập hệ thống quan trắc tự động liên tỉnh trước năm 2030.

Nghiên cứu mở ra bước tiến quan trọng trong việc chuyển đổi từ quản lý tài nguyên nước bị động sang quản trị chủ động dựa trên mô hình toán và kịch bản dự báo. Các cơ quan ban ngành và chính quyền địa phương cần khẩn trương đưa bộ chỉ số mô hình vào quy hoạch phân vùng môi trường lưu vực nhằm đảm bảo an ninh nguồn nước bền vững cho tương lai.