Tổng quan nghiên cứu

Biến đổi khí hậu toàn cầu đang gia tăng nhanh chóng với tần suất và cường độ các hiện tượng thời tiết cực đoan ngày càng phức tạp. Tại Việt Nam, các đợt thiên tai cực đoan hàng năm gây thiệt hại ước tính khoảng 1% đến 1,5% GDP, đe dọa nghiêm trọng đến đời sống dân sinh và an ninh lương thực. Trong bối cảnh đó, bài toán dự báo khí hậu hạn mùa với độ dẫn trước từ 1 đến 6 tháng đóng vai trò then chốt trong công tác chủ động phòng tránh rủi ro. Tuy nhiên, các phương pháp dự báo thống kê truyền thống bộc lộ nhiều hạn chế do chỉ dựa trên chuỗi số liệu quá khứ, không mô tả đầy đủ các quá trình vật lý khí quyển và khó nắm bắt các pha chuyển đổi thời tiết đột biến.

Luận văn tập trung giải quyết vấn đề ứng dụng mô hình số trị động lực khu vực tiên tiến nhằm thử nghiệm dự báo hạn mùa các trường khí tượng trung bình tháng và một số chỉ số khí hậu cực đoan tại Việt Nam. Mục tiêu cụ thể là xây dựng quy trình công nghệ dự báo hạn mùa bằng mô hình clWRF, đánh giá sai số mô phỏng và kiểm chứng hiệu quả của phương pháp hiệu chỉnh thống kê.

Phạm vi nghiên cứu bao trùm toàn bộ lãnh thổ Việt Nam và khu vực Biển Đông, với miền tính 144x130 điểm lưới ở độ phân giải 36km. Chuỗi số liệu đánh giá bao gồm 20 năm mô phỏng khí hậu quá khứ giai đoạn 1981-2000 cùng 58 trường hợp chạy thử nghiệm dự báo trong giai đoạn từ tháng 1 năm 2012 đến tháng 5 năm 2013 tại 77 trạm quan trắc. Kết quả nghiên cứu có ý nghĩa thực tiễn to lớn khi giúp giảm sai số tuyệt đối trung bình nhiệt độ từ trên 3,5 độ C xuống dưới 1,2 độ C, thiết lập tiền đề khoa học vững chắc để nâng cao độ tin cậy của các bản tin cảnh báo thiên tai sớm tại Việt Nam.

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 trên nền tảng lý thuyết hoàn lưu chung khí quyển và kỹ thuật hạ thấp quy mô động lực từ các mô hình toàn cầu xuống quy mô khu vực. Mô hình số trị được lựa chọn là phiên bản cải tiến clWRF phát triển từ mô hình WRF 3.1. Điểm ưu việt của clWRF là tích hợp các mô-đun phục vụ nghiên cứu khí hậu dài hạn, đặc biệt là sơ đồ bức xạ sóng dài và sóng ngắn CAM có khả năng cập nhật nồng độ khí nhà kính theo các kịch bản phát thải SRES và RCP cùng sự biến đổi ôzôn theo từng tháng.

Bên cạnh đó, khung lý thuyết sử dụng lý thuyết giá trị cực trị dựa trên hàm phân bố xác suất và hàm mật độ tích lũy để xác định các ngưỡng phân vị chuẩn. Bốn khái niệm chính được định nghĩa chặt chẽ bao gồm:

  • Ngưỡng phân vị cực đoan nhiệt độ: Phân vị thứ 10 và thứ 90 của nhiệt độ cực tiểu ngày (TN10, TN90) và nhiệt độ cực đại ngày (TX10, TX90).
  • Ngưỡng phân vị cực đoan lượng mưa: Phân vị thứ 95 và thứ 99 của lượng mưa ngày (R95, R99).
  • Sai số hệ thống mô hình: Độ lệch cố định do cấu trúc tham số hóa vật lý và địa hình gây ra trong quá trình tích phân dài hạn.
  • Hạ thấp quy mô động lực: Quá trình lồng ghép điều kiện biên quy mô lớn từ mô hình toàn cầu vào mô hình khu vực độ phân giải cao.

Hệ thống vật lý trong mô hình bao gồm sơ đồ vi vật lý băng 3 lớp WSM, mô hình bề mặt đất Noah và sơ đồ tham số hóa đối lưu Kain-Fritsch.

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

Cỡ mẫu nghiên cứu được thiết lập trên mạng lưới 77 trạm khí tượng bề mặt đại diện cho 7 vùng khí hậu đặc trưng của Việt Nam: Tây Bắc (B1), Đông Bắc (B2), Đồng bằng Bắc Bộ (B3), Bắc Trung Bộ (B4), Nam Trung Bộ (N1), Tây Nguyên (N2) và Nam Bộ (N3). Đối với bài toán cực đoan, cỡ mẫu gồm 67 trạm có chuỗi số liệu ngày liên tục và đồng nhất. Phương pháp chọn mẫu phân tầng theo vùng địa lý đảm bảo phản ánh đầy đủ tính đa dạng của địa hình đồi núi, đồng bằng và ven biển.

Nguồn dữ liệu đầu vào bao gồm:

  • Số liệu tái phân tích khí hậu CFSR độ phân giải 0,5 độ kinh vĩ từ NCEP giai đoạn 1980-2000 (với năm 1980 dùng làm thời gian khởi động mô hình - spin-up).
  • Dữ liệu dự báo toàn cầu thời gian thực CFSv1 tại hạn dự báo 6 tháng với chu kỳ cập nhật 7 ngày một lần.
  • Bộ số liệu quan trắc thực tế 20 năm (1981-2000) và chuỗi quan trắc 17 tháng thực nghiệm (tháng 1/2012 đến tháng 5/2013).

Lý do lựa chọn phương pháp phân tích kết hợp giữa hạ thấp quy mô động lực và hiệu chỉnh thống kê delta là vì các mô hình động lực luôn tồn tại độ lệch khí hậu riêng. Phương pháp hiệu chỉnh dựa trên độ lệch trung bình 20 năm giữa mô phỏng và quan trắc giúp loại bỏ phần lớn sai số hệ thống mà không làm mất đi các tín hiệu động lực quy mô lớn. Các chỉ số thống kê định lượng gồm sai số trung bình (ME), sai số tuyệt đối trung bình (MAE) và sai số bình phương trung bình quân phương (RMSE) được tính toán chi tiết cho từng hạn dự báo từ 1 đến 6 tháng.

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

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

Thứ nhất, mô hình clWRF nguyên bản tồn tại sai số hệ thống thiên lạnh rõ rệt trên phạm vi toàn quốc. Đối với nhiệt độ trung bình bề mặt mực 2m (T2m) và nhiệt độ cực tiểu trung bình tháng (Tm), giá trị dự báo chưa hiệu chỉnh luôn thấp hơn quan trắc từ 3,0 độ C đến 5,0 độ C. Sai số ME dao động trong khoảng -5,0 độ C đến -2,0 độ C, trong khi sai số MAE và RMSE trước hiệu chỉnh duy trì ở mức cao, phổ biến từ 2,5 độ C đến 3,8 độ C tại hầu hết các hạn dự báo.

Thứ hai, phương pháp hiệu chỉnh thống kê đã cải thiện vượt bậc chất lượng dự báo với tỷ lệ giảm sai số đạt trên 60%. Sau khi hiệu chỉnh, chênh lệch nhiệt độ T2m giữa mô hình và thực tế giảm xuống chỉ còn khoảng 1,0 độ C đến 1,5 độ C. Giá trị ME sau hiệu chỉnh dao động quanh mức -2,5 độ C đến +0,5 độ C (phổ biến từ -1,0 độ C đến -1,5 độ C), và RMSE giảm mạnh xuống khoảng 0,8 độ C đến 1,8 độ C.

Thứ ba, độ chính xác của mô hình thể hiện sự phân hóa sâu sắc theo mùa và vùng khí hậu:

  • Về mặt thời gian: Mô hình dự báo tốt nhất trong các tháng chính hè (tháng 6, 7, 8 với nền nhiệt ổn định khoảng 27,0 độ C đến 28,5 độ C) và chính đông. Ngược lại, sai số lớn nhất xuất hiện ở các tháng chuyển mùa (tháng 4 và tháng 10), với MAE sau hiệu chỉnh vẫn ở mức 2,0 độ C đến 2,5 độ C.
  • Về mặt không gian: Khu vực Đồng bằng Bắc Bộ (B3) và Bắc Trung Bộ (B4) có sai số mùa đông trước hiệu chỉnh lên tới 5,0 độ C, sau hiệu chỉnh giảm còn 2,0 độ C đến 2,5 độ C. Khu vực Nam Bộ (N3) thể hiện tính ổn định cao nhất, sai số ME chỉ dao động trong khoảng -3,0 độ C đến +1,0 độ C và ít biến động theo hạn dự báo.

Thứ tư, đối với các chỉ số cực đoan, mô hình tái tạo tương đối tốt xu thế không gian của các ngày có nhiệt độ cực đại vượt ngưỡng phân vị 90% (TX90) và nhiệt độ cực tiểu dưới ngưỡng phân vị 10% (TN10). Dự báo lượng mưa tháng và mưa cực đoan (R95, R99) nắm bắt được quy luật phân bố chung nhưng có mức độ dao động sai số lớn hơn so với nhiệt độ.

Thảo luận kết quả

Nguyên nhân vật lý dẫn đến hiện tượng thiên lạnh của mô hình clWRF bắt nguồn từ sự tương tác phức tạp giữa sơ đồ bức xạ CAM, mô hình đất Noah và sơ đồ đối lưu Kain-Fritsch. Sơ đồ tham số hóa đối lưu có xu hướng kích hoạt mây nhiều hơn thực tế tại các vùng nhiệt đới ẩm, làm giảm lượng bức xạ mặt trời xuống bề mặt đất. Đồng thời, độ ẩm đất bị ước tính cao dẫn đến gia tăng bốc thoát hơi nước và làm giảm nhiệt độ bề mặt. Kết quả này hoàn toàn tương đồng với các nghiên cứu khí hậu khu vực bằng mô hình RegCM và MM5 tại Đông Nam Á của các nhóm tác giả trong nước và quốc tế.

Khi biểu diễn trên biểu đồ ma trận hạn dự báo từ 1 đến 6 tháng, sai số MAE của nhiệt độ biến thiên rất ít theo độ dẫn thời gian mà phụ thuộc chủ yếu vào tháng đích được dự báo. Điều này chứng minh rằng việc dự báo các tháng mùa hè từ hạn 6 tháng trước (ví dụ đứng từ tháng 1 dự báo cho tháng 7) vẫn đạt độ chính xác xấp xỉ so với hạn 1 tháng (đứng từ tháng 6 dự báo cho tháng 7). Trên các bảng tổng hợp số liệu 17 tháng thử nghiệm, đồ thị so sánh giữa đường quan trắc, đường dự báo thô và đường dự báo hiệu chỉnh cho thấy đường hiệu chỉnh bám rất sát biến trình năm của nhiệt độ thực tế, khẳng định tính khả thi của quy trình công nghệ trong nghiệp vụ.

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

Nhằm hoàn thiện và đưa công nghệ dự báo hạn mùa số trị vào ứng dụng thực tiễn đạt hiệu quả cao, bốn giải pháp trọng tâm được đề xuất như sau:

  • Tinh chỉnh và thử nghiệm đa dạng hóa các sơ đồ vật lý mô hình: Trung tâm Dự báo Khí tượng Thủy văn Quốc gia phối hợp với các viện nghiên cứu tiến hành thử nghiệm kết hợp sơ đồ vi vật lý WDM6 với sơ đồ tham số hóa đối lưu Betts-Miller-Janjic hoặc Grell-Freitas. Mục tiêu giảm sai số thô ban đầu của nhiệt độ và lượng mưa từ 15% đến 20%, thực hiện trong giai đoạn 2026-2027.
  • Ứng dụng các thuật toán hiệu chỉnh sai số phi tuyến tính và học máy: Nhóm nghiên cứu mô hình hóa khí hậu triển khai các kỹ thuật hiệu chỉnh phân vị đa chiều (Quantile Mapping) kết hợp thuật toán Random Forest để xử lý sai số cục bộ trong các tháng chuyển tiếp mùa (tháng 4 và tháng 10). Mục tiêu đưa sai số MAE của các tháng chuyển mùa xuống dưới 1,2 độ C trước quý IV năm 2027.
  • Nâng cao độ phân giải không gian miền tính xuống quy mô 10-12km: Bộ Tài nguyên và Môi trường phê duyệt dự án đầu tư nâng cấp hạ tầng tính toán hiệu năng cao (HPC), nâng độ phân giải lưới từ 36km xuống 10km. Mục tiêu tăng độ chính xác dự báo lượng mưa cực trị (R95, R99) tại các khu vực địa hình chia cắt phức tạp thêm 25% trong lộ trình 2027-2028.
  • Xây dựng hệ thống dự báo tổ hợp đa mô hình: Các cơ quan chuyên trách thiết lập hệ thống tích hợp kết quả từ clWRF, RegCM4 và hệ thống dự báo toàn cầu CFSv2. Mục tiêu đạt độ tin cậy trên 80% trong các bản tin cảnh báo sớm các đợt nắng nóng diện rộng và không khí lạnh kỷ lục trước từ 2 đến 3 tháng, vận hành thử nghiệm từ năm 2028.

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

Nội dung luận văn mang lại giá trị học thuật và ứng dụng chuyên sâu cho bốn nhóm đối tượng chính:

  • Cán bộ và chuyên viên dự báo khí tượng thủy văn: Khai thác chi tiết quy trình công nghệ tích hợp giữa mô hình động lực clWRF và phương pháp hiệu chỉnh delta 20 năm để ứng dụng trực tiếp vào nghiệp vụ dự báo mùa tại các Đài Khí tượng Thủy văn khu vực.
  • Giảng viên và học viên cao học, nghiên cứu sinh ngành Khoa học Trái đất: Sử dụng khung lý thuyết, phương pháp thiết lập miền tính 36km và hệ thống đánh giá sai số 77 trạm làm tài liệu tham khảo chuẩn mực trong giảng dạy và nghiên cứu mô phỏng khí hậu.
  • Cơ quan quản lý nhà nước về phòng chống thiên tai: Sử dụng các phát hiện về chu kỳ xuất hiện và quy luật phân bố của các hiện tượng cực đoan để xây dựng phương án ứng phó rủi ro thiên tai cấp quốc gia và cấp tỉnh trước từ 3 đến 6 tháng.
  • Chuyên gia quy hoạch nông nghiệp và quản lý tài nguyên nước: Tham khảo các phân tích dị thường nhiệt ẩm hạn mùa để lập kế hoạch điều tiết lưu lượng hồ chứa thủy điện và bố trí cơ cấu mùa vụ nông nghiệp tối ưu cho 7 vùng sinh thái.

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

Mô hình clWRF có ưu điểm gì vượt trội so với mô hình WRF tiêu chuẩn trong dự báo khí hậu? Mô hình clWRF tích hợp sơ đồ bức xạ CAM tiên tiến cho phép cập nhật biến trình ôzôn hàng tháng và nồng độ khí nhà kính theo các kịch bản SRES/RCP. Điều này giúp mô hình duy trì tính ổn định trong các phép tích phân thời gian dài hạn từ 6 tháng đến hàng thập kỷ mà phiên bản WRF thông thường không hỗ trợ tối ưu.

Tại sao mô hình clWRF lại xuất hiện sai số thiên lạnh khi áp dụng tại Việt Nam? Sai số thiên lạnh phát sinh do sự kết hợp giữa sơ đồ đối lưu Kain-Fritsch và mô hình đất Noah ước tính độ che phủ mây và độ ẩm đất cao hơn thực tế tại vùng nhiệt đới ẩm. Hiện tượng này làm giảm bức xạ sóng ngắn tới bề mặt, khiến nhiệt độ dự báo thô thấp hơn quan trắc từ 3,0 độ C đến 5,0 độ C.

Phương pháp hiệu chỉnh thống kê trong nghiên cứu vận hành theo nguyên lý nào? Nghiên cứu sử dụng số liệu mô phỏng khí hậu 20 năm (1981-2000) trừ đi số liệu quan trắc cùng thời kỳ tại 77 trạm để xác định độ lệch chuẩn cố định cho từng tháng. Giá trị chênh lệch này được cộng trực tiếp vào kết quả dự báo thời gian thực, giúp giảm sai số MAE của nhiệt độ xuống dưới 1,5 độ C.

Độ chính xác của dự báo thay đổi như thế nào từ hạn 1 tháng đến hạn 6 tháng? Kết quả thực nghiệm cho thấy chất lượng dự báo nhiệt độ trung bình tháng ít bị suy giảm theo độ dài hạn dự báo. Sai số phụ thuộc chủ yếu vào đặc điểm mùa của tháng đích (tháng chính hè và chính đông dự báo chính xác hơn tháng chuyển mùa) hơn là khoảng cách thời gian đứng làm dự báo.

Khả năng cảnh báo các hiện tượng thời tiết cực đoan của mô hình đạt mức độ nào? Mô hình mô phỏng tương đối chính xác hình thế phân bố không gian của các đợt nắng nóng diện rộng (TX90) và các đợt rét hại (TN10). Tuy nhiên, đối với cực đoan lượng mưa (R95, R99), mô hình mới nắm bắt được xu thế chung chứ chưa định lượng hoàn toàn chính xác tâm mưa lớn cục bộ do hạn chế độ phân giải 36km.

Kết luận

  • Luận văn đã xây dựng thành công quy trình thử nghiệm dự báo hạn mùa từ 1 đến 6 tháng tại Việt Nam bằng mô hình khí hậu động lực khu vực clWRF trên lưới phân giải 36km.
  • Thiết lập hoàn chỉnh bộ cơ sở dữ liệu mô phỏng khí hậu 20 năm (1981-2000) tại 77 trạm quan trắc, cung cấp nguồn dữ liệu chuẩn cho việc xác định sai số hệ thống và ngưỡng cực trị.
  • Chứng minh phương pháp hiệu chỉnh thống kê delta giúp loại bỏ triệt để sai số thiên lạnh, nâng cao độ chính xác dự báo nhiệt độ với mức giảm sai số toàn diện trên 60%.
  • Làm rõ tính quy luật của sai số theo thời gian và không gian, trong đó các tháng chính hè, chính đông và khu vực Nam Bộ đạt độ tin cậy dự báo cao nhất.
  • Đánh giá khả thi năng lực dự báo các chỉ số cực đoan nhiệt độ và lượng mưa hạn mùa, tạo bước đột phá trong công tác cảnh báo thiên tai sớm.

Đóng góp lớn nhất của công trình là đặt nền móng thực nghiệm vững chắc cho bài toán dự báo mùa số trị kết hợp hiệu chỉnh thống kê tại Việt Nam. Lộ trình phát triển tiếp theo giai đoạn 2026-2028 tập trung vào việc nâng cao độ phân giải không gian lên 10km và tích hợp trí tuệ nhân tạo. Kính mời các nhà khoa học, chuyên gia khí tượng và cơ quan quản lý tiếp tục khai thác, phát triển nguồn tư liệu quý giá này nhằm nâng cao năng lực chủ động ứng phó biến đổi khí hậu quốc gia.