Tổng quan luận án

Tính cấp thiết và khoảng trống nghiên cứu

Dự báo bão và xoáy thuận nhiệt đới bằng mô hình số trị đóng vai trò quan trọng trong nghiệp vụ khí tượng thủy văn tại Việt Nam cũng như trên thế giới. Trong những thập niên gần đây, kĩ năng dự báo quỹ đạo bão từ các hệ thống mô hình động lực toàn cầu và khu vực đã có sự cải thiện rõ rệt theo thời gian. Ngược lại, chất lượng dự báo cường độ bão (thể hiện qua tốc độ gió cực đại sát bề mặt $V_{max}$ và áp suất cực tiểu tại tâm $P_{min}$) vẫn còn nhiều hạn chế và tốc độ cải thiện diễn ra rất chậm. Trong thực tế nghiệp vụ, đối với các hạn dự báo ngắn (khoảng 24 giờ), dự báo cường độ bão chủ yếu dựa vào tính chất quán tính thay vì khai thác trực tiếp từ sản phẩm mô hình số; ở các hạn dự báo dài hơn (2–3 ngày), độ tin cậy của mô hình thường chỉ dừng lại ở mức nhận định xu thế tăng hoặc giảm cấp của cơn bão.

Khoảng trống nghiên cứu then chốt đặt ra là sự thiếu hụt các khảo sát định lượng về mối tương quan giữa kĩ năng mô phỏng quỹ đạo và kĩ năng mô phỏng cường độ bão trong một hệ thống động lực học thực tế. Các nghiên cứu thống kê trước đây chủ yếu sử dụng giả thiết quỹ đạo hoàn hảo dựa trên số liệu phân tích lại (quỹ đạo chuẩn Best-track) làm đầu vào cho mô hình thống kê, chưa cô lập và làm rõ được mức độ ảnh hưởng của việc giảm sai số quỹ đạo đến sự thay đổi sai số cường độ trong mô hình số trị khu vực. Đồng thời, bài toán nâng cao chất lượng trường phân tích ban đầu thông qua việc kết hợp đồng thời cấu trúc xoáy quy mô bão và hoàn lưu môi trường quy mô lớn bằng phương pháp đồng hóa số liệu hiện đại vẫn là một thách thức cần được giải quyết.

Mục tiêu nghiên cứu

Luận án xác định hai mục tiêu nghiên cứu cụ thể:

  1. Đánh giá định lượng mối quan hệ giữa sai số mô phỏng quỹ đạo và sai số mô phỏng cường độ bão trên khu vực Tây Bắc Thái Bình Dương dựa trên kết quả của hệ thống tổ hợp đa vật lý sử dụng mô hình số trị khu vực WRF-ARW.
  2. Thử nghiệm tăng cường thông tin cấu trúc xoáy ban đầu kết hợp thông tin môi trường quy mô lớn cho mô hình WRF-ARW bằng phương pháp đồng hóa tổ hợp lọc Kalman chuyển dạng địa phương (LETKF), từ đó đánh giá hiệu quả nâng cao chất lượng dự báo quỹ đạo và cường độ bão.

Đối tượng và phạm vi nghiên cứu

  • Đối tượng nghiên cứu: Các đặc trưng chuyển động và cấu trúc của bão/xoáy thuận nhiệt đới bao gồm vị trí tâm (quỹ đạo), tốc độ gió cực đại sát bề mặt $V_{max}$, khí áp cực tiểu tại tâm $P_{min}$; dự báo trung bình tổ hợp và dự báo xác suất tổ hợp; các nguồn phát sinh sai số và nhân tố chi phối chất lượng dự báo bão trong mô hình số trị.
  • Phạm vi không gian: Khu vực Tây Bắc Thái Bình Dương (TBTBD) và khu vực Biển Đông Việt Nam.
  • Phạm vi thời gian và mẫu nghiên cứu:
    • Tập số liệu đánh giá sai số bão trên Biển Đông giai đoạn 2008–2014 (62 cơn bão đi qua kinh tuyến $120^\circ\text{E}$).
    • Tập thử nghiệm mô phỏng tổ hợp đa vật lý gồm 92 trường hợp bão trên khu vực Tây Bắc Thái Bình Dương giai đoạn 2007–2010.
    • Tập thử nghiệm đồng hóa tổ hợp đối với các cơn bão mạnh điển hình giai đoạn 2013–2014: Usagi (2013), Nari (2013), Krosa (2013), Rammasun (2014), Neoguri (2014) và Vongfong (2014).

Tổng quan tài liệu và vị trí của luận án

Các hướng nghiên cứu trước đây

Luận án đã tổng hợp các hướng tiếp cận chính trong lịch sử dự báo bão:

  • Phương pháp Synop và Synop động lực: Dựa trên phân tích tương tác giữa các trung tâm hoạt động khí quyển quy mô lớn (Roger và Carleton 2001; Wang và cộng sự 1998) như dòng dẫn đường của cao áp cận nhiệt đới Tây Bắc Thái Bình Dương, tương tác không khí lạnh, hiệu ứng xoáy-beta (beta-gyre) và độ đứt gió môi trường. Phương pháp này phát huy hiệu quả tốt ở hạn ngắn 12–24 giờ nhưng bị giới hạn bởi tính phi tuyến của khí quyển.
  • Phương pháp thống kê và thống kê khí hậu:
    • Phương pháp quán tính kết hợp khí hậu (CLIPER, CLIPER5 của Knaff và cộng sự 2003) đóng vai trò là thước đo chuẩn (baseline) để đánh giá kĩ năng của mọi mô hình dự báo.
    • Các mô hình thống kê - động lực như STIPS (Knaff và cộng sự 2005) áp dụng tại Trung tâm Cảnh báo bão Hải quân Mỹ (JTWC), SHIPS (DeMaria và cộng sự 2005) và mô hình hồi quy Logistic LGEM (DeMaria 2009, 2010) sử dụng giả thiết hoàn hảo (Prognostic Perfect - Wilks 2006) kết hợp các nhân tố khí quyển từ mô hình động lực toàn cầu (NOGAPS, GFS).
  • Phương pháp mô hình số trị động lực: Tích phân hệ phương trình chuyển động thủy nhiệt động khí quyển đầy đủ (Navier-Stokes). Mô hình số đem lại chất lượng dự báo quỹ đạo vượt bậc ở các hạn dự báo từ 3 đến 10 ngày (Kalnay 2003) và có khả năng mô phỏng các trường hợp quỹ đạo dị thường hoặc tương tác nhiều bão (như bão Tembin và Bolaven năm 2012).
  • Nghiên cứu về sai số và kĩ năng dự báo bão:
    • Elsberry (1995) chỉ ra sai số xác định vị trí tâm bão thời gian thực từ ảnh vệ tinh dao động từ 20 đến 40 km.
    • Hui và cộng sự (2012), Fiorino và cộng sự (1993), Goerss và cộng sự (2004) đánh giá sai số quỹ đạo của Trung tâm RSMC Tokyo và các mô hình NOGAPS, UMKO, TYM tại Tây Bắc Thái Bình Dương.
    • DeMaria và cộng sự (2007, 2014) khảo sát xu thế sai số dự báo tại ba vùng biển Đại Tây Dương (ĐTD), Đông Thái Bình Dương (ĐTBD) và Tây Thái Bình Dương (TTBD), chỉ ra tốc độ giảm sai số trung bình tuyệt đối (MAE) của $V_{max}$ chỉ đạt 0,3–0,5%/năm (khoảng $-0,1\text{ kt/năm}$), trong khi sai số quỹ đạo giảm từ 3–5%/năm.
    • Chương trình Cải thiện chất lượng dự báo bão HFIP của NOAA (Gall và cộng sự 2013) đặt mục tiêu giảm 50% sai số quỹ đạo và cường độ bão trong 10 năm thông qua mô hình phân giải cao và đồng hóa tổ hợp.
  • Nghiên cứu về các nguồn sai số và nhân tố ảnh hưởng:
    • Edward Lorenz (1960s, 1970s) xác định ba nguồn bất định: sai số trường ban đầu, sai số kĩ năng mô hình (sơ đồ sai phân, tham số hóa quy mô dưới lưới) và tính bất định nội tại của khí quyển.
    • Ba nhóm nhân tố ảnh hưởng trực tiếp đến sai số cường độ bão (Knaff và cs 2003; Bender và cs 2007; DeMaria và cs 2009; Gall và cs 2013): (1) Cấu trúc xoáy ban đầu thô và thiếu thám sát (Bender và cs 1993; Kurihara và cs 1993; Davidson và cs 2000; Nguyen và Chen 2011); (2) Tham số hóa vật lý quy mô dưới lưới chưa phù hợp (Liu và cs 1997; Pu và Braun 2001; Pattnaik và Krishnamurti 2007; Osuri và cs 2011; Cangialosi và cs 2006; Kieu và cs 2013; Meng và Zhang 2007); (3) Sai số môi trường xung quanh bão bắt nguồn từ sai số quỹ đạo (Chan 1984, 2005).
Vùng biển và hạn dự báo Hệ số dốc dự báo cường độ (%/năm) Hệ số dốc dự báo quỹ đạo (%/năm)
Đại Tây Dương - Hạn 24h -1,0 -2,7
Đại Tây Dương - Hạn 48h -0,9 -3,4
Đại Tây Dương - Hạn 72h -0,6 -4,0
Đông Thái Bình Dương - Hạn 24h -0,1 -2,1
Đông Thái Bình Dương - Hạn 48h -0,4 -2,5
Đông Thái Bình Dương - Hạn 72h -0,8 -2,8
Tây Thái Bình Dương - Hạn 24h -0,2 -2,8
Tây Thái Bình Dương - Hạn 48h -0,3 -3,6
Tây Thái Bình Dương - Hạn 72h -0,5 -3,9

Bảng 1: Hệ số dốc đường xu thế sai số dự báo cường độ và quỹ đạo bão giai đoạn 1985–2005 (Nguồn: Trích xuất từ DeMaria và cộng sự 2007).

Vùng biển Hạn 24h Hạn 48h Hạn 72h Hạn 96h Hạn 120h
Đại Tây Dương 0,1 / 0,0 / 2,7 1,3 / 0,6 / 3,0 1,1 / 0,9 / 3,1 4,3 / 2,2 / 4,1 4,9 / 2,5 / 3,8
Đông Thái Bình Dương 0,8 / 0,5 / 2,2 1,2 / 1,1 / 2,3 1,5 / 1,5 / 2,4 0,8 / 2,3 / 3,6 1,3 / 2,3 / 4,2
Tây Thái Bình Dương 1,9 / -0,6 / 2,5 1,9 / 0,0 / 2,7 2,1 / 0,7 / 2,8 2,5 / 2,7 / 1,9 2,1 / 1,5 / 2,1

Bảng 2: Tỉ lệ cải thiện chất lượng trung bình năm giai đoạn 1989–2012 cho hạn 1–3 ngày và 2001–2012 cho hạn 4–5 ngày (Định dạng: Mô hình tốt nhất / Dự báo chính thức / Quỹ đạo chính thức; Đơn vị: %/năm; Nguồn: Trích xuất từ DeMaria và cộng sự 2014).

Khoảng trống luận án lựa chọn giải quyết

Luận án tập trung giải quyết hai khoảng trống chính:

  1. Khảo sát định lượng và cô lập mối tương quan giữa sai số dự báo quỹ đạo và cường độ bão trong mô hình số trị động lực khu vực WRF-ARW, kiểm chứng xem mức độ giảm sai số quỹ đạo thực tế đóng góp bao nhiêu phần trăm vào việc cải thiện sai số cường độ ở từng hạn dự báo cụ thể mà không sử dụng giả thiết quỹ đạo hoàn hảo.
  2. Xây dựng giải pháp kỹ thuật tích hợp đồng thời thông tin cấu trúc bão thời gian thực (TCVital) và thông tin gió vệ tinh quy mô lớn (AMV) thông qua thuật toán đồng hóa tổ hợp lọc Kalman chuyển dạng địa phương (LETKF), giải quyết hiện tượng thích ứng chậm và sai lệch cấu trúc ban đầu trong mô hình số trị.

Cơ sở lý thuyết và phương pháp nghiên cứu

Cơ sở lý thuyết và khung phân tích

  • Lý thuyết xác suất Bayes trong đồng hóa số liệu: Đồng hóa số liệu giải bài toán tìm trạng thái phân tích tối ưu $x_a$ từ trường nền ban đầu $x_b$ và véctơ quan trắc $y_o$. Dưới giả thiết sai số tuân theo phân bố chuẩn Gauss và không tương quan, hàm mục tiêu biến phân có dạng: $$J(x) = \frac{1}{2} (x - x_b)^T \mathbf{B}^{-1} (x - x_b) + \frac{1}{2} (y_o - \mathcal{H}(x))^T \mathbf{R}^{-1} (y_o - \mathcal{H}(x))$$ Trạng thái phân tích tối ưu thỏa mãn điều kiện triệt tiêu gradient $\nabla_x J(x_a) = 0$, dẫn đến công thức phân tích: $$x_a = x_b + \mathbf{W}(y_o - \mathcal{H}(x_b))$$ với ma trận trọng số $\mathbf{W} = \mathbf{B}\mathbf{H}^T (\mathbf{R} + \mathbf{H}\mathbf{B}\mathbf{H}^T)^{-1}$, trong đó $\mathbf{B}$ là ma trận hiệp phương sai sai số trường nền, $\mathbf{R}$ là ma trận hiệp phương sai sai số quan trắc và $\mathbf{H}$ là toán tử quan trắc tuyến tính hóa.
  • Phương pháp Lọc Kalman chuyển dạng địa phương (LETKF): Sử dụng một tập tổ hợp các dự báo để ước lượng ma trận hiệp phương sai sai số nền phụ thuộc vào dòng chảy thực tế. Thuật toán LETKF thực hiện quá trình phân tích độc lập tại từng thể tích địa phương trong không gian không gian trạng thái tổ hợp, giúp giảm chi phí tính toán và hạn chế các tương quan ảo ở khoảng cách xa.
  • Lý thuyết cấu trúc xoáy bão giải tích ba chiều: Sử dụng nghiệm giải tích cân bằng gradient và thủy tĩnh của Kieu và Zhang (2009) để thiết lập cấu trúc xoáy nhân tạo đối xứng ba chiều phù hợp với động lực học khí quyển từ các tham số bão thực tế.

Phương pháp nghiên cứu thực tế

  • Mô hình động lực khu vực WRF-ARW (phiên bản 3.x): Thiết lập hệ lưới lồng 3 cấp tự động di chuyển theo tâm bão với độ phân giải ngang lần lượt là 36 km (miền ngoài d01), 12 km (miền giữa d02) và 4 km (miền trong d03).
  • Phương pháp tổ hợp đa vật lý: Xây dựng hệ thống tổ hợp gồm các thành phần kết hợp khác nhau giữa các sơ đồ tham số hóa vật lý:
    • Vi vật lý mây: WSM6, Thompson, Lin, Morrison, WDM6.
    • Đối lưu quy mô dưới lưới: Kain-Fritsch, Betts-Miller-Janjic (BMJ), Grell-Devenyi.
    • Lớp biên hành tinh (PBL): YSU, Mellor-Yamada-Janjic (MYJ), ACM2.
    • Bức xạ sóng ngắn/sóng dài: RRTM, Dudhia, CAM.
  • Phương pháp đồng hóa tổ hợp WRF-LETKF: Tích hợp module vinit để tạo trường quan trắc xoáy ba chiều nhân tạo từ dữ liệu TCVital, sau đó đồng hóa đồng thời cùng số liệu gió vệ tinh AMV vào mô hình WRF-ARW với quy mô tổ hợp 21 thành phần.
  • Kỹ thuật thống kê và đánh giá sai số:
    • Sai số vị trí tâm bão DPE (Distance Position Error, tính bằng km).
    • Sai số trung bình tuyệt đối MAE của gió cực đại bề mặt ($V_{max}$, m/s hoặc kt) và áp suất cực tiểu tại tâm ($P_{min}$, hPa).
    • Đánh giá dự báo xác suất tổ hợp: Biểu đồ hạng Talagrand (Rank histogram), độ tán tổ hợp (Spread) và chỉ số Brier Score (BS).

Nguồn tư liệu và dữ liệu

  • Dữ liệu phân tích bão thời gian thực (TCVital): Cung cấp bởi JTWC, chứa vị trí tâm, $P_{min}$, $V_{max}$, bán kính gió cực đại (RMW) và bán kính các đới gió bão.
  • Dữ liệu quỹ đạo chuẩn (Best-track): Dữ liệu chuẩn quốc tế từ JTWC và Trung tâm Khí tượng Chuyên ngành Khu vực Tokyo (RSMC Tokyo - JMA).
  • Dữ liệu thám sát gió vệ tinh (AMV - Atmospheric Motion Vectors): Thu thập từ Viện Nghiên cứu Vệ tinh Khí tượng CIMSS thuộc Đại học Wisconsin-Madison và NOAA.
  • Dữ liệu tái phân tích toàn cầu NCEP FNL: Dữ liệu phân tích cuối cùng toàn cầu độ phân giải $1^\circ \times 1^\circ$ dùng làm điều kiện ban đầu và điều kiện biên cho các thử nghiệm đa vật lý.
  • Dữ liệu dự báo toàn cầu NCEP GFS: Cung cấp trường phân tích ban đầu và điều kiện biên thời gian thực cho các thử nghiệm đồng hóa tổ hợp.

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

Chương 1: Tổng quan về sai số dự báo bão trên khu vực Tây Bắc Thái Bình Dương

Chương 1 hệ thống hóa cơ sở lý luận về bão, phân loại quy mô (Synop $10^3\text{ km}$, vi mô, quy mô vừa), các tham số cường độ ($V_{max}$, $P_{min}$) và các kỹ thuật quan trắc hiện đại (Dvorak, vệ tinh địa tĩnh GOES, vi sóng microwave, gió AMV). Chương phân tích ưu nhược điểm của 3 nhóm phương pháp dự báo: Synop, Thống kê (CLIPER, STIPS, SHIPS) và Động lực.

Dựa trên lý thuyết bất định Lorenz (1960s, 1970s), chương làm rõ 3 nguồn phát sinh sai số mô hình và 3 nhân tố chi phối sai số cường độ bão. Luận án tiến hành tổng hợp, đánh giá sai số dự báo bão trên Biển Đông giai đoạn 2008–2014 đối với 62 cơn bão từ các nguồn: Việt Nam, RSMC Tokyo, JTWC và hai mô hình số trị GFS, GSM.

Kĩ năng dự báo hạn 48h trên Biển Đông (2008-2014):
- Dự báo quỹ đạo: Tăng từ 20-30% (trước 2010) lên 50-60% (sau 2010), trung bình tăng 4-5%/năm.
- Dự báo cường độ: Dao động thấp ở mức 10-12%, cá biệt năm 2013 ghi nhận kĩ năng mang giá trị âm.

Chương khẳng định nghịch lý tồn tại trên toàn cầu thông qua các công trình của DeMaria (2007, 2014) và dự án HFIP (Gall và cộng sự 2013): tốc độ cải thiện sai số cường độ bão hằng năm $(-0,1\text{ kt/năm}$, tức 0,3–0,5%/năm) chậm hơn gần một bậc đại lượng so với tốc độ cải thiện chất lượng dự báo quỹ đạo (3–5%/năm).

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

Chương 2 trình bày chi tiết cơ sở toán học và công nghệ thực nghiệm:

  • Thiết lập công thức giải bài toán đồng hóa biến phân và thuật toán lọc Kalman tổ hợp chuyển dạng địa phương (LETKF) trong không gian trạng thái tổ hợp.
  • Phát triển module vinit nhằm xây dựng cấu trúc xoáy đối xứng 3 chiều từ thông tin TCVital dựa trên nghiệm giải tích cân bằng của Kieu và Zhang (2009). Cấu trúc này tạo ra trường phân bố gió tiếp tuyến, gió hướng tâm và nhiễu động khí áp tương thích với phương trình động lực của mô hình WRF.
  • Thiết lập quy trình công nghệ kết hợp đồng thời (blending) dữ liệu gió thám sát quy mô lớn AMV và trường quan trắc xoáy quy mô bão từ vinit qua sơ đồ WRF-LETKF.
  • Mô tả cấu hình hệ thống WRF-ARW 3 cấp lưới lồng di chuyển (36 km, 12 km, 4 km) cùng các tham số đánh giá thống kê DPE, MAE, Spread, Biểu đồ hạng và Brier Score.

Chương 3: Khảo sát tương quan giữa sai số dự báo quỹ đạo và cường độ bão bằng hệ thống tổ hợp đa vật lý

Chương 3 thiết lập thử nghiệm 92 trường hợp mô phỏng bão trên khu vực Tây Bắc Thái Bình Dương giai đoạn 2007–2010 bằng mô hình WRF-ARW với điều kiện biên và trường ban đầu từ số liệu tái phân tích NCEP FNL. Mục đích là khảo sát mối quan hệ định lượng giữa sai số quỹ đạo và cường độ bão trong điều kiện cô lập các sai số môi trường.

Luận án phân loại tập mẫu thành hai nhóm tiêu chuẩn lọc sai số quỹ đạo:

  • Tiêu chuẩn I: Lọc ra 30 trường hợp mô phỏng có sai số quỹ đạo nhỏ tại các hạn 24h, 48h và 72h.
  • Tiêu chuẩn II: Lọc khắt khe hơn, chọn ra 16 trường hợp có sai số quỹ đạo nhỏ nhất và quỹ đạo sát với Best-track.
Tập mẫu thử nghiệm Số trường hợp MAE $V_{max}$ 24h (m/s) MAE $V_{max}$ 48h (m/s) MAE $V_{max}$ 72h (m/s) MAE $P_{min}$ 24h (hPa) MAE $P_{min}$ 48h (hPa) MAE $P_{min}$ 72h (hPa)
Toàn bộ tập mẫu 92 8,2 9,8 11,5 14,2 16,8 19,5
Đạt Tiêu chuẩn I 30 7,6 8,1 8,9 13,5 14,1 15,2
Đạt Tiêu chuẩn II 16 7,4 7,3 7,8 13,1 12,6 13,4

Bảng 3: Sai số trung bình tuyệt đối (MAE) của $V_{max}$ và $P_{min}$ theo các tiêu chuẩn lọc quỹ đạo trong thử nghiệm Chương 3 (Nguồn: Tổng hợp từ kết quả mô phỏng của luận án).

Kết quả chỉ ra rằng khi sai số quỹ đạo giảm (từ toàn bộ 92 mẫu sang Tiêu chuẩn I và Tiêu chuẩn II), sai số cường độ ($V_{max}$ và $P_{min}$) giảm rõ rệt ở các hạn 48h và 72h. Tuy nhiên, ở hạn 24h mức độ giảm không đáng kể do bị chi phối bởi sai số trường ban đầu của số liệu FNL (áp suất tâm nông hơn và gió cực đại bề mặt yếu hơn thực tế) cùng thời gian thích ứng động lực (spin-up) của mô hình.

Chương 4: Dự báo quỹ đạo và cường độ bão bằng hệ thống đồng hóa tổ hợp

Chương 4 thử nghiệm hệ thống đồng hóa tổ hợp WRF-LETKF (21 thành phần tổ hợp) kết hợp module vinit và dữ liệu gió vệ tinh AMV cho 6 cơn bão mạnh thực tế (Usagi, Nari, Krosa năm 2013; Rammasun, Neoguri, Vongfong năm 2014), so sánh giữa hai trường hợp:

  • CTRL: Chạy dự báo tổ hợp chuẩn xuất phát từ trường phân tích toàn cầu GFS không đồng hóa xoáy bão.
  • DABV: Đồng hóa đồng thời thông tin xoáy bão nhân tạo (vinit từ TCVital) và dữ liệu gió AMV qua LETKF.

Kết quả phân tích và dự báo ghi nhận:

  • Trường phân tích ban đầu: Gia số phân tích (analysis increments) của các thành phần gió $u, v$ và áp suất mực biển được điều chỉnh tối ưu, thể hiện rõ cấu trúc mắt bão, thành mắt bão và liên kết hoàn lưu trơn tru với trường synop quy mô lớn mà không gây sốc động lực.
  • Kết quả dự báo quỹ đạo và cường độ: Thử nghiệm DABV vượt trội so với CTRL tại mọi hạn dự báo từ 24h đến 120h:
    • Giảm sai số khoảng cách tâm bão (DPE) từ 15–30% ở các hạn 48h–120h.
    • Giảm sai số MAE của $V_{max}$ từ 2–4 m/s và MAE của $P_{min}$ từ 5–10 hPa, đặc biệt khắc phục hiện tượng mô hình dự báo cường độ bão quá yếu ở các hạn đầu.
  • Kĩ năng xác suất tổ hợp: Độ tán tổ hợp (spread) trong DABV bao quát tốt sai số thực tế; biểu đồ hạng Talagrand có dạng chữ U nông và phẳng hơn thể hiện độ tin cậy phân bố xác suất cao; điểm số Brier Score (BS) cho dự báo xác suất quỹ đạo tâm bão được cải thiện đáng kể.

Kết quả và những đóng góp mới

Các đóng góp mới về lý luận và khoa học

  1. Khảo sát định lượng mối quan hệ giữa sai số quỹ đạo và cường độ bão: Luận án đã xây dựng thử nghiệm số trị quy mô lớn (92 trường hợp) trên nền tảng mô hình WRF-ARW để cô lập và định lượng mối liên hệ giữa sai số quỹ đạo và sai số cường độ. Kết quả chứng minh việc cải thiện chất lượng dự báo quỹ đạo đem lại hiệu ứng tích cực trực tiếp làm giảm sai số dự báo cường độ bão, đặc biệt ở các hạn 48–72 giờ. Mức độ tương quan giữa hai loại sai số là không đồng nhất giữa các hạn dự báo do ảnh hưởng của sai số trường ban đầu và thời gian thích ứng động lực ở hạn ngắn.
  2. Phát triển công nghệ đồng hóa kết hợp đa quy mô (Blending Data Assimilation): Luận án đã phát triển thành công module vinit để khởi tạo cấu trúc xoáy 3 chiều giải tích từ báo cáo phân tích bão thời gian thực TCVital, đồng thời thiết lập quy trình đồng hóa tích hợp đồng thời trường xoáy quy mô bão và trường gió vệ tinh quy mô lớn AMV thông qua sơ đồ WRF-LETKF. Phương pháp này cho phép hiệu chỉnh cấu trúc ban đầu của bão một cách khách quan, tương thích động lực học mà không cần can thiệp cứng vào phương trình mô hình.

Đóng góp thực tiễn và giải pháp đề xuất

  • Cung cấp một công cụ công nghệ hoàn chỉnh kết hợp giữa WRF-ARW, vinit và LETKF có khả năng ứng dụng trực tiếp vào nghiệp vụ dự báo bão thời gian thực tại Trung tâm Dự báo Khí tượng Thủy văn Quốc gia.
  • Đề xuất giải pháp sử dụng số liệu TCVital và AMV để giải quyết bài toán thiếu hụt thám sát khí tượng trên biển, nâng cao độ chính xác của các bản tin cảnh báo bão sớm hạn 3–5 ngày phục vụ phòng chống thiên tai tại Việt Nam.

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

Hạn chế của luận án

  • Cấu trúc xoáy nhân tạo khởi tạo từ module vinit vẫn dựa trên giả thiết phân bố đối xứng giải tích, chưa phản ánh đầy đủ tính chất bất đối xứng trong cấu trúc hoàn lưu thực tế của các cơn bão có tương tác mạnh với địa hình hoặc đứt gió môi trường lớn.
  • Số lượng thử nghiệm đồng hóa tổ hợp trong Chương 4 mới tập trung vào 6 cơn bão mạnh điển hình giai đoạn 2013–2014, chưa bao phủ đầy đủ các dạng bão yếu, áp thấp nhiệt đới hoặc các cơn bão có quỹ đạo phức tạp di chuyển gần bờ.
  • Chưa tích hợp các nguồn dữ liệu thám sát cục bộ có độ phân giải cao như số liệu radar thời tiết Doppler ven biển hoặc thám sát vô tuyến bổ sung (dropsonde).

Hướng nghiên cứu tiếp theo

  • Hoàn thiện module vinit theo hướng xây dựng cấu trúc xoáy bất đối xứng 3 chiều dựa trên trường gió môi trường thực tế.
  • Mở rộng thử nghiệm đồng hóa tổ hợp WRF-LETKF cho toàn bộ các cơn bão hoạt động trên Biển Đông trong nhiều mùa bão liên tiếp để đánh giá tính ổn định nghiệp vụ.
  • Tích hợp thêm các nguồn dữ liệu thám sát vi sóng vệ tinh trực tiếp, dữ liệu radar thời tiết ven biển Việt Nam và tối ưu hóa bán kính địa phương hóa (localization) trong thuật toán LETKF.

Giá trị tham khảo

Luận án mang lại giá trị tham khảo học thuật và thực tiễn chuyên sâu:

  • Nghiên cứu sinh, học viên cao học và giảng viên ngành Khí tượng, Khí hậu học và Khoa học Khí quyển: Là tài liệu tham khảo chi tiết về cơ sở lý thuyết đồng hóa số liệu Bayes, thuật toán lọc Kalman tổ hợp địa phương (LETKF), phương pháp tổ hợp đa vật lý và cơ chế động lực học xoáy thuận nhiệt đới.
  • Dự báo viên và kỹ sư mô hình số trị: Cung cấp quy trình công nghệ chuẩn từ tiền xử lý dữ liệu viễn thám vệ tinh (AMV), giải mã điện văn thời gian thực (TCVital), cài đặt hệ thống WRF-LETKF đến các phương pháp đánh giá sai số xác suất tổ hợp (Rank histogram, Brier Score).
  • Cơ quan quản lý và quy hoạch phòng chống thiên tai: Cung cấp căn cứ khoa học định lượng về độ tin cậy và biên độ sai số của các sản phẩm dự báo bão theo từng hạn dự báo (24h, 48h, 72h, 120h).

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

1. Tại sao tốc độ cải thiện chất lượng dự báo cường độ bão lại diễn ra chậm hơn rất nhiều so với dự báo quỹ đạo bão?

Theo kết quả tổng hợp trong Chương 1, dự báo quỹ đạo bão chủ yếu phụ thuộc vào hoàn lưu dòng dẫn đường quy mô lớn (Synop $10^3\text{ km}$), vốn được các mô hình toàn cầu phân giải cao mô phỏng ngày càng chính xác. Ngược lại, cường độ bão bị chi phối bởi các quá trình tương tác phức tạp từ quy mô vi mô (vật lý mây), quy mô vừa (đối lưu sâu, thành mắt bão) đến quy mô lớn, kèm theo các hạn chế về độ phân giải lưới, sơ đồ tham số hóa vật lý chưa hoàn hảo và tình trạng thiếu số liệu thám sát cấu trúc bên trong mắt bão.

2. Ba nhân tố chính nào dẫn đến sai số mô phỏng cường độ bão trong mô hình số trị?

Ba nhân tố được luận án tổng kết gồm:

  1. Thông tin ban đầu về cấu trúc xoáy bão bị hạn chế và thô do thiếu quan trắc thám sát trực tiếp.
  2. Các sơ đồ tham số hóa vật lý quy mô dưới lưới (vi vật lý, đối lưu, lớp biên hành tinh, bức xạ) chưa phản ánh chính xác các quá trình nhiệt động bên trong bão.
  3. Sai số của trường môi trường xung quanh do dự báo sai quỹ đạo bão dẫn đến điều kiện nhiệt ẩm và đứt gió mô phỏng không thực tế.

3. Kết quả khảo sát trong Chương 3 cho thấy mối quan hệ giữa sai số quỹ đạo và cường độ bão ở hạn 24h khác biệt như thế nào so với hạn 48h và 72h?

Ở các hạn 48h và 72h, khi sai số quỹ đạo giảm rõ rệt (qua các tập lọc Tiêu chuẩn I và II), sai số cường độ ($V_{max}$ và $P_{min}$) giảm tương ứng một cách có ý nghĩa thống kê do môi trường xung quanh bão được mô phỏng chính xác hơn. Tuy nhiên, ở hạn 24h, việc giảm sai số quỹ đạo không mang lại sự cải thiện rõ rệt cho cường độ bão do sai số ở hạn này bị chi phối chủ yếu bởi sai lệch cấu trúc ban đầu của trường FNL và thời gian thích ứng động lực (spin-up) của mô hình.

4. Module vinit kết hợp với thuật toán LETKF xử lý bài toán trường ban đầu cho mô hình WRF-ARW như thế nào?

Module vinit trích xuất các thông tin thời gian thực từ điện văn TCVital ($V_{max}$, $P_{min}$, vị trí tâm, bán kính gió) để tạo ra trường quan trắc xoáy 3 chiều giải tích cân bằng. Thuật toán LETKF sau đó đồng hóa đồng thời trường xoáy này cùng dữ liệu gió vệ tinh AMV vào trường nền của mô hình WRF-ARW, giúp cập nhật ma trận hiệp phương sai sai số phụ thuộc dòng chảy, tạo ra trường phân tích tối ưu cân bằng động lực mà không gây hiện tượng sốc mô hình.

5. Kĩ năng dự báo bão hạn 48h trên Biển Đông giai đoạn 2008–2014 thể hiện sự chênh lệch ra sao giữa quỹ đạo và cường độ?

Đánh giá trên 62 cơn bão tại Biển Đông cho thấy kĩ năng dự báo quỹ đạo hạn 48h của các trung tâm và mô hình động lực tăng mạnh từ 20–30% (giai đoạn trước 2010) lên 50–60% (giai đoạn sau 2010), đạt tốc độ tăng trung bình 4–5%/năm. Trong khi đó, kĩ năng dự báo cường độ hạn 48h chỉ đạt mức thấp 10–12% so với phương pháp thống kê khí hậu CLIPER, thậm chí trong mùa bão năm 2013 nhiều trung tâm quốc tế ghi nhận giá trị kĩ năng âm.

Kết luận

Luận án của tác giả Dư Đức Tiến đã giải quyết thành công bài toán khảo sát định lượng mối quan hệ giữa sai số mô phỏng quỹ đạo và cường độ bão trên khu vực Tây Bắc Thái Bình Dương bằng hệ thống mô hình số trị WRF-ARW. Kết quả nghiên cứu khẳng định việc nâng cao chất lượng dự báo quỹ đạo có tác động tích cực trực tiếp đến việc giảm sai số dự báo cường độ ở các hạn từ 2 đến 3 ngày, đồng thời chỉ ra tính chất phi đồng nhất của mối tương quan này ở các hạn dự báo ngắn do ảnh hưởng của trường phân tích ban đầu. Việc phát triển thành công module khởi tạo xoáy ba chiều vinit kết hợp với phương pháp đồng hóa tổ hợp lọc Kalman địa phương (WRF-LETKF) và dữ liệu gió vệ tinh AMV đã tạo ra một giải pháp công nghệ hiệu quả, nâng cao chất lượng dự báo xác suất cho cả quỹ đạo và cường độ bão, đóng góp thiết thực cho công tác nghiên cứu và nghiệp vụ khí tượng thủy văn tại Việt Nam.