Tổng quan nghiên cứu

Hạn hán là loại hình thiên tai nghiêm trọng thứ 3 tại Việt Nam, chỉ sau bão và lũ lụt, gây ra những tổn thất nặng nề cho kinh tế - xã hội và tài nguyên môi trường. Thực tế lịch sử ghi nhận đợt hạn hán khốc liệt năm 1997-1998 do El Nino kéo dài 15 tháng đã khiến hơn 843.000 ha lúa bị thiệt hại và gần 3.000.000 người dân rơi vào tình trạng thiếu nước sinh hoạt trầm trọng. Dưới tác động của biến đổi khí hậu toàn cầu, tính chất và quy luật của hạn hán ngày càng diễn biến phức tạp, đòi hỏi các công cụ mô phỏng động lực hiện đại thay vì chỉ dựa vào chuỗi số liệu quan trắc quá khứ.

Vấn đề nghiên cứu trọng tâm của luận văn là đánh giá và dự tính sự biến đổi hạn hán tại 7 vùng khí hậu đặc trưng của Việt Nam thông qua mô hình khí hậu khu vực. Mục tiêu cụ thể là kiểm định năng lực mô phỏng khí hậu của mô hình RegCM3 trong thời kỳ chuẩn 1970-1999 (30 năm), từ đó đưa ra các kịch bản dự tính hạn hán chi tiết cho giai đoạn 2011-2050 (40 năm) theo hai kịch bản phát thải khí nhà kính A1B và A2. Phạm vi không gian bao phủ toàn bộ lãnh thổ Việt Nam từ Tây Bắc đến Nam Bộ với độ phân giải mô hình 36 km. Kết quả nghiên cứu mang ý nghĩa thực tiễn to lớn trong việc thiết lập cơ sở dữ liệu định lượng, giúp các cơ quan quản lý chủ động xây dựng phương án phòng ngừa rủi ro thiên tai và điều tiết tài nguyên nước bền vững.

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 xây dựng trên nền tảng lý thuyết hạn khí tượng học và phương pháp mô hình hóa động lực khí hậu khu vực theo báo cáo đánh giá của Ban Liên chính phủ về Biến đổi Khí hậu. Bốn khái niệm học thuật cốt lõi được ứng dụng xuyên suốt bao gồm:

  • Hạn khí tượng: Trạng thái thiếu hụt lượng giáng thủy tự nhiên trong một khoảng thời gian kéo dài so với giá trị trung bình nhiều năm của cùng thời kỳ.
  • Chỉ số khô cằn De Martonne (J): Được xác định theo biểu thức tương quan giữa lượng mưa tháng và nhiệt độ trung bình tháng, phục vụ việc phân loại và nhận diện các tháng hạn hán trong năm.
  • Chỉ số Ped: Chỉ số đánh giá dị thường khí hậu dựa trên độ lệch chuẩn hóa của nhiệt độ và lượng giáng thủy, cho phép định lượng mức độ khắc nghiệt và xu thế hạn hán trên quy mô năm.
  • Kịch bản phát thải khí nhà kính SRES: Khung kịch bản phát thải toàn cầu, trong đó kịch bản A1B đại diện cho xu hướng cân bằng công nghệ - năng lượng và kịch bản A2 đại diện cho xu hướng phát thải cao với tốc độ tăng dân số nhanh.

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

Nghiên cứu thu thập và xử lý bộ số liệu quan trắc khí tượng bề mặt gồm nhiệt độ trung bình tháng và lượng mưa tích lũy tháng tại 50 trạm quan trắc tiêu chuẩn thuộc 7 vùng khí hậu Việt Nam trong thời kỳ chuẩn 1970-1999. Cỡ mẫu 50 trạm được lựa chọn theo phương pháp phân tầng địa lý đại diện, bao phủ đầy đủ các điều kiện địa hình phức tạp từ vùng núi cao Tây Bắc đến đồng bằng ven biển Nam Bộ, đảm bảo tính liên tục và độ tin cậy chuẩn hóa của chuỗi dữ liệu 30 năm.

Phương pháp phân tích chủ đạo là hạ quy mô động lực bằng mô hình khí hậu khu vực RegCM3 của Trung tâm Quốc gia Nghiên cứu Khí quyển Hoa Kỳ, lồng ghép điều kiện biên từ mô hình toàn cầu CCSM3.0. Miền tính được thiết lập trên lưới 144x105 nút với độ phân giải ngang 36 km, sử dụng sơ đồ tham số hóa đối lưu Grell-Arakawa và mô hình bề mặt BATS. Để đánh giá sai số mô hình so với quan trắc, ba chỉ số thống kê chuẩn được áp dụng là sai số trung bình (ME), sai số trung bình tuyệt đối (MAE) và sai số quân phương (RMSE). Tiến trình nghiên cứu phân tách rõ rệt giữa giai đoạn kiểm định mô hình (1970-1999) và giai đoạn dự tính tương lai 2011-2050, được chia thành hai thời đoạn đánh giá 2011-2030 và 2031-2050.

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

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

Quá trình kiểm định và mô phỏng số liệu cho thấy bốn phát hiện khoa học quan trọng:

  • Khả năng mô phỏng nhiệt độ và lượng mưa của RegCM3: Mô hình có xu thế mô phỏng nhiệt độ thiên thấp (ME âm) trong mùa đông và đầu mùa xuân trên các vùng phía Bắc. Sai số ME âm lớn nhất ghi nhận tại Tây Bắc vào tháng 3 đạt -5,8°C với MAE là 5,89°C; sai số RMSE cao nhất xuất hiện tại Bắc Trung Bộ vào tháng 3 với giá trị 6,67°C. Lượng mưa mô phỏng có xu hướng thấp hơn quan trắc trong các tháng mùa mưa, với ME âm từ -250 mm đến -316,07 mm tại Bắc Trung Bộ vào tháng 10.
  • Đặc trưng mùa hạn theo chỉ số J: Kết quả mô phỏng xác định chính xác thời kỳ khô hạn tại vùng Tây Bắc từ tháng 10 đến tháng 4 với giá trị J nhỏ nhất là 6,68 mm/°C vào tháng 11 (quan trắc là 9,47 mm/°C vào tháng 12); vùng Đông Bắc bị hạn từ tháng 10 đến tháng 3 với J nhỏ nhất đạt 8,54 mm/°C. Tại Nam Bộ, hạn hán khốc liệt tập trung từ tháng 12 đến tháng 4 năm sau với J quan trắc nhỏ nhất đạt 1,76 mm/°C vào tháng 2.
  • Xu thế biến đổi hạn theo chỉ số Ped: Trong thời kỳ chuẩn 1970-1999, chỉ số Ped quan trắc tại 4 vùng khí hậu phía Bắc đều thể hiện xu thế tăng tuyến tính, phản ánh tần suất hạn hán gia tăng. Mức độ trùng khớp giữa mô hình và thực tế đạt 57% tại Bắc Trung Bộ (14/30 năm khớp hoàn toàn giữa các pha hạn và ẩm) và 47% tại Tây Bắc. Các năm hạn nặng tính toán được đều trùng khớp với các năm xuất hiện hiện tượng El Nino mạnh như 1987, 1991, 1998.
  • Dự tính hạn hán tương lai (2011-2050): Dưới tác động của kịch bản A1B và A2, nền nhiệt độ trung bình toàn quốc dự tính tăng từ 1,2°C đến 2,5°C so với thời kỳ chuẩn. Sự suy giảm lượng mưa mùa khô từ 10% đến 25% kết hợp bốc hơi mạnh làm tăng đáng kể tần suất và cường độ hạn hán tại Tây Nguyên và Nam Trung Bộ.

Thảo luận kết quả

Sự sai lệch giữa mô phỏng của RegCM3 và quan trắc xuất phát từ địa hình chia cắt phức tạp của dãy Hoàng Liên Sơn và dãy Trường Sơn, nơi các quá trình hoàn lưu vi khí hậu quy mô nhỏ chưa được mô hình hóa toàn diện ở độ phân giải 36 km. Ngoài ra, mô hình xuất hiện hiện tượng trễ pha mùa mưa từ 1 đến 2 tháng tại khu vực duyên hải miền Trung, dẫn đến việc mô phỏng đỉnh hạn trễ hơn so với thực tế.

Khi so sánh với các công trình quốc tế của Sheffield & Wood (2008) và Kenneth Strzepek (2010), xu thế gia tăng hạn hán theo kịch bản A2 hoàn toàn tương thích với các dự báo toàn cầu về sự suy giảm độ ẩm đất và gia tăng các đợt khô hạn cực đoan. Về mặt ứng dụng trực quan, toàn bộ kết quả có thể được thể hiện hiệu quả qua biểu đồ chuỗi thời gian so sánh đường xu thế tuyến tính của chỉ số Ped, kết hợp bảng ma trận phân loại chỉ số J theo 12 tháng, tạo điều kiện thuận lợi cho các cơ quan địa phương theo dõi ngưỡng báo động hạn.

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

  • Nâng cấp hệ thống quan trắc và tích hợp chỉ số cảnh báo sớm: Cục Khí tượng Thủy văn cần chủ trì tích hợp bộ chỉ số J và Ped vào quy trình dự báo nghiệp vụ trước năm 2026, nâng độ phân giải mô hình cảnh báo xuống mức 10-15 km nhằm giảm sai số RMSE nhiệt độ xuống dưới 1,5°C và cảnh báo sớm hạn hán trước từ 30 đến 60 ngày.
  • Chuyển đổi cơ cấu cây trồng và mùa vụ nông nghiệp thích ứng: Bộ Nông nghiệp và Phát triển Nông thôn phối hợp cùng các tỉnh Nam Trung Bộ và Tây Nguyên chuyển đổi ít nhất 40.000 ha đất lúa kém hiệu quả sang cây trồng cạn chịu hạn trước năm 2028, giúp tiết kiệm từ 25% đến 35% lượng nước tưới tiêu trong các tháng mùa cạn.
  • Hoàn thiện quy trình vận hành liên hồ chứa đa mục tiêu: Tập đoàn Điện lực Việt Nam và các công ty quản lý khai thác công trình thủy lợi cần phê duyệt quy chế phối hợp xả nước mùa kiệt trước năm 2025, duy trì dung tích trữ hữu ích trên 80% tại các hồ chứa trọng điểm trong các tháng cao điểm khô hạn từ tháng 1 đến tháng 4.
  • Xây dựng bản đồ phân vùng rủi ro hạn hán cấp huyện: Ủy ban Nhân dân các tỉnh, thành phố hoàn thành việc tích hợp bản đồ nguy cơ hạn hán chi tiết theo kịch bản A1B và A2 vào quy hoạch phát triển kinh tế - xã hội giai đoạn 2025-2030, hướng đến mục tiêu giảm thiểu 50% thiệt hại kinh tế do khô hạn gây ra hàng năm.

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

  • Cơ quan quản lý nhà nước về khí tượng thủy văn và biến đổi khí hậu: Sử dụng bộ phương pháp đánh giá sai số mô hình và hệ thống chỉ số Ped, J để chuẩn hóa công tác dự báo hạn hán hạn dài trên 7 vùng khí hậu.
  • Đơn vị quản lý tài nguyên nước, thủy lợi và phòng chống thiên tai: Khai thác chi tiết khung thời gian xuất hiện hạn hán theo tháng để lập kế hoạch tích trữ nước và xây dựng phương án cấp nước sinh hoạt khẩn cấp cho các vùng khô hạn.
  • Viện nghiên cứu nông nghiệp và nhà hoạch định quy hoạch nông thôn: Ứng dụng số liệu dự tính lượng mưa và nhiệt độ giai đoạn 2011-2050 để điều chỉnh lịch thời vụ, lựa chọn bộ giống cây trồng thích ứng tại các điểm nóng hạn hán miền Trung và Tây Nguyên.
  • Giảng viên, nghiên cứu sinh và học viên cao học chuyên ngành Khí tượng, Biến đổi Khí hậu: Tham khảo quy trình kỹ thuật hạ quy mô động lực bằng mô hình RegCM3, phương pháp ghép nối số liệu toàn cầu CCSM3.0 và kỹ thuật xử lý dữ liệu khí hậu đa biến.

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

  • Tại sao nghiên cứu lại lựa chọn hai chỉ số Ped và J để phân tích hạn hán tại Việt Nam? Chỉ số J của De Martonne phản ánh trực quan mức độ khô cằn theo tháng qua tỷ số giữa lượng mưa và nhiệt độ. Chỉ số Ped chuẩn hóa độ lệch của cả hai yếu tố trên thang năm, giúp định lượng chính xác các năm hạn khắc nghiệt và xu thế biến đổi dài hạn dựa trên chuỗi số liệu 50 trạm khí tượng.

  • Mô hình khí hậu khu vực RegCM3 có ưu thế gì so với mô hình toàn cầu trong nghiên cứu này? RegCM3 hạ quy mô không gian xuống lưới phân giải 36 km với 144x105 nút, mô tả chi tiết hơn đặc điểm địa hình đồi núi phức tạp của Việt Nam so với lưới thô hàng trăm km của mô hình toàn cầu CCSM3.0, từ đó mô phỏng hoàn lưu gió mùa và mưa cục bộ chính xác hơn.

  • Đợt hạn hán lịch sử năm 1997-1998 đã gây ra những tổn thất cụ thể nào tại Việt Nam? Do tác động của hiện tượng El Nino kéo dài 15 tháng, đợt hạn 1997-1998 đã làm mất trắng 254.740 ha lúa vụ đông xuân, ảnh hưởng 435.320 ha lúa hè thu và khiến gần 3.000.000 người dân trên toàn quốc lâm vào cảnh thiếu nước sinh hoạt.

  • Nguyên nhân nào dẫn đến sự sai số giữa kết quả mô phỏng của RegCM3 và số liệu quan trắc? Sai số xuất phát từ ảnh hưởng của địa hình chia cắt mạnh và sự biến động phức tạp của gió mùa Đông Bắc. Mô hình có xu hướng mô phỏng nhiệt độ thiên thấp trong mùa đông (ME đạt -5,8°C tại Tây Bắc) và đánh giá lượng mưa mùa mưa thấp hơn thực tế khoảng 200-300 mm.

  • Kịch bản phát thải A1B và A2 đưa ra những dự báo cốt lõi nào về hạn hán giai đoạn 2011-2050? Cả hai kịch bản đều chỉ ra xu thế tăng nhiệt độ từ 1,2°C đến 2,5°C trên cả nước. Lượng mưa trong các tháng mùa khô giảm từ 10% đến 25% tại Nam Trung Bộ và Tây Nguyên, làm gia tăng tần suất các đợt hạn nặng và kéo dài thời gian khô hạn.

Kết luận

  • Tổng kết toàn diện đặc trưng hạn hán 30 năm (1970-1999) trên 7 vùng khí hậu qua chuỗi dữ liệu 50 trạm quan trắc tiêu chuẩn.
  • Kiểm chứng thành công khả năng ứng dụng của mô hình khu vực RegCM3 độ phân giải 36 km kết hợp chỉ số J và Ped trong nghiên cứu hạn hán.
  • Dự tính bức tranh hạn hán chi tiết giai đoạn 2011-2050 theo kịch bản A1B và A2, khẳng định nguy cơ gia tăng khô hạn tại miền Trung và Tây Nguyên.
  • Cung cấp luận cứ khoa học định lượng vững chắc giúp các bộ, ngành hoàn thiện chiến lược thích ứng biến đổi khí hậu trước năm 2030.
  • Khuyến nghị các nghiên cứu tiếp theo mở rộng phương pháp mô hình hóa tổ hợp với độ phân giải cao dưới 10 km để tối ưu hóa dự báo cục bộ.