Tổng quan nghiên cứu

Hạn hán là một trong những hiểm họa tự nhiên khốc liệt nhất, đứng thứ ba trong các loại hình thiên tai phổ biến tại Việt Nam sau bão và lũ lụt. Thống kê toàn cầu giai đoạn 1900-2013 ghi nhận 642 đợt hạn hán lớn, cướp đi sinh mạng của khoảng 12 triệu người, ảnh hưởng trực tiếp tới hơn 2 tỷ người và gây tổn thất kinh tế ước tính vượt 135 tỷ USD. Tại Việt Nam, sự suy giảm lượng mưa kết hợp với nền nhiệt độ gia tăng bất thường đã làm suy kiệt dòng chảy sông ngòi từ 0,5 đến 1,5m dưới mức trung bình nhiều năm, đồng thời đẩy mặn xâm nhập sâu 15-20km vào nội đồng tại các khu vực duyên hải và đồng bằng.

Vấn đề nghiên cứu trọng tâm của luận văn tập trung vào việc phân tích đặc điểm, cơ chế tác động của các nhân tố hoàn lưu khí quyển quy mô lớn đối với điều kiện hạn hán trên lãnh thổ Việt Nam. Mục tiêu cụ thể là lượng hóa tần suất và mức độ khô hạn, xác định vai trò của chu trình ENSO, áp cao cận nhiệt đới Tây Thái Bình Dương, áp thấp nóng phía tây và áp cao lục địa Siberi, từ đó xây dựng cơ sở khoa học cho công tác dự báo sớm và dự tính xu thế hạn trong tương lai.

Phạm vi không gian nghiên cứu bao quát toàn bộ 7 phân vùng khí hậu cùng hệ thống các trạm đảo tiêu biểu của Việt Nam, với chuỗi số liệu quan trắc thực tế kéo dài 34 năm (1981-2014) và dự tính khí hậu giai đoạn 2016-2027. Kết quả nghiên cứu mang ý nghĩa thực tiễn to lớn trong việc giảm thiểu nguy cơ thâm hụt 30-70% lượng mưa tự nhiên, hỗ trợ hoạch định chiến lược an ninh nguồn nước và phòng chống thiên tai cấp quốc gia.

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 vấn đề dựa trên nền tảng lý thuyết hoàn lưu khí quyển quy mô lớn và sự tương tác biển - khí quyển. Cơ chế hình thành thời tiết khô hạn được xem xét thông qua sự biến động của hoàn lưu Walker trong chu trình ENSO (El Nino - Dao động Nam), áp cao cận nhiệt đới Tây Thái Bình Dương ở mực 850mb, áp thấp nóng phía tây (áp thấp Nam Á hay áp thấp Ấn Miến), hoàn lưu áp cao lạnh lục địa Siberi và nhánh dòng giáng của hoàn lưu Hadley.

Về mặt khái niệm, công trình kế thừa cách phân loại của Wilhite và Glantz (1985), chia hạn hán thành 4 cấp độ: hạn khí tượng (thiếu hụt giáng thủy kéo dài), hạn nông nghiệp (thiếu hụt độ ẩm đất nuôi dưỡng cây trồng), hạn thủy văn (suy kiệt dòng chảy mặt và nước ngầm) và hạn kinh tế - xã hội (mất cân bằng cung cầu nước phục vụ sản xuất và sinh hoạt).

Để định lượng mức độ khô hạn, luận văn áp dụng 3 mô hình chỉ số hạn khí tượng tiêu biểu: Chỉ số chuẩn hóa lượng mưa (SPI - Standardized Precipitation Index theo McKee và cộng sự, 1993), Chỉ số khô hạn De Martonne (chỉ số J, 1926) và Chỉ số khô cằn khí quyển Ped (Pedey, 1975). Các chỉ số này cho phép lượng hóa chính xác điều kiện thiếu hụt lượng mưa đơn lẻ cũng như sự kết hợp tương tác giữa nhiệt độ và bốc thoát hơi nước.

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

Nguồn dữ liệu của luận văn được thu thập từ mạng lưới 55 trạm khí tượng bề mặt đại diện cho 7 vùng khí hậu (Tây Bắc, Đông Bắc, Đồng bằng Bắc Bộ, Bắc Trung Bộ, Nam Trung Bộ, Tây Nguyên, Nam Bộ) và 5 trạm đảo tiền tiêu (Bạch Long Vỹ, Hoàng Sa, Trường Sa, Cồn Cỏ, Côn Đảo) từ Trung tâm Dự báo Khí tượng Thủy văn Quốc gia. Cỡ mẫu nghiên cứu gồm chuỗi số liệu nhiệt độ trung bình tháng và tổng lượng mưa tháng liên tục 34 năm (1981-2014) với tổng cộng 408 bước thời gian tháng. Dữ liệu nhiệt độ mặt nước biển (SST) và chỉ số Oceanic Niño Index (ONI) được trích xuất từ Cục Quản lý Đại dương và Khí quyển Quốc gia Hoa Kỳ (NOAA).

Phương pháp chọn mẫu áp dụng nguyên tắc phân vùng khí hậu chuẩn của Việt Nam nhằm phản ánh đầy đủ tính đa dạng địa hình và vi khí hậu. Luận văn sử dụng phương pháp thống kê khí hậu, tính toán chuẩn sai, phân tích xác suất tần suất xuất hiện hạn (%) theo 4 cấp độ (không hạn, hạn nhẹ, hạn vừa, hạn nặng) và phân tích tương quan chuỗi thời gian có độ trễ.

Lý do lựa chọn đồng thời 3 chỉ số J, SPI và Ped là nhằm phát huy ưu điểm bổ trợ lẫn nhau: chỉ số SPI tối ưu cho việc theo dõi lượng mưa độc lập, chỉ số J nhạy bén trong việc xác định thời điểm bắt đầu hạn, còn chỉ số Ped phản ánh hiệu quả tác động kép của nhiệt độ cao và giáng thủy suy giảm. Ngoài ra, mô hình khí hậu NCAR-CCSM được ứng dụng để dự tính mức độ hạn hán tương lai (2016-2027) theo các kịch bản phát thải RCP 2.6, RCP 4.5 và RCP 8.5.

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

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

Phân tích định lượng chuỗi số liệu 1981-2014 cho thấy tần suất xuất hiện hạn hán có sự phân hóa mạnh mẽ theo từng chỉ số và không gian địa lý. Chỉ số J luôn cho tần suất xuất hiện hạn cao nhất tại 4 vùng phía Bắc (40-45%) và vượt mức 50% tại 3 vùng phía Nam, trong khi chỉ số SPI duy trì ổn định ở mức 35-40% và chỉ số Ped đạt khoảng 25% (tỷ lệ không hạn đạt 73-76%).

Về phân bố cấp độ hạn, các vùng khí hậu phía Nam ghi nhận tỷ lệ hạn nặng cao hơn đáng kể so với phía Bắc. Trong đó, vùng Tây Nguyên và Nam Bộ chịu áp lực hạn nặng lớn nhất, còn vùng Nam Trung Bộ tập trung chủ yếu ở ngƣỡng hạn nhẹ và hạn vừa. Điểm đặc thù diễn ra tại Bắc Trung Bộ khi tần suất hạn vừa và hạn nặng xác định theo chỉ số SPI và Ped lại vượt trội hơn so với chỉ số J.

Biến trình hạn theo 12 tháng thể hiện tính quy luật mùa sâu sắc. Tại Tây Bắc, Đông Bắc và Đồng bằng Bắc Bộ, hạn tập trung cao độ từ tháng 11 đến tháng 3-4 năm sau với tần suất trên 70% đến 94% theo chỉ số J, trong khi từ tháng 6 đến tháng 8 tần suất hạn gần như bằng 0%. Vùng Bắc Trung Bộ có hai cực đại hạn vào tháng 2-3 (88%) và tháng 6-7 (50%), hoàn toàn không có hạn trong tháng 9-10. Tại Tây Nguyên và Nam Bộ, hạn hán mùa khô khốc liệt đạt đỉnh điểm trong giai đoạn tháng 1-3 với tần suất hạn theo chỉ số J lên tới 80-100%.

Tại các trạm đảo, trạm Bạch Long Vỹ và Hoàng Sa ghi nhận tần suất hạn theo chỉ số J từ 55% đến 65%, trong khi các trạm đảo phía Nam (Trường Sa, Côn Đảo) dao động 40-45%. Hoạt động của ENSO thể hiện tác động cực đoan: trong đợt El Nino mạnh 1997-1998, lượng mưa thiếu hụt 30-70% so với trung bình nhiều năm, nền nhiệt độ tăng 1-3°C trên diện rộng; ngược lại, giai đoạn La Nina 1988-1989 mang lại lượng mưa dồi dào và làm suy giảm rõ rệt tần suất khô hạn.

Thảo luận kết quả

Cơ chế vật lý đằng sau các phát hiện trên được lý giải bởi sự biến động đồng thời của các hệ thống hoàn lưu khí quyển. Trong các năm El Nino, hoàn lưu Walker hoạt động suy yếu khiến trung tâm đối lưu sâu dịch chuyển về phía đông, làm giảm số lượng bão trên Biển Đông và khiến mùa mưa tại Việt Nam kết thúc sớm khoảng 1 tháng. Đồng thời, sự mở rộng bất thường của áp thấp nóng phía tây kết hợp với trường gió giáng ở tầng cao từ áp cao cận nhiệt đới Tây Thái Bình Dương đã triệt tiêu sự hình thành mây, đẩy nhiệt độ mặt đất lên cao và kéo dài chuỗi ngày không mưa.

Sự khác biệt tại Bắc Trung Bộ vào các tháng mùa hè bắt nguồn từ hiệu ứng gió Phơn khô nóng khi đới gió Tây Nam vượt qua dãy Trường Sơn. Trái lại, vào mùa đông, đới gió mùa Đông Bắc đem lại mưa lớn ở sườn đón gió ven biển miền Trung nhưng lại trở nên khô kiệt khi tràn sang khu vực Tây Nguyên và cực Nam Trung Bộ, giải thích vì sao tần suất hạn tại Tây Nguyên đạt gần như tuyệt đối (80-100%) vào đầu năm.

Các kết quả này củng cố các nghiên cứu trước đây về độ trễ 4 tháng giữa nhiệt độ bề mặt nước biển khu vực Nino 3.4 và phản ứng nhiệt độ tại Việt Nam. Toàn bộ diễn biến phức tạp trên được trực quan hóa mạch lạc trong luận văn thông qua các biểu đồ đường tần suất hạn 12 tháng, bản đồ phân bố chuẩn sai nhiệt độ bề mặt nước biển (SST) trên Biển Đông qua các thập kỷ và hệ thống bảng ma trận so sánh cấp độ hạn giữa 3 chỉ số J, SPI và Ped.

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

Thứ nhất, nâng cấp hệ thống cảnh báo sớm và giám sát hạn hán đa chỉ số: Trung tâm Dự báo Khí tượng Thủy văn Quốc gia cần chủ trì tích hợp bộ ba chỉ số J, SPI và Ped vào quy trình dự báo nghiệp vụ hàng tháng. Khai thác mối quan hệ tương quan có độ trễ 1-4 tháng giữa chỉ số ONI, chuẩn sai SST Biển Đông và hoàn lưu áp cao cận nhiệt đới nhằm nâng cao độ tin cậy dự báo hạn mùa đạt trên 85% trong giai đoạn 2026-2028.

Thứ hai, tối ưu hóa quy trình điều tiết liên hồ chứa và trữ nước thượng nguồn: Bộ Nông nghiệp và Phát triển Nông thôn phối hợp cùng Tập đoàn Điện lực Việt Nam và các công ty thủy điện rà soát quy trình vận hành liên hồ chứa trên lưu vực sông Hồng, sông Mã, sông Vu Gia - Thu Bồn và sông Ba. Mục tiêu là chủ động tích nước sớm trước khi mùa lũ kết thúc 15-20 ngày, giảm thiểu 40% nguy cơ suy kiệt mực nước hồ chứa về mực nước chết trong các tháng cao điểm mùa khô (từ tháng 1 đến tháng 4 hàng năm).

Thứ ba, chuyển đổi cơ cấu cây trồng và mùa vụ thích ứng hạn tại các vùng trọng điểm: Sở Nông nghiệp và Phát triển Nông thôn các tỉnh Nam Trung Bộ, Tây Nguyên và Nam Bộ cần cơ cấu lại lịch mùa vụ, ưu tiên mở rộng các giống cây trồng ngắn ngày có khả năng chống chịu khô hạn cao. Áp dụng kỹ thuật tưới tiết kiệm nước nhằm cắt giảm 25-30% nhu cầu nước tưới tiêu trong các năm xuất hiện hiện tượng El Nino giai đoạn 2026-2030.

Thứ tư, hoàn thiện công trình ngăn mặn và trữ ngọt thích ứng hạn thủy văn: Ủy ban Nhân dân các tỉnh ven biển Bắc Trung Bộ và Đồng bằng sông Cửu Long cần triển khai nạo vét hệ thống kênh mương, xây dựng đập tạm ngăn mặn trước mùa khô, khống chế ranh mặn 4g/l không xâm nhập sâu quá 15km vào các vùng sản xuất nông nghiệp và vùng cấp nước sinh hoạt trọng điểm.

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

Nhóm chuyên gia và dự báo viên khí tượng thủy văn: Nghiên cứu cung cấp phương pháp luận chuẩn xác trong việc tính toán, so sánh độ nhạy của 3 chỉ số J, SPI, Ped kết hợp phân tích các chỉ số khí quyển - đại dương quy mô lớn (ENSO, SST, ONI). Đây là tài liệu nghiệp vụ giá trị để cải tiến mô hình dự báo mùa và cảnh báo sớm thiên tai hạn hán tại các đài khu vực.

Nhóm cán bộ quản lý tài nguyên nước và công trình thủy lợi: Luận văn cung cấp số liệu định lượng chi tiết về tần suất hạn theo từng tháng và chu kỳ thiếu hụt dòng chảy 34 năm trên 7 vùng khí hậu. Nguồn dữ liệu này là cơ sở then chốt để lập quy hoạch khai thác nguồn nước, xây dựng kế hoạch phân bổ nước tưới và điều hòa dung tích các hồ chứa lớn.

Nhóm nhà hoạch định chính sách nông nghiệp và phát triển nông thôn: Các phân tích về phân hóa không gian hạn hán giữa các miền giúp các nhà quản lý xây dựng chiến lược chuyển đổi cơ cấu cây trồng, bố trí mùa vụ hợp lý cho từng địa phương nhằm giảm thiểu thiệt hại kinh tế do hạn hán nông nghiệp gây ra.

Nhóm giảng viên, học viên cao học và sinh viên chuyên ngành khoa học Trái đất: Công trình là tài liệu học thuật mẫu mực về cấu trúc nghiên cứu khí hậu học, kết hợp hài hòa giữa khai thác chuỗi số liệu quan trắc dài hạn 1981-2014 và ứng dụng mô hình khí hậu toàn cầu (NCAR-CCSM) dự tính theo các kịch bản biến đổi khí hậu RCP.

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

Luận văn sử dụng những chỉ số nào để đánh giá và lượng hóa hạn hán tại Việt Nam? Nghiên cứu ứng dụng đồng thời 3 chỉ số gồm: Chỉ số chuẩn hóa lượng mưa (SPI) phân tích xác suất thiếu hụt mưa; Chỉ số khô hạn De Martonne (J) đánh giá cân bằng nhiệt - ẩm tháng; và Chỉ số khô cằn khí quyển Ped tính toán độ lệch chuẩn của cả nhiệt độ và giáng thủy trên chuỗi số liệu 1981-2014.

Hiện tượng El Nino tác động cụ thể ra sao đến tình hình hạn hán tại Việt Nam? Trong các năm El Nino mạnh điển hình như 1997-1998, hoàn lưu Walker suy yếu khiến mùa mưa kết thúc sớm 1 tháng, tổng lượng mưa thâm hụt 30-70% so với trung bình nhiều năm, nền nhiệt độ tăng cao từ 1-3°C và xâm nhập mặn lấn sâu 15-20km vào đất liền.

Vùng khí hậu nào tại Việt Nam chịu rủi ro hạn hán nghiêm trọng nhất? Khu vực Tây Nguyên và Nam Bộ chịu rủi ro hạn nặng cao nhất cả nước trong các tháng mùa khô (tháng 1 đến tháng 3), với tần suất hạn xác định theo chỉ số J đạt từ 80% đến 100% do gió mùa mùa đông biến tính khô kiệt sau khi vượt dãy Trường Sơn.

Vì sao vùng khí hậu Bắc Trung Bộ lại có quy luật hạn hán khác biệt so với các vùng phía Bắc còn lại? Bắc Trung Bộ xuất hiện hai cực đại hạn vào tháng 2-3 (88%) và tháng 6-7 (50%). Đỉnh hạn mùa hè tại đây hình thành do đới gió Tây Nam vượt dãy Trường Sơn gây ra hiệu ứng Phơn khô nóng gay gắt kết hợp cùng sự phát triển mạnh mẽ của áp thấp nóng phía tây.

Tại sao cần kết hợp nhiều chỉ số hạn thay vì sử dụng một chỉ số đơn lẻ? Mỗi chỉ số phản ánh một khía cạnh vật lý riêng biệt: chỉ số SPI đánh giá thuần túy lượng giáng thủy, chỉ số J nhạy bén với thời điểm khởi phát hạn, còn chỉ số Ped phản ánh tác động đồng thời của nhiệt độ cao và ít mưa. Việc kết hợp giúp tránh đánh giá sai lệch điều kiện khô hạn.

Kết luận

  • Luận văn đã thiết lập bức tranh toàn diện và có hệ thống về tần suất, mức độ phân hóa hạn hán trong 34 năm (1981-2014) tại 7 vùng khí hậu và 5 trạm đảo tiêu biểu của Việt Nam qua bộ 3 chỉ số J, SPI và Ped.
  • Làm sáng tỏ cơ chế vật lý chi phối hạn hán từ các hoàn lưu khí quyển quy mô lớn, bao gồm chu trình ENSO, áp cao cận nhiệt đới Tây Thái Bình Dương ở mực 850mb, áp thấp nóng phía tây và áp cao lạnh lục địa Siberi.
  • Xác định quy luật phân mùa hạn rõ rệt: Tây Nguyên và Nam Bộ chịu tần suất hạn mùa khô cao nhất (80-100% vào tháng 1-3), trong khi Bắc Trung Bộ ghi nhận đỉnh hạn mùa hè do hiệu ứng Phơn và miền Bắc tập trung hạn trong các tháng chính đông.
  • Chứng minh tính ứng dụng cao của mối tương quan có độ trễ 4 tháng giữa nhiệt độ bề mặt nước biển (SST) vùng Nino 3.4, chỉ số ONI và nhiệt độ đất liền trong công tác thiết lập hệ thống cảnh báo sớm hạn hán.
  • Cung cấp cơ sở khoa học định lượng phục vụ công tác dự tính xu thế khô hạn giai đoạn 2016-2027 theo các kịch bản phát thải RCP 2.6, RCP 4.5 và RCP 8.5 từ mô hình NCAR-CCSM.

Trong định hướng phát triển giai đoạn 2026-2030, các nghiên cứu tiếp theo cần mở rộng tích hợp dữ liệu viễn thám độ ẩm đất tầng sâu và ảnh vệ tinh nhằm nâng cao độ phân giải không gian trong giám sát hạn hán thời gian thực. Các cơ quan quản lý và viện nghiên cứu chuyên ngành được khuyến khích khai thác triệt để những kết quả định lượng của luận văn này nhằm nâng cao năng lực ứng phó chủ động với thiên tai và biến đổi khí hậu tại Việt Nam.