Bối cảnh và vấn đề nghiên cứu

Trong công tác dự báo thời tiết bằng phương pháp số trị, mô hình nghiên cứu và dự báo thời tiết WRF (Weather Research and Forecasting) là một công cụ hiện đại được ứng dụng rộng rãi tại các cơ quan nghiên cứu và cơ quan khí tượng nghiệp vụ. Trong mô hình WRF, các quá trình tham số hóa vật lý—đặc biệt là các sơ đồ vi vật lý mô phỏng quá trình tạo mây, giáng thủy và bốc hơi—đóng vai trò quyết định đến độ chính xác của các bản tin dự báo. Các nghiên cứu thực nghiệm đã chỉ ra rằng ngay cả khi sử dụng cùng một điều kiện ban đầu, việc thay đổi sơ đồ vi vật lý có thể dẫn đến những sai khác đáng kể trong kết quả mô phỏng các hiện tượng thời tiết như mưa lớn, xoáy thuận nhiệt đới và trường nhiệt độ bề mặt.

Mặc dù giải pháp kết hợp nhiều sơ đồ vi vật lý để thiết lập hệ thống dự báo tổ hợp mang lại hiệu quả cao, việc triển khai hệ thống này đòi hỏi năng lực tính toán và nguồn kinh phí rất lớn, tạo nên rào cản kỹ thuật đối với các đài khí tượng thủy văn khu vực. Do đó, việc xác định một sơ đồ vi vật lý đơn lẻ tối ưu để ứng dụng trực tiếp vào quy trình dự báo nghiệp vụ là một nhu cầu cấp thiết. Xuất phát từ thực tế đó, đồ án tốt nghiệp kỹ sư ngành Khí tượng học với đề tài "Lựa chọn sơ đồ vi vật lý tối ưu trong mô hình WRF phục vụ dự báo trường nhiệt độ ở khu vực Nam Bộ trong giai đoạn bùng nổ gió mùa" đã được thực hiện nhằm giải quyết bài toán dự báo trường nhiệt độ tại khu vực Nam Bộ.

Nhiệm vụ nghiên cứu của đồ án được xác định cụ thể gồm 5 nội dung chính:

  1. Sơ lược về hoàn lưu chung khí quyển trong giai đoạn bùng nổ gió mùa.
  2. Tìm hiểu về mô hình WRF và 6 sơ đồ vi vật lý tham số hóa mây - giáng thủy.
  3. Khảo sát từng sơ đồ vi vật lý trong việc dự báo trường nhiệt độ ở khu vực Nam Bộ trong giai đoạn bùng nổ gió mùa.
  4. Phân tích kết quả trường nhiệt độ thu được qua 6 sơ đồ vi vật lý tại khu vực Nam Bộ.
  5. Đánh giá kết quả sai số dự báo của 6 sơ đồ vi vật lý so với số liệu quan trắc trong hạn dự báo 3 ngày.

Đối tượng nghiên cứu của đồ án là trường nhiệt độ bề mặt và cấu trúc hoàn lưu khí quyển quy mô lớn liên quan tại khu vực Nam Bộ. Phạm vi không gian nghiên cứu bao trùm khu vực Nam Bộ và lân cận thông qua miền tính lồng ghép. Phạm vi thời gian nghiên cứu tập trung vào giai đoạn bùng nổ gió mùa chuyển mùa từ cuối tháng 4 đến đầu tháng 5 năm 2009 (từ ngày 29/04/2009 đến ngày 02/05/2009).

Cơ sở lý thuyết và phương pháp

Đồ án xây dựng trên nền tảng lý thuyết động lực học khí quyển vùng nhiệt đới và hệ thống mô hình dự báo số trị WRF phiên bản ARW (Advanced Research WRF). Mô hình ARW là hệ phương trình khí quyển hoàn toàn chịu nén, phi thủy tĩnh, sử dụng hệ tọa độ ngang lưới so le Arakawa C và hệ tọa độ thẳng đứng theo địa hình (địa thế vị quy chuẩn). Việc tích phân số trị theo thời gian sử dụng sơ đồ sai phân Runge-Kutta bậc 3, kết hợp với các bước thời gian phân chia nhỏ cho sóng âm thanh và sóng trọng trường.

Quy trình chuẩn bị dữ liệu đầu vào và khởi tạo mô hình được vận hành thông qua hệ thống tiền xử lý WPS (WRF Preprocessing System) bao gồm 3 chương trình chính: geogrid.exe (xác định phạm vi miền tính, nội suy dữ liệu tĩnh về địa hình, loại đất và sử dụng đất theo chuẩn USGS/MODIS), ungrib.exe (giải nén dữ liệu khí tượng định dạng GRIB), và metgrib.exe (nội suy các trường khí tượng lên lưới tính toán). Hệ thống lõi wrf.exe nhận điều kiện ban đầu và điều kiện biên từ chương trình real.exe để tiến hành tích phân phương trình dự báo.

Nghiên cứu tiến hành khảo sát và đánh giá so sánh 6 sơ đồ vi vật lý được tích hợp trong mô hình WRF:

  • Sơ đồ Kessler (mp_physics = 1): Sơ đồ vi vật lý mây ẩm đơn giản, dự đoán 3 biến thủy vĩ (hơi nước, nước mây, nước mưa), không xét pha băng.
  • Sơ đồ Lin/Purdue Lin (mp_physics = 2): Sơ đồ vi vật lý tinh vi gồm 6 biến thủy vĩ (hơi nước, nước mây, nước mưa, mây băng, tuyết và graupel), tính toán đầy đủ các quá trình pha băng và pha hỗn hợp.
  • Sơ đồ WSM3 (mp_physics = 3): Sơ đồ 3 cấp WSM đơn giản có xét pha băng chẩn đoán theo nhiệt độ đóng băng, dự đoán 3 biến thủy vĩ.
  • Sơ đồ WSM5 (mp_physics = 4): Sơ đồ 5 cấp WSM dự đoán hơi nước, nước mây, nước mưa, mây băng và tuyết; cho phép tồn tại nước siêu lạnh.
  • Sơ đồ Ferrier/Eta (mp_physics = 5): Sơ đồ vận hành nghiệp vụ của NCEP, dự đoán biến tổng phần ngưng tụ kết hợp với chẩn đoán phân chia pha rắn và lỏng.
  • Sơ đồ WSM6 (mp_physics = 6): Sơ đồ 6 cấp WSM mở rộng từ WSM5, bổ sung biến graupel và tham số hóa vận tốc rơi riêng biệt cho hỗn hợp pha rắn.
Sơ đồ vi vật lý Ký hiệu thử nghiệm Số lượng biến Quá trình pha băng Quá trình pha hỗn hợp
Kessler TN1 3 Không (N) Không (N)
Purdue Lin TN2 6 Có (Y) Có (Y)
WSM3 TN3 3 Có (Y) Không (N)
WSM5 TN4 5 Có (Y) Không (N)
Eta (Ferrier) TN5 2 Có (Y) Có (Y)
WSM6 TN6 6 Có (Y) Có (Y)

Phương pháp đánh giá định lượng sử dụng các chỉ số thống kê sai số chuẩn bao gồm:

  • Sai số tuyệt đối trung bình (MAE - Mean Absolute Error): $$\text{MAE} = \frac{1}{n} \sum_{i=1}^{n} |F_i - O_i|$$
  • Sai số quân phương (RMSE - Root Mean Squared Error): $$\text{RMSE} = \sqrt{\frac{1}{n} \sum_{i=1}^{n} (F_i - O_i)^2}$$
  • Sai số căn quân phương năng lượng trung bình thể tích (EME - Volume-averaged Energy Root Mean Squared Errors): $$\text{EME} = \left[ \frac{1}{2} \left( U'^2 + V'^2 + \frac{C_p}{\bar{T}} T'^2 \right) \right]^{1/2}$$ Trong đó $F_i$ là giá trị dự báo, $O_i$ là giá trị quan trắc, $U', V', T'$ là độ lệch giữa dự báo và thực tế của các thành phần gió vĩ tuyến, gió kinh tuyến và nhiệt độ, $C_p$ là nhiệt dung đẳng áp, $\bar{T}$ là nhiệt độ trung bình.

Nguồn dữ liệu sử dụng bao gồm: dữ liệu tái phân tích khí quyển từ Trung tâm Quốc gia Dự báo Môi trường Hoa Kỳ (NCEP) đóng vai trò số liệu phân tích kiểm chứng (ANA), dữ liệu điều kiện biên toàn cầu GFS cập nhật định kỳ 6 giờ/lần, và số liệu thực đo nhiệt độ tại mạng lưới các trạm khí tượng bề mặt khu vực Nam Bộ.

Nội dung chính theo từng chương

Chương 1: Tổng quan về mô hình WRF và các sơ đồ vi vật lý

Chương 1 hệ thống hóa cấu trúc tổng thể và nguyên lý vận hành của hệ thống mô hình WRF (WRFV3) cùng hệ thống tiền xử lý WPS. Tác giả phân tích chi tiết cơ chế động lực học phi thủy tĩnh, phương pháp chia lưới Arakawa C và quy trình tích hợp các sơ đồ tham số hóa vật lý bức xạ, lớp biên hành tinh (PBL), đối lưu mây tích và vi vật lý.

Nội dung trọng tâm của chương là phân loại và làm rõ bản chất nhiệt động lực học của 6 sơ đồ vi vật lý tham số hóa mây và giáng thủy (Kessler, Lin, WSM3, WSM5, Eta, WSM6). Tác giả làm rõ sự khác biệt giữa các sơ đồ đơn giản không chứa pha rắn (Kessler) và các sơ đồ phức tạp mô phỏng chi tiết 6 dạng hạt thủy khí cùng pha hỗn hợp (Lin, WSM6), chỉ ra tầm ảnh hưởng của việc mô phỏng graupel và nước siêu lạnh đối với cân bằng bức xạ nhiệt và phân bố nhiệt độ lớp biên bề mặt.

Chương 2: Đặc điểm thời tiết giai đoạn bùng nổ gió mùa tại Nam Bộ và phương pháp đánh giá

Chương 2 phân tích quy luật khí hậu và hình thái synop trong giai đoạn chuyển mùa và bùng nổ gió mùa Tây Nam tại Nam Bộ. Hoàn lưu khí quyển Nam Bộ chịu sự khống chế sâu sắc của chu kỳ gió mùa nội chí tuyến: mùa khô thịnh hành gió đông bắc do tác động của áp cao lục địa và đới tín phong; mùa mưa thịnh hành gió tây nam do sự mở rộng của áp thấp nóng Ấn Miến và dòng khí vượt xích đạo từ Nam bán cầu.

Giai đoạn bùng nổ gió mùa (thường từ cuối tháng 4 đến đầu tháng 5) là thời điểm chuyển pha đột ngột với sự tranh chấp quyết liệt giữa lưỡi áp cao cận nhiệt Tây Thái Bình Dương và rãnh thấp kinh tuyến vịnh Bengal, đi kèm sự dịch chuyển lên phía bắc của dòng xiết cận chí tuyến trên mực 200hPa. Chương 2 cũng trình bày chi tiết cơ sở toán học của các chỉ tiêu thống kê đánh giá dự báo (MAE, RMSE) và công thức sai số năng lượng tổng hợp EME nhằm phục vụ việc kiểm định kết quả mô phỏng.

Chương 3: Kết quả thí nghiệm

Chương 3 trình bày toàn bộ kết quả thực nghiệm dự báo số trị cho đợt bùng nổ gió mùa từ 29/04/2009 đến 02/05/2009. Cấu hình thử nghiệm được thiết lập trên hệ tọa độ Mercator với tâm miền tính tại $12^\circ\text{N} - 106.5^\circ\text{E}$, cấu trúc 2 miền lưới lồng nhau:

  • Miền lưới 1 (D01): $103 \times 103$ điểm lưới, độ phân giải không gian $54\text{ km}$, bao phủ phạm vi $13^\circ\text{S} - 37^\circ\text{N}$ và $81.5^\circ\text{E}$.
  • Miền lưới 2 (D02): $91 \times 91$ điểm lưới, độ phân giải không gian $18\text{ km}$, bao trùm khu vực Nam Bộ và lân cận. Đỉnh mô hình đặt ở mực áp suất $10\text{ hPa}$.

Tác giả tiến hành phân tích chi tiết kết quả mô phỏng trên các trường khí tượng:

  • Hoàn lưu quy mô lớn mực 850hPa: Ngày 30/04/2009, số liệu phân tích ANA phản ánh lưỡi áp cao cận nhiệt bao trùm miền Bắc và vùng áp thấp yếu tại Biển Đông ($12^\circ\text{N} - 118^\circ\text{E}$). Thử nghiệm TN2 (Lin) và TN5 (Eta) mô phỏng chính xác cấu trúc áp cao cận nhiệt; TN2 và TN6 (WSM6) tái hiện đúng vị trí yên khí áp tại $4^\circ\text{N} - 106^\circ\text{E}$. Đến ngày 01/05/2009, TN2 tiếp tục mô phỏng sát nhất sự suy yếu của yên khí áp và dịch chuyển của tâm áp thấp Biển Đông.
  • Hoàn lưu mực 500hPa và 200hPa: Ở mực 500hPa, TN1 và TN2 tái hiện tốt rãnh gió tây phía bắc và dòng gió vượt xích đạo tây nam. Ở mực 200hPa, các sơ đồ TN1, TN2, TN5 mô phỏng tương đồng đới gió tây với vận tốc dòng xiết đạt 50–80 m/s quanh vĩ tuyến $30^\circ\text{N}$.
  • Cấu trúc thẳng đứng của gió qua trạm Tân Sơn Hòa: Mặt cắt gió vĩ tuyến ($U$) dọc kinh tuyến $106.7^\circ\text{E}$ và gió kinh tuyến ($V$) dọc vĩ tuyến $10.8^\circ\text{N}$ cho thấy TN5 và TN6 mô phỏng chuẩn xác đới gió tây đạt 8 m/s tại các mực 900hPa và 400hPa, trong khi TN1 có xu hướng khuếch đại cường độ gió so với quan trắc.
  • Sai số tốc độ gió và địa thế vị mực 500hPa: Tại khu vực Nam Bộ, TN3 và TN5 đạt sai số tốc độ gió nhỏ nhất ngày 30/04, trong khi TN2, TN4, TN5 duy trì sai số nhỏ nhất ($\pm 1\text{ m/s}$) vào ngày 01/05. Đối với sai số độ cao địa thế vị, TN2, TN4, TN5 đạt độ lệch thấp nhất trong khoảng $\pm 100\text{ m}^2/\text{s}^2$.

Kết quả và đóng góp

Qua việc so sánh định tính các bản đồ synop và định lượng các chỉ số sai số đối với trường gió, địa thế vị và trường nhiệt độ bề mặt 2m, đồ án rút ra các kết quả đánh giá thực nghiệm tổng hợp giữa 6 phương án thử nghiệm:

Sơ đồ thử nghiệm Mô phỏng hoàn lưu 850hPa Mô phỏng dòng xiết 200hPa Sai số vận tốc gió 500hPa Sai số địa thế vị 500hPa Hiệu quả dự báo trường nhiệt độ
TN1 (Kessler) Xuất hiện sai lệch lớn ở vùng biển Mô phỏng được đới gió tây Cao hơn thực tế Thấp hơn thực tế Sai số lớn
TN2 (Lin) Tương đồng cao với ANA Tương đồng cao với ANA Nhỏ ($\pm 1\text{ m/s}$) Nhỏ ($\pm 100\text{ m}^2/\text{s}^2$) Tối ưu cao
TN3 (WSM3) Vùng thấp Biển Đông bị yếu Tương đồng trung bình Trung bình Lệch âm so với thực tế Trung bình
TN4 (WSM5) Vùng thấp Biển Đông bị yếu Tương đồng trung bình Nhỏ vào hạn 24h Nhỏ ($\pm 100\text{ m}^2/\text{s}^2$) Trung bình
TN5 (Eta) Tương đồng cao với ANA Tương đồng cao với ANA Nhỏ ($\pm 1\text{ m/s}$) Nhỏ ($\pm 100\text{ m}^2/\text{s}^2$) Tối ưu cao
TN6 (WSM6) Tái hiện được yên khí áp xích đạo Tương đồng trung bình Trung bình Nhỏ ($\pm 100\text{ m}^2/\text{s}^2$) Tối ưu cao

Kết quả nghiên cứu kết luận rằng trong giai đoạn bùng nổ gió mùa năm 2009 tại Nam Bộ, các sơ đồ vi vật lý Lin (mp_physics = 2), Eta (mp_physics = 5) và WSM6 (mp_physics = 6) cho chất lượng dự báo trường nhiệt độ và các trường hoàn lưu khí quyển liên quan đạt độ chính xác cao nhất so với các sơ đồ còn lại.

Về mặt đóng góp khoa học và thực tiễn, đồ án cung cấp tài liệu kiểm nghiệm thực tế đầu tiên về việc đánh giá độ nhạy của các sơ đồ vi vật lý chuyên biệt cho bài toán dự báo trường nhiệt độ tại Nam Bộ trong thời kỳ bùng nổ gió mùa. Nghiên cứu đưa ra khuyến nghị thực tế cho các đài dự báo khí tượng thủy văn khu vực trong việc lựa chọn cấu hình vi vật lý tối ưu (Lin, Eta hoặc WSM6) khi vận hành mô hình dự báo đơn lẻ WRF nhằm tiết kiệm chi phí tính toán mà vẫn đảm bảo độ tin cậy của bản tin thời tiết.

Hạn chế và hướng nghiên cứu tiếp

Đồ án tồn tại một số hạn chế gắn liền với phạm vi và điều kiện nghiên cứu:

  • Dữ liệu thực nghiệm mới chỉ tập trung vào một đợt bùng nổ gió mùa cụ thể (đợt chuyển mùa cuối tháng 4 - đầu tháng 5 năm 2009), chưa thực hiện khảo nghiệm trên chuỗi số liệu nhiều năm để tăng tính khái quát hóa thống kê.
  • Nghiên cứu tập trung đánh giá độ nhạy của tham số hóa vi vật lý mà chưa khảo sát tương tác đồng thời giữa vi vật lý với các sơ đồ tham số hóa lớp biên hành tinh (PBL) và tham số hóa đối lưu mây tích khác nhau.

Hướng nghiên cứu tiếp theo được định hướng:

  • Mở rộng phạm vi thử nghiệm cho nhiều đợt bùng nổ gió mùa trong các năm có điều kiện khí hậu đối nghịch (El Niño và La Niña).
  • Nghiên cứu phối hợp đa tham số (kết hợp sơ đồ vi vật lý Lin/WSM6 với các sơ đồ đối lưu Kain-Fritsch, Betts-Miller-Janjic và các sơ đồ PBL khác nhau) để xây dựng bộ tổ hợp vật lý tối ưu cho khu vực Nam Bộ.

Giá trị tham khảo

Đồ án tốt nghiệp là tài liệu tham khảo hữu ích đối với:

  • Sinh viên, học viên cao học và cán bộ nghiên cứu chuyên ngành Khí tượng học, Thủy văn học và Khoa học Môi trường trong việc tiếp cận phương pháp cài đặt, thiết lập miền tính và cấu hình các gói vật lý trong mô hình WRF.
  • Các dự báo viên tại Đài Khí tượng Thủy văn khu vực Nam Bộ và các trung tâm dự báo thời tiết trong việc lựa chọn sơ đồ tham số hóa vi vật lý tối ưu khi lập quy trình công nghệ dự báo trường nhiệt độ và hoàn lưu gió mùa.
  • Các nhà nghiên cứu mô hình hóa khí quyển tham khảo cấu trúc thiết lập miền tính lồng Mercator ($54\text{ km} / 18\text{ km}$) và phương pháp đánh giá sai số năng lượng thể tích EME kết hợp đối chiếu trạm thực tế.

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

1. Mục tiêu cốt lõi của đồ án là gì?
Đồ án nhằm khảo sát, so sánh và lựa chọn sơ đồ vi vật lý tối ưu nhất trong số 6 sơ đồ tích hợp trong mô hình WRF để phục vụ công tác dự báo trường nhiệt độ tại khu vực Nam Bộ trong giai đoạn bùng nổ gió mùa.

2. Sáu sơ đồ vi vật lý được đưa vào thử nghiệm gồm những sơ đồ nào?
Sáu sơ đồ được thử nghiệm bao gồm: Kessler (TN1), Purdue Lin (TN2), WSM3 (TN3), WSM5 (TN4), Ferrier/Eta (TN5), và WSM6 (TN6).

3. Miền tính và độ phân giải không gian của mô hình WRF trong đồ án được thiết lập ra sao?
Mô hình sử dụng phép chiếu Mercator với tâm miền tính tại $12^\circ\text{N} - 106.5^\circ\text{E}$, đỉnh mô hình $10\text{ hPa}$, gồm 2 miền lưới lồng nhau: miền ngoài (D01) có $103 \times 103$ điểm lưới với độ phân giải $54\text{ km}$; miền trong (D02) có $91 \times 91$ điểm lưới với độ phân giải $18\text{ km}$ bao trùm Nam Bộ. Điều kiện biên được cập nhật 6 giờ/lần từ dữ liệu GFS.

4. Tác giả sử dụng các phương pháp và chỉ số nào để đánh giá sai số dự báo?
Tác giả sử dụng các chỉ số thống kê sai số tuyệt đối trung bình (MAE), sai số quân phương (RMSE), sai số căn quân phương năng lượng trung bình thể tích (EME), kết hợp phương pháp đối chiếu hình thái synop với số liệu tái phân tích NCEP (ANA) và so sánh trực tiếp với số liệu đo đạc tại các trạm khí tượng bề mặt Nam Bộ trong hạn dự báo 3 ngày.

5. Kết quả thực nghiệm xác định những sơ đồ vi vật lý nào đạt hiệu quả dự báo tốt nhất?
Kết quả phân tích hoàn lưu và trường nhiệt độ khẳng định ba sơ đồ vi vật lý Lin (TN2), Eta (TN5) và WSM6 (TN6) cho kết quả mô phỏng trường nhiệt độ, trường gió và độ cao địa thế vị chính xác nhất, có sai số nhỏ nhất so với số liệu quan trắc thực tế.

Kết luận

Đồ án tốt nghiệp của tác giả Âu Thị Thanh Tâm đã thực hiện thành công việc đánh giá so sánh 6 sơ đồ vi vật lý trong mô hình WRF đối với trường nhiệt độ khu vực Nam Bộ trong đợt bùng nổ gió mùa năm 2009. Thông qua quy trình mô phỏng số trị bài bản và các phương pháp thống kê sai số chặt chẽ, nghiên cứu đã chỉ ra các sơ đồ Lin, Eta và WSM6 mang lại hiệu quả dự báo vượt trội về cả trường hoàn lưu quy mô lớn lẫn trường nhiệt độ bề mặt. Kết quả này cung cấp luận cứ kỹ thuật quan trọng giúp các cơ quan dự báo khí tượng thủy văn khu vực tối ưu hóa cấu hình mô hình số trị đơn lẻ trong điều kiện tài nguyên tính toán còn hạn chế.