Tổng quan nghiên cứu

Biến đổi khí hậu toàn cầu đang diễn biến phức tạp và tác động sâu sắc đến mọi mặt đời sống kinh tế - xã hội. Theo báo cáo của Ban Liên chính phủ về Biến đổi Khí hậu, nhiệt độ trung bình bề mặt Trái Đất đã gia tăng khoảng 0,74°C trong giai đoạn 1906–2005, trong đó tốc độ ấm lên ở 50 năm gần đây cao gần gấp đôi giai đoạn trước. Tại Việt Nam, chuỗi số liệu quan trắc 50 năm từ 1958 đến 2007 cho thấy nền nhiệt độ trung bình năm đã tăng từ 0,5°C đến 0,7°C, đi kèm sự gia tăng về tần suất và cường độ của các hiện tượng thời tiết nguy hiểm. Vấn đề nghiên cứu trọng tâm là đánh giá một cách định lượng tác động của biến đổi khí hậu đến các cực trị nhiệt độ và lượng mưa – những yếu tố trực tiếp chi phối rủi ro thiên tai như nắng nóng gay gắt, hạn hán, lũ lụt và rét đậm, rét hại.

Mục tiêu cụ thể của công trình là kiểm chứng năng lực mô phỏng khí hậu của mô hình động lực khu vực CCAM cho thời kỳ chuẩn 1980–1999, đồng thời xây dựng bức tranh dự tính chi tiết về sự biến đổi của các yếu tố cực trị nhiệt độ và lượng mưa trong thế kỷ XXI. Phạm vi không gian nghiên cứu bao quát toàn bộ lãnh thổ Việt Nam với 7 vùng khí hậu đặc trưng, trải dài qua 3 giai đoạn tương lai gồm 2011–2040, 2041–2070 và 2071–2100 theo hai kịch bản phát thải khí nhà kính trung bình A1B và phát thải cao A2. Nghiên cứu mang ý nghĩa khoa học và thực tiễn to lớn khi cung cấp bộ chỉ số cực trị ở độ phân giải cao 25 km, tạo tiền đề định lượng vững chắc giúp các nhà quản lý xây dựng kế hoạch hành động ứng phó thiên tai, giảm thiểu thiệt hại kinh tế ước tính hàng nghìn tỷ đồng mỗi năm do khí hậu cực đoan gây ra.

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 nền tảng lý thuyết biến đổi khí hậu theo Công ước khung của Liên Hợp Quốc về Biến đổi Khí hậu (UNFCCC 1992) 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 2001, 2007). Khung lý thuyết phân định rõ giữa dao động khí hậu tự nhiên và biến đổi khí hậu nhân vi do nồng độ khí nhà kính gia tăng. Về mặt cực trị, nghiên cứu tiếp cận định nghĩa hiện tượng khí hậu cực đoan theo IPCC và các công trình khí tượng học hiện đại: hiện tượng cực đoan là trạng thái thời tiết hiếm gặp với xác suất xuất hiện dưới 10%, có cường độ lớn và tính chất khắc nghiệt đe dọa sinh thái và nhân sinh.

Hệ thống chỉ số nghiên cứu kế thừa bộ 27 chỉ số khí hậu tiêu chuẩn do nhóm chuyên gia ETCCDI thuộc Trung tâm Hadley và IPCC thiết lập, kết hợp cùng các quy chuẩn khí tượng tác nghiệp tại Việt Nam. Bốn nhóm khái niệm cốt lõi được định lượng bao gồm:

  1. Nhiệt độ cực đại ngày (Tx) và số ngày nắng nóng trong tháng có nhiệt độ cực đại vượt ngưỡng 35°C (Tx35).
  2. Nhiệt độ cực tiểu ngày (Tm) và số ngày trong tháng có nhiệt độ cực tiểu dưới ngưỡng 15°C (Tm15) phản ánh điều kiện rét đậm, rét hại.
  3. Lượng mưa ngày cực đại (Rx) và số ngày mưa lớn trong tháng đạt lượng mưa từ 50 mm trở lên (R50).
  4. Lượng mưa tích lũy 5 ngày liên tiếp lớn nhất trong tháng (Rx5day) phản ánh nguy cơ ngập lụt diện rộng.

Về kịch bản phát thải, luận văn áp dụng hai họ kịch bản SRES: kịch bản A1B thể hiện 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 với tốc độ tăng trưởng kinh tế nhanh; kịch bản A2 mô tả viễn cảnh thế giới phân mảnh với mức phát thải khí nhà kính ở ngưỡng rất cao.

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

Nghiên cứu tích hợp nguồn dữ liệu đa tầng với độ tin cậy cao. Cỡ mẫu thực nghiệm bao gồm chuỗi số liệu quan trắc hàng ngày liên tục 20 năm trong thời kỳ chuẩn 1980–1999 tại 73 trạm khí tượng bề mặt tiêu chuẩn. Phương pháp chọn mẫu là chọn mẫu có chủ đích đại diện, phân bố đồng đều trên cả 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 và Đồng bằng Nam Bộ, đảm bảo loại bỏ hoàn toàn các sai số thô trong quá trình đo đạc. Nguồn dữ liệu kiểm chứng độc lập bao gồm bộ số liệu tái phân tích toàn cầu ERA-40 của Trung tâm Dự báo Thời tiết Hạn vừa Châu Âu tại mực áp suất 850 hPa và bộ dữ liệu phân tích mưa, nhiệt dạng lưới APHRODITE độ phân giải 0,5° x 0,5° của Nhật Bản.

Công cụ mô phỏng chính là mô hình hoàn lưu khí quyển bảo giác lập phương CCAM (Cubic Conformal Atmospheric Model) do tổ chức CSIRO (Úc) phát triển. Lý do lựa chọn CCAM xuất phát từ ưu thế vượt trội của lưới bảo giác lập phương kết hợp phương pháp tích phân bán Lagrange (semi-Lagrangian), giúp giải quyết triệt để vấn đề hội tụ kinh tuyến và biến dạng lưới tại các vùng cực mà mô hình kinh vĩ độ truyền thống gặp phải. Cấu hình miền tính bao gồm 121x121 điểm lưới, tâm đặt tại tọa độ 11,5°B - 107,5°Đ, bao phủ từ 0° đến 30°B và 90° đến 120°Đ với 18 mực sigma thẳng đứng và độ phân giải ngang đạt 25 km. Timeline tích phân mô phỏng được thực hiện liên tục từ ngày 01/01/2000 đến 31/12/2100 cho cả hai kịch bản A1B và A2, sau khi đã hoàn thành bước mô phỏng kiểm định cho giai đoạn chuẩn 1980–1999.

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à phân tích mô phỏng đem lại bốn phát hiện khoa học quan trọng:

Thứ nhất, trong thời kỳ chuẩn 1980–1999, mô hình CCAM tái tạo rất tốt hoàn lưu gió và trường nhiệt độ tại mực 850 hPa so với dữ liệu tái phân tích ERA-40. Trường nhiệt độ mô phỏng phản ánh đúng sự phân hóa theo vĩ độ và địa hình, dù tồn tại độ lệch thiên dương từ 1°C đến 3°C trên khu vực đất liền. Tốc độ gió được mô phỏng sát thực tế với sai số nhỏ chỉ khoảng 1 đến 2 m/s, trong đó gió mùa mùa hè đạt mức độ tương đồng cao nhất giữa mô hình và thực tế.

Thứ hai, khi đối soát với 73 trạm quan trắc mặt đất, CCAM tái hiện chính xác biến trình năm của nhiệt độ cực tiểu trung bình tháng tại 6 trên 7 vùng khí hậu với sai số trung bình chỉ từ 1°C đến 2°C. Nhiệt độ cực đại cũng đạt độ chuẩn xác rất cao tại Nam Trung Bộ, Tây Nguyên và Nam Bộ với độ lệch dưới 2°C. Riêng khu vực Tây Bắc ghi nhận độ lệch nhiệt độ cực tiểu lớn hơn vào tháng 1 (lên đến khoảng 9°C) do ảnh hưởng của địa hình núi cao hiểm trở làm sai lệch phép nội suy lưới.

Thứ ba, về trường mưa, mô hình CCAM nắm bắt rất tốt quy luật phân bố mùa và thời điểm cực đại lượng mưa trên toàn lãnh thổ. Tuy nhiên, mô hình có xu thế mô phỏng thấp hơn lượng mưa thực tế tại các tháng cao điểm mưa bão ở miền Trung. Điển hình tại Bắc Trung Bộ trong tháng 10, lượng mưa mô phỏng thấp hơn giá trị quan trắc khoảng 380 mm, cho thấy hiệu ứng địa hình chắn gió chưa được mô hình hóa triệt để ở pha đối lưu mạnh.

Thứ tư, kết quả dự tính khí hậu thế kỷ XXI chỉ ra xu thế ấm lên áp đảo. Số ngày nắng nóng cực đoan (Tx35) gia tăng mạnh mẽ qua 3 giai đoạn 2011–2040, 2041–2070 và 2071–2100 dưới cả hai kịch bản phát thải A1B và A2. Ngược lại, số ngày có nhiệt độ cực tiểu thấp (Tm15) suy giảm rõ rệt trên diện rộng, đặc biệt ở Bắc Bộ và Bắc Trung Bộ, báo hiệu các đợt rét đậm mùa đông sẽ giảm về số lượng nhưng gia tăng tính chất bất thường.

Thảo luận kết quả

Nguyên nhân dẫn đến sự sai khác cục bộ giữa mô phỏng và quan trắc xuất phát từ sự tương tác phức tạp giữa địa hình dốc đứng của dãy Trường Sơn với các hệ thống khí quyển quy mô lớn. Hiệu ứng phơn khô nóng ở sườn đông Trường Sơn trong mùa hè được mô hình RCM-CCAM khuếch đại rõ nét hơn so với mô hình toàn cầu, dẫn đến hiện tượng nhiệt độ bề mặt thiên cao nhẹ.

So sánh với các nghiên cứu trước đây của IPCC (2007) và các công trình nghiên cứu khí hậu trong nước, xu thế biến đổi cực trị trong luận văn hoàn toàn tương đồng. Tốc độ gia tăng của nhiệt độ cực tiểu nhanh hơn nhiệt độ cực đại (thực tế lịch sử tại Việt Nam ghi nhận nhiệt độ cực tiểu tăng khoảng 0,9°C/thập kỷ trong khi cực đại chỉ tăng khoảng 0,1°C/thập kỷ), dẫn đến sự thu hẹp biên độ dao động nhiệt độ ngày đêm.

Để trực quan hóa các kết quả này, dữ liệu biến trình 12 tháng tại 7 vùng khí hậu được thể hiện tối ưu qua các biểu đồ đường kép so sánh chuỗi số liệu quan trắc trạm với kết quả CCAM. Đồng thời, cấu trúc bảng số liệu tổng hợp chuẩn sai nhiệt độ cực trị và bản đồ phân bố không gian 2D với thang màu chi tiết giúp làm nổi bật các tâm điểm gia tăng số ngày nắng nóng Tx35 và biến động mưa lớn R50 trên từng khu vực địa lý.

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

Dựa trên các kết quả dự tính cực trị khí hậu, luận văn đưa ra 4 nhóm giải pháp chiến lược:

Thứ nhất, nâng cấp và tích hợp hệ thống mô hình dự báo cực trị khí hậu phân giải cao. Cục Khí tượng Thủy văn phối hợp cùng Trung tâm Dự báo Khí tượng Thủy văn Quốc gia cần đưa mô hình CCAM độ phân giải 10–25 km vào quy trình nghiệp vụ dự báo hạn mùa, phấn đấu nâng cao độ chính xác cảnh báo các đợt sóng nhiệt gay gắt và mưa lớn cực đoan lên trên 85% trong giai đoạn 2026–2030.

Thứ hai, điều chỉnh và chuyển đổi cơ cấu mùa vụ, giống cây trồng thích ứng với nền nhiệt gia tăng. Bộ Nông nghiệp và Phát triển Nông thôn cùng chính quyền các tỉnh Tây Nguyên và Nam Bộ cần chủ động đưa vào canh tác các giống lúa, hoa màu có khả năng chịu nhiệt và hạn hán cao, thích ứng với ngưỡng tăng nhiệt độ 1,5°C - 2,5°C, hoàn thành việc chuyển đổi cơ cấu cây trồng trước năm 2030 nhằm giảm thiểu 20-30% rủi ro thiệt hại nông nghiệp.

Thứ ba, gia cố và nâng cấp hạ tầng thoát nước, đê điều ứng phó với mưa cực đoan và lũ quét. Bộ Xây dựng và Ủy ban nhân dân các tỉnh ven biển miền Trung cần rà soát, nâng chuẩn thiết kế hệ thống tiêu thoát nước đô thị và hồ chứa thủy điện, bảo đảm năng lực ứng phó với lượng mưa ngày cực đại vượt trên 200 mm/ngày, hoàn thành các công trình thoát lũ trọng điểm trước năm 2035.

Thứ tư, xây dựng và ban hành bản đồ phân vùng rủi ro thiên tai khí hậu cực đoan chi tiết cấp huyện. Viện Khoa học Khí tượng Thủy văn và Môi trường cần chủ trì hoàn thiện bộ cơ sở dữ liệu số hóa về nguy cơ nắng nóng và hạn hán trước năm 2028, làm căn cứ bắt buộc khi thẩm định các quy hoạch phát triển kinh tế - xã hội vùng.

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

Luận văn là tài liệu tham khảo học thuật và ứng dụng giá trị cao cho 4 nhóm đối tượng chính:

Thứ nhất, học viên cao học, nghiên cứu sinh và sinh viên chuyên ngành Khí tượng, Thủy văn, Biến đổi khí hậu: Tài liệu cung cấp quy trình chuẩn mực về phương pháp cài đặt, hiệu chỉnh và chạy mô hình khí hậu khu vực CCAM, phương pháp xử lý dữ liệu lưới tái phân tích ERA-40 và cách tính toán 27 chỉ số cực trị khí hậu.

Thứ hai, dự báo viên và cán bộ kỹ thuật tại các cơ quan khí tượng thủy văn: Nắm bắt bản chất sai số hệ thống của mô hình động lực (độ lệch thiên dương nhiệt độ, mức độ hụt mưa mùa thu ở miền Trung) để hiệu chỉnh sai số thống kê trong các bản tin dự báo khí hậu hạn dài.

Thứ ba, nhà hoạch định chính sách và cán bộ quản lý môi trường tại các Bộ, ngành, địa phương: Tiếp cận bộ số liệu kịch bản biến đổi cực trị khí hậu A1B và A2 để xây dựng Kế hoạch hành động ứng phó với biến đổi khí hậu cấp tỉnh và quy hoạch không gian lãnh thổ bền vững.

Thứ tư, kỹ sư thiết kế công trình thủy lợi, giao thông, xây dựng đô thị: Khai thác số liệu dự tính về lượng mưa ngày cực đại Rx và lượng mưa 5 ngày lớn nhất Rx5day để tính toán khẩu độ công trình thoát nước, đập thủy điện và cao trình nền đường, giảm thiểu nguy cơ sạt lở ngập úng.

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

Ưu điểm lớn nhất của mô hình khí quyển bảo giác lập phương CCAM là gì?
Mô hình CCAM sử dụng lưới bảo giác chiếu từ khối lập phương lên mặt cầu, giải quyết triệt để bài toán biến dạng lưới và bước thời gian tích phân cực ngắn ở hai cực của mô hình kinh vĩ độ truyền thống. Điều này cho phép mô hình tích hợp phương pháp bán Lagrange để chi tiết hóa động lực phân giải cao 25 km cho khu vực Việt Nam một cách ổn định và tiết kiệm tài nguyên tính toán.

Các chỉ số Tx35 và Tm15 mang ý nghĩa gì trong nghiên cứu khí hậu cực đoan?
Chỉ số Tx35 đại diện cho số ngày trong tháng có nhiệt độ cao nhất vượt 35°C, phản ánh tần suất xuất hiện các đợt nắng nóng gay gắt ảnh hưởng trực tiếp đến sức khỏe và lưới điện. Ngược lại, chỉ số Tm15 đo lường số ngày có nhiệt độ thấp nhất dưới 15°C, là chỉ tiêu nhận diện các đợt rét đậm, rét hại gây thiệt hại nặng nề cho nông nghiệp và chăn nuôi tại miền Bắc.

Mô hình CCAM mô phỏng nhiệt độ và lượng mưa thời kỳ chuẩn 1980–1999 chính xác đến mức nào?
Mô hình tái tạo rất chuẩn xác biến trình năm của nhiệt độ cực tiểu tại 6/7 vùng khí hậu với sai số chỉ 1°C - 2°C so với trạm quan trắc. Về lượng mưa, mô hình phản ánh đúng mùa mưa ẩm nhưng có xu hướng ước tính thấp hơn thực tế khoảng 380 mm trong các tháng mưa đỉnh điểm tháng 10 tại khu vực miền Trung.

Kịch bản phát thải A1B và A2 phản ánh những viễn cảnh nào trong tương lai?
Kịch bản A1B mô tả thế giới phát triển công nghệ vượt bậc với sự cân bằng giữa năng lượng hóa thạch và năng lượng tái tạo, đại diện cho mức phát thải trung bình. Kịch bản A2 phản ánh thế giới bảo hộ khu vực, gia tăng dân số liên tục và chuyển dịch công nghệ chậm, đại diện cho kịch bản phát thải cao với mức tăng nhiệt độ và cực trị khí hậu khốc liệt nhất vào cuối thế kỷ XXI.

Khu vực nào ở Việt Nam ghi nhận sự gia tăng nhiệt độ mạnh nhất trong các thập kỷ qua?
Theo số liệu quan trắc 50 năm, khu vực Tây Bắc, Đông Bắc và Đồng bằng Bắc Bộ có tốc độ tăng nhiệt độ mùa đông nhanh nhất cả nước, đạt mức tăng từ 1,3°C đến 1,5°C/50 năm. Tốc độ ấm lên của nhiệt độ cực tiểu tại các khu vực này diễn ra nhanh vượt trội so với nhiệt độ cực đại.

Kết luận

Luận văn đã thực hiện trọn vẹn mục tiêu đánh giá tác động của biến đổi khí hậu đến các cực trị nhiệt độ và lượng mưa tại Việt Nam thông qua mô hình động lực CCAM. Năm đóng góp khoa học then chốt bao gồm:

  • Khẳng định độ tin cậy của mô hình CCAM trong việc mô phỏng hoàn lưu khí quyển mực 850 hPa và trường nhiệt ẩm tại Việt Nam với sai số nhiệt độ bề mặt chỉ từ 1°C đến 3°C.
  • Chứng minh khả năng tái tạo chuẩn xác biến trình năm của nhiệt độ cực trị tại 73 trạm quan trắc trên 7 vùng khí hậu, đặc biệt là nhiệt độ cực tiểu với độ lệch chỉ 1°C - 2°C.
  • Chỉ ra đặc tính suy giảm có hệ thống của lượng mưa mô phỏng trong các tháng mưa lũ lớn ở miền Trung, cung cấp cơ sở để hiệu chỉnh mô hình số trị.
  • Dự tính định lượng sự gia tăng vượt bậc của số ngày nắng nóng Tx35 và sự suy giảm của số ngày rét Tm15 trong 3 giai đoạn 2011–2040, 2041–2070 và 2071–2100 theo kịch bản A1B và A2.
  • Cung cấp bộ cơ sở dữ liệu số hóa dạng lưới độ phân giải 25 km phục vụ trực tiếp công tác xây dựng kịch bản biến đổi khí hậu quốc gia.

Bước phát triển tiếp theo trong giai đoạn 2026–2030 cần tập trung nâng cao độ phân giải mô hình xuống dưới 10 km và tích hợp các sơ đồ tham số hóa mây tiên tiến. Các nhà nghiên cứu, kỹ sư và cơ quan quản lý hãy chủ động khai thác ngay các dữ liệu dự tính cực trị của luận văn để tối ưu hóa quy hoạch hạ tầng và nâng cao năng lực thích ứng với biến đổi khí hậu bền vững.