Tổng quan nghiên cứu
Biến đổi khí hậu đang là thách thức an ninh phi truyền thống nghiêm trọng hàng đầu đối với toàn cầu và Việt Nam. Theo số liệu quan trắc trong giai đoạn 50 năm (1958-2007), nhiệt độ trung bình năm tại Việt Nam đã gia tăng khoảng 0,5°C đến 0,7°C, trong đó tốc độ tăng nhiệt của mùa đông diễn ra nhanh hơn mùa hè và các khu vực phía Bắc có mức tăng nhiệt độ cao hơn đáng kể so với phía Nam. Điển hình tại 3 đô thị lớn gồm Hà Nội, Đà Nẵng và Thành phố Hồ Chí Minh, nhiệt độ trung bình thập kỷ 1991-2000 so với giai đoạn 1931-1940 đã tăng lần lượt là 0,8°C; 0,4°C và 0,6°C. Bước sang năm 2007, mức chênh lệch nhiệt độ trung bình tại cả ba trạm này tiếp tục tăng cao hơn từ 0,8°C đến 1,3°C so với nền khí hậu 1931-1940. Song song với xu thế nhiệt, chế độ mưa và các hiện tượng thời tiết cực đoan như bão, hạn hán, lũ lụt cũng ghi nhận những biến động bất thường và phân hóa sâu sắc trên phạm vi cả nước.
Thực tiễn nghiên cứu khí hậu đặt ra yêu cầu cấp thiết về việc xây dựng các công cụ định lượng tổng hợp thay vì chỉ đánh giá riêng lẻ từng thông số nhiệt độ hay lượng mưa. Luận văn Thạc sĩ Khoa học chuyên ngành Khí tượng và Khí hậu học (Mã số: 60 44 87) của tác giả Trịnh Hoàng Dương, dưới sự hướng dẫn khoa học của TS. Nguyễn Thị Hiền Thuận tại Trường Đại học Khoa học Tự nhiên - Đại học Quốc gia Hà Nội, tập trung giải quyết bài toán: Đánh giá định lượng biến đổi khí hậu ở Việt Nam thông qua việc ứng dụng hai chỉ số tích hợp hiện đại là Chỉ số biến đổi khí hậu khu vực (Regional Climate Change Index - RCCI) và Chỉ số biến đổi khí hậu (Climate Change Index - CCI).
Mục tiêu cốt lõi của đề tài là khảo sát, phân tích chuỗi số liệu quan trắc lịch sử trong giai đoạn 1961-2000 trên 7 vùng khí hậu đặc trưng của Việt Nam và ước tính mức độ biến đổi khí hậu tương lai cho các giai đoạn 2041-2070 và 2071-2100 dựa trên mô hình khí hậu khu vực PRECIS độ phân giải cao 25 km theo các kịch bản phát thải của Ban Liên chính phủ về Biến đổi Khí hậu (IPCC). Nghiên cứu mang lại ý nghĩa khoa học và thực tiễn to lớn, cung cấp cơ sở dữ liệu định lượng phục vụ công tác hoạch định chính sách, quy hoạch phát triển kinh tế - xã hội và xây dựng các chiến lược thích ứng với rủi ro khí hậu cho từng vùng lãnh thổ.
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 vận dụng khung lý thuyết hiện đại về biến đổi khí hậu toàn cầu được thiết lập bởi IPCC và các lý thuyết định lượng điểm nóng khí hậu khu vực:
- Lý thuyết điểm nóng khí hậu và Chỉ số RCCI của Filippo Giorgi (2006): Lý thuyết này xác định các vùng lãnh thổ có phản ứng nhạy cảm và mãnh liệt nhất trước tác động của biến đổi khí hậu toàn cầu thông qua việc tích hợp đồng thời sự biến đổi của giá trị trung bình và sự biến đổi trong dao động liên năm của nhiệt độ và lượng mưa.
- Lý thuyết đánh giá hiện tượng cực đoan và Chỉ số CCI của Michele B. Baettig và cộng sự (2007): Khung tiếp cận này tập trung vào tần suất xuất hiện của các năm cực đoan (năm nóng nhất, năm khô nhất, năm ẩm nhất) vượt qua các ngưỡng phân vị thống kê lịch sử (phân vị 95% cho nhiệt độ và lượng mưa cực đại, phân vị 5% cho lượng mưa cực tiểu).
- Khung lý thuyết phơi nhiễm khí hậu - xã hội (Socio-climatic Exposure) của Noah S. Diffenbaugh và cộng sự (2007, 2008): Tích hợp chỉ số biến đổi khí hậu vật lý với các chỉ số phát triển con người (HDI), mật độ dân số và chỉ số nghèo đói để xác định mức độ tổn thương tổng hợp.
Các khái niệm chuyên ngành cốt lõi được chuẩn hóa trong nghiên cứu bao gồm:
- Biến đổi khí hậu (Climate Change): Sự biến đổi trạng thái trung bình dài hạn của các yếu tố khí hậu qua các thập kỷ.
- Biến động khí hậu (Climate Variability): Mức độ dao động tự nhiên quanh giá trị trung bình giữa các năm, được đo lường bằng độ lệch chuẩn ($\Delta\sigma T$) đối với nhiệt độ và hệ số biến động ($\Delta\sigma P$) đối với lượng mưa.
- Thời kỳ tham chiếu (Baseline Period): Giai đoạn khí hậu nền chuẩn (1961-1980 hoặc 1961-1990) làm mốc so sánh.
- Phân vị cực đoan (Percentile Extremes): Ngưỡng xác suất 5% và 95% nhằm nhận diện các dị thường thời tiết khắc nghiệt.
Phương pháp nghiên cứu
Nguồn dữ liệu và cỡ mẫu nghiên cứu được thu thập toàn diện từ mạng lưới trạm quan trắc bề mặt quốc gia và dữ liệu mô phỏng số trị:
- Dữ liệu quan trắc lịch sử: Thu thập từ chuỗi trạm khí tượng bề mặt tiêu chuẩn đại diện cho 7 vùng khí hậu: Tây Bắc (B1), Đông Bắc (B2), Đồng bằng Bắc Bộ (B3), Trung Trung Bộ (B4), Nam Trung Bộ (N1), Tây Nguyên (N2) và Nam Bộ (N3) trong giai đoạn 40 năm từ 1961 đến 2000. Cỡ mẫu bao gồm đầy đủ dữ liệu nhiệt độ trung bình, nhiệt độ tối cao, nhiệt độ tối thấp và lượng mưa ngày, tháng, mùa. Dữ liệu đã trải qua kiểm định sai số thô và kiểm tra tính đồng nhất nghiêm ngặt.
- Dữ liệu mô phỏng khí hậu tương lai: Trích xuất từ mô hình khí hậu khu vực PRECIS (Providing Regional Climates for Impacts Studies) do Trung tâm Hadley (Vương quốc Anh) phát triển, với miền tính bao phủ kinh độ 93°E - 120°E, vĩ độ 4°N - 36°N, kích thước lưới $140 \times 160$ nút, 19 tầng thẳng đứng, sử dụng lưới xen kẽ Arakawa B và độ phân giải ngang 25 km $\times$ 25 km. Các kịch bản phát thải khí nhà kính được áp dụng gồm B1, B2, A1B và A2.
Phương pháp xử lý và phân tích số liệu:
- Phương pháp tính toán chỉ số RCCI: Tích hợp 8 biến chỉ thị gồm 4 biến thay đổi giá trị trung bình ($\Delta T$, $\Delta P%$) và 4 biến thay đổi độ biến động ($\Delta\sigma T%$, $\Delta\sigma P%$) cho 2 mùa trong năm (mùa MAMJJA gồm tháng 3 đến tháng 8 và mùa SONDJF gồm tháng 9 đến tháng 2 năm sau). Các biến được phân hạng và gán trọng số nguyên $n$ (từ 0 đến 4) hoặc chuẩn hóa qua phương pháp khoảng cách hình học gồm Khoảng cách Tiêu chuẩn Euclidean (Standard Euclidean Distance - SED) và Hệ số bình phương khoảng cách cung không đồng dạng (Squared Chord Distance - SCD).
- Phương pháp tính toán chỉ số CCI: Xác định tần suất các sự kiện mùa và năm vượt ngưỡng phân vị 95th (cực nóng, cực ẩm) và thấp hơn phân vị 5th (cực khô) của thời kỳ nền 1961-1980/1961-1990 so với giai đoạn kiểm nghiệm 1981-2000 cũng như các giai đoạn tương lai 2041-2070 và 2071-2100.
- Phương pháp hồi quy tuyến tính và thống kê đa biến: Đánh giá xu thế biến đổi dài hạn qua hệ số góc $A_1$, phân tích ma trận hệ số tương quan và tỷ lệ đóng góp của từng chỉ số thành phần vào chỉ số tổng hợp.
Lý do lựa chọn các phương pháp này xuất phát từ khả năng loại bỏ tính chủ quan của các phép so sánh đơn lẻ, đồng thời tích hợp đa chiều cả cường độ biến đổi nhiệt - ẩm và tính bất định của các hiện tượng cực đoan.
Kết quả nghiên cứu và thảo luận
Những phát hiện chính
Quá trình tính toán và phân tích chuỗi số liệu quan trắc cùng dữ liệu mô phỏng PRECIS đã mang lại 4 phát hiện quan trọng:
- Sự gia tăng nhiệt độ rõ rệt và bất đối xứng giữa các mùa: Trong giai đoạn 1981-2000 so với giai đoạn tham chiếu 1961-1980, nhiệt độ trung bình mùa lạnh (DJF) và mùa nóng (JJA) tại toàn bộ 7 vùng khí hậu đều có xu hướng tăng với mức chênh lệch trung bình từ 0,3°C đến 0,8°C. Mức tăng nhiệt mùa đông tại các vùng phía Bắc (B1, B2, B3) cao hơn mùa hè từ 0,2°C đến 0,4°C. Hệ số hồi quy tuyến tính $A_1$ của nhiệt độ tại hầu hết các trạm vùng Nam Bộ (N3) và Bắc Bộ đều mang giá trị dương rõ rệt.
- Sự phân cực lượng mưa và gia tăng biến động liên năm: Biến đổi lượng mưa thể hiện sự phân hóa không gian sâu sắc. Lượng mưa mùa khô có xu hướng giảm từ 5% đến 12% ở các vùng khí hậu phía Bắc nhưng lượng mưa mùa mưa lại tăng từ 6% đến 15% tại các vùng ven biển Trung Bộ và Nam Bộ. Hệ số biến động lượng mưa ($\Delta\sigma P%$) tăng vượt mức 15% ở nhiều trạm quan trắc, cho thấy tính bất định của chu kỳ mưa ngày càng khốc liệt.
- Nhận diện điểm nóng khí hậu qua chỉ số RCCI: Kết quả tổng hợp chỉ số RCCI lịch sử cho thấy vùng Tây Bắc (B1) và Đông Bắc (B2) là hai khu vực có giá trị chỉ số cao nhất cả nước. Tỷ lệ đóng góp của 4 chỉ số thành phần nhiệt độ vào tổng chỉ số RCCI chiếm ưu thế vượt trội với mức từ 65% đến 72%, trong khi các thành phần biến động lượng mưa đóng góp từ 28% đến 35%.
- Bùng nổ các sự kiện cực đoan trong thế kỷ 21 theo chỉ số CCI: Dự tính khí hậu giai đoạn tương lai 2071-2100 theo kịch bản phát thải cao (A2) cho thấy chỉ số CCI tăng gấp 1,5 đến 2,2 lần so với thời kỳ nền. Số năm cực nóng (vượt ngưỡng phân vị 95th) xuất hiện với tần suất áp đảo, chiếm từ 80% đến 95% số năm trong giai đoạn cuối thế kỷ 21 trên phạm vi cả 7 vùng khí hậu, đặc biệt nghiêm trọng tại Nam Trung Bộ (N1) và Tây Nguyên (N2).
Thảo luận kết quả
Nguyên nhân cốt lõi dẫn đến sự phân hóa của các chỉ số khí hậu bắt nguồn từ sự tương tác phức tạp giữa hoàn lưu gió mùa Đông Á - Nam Á và địa hình đón gió hoặc khuất gió của các dãy núi lớn (dãy Hoàng Liên Sơn, Trường Sơn). Vùng núi phía Bắc chịu tác động trực tiếp của không khí lạnh lục địa biến tính kết hợp hiệu ứng nóng lên toàn cầu làm biên độ nhiệt độ mùa đông dịch chuyển mạnh. Trong khi đó, các vùng phía Nam chịu chi phối mạnh mẽ bởi gió mùa Tây Nam và hoàn lưu đối lưu nhiệt đới, dẫn đến sự gia tăng đột biến của các đợt mưa cực đoan và kéo dài thời gian hạn hán cục bộ.
So sánh với các nghiên cứu quốc tế, kết quả của luận văn hoàn toàn tương đồng với các phát hiện của Filippo Giorgi (2006) đối với khu vực Đông Á và nghiên cứu của Michele B. Baettig (2007) trên quy mô toàn cầu. Nghiên cứu khẳng định Việt Nam là một trong những điểm nóng biến đổi khí hậu nhạy cảm nhất khu vực nhiệt đới gió mùa.
Trong thực tế phân tích, toàn bộ các kết quả này có thể được trực quan hóa tối ưu thông qua:
- Bản đồ không gian GIS: Thể hiện phân bố giá trị RCCI và CCI trên 7 vùng khí hậu dưới dạng thang màu phân cấp, giúp định vị nhanh các khu vực rủi ro cao.
- Biểu đồ phân vị hộp (Boxplot) và ma trận tương quan: Minh họa sự dịch chuyển của hàm mật độ xác suất nhiệt độ/lượng mưa từ thời kỳ nền sang thế kỷ 21, làm nổi bật tỷ lệ các năm vượt quá ngưỡng 95th và 5th.
- Bảng tổng hợp trọng số đa biến: Trình bày chi tiết đóng góp phần trăm của 8 chỉ số thành phần (ΔT, ΔP, ΔσT, ΔσP cho hai mùa) đối với từng trạm và từng tiểu vùng sinh thái.
Đề xuất và khuyến nghị
Dựa trên kết quả định lượng của bộ chỉ số RCCI và CCI, luận văn đề xuất 4 nhóm giải pháp chiến lược:
- Tích hợp chỉ số RCCI và CCI vào quy hoạch phát triển kinh tế - xã hội vùng:
- Hành động: Bộ Kế hoạch và Đầu tư chủ trì, phối hợp cùng Bộ Tài nguyên và Môi trường ban hành quy chuẩn kỹ thuật bắt buộc lồng ghép hai chỉ số biến đổi khí hậu tích hợp vào báo cáo đánh giá môi trường chiến lược (ĐMC) và quy hoạch phát triển của 7 vùng kinh tế trọng điểm.
- Target metric: 100% các dự án hạ tầng lớn và quy hoạch sử dụng đất cấp tỉnh hoàn thành đánh giá định lượng rủi ro khí hậu trước khi phê duyệt.
- Timeline: Giai đoạn 2026-2028.
- Nâng cấp độ phân giải mạng lưới quan trắc và hệ thống dự báo số trị:
- Hành động: Tổng cục Khí tượng Thủy văn đầu tư nâng cấp mô hình PRECIS kết hợp các mô hình động lực thế hệ mới (như RegCM4, WRF) đạt độ phân giải ngang dưới 10 km $\times$ 10 km, đồng thời triển khai thêm ít nhất 50 trạm đo mưa và nhiệt độ tự động tại các khu vực điểm nóng Tây Bắc và Bắc Trung Bộ.
- Target metric: Giảm sai số mô phỏng khí hậu cực đoan xuống dưới 10% và tăng độ tin cậy của kịch bản thời gian thực.
- Timeline: Giai đoạn 2026-2030.
- Tái cơ cấu mùa vụ nông nghiệp và công trình điều tiết tài nguyên nước:
- Hành động: Bộ Nông nghiệp và Phát triển Nông thôn chỉ đạo các địa phương vùng Nam Trung Bộ, Tây Nguyên và Đồng bằng sông Cửu Long chuyển đổi tối thiểu 30% diện tích lúa chịu hạn sang cây trồng cạn có giá trị kinh tế cao, đồng thời hiện đại hóa hệ thống hồ chứa phân tán.
- Target metric: Nâng tỷ lệ an toàn cấp nước nông nghiệp mùa khô đạt trên 85% và cắt giảm 20% thiệt hại do hạn mặn gây ra.
- Timeline: Hoàn thành trước năm 2030.
- Xây dựng bộ chỉ số phơi nhiễm khí hậu - xã hội cấp địa phương:
- Hành động: Các Viện nghiên cứu phối hợp với Ủy ban nhân dân các tỉnh xây dựng bản đồ nguy cơ phơi nhiễm dựa trên công thức tích hợp chỉ số CCI với mật độ dân số và chỉ số nghèo đói tại cấp huyện, xã.
- Target metric: Thiết lập hệ thống cơ sở dữ liệu số hóa giám sát tổn thương khí hậu cho 63 tỉnh, thành phố trực thuộc Trung ương.
- Timeline: Ban hành và thử nghiệm diện rộng trong giai đoạn 2026-2027.
Đối tượng nên tham khảo luận văn
Nội dung và phương pháp luận của luận văn mang lại giá trị thực tiễn cho 4 nhóm đối tượng trọng tâm:
- Cán bộ hoạch định chính sách và cơ quan quản lý nhà nước (Bộ Tài nguyên và Môi trường, Cục Biến đổi khí hậu):
- Lợi ích: Nắm bắt phương pháp luận định lượng tổng hợp để thay thế các chỉ tiêu đơn lẻ truyền thống.
- Use case: Ứng dụng công thức tính RCCI và CCI làm căn cứ khoa học phân bổ nguồn vốn đầu tư công, xây dựng Kế hoạch quốc gia thích ứng với biến đổi khí hậu (NAP) và Đóng góp do quốc gia tự quyết định (NDC).
- Chuyên gia mô hình hóa khí hậu và các nhà nghiên cứu khí tượng thủy văn:
- Lợi ích: Tiếp cận quy trình chi tiết hóa động lực bằng mô hình PRECIS (lưới Arakawa B, 19 tầng áp suất, độ phân giải 25 km) kết hợp các kỹ thuật chuẩn hóa hình học (SED, SCD).
- Use case: Làm tài liệu tham chiếu để hiệu chỉnh tham số vật lý cho các mô hình khí hậu khu vực thế hệ mới và mở rộng tính toán cho các kịch bản phát thải cập nhật.
- Kỹ sư quy hoạch thủy lợi, nông nghiệp và tài nguyên nước:
- Lợi ích: Khai thác bộ số liệu về hệ số biến động lượng mưa ($\Delta\sigma P%$) và tần suất năm khô/năm ẩm vượt ngưỡng 5th/95th.
- Use case: Tính toán dung tích thiết kế công trình hồ chứa, đê điều, cống ngăn mặn và tối ưu hóa biểu đồ điều tiết nước thích ứng các chu kỳ khô hạn kéo dài.
- Giảng viên, nghiên cứu sinh, học viên cao học và sinh viên các khối ngành Khoa học Trái đất, Môi trường:
- Lợi ích: Tài liệu học thuật chuẩn mực về phương pháp xử lý chuỗi số liệu khí hậu dài hạn, thống kê phân vị và phân tích độ bất định trong dự tính khí hậu tương lai.
- Use case: Sử dụng làm giáo trình tham khảo, khung phương pháp luận cho các đề tài khóa luận, luận văn thạc sĩ và bài báo khoa học quốc tế.
Câu hỏi thường gặp
Chỉ số RCCI khác biệt như thế nào so với các phương pháp đánh giá biến đổi khí hậu truyền thống? Các phương pháp truyền thống thường chỉ đánh giá mức chênh lệch trung bình của nhiệt độ hoặc lượng mưa giữa hai thời kỳ. Chỉ số RCCI của Filippo Giorgi tích hợp đồng thời 8 biến chỉ thị bao gồm cả giá trị trung bình lẫn độ biến động liên năm ($\Delta\sigma T$, $\Delta\sigma P$) cho 2 mùa trong năm, giúp nhận diện chính xác các điểm nóng khí hậu nhạy cảm nhất.
Tại sao luận văn lại lựa chọn mô hình PRECIS để mô phỏng khí hậu tương lai tại Việt Nam? PRECIS là mô hình khí hậu khu vực độ phân giải cao (25 km $\times$ 25 km) do Trung tâm Hadley phát triển, sử dụng 19 tầng khí quyển và lưới Arakawa B. Mô hình có khả năng mô phỏng vượt trội ảnh hưởng của địa hình phức tạp và hoàn lưu gió mùa tại Việt Nam, khắc phục hoàn toàn nhược điểm độ phân giải thô của các mô hình hoàn lưu chung toàn cầu (GCM).
Vùng khí hậu nào ở Việt Nam chịu mức độ biến đổi nhiệt độ và lượng mưa khốc liệt nhất? Theo kết quả tính toán chỉ số RCCI lịch sử và mô phỏng tương lai, khu vực Tây Bắc (B1) và Đông Bắc (B2) ghi nhận mức biến đổi nhiệt độ mạnh nhất với tỷ lệ đóng góp nhiệt vượt 70%. Trong khi đó, khu vực Nam Trung Bộ (N1) và Tây Nguyên (N2) chịu ảnh hưởng nghiêm trọng nhất về tần suất xuất hiện các năm cực nóng và cực khô theo chỉ số CCI.
Các ngưỡng phân vị 95th và 5th trong chỉ số CCI có ý nghĩa thực tiễn gì? Phân vị 95th đại diện cho các giá trị cực đại (chỉ xuất hiện với xác suất 5% trong lịch sử, tương ứng các năm cực nóng hoặc mưa đặc biệt lớn), trong khi phân vị 5th đại diện cho ngưỡng cực tiểu (năm hạn hán gay gắt nhất). Khi tần suất vượt ngưỡng 95th trong tương lai đạt trên 80%, điều đó báo hiệu các hiện tượng thời tiết cực đoan đã trở thành trạng thái phổ biến mới.
Làm thế nào để kết hợp chỉ số CCI vào việc đánh giá tính dễ bị tổn thương kinh tế - xã hội? Dựa trên mô hình của Noah S. Diffenbaugh, chỉ số CCI đóng vai trò là biến số phơi nhiễm vật lý. Khi nhân chỉ số CCI với tổ hợp các chỉ số trọng số về dân số, tỷ lệ nghèo đói và chỉ số phát triển con người (HDI), các nhà quản lý sẽ thu được bản đồ định lượng chính xác các khu vực dân cư dễ bị tổn thương nhất để ưu tiên cứu trợ và đầu tư công.
Kết luận
- Luận văn đã thiết lập thành công quy trình đánh giá định lượng biến đổi khí hậu tại Việt Nam thông qua việc ứng dụng hai chỉ số tích hợp hiện đại RCCI và CCI trên chuỗi số liệu quan trắc lịch sử 1961-2000 và mô phỏng tương lai đến năm 2100.
- Kết quả nghiên cứu khẳng định sự gia tăng nhiệt độ trung bình từ 0,3°C đến 0,8°C trên cả nước trong giai đoạn 1981-2000 so với nền 1961-1980, với mùa đông tại miền Bắc có mức tăng nhiệt vượt trội so với mùa hè.
- Mô phỏng tương lai bằng mô hình PRECIS độ phân giải 25 km theo các kịch bản phát thải chỉ ra rằng vào cuối thế kỷ 21, tần suất các năm cực nóng vượt ngưỡng phân vị 95th sẽ chiếm từ 80% đến 95% tổng số năm mô phỏng.
- Tây Bắc, Đông Bắc và Nam Trung Bộ được định lượng là các điểm nóng khí hậu trọng điểm, đòi hỏi các chính sách ưu tiên thích ứng và phân bổ ngân sách giảm nhẹ thiên tai ngay trong giai đoạn 2026-2030.
- Về lộ trình tiếp theo, các cơ quan quản lý cần sớm chuẩn hóa việc đưa bộ chỉ số biến đổi khí hậu tích hợp vào khung pháp lý quy hoạch không gian và phát triển kinh tế vùng nhằm bảo đảm phát triển bền vững quốc gia.
Quý độc giả, các nhà nghiên cứu và cơ quan quản lý quan tâm đến phương pháp luận định lượng khí hậu chi tiết có thể tra cứu toàn văn luận văn thạc sĩ của tác giả Trịnh Hoàng Dương tại Trung tâm Thông tin - Thư viện Đại học Quốc gia Hà Nội để khai thác trọn vẹn hệ thống dữ liệu và mô hình tính toán.