Tổng quan nghiên cứu

Biến đổi khí hậu đang trở thành thách thức môi trường cấp bách toàn cầu trong thế kỷ 21, làm biến động sâu sắc chu trình thủy văn và gia tăng tần suất thiên tai cực đoan. Theo ước tính của Bộ Tài nguyên và Môi trường, trong vòng 50 năm qua, nhiệt độ trung bình tại Việt Nam đã tăng từ 0,5°C đến 0,7°C và mực nước biển dâng khoảng 20 cm. Các dự báo khí hậu chỉ ra rằng đến năm 2100, nhiệt độ trung bình toàn quốc có thể tăng thêm 2°C đến 3°C, dẫn tới nguy cơ mất an ninh nguồn nước nghiêm trọng.

Lưu vực sông Cầu là một trong những hệ thống sông huyết mạch của vùng kinh tế trọng điểm phía Bắc, có diện tích tự nhiên 6.030 km² cùng mạng lưới sông ngòi phân nhánh với tổng chiều dài khoảng 1.600 km. Dòng chảy sông Cầu giữ vai trò then chốt trong việc cung cấp nguồn nước sinh hoạt, sản xuất công nghiệp và nông nghiệp cho hàng triệu người dân trải rộng trên địa bàn 6 tỉnh thành gồm Bắc Kạn, Thái Nguyên, Vĩnh Phúc, Bắc Giang, Bắc Ninh và Hà Nội. Tuy nhiên, tình trạng biến đổi thời tiết và các dao động khí hậu bất thường đã khiến chế độ thủy văn trên lưu vực biến động phức tạp, gây ra lũ quét dữ dội vào mùa mưa và hạn hán, cạn kiệt nguồn nước nghiêm trọng trong mùa khô.

Mục tiêu cốt lõi của nghiên cứu là đánh giá định tính và định lượng sự biến đổi của các yếu tố khí hậu cơ bản như nhiệt độ, lượng mưa, đồng thời làm rõ mức độ tác động của biến đổi khí hậu đến các đặc trưng dòng chảy mặt (lưu lượng dòng chảy năm, dòng chảy lũ, dòng chảy kiệt và mực nước đỉnh lũ). Nghiên cứu được triển khai trên phạm vi toàn bộ lưu vực sông Cầu, phân tích chuỗi dữ liệu thực đo lịch sử giai đoạn 1961–2013 và tính toán dự phóng dòng chảy tương lai giai đoạn 2014–2033 theo kịch bản phát thải trung bình A1B. Kết quả luận văn cung cấp cơ sở khoa học và hệ số liệu định lượng tin cậy, giúp nâng cao năng lực ứng phó thiên tai và quản lý bền vững tài nguyên nước lưu vực.

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 cân bằng nước lưu vực sông và phương pháp tiếp cận hệ thống động lực tích hợp tự nhiên - kinh tế - xã hội. Theo định nghĩa chuẩn hóa của Ủy ban Liên chính phủ về Biến đổi Khí hậu (IPCC), biến đổi khí hậu là sự biến đổi trạng thái của hệ thống khí hậu được duy trì trong một khoảng thời gian đủ dài, điển hình là qua nhiều thập kỷ. Các khái niệm học thuật chủ chốt được vận dụng xuyên suốt gồm: kịch bản phát thải khí nhà kính (SRES), dao động chuẩn sai khí hậu, lưu lượng dòng chảy mặt ($Q$), module dòng chảy ($M$) và tổng lượng dòng chảy ($W$).

Khung phân tích kế thừa trực tiếp các định hướng từ Chiến lược quốc gia về biến đổi khí hậu và Luật Tài nguyên nước số 17/2012/QH13. Mối quan hệ giữa sự gia tăng nồng độ khí nhà kính và dao động thủy văn được phân tích qua 4 họ kịch bản phát thải gốc gồm A1, A2, B1 và B2. Trong đó, kịch bản trung bình A1B được chọn làm trục tham chiếu chính nhằm mô phỏng trạng thái cân bằng giữa các nguồn năng lượng hóa thạch và năng lượng tái tạo trong giai đoạn phát triển kinh tế tương lai.

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

Nguồn dữ liệu của luận văn được thu thập đồng bộ từ chuỗi quan trắc khí tượng thủy văn 40 năm liên tục (1974–2013) tại 8 trạm khí tượng tiêu biểu (Bắc Kạn, Định Hóa, Thái Nguyên, Tam Đảo, Vĩnh Yên, Bắc Ninh...) và các trạm thủy văn chủ lực như Gia Bảy, Chợ Mới, Thác Bưởi. Cỡ mẫu số liệu gồm hơn 14.600 ngày quan trắc kết hợp cùng số liệu lưới khí hậu tái phân tích toàn cầu. Phương pháp chọn mẫu là chọn mẫu cụm có chủ đích toàn bộ các trạm chuẩn quốc gia đại diện cho cả 3 vùng địa hình: thượng lưu miền núi, trung lưu đồi bán sơn địa và hạ lưu đồng bằng phù sa.

Quy trình xử lý và phân tích số liệu sử dụng 3 phương pháp trụ cột:

  1. Phương pháp thống kê toán học: Ứng dụng phương trình xu thế tuyến tính $y = a_1x + b$ để xác định hệ số góc $a_1$ phản ánh tốc độ biến động hàng năm của nhiệt độ trung bình, lượng mưa và lưu lượng dòng chảy qua hai thời đoạn đối chiếu 1974–1993 và 1994–2013.
  2. Phương pháp công nghệ thông tin và GIS: Sử dụng phần mềm mã nguồn mở OpenGrADS để xử lý tệp dữ liệu khí hậu dạng lưới và hệ thống ArcGIS để số hóa bản đồ phân vùng lưu vực sông Cầu thành các tiểu lưu vực phụ thuộc.
  3. Phương pháp mô hình hóa thủy văn: Vận dụng mô hình dòng chảy mưa NAM (thuộc bộ mô hình MIKE 11) được chuẩn hóa theo tiêu chuẩn của Tổ chức Khí tượng Thế giới (WMO). Lý do lựa chọn mô hình NAM là nhờ cấu trúc mô hình bán phân bố dạng bể chứa tối ưu, có khả năng mô phỏng liên tục dòng chảy mặt, dòng chảy ngầm và độ ẩm đất với số lượng thông số hiệu chỉnh hợp lý, đem lại hệ số tương quan Nash-Sutcliffe đạt trên 0,85 trong quá trình kiểm định thực tế.

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

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

Nghiên cứu đã chỉ ra những thay đổi rõ nét của khí hậu và thủy văn trên lưu vực sông Cầu qua chuỗi số liệu thực tế:

  • Xu thế gia tăng nhiệt độ không khí: Nhiệt độ trung bình năm trên toàn lưu vực tăng với tốc độ từ 0,2°C đến 0,3°C/thập kỷ trong giai đoạn 1974–2013. Đặc biệt, nhiệt độ mùa đông (tháng 1 và tháng 2) có mức tăng mạnh nhất với hệ số góc $a_1$ đạt khoảng +0,035°C/năm. Nhiệt độ cực tiểu ($T_m$) tại trạm Bắc Kạn và Thái Nguyên tăng trung bình 0,5°C/thập kỷ, dẫn đến số ngày rét đậm và rét hại giảm từ 3 đến 5 ngày/năm so với thời kỳ 1961–1990.
  • Biến động phân bố lượng mưa: Lượng mưa trung bình toàn lưu vực đạt 1.700 mm/năm nhưng phân bố không đồng đều theo không gian và thời gian. Vùng mưa lớn nhất tập trung tại Tam Đảo với 2.495 mm/năm và Thái Nguyên đạt 2.074 mm/năm, trong khi thượng nguồn Bắc Kạn chỉ đạt 1.533 mm/năm. Trong thời đoạn 1994–2013, số ngày có mưa lớn dị thường ($R > 50\text{ mm/ngày}$) tăng từ 12% đến 18% tại các trạm trung lưu và thượng lưu.
  • Tác động đến dòng chảy mùa lũ: Dòng chảy mùa lũ (từ tháng 6 đến tháng 10) chiếm từ 75% đến 82% tổng lượng dòng chảy năm. Theo kịch bản A1B, lưu lượng dòng chảy mùa lũ giai đoạn 2014–2033 dự kiến tăng từ 4,2% đến 8,6% so với thời kỳ nền. Mực nước đỉnh lũ lớn nhất ($h_{max}$) tại trạm Gia Bảy có xu thế nâng cao với tốc độ tăng từ 0,8 cm đến 1,5 cm/năm, làm gia tăng nguy cơ ngập úng đô thị tại thành phố Thái Nguyên.
  • Suy giảm dòng chảy mùa cạn: Ngược lại với mùa lũ, dòng chảy mùa kiệt (tháng 12 đến tháng 4 năm sau) chỉ chiếm từ 18% đến 25% lượng nước cả năm và đang có xu thế suy giảm rõ rệt. Lưu lượng dòng chảy mùa cạn trong giai đoạn 1994–2013 giảm khoảng 5,3% so với giai đoạn 1974–1993. Dự tính đến năm 2033, lưu lượng dòng chảy 3 tháng cạn nhất tại các tiểu lưu vực thượng nguồn có thể tiếp tục sụt giảm từ 3,8% đến 7,5%.

Thảo luận kết quả

Sự biến động tương phản giữa dòng chảy mùa lũ và mùa cạn phản ánh rõ tác động kép của hiện tượng nóng lên toàn cầu và suy giảm khả năng điều tiết tự nhiên của thảm phủ lưu vực. Việc gia tăng nhiệt độ bề mặt thúc đẩy cường độ bốc thoát hơi nước tiềm năng, kết hợp với các đợt mưa rào tập trung cường độ lớn làm cho quá trình hình thành dòng chảy mặt diễn ra nhanh hơn. Địa hình thượng lưu sông Cầu có độ dốc đáy sông trên 10‰ cùng lòng sông hẹp, khiến lũ tập trung nhanh, biên độ dao động mực nước lớn và thời gian truyền lũ bị rút ngắn từ 2 đến 4 giờ.

Về phương thức biểu diễn số liệu, các phát hiện nghiên cứu được hệ thống hóa trực quan thông qua bảng đối sánh lưu lượng trung bình tháng và đồ thị đường cong phân bố tích lũy thời gian dòng chảy. Việc thể hiện các ma trận số liệu lưu lượng ($Q$), module dòng chảy ($M$) và tổng lượng nước ($W$) dưới dạng biểu đồ thanh so sánh giữa thời đoạn 1974–1993 và 1994–2013 giúp nhận diện trực tiếp chu kỳ suy giảm của dòng chảy kiệt.

So sánh với các nghiên cứu thủy văn trên hệ thống sông Hồng - Thái Bình của Viện Khoa học Khí tượng Thủy văn và Môi trường, quy luật gia tăng dòng chảy mùa lũ từ 3% đến 10% và giảm dòng chảy mùa cạn từ 4% đến 12% hoàn toàn nhất quán. Tuy nhiên, lưu vực sông Cầu có tính chất cực đoan cục bộ cao hơn do sự suy giảm diện tích rừng tự nhiên tại vùng đệm và hoạt động khai thác cát, khoáng sản lòng sông làm thay đổi hình thái lòng dẫn.

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

Nhằm chủ động ứng phó với diễn biến biến đổi khí hậu và bảo vệ an ninh tài nguyên nước mặt trên lưu vực sông Cầu, luận văn đề xuất 4 nhóm giải pháp chiến lược:

  1. Hiện đại hóa mạng lưới quan trắc khí tượng thủy văn tự động: Sở Tài nguyên và Môi trường các tỉnh Bắc Kạn, Thái Nguyên, Bắc Giang chủ trì lắp đặt bổ sung 15 trạm đo mưa tự động và 6 trạm đo mực nước tự động tại các nhánh sông chính trước năm 2026, nâng cao tỷ lệ dự báo sớm lũ quét và ngập lụt chính xác đạt trên 90%.
  2. Xây dựng quy trình vận hành điều tiết liên hồ chứa linh hoạt: Bộ Nông nghiệp và Phát triển nông thôn phối hợp với Ủy ban nhân dân các tỉnh lưu vực hoàn thiện quy chế xả nước liên hồ chứa thủy lợi trước năm 2027, duy trì xả dòng chảy môi trường tối thiểu trong các tháng mùa kiệt nhằm bù đắp ít nhất 15% lượng nước thiếu hụt cho vùng hạ du.
  3. Khôi phục và nâng cao chất lượng rừng phòng hộ đầu nguồn: Ủy ban nhân dân tỉnh Bắc Kạn và Thái Nguyên triển khai chương trình trồng mới và làm giàu 25.000 ha rừng tự nhiên giai đoạn 2025–2030, phấn đấu nâng tỷ lệ che phủ rừng đầu nguồn lên trên 65% để tăng dung tích trữ nước tự nhiên và giảm hệ số dòng chảy lũ từ 5% đến 8%.
  4. Chuyển đổi cơ cấu cây trồng và quy hoạch sử dụng nước tiết kiệm: Sở Nông nghiệp và Phát triển nông thôn các tỉnh Bắc Ninh, Bắc Giang chuyển đổi tối thiểu 20% diện tích canh tác lúa kém hiệu quả sang cây trồng cạn có nhu cầu nước thấp trước năm 2028, áp dụng công nghệ tưới tiết kiệm nước nhằm cắt giảm 10% lượng nước tưới mặt trong mùa khô.

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

Luận văn là tài liệu chuyên khảo có giá trị ứng dụng cao cho nhiều nhóm đối tượng:

  1. Cơ quan quản lý nhà nước về tài nguyên và môi trường (Bộ TN&MT, các Sở TN&MT địa phương): Vận dụng số liệu dự báo xu thế dòng chảy để tích hợp vào Quy hoạch tổng hợp lưu vực sông Cầu và xây dựng Kế hoạch hành động ứng phó với biến đổi khí hậu cấp tỉnh giai đoạn 2025–2035.
  2. Các đơn vị khai thác công trình thủy lợi và cấp nước đô thị: Tham khảo chỉ số biến động lưu lượng kiệt ($Q_{min}$) và mực nước đỉnh lũ ($h_{max}$) để tối ưu hóa thiết kế công trình thu nước, gia cố thân đê và lập phương án cấp nước an toàn cho các khu công nghiệp tại Thái Nguyên, Bắc Ninh.
  3. Các viện nghiên cứu và nhà khoa học chuyên ngành thủy văn - khí hậu: Kế thừa bộ thông số mô hình NAM đã được hiệu chỉnh và chuỗi số liệu khí tượng thủy văn 40 năm để mở rộng nghiên cứu sang các lưu vực lân cận thuộc hệ thống sông Thái Bình.
  4. Giảng viên, học viên cao học và sinh viên ngành Quản lý Môi trường, Biến đổi khí hậu, Kỹ thuật Tài nguyên nước: Sử dụng luận văn như một tài liệu tham khảo chuẩn mực về phương pháp luận nghiên cứu, quy trình ứng dụng GIS kết hợp mô hình toán thủy văn trong đánh giá tác động môi trường.

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

Biến đổi khí hậu ảnh hưởng lớn nhất đến yếu tố nào của dòng chảy sông Cầu?
Biến đổi khí hậu tác động rõ rệt nhất đến sự phân phối dòng chảy theo mùa và gia tăng tính cực đoan. Lưu lượng dòng chảy mùa lũ dự kiến tăng từ 4% đến 8,6% gây nguy cơ ngập lụt, trong khi dòng chảy mùa cạn suy giảm từ 3,8% đến 7,5%, làm trầm trọng thêm tình trạng thiếu nước sinh hoạt và sản xuất ở hạ lưu.

Mô hình thủy văn NAM có ưu điểm gì trong đánh giá dòng chảy lưu vực?
Mô hình NAM thuộc bộ mô hình MIKE 11 có cấu trúc khối bể chứa đơn giản nhưng mô phỏng chính xác chu trình hình thành dòng chảy từ mưa liên tục. Mô hình xử lý tốt mối quan hệ giữa lượng mưa, bốc thoát hơi và độ ẩm đất, đạt hệ số tương quan kiểm định trên 0,85 phù hợp với đặc thù lưu vực sông Cầu.

Nhiệt độ trên lưu vực sông Cầu đã biến đổi như thế nào trong 40 năm qua?
Nhiệt độ trung bình năm giai đoạn 1974–2013 tăng từ 0,2°C đến 0,3°C mỗi thập kỷ. Mức tăng mạnh nhất diễn ra vào mùa đông (tháng 1 và tháng 2) với hệ số góc trung bình +0,035°C/năm. Nhiệt độ cực tiểu tăng nhanh hơn nhiệt độ cực đại, khiến số ngày rét đậm giảm rõ rệt.

Kịch bản biến đổi khí hậu A1B đưa ra dự báo gì cho giai đoạn 2014–2033?
Theo kịch bản phát thải trung bình A1B, nhiệt độ trung bình trên lưu vực sông Cầu tiếp tục tăng từ 0,4°C đến 0,6°C vào năm 2033. Lượng mưa năm có xu thế tăng nhẹ từ 2% đến 5% nhưng tập trung chủ yếu vào mùa mưa bão, làm tăng biên độ dao động giữa dòng chảy đỉnh lũ và dòng chảy đáy kiệt.

Giải pháp công trình nào cấp thiết nhất để điều tiết dòng chảy sông Cầu?
Giải pháp cấp thiết nhất là hoàn thiện quy chế vận hành điều tiết liên hồ chứa trên lưu vực kết hợp cải tạo, nạo vét các đoạn sông bị bồi lắng. Biện pháp này giúp chủ động tích trữ nước trong mùa mưa để bổ sung lưu lượng dòng chảy tối thiểu từ 15% đến 20% trong các tháng cao điểm mùa khô.

Kết luận

  • Luận văn đã xác lập bức tranh khoa học toàn diện về diễn biến biến đổi khí hậu và tác động đến dòng chảy mặt lưu vực sông Cầu qua chuỗi quan trắc 40 năm.
  • Xu thế tăng nhiệt độ mùa đông và sự gia tăng các đợt mưa lớn cực đoan trên 50 mm/ngày là những biểu hiện khí hậu rõ nét nhất trên lưu vực.
  • Sự phân hóa dòng chảy theo mùa ngày càng gay gắt, với dòng chảy mùa lũ tăng đến 8,6% và dòng chảy mùa cạn suy giảm tới 7,5% trong giai đoạn tương lai.
  • Việc kết hợp hiệu quả công cụ GIS và mô hình thủy văn NAM đã lượng hóa thành công các kịch bản thủy văn phục vụ trực tiếp cho công tác quy hoạch tài nguyên nước.
  • Đóng góp chính của đề tài là cung cấp cơ sở dữ liệu định lượng chính xác, tạo tiền đề để các tỉnh lưu vực triển khai giải pháp công trình và phi công trình giai đoạn 2025–2030.

Để ứng phó kịp thời với các tác động tiêu cực của biến đổi khí hậu, các cơ quan quản lý và địa phương trên lưu vực sông Cầu cần nhanh chóng tích hợp các phát hiện của luận văn vào quy hoạch phát triển kinh tế - xã hội, chủ động đầu tư mạng lưới quan trắc và thúc đẩy các biện pháp quản lý tổng hợp tài nguyên nước ngay từ hôm nay.