Tổng quan nghiên cứu

Khu vực Nam Trung Bộ là một trong những trọng điểm chịu ảnh hưởng nặng nề nhất của thiên tai khô hạn tại Việt Nam, nơi có lượng mưa trung bình năm vùng đồng bằng chỉ đạt khoảng 720 mm nhưng lượng bốc hơi tiềm năng lên tới 1.860 mm, gấp gần 2,6 lần lượng giáng thủy. Dưới tác động sâu sắc của biến đổi khí hậu toàn cầu, trong 50 năm qua nhiệt độ trung bình trên toàn lãnh thổ Việt Nam đã tăng khoảng 0,5 đến 0,7°C và mực nước biển dâng xấp xỉ 20 cm, làm gia tăng tính bất ổn định của các hiện tượng thời tiết cực đoan. Nghiên cứu này được thực hiện nhằm mục tiêu đánh giá định lượng hiện trạng, tính chất và quy luật diễn biến của hạn hán trong quá khứ giai đoạn 1961 - 2012, đồng thời ứng dụng mô hình khí hậu khu vực độ phân giải cao để dự tính mức độ hạn hán trong các thời kỳ tương lai của thế kỷ 21.

Phạm vi không gian nghiên cứu bao quát 5 tỉnh duyên hải Nam Trung Bộ gồm Bình Định, Phú Yên, Khánh Hòa, Ninh Thuận và Bình Thuận với tổng diện tích tự nhiên trên 21.000 km². Luận văn tập trung giải quyết bài toán cốt lõi về việc định lượng hóa mức độ khô hạn theo nhiều thang thời gian khác nhau (1, 3, 6 và 12 tháng) dựa trên chỉ số chuẩn hóa lượng mưa. Kết quả nghiên cứu mang ý nghĩa thực tiễn to lớn trong việc cung cấp cơ sở khoa học chính xác phục vụ quy hoạch, phân bổ nguồn tài nguyên nước cho hơn 400.000 ha đất sản xuất nông nghiệp, hỗ trợ các nhà quản lý xây dựng phương án ứng phó khẩn cấp và giảm thiểu thiệt hại kinh tế hàng trăm tỷ đồng do hạn hán gây ra mỗi năm.

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 dựa trên khung lý thuyết chuẩn hóa của Tổ chức Khí tượng Thế giới (WMO) và các báo cáo đánh giá của Ban Liên chính phủ về Biến đổi Khí hậu (IPCC AR4, SREX). Khung lý thuyết phân chia hạn hán thành bốn cấp độ diễn tiến: hạn khí tượng (thiếu hụt giáng thủy kéo dài), hạn nông nghiệp (suy giảm độ ẩm trong đất và mất cân bằng bốc thoát hơi), hạn thủy văn (suy kiệt dòng chảy mặt và mực nước ngầm) và hạn kinh tế - xã hội (ảnh hưởng trực tiếp đến đời sống dân sinh, an ninh lương thực).

Mô hình nghiên cứu áp dụng lý thuyết cân bằng ẩm khí quyển - mặt đệm kết hợp với hệ thống các chỉ số hạn chuyên dụng. Trong đó, chỉ số chuẩn hóa lượng mưa (SPI) do McKee đề xuất năm 1993 được lựa chọn làm công cụ hạt nhân nhờ tính linh hoạt trên nhiều quy mô thời gian. Thang phân cấp hạn được xác định chặt chẽ: bắt đầu hạn (-0,5 < SPI ≤ 0,25), hạn vừa (-1,0 < SPI ≤ -0,5), hạn nặng (-1,5 < SPI ≤ -1,0) và hạn rất nặng đến đặc biệt nghiêm trọng (SPI ≤ -1,5). Nghiên cứu tích hợp lý thuyết biến đổi khí hậu theo kịch bản phát thải trung bình A1B để đánh giá mối liên hệ giữa sự nóng lên toàn cầu và chu kỳ dao động khí hậu El Nino - Dao động Nam (ENSO).

   [Phát thải KNK & Biến đổi khí hậu]
   [Mô hình Khí hậu Toàn cầu HadCM3]
   [Mô hình Khí hậu Khu vực PRECIS (25x25 km)]
   [Chuỗi số liệu mưa mô phỏng (A1B: 2020 - 2099)]
      [Tính toán Chỉ số SPI]          [Số liệu Quan trắc 1961-2012]
      (1, 3, 6, 12 tháng)             (8 trạm khí tượng tiêu biểu)
                     [Phân tích Xu thế Tuyến tính]
                     [Bản đồ & Kịch bản Hạn hán]

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

Dữ liệu quan trắc lịch sử được thu thập liên tục trong 52 năm (1961 - 2012) tại 8 trạm khí tượng tiêu chuẩn phủ đều khu vực nghiên cứu gồm: Quy Nhơn, Sơn Hòa, Tuy Hòa, Nha Trang, Cam Ranh, Phan Rang, Phan Thiết và Hàm Tân. Cỡ mẫu 8 trạm với chuỗi thời gian hơn 50 năm đáp ứng hoàn toàn yêu cầu về tính đồng nhất không gian và độ tin cậy thống kê khí hậu. Phương pháp chọn mẫu tập trung vào các điểm nút đại diện cho cả vùng đồng bằng ven biển và vùng đồi núi thấp chuyển tiếp, phản ánh trọn vẹn sự phân hóa địa hình từ vĩ độ 10°41'B đến 13°46'B.

Phương pháp phân tích xu thế tuyến tính bằng kỹ thuật bình phương tối thiểu được sử dụng để xác định hệ số góc biến đổi lượng mưa và chỉ số SPI theo thời gian. Để dự tính tương lai, luận văn sử dụng mô hình động lực khí hậu khu vực PRECIS do Trung tâm Hadley phát triển, chạy với độ phân giải ngang 25x25 km, cấu trúc 19 tầng khí áp theo hệ tọa độ lai. Nghiên cứu thực hiện mô phỏng tổ hợp từ 5 phương án thành phần (Q0, Q3, Q10, Q11, Q13) thuộc kịch bản phát thải trung bình A1B cho 4 thời kỳ tương lai: 2020 - 2039, 2040 - 2059, 2060 - 2079 và 2080 - 2099 so với thời kỳ nền 1980 - 1999.

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

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

Thứ nhất, phân tích chuỗi số liệu 1961 - 2012 cho thấy có sự phân hóa rõ rệt về xu thế hạn hán theo không gian. Tại 6/7 trạm quan trắc (Tuy Hòa, Sơn Hòa, Nha Trang, Cam Ranh, Phan Rang, Phan Thiết), chỉ số SPI có xu thế tuyến tính tăng nhẹ với hệ số góc dương từ 0,0112 đến 0,0567, phản ánh lượng mưa năm có xu hướng gia tăng nhẹ từ 5 đến 20% trong nửa thế kỷ. Ngược lại, trạm Hàm Tân (Bình Thuận) là điểm duy nhất thể hiện xu thế giảm tuyến tính của chỉ số SPI với hệ số b1 = -0,0181, báo hiệu nguy cơ khô hạn gia tăng cục bộ tại khu vực cực Nam.

Thứ hai, mức độ khắc nghiệt của các đợt hạn nặng có xu hướng gia tăng mạnh mẽ và gắn liền mật thiết với chu kỳ hoạt động của hiện tượng El Nino. Trong các năm 1983, 1997 - 1998, 2004 - 2005, chỉ số SPI tại hầu hết các trạm đều chạm ngưỡng cực tiểu từ -1,5 đến -2,5 (hạn rất nặng đến đặc biệt nghiêm trọng). Điển hình trong đợt hạn 2004 - 2005 tại Ninh Thuận, lượng mưa 4 tháng mùa khô chỉ đạt 41% mức trung bình nhiều năm, khiến dung tích trữ của các hồ chứa giảm hơn 50% so với thiết kế và hồ Tân Giang xuống dưới mực nước chết.

Thứ ba, phân bố hạn theo thời gian thể hiện tính mùa sâu sắc. Mùa khô kéo dài từ tháng 1 đến tháng 8 ghi nhận tần suất xuất hiện hạn đặc biệt cao, trong đó các tháng 1 đến 4 có tần suất hạn đạt tới 100% tại Ninh Thuận và Bình Thuận. Trái lại, trong 3 tháng mùa mưa chính vụ (tháng 9 đến tháng 11), tần suất hạn tại các tỉnh phía Bắc gần như bằng 0%.

Thứ tư, kết quả dự tính khí hậu tương lai từ mô hình PRECIS theo kịch bản A1B chỉ ra rằng số tháng hạn quy mô 1 tháng, 3 tháng và 6 tháng tiếp tục duy trì ở mức cao trong hai thập kỷ 2020 - 2039 và 2040 - 2059. Các đợt hạn kéo dài liên tục từ 6 đến 12 tháng có xu hướng gia tăng về tần suất xuất hiện, đặc biệt tại khu vực cực Nam Trung Bộ.

Tần suất và Mức độ Hạn hán tại các Trạm Khí tượng Đại diện (1961 - 2012)
Trạm Khí tượng     Tỉnh/Thành    Xu thế SPI (b1)   SPI Cực tiểu (Năm El Nino)
Tuy Hòa            Phú Yên       +0,0315 (Giảm)    -2,00 (Năm 2005)
Sơn Hòa            Phú Yên       +0,0112 (Giảm)    -1,50 (Năm 1983)
Nha Trang          Khánh Hòa     +0,0425 (Giảm)    -1,80 (Năm 1967)
Cam Ranh           Khánh Hòa     +0,0505 (Giảm)    -1,50 (Năm 2005)
Phan Rang          Ninh Thuận    +0,0567 (Giảm)    -2,50 (Năm 1988)
Phan Thiết         Bình Thuận    +0,0137 (Giảm)    -2,00 (Năm 1977)
Hàm Tân            Bình Thuận    -0,0181 (Tăng)    -2,50 (Năm 2005)

Thảo luận kết quả

Diễn biến hạn hán phức tạp tại Nam Trung Bộ bắt nguồn từ sự tương tác giữa địa hình khuất gió và hoàn lưu khí quyển. Các dãy núi cao từ 1.200 m đến 2.000 m ở phía Tây tạo bức tường chắn gió mùa Tây Nam, gây hiệu ứng phơn khô nóng gay gắt trên diện rộng. Kết quả này tương đồng với các công bố của Viện Khoa học Khí tượng Thủy văn và Biến đổi khí hậu (IMHEN, 2010), khẳng định rằng dù tổng lượng mưa năm có xu thế tăng nhẹ, nhưng sự phân bố mưa ngày càng cực đoan: mưa tập trung cường suất lớn trong thời gian ngắn gây lũ lụt, xen kẽ với các đợt khô hạn kéo dài khốc liệt.

Các dữ liệu nghiên cứu được biểu diễn trực quan qua hệ thống đồ thị đường thể hiện biến thiên chỉ số SPI hàng năm và bản đồ phân bố không gian theo ô lưới 25x25 km của mô hình PRECIS. Các biểu đồ cột phân cấp mức độ hạn sử dụng dải màu trực quan (cam cho bắt đầu hạn, đỏ nhạt cho hạn vừa, đỏ sẫm cho hạn nặng và xanh cho điều kiện ẩm ướt) giúp nhận diện rõ nét tính chu kỳ của các đợt hạn lịch sử và các giai đoạn tương lai.

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

Dựa trên các kết quả tính toán định lượng, luận văn đề xuất 4 nhóm giải pháp chiến lược nhằm thích ứng hiệu quả với biến đổi khí hậu và giảm thiểu thiệt hại do hạn hán:

                          CÁC NHÓM GIẢI PHÁP ỨNG PHÓ HẠN HÁN
[Công trình Thủy lợi] [Chuyển đổi Nông nghiệp] [Cảnh báo Sớm]     [Phục hồi Rừng]
- Nâng dung tích hồ    - Tưới tiết kiệm 50.000 ha - Dự báo SPI đa kỳ  - Độ che phủ > 55%
- Kết nối liên lưu vực - Giảm 30 - 40% nước tưới - Độ tin cậy > 85%   - Bảo vệ nguồn sinh thủy
(2025 - 2030)          (2024 - 2028)              (2025 - 2026)       (2025 - 2030)
  1. Nâng cấp và hoàn thiện hệ thống công trình thủy lợi liên lưu vực: Đẩy mạnh đầu tư xây dựng các tuyến kênh dẫn liên tỉnh, nâng cấp dung tích trữ của các hồ chứa lớn thêm 20 đến 30% tại Ninh Thuận và Bình Thuận. Chủ thể thực hiện là Bộ Nông nghiệp và Phát triển Nông thôn phối hợp với UBND các tỉnh, hoàn thành giai đoạn 2025 - 2030.
  2. Chuyển đổi cơ cấu cây trồng và ứng dụng công nghệ tưới tiết kiệm: Triển khai chuyển đổi hơn 50.000 ha đất trồng lúa kém hiệu quả sang các loại cây trồng cạn có khả năng chịu hạn cao (ngô lai, cây ăn quả, nho, táo); áp dụng công nghệ tưới nhỏ giọt và tưới phun mưa tự động giúp cắt giảm 30 đến 40% lượng nước tiêu thụ. Đơn vị chủ trì là Sở Nông nghiệp và Phát triển Nông thôn các địa phương, thực hiện trong giai đoạn 2024 - 2028.
  3. Hiện đại hóa công nghệ giám sát và cảnh báo sớm hạn hán: Ứng dụng hệ thống dự báo số mục hạn dựa trên chỉ số SPI tích hợp dữ liệu viễn thám và mô hình động lực độ phân giải cao, nâng cao độ chính xác cảnh báo trước 3 đến 6 tháng đạt trên 85%. Chủ thể thực hiện là Trung tâm Dự báo Khí tượng Thủy văn Quốc gia, hoàn thành trước năm 2026.
  4. Bảo vệ và phục hồi rừng phòng hộ đầu nguồn: Đẩy mạnh trồng rừng giữ nước tại các lưu vực sông Ba, sông Cái và sông La Ngà, phấn đấu nâng tỷ lệ che phủ rừng toàn vùng lên trên 55% nhằm tăng cường khả năng giữ ẩm tự nhiên và chống hoang mạc hóa đất đai. Chi cục Kiểm lâm các tỉnh phối hợp cùng cộng đồng địa phương triển khai liên tục từ năm 2025 đến 2030.

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

Luận văn cung cấp nguồn tư liệu học thuật và thực tiễn giá trị cho 4 nhóm đối tượng chính:

  • Cơ quan quản lý nhà nước và nhà hoạch định chính sách: Các chuyên viên tại Sở Tài nguyên và Môi trường, Sở Nông nghiệp và Phát triển Nông thôn sử dụng các bản đồ phân vùng rủi ro hạn hán và số liệu dự tính mô hình PRECIS để xây dựng kế hoạch phân bổ nguồn nước và chiến lược thích ứng biến đổi khí hậu cấp tỉnh.
  • Các nhà nghiên cứu và viện chuyên ngành: Các nhà khoa học khí tượng thủy văn, địa lý tài nguyên khai thác phương pháp luận tính toán chỉ số SPI đa quy mô và quy trình tổ hợp mô hình khí hậu khu vực để mở rộng cho các lưu vực sông khác trên cả nước.
  • Doanh nghiệp nông nghiệp và hợp tác xã địa phương: Ban lãnh đạo các hợp tác xã nông nghiệp tham khảo thông tin chu kỳ mùa khô và quy luật xuất hiện hạn 1 - 3 tháng để bố trí lịch thời vụ chính xác, chủ động đầu tư hệ thống tích trữ nước và kỹ thuật tưới phù hợp.
  • Giảng viên, nghiên cứu sinh và sinh viên đại học: Giảng viên và học viên cao học các ngành Khí tượng, Thủy văn, Biến đổi Khí hậu và Quản lý Tài nguyên Nước sử dụng tài liệu làm dẫn chứng học thuật và khung tham chiếu phương pháp nghiên cứu xử lý chuỗi số liệu khí hậu dài hạn.

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

Chỉ số SPI có ưu điểm gì vượt trội trong đánh giá hạn hán so với các chỉ số truyền thống?

Chỉ số SPI chỉ yêu cầu đầu vào duy nhất là chuỗi số liệu lượng mưa liên tục, cho phép tính toán linh hoạt trên nhiều thang thời gian khác nhau (1, 3, 6, 12, 24 tháng). Điều này giúp đánh giá chính xác cả hạn khí tượng ngắn hạn lẫn hạn nông nghiệp và thủy văn dài hạn mà không bị phụ thuộc vào sự phức tạp của các thông số độ ẩm đất hay bức xạ.

Vì sao lượng mưa năm tại Nam Trung Bộ có xu thế tăng nhưng hạn hán mùa khô vẫn diễn biến khốc liệt?

Mặc dù lượng mưa năm tăng từ 0 đến 40% tại một số trạm, nhưng lượng mưa phân bố không đều và tập trung tới 85 - 90% vào các tháng mùa mưa ngắn ngủi (tháng 9 đến tháng 11). Trong mùa khô kéo dài 6 đến 8 tháng, lượng mưa thiếu hụt nghiêm trọng kết hợp với nhiệt độ cao làm bốc hơi tới 1.860 mm/năm khiến khô hạn thêm trầm trọng.

Mô hình PRECIS đóng vai trò gì trong việc dự tính kịch bản hạn hán tương lai?

PRECIS là mô hình khí hậu khu vực có độ phân giải cao (25x25 km), mô tả chi tiết hiệu ứng địa hình phức tạp và tương tác khí quyển tại Nam Trung Bộ. Luận văn sử dụng phương pháp tổ hợp 5 phương án khác nhau của PRECIS dưới kịch bản A1B, giúp giảm thiểu độ bất định và tăng tính tin cậy cho các dự báo hạn hán đến năm 2099.

Hiện tượng El Nino tác động cụ thể như thế nào đến tình hình hạn hán tại Nam Trung Bộ?

Trong những năm El Nino hoạt động mạnh (như 1982 - 1983, 1997 - 1998, 2004 - 2005), lượng mưa mùa khô giảm từ 30 đến 60% so với trung bình nhiều năm, bức xạ mặt trời và nhiệt độ tăng từ 0,5 đến 1,5°C. Điều này thúc đẩy tốc độ bốc hơi nước cực nhanh, đẩy chỉ số SPI xuống mức âm sâu dưới -2,0, gây ra hạn hán trên diện rộng.

Các tỉnh nào tại khu vực Nam Trung Bộ có nguy cơ hạn hán cao nhất trong tương lai?

Ninh Thuận và Bình Thuận là hai tỉnh chịu rủi ro hạn hán nghiêm trọng nhất toàn vùng. Kết quả phân tích cho thấy trạm Hàm Tân (Bình Thuận) có xu thế hạn tăng liên tục (b1 = -0,0181), đồng thời tần suất hạn các tháng mùa khô tại khu vực này thường xuyên chạm ngưỡng 100%, đòi hỏi các giải pháp can thiệp khẩn cấp.

Kết luận

  • Luận văn đã xác lập bức tranh toàn diện về đặc trưng hạn hán Nam Trung Bộ giai đoạn 1961 - 2012 qua chỉ số SPI, khẳng định tính bất thường của hạn hán gắn liền với các đợt El Nino lịch sử.
  • Nghiên cứu phát hiện sự đối lập không gian rõ rệt: trong khi phần lớn khu vực có xu thế tăng nhẹ chỉ số SPI thì vùng cực Nam (Hàm Tân, Bình Thuận) lại có xu thế gia tăng mức độ khô hạn.
  • Ứng dụng thành công mô hình khí hậu khu vực PRECIS độ phân giải 25x25 km để mô phỏng và dự tính chi tiết các cấp độ hạn 1, 3, 6, 12 tháng theo kịch bản phát thải trung bình A1B đến năm 2099.
  • Đóng góp cơ sở dữ liệu khoa học chuẩn xác phục vụ trực tiếp công tác điều hành thủy lợi, chuyển đổi cơ cấu sản xuất nông nghiệp và quy hoạch tài nguyên nước bền vững cho 5 tỉnh duyên hải miền Trung.
  • Trong giai đoạn tiếp theo, các địa phương cần đẩy mạnh triển khai đồng bộ 4 nhóm giải pháp công trình và phi công trình trước năm 2030 nhằm bảo đảm an ninh nguồn nước và phát triển kinh tế - xã hội vững chắc.