Tổng quan nghiên cứu

Lưu vực sông Phan – Cà Lồ thuộc tỉnh Vĩnh Phúc có tổng diện tích tự nhiên khoảng 710 km², chiếm gần 2/3 diện tích toàn tỉnh, giữ vai trò huyết mạch trong việc cung cấp nguồn nước tưới tiêu cho 38.200 ha đất canh tác nông nghiệp và phục vụ đời sống dân sinh cho hơn 823.529 người dân với mật độ dân số cao 1.160 người/km². Tuy nhiên, dưới tác động ngày càng gia tăng của biến đổi khí hậu toàn cầu và áp lực phát triển kinh tế với 2.365 ha đất công nghiệp, tình trạng khô hạn và thiếu hụt nguồn nước mặt trong mùa khô đang diễn ra hết sức phức tạp. Lượng mưa tại lưu vực phân bố không đều khi có tới 80% tổng lượng mưa tập trung trong mùa mưa từ tháng 5 đến tháng 10, trong khi mùa khô kéo dài từ tháng 11 đến tháng 4 năm sau chỉ chiếm 20% lượng mưa cả năm.

Thực tế quan trắc cho thấy vào các năm hạn hán gay gắt như vụ đông xuân 2010, toàn vùng nghiên cứu thiếu hụt tới 40 triệu m³ nước tưới so với tổng nhu cầu 218 triệu m³. Tương tự, vào vụ đông xuân 2015-2016, lượng nước thiếu hụt hơn 5 triệu m³ đã khiến dung tích hữu ích của các hồ chứa trọng điểm sụt giảm nghiêm trọng, trong đó hồ Xạ Hương thiếu hụt 4,35 triệu m³ và hồ Thanh Lanh thiếu hụt 1,66 triệu m³.

Xuất phát từ thực tiễn trên, luận văn tập trung vào mục tiêu xác định và ứng dụng bộ chỉ số đánh giá khô hạn chuẩn hóa lượng mưa SPI và chỉ số cân bằng ẩm Kappa, kết hợp công cụ Hệ thông tin địa lý GIS để giám sát diễn biến thiếu hụt nguồn nước theo không gian và thời gian. Phạm vi nghiên cứu bao trùm 7 đơn vị hành chính của lưu vực sông Phan – Cà Lồ với chuỗi số liệu quan trắc 30 năm từ 1985 đến 2015. Kết quả nghiên cứu mang ý nghĩa khoa học và thực tiễn sâu sắc, cung cấp cơ sở định lượng giúp các cơ quan quản lý nâng cao độ chính xác trong dự báo thiếu nước lên trên 85% và tối ưu hóa quy trình vận hành hệ thống công trình thủy lợi.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Luận văn tiếp cận toàn diện dựa trên lý thuyết phân loại hạn hán của Tổ chức Khí tượng Thế giới WMO, bao gồm 4 nhóm chính: hạn khí tượng, hạn nông nghiệp, hạn thủy văn và hạn kinh tế - xã hội. Trong đó, hạn khí tượng được xác định là nguyên nhân khởi nguồn trực tiếp dẫn đến sự thiếu hụt độ ẩm trong đất và cạn kiệt dòng chảy mặt.

Mô hình nghiên cứu trọng tâm là Chỉ số chuẩn hóa lượng mưa SPI do McKee và các cộng sự đề xuất năm 1993. Chỉ số SPI được xây dựng trên cơ sở chuẩn hóa chuỗi phân bố xác suất lượng mưa thực tế về dạng phân bố chuẩn tắc Gauss thông qua hàm mật độ xác suất Gamma hai tham số gồm tham số hình dạng alpha và tham số tỉ lệ beta. Chỉ số này cho phép đánh giá mức độ ẩm ướt hoặc khô hạn theo các thang thời gian linh hoạt gồm SPI1 (hạn theo tháng), SPI3 (hạn theo mùa vụ 3 tháng) và SPI6 (hạn theo chu kỳ 6 tháng) với thang phân cấp từ cực kỳ ẩm ướt (SPI ≥ 2,0) đến hạn cực nặng (SPI ≤ -2,0).

Bên cạnh đó, nghiên cứu ứng dụng chỉ số cân bằng ẩm Kappa nhằm đánh giá tương quan giữa lượng mưa thực tế và lượng bốc hơi tiềm năng. Chỉ số này đặc biệt có ý nghĩa khi so sánh giữa các tiểu vùng sinh thái có sự chênh lệch lớn về điều kiện vi khí hậu, chẳng hạn như trạm Tam Đảo có khí hậu núi ẩm với lượng bốc hơi bình quân thấp và trạm Vĩnh Yên thuộc vùng trung du có lượng bốc hơi trung bình năm lên tới 1.040 mm.

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

Nghiên cứu sử dụng cỡ mẫu dữ liệu khí tượng thủy văn hoàn chỉnh trong chuỗi thời gian 30 năm liên tục từ năm 1985 đến năm 2015. Dữ liệu được thu thập từ 10 trạm quan trắc khí tượng đại diện phân bố đều trong và lân cận lưu vực, bao gồm: Tam Đảo, Vĩnh Phúc (Vĩnh Yên), Phú Thọ, Việt Trì, Kỳ Phú, Đông Anh, Lâm Thao, Phù Ninh, Thanh Thủy và Sơn Tây.

Phương pháp chọn mẫu tuân thủ nghiêm ngặt tiêu chuẩn khí tượng thế giới, ưu tiên các trạm có độ dài quan trắc trên 30 năm, độ tin cậy cao và tính đại diện cho ba vùng địa hình đặc thù: vùng núi cao (độ cao trên 300 m), vùng đồi gò trung du và vùng đồng bằng bãi bồi ven sông.

Lý do lựa chọn phương pháp phân tích kết hợp giữa thống kê toán học xác suất và Hệ thông tin địa lý GIS là vì chỉ số SPI có khả năng loại bỏ ảnh hưởng của tính mùa vụ và sai lệch địa hình, đồng thời thang thời gian SPI3 tương thích hoàn hảo với chu kỳ canh tác mùa vụ nông nghiệp tại địa phương. Công cụ GIS cho phép thực hiện nội suy không gian để chuyển đổi các điểm dữ liệu rời rạc thành bản đồ phân vùng hạn hán trực quan, định lượng chính xác diện tích và mức độ thiếu hụt nguồn nước trên toàn bộ 710 km² lưu vực.

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

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

Thứ nhất, kết quả tính toán chuỗi chỉ số SPI qua 30 năm đã nhận diện rõ quy luật hạn hán theo thời gian. Mùa khô kéo dài từ tháng 11 đến tháng 3 năm sau là giai đoạn xuất hiện hạn hán gay gắt nhất, chiếm trên 80% tổng số đợt thiếu hụt nguồn nước mặt. Các đợt hạn cực đoan nhất trong lịch sử được ghi nhận vào các năm 1988, 1998, 2003, 2005, 2010 và 2015. Trong các giai đoạn này, chỉ số SPI3 đồng loạt giảm sâu xuống dưới ngưỡng -1,5 (hạn nặng) và thậm chí dưới -2,0 (hạn cực nặng), phản ánh mức thiếu hụt lượng mưa trên 50% so với trung bình nhiều năm trong 3 tháng liên tiếp.

Thứ hai, nghiên cứu phát hiện sự phân hóa không gian khô hạn rất rõ rệt giữa các tiểu vùng địa hình. Vùng có nguy cơ hạn hán cao nhất chiếm khoảng 30% đến 35% diện tích lưu vực, tập trung chủ yếu tại các xã chân núi Tam Đảo thuộc các huyện Tam Dương, Tam Đảo, Bình Xuyên và các vùng bãi bồi ngoài đê thuộc huyện Vĩnh Tường, Yên Lạc. Cụ thể, khu vực kênh Hữu thuộc hệ thống Liễn Sơn có khoảng 670 ha và vùng bãi sông Phó Đáy có khoảng 500 ha thường xuyên rơi vào tình trạng thiếu nước tưới nghiêm trọng do địa hình chia cắt và mực nước sông bị hạ thấp.

Thứ ba, kết quả so sánh giữa chỉ số SPI và chỉ số Kappa khẳng định tính đồng thuận cao trong việc nhận dạng các năm hạn hán lịch sử. Chỉ số SPI3 phản ánh chính xác 85% các đợt hạn thực tế tác động đến sản xuất lúa vụ đông xuân. Trong khi đó, chỉ số Kappa thể hiện rõ rệt sự gia tăng độ căng thẳng bốc thoát hơi nước trong các tháng mùa hè tại khu vực trung du Vĩnh Yên (lượng bốc hơi trung bình tháng đạt 107,58 mm) so với vùng núi cao Tam Đảo (lượng bốc hơi trung bình tháng đạt 50,2 mm).

Thảo luận kết quả

Nguyên nhân cốt lõi dẫn đến tình trạng khô hạn nghiêm trọng tại lưu vực sông Phan – Cà Lồ là sự kết hợp giữa hiện tượng khí hậu cực đoan El Nino và sự phân bố địa hình dốc theo hướng Tây Bắc – Đông Nam. Trong các năm El Nino hoạt động mạnh như 1997-1998 và 2014-2016, nền nhiệt độ mùa hè tăng vọt lên mức 38-42°C khiến lượng bốc hơi tăng mạnh, đồng thời mùa mưa kết thúc sớm làm dòng chảy kiệt sụt giảm từ 60% đến 90%. Ngoài ra, mạng lưới sông ngòi quanh co với hệ số uốn khúc K xấp xỉ 1,8 cùng với hơn 102 công trình cầu, cống ngăn trên trục chính đã làm giảm vận tốc và dung tích điều tiết tự nhiên của dòng chảy.

Các tập dữ liệu nghiên cứu được trực quan hóa chi tiết thông qua các biểu đồ đường thể hiện xu thế biến đổi lượng mưa mùa khô giai đoạn 1985-2015, bảng thống kê phân loại tần suất hạn hán theo các ngưỡng SPI và chuỗi bản đồ số GIS phân vùng diện tích hạn SPI3 cho các năm điển hình 1988, 1995, 2005, 2010 và 2015. Bản đồ GIS đã mô tả sinh động sự lan rộng của các vùng hạn từ sườn núi Tam Đảo lan tỏa xuống khu vực hạ lưu sông Cà Lồ.

So sánh với các công trình nghiên cứu của các tác giả tại vùng đồng bằng sông Hồng và lưu vực sông Cái, kết quả của luận văn hoàn toàn tương đồng về tính chu kỳ và mức độ nhạy cảm của chỉ số SPI đối với hạn nông nghiệp. Điều này khẳng định việc tích hợp chỉ số SPI đa thời đoạn và công nghệ GIS là giải pháp khoa học tối ưu, cho phép giám sát toàn diện tình trạng thiếu hụt nguồn nước trên các lưu vực sông vừa và nhỏ.

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

Nâng cấp và cải tạo hệ thống hạ tầng công trình thủy lợi: Tiến hành nạo vét toàn bộ 86 km trục chính sông Phan – Cà Lồ và sửa chữa, hiện đại hóa 76 công trình điều tiết nước trên tuyến sông Phan; mục tiêu đảm bảo duy trì lưu lượng dòng chảy kiệt tối thiểu 30 m³/s trong mùa khô; lộ trình triển khai từ năm 2026 đến năm 2028; chủ thể thực hiện là Sở Nông nghiệp và Phát triển nông thôn tỉnh Vĩnh Phúc phối hợp cùng Chi cục Thủy lợi.

Hoàn thiện quy trình vận hành liên hồ chứa thông minh: Xây dựng và ban hành quy chế phối hợp tích nước, xả nước hợp lý cho các hồ chứa lớn gồm Đại Lải, Xạ Hương, Thanh Lanh và Làng Hà; mục tiêu đảm bảo tích trữ đạt 100% dung tích thiết kế trước tháng 11 hàng năm, bù đắp lượng thiếu hụt từ 5 đến 40 triệu m³ nước tưới cho vụ đông xuân; lộ trình thực hiện định kỳ hàng năm từ năm 2026 đến năm 2030; chủ thể thực hiện là Công ty TNHH MTV Thủy lợi Liễn Sơn và Công ty TNHH MTV Thủy lợi Tam Đảo.

Xây dựng hệ thống giám sát và cảnh báo hạn hán thời gian thực: Thiết lập phần mềm chuyên dụng tích hợp dữ liệu quan trắc tự động, chỉ số SPI3, SPI6 và công nghệ GIS; mục tiêu nâng cao độ chính xác cảnh báo sớm trước 15-30 ngày đạt trên 90%; lộ trình hoàn thành trong giai đoạn 2026-2027; chủ thể thực hiện là Trung tâm Dự báo Khí tượng Thủy văn tỉnh phối hợp với các viện nghiên cứu chuyên ngành.

Chuyển đổi cơ cấu cây trồng và nhân rộng công nghệ tưới tiết kiệm: Rà soát và chuyển đổi khoảng 1.500 ha đến 2.000 ha diện tích lúa thường xuyên bị hạn tại vùng kênh Hữu và vùng bãi sông Phó Đáy sang trồng các loại cây chịu hạn có giá trị kinh tế cao như thanh long, cây ăn quả và ngô lai; mục tiêu cắt giảm từ 20% đến 30% định mức nước tưới; lộ trình thực hiện từ năm 2026 đến năm 2029; chủ thể thực hiện là Ủy ban nhân dân các huyện Tam Dương, Tam Đảo, Vĩnh Tường, Yên Lạc cùng các hợp tác xã nông nghiệp.

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

Cán bộ quản lý và hoạch định chính sách tại các cơ quan quản lý tài nguyên nước: Lãnh đạo và chuyên viên thuộc Sở Nông nghiệp và Phát triển nông thôn, Sở Tài nguyên và Môi trường tỉnh Vĩnh Phúc có thể ứng dụng các bản đồ phân vùng hạn hán GIS để xây dựng quy hoạch phân bổ tài nguyên nước hợp lý cho 38.200 ha đất nông nghiệp và các khu công nghiệp.

Kỹ sư vận hành hệ thống công trình thủy lợi: Đội ngũ kỹ sư tại các công ty khai thác công trình thủy lợi (Liễn Sơn, Tam Đảo, Phúc Yên) có thể áp dụng các ngưỡng giá trị chỉ số SPI3 và Kappa để chủ động lập kế hoạch tích trữ nước tại các hồ đập và điều hành trạm bơm trong các thời kỳ khô hạn cực đoan.

Giảng viên, nhà khoa học và học viên cao học chuyên ngành Thủy văn và Kỹ thuật tài nguyên nước: Các nhà nghiên cứu có thể kế thừa phương pháp chuẩn hóa xác suất Gamma trên chuỗi số liệu 30 năm và quy trình tích hợp GIS để phát triển các mô hình dự báo hạn hán chuyên sâu cho các lưu vực sông khác trên toàn quốc.

Chính quyền địa phương và cán bộ khuyến nông cơ sở: Ủy ban nhân dân các cấp và cán bộ nông nghiệp tại 7 huyện, thành phố vùng nghiên cứu có thể sử dụng kết quả đánh giá vùng hạn trọng điểm để hướng dẫn nông dân điều chỉnh lịch thời vụ và bố trí cây trồng chịu hạn thích ứng hiệu quả với thiên tai.

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

Tại sao luận văn lại lựa chọn chỉ số chuẩn hóa lượng mưa SPI làm công cụ chủ đạo để giám sát khô hạn? Chỉ số SPI được Tổ chức Khí tượng Thế giới WMO khuyến nghị sử dụng rộng rãi nhờ tính chuẩn hóa cao, cho phép so sánh khách quan tình trạng khô hạn giữa các vùng khí hậu khác nhau. Nghiên cứu đã chứng minh chỉ số SPI có khả năng tính toán linh hoạt trên chuỗi số liệu 30 năm tại 10 trạm đo, phản ánh chính xác tình trạng thiếu hụt lượng mưa mà không đòi hỏi quá nhiều thông số phức tạp.

Chỉ số SPI3 có ý nghĩa thực tiễn như thế nào đối với sản xuất nông nghiệp vụ đông xuân? Chỉ số SPI3 phản ánh lượng mưa tích lũy trong khoảng thời gian 3 tháng liên tiếp, trùng khớp hoàn toàn với chu kỳ sinh trưởng của cây trồng một mùa vụ. Khi giá trị SPI3 giảm xuống dưới -1,5 trong các năm 2010 và 2015, thực tế ghi nhận toàn tỉnh Vĩnh Phúc thiếu hụt từ 5 đến 40 triệu m³ nước tưới, gây khô hạn cục bộ cho hàng nghìn hecta lúa vụ đông xuân.

Những khu vực nào trong lưu vực sông Phan – Cà Lồ chịu ảnh hưởng hạn hán nặng nề nhất? Khu vực chịu rủi ro hạn hán cao nhất tập trung tại các xã vùng chân núi Tam Đảo thuộc huyện Tam Dương, Tam Đảo, Bình Xuyên do địa hình dốc chia cắt, và vùng kênh Hữu (khoảng 670 ha) cùng bãi bồi sông Phó Đáy (khoảng 500 ha) thuộc huyện Vĩnh Tường, Yên Lạc do nằm xa công trình đầu mối và nguồn nước mặt bị suy giảm.

Điểm khác biệt cơ bản giữa kết quả đánh giá của chỉ số SPI và chỉ số Kappa là gì? Chỉ số SPI tập trung phản ánh mức độ thiếu hụt lượng mưa khí tượng theo phân bố xác suất Gamma, trong khi chỉ số Kappa tích hợp thêm yếu tố lượng bốc hơi tiềm năng. Tại trạm Vĩnh Yên với lượng bốc hơi hàng năm lên tới 1.040 mm, chỉ số Kappa thể hiện rõ nét hơn mức độ căng thẳng ẩm độ so với trạm Tam Đảo có khí hậu vùng núi mát mẻ.

Công nghệ GIS đóng góp vai trò gì trong việc nâng cao chất lượng đánh giá hạn hán của luận văn? Hệ thống thông tin địa lý GIS đóng vai trò nội suy không gian các giá trị SPI từ 10 trạm quan trắc để xây dựng bản đồ chuyên đề phân vùng hạn hán các năm điển hình 1988, 1995, 2005, 2010 và 2015. Bản đồ GIS giúp trực quan hóa ranh giới các vùng hạn, xác định chính xác tỷ lệ diện tích bị ảnh hưởng trên tổng số 710 km² lưu vực.

Kết luận

  • Xây dựng thành công quy trình ứng dụng bộ chỉ số hạn SPI và Kappa trên chuỗi số liệu quan trắc 30 năm (1985–2015) từ 10 trạm khí tượng thủy văn đại diện cho toàn lưu vực.
  • Nhận diện chính xác quy luật phân hóa hạn hán theo thời gian với trọng tâm vào các tháng mùa khô từ tháng 11 đến tháng 3, làm rõ mức độ thiếu hụt nước từ 5 đến 40 triệu m³ trong các năm hạn cực đoan 2010 và 2015.
  • Ứng dụng công nghệ GIS thành lập hệ thống bản đồ số phân vùng khô hạn cho 710 km² diện tích lưu vực sông Phan – Cà Lồ, chỉ ra các vùng trọng điểm hạn hán tại Tam Đảo, Tam Dương, Bình Xuyên và Vĩnh Tường.
  • Đề xuất hệ thống 4 nhóm giải pháp công trình và phi công trình đồng bộ, tạo nền tảng vững chắc để đảm bảo an ninh nguồn nước cho 38.200 ha đất canh tác và hơn 823.529 người dân.
  • Đóng góp luận cứ khoa học quan trọng phục vụ công tác quy hoạch tài nguyên nước, quản lý rủi ro thiên tai và hiện đại hóa quy trình vận hành 102 công trình thủy lợi trên hệ thống sông Phan – Cà Lồ.

Trong giai đoạn 2026–2030, các nghiên cứu tiếp theo cần tích hợp thêm ảnh vệ tinh viễn thám và dữ liệu dự báo khí hậu toàn cầu để mở rộng khả năng cảnh báo sớm hạn hán theo thời gian thực. Đề nghị các cơ quan quản lý nhà nước và doanh nghiệp khai thác thủy lợi tại tỉnh Vĩnh Phúc nhanh chóng tiếp nhận, ứng dụng kết quả luận văn vào công tác điều hành cấp nước và ứng phó biến đổi khí hậu tại địa phương.