Tổng quan nghiên cứu

Biến đổi khí hậu đang tạo ra những sức ép nặng nề lên an ninh nguồn nước toàn cầu và khu vực. Tại Việt Nam, các ghi nhận lịch sử trong 50 năm qua cho thấy nhiệt độ trung bình năm đã tăng từ 0,5 đến 0,7°C, trong khi mực nước biển dâng khoảng 20 cm, làm gia tăng tần suất của các hiện tượng thời tiết cực đoan. Lưu vực sông Lô đóng vai trò huyết mạch đối với vùng kinh tế trọng điểm Bắc Bộ với tổng diện tích lưu vực đạt 37.878 km², trong đó phần diện tích thuộc lãnh thổ Việt Nam là 22.629 km² (chiếm 6,84% diện tích cả nước) trải dài trên địa bàn của 8 tỉnh miền núi và trung du. Đây là nơi sinh sống của khoảng 2,94 triệu người với 78,5% lực lượng lao động phụ thuộc trực tiếp vào ngành nông nghiệp – một lĩnh vực cực kỳ nhạy cảm với các biến động thủy văn.

Nghiên cứu tập trung giải quyết bài toán cốt lõi: Đánh giá định lượng mức độ tác động của biến đổi khí hậu đến các đặc trưng tài nguyên nước mặt trên lưu vực sông Lô theo chuỗi không gian và thời gian. Mục tiêu cụ thể là xác định xu thế biến động của dòng chảy đến, nhu cầu dùng nước, cân bằng nước hệ thống, đỉnh lũ thiết kế và mức độ hạn hán theo các kịch bản phát thải khí nhà kính khác nhau (A2, B2, B1). Dựa trên chuỗi số liệu quan trắc thực đo giai đoạn 1970–2010 và các mô hình dự tính khí hậu tương lai đến cuối thế kỷ 21, công trình cung cấp cơ sở khoa học định lượng chuẩn xác, hỗ trợ đắc lực cho công tác quy hoạch thủy lợi, phòng chống thiên tai và quản lý bền vững nguồn nước mặt cho toàn bộ khu vực Đông Bắc Bộ.

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 chu trình thủy văn, động lực học dòng chảy trong hệ thống sông ngòi và các lý thuyết về kịch bản phát thải khí nhà kính toàn cầu. Trọng tâm lý thuyết xoay quanh báo cáo đặc biệt về kịch bản phát thải (SRES) của Ban Liên chính phủ về Biến đổi Khí hậu (IPCC) với ba kịch bản phát thải tiêu biểu: A2 (phát thải cao, mô phỏng thế giới phát triển không đồng nhất), B2 (phát thải trung bình, chú trọng tính bền vững theo vùng) và B1 (phát thải thấp, hướng tới nền kinh tế sạch và sử dụng tài nguyên hiệu quả).

Để mô phỏng động lực dòng chảy hở trong mạng lưới sông tự nhiên, nghiên cứu ứng dụng hệ phương trình Saint-Venant 1 chiều gồm phương trình liên tục và phương trình chuyển động, kết hợp công thức Manning để tính toán sức cản lòng dẫn. Ngoài ra, lý thuyết cân bằng nước hệ thống được tích hợp nhằm mô hình hóa tương tác phức tạp giữa nguồn nước đến tự nhiên, nhu cầu khai thác sử dụng nước của các ngành kinh tế - xã hội và khả năng điều tiết của các công trình hồ chứa lớn trên lưu vực.

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

Luận văn kết hợp đồng bộ các phương pháp thống kê, viễn thám GIS và chuỗi mô hình toán chuyên ngành. Cỡ mẫu dữ liệu đầu vào bao gồm chuỗi số liệu khí tượng quan trắc đồng bộ giai đoạn 1960–2012 từ 16 trạm khí tượng tiêu biểu (Hà Giang, Bắc Quang, Chiêm Hóa, Phú Hộ, Việt Trì...) và số liệu thủy văn liên tục của 7 trạm chính trên hệ thống sông Lô, sông Gâm, sông Chảy. Phương pháp chọn mẫu theo mạng lưới trạm đại diện phân bố đều từ thượng nguồn đến hạ lưu, bao phủ đầy đủ các tiểu vùng địa hình từ vùng núi cao trên 2.000 m đến vùng đồng bằng đệ tam. Chuỗi số liệu từ năm 1980 đến năm 1990 được dùng để hiệu chỉnh bộ thông số, và chuỗi 1991–2000 được sử dụng để kiểm định mô hình.

Bộ công cụ phân tích bao gồm:

  • Mô hình thủy văn NAM: Mô phỏng quá trình mưa - dòng chảy thông qua cấu trúc 5 bể chứa ẩm thẳng đứng (bể chứa tuyết, bể chứa bề mặt, bể chứa rễ cây, bể ngầm tầng trên và bể ngầm tầng dưới). Toàn bộ lưu vực được chia thành 7 tiểu lưu vực tính toán.
  • Mô hình thủy lực MIKE 11: Ứng dụng sơ đồ sai phân hữu hạn 6 điểm ẩn Abbott để giải hệ phương trình Saint-Venant, mô phỏng chế độ truyền lũ và diễn biến mực nước trên toàn mạng lưới sông.
  • Mô hình MIKE BASIN: Phân tích cân bằng nước và phân bổ nguồn nước cho 7 tiểu vùng sử dụng nước theo từng tháng.
  • Mô hình PRECIS và công cụ chi tiết hóa thống kê (SDSM, SIMCLIM): Xây dựng các kịch bản khí hậu chi tiết với độ phân giải lưới 25x25 km.

Lý do lựa chọn chuỗi mô hình này xuất phát từ khả năng tái hiện chân thực các đặc tính phi tuyến của dòng chảy lưu vực miền núi có địa hình phân cắt mạnh và hệ thống hồ chứa điều tiết phức tạp. Kết quả kiểm định mô hình NAM cho thấy chỉ số hiệu quả Nash-Sutcliffe (NASH) tại tất cả các trạm chính đều đạt độ tin cậy rất cao, dao động từ 0,70 đến 0,93 (trạm Vụ Quang đạt NASH hiệu chỉnh 0,90 và kiểm định 0,85; trạm Hàm Yên đạt NASH kiểm định 0,93).

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

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

Phân tích chuỗi số liệu thực đo giai đoạn 1970–2010 đã chỉ rõ xu thế nhiệt độ trung bình năm trên lưu vực sông Lô tăng khoảng 0,95°C, trong khi lượng mưa năm có xu thế giảm khoảng 18,46%. Phân bố lượng mưa phân hóa sâu sắc khi 80–85% tổng lượng mưa năm chỉ tập trung trong 6 tháng mùa mưa, khiến dòng chảy kiệt suy giảm rõ rệt.

Dự tính theo các kịch bản biến đổi khí hậu đến cuối thế kỷ 21 (giai đoạn 2080–2099) so với thời kỳ nền (1980–1999) cho thấy:

  • Về nhiệt độ: Nhiệt độ trung bình năm tăng 1,7°C theo kịch bản B1; 2,3°C theo kịch bản B2 và tăng tới 2,8°C theo kịch bản A2. Mức tăng nhiệt độ vào mùa mưa cao hơn mùa khô khoảng 0,5°C trên hầu khắp các trạm quan trắc.
  • Về lượng mưa: Tổng lượng mưa năm có xu hướng tăng từ 2,9% (kịch bản B1), 4,1% (kịch bản B2) đến 4,9% (kịch bản A2). Tuy nhiên, lượng mưa mùa khô giảm nghiêm trọng, đặc biệt trong tháng 1 mức giảm lên đến 15,8% (B1), 17,5% (B2) và 18,7% (A2). Ngược lại, lượng mưa các tháng đỉnh mùa lũ (tháng 7 và 8) tăng mạnh từ 6% đến trên 10%.
  • Về chế độ dòng chảy và lũ: Dòng chảy trung bình năm tại trạm Vụ Quang có xu thế tăng nhẹ, nhưng dòng chảy mùa lũ tăng từ 5% đến 12%, đẩy lưu lượng đỉnh lũ thiết kế tần suất 1% và 5% tăng đáng kể. Mực nước lớn nhất tại trạm Vụ Quang tăng từ 0,35 m đến 0,78 m dưới tác động của kịch bản A2.
  • Về hạn hán và cân bằng nước: Chỉ số hạn hán ($K_{hạn}$) tính toán trong các tháng mùa kiệt (tháng 11 đến tháng 3) có xu hướng gia tăng mức độ khốc liệt. Lượng nước thiếu hụt trên lưu vực mở rộng mạnh ở các khu vực hạ du, gây áp lực thiếu nước nghiêm trọng cho sản xuất nông nghiệp từ 15% đến 28% so với nhu cầu thực tế.

Thảo luận kết quả

Sự suy giảm nghiêm trọng của dòng chảy mùa khô kết hợp với mức bốc hơi tiềm năng tăng cao (trung bình toàn lưu vực hiện tại đạt 765 mm/năm, tại một số nơi như Phú Thọ và Việt Trì đạt tới 977 mm/năm) là nguyên nhân chính dẫn đến tình trạng mất cân bằng nước nghiêm trọng vào các tháng đầu năm. Dữ liệu mô phỏng khi được thể hiện qua các đồ thị thủy văn mùa và biểu đồ đường quá trình lưu lượng cho thấy sự tương phản gay gắt giữa hai thái cực lũ và kiệt. Đỉnh lũ chuyển dịch cao hơn và dốc hơn, trong khi đáy kiệt bị kéo dài và hạ thấp sâu hơn.

Bản đồ phân vùng hạn hán trên lưu vực sông Lô qua các giai đoạn 2050 và 2080 thể hiện rõ diện tích vùng bị khô hạn cấp độ nặng lan rộng tại các khu vực thung lũng sông Gâm và vùng đệm trung lưu. So sánh với các báo cáo đánh giá của Bộ Tài nguyên và Môi trường cũng như các nghiên cứu trên hệ thống lưu vực sông Hồng - Thái Bình, kết quả nghiên cứu này hoàn toàn đồng nhất về xu hướng cực đoan hóa của tài nguyên nước. Điểm khác biệt mang tính đóng góp cao của luận văn là đã định lượng chi tiết tương tác thủy lực mạng sông có xét đến dung tích điều tiết của hai hồ chứa lớn là Thác Bà và Tuyên Quang, từ đó phản ánh sát thực tế khả năng cấp nước và cắt lũ của hệ thống công trình hiện hữu.

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

Dựa trên các kết quả tính toán định lượng về diễn biến tài nguyên nước, luận văn đề xuất 4 nhóm giải pháp chiến lược:

  • Tối ưu hóa quy trình vận hành liên hồ chứa: Bộ Nông nghiệp và Phát triển nông thôn cùng Tập đoàn Điện lực Việt Nam cần phối hợp rà soát, điều chỉnh quy trình vận hành liên hồ Tuyên Quang và Thác Bà trước năm 2028. Cần nâng cao dung tích phòng lũ chủ động thêm 10–15% vào các tháng cao điểm lũ (tháng 7–8) và thực hiện tích nước sớm hơn vào cuối mùa lũ để bảo đảm xả dòng chảy duy trì môi trường tối thiểu không dưới 45 m³/s trong các tháng kiệt nhất (tháng 1–3).
  • Hiện đại hóa hệ thống quan trắc và cảnh báo sớm: Ủy ban nhân dân các tỉnh thuộc lưu vực phối hợp với Trung tâm Dự báo Khí tượng Thủy văn Quốc gia triển khai lắp đặt bổ sung 15 trạm đo mưa tự động và 4 trạm quan trắc mực nước tự động tại các vị trí đầu nguồn xung yếu (như Bắc Quang, Hoàng Su Phì, thượng lưu sông Con và sông Miện) giai đoạn 2025–2030, nhằm tăng thời gian dự báo lũ sớm lên tối thiểu 6 đến 12 giờ.
  • Chuyển đổi cơ cấu kinh tế nông nghiệp và tưới tiết kiệm: Sở Nông nghiệp và Phát triển nông thôn các tỉnh Tuyên Quang, Hà Giang, Phú Thọ chỉ đạo chuyển đổi ít nhất 20% diện tích đất lúa nước có nguy cơ thiếu nước cao sang cây trồng cạn có giá trị kinh tế cao, đồng thời áp dụng công nghệ tưới tiết kiệm nước nhằm cắt giảm 15% tổng nhu cầu nước tưới nông nghiệp vào năm 2035.
  • Củng cố hạ tầng đê điều và thoát lũ vùng hạ du: Cục Quản lý Đê điều và Phòng chống Thiên tai tiến hành nâng cấp, gia cố các tuyến đê xung yếu tại Phú Thọ, Vĩnh Phúc đáp ứng mực nước đỉnh lũ thiết kế tăng thêm 0,5–0,8 m theo kịch bản biến đổi khí hậu A2, hoàn thành trước năm 2040.

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

  • Cơ quan quản lý nhà nước về tài nguyên nước và khí tượng thủy văn: Cục Quản lý Tài nguyên nước, Viện Khoa học Khí tượng Thủy văn và Biến đổi khí hậu, Sở Tài nguyên và Môi trường 8 tỉnh thuộc lưu vực sử dụng số liệu làm căn cứ lập quy hoạch phân bổ tài nguyên nước và kế hoạch hành động ứng phó biến đổi khí hậu cấp tỉnh.
  • Các đơn vị quản lý, vận hành hồ chứa và thủy điện: Ban quản lý hồ thủy điện Tuyên Quang, thủy điện Thác Bà và các công trình thủy lợi khai thác thông tin về lưu lượng dòng chảy đến để xây dựng phương án tích nước, xả lũ và phát điện tối ưu.
  • Giảng viên, nghiên cứu sinh và sinh viên ngành Thủy văn học, Kỹ thuật tài nguyên nước: Sử dụng tài liệu như một nghiên cứu điển hình về phương pháp luận kết hợp chuỗi mô hình NAM - MIKE 11 - MIKE BASIN và kỹ thuật downscaling khí hậu PRECIS.
  • Các tổ chức phát triển và chuyên gia quy hoạch phát triển vùng: Tham khảo các phân tích không gian về cân bằng nước và mức độ dễ bị tổn thương để lập kế hoạch đầu tư các dự án thích ứng biến đổi khí hậu dựa vào hệ sinh thái lưu vực sông.

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

Mô hình thủy văn NAM có độ chính xác như thế nào khi áp dụng cho lưu vực sông Lô?
Mô hình NAM đạt độ tin cậy rất cao trên lưu vực sông Lô với hệ số hiệu quả Nash-Sutcliffe (NASH) dao động từ 0,74 đến 0,90 trong giai đoạn hiệu chỉnh (1980–1990) và đạt từ 0,70 đến 0,93 trong giai đoạn kiểm định (1991–2000) tại 6 trạm thủy văn chính, đáp ứng đầy đủ tiêu chuẩn khoa học để mô phỏng dòng chảy theo các kịch bản khí hậu tương lai.

Biến đổi khí hậu làm suy giảm dòng chảy mùa kiệt lưu vực sông Lô ở mức độ nào?
Dưới tác động của biến đổi khí hậu, lượng mưa các tháng mùa kiệt có xu thế giảm sâu, đặc biệt trong tháng 1 giảm từ 15,8% đến 18,7%. Điều này khiến dòng chảy kiệt các tháng 1, 2, 3 suy giảm từ 8% đến 14%, đẩy chỉ số hạn hán $K_{hạn}$ lên mức báo động và làm thiếu hụt nguồn nước tưới nghiêm trọng tại vùng trung và hạ lưu.

Kịch bản phát thải nào gây ra tác động tiêu cực nhất đến tài nguyên nước?
Kịch bản phát thải cao A2 gây ra tác động bất lợi lớn nhất khi đẩy nhiệt độ trung bình năm tăng tới 2,8°C vào cuối thế kỷ 21, làm tăng mạnh bốc hơi tiềm năng, khiến chênh lệch dòng chảy giữa mùa lũ và mùa kiệt trở nên sâu sắc nhất, đồng thời đẩy mực nước đỉnh lũ tại trạm Vụ Quang tăng thêm khoảng 0,78 m.

Tại sao lượng mưa năm dự báo tăng nhưng nguy cơ hạn hán cục bộ vẫn trầm trọng hơn?
Mặc dù tổng lượng mưa năm tăng từ 2,9% đến 4,9%, lượng mưa gia tăng chủ yếu dồn vào các tháng cao điểm mùa lũ (tháng 7 và 8), trong khi các tháng mùa khô lại giảm sâu. Đi đôi với nhiệt độ tăng cao làm tăng lượng bốc thoát hơi nước của cây trồng, sự phân bố mưa lệch mùa đã làm gia tăng tần suất và cường độ hạn hán mùa khô.

Hồ chứa Thác Bà và Tuyên Quang giữ vai trò gì trong việc ứng phó với biến đổi khí hậu?
Hai hồ chứa lớn này đóng vai trò then chốt trong việc cắt giảm đỉnh lũ mùa mưa cho vùng hạ du và điều tiết bổ sung dòng chảy kiệt trong các tháng đầu năm. Theo mô hình MIKE BASIN, việc vận hành hợp lý dung tích hữu ích của hai hồ giúp bù đắp từ 30% đến 45% lượng nước thiếu hụt cho các tiểu vùng sử dụng nước hạ lưu.

Kết luận

  • Luận văn đã định lượng hóa thành công tác động của biến đổi khí hậu đến tài nguyên nước lưu vực sông Lô, chỉ ra mức tăng nhiệt độ từ 1,7°C đến 2,8°C và lượng mưa năm tăng từ 2,9% đến 4,9% theo các kịch bản B1, B2, A2 vào giai đoạn 2080–2099.
  • Dòng chảy mặt có sự phân hóa mùa sâu sắc: dòng chảy mùa lũ tăng từ 5% đến 12% làm gia tăng áp lực ngập lụt, trong khi dòng chảy mùa kiệt suy giảm mạnh, đẩy cao nguy cơ hạn hán và thiếu hụt nước tại hạ lưu.
  • Phương pháp luận tích hợp chuỗi mô hình tiên tiến NAM, MIKE 11 và MIKE BASIN đã chứng minh tính ưu việt với độ chính xác cao (NASH > 0,70–0,93), tạo khung công cụ tiêu chuẩn cho các nghiên cứu thủy văn lưu vực lớn.
  • Các kết quả tính toán cân bằng nước cung cấp luận cứ thực tiễn trực tiếp cho việc điều chỉnh quy trình vận hành liên hồ chứa Thác Bà - Tuyên Quang và định hướng chuyển dịch cơ cấu mùa vụ nông nghiệp vùng Đông Bắc.
  • Để ứng phó hiệu quả, các cơ quan ban ngành cần khẩn trương hoàn thiện thể chế quản lý tổng hợp tài nguyên nước theo lưu vực sông, tăng cường đầu tư mạng lưới quan trắc tự động và triển khai các công trình phòng chống thiên tai thích ứng linh hoạt với các kịch bản khí hậu trong tương lai.