Giới thiệu dự án
Dự báo thời tiết số trị (Numerical Weather Prediction - NWP) đóng vai trò xương sống trong công tác khí tượng thủy văn hiện đại. Tại khu vực Nam Bộ, giai đoạn chuyển mùa từ mùa khô sang mùa mưa – đặc biệt là thời kỳ bùng nổ gió mùa Tây Nam (Southwest Monsoon Onset) vào cuối tháng 4 đầu tháng 5 – thường đi kèm các hiện tượng thời tiết cực đoan như dông nhiệt, lốc xoáy, mưa lớn cục bộ và sự biến động nhiệt độ bề mặt gay gắt. Các hiện tượng này ảnh hưởng trực tiếp đến hơn 20 triệu dân cư, ngành nông nghiệp lúa nước và các hoạt động kinh tế biển.
Tuy nhiên, việc mô phỏng chính xác trường nhiệt độ ($T_{2m}$ và nhiệt độ các tầng đối lưu 850hPa, 500hPa, 200hPa) trong thời kỳ bùng nổ gió mùa gặp nhiều khó khăn do tính phi tuyến mạnh mẽ của hoàn lưu quy mô vừa (mesoscale) và quá trình chuyển pha nhiệt động lực học trong các khối mây đối lưu sâu. Mô hình WRF (Weather Research and Forecasting) phiên bản ARW (Advanced Research WRF) cung cấp nhiều sơ đồ tham số hóa vi vật lý (Microphysics Parameterization Schemes), nhưng mỗi sơ đồ lại phản ánh các cơ chế ngưng kết, đóng băng và giải phóng ẩn nhiệt khác nhau.
Đồ án tập trung giải quyết bài toán: Xác định sơ đồ vi vật lý tối ưu nhất trong 6 sơ đồ điển hình của WRF nhằm nâng cao độ chính xác dự báo trường nhiệt độ hạn 72 giờ tại khu vực Nam Bộ.
Mục tiêu cụ thể của đồ án:
- Phân tích đặc điểm hoàn lưu khí quyển quy mô lớn trong giai đoạn bùng nổ gió mùa Tây Nam năm 2009 tại khu vực Nam Bộ và Biển Đông.
- Thiết lập cấu hình hệ thống mô hình WRF-ARW với 2 miền lưới lồng ghép (horizontal resolutions: 54 km và 18 km).
- Triển khai 6 kịch bản thử nghiệm vi vật lý độc lập (
mp_physics = 1, 2, 3, 4, 5, 6 tương ứng Kessler, Lin, WSM3, WSM5, Eta Ferrier, WSM6).
- Đánh giá kiểm chứng trường nhiệt độ dự báo và cấu trúc dòng thẳng đứng với dữ liệu phân tích tái phân tích toàn cầu NCEP Final Analysis (FNL/GFS) và chuỗi số liệu thực đo tại các trạm khí tượng thủy văn Nam Bộ.
- Đề xuất sơ đồ đơn lẻ tối ưu hóa tài nguyên tính toán cho các đài khí tượng thủy văn khu vực chưa đủ năng lực vận hành hệ thống dự báo tổ hợp (Ensemble Prediction System - EPS).
Phạm vi nghiên cứu giới hạn trong đợt chuyển mùa lịch sử từ ngày 29/04/2009 đến 02/05/2009 với hạn dự báo 72 giờ (3 ngày).
Phân tích và thiết kế giải pháp
Phân tích hiện trạng
Tại các trung tâm dự báo khu vực, việc chạy dự báo tổ hợp đa sơ đồ (multi-physics ensemble) đòi hỏi hệ thống siêu máy tính (HPC cluster) với hàng trăm core CPU và bộ nhớ khổng lồ. Trong điều kiện tài nguyên máy chủ hạn chế, việc chọn ra một cấu hình vi vật lý "chuẩn" (deterministic single-run) có độ tin cậy cao là ưu tiên hàng đầu.
| Sơ đồ vi vật lý |
Số lượng biến pha ẩm |
Quá trình pha băng |
Quá trình hỗn hợp pha |
Ưu điểm chính |
Hạn chế lớn |
| Kessler (TN1) |
3 |
Không (N) |
Không (N) |
Tốc độ tính toán siêu nhanh, tài nguyên thấp |
Bỏ qua hoàn toàn pha băng, sai số nhiệt độ tầng đối lưu lớn |
| Purdue Lin (TN2) |
6 |
Có (Y) |
Có (Y) |
Mô phỏng chi tiết 6 dạng hạt ẩm kể cả mưa đá (graupel) |
Tốn tài nguyên, tính toán bão hòa hơi nước phức tạp |
| WSM3 (TN3) |
3 |
Có (Y) |
Không (N) |
Hiệu quả tính toán cho mây băng đơn giản |
Thiếu nước siêu lạnh, chuyển pha nhiệt độ đột ngột |
| WSM5 (TN4) |
5 |
Có (Y) |
Không (N) |
Tách biệt băng mây và tuyết, cho phép nước siêu lạnh |
Chưa mô phỏng được graupel, sai số ẩn nhiệt đối lưu mạnh |
| Ferrier / Eta (TN5) |
2 |
Có (Y) |
Có (Y) |
Tối ưu hóa bộ nhớ cho mô hình nghiệp vụ NCEP |
Chẩn đoán phân đoạn pha băng gộp, hạn chế phân giải cao |
| WSM6 (TN6) |
6 |
Có (Y) |
Có (Y) |
Tích hợp graupel và tốc độ rơi hạt hỗn hợp mới |
Cần năng lực tính toán trung bình-cao |
Ma trận yêu cầu hệ thống theo mô hình MoSCoW:
- Must have: Mô phỏng chính xác trường dòng gió 850hPa, 500hPa, 200hPa và nhiệt độ bề mặt ($T_{2m}$); xuất file NetCDF chuẩn
wrfout.
- Should have: Tích hợp đầy đủ quá trình giải phóng ẩn nhiệt ngưng kết ($L_v$), ẩn nhiệt đông đặc ($L_f$) và ẩn nhiệt thăng hoa ($L_s$).
- Could have: Đánh giá sai số năng lượng tổng thể thông qua chỉ số $EME$.
- Won't have (trong phạm vi đề tài): Chạy đồng hóa dữ liệu vệ tinh Radar/Doppler thời gian thực (Data Assimilation).
Thiết kế hệ thống
Kiến trúc thử nghiệm dự báo số trị sử dụng mô hình WRF-ARW Version 3 bao gồm các module chính:
+-------------------------------------------------------------------------+
| WRF PRE-PROCESSING SYSTEM (WPS) |
| +--------------------+ +---------------------+ +----------------+ |
| | geogrid.exe (D01, |-->| ungrib.exe (GRIB1/ |-->| metgrib.exe | |
| | D02 Domain Setup) | | GRIB2 GFS Parsing) | | (Interpolate) | |
| +--------------------+ +---------------------+ +-------+--------+ |
+-------------------------------------------------------------|-----------+
v
+-------------------------------------------------------------------------+
| WRF MODEL INTEGRATION ENGINE |
| +-------------------------------------------------------------------+ |
| | real.exe (Boundary & Initial Condition Generation: wrfbdy/wrfinput)| |
| +-----------------------------------+-------------------------------+ |
| v |
| +-------------------------------------------------------------------+ |
| | wrf.exe (Non-hydrostatic, Time-split RK3 Integration, C-Grid) | |
| | * Dynamics: Fully compressible non-hydrostatic | |
| | * Microphysics: Options 1 to 6 (Kessler -> WSM6) | |
| | * Cumulus: Kain-Fritsch / Betts-Miller-Janjic | |
| | * Radiation: RRTM Longwave / Dudhia Shortwave | |
| | * PBL: Yonsei University (YSU) Scheme | |
| +-----------------------------------+-------------------------------+ |
+--------------------------------------|----------------------------------+
v
+-------------------------------------------------------------------------+
| POST-PROCESSING & STATISTICAL METRICS |
| +-------------------------------------------------------------------+ |
| | NCL / GrADS Scripts -> Extract Variables -> MAE, RMSE, EME Engine | |
| +-------------------------------------------------------------------+ |
+-------------------------------------------------------------------------+
Thông số cấu hình không gian và động lực học:
- Tâm miền tính: $12.0^\circ\text{N} - 106.5^\circ\text{E}$ (phép chiếu Mercator).
- Miền lưới 1 (D01 - Vùng rộng lớn): Kích thước $103 \times 103$ điểm lưới, độ phân giải ngang $54\text{ km}$, bao phủ từ $13^\circ\text{S} - 37^\circ\text{N}$ và $81.5^\circ\text{E} - 131.5^\circ\text{E}$.
- Miền lưới 2 (D02 - Khu vực Nam Bộ lồng ghép): Kích thước $91 \times 91$ điểm lưới, độ phân giải ngang $18\text{ km}$.
- Cấu trúc thẳng đứng: 28 tầng $\eta$ địa hình theo sau (terrain-following hydrostatic pressure coordinate), đỉnh mô hình tại $10\text{ hPa}$.
- Dữ liệu đầu vào: Toàn cầu NCEP GFS $0.5^\circ \times 0.5^\circ$, cập nhật biên 6 giờ/lần.
Methodology
Quy trình nghiên cứu áp dụng phương pháp thực nghiệm kiểm chứng số trị kết hợp phân tích thống kê khí tượng:
- Thiết lập chuẩn hóa môi trường mô phỏng: Cố định toàn bộ các tham số vật lý bức xạ, lớp biên khí quyển (PBL), đối lưu mây tích, mặt đệm đất; chỉ thay đổi duy nhất tham số vi vật lý
mp_physics.
- Kế hoạch thực hiện:
- Giai đoạn 1: Xử lý dữ liệu nền địa hình (USGS 30s) và dữ liệu khí tượng GFS cho đợt bùng nổ gió mùa 29/04/2009 - 02/05/2009.
- Giai đoạn 2: Khởi tạo 6 mô phỏng WRF tương ứng với TN1, TN2, TN3, TN4, TN5, TN6.
- Giai đoạn 3: Trích xuất trường nhiệt độ ($T_{2m}$, $T_{850}$, $T_{500}$), độ cao địa thế vị ($H_{500}$), thành phần gió kinh hướng ($V$) và vĩ hướng ($U$).
- Giai đoạn 4: Tính toán sai số định lượng so với số liệu tái phân tích ANA (NCEP) và trạm quan trắc bề mặt Tân Sơn Hòa, Cần Thơ, Cà Mau, Vũng Tàu.
Implementation và kết quả
Development process
Cấu hình chi tiết trong tệp điều khiển namelist.input của WRF cho 6 thử nghiệm:
&domains
time_step = 180,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 2,
e_we = 103, 91,
e_sn = 103, 91,
e_vert = 28, 28,
dx = 54000, 18000,
dy = 54000, 18000,
grid_id = 1, 2,
parent_id = 0, 1,
i_parent_start = 1, 35,
j_parent_start = 1, 35,
parent_grid_ratio = 1, 3,
parent_time_step_ratio = 1, 3,
/
&physics
mp_physics = 6, 6, ! Thay doi tu 1 -> 6 cho TN1 -> TN6
ra_lw_physics = 1, 1, ! RRTM scheme
ra_sw_physics = 1, 1, ! Dudhia scheme
radt = 30, 10,
sf_sfclay_physics = 1, 1, ! Monin-Obukhov
sf_surface_physics = 2, 2, ! Noah LSM
bl_pbl_physics = 1, 1, ! YSU scheme
bldt = 0, 0,
cu_physics = 1, 1, ! Kain-Fritsch scheme
cudt = 5, 5,
isfflx = 1,
ifsnow = 0,
icloud = 1,
/
Các công thức toán học đánh giá sai số được lập trình tự động hóa:
-
Sai số tuyệt đối trung bình (Mean Absolute Error - MAE):
$$\text{MAE} = \frac{1}{n} \sum_{i=1}^n |F_i - O_i|$$
-
Sai số căn quân phương (Root Mean Square Error - RMSE):
$$\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 (Volume-averaged Energy Root Mean Squared Error - EME):
$$\text{EME} = \left[ \frac{1}{2} \overline{\left( U'^2 + V'^2 + \frac{C_p}{\bar{T}} T'^2 \right)} \right]^{1/2}$$
Trong đó: $U', V', T'$ lần lượt là độ lệch sai số giữa giá trị dự báo và phân tích của thành phần gió vĩ hướng, kinh hướng và nhiệt độ; $C_p = 1004.5\text{ J}/(\text{kg}\cdot\text{K})$ là nhiệt dung đẳng áp của không khí khô; $\bar{T}$ là nhiệt độ trung bình toàn miền tính.
Testing và validation
Quá trình mô phỏng được tiến hành liên tục trong 72 giờ từ 00:00 UTC ngày 29/04/2009 đến 00:00 UTC ngày 02/05/2009.
1. Đánh giá hoàn lưu quy mô lớn và cấu trúc thẳng đứng
- Mực 850hPa: Ngày 30/04/2009, khu vực Nam Bộ nằm trong vùng tranh chấp giữa rìa áp cao cận nhiệt đới Tây Thái Bình Dương ở phía Đông và rãnh áp thấp Vịnh Bengal phát triển ở phía Tây. Sơ đồ Lin (TN2), Eta (TN5) và WSM6 (TN6) mô phỏng chính xác vị trí yên khí áp tại tọa độ $4.0^\circ\text{N} - 106.0^\circ\text{E}$, khớp hoàn toàn với số liệu phân tích ANA của NCEP. Ngược lại, Kessler (TN1) làm rãnh thấp suy yếu quá nhanh.
- Mực 500hPa: Xuất hiện rãnh gió Tây ở phía Bắc và dòng gió Tây Nam vượt xích đạo thổi thẳng vào Nam Bộ. TN2, TN5 và TN6 tái tạo xuất sắc trường dòng gió với sai số tốc độ gió tại Nam Bộ nhỏ nhất, dao động trong khoảng $\pm 1.0\text{ m/s}$.
- Mặt cắt thẳng đứng gió vĩ tuyến ($U$) dọc theo kinh tuyến $106.7^\circ\text{E}$ qua trạm Tân Sơn Hòa: TN5 và TN6 mô phỏng chính xác sự phân tầng động lực: đới gió Tây tầng thấp ($900 - 400\text{ hPa}$) đạt vận tốc $8\text{ m/s}$ tại $10^\circ\text{N}$, và dòng xiết gió Tây trên cao ($200\text{ hPa}$) đạt vận tốc cực đại $36\text{ m/s}$ ở phía Bắc $20^\circ\text{N}$.
Độ cao (hPa)
0°N (Xích đạo) 10°N (Nam Bộ) 20°N
2. Sai số độ cao địa thế vị mực 500hPa ($H_{500}$)
- TN2, TN5 và TN6 cho sai số địa thế vị nhỏ nhất tại khu vực Nam Bộ, nằm trong biên độ tối ưu $\pm 100\text{ m}^2/\text{s}^2$ sau 48 giờ dự báo.
- TN1 cho giá trị địa thế vị thấp hơn thực tế (âm sâu), trong khi TN3 và TN4 có xu hướng đánh giá cao cường độ áp cao cận nhiệt đới.
Kết quả đạt được
Bảng tổng hợp sai số nhiệt độ tuyệt đối trung bình ($\text{MAE}{T{2m}}$, đơn vị: $^\circ\text{C}$) tại các trạm quan trắc Nam Bộ trong 3 ngày dự báo:
| Thử nghiệm |
Sơ đồ vi vật lý |
MAE Ngày 1 (29/04/2009) |
MAE Ngày 2 (30/04/2009) |
MAE Ngày 3 (01/05/2009) |
MAE Trung bình 72h |
| TN1 |
Kessler |
$1.82^\circ\text{C}$ |
$2.15^\circ\text{C}$ |
$2.48^\circ\text{C}$ |
$2.15^\circ\text{C}$ |
| TN2 |
Purdue Lin |
$1.12^\circ\text{C}$ |
$1.28^\circ\text{C}$ |
$1.35^\circ\text{C}$ |
$1.25^\circ\text{C}$ |
| TN3 |
WSM3 |
$1.65^\circ\text{C}$ |
$1.89^\circ\text{C}$ |
$2.10^\circ\text{C}$ |
$1.88^\circ\text{C}$ |
| TN4 |
WSM5 |
$1.42^\circ\text{C}$ |
$1.56^\circ\text{C}$ |
$1.72^\circ\text{C}$ |
$1.57^\circ\text{C}$ |
| TN5 |
Ferrier Eta |
$1.15^\circ\text{C}$ |
$1.30^\circ\text{C}$ |
$1.39^\circ\text{C}$ |
$1.28^\circ\text{C}$ |
| TN6 |
WSM6 |
$1.08^\circ\text{C}$ |
$1.22^\circ\text{C}$ |
$1.31^\circ\text{C}$ |
$1.20^\circ\text{C}$ |
Biểu đồ so sánh sai số MAE nhiệt độ trung bình 72h (°C)
2.15°C 1.25°C 1.88°C 1.57°C 1.28°C 1.28°C 1.20°C
Kết luận thực nghiệm:
- Sơ đồ WSM6 (TN6) đạt độ chính xác cao nhất với sai số $\text{MAE} = 1.20^\circ\text{C}$, cải thiện $44.2%$ độ chính xác so với sơ đồ cơ sở Kessler (TN1).
- Sơ đồ Purdue Lin (TN2) và Ferrier Eta (TN5) bám sát phía sau với sai số lần lượt là $1.25^\circ\text{C}$ và $1.28^\circ\text{C}$.
- Sự vượt trội của WSM6 và Lin đến từ việc tính toán đầy đủ thành phần graupel và tương tác hỗn hợp pha (mixed-phase), giúp tái lập chính xác cơ chế giải phóng ẩn nhiệt khi bùng nổ gió mùa Tây Nam.
Đổi mới và đóng góp
- Lấp đầy khoảng trống nghiên cứu cục bộ: Đây là công trình đầu tiên thực hiện đánh giá đối chứng chi tiết và có hệ thống 6 sơ đồ vi vật lý độc lập của WRF chuyên biệt cho bài toán trường nhiệt độ Nam Bộ trong thời kỳ bùng nổ gió mùa.
- Định lượng hóa hiệu năng vật lý pha ẩm: Chứng minh bằng thực nghiệm số trị rằng việc bỏ qua pha băng (như trong Kessler) dẫn đến sai số nhiệt độ bề mặt vượt ngưỡng $2.4^\circ\text{C}$ ở hạn dự báo 72h; trong khi việc đưa hạt graupel vào (WSM6) giúp hạ sai số xuống mức chấp nhận được trong nghiệp vụ ($\sim 1.2^\circ\text{C}$).
- Bộ tiêu chí đánh giá đa chiều: Không chỉ dừng lại ở các sai số thống kê truyền thống tại trạm bề mặt (MAE, RMSE), nghiên cứu áp dụng công thức sai số năng lượng tổng thể EME và mặt cắt cấu trúc thẳng đứng qua kinh vĩ độ trạm cao không Tân Sơn Hòa.
- Khuyến nghị cấu hình tối ưu chi phí - hiệu năng: Xác lập sơ đồ WSM6 và Eta là cấu hình chuẩn cho các đài dự báo vừa và nhỏ, giúp đạt độ chính xác tương đương các hệ thống dự báo lớn mà không cần đầu tư cụm máy chủ phân tán đắt đỏ.
Ứng dụng thực tế và triển khai
Tình huống ứng dụng thực tế
- Cảnh báo nắng nóng gay gắt và dông sét chuyển mùa: Cung cấp thông tin dự báo trường nhiệt độ $T_{2m}$ chi tiết đến cấp huyện/thị xã tại Nam Bộ trước các đợt mưa rào đầu mùa.
- Phục vụ nông nghiệp Đồng bằng sông Cửu Long: Cảnh báo kịp thời nhiệt độ cao kết hợp bốc thoát hơi nước mạnh cuối mùa khô, hỗ trợ kế hoạch xuống giống lúa Hè Thu.
- Vận hành hệ thống điện và năng lượng tái tạo: Dự báo phụ tải điện sinh hoạt tại TP. Hồ Chí Minh dựa trên biên độ dao động nhiệt ngày và đêm.
Yêu cầu triển khai hệ thống (System Requirements)
- Hệ điều hành: Linux CentOS 7+ hoặc Ubuntu Server 20.04/22.04 LTS (x86_64).
- Trình biên dịch & Thư viện: GNU Fortran / C Compiler (
gcc, gfortran v9.x trở lên), OpenMPI v4.x, NetCDF-Fortran v4.5+, HDF5, JasPer, LibPNG, zlib.
- Phần cứng đề xuất:
- CPU: Tối thiểu 16 Core (Khuyến nghị 32 - 64 Core Intel Xeon hoặc AMD EPYC).
- RAM: 64 GB ECC DDR4/DDR5.
- Dung lượng lưu trữ: 2 TB SSD NVMe (ghi đọc dữ liệu I/O
wrfout tốc độ cao).
- Quy trình tự động hóa (Automation Script Pipeline):
#!/bin/bash
# Pipeline tu dong hoa chay mo hinh WRF voi so do WSM6
set -e
DATE=$(date -u +%Y%m%d)
DATADIR="/data/nwp/gfs/${DATE}00"
WRF_RUN="/opt/models/wrf/run"
echo "=== [1] Tai du lieu GFS 0.5 do tu NCEP ==="
# Download GFS GRIB2 files via curl/wget...
echo "=== [2] Chay WPS Preprocessing ==="
cd /opt/models/wps
./ungrib.exe
./metgrib.exe
echo "=== [3] Chay WRF ARW Solver (WSM6 Scheme) ==="
cd ${WRF_RUN}
mpirun -np 32 ./real.exe
mpirun -np 32 ./wrf.exe
echo "=== [4] Hau xu ly va trich xuat truong nhiet do ==="
python3 /opt/scripts/extract_temperature_metrics.py --input wrfout_d02* --output /var/www/html/forecast/
echo "Dự báo hoàn tất thành công!"
Hạn chế và hướng phát triển
Hạn chế kỹ thuật
- Chuỗi số liệu thử nghiệm còn hẹp: Nghiên cứu tập trung chuyên sâu vào đợt bùng nổ gió mùa năm 2009; cần mở rộng kiểm chứng trên chuỗi 10 - 20 năm để đạt tính đại diện khí hậu dài hạn.
- Chưa tích hợp đồng hóa dữ liệu địa phương (Data Assimilation): Trường ban đầu phụ thuộc hoàn toàn vào phân tích toàn cầu GFS, chưa khai thác dữ liệu quan trắc bề mặt tự động (AWS) và Radar Nhà Bè.
- Độ phân giải 18 km còn tương đối thô: Chưa phản ánh đầy đủ hiệu ứng đảo nhiệt đô thị (Urban Heat Island - UHI) tại các siêu đô thị như TP. Hồ Chí Minh.
Hướng phát triển tiếp theo
- Tăng độ phân giải miền lồng: Thiết lập miền lưới D03 với độ phân giải siêu cao $3\text{ km} - 1\text{ km}$ dạng mây phân giải trực tiếp (cloud-resolving scale), tắt sơ đồ đối lưu mây tích.
- Ứng dụng WRF-DA (Data Assimilation): Đồng hóa số liệu phản hồi vô tuyến Radar và vệ tinh địa tĩnh Himawari-8/9 bằng thuật toán 3DVAR/4DVAR.
- Kết hợp mô hình đô thị (WRF-Urban / BEP): Mô phỏng chi tiết tương tác bức xạ và bề mặt bê tông hóa tại khu vực nội thành Nam Bộ.
Đối tượng hưởng lợi
- Sinh viên & Kỹ sư ngành Khí tượng, Khoa học Môi trường: Nắm bắt quy trình khoa học chuẩn từ tiền xử lý WPS, cấu hình tham số vật lý mây đến kỹ thuật hậu xử lý sai số thống kê.
- Dự báo viên tại các Đài Khí tượng Thủy văn Khu vực: Có được cơ sở khoa học để thiết lập mô hình chạy nghiệp vụ hàng ngày với chi phí phần cứng tối thiểu mà vẫn đảm bảo độ chính xác nhiệt độ.
- Ngành Nông nghiệp & Thủy sản Nam Bộ: Nhận bản tin cảnh báo nhiệt độ chuyển mùa chuẩn xác hơn, giảm thiểu hiện tượng sốc nhiệt trên cây trồng và thủy sản nuôi trồng.
Câu hỏi thường gặp
1. Tại sao sơ đồ Kessler lại có sai số dự báo nhiệt độ lớn nhất?
Sơ đồ Kessler là sơ đồ vi vật lý thuần mây ấm đơn giản (warm-rain scheme), chỉ xem xét 3 biến pha ẩm ($q_v, q_c, q_r$) và hoàn toàn loại bỏ pha băng. Trong khí quyển nhiệt đới Nam Bộ vào thời kỳ chuyển mùa, mây đối lưu phát triển rất mạnh lên đến độ cao $12 - 16\text{ km}$ (vượt xa tầng đông kết $0^\circ\text{C}$). Việc không tính toán quá trình đóng băng và giải phóng ẩn nhiệt đông đặc dẫn đến ước lượng sai hoàn toàn cấu trúc nhiệt động lực học và nhiệt độ tầng mặt.
2. Sự khác biệt cốt lõi giữa WSM5 và WSM6 trong mô hình WRF là gì?
WSM5 (5-class) chỉ tính toán 5 dạng hạt: hơi nước, nước mây, nước mưa, băng mây và tuyết. Sơ đồ WSM6 (6-class) bổ sung thêm dạng hạt ẩm thứ 6 là mưa đá/graupel. Trong các đám mây đối lưu nhiệt đới, graupel đóng vai trò quyết định trong việc hình thành dòng giáng lạnh (cold pool downdraft) và làm thay đổi mạnh mẽ trường nhiệt độ mặt đất thông qua quá trình tan chảy và bốc hơi dưới chân mây.
3. Tại sao sơ đồ Ferrier (New Eta) lại đạt kết quả tốt dù số biến trạng thái ít?
Sơ đồ Ferrier Eta được thiết kế đặc thù cho các mô hình nghiệp vụ của NCEP (như mô hình Eta, NAM). Nó sử dụng kỹ thuật chẩn đoán hỗn hợp pha thông minh (diagnostic mixed-phase fraction) từ một biến ngưng kết tổng thể, giúp tiết kiệm bộ nhớ RAM và xung nhịp CPU nhưng vẫn phản ánh chính xác phân bố ẩn nhiệt trong quá trình thăng hoa và đông kết.
4. Cần cấu hình tài nguyên phần cứng thế nào để chạy ổn định mô hình này?
Để dự báo nghiệp vụ hạn 72 giờ cho khu vực Nam Bộ với 2 miền lưới ($54\text{ km} - 18\text{ km}$), cấu hình tối thiểu là máy trạm 16 Core CPU (Intel Xeon hoặc AMD EPYC), 32 GB RAM, hệ điều hành Linux 64-bit. Thời gian tích phân mô hình mất khoảng 25 - 40 phút cho mỗi phiên chạy.
5. Chỉ số sai số năng lượng $EME$ có ý nghĩa gì hơn so với RMSE thông thường?
RMSE chỉ đánh giá sai số cục bộ trên từng biến số riêng lẻ ($T$ hoặc $U$ hoặc $V$) mà không phản ánh được sự tương tác động lực học giữa trường nhiệt và trường gió. Chỉ số $EME$ tích hợp sai số của cả thành phần động năng động lực ($U'^2 + V'^2$) và thế năng nhiệt động ($\frac{C_p}{\bar{T}}T'^2$) trên toàn bộ thể tích khí quyển 3 chiều, giúp đánh giá độ ổn định tổng thể của toàn hệ thống mô phỏng.
Kết luận
Đồ án tốt nghiệp đã giải quyết thành công bài toán tối ưu hóa cấu hình vi vật lý trong mô hình số trị 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 Tây Nam năm 2009. Thông qua việc phân tích đối chứng 6 sơ đồ vi vật lý độc lập (Kessler, Lin, WSM3, WSM5, Eta, WSM6) với số liệu phân tích chuẩn NCEP và số liệu thực đo mặt đất:
- WSM6 và Purdue Lin được xác định là hai sơ đồ vi vật lý tối ưu nhất, mang lại trường dự báo nhiệt độ bề mặt và nhiệt độ các tầng cao với sai số thấp nhất ($\text{MAE} \approx 1.20^\circ\text{C}$ sau 72h), phản ánh chân thực hoàn lưu gió mùa và cấu trúc dòng thẳng đứng.
- Sơ đồ Ferrier Eta là lựa chọn thay thế xuất sắc khi bị giới hạn về tài nguyên tính toán.
- Kết quả nghiên cứu cung cấp cơ sở khoa học và kỹ thuật trực tiếp để chuẩn hóa quy trình dự báo thời tiết hạn ngắn tại các Đài Khí tượng Thủy văn khu vực phía Nam.
