Tổng quan nghiên cứu

Biến đổi khí hậu toàn cầu đang thúc đẩy sự gia tăng nghiêm trọng cả về tần suất, cường độ và thời gian kéo dài của các hiện tượng thời tiết cực đoan. Theo các báo cáo khí hậu quốc tế, các đợt nắng nóng cực đoan chiếm tới 75% tổng số ca tử vong liên quan đến thiên tai thời tiết đối với các nhóm dân số dễ bị tổn thương. Tại Việt Nam, Kịch bản Biến đổi Khí hậu và Nước biển dâng công bố năm 2016 đã chỉ ra rằng số ngày nắng nóng có xu thế tăng từ 20 đến 40 ngày trên cả nước, trong đó dải đất Bắc Trung Bộ là khu vực ghi nhận tốc độ gia tăng nhanh nhất và khắc nghiệt nhất. Các kỷ lục nhiệt độ lịch sử liên tiếp được thiết lập tại khu vực này, tiêu biểu như 42,7°C tại Tương Dương và 43,4°C tại Hương Khê, gây ảnh hưởng nặng nề đến sản xuất nông nghiệp, an ninh năng lượng và sức khỏe cộng đồng.

Xuất phát từ thực tiễn trên, đề tài tập trung nghiên cứu và đánh giá năng lực dự báo hạn mùa từ 1 đến 5 tháng đối với các chỉ số cực đoan nhiệt độ mùa hè (tháng 6, 7 và 8) trên khu vực Bắc Trung Bộ bằng mô hình khí hậu khu vực RegCM phiên bản 4.6. Nghiên cứu khai thác chuỗi số liệu dự báo lại dài hạn trong 30 năm (giai đoạn 1983-2012) từ mô hình toàn cầu CFS của NCEP (Hoa Kỳ) kết hợp bộ số liệu quan trắc tại 25 trạm khí tượng phân bố khắp 6 tỉnh Bắc Trung Bộ. Đồng thời, đề tài tiến hành thử nghiệm phương pháp hiệu chỉnh sai số hệ thống cho giai đoạn 2013-2014. Kết quả nghiên cứu tạo tiền đề khoa học vững chắc nhằm nâng cao độ tin cậy của các bản tin cảnh báo sớm, góp phần giảm thiểu rủi ro thiên tai khí hậu cho khu vực miền Trung.

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 thủy nhiệt động lực học khí quyển và phương pháp hạ quy mô động lực (dynamical downscaling). Mô hình khí hậu khu vực RegCM4.6 giải hệ phương trình nguyên thủy thủy tĩnh nén trong hệ tọa độ sigma-p, tích hợp các sơ đồ tham số hóa vật lý hiện đại nhằm mô phỏng chi tiết các quá trình trao đổi năng lượng, bức xạ và nhiệt ẩm.

Đồng thời, nghiên cứu vận dụng lý thuyết nhiệt động học địa phương để mô tả hiệu ứng phơn (Foehn) – nhân tố cốt lõi gây ra các đợt khô nóng khốc liệt tại Bắc Trung Bộ khi khối không khí gió mùa Tây Nam vượt qua dãy Trường Sơn Bắc bị mất ẩm và đốt nóng đoạn nhiệt khô với gradien nhiệt độ giảm theo độ cao đạt δd = -0,98°C/100m.

Khung đánh giá chỉ số cực đoan bám sát hướng dẫn của Tổ chức Khí tượng Thế giới (WMO) và nhóm chuyên gia ETCCDI, tập trung vào 3 nhóm khái niệm cốt lõi:

  • Các yếu tố cực trị nhiệt độ: Nhiệt độ tối cao trung bình tháng (Tx) và Nhiệt độ tối cao tuyệt đối (TXx).
  • Chỉ số tần suất cực đoan: Phần trăm số ngày trong tháng có nhiệt độ cực đại vượt ngưỡng phân vị 90% (TX90p).
  • Hiện tượng cực đoan nắng nóng: Số ngày nắng nóng (SNNN), nắng nóng diện rộng (SNNNDR), nắng nóng gay gắt diện rộng (SNNNGGDR) và số đợt nắng nóng kéo dài được xác định theo cả ngưỡng nhiệt cố định (35°C, 37°C) lẫn ngưỡng phân vị thực nghiệm tại từng trạm (P35, P37).

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

Nguồn dữ liệu đầu vào sử dụng sản phẩm dự báo lại toàn cầu CFSv2 với độ phân giải lưới 1° x 1° kinh vĩ, gồm 18 tầng áp suất theo phương thẳng đứng và nhiệt độ mặt nước biển OISST cập nhật liên tục từ NOAA. Chuỗi số liệu thực đo nhiệt độ tối cao ngày được thu thập đồng bộ tại 25 trạm quan trắc khí tượng trải dài từ Thanh Hóa đến Thừa Thiên Huế trong giai đoạn 32 năm (1983-2014). Cỡ mẫu nghiên cứu bao gồm hơn 73.600 quan trắc trạm-ngày, đảm bảo độ đại diện thống kê cao và phản ánh đầy đủ các pha dao động khí hậu quy mô lớn như ENSO.

Quy trình mô phỏng thiết lập hệ thống 2 miền tính lồng nhau theo phương thức offline: miền ngoài có độ phân giải 60 km (119 x 148 nút lưới) bao trùm các hoàn lưu chi phối khu vực châu Á, miền trong có độ phân giải 20 km (139 x 139 nút lưới) bao trọn lãnh thổ Việt Nam. Cấu hình vật lý gồm sơ đồ đối lưu Tiedtke, sơ đồ lớp biên Holtslag, sơ đồ bức xạ CCM3, sơ đồ giáng thủy SUBEX và mô hình bề mặt đất BATS1e.

Số liệu đầu ra của mô hình được nội suy về tọa độ 25 trạm và hiệu chỉnh độ lệch địa hình thực tế bằng công thức gradien đoạn nhiệt khô. Đánh giá kỹ năng dự báo sử dụng tập hợp các chỉ số thống kê toán học: Sai số trung bình (ME), Sai số tuyệt đối trung bình (MAE), Sai số quân phương (RMSE), Hệ số tương quan Pearson (r) và phân bố mật độ xác suất thực nghiệm. Cuối cùng, phương pháp hiệu chỉnh sai số hệ thống (BSC) được áp dụng độc lập cho giai đoạn thử nghiệm 2013-2014 nhằm cô lập và triệt tiêu các độ lệch hệ thống phát sinh từ cấu trúc số trị.

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

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

Quá trình phân tích chuỗi dữ liệu 30 năm (1983-2012) đã làm sáng tỏ nhiều đặc tính dự báo quan trọng của mô hình RegCM4.6 đối với chế độ nhiệt mùa hè tại Bắc Trung Bộ:

  • Năng lực mô phỏng nhiệt độ tối cao trung bình tháng (Tx): Độ chính xác của mô hình có sự phân hóa rõ nét theo từng tháng mùa hè. Trong tháng 6, mô hình đạt chất lượng dự báo tốt nhất với sai số MAE dao động từ 1,5°C đến 2,0°C và ME trung bình vùng xấp xỉ -0,5°C đến 1,1°C. Sang tháng 7, sai số MAE tăng lên mức 1,9°C đến 3,1°C đi kèm xu thế dự báo thiên cao với ME dương từ 1,6°C đến 2,9°C. Tháng 8 ghi nhận sai số lớn nhất khi MAE đạt từ 2,8°C đến 4,4°C và ME thiên cao từ 2,7°C đến 4,3°C trên toàn khu vực.
  • Biến động sai số theo hạn dự báo: Kết quả cho thấy xu hướng khác biệt so với dự báo thời tiết thông thường. Sai số MAE và RMSE của Tx trong các tháng 7 và 8 có xu hướng giảm khi hạn dự báo tăng từ 1 tháng lên 4-5 tháng. Cụ thể trong tháng 8, sai số MAE trung bình toàn vùng giảm từ 4,4°C ở hạn 1 tháng xuống còn 2,8°C ở hạn 5 tháng.
  • Tương quan chuỗi thời gian: Hệ số tương quan r giữa quan trắc và mô hình cho toàn khu vực đạt giá trị từ 0,3 đến 0,6. Tháng 7 đạt tương quan cao nhất (r = 0,4 - 0,6), tiếp đến là tháng 8 (r = 0,4 - 0,5) và tháng 6 (r = 0,3 - 0,4). Mức độ tương quan này phản ánh khả năng nắm bắt tương đối đồng nhất xu thế không gian của hoàn lưu nhiệt độ mùa hè.
  • Hiệu quả hiệu chỉnh độ cao địa hình: Việc áp dụng gradien đoạn nhiệt khô (-0,98°C/100m) đã cải thiện trực tiếp số liệu đầu vào tại các trạm vùng núi cao. Điển hình tại trạm Tương Dương với chênh lệch độ cao mô hình so với thực tế là 476,8m, giá trị Tx được hiệu chỉnh bù nhiệt lên tới 4,67°C; trạm Hồi Xuân chênh lệch 300,4m được hiệu chỉnh 2,94°C.

Thảo luận kết quả

Nguyên nhân dẫn đến sai số thiên cao đáng kể trong tháng 7 và tháng 8 bắt nguồn từ sự tương tác phi tuyến phức tạp giữa sơ đồ bức xạ CCM3, sơ đồ đất BATS1e và hiệu ứng gió phơn địa phương. Khi mô phỏng bề mặt đệm trong điều kiện độ ẩm đất suy giảm vào nửa cuối mùa hè, mô hình có xu hướng đánh giá quá mức thông lượng nhiệt hiện và tích tụ nhiệt lượng tại lớp biên khí quyển, dẫn đến nhiệt độ cực đại bị đẩy lên cao hơn thực tế từ 2°C đến 5°C tại vùng đồng bằng ven biển Thanh Hóa và Nghệ An.

Kết quả này hoàn toàn tương đồng với các công bố khoa học quốc tế trong khuôn khổ dự án CORDEX Đông Nam Á khi sử dụng mô hình RegCM4 kết hợp dữ liệu tái phân tích, vốn cũng chỉ ra độ nhạy cao của sơ đồ đối lưu đối với trường nhiệt cực trị. Để minh họa trực quan, các trường dữ liệu sai số được thể hiện hiệu quả nhất thông qua hệ thống bản đồ phân bố không gian ME, đồ thị tần suất tích lũy sai số MAE và đồ thị phân tán điểm (scatter plot) kết hợp đường hồi quy thực nghiệm. Việc đối sánh giữa phân bố tần suất lý thuyết và quan trắc giúp phát hiện chính xác các khoảng giá trị mà mô hình bị lệch pha, làm căn cứ xác lập hàm truyền hiệu chỉnh tối ưu.

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

Nhằm hoàn thiện quy trình nghiệp vụ và ứng dụng hiệu quả sản phẩm mô hình số trị vào công tác dự báo khí hậu phục vụ phát triển kinh tế - xã hội, các giải pháp trọng tâm được đề xuất bao gồm:

  • Chuẩn hóa quy trình hiệu chỉnh sai số hệ thống (BSC) đa biến: Trung tâm Dự báo Khí tượng Thủy văn Quốc gia cần tích hợp mô-đun hiệu chỉnh BSC tự động vào chuỗi công nghệ dự báo hạn mùa, phấn đấu giảm sai số MAE của dự báo nhiệt độ tối cao xuống dưới 1,2°C tại 100% các trạm quan trắc trọng điểm trước năm 2025.
  • Nâng cấp độ phân giải và chuyển đổi sang lõi phi thủy tĩnh: Viện Khoa học Khí tượng Thủy văn và Biến đổi khí hậu chủ trì triển khai thử nghiệm phiên bản mô hình RegCM thế hệ mới với lưới tính độ phân giải cao ≤5 km cho khu vực miền Trung, hoàn thành đánh giá tham số hóa vật lý lớp biên trong giai đoạn 2024-2026.
  • Áp dụng bộ chỉ số ngưỡng cực đoan động dựa trên phân vị (P35, P37): Đài Khí tượng Thủy văn khu vực Bắc Trung Bộ cần chuyển đổi linh hoạt từ ngưỡng nhiệt tuyệt đối sang ngưỡng bách phân vị trạm trong xác định số ngày nắng nóng gay gắt, giúp tăng độ chính xác phát hiện các đợt sóng nhiệt cục bộ thêm ít nhất 20% trong giai đoạn 2024-2025.
  • Xây dựng hệ thống cảnh báo sớm rủi ro nhiệt nội mùa đến hạn mùa (S2S): Ban Chỉ đạo Quốc gia về Phòng chống Thiên tai phối hợp cùng các bộ, ngành liên quan ứng dụng sản phẩm dự báo hạn 1-3 tháng của mô hình để lập kế hoạch điều tiết nguồn nước tưới tiêu và điều độ phụ tải điện lưới quốc gia từ 30 đến 60 ngày trước các đợt cao điểm mùa hè.

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

Nội dung và kết quả thực nghiệm của công trình mang lại giá trị thực tiễn và học thuật cho 4 nhóm đối tượng chính:

  • Chuyên viên dự báo khí tượng thủy văn: Nắm vững phương pháp thiết lập mô hình khí hậu khu vực, kỹ thuật lồng ghép miền tính offline và thuật toán xử lý dữ liệu dự báo lại từ mô hình toàn cầu CFSv2 phục vụ công tác nghiệp vụ hàng ngày.
  • Học viên cao học và nghiên cứu sinh chuyên ngành Khí tượng, Biến đổi khí hậu: Khai thác làm tài liệu tham khảo chuẩn mực về phương pháp luận đánh giá chỉ số khí hậu cực đoan ETCCDI, kỹ thuật xử lý sai số địa hình và các công thức hiệu chỉnh thống kê sau mô hình.
  • Cơ quan quản lý nhà nước và hoạch định chính sách phòng chống thiên tai: Sử dụng các dẫn chứng định lượng về quy luật nhiệt độ cực trị Bắc Trung Bộ để xây dựng kế hoạch hành động ứng phó nắng nóng và hạn hán cấp tỉnh với tầm nhìn trung và dài hạn.
  • Doanh nghiệp trong lĩnh vực nông nghiệp, năng lượng và bảo hiểm: Ứng dụng các kịch bản dự báo hạn mùa nhiệt độ để tối ưu hóa lịch thời vụ gieo trồng, quản lý rủi ro cháy rừng và dự báo nhu cầu tiêu thụ điện năng trong các tháng cao điểm mùa khô.

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

1. Mô hình RegCM phiên bản 4.6 có những cải tiến kỹ thuật nổi bật nào so với các phiên bản tiền nhiệm?
RegCM4.6 nâng cấp toàn diện cấu trúc mã nguồn, tối ưu hóa khả năng chạy song song và bổ sung các sơ đồ đối lưu tiên tiến như Tiedtke, Kain-Fritsch bên cạnh sơ đồ hỗn hợp Grell-Emanuel. Đặc biệt, phiên bản này hỗ trợ tích hợp trực tiếp điều kiện biên phân giải cao từ hệ thống CFSv2 và cung cấp tùy chọn mô hình đất CLM4, nâng cao độ chính xác mô phỏng tương tác đất - khí quyển.

2. Vì sao mô hình RegCM4.6 lại có xu thế dự báo nhiệt độ tối cao thiên cao rõ rệt trong tháng 7 và tháng 8?
Nguyên nhân chủ yếu do sự kết hợp giữa hiệu ứng đốt nóng đoạn nhiệt phơn qua sườn đông Trường Sơn và hiện tượng phản hồi khô hạn của mô hình đất BATS1e. Trong các tháng giữa và cuối hè, mô hình tính toán lượng bốc thoát hơi nước thực tế thấp hơn thực tế, làm giảm mạnh thông lượng nhiệt ẩn và tăng thông lượng nhiệt hiện, khiến giá trị Tx dự báo cao hơn quan trắc từ 2,0°C đến 4,0°C.

3. Việc hiệu chỉnh độ cao địa hình theo đoạn nhiệt khô (-0,98°C/100m) mang lại ý nghĩa gì cho kết quả dự báo?
Địa hình mô hình ở lưới 20 km thường bị làm trơn đáng kể so với thực tế, tạo ra chênh lệch độ cao trạm lên tới hàng trăm mét. Áp dụng gradien đoạn nhiệt khô đặc trưng cho điều kiện gió phơn giúp loại bỏ sai số cơ học do độ cao, bù nhiệt chính xác từ 0,5°C đến 4,67°C tại các trạm miền núi như Tương Dương, Hồi Xuân, đảm bảo tính đồng nhất trước khi đánh giá kỹ năng mô hình.

4. Tại sao sai số MAE của dự báo nhiệt độ tối cao ở hạn 4-5 tháng lại có xu hướng thấp hơn hạn 1 tháng?
Hiện tượng này bắt nguồn từ quá trình điều chỉnh động lực học (spin-up) của mô hình khí hậu toàn cầu CFSv2. Ở các hạn dự báo ngắn (1 tháng), các dao động nhiễu động ban đầu chưa đạt trạng thái cân bằng thống kê hoàn toàn; trong khi ở hạn 4-5 tháng, các hoàn lưu quy mô lớn đạt độ ổn định cao hơn, giúp trường biên truyền vào RegCM4.6 ít bị đột biến cực đoan cục bộ.

5. Hiệu chỉnh sai số hệ thống (BSC) giúp cải thiện độ chính xác dự báo cụ thể ra sao trong thực nghiệm?
Phương pháp BSC sử dụng trường sai số ME và MAE được xây dựng từ chuỗi chuẩn 30 năm (1983-2012) để trừ khử độ lệch cố định cho từng trạm và từng hạn dự báo. Trong giai đoạn thử nghiệm 2013-2014, BSC đã thu hẹp đáng kể khoảng phân tán sai số, triệt tiêu độ lệch thiên cao của Tx tháng 7 và tháng 8, đưa đường hồi quy thực nghiệm tiệm cận sát đường hồi quy lý tưởng.

Kết luận

  • Luận văn đã thực hiện đánh giá toàn diện, có hệ thống năng lực dự báo hạn mùa từ 1 đến 5 tháng của mô hình RegCM4.6 đối với các cực đoan nhiệt độ mùa hè trên khu vực Bắc Trung Bộ dựa trên chuỗi số liệu 30 năm (1983-2012).
  • Xác lập bức tranh chi tiết về sai số mô hình với MAE dự báo Tx dao động từ 1,5°C đến 4,4°C; đồng thời phát hiện quy luật dự báo thiên thấp trong tháng 6 (ở hạn 3-5 tháng) và thiên cao rõ rệt trong tháng 7 và tháng 8 trên toàn vùng.
  • Khẳng định vai trò thiết yếu của việc hiệu chỉnh sai số địa hình bằng gradien đoạn nhiệt khô (-0,98°C/100m), giúp loại bỏ sai lệch nhiệt độ cơ học lên tới 4,67°C tại các trạm địa hình phức tạp.
  • Chứng minh tính hiệu quả vượt trội của phương pháp hiệu chỉnh sai số hệ thống (BSC) trong việc làm giảm sai số dự báo các chỉ số cực trị Tx, TXx và các đợt nắng nóng diện rộng trong giai đoạn kiểm định 2013-2014.
  • Đóng góp nguồn cơ sở dữ liệu khí hậu mô hình chuẩn hóa, mở ra hướng nghiên cứu ứng dụng dự báo tổ hợp đa mô hình kết hợp hạ quy mô thống kê nhằm phục vụ công tác cảnh báo sớm thiên tai tại Việt Nam.

Trong giai đoạn tiếp theo (2024-2026), các đơn vị nghiên cứu cần tiếp tục mở rộng thử nghiệm mô hình đối với các hiện tượng cực đoan lượng mưa và tích hợp trí tuệ nhân tạo (AI) vào khâu hậu xử lý sai số. Quý độc giả, các nhà nghiên cứu và cơ quan chuyên môn quan tâm có thể khai thác trực tiếp phương pháp luận và bộ tham số thực nghiệm của luận văn để phục vụ công tác nghiên cứu và triển khai ứng dụng thực tế.