Tổng quan nghiên cứu

Biến đổi khí hậu toàn cầu đang diễn ra với tốc độ ngày càng phức tạp và khó lường. Theo báo cáo đánh giá lần thứ năm của Ban Liên chính phủ về Biến đổi Khí hậu, nhiệt độ bề mặt trái đất trong giai đoạn 1901 - 2012 đã tăng 0,89°C, trong đó giai đoạn 1951 - 2012 tăng tới 0,72°C. Tại Việt Nam, số liệu quan trắc trong 50 năm qua cho thấy nhiệt độ trung bình năm đã tăng từ 0,5°C đến 0,7°C, mực nước biển dâng khoảng 20 cm, kéo theo sự gia tăng rõ rệt về tần suất và cường độ của các thiên tai nguy hiểm. Tuy nhiên, các nghiên cứu trước đây tại Việt Nam thường chỉ tiếp cận các yếu tố khí hậu trung bình hoặc đánh giá riêng lẻ từng biến khí tượng mà chưa xây dựng được một hệ thống chỉ số cực đoan đồng bộ, có tính tương thích quốc tế.

Luận văn thạc sĩ khoa học chuyên ngành Khí tượng và khí hậu học với đề tài "Đánh giá biến đổi khí hậu ở Việt Nam dựa trên chỉ số khí hậu cực đoan của IPCC" của tác giả Nguyễn Đình Dũng, dưới sự hướng dẫn khoa học của Giáo sư Phan Văn Tân tại Trường Đại học Khoa học Tự nhiên - Đại học Quốc gia Hà Nội, đã giải quyết căn bản khoảng trống học thuật này. Mục tiêu nghiên cứu cụ thể của luận văn bao gồm: thứ nhất, đánh giá quy luật và xu thế biến đổi của các cực đoan khí hậu trên quy mô toàn quốc trong giai đoạn lịch sử 1961 - 2010 thông qua hệ thống chỉ số chuẩn hóa về nhiệt độ và lượng mưa; thứ hai, xác định mức độ biến đổi của các hiện tượng cực đoan khí hậu trong tương lai ở hai thời kỳ giữa thế kỷ 21 (2040 - 2059) và cuối thế kỷ 21 (2080 - 2099) theo kịch bản phát thải trung bình A1B bằng mô hình khí hậu khu vực PRECIS.

Phạm vi nghiên cứu bao quát mạng lưới 92 trạm khí tượng bề mặt và 96 trạm đo mưa trên khắp 7 vùng khí hậu của Việt Nam. Ý nghĩa thực tiễn của công trình thể hiện ở việc lượng hóa chính xác các metric rủi ro khí hậu như tốc độ biến đổi nhiệt độ cực trị (°C/10 năm), tần suất ngày nắng nóng (ngày/10 năm) và lượng mưa ngày cực đại (mm/50 năm). Kết quả này là cơ sở khoa học trực tiếp phục vụ công tác quy hoạch phát triển kinh tế - xã hội, quản lý rủi ro thiên tai và xây dựng chiến lược thích ứng với biến đổi khí hậu 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 vận dụng khung lý thuyết chuẩn mực của Nhóm Chuyên gia về Phát hiện, Theo dõi và Chỉ số hóa Biến đổi Khí hậu (ETCCDI) trực thuộc Tổ chức Khí tượng Thế giới (WMO) và Ủy ban Khí hậu Toàn cầu (CLIVAR/JCOMM), tài liệu kỹ thuật WCDMP-No.72. Bộ khung này định nghĩa 27 chỉ số cực đoan khí hậu tiêu chuẩn toàn cầu, dựa trên nguyên lý phân bố xác suất cực trị và các ngưỡng phân vị phi tham số. Bên cạnh đó, luận văn áp dụng lý thuyết mô hình hóa khí hậu khu vực động lực học thông qua hệ thống mô hình PRECIS (Providing REgional Climates for Impacts Studies) do Trung tâm Hadley (Vương quốc Anh) phát triển, giúp chi tiết hóa các kịch bản toàn cầu xuống độ phân giải ngang 25 km x 25 km.

Các khái niệm then chốt được chuẩn hóa trong nghiên cứu gồm:

  • Chỉ số cực đoan khí hậu (Climate Extreme Indices - CEI): Các đại lượng thống kê dẫn xuất từ chuỗi quan trắc ngày nhằm định lượng hóa tần suất, cường độ và thời gian kéo dài của các trạng thái thời tiết khắc nghiệt.
  • Hiện tượng khí hậu cực trị (Extreme Climate Events - ECE): Trạng thái thời tiết có tần suất xuất hiện thống kê nhỏ hơn 10% trong phân bố xác suất nhiều năm tại một khu vực cụ thể.
  • Ngưỡng phân vị (Percentiles): Các mốc phân vị 10%, 90%, 95%, 99% được xác định từ hàm tích lũy xác suất trong thời kỳ chuẩn 1971 - 2000 nhằm phân loại ngày lạnh, đêm lạnh, ngày ấm, đêm nóng và lượng mưa rất lớn.
  • Kịch bản phát thải SRES A1B: Kịch bản phát thải khí nhà kính trung bình, mô tả một thế giới tương lai tăng trưởng kinh tế nhanh với sự cân bằng giữa các nguồn năng lượng hóa thạch và phi hóa thạch.

Bộ chỉ số được tác giả tinh chỉnh để phù hợp với điều kiện nhiệt đới gió mùa tại Việt Nam, bao gồm 16 chỉ số nhiệt độ (bổ sung ngưỡng nắng nóng SU35, SU37, đêm nhiệt đới TR25, TR30) và 11 chỉ số lượng mưa (bổ sung ngưỡng mưa lớn R50, R100).

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

Nguồn dữ liệu quan trắc lịch sử được thu thập từ Trung tâm Tư liệu Khí tượng Thủy văn Quốc gia với chuỗi số liệu nhiệt độ tối cao ngày (Tx), nhiệt độ tối thấp ngày (Tn) tại 92 trạm và lượng mưa ngày (R) tại 96 trạm (bổ sung 4 trạm Mộc Châu, Bắc Giang, Tĩnh Gia, Phan Thiết) trong 50 năm liên tục từ năm 1961 đến năm 2010. Cỡ mẫu này đại diện toàn diện cho 7 vùng khí hậu đặc trưng từ Tây Bắc Bộ, Đông Bắc Bộ, Đồng bằng Bắc Bộ, Bắc Trung Bộ, Nam Trung Bộ, Tây Nguyên đến Nam Bộ. Phương pháp chọn mẫu tập trung vào các trạm có độ dài chuỗi quan trắc trên 40 năm, tỷ lệ số liệu khuyết thiếu dưới 5%, đảm bảo tính đại diện địa lý và tính đồng nhất khí hậu.

Quy trình xử lý số liệu được thực hiện nghiêm ngặt qua phần mềm RHtestsV3 nhằm kiểm soát chất lượng, phát hiện và loại bỏ sai số thô, xử lý mã khuyết (-99.9) và kiểm định tính đồng nhất của chuỗi thời gian. Phương pháp phân tích xu thế biến đổi sử dụng mô hình hồi quy tuyến tính bằng phương pháp bình phương tối thiểu (Ordinary Least Squares - OLS) với phương trình dạng xt = b0 + b1*t, trong đó b1 là tốc độ biến đổi và b0 là gốc xu thế. Việc tính toán các chỉ số phân vị được tự động hóa qua phần mềm chuyên dụng RClimDex dựa trên nền tảng ngôn ngữ R với thời kỳ cơ sở chuẩn 1971 - 2000. Lý do lựa chọn phương pháp này là đảm bảo độ tin cậy thống kê cao, tuân thủ đúng chuẩn mực của WMO, giúp kết quả nghiên cứu của Việt Nam có thể so sánh trực tiếp với các báo cáo toàn cầu của IPCC.

Đối với dữ liệu tương lai, nghiên cứu khai thác sản phẩm dự tính từ mô hình PRECIS do Viện Khoa học Khí tượng Thủy văn và Biến đổi khí hậu cung cấp, sử dụng tổ hợp đa thành phần gồm 5 phương án (Q0, Q3, Q10, Q11, Q13) từ mô hình toàn cầu HadCM3 nhằm giảm thiểu độ bất định của kịch bản phát thải A1B trong hai giai đoạn 2040 - 2059 và 2080 - 2099.

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

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

Kết quả tính toán chuỗi số liệu quan trắc 50 năm (1961 - 2010) tại mạng lưới trạm khí tượng đã làm rõ những biến đổi sâu sắc của khí hậu cực đoan tại Việt Nam:

  • Xu thế gia tăng vượt trội của các hiện tượng nhiệt độ cực cao: Chỉ số số ngày nắng nóng (SU35 với Tx > 35°C) tăng tại 71 trên tổng số 92 trạm quan trắc (chiếm tỷ lệ 77,17%), trong đó có 34 trạm ghi nhận tốc độ tăng trên 1,0 ngày/10 năm. Mức tăng mạnh nhất xấp xỉ 8,0 ngày/10 năm tại Châu Đốc (An Giang) và Tuyên Hóa (Quảng Bình). Chỉ số nhiệt độ tối cao tuyệt đối (TXx) tăng tại 68/92 trạm (chiếm 73,91%), mức tăng kỷ lục đạt 0,74°C/10 năm tại Càng Long (Trà Vinh) và 0,70°C/10 năm tại đảo Trường Sa.

  • Hiện tượng ấm lên rõ rệt của nhiệt độ cực tiểu: Nhiệt độ tối thấp tuyệt đối (TNn) có xu thế tăng tại 87/92 trạm (chiếm 94,57%), với 59 trạm tăng trên 0,1°C/10 năm. Nhiều trạm ghi nhận mức tăng vượt 1,0°C/10 năm như M'Drắk (1,56°C/10 năm), Kon Tum (1,43°C/10 năm), Châu Đốc (1,36°C/10 năm), Phan Rang (1,15°C/10 năm), Điện Biên và Đắk Nông (1,04°C/10 năm). Số ngày ấm (TX90p) tăng tại 53/92 trạm (57,61%), tăng cao nhất tại các trạm hải đảo như Trường Sa (+15,14 ngày/10 năm) và Phú Quốc (+8,66 ngày/10 năm). Trong khi đó, số đêm mát (Tn90p) và đêm lạnh (TN10p) giảm mạnh, điển hình tại Vũng Tàu giảm 7,72 ngày/10 năm và Cần Thơ giảm 6,66 ngày/10 năm.

  • Phân hóa không gian sâu sắc của lượng mưa cực đoan: Lượng mưa ngày lớn nhất tháng (Rx1day) và lượng mưa 5 ngày lớn nhất tháng (Rx5day) có xu hướng tăng mạnh từ 5% đến 20% tại các vùng Nam Trung Bộ, Tây Nguyên và Nam Bộ. Ngược lại, lượng mưa cực đoan có xu hướng giảm nhẹ từ 5% đến 13% tại một số khu vực Đồng bằng và Đông Bắc Bộ trong mùa mưa. Số ngày mưa lớn (R50 với lượng mưa từ 50 mm trở lên) và mưa rất lớn (R100 từ 100 mm trở lên) tăng đáng kể ở dải ven biển Trung Bộ.

  • Dự tính tương lai theo mô hình PRECIS (2040 - 2059 và 2080 - 2099): Nhiệt độ tối cao và tối thấp tuyệt đối tiếp tục tăng từ 1,5°C đến trên 3,5°C trên phạm vi cả nước vào cuối thế kỷ 21. Số ngày nắng nóng gay gắt (SU37) dự kiến gia tăng từ 20 đến 45 ngày mỗi năm. Các chỉ số mưa cực đoan (Rx1day, Rx5day, R95p) tăng từ 10% đến hơn 30% ở hầu hết các vùng, trong khi số ngày khô hạn liên tục (CDD) gia tăng cục bộ tại một số khu vực Nam Trung Bộ.

Thảo luận kết quả

Sự gia tăng mạnh mẽ của các chỉ số nhiệt độ cực đoan phản ánh xu thế nóng lên toàn cầu dưới tác động của gia tăng nồng độ khí nhà kính. Đáng chú ý, tốc độ tăng của nhiệt độ tối thấp tuyệt đối (TNn) cao hơn rõ rệt so với nhiệt độ tối cao tuyệt đối (TXx), dẫn tới hiện tượng thu hẹp biên độ nhiệt độ ngày đêm (DTR) trên diện rộng. Cơ chế vật lý đằng sau hiện tượng này là sự gia tăng hơi nước và mây trong khí quyển nhiệt đới, giữ lại bức xạ sóng dài ban đêm hiệu quả hơn.

Sự phân hóa của các chỉ số lượng mưa cực đoan giữa miền Bắc và miền Nam hoàn toàn tương thích với các báo cáo AR4, AR5 của IPCC và các nghiên cứu của Bộ Tài nguyên và Môi trường công bố năm 2012. Hoạt động bất thường của hoàn lưu gió mùa mùa hè kết hợp với các đợt ENSO (El Niño và La Niña) đã làm gia tăng các trận mưa cường suất lớn tại Trung Bộ và Nam Bộ, trong khi kéo dài chuỗi ngày khô hạn tại một số thời điểm trong năm.

Để trực quan hóa các phát hiện này, dữ liệu nghiên cứu được trình bày tối ưu qua hệ thống bản đồ phân bố dạng điểm (point distribution maps) thể hiện độ dốc xu thế b1 trên nền địa hình không nội suy, kết hợp với các biểu đồ chuỗi thời gian (time-series plots) đường hồi quy và biểu đồ phân bố tần suất tích lũy. Cách biểu diễn này giúp các nhà nghiên cứu và quản lý dễ dàng so sánh mức độ biến đổi cực đoan giữa các trạm ven biển, hải đảo và lục địa.

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

Dựa trên các kết luận định lượng về xu thế biến đổi khí hậu cực đoan, luận văn cung cấp cơ sở để đưa ra các đề xuất và khuyến nghị hành động thiết thực sau:

  • Nâng cấp hệ thống quan trắc và cảnh báo sớm thiên tai khí tượng thủy văn: Hiện đại hóa 100% mạng lưới trạm khí tượng tự động, tích hợp radar thời tiết độ phân giải cao nhằm nâng cao độ chính xác dự báo mưa cực đoan định lượng thêm 25 - 30%, cảnh báo sớm các đợt nắng nóng gay gắt trước 3 đến 5 ngày. Đơn vị thực hiện: Cục Khí tượng Thủy văn phối hợp cùng Viện Khoa học Khí tượng Thủy văn và Biến đổi khí hậu, hoàn thành trong giai đoạn 2025 - 2030.

  • Rà soát và điều chỉnh quy chuẩn kỹ thuật công trình hạ tầng và thoát nước đô thị: Cập nhật chỉ số lượng mưa ngày lớn nhất (Rx1day) và lượng mưa 5 ngày lớn nhất (Rx5day) vào tiêu chuẩn thiết kế hệ thống tiêu thoát nước, nâng hệ số an toàn lưu lượng xả lũ của đê điều và hồ chứa thêm 15 - 20% tại các đô thị trọng điểm như Hà Nội, Thành phố Hồ Chí Minh, Đà Nẵng và Cần Thơ. Đơn vị chủ trì: Bộ Xây dựng phối hợp cùng Bộ Nông nghiệp và Phát triển Nông thôn, triển khai trong giai đoạn 2025 - 2028.

  • Chuyển đổi cơ cấu cây trồng và mùa vụ nông nghiệp thích ứng với nhiệt độ cao: Chuyển đổi tối thiểu 20% diện tích đất nông nghiệp thường xuyên bị hạn hán tại Nam Trung Bộ và Tây Nguyên sang các giống cây trồng chịu nhiệt, tiết kiệm nước; áp dụng công nghệ tưới nhỏ giọt nhằm giảm 30% lượng nước tiêu thụ trong các giai đoạn số ngày khô liên tục (CDD) kéo dài. Đơn vị thực hiện: Bộ Nông nghiệp và Phát triển Nông thôn cùng Sở Nông nghiệp các địa phương, thời gian từ 2025 đến 2027.

  • Xây dựng kế hoạch hành động ứng phó biến đổi khí hậu dựa trên bộ chỉ số cực đoan cấp tỉnh: Ban hành tài liệu hướng dẫn kỹ thuật áp dụng bộ 27 chỉ số ETCCDI vào đánh giá rủi ro khí hậu định kỳ 5 năm cho 63 tỉnh, thành phố trực thuộc trung ương, ưu tiên xây dựng bản đồ rủi ro ngập lụt và sóng nhiệt cho 28 tỉnh ven biển. Đơn vị thực hiện: Ủy ban nhân dân các tỉnh, thành phố phối hợp với Bộ Tài nguyên và Môi trường, hoàn thành trước năm 2027.

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

Công trình nghiên cứu mang lại giá trị thực tiễn và học thuật to lớn cho nhiều nhóm độc giả chuyên biệt:

  • Cơ quan quản lý nhà nước và nhà hoạch định chính sách: Bộ Tài nguyên và Môi trường, Bộ Kế hoạch và Đầu tư, Cục Biến đổi khí hậu có thể trực tiếp sử dụng số liệu định lượng 50 năm để cập nhật Kịch bản Biến đổi Khí hậu Quốc gia, xây dựng Chiến lược quốc gia về Biến đổi khí hậu và Kế hoạch thích ứng quốc gia (NAP).

  • Giảng viên, nghiên cứu sinh và sinh viên ngành Khí tượng, Thủy văn, Biến đổi khí hậu: Kế thừa khung phương pháp luận chuẩn mực của ETCCDI, quy trình làm sạch dữ liệu bằng RHtestsV3, phương pháp phân tích chỉ số qua RClimDex và quy trình chạy mô hình động lực học PRECIS để phát triển các đề tài nghiên cứu chuyên sâu.

  • Kỹ sư quy hoạch đô thị, công trình thủy lợi và giao thông vận tải: Ứng dụng các thông số tính toán cực trị mưa (Rx1day, Rx5day, R50, R100) và nhiệt độ cực trị (TXx, TNn) để tính toán tải trọng nhiệt của vật liệu, thiết kế khẩu độ cầu cống, nâng cấp đê biển và tính toán dung tích phòng lũ cho các hồ chứa thủy điện, thủy lợi.

  • Các tổ chức phi chính phủ, doanh nghiệp bảo hiểm và phát triển nông nghiệp: Sử dụng ma trận rủi ro khí hậu cực đoan và số liệu phân vị 90%, 95% để thiết kế các gói bảo hiểm chỉ số thời tiết nông nghiệp, xây dựng các dự án tài trợ thích ứng biến đổi khí hậu dựa vào cộng đồng tại các điểm nóng khí hậu như Tây Nguyên và Đồng bằng sông Cửu Long.

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

  • Bộ chỉ số cực đoan khí hậu của IPCC có ưu điểm gì vượt trội so với cách tiếp cận truyền thống? Bộ chỉ số của IPCC do nhóm ETCCDI thiết lập cung cấp hệ thống 27 chỉ số chuẩn hóa từ số liệu nhiệt độ và lượng mưa ngày. Phương pháp này lượng hóa chính xác các sự kiện hiếm thông qua các ngưỡng phân vị thống kê phi tham số (10%, 90%, 95%), khắc phục hoàn toàn nhược điểm của việc chỉ đánh giá giá trị trung bình tháng hoặc năm vốn thường làm mờ đi các cực trị nguy hiểm.

  • Vì sao tác giả cần hiệu chỉnh các chỉ số chuẩn như SU25, R10 khi áp dụng tại Việt Nam? Chỉ số quốc tế SU25 (nhiệt độ tối cao ngày trên 25°C) hay R10 (lượng mưa từ 10 mm/ngày) được xây dựng cho vùng ôn đới, nơi nhiệt độ và lượng mưa thấp. Tại vùng nhiệt đới gió mùa như Việt Nam, các ngưỡng này diễn ra gần như quanh năm. Việc điều chỉnh thành SU35, SU37 (nắng nóng và nắng nóng gay gắt) cùng R50, R100 (mưa lớn và mưa rất lớn) phản ánh chính xác ngưỡng thiên tai thực tế.

  • Nguyên nhân nào khiến nhiệt độ tối thấp tuyệt đối (TNn) tăng mạnh hơn nhiệt độ tối cao tuyệt đối (TXx)? Nhiệt độ tối thấp tuyệt đối tăng nhanh hơn là do sự gia tăng nồng độ khí nhà kính và độ ẩm trong khí quyển, làm tăng hiệu ứng giữ nhiệt bức xạ sóng dài vào ban đêm. Kết quả là ban đêm ấm lên nhanh hơn ban ngày, làm giảm biên độ dao động nhiệt ngày (DTR) trên 94,57% số trạm quan trắc toàn quốc.

  • Mô hình PRECIS đóng góp vai trò gì trong việc dự tính các chỉ số khí hậu tương lai? Mô hình khí hậu khu vực PRECIS giải quyết bài toán chi tiết hóa động lực học từ mô hình hoàn lưu chung toàn cầu HadCM3 xuống độ phân giải cao 25 km x 25 km. Phương pháp tổ hợp 5 kịch bản thành phần (Q0, Q3, Q10, Q11, Q13) giúp kiểm soát và đánh giá độ bất định của dự tính khí hậu cực đoan trong thế kỷ 21.

  • Dữ liệu lượng mưa cực đoan trong luận văn hỗ trợ như thế nào cho công tác quy hoạch đô thị? Dữ liệu về các chỉ số Rx1day, Rx5day và lượng mưa ngày rất ẩm R95p cung cấp cơ sở để kỹ sư xác định chính xác cường độ mưa tính toán theo chu kỳ lặp lại 10 năm, 50 năm hoặc 100 năm. Nhờ đó, kích thước cống ngầm và hồ điều hòa đô thị được thiết kế tối ưu, giảm thiểu nguy cơ ngập lụt lịch sử.

Kết luận

  • Luận văn đã chuẩn hóa thành công chuỗi số liệu quan trắc 50 năm (1961 - 2010) tại 92 trạm nhiệt độ và 96 trạm đo mưa trên toàn lãnh thổ Việt Nam bằng phần mềm kiểm soát chất lượng RHtestsV3.
  • Thiết lập bức tranh toàn diện về xu thế ấm lên rõ rệt với 94,57% số trạm tăng nhiệt độ tối thấp tuyệt đối và 77,17% số trạm tăng số ngày nắng nóng SU35, trong đó mức tăng kỷ lục đạt xấp xỉ 8 ngày/10 năm.
  • Xác định sự phân hóa sâu sắc của lượng mưa cực đoan với xu thế tăng mạnh từ 5% đến 20% tại Nam Trung Bộ và Nam Bộ, đồng thời gia tăng tần suất các trận mưa lớn vượt ngưỡng R50, R100.
  • Hoàn thành dự tính khí hậu cực đoan cho hai thời kỳ 2040 - 2059 và 2080 - 2099 bằng tổ hợp 5 phương án mô hình PRECIS, chỉ ra mức tăng nhiệt độ cực trị từ 1,5°C đến 3,5°C và lượng mưa cực đoan tăng từ 10% đến 30%.
  • Đóng góp bộ công cụ và hệ thống chỉ số cực đoan IPCC được hiệu chuẩn hoàn hảo cho điều kiện nhiệt đới gió mùa Việt Nam, phục vụ công tác giám sát khí hậu theo chuẩn quốc tế.

Trong giai đoạn tiếp theo từ năm 2025 đến năm 2030, các cơ quan nghiên cứu cần tiếp tục cập nhật chuỗi số liệu quan trắc đến hiện tại, mở rộng đánh giá theo các kịch bản nồng độ khí nhà kính mới (RCP và SSP) với độ phân giải cao hơn. Hãy khai thác ngay kết quả nghiên cứu này để tích hợp vào các quy hoạch phát triển bền vững và chiến lược phòng chống thiên tai tại địa phương của bạn.