Tổng quan về luận án

Xoáy thuận nhiệt đới (XTNĐ) hình thành và hoạt động trên vùng biển Tây Bắc Thái Bình Dương (TBTBD) chiếm khoảng 38% tổng số cơn bão toàn cầu với tần suất xấp xỉ 30 cơn/năm, tạo ra áp lực thiên tai thường trực đối với dải duyên hải Việt Nam. Theo thống kê của Tổ chức Khí tượng Thế giới (WMO), "trong 50 năm qua, có 1.942 trường hợp thiên tai do XTNĐ gây ra, làm gần 800.000 người thiệt mạng và gây thiệt hại kinh tế 407,6 tỷ USD, trung bình mỗi ngày có 43 người chết và thiệt hại 78 triệu USD". Tại Việt Nam, công tác dự báo quỹ đạo và cường độ bão khi tiếp cận vùng bờ biển luôn đối mặt với sai số lớn do các quá trình tương tác phi tuyến phức tạp giữa hoàn lưu xoáy bão quy mô vừa với địa hình chắn bờ biển (đặc biệt là dãy núi Trường Sơn và địa hình đồi núi Bắc Bộ) và các đợt xâm nhập của không khí lạnh (KKL) cực đới biến tính từ Áp cao Siberia trong thời kỳ chuyển mùa (từ tháng 9 đến tháng 12).

Khoảng trống nghiên cứu (research gap) trọng tâm nằm ở chỗ các nghiên cứu trước đây tại Việt Nam chủ yếu tiếp cận theo phương pháp thống kê synốp định tính hoặc áp dụng các mô hình chính áp đơn giản (WBAR, BoLAM), chưa định lượng hóa chi tiết cơ chế biến đổi ba chiều của cấu trúc mắt bão, thành mắt bão, dải mây đối lưu xoắn và phân bố trường mưa, trường gió khi bão chịu tác động đồng thời của địa hình gồ ghề và sự thâm nhập nêm lạnh tầng thấp. Luận án tiến sĩ chuyên ngành Khí tượng và Khí hậu học (Mã số: 9440222) của NCS. Nguyễn Bình Phong với đề tài "Nghiên cứu ảnh hưởng của địa hình và không khí lạnh đến cấu trúc của bão ở vùng duyên hải Việt Nam" dưới sự hướng dẫn của PGS. Nguyễn Văn Thắng và TS. Nguyễn Văn Hiệp tại Viện Khoa học Khí tượng Thủy văn và Biến đổi khí hậu (2022) đã giải quyết triệt để bài toán này.

Hệ thống câu hỏi nghiên cứu và giả thuyết khoa học được thiết lập chặt chẽ:

  1. Câu hỏi nghiên cứu 1 (RQ1): Địa hình vùng duyên hải Việt Nam làm biến dạng cấu trúc đối xứng của các trường khí tượng (áp, gió, ẩm, phản hồi vô tuyến) và làm lệch quỹ đạo của xoáy bão như thế nào trước và sau khi đổ bộ?
    • Giả thuyết 1 (H1): Hiệu ứng ma sát và nâng cưỡng bức của sườn đón gió gây bất đối xứng mạnh trong trường gió mực 10m và tập trung tỉ số xáo trộn nước mưa ở phía trước tâm bão, đồng thời lực cản địa hình làm chậm vận tốc di chuyển và gây hiệu ứng lệch hướng quỹ đạo về phía nam.
  2. Câu hỏi nghiên cứu 2 (RQ2): Không khí lạnh xâm nhập vào hoàn lưu bão trên Biển Đông theo những cơ chế nhiệt động lực học nào và làm thay đổi cấu trúc đối lưu, cường độ cực tiểu của khí áp tâm bão ra sao?
    • Giả thuyết 2 (H2): Sự xâm nhập trực tiếp của KKL khô lạnh tạo thành nêm lạnh ở tầng biên làm triệt tiêu nguồn cung cấp ẩm và phá vỡ cấu trúc lõi nóng; ngược lại, tác động gián tiếp từ xa làm gia tăng gradient khí áp nằm ngang, mở rộng trường gió mạnh ở rìa phía bắc.
  3. Câu hỏi nghiên cứu 3 (RQ3): Phương pháp ban đầu hóa xoáy động lực học kết hợp dịch chuyển xoáy cải thiện chất lượng trường ban đầu trong mô hình số trị quy mô vừa WRF-ARW ở mức độ nào?
    • Giả thuyết 3 (H3): Tinh chỉnh cấu trúc xoáy ban đầu bằng sơ đồ ban đầu hóa xoáy chuyên biệt giúp tái lập chính xác gradient khí áp xuyên tâm và dị thường nhiệt độ thẳng đứng, giảm thiểu sai số dự báo vị trí và cường độ trong 24 giờ đầu.

Khung lý thuyết của luận án được xây dựng dựa trên động lực học khí quyển phi thủy tĩnh, lý thuyết lớp biên Ekman, cơ chế bất ổn định có điều kiện của khí quyển (Conditional Instability of the Second Kind - CISK), lý thuyết cân bằng gió nhiệt và quá trình truyền nhiệt - ẩm đối lưu sâu. Phạm vi nghiên cứu bao gồm toàn bộ khu vực Biển Đông và vùng duyên hải Việt Nam, khảo sát chuỗi số liệu thực tế giai đoạn 2014-2018 với tập mẫu 18 cơn bão nhiệt đới điển hình và thực hiện tổ hợp 228 thí nghiệm mô phỏng số trị trên mô hình WRF.

Literature Review và Positioning

Lịch sử nghiên cứu cấu trúc bão quốc tế khởi nguồn từ các phép đo thám sát máy bay đầu tiên của Duckworth và O'Hair năm 1943, tiếp nối bằng mô hình mặt cắt thẳng đứng kinh điển của Palmén (1948) và Riehl (1948, 1950, 1954). Về mặt nhiệt động lực học, Kleinschmidt (1951) đã đặt nền móng với công thức biểu diễn vận tốc gió cực đại $V$ là hàm số của năng lượng nhiệt đại dương:

$$V = \sqrt{2E\left(\frac{q^2}{1 - q^2}\right)}$$

Trong đó $E$ là năng lượng nhiệt đại dương trên một đơn vị khối lượng và $q$ là hệ số thực nghiệm. Tiếp sau đó, Miller (1958) phát triển lý thuyết áp suất cực tiểu tại tâm bão dựa trên phương trình thủy tĩnh và tích phân profile nhiệt độ thẳng đứng. Trong kỷ nguyên hiện đại, Hawkins và Imbembo (1976), Holland và Merrill (1984), Emanuel và Rotunno (2011), Chavas và Emanuel (2014) đã hoàn thiện lý thuyết cấu trúc đối xứng của trường gió phương vị và profile áp suất bề mặt tuân theo hàm phân bố log-normal.

Về ảnh hưởng của địa hình đến bão, các nghiên cứu quốc tế đã ghi nhận sự chuyển hướng và suy giảm cấu trúc bão khi tiếp cận các địa hình đồi núi phức tạp. Yeh và Elsberry (1993) khi khảo sát 103 cơn bão tiếp cận dãy núi Trung ương Đài Loan (CMR) đã chỉ ra hiện tượng suy yếu cường độ và giảm tốc độ di chuyển. Jian và Wu (2008) sử dụng mô hình WRF độ phân giải 4 km nghiên cứu siêu bão Haitang (2005) đã chứng minh cơ chế hình thành chuyển động thắt nút (looping track) và lệch hướng nam do sự hội tụ gió mực thấp bị cưỡng bức giữa xoáy bão và địa hình núi cao. Huang và Yang (2014) khi mô phỏng bão Morakot (2009) khẳng định địa hình làm giảm hiệu quả mưa từ 15-20% sau khi bão vượt qua sườn núi khuất gió.

Về tương tác giữa bão và không khí lạnh, tồn tại cuộc tranh luận học thuật sâu sắc giữa hai trường phái:

  • Quan điểm tăng cường (Intensification View): Fett (1968), Li Tianfu (2003) và đặc biệt là Han Ying et al. (2008) chứng minh rằng khi gradient nhiệt độ nằm ngang gia tăng do sự tiếp cận của khối KKL ở rìa bắc, lực đẩy vào trung tâm sẽ làm tăng cường độ hội tụ tầng thấp, gia tăng vận chuyển năng lượng và nâng cao cường độ xoáy bão.
  • Quan điểm suy yếu và triệt tiêu (Decay / Extratropical Transition View): Riehl (1963b), Palmén và Newton (1969), Lê Văn Thảo (1996), Trần Đình Bá (1979, 1997), Nguyễn Văn Thắng (2011) lại khẳng định sự xâm nhập của khối không khí cực đới khô lạnh đóng vai trò như một "nêm lạnh" chui xuống phần đáy cột không khí, cắt đứt hoàn lưu ẩm bốc hơi từ đại dương, làm tiêu tán lõi nóng và khiến bão suy yếu nhanh chóng hoặc chuyển hóa thành xoáy thuận ngoại nhiệt đới.

Tại Việt Nam, các công trình của Phan Văn Tân và Bùi Hoàng Hải (2004), Phan Văn Tân và Nguyễn Lê Dũng (2008), Hoàng Đức Cường và cs (2010), Trần Tân Tiến và cs (2009, 2010), Dư Đức Tiến và cs (2016) đã bước đầu ứng dụng mô hình MM5, WRF, đồng hóa dữ liệu 3D-VAR và tổ hợp lọc Kalman (LetKF) để dự báo quỹ đạo. Tuy nhiên, việc bóc tách định lượng mức độ ảnh hưởng độc lập và tương hỗ của độ cao địa hình ven biển và các đợt KKL đến cấu trúc nhiệt động học ba chiều của bão vẫn chưa được phân tích toàn diện. Luận án của NCS. Nguyễn Bình Phong định vị chính xác vào khoảng trống này, kết nối phương pháp mô phỏng số độ phân giải cao với các kịch bản thực nghiệm thay đổi độ cao địa hình để cung cấp bằng chứng động lực học xác thực.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án mở rộng lý thuyết nhiệt động lực học bão của Kleinschmidt (1951) và Emanuel và Rotunno (2011) trong điều kiện biên bất đồng nhất bằng cách xác lập mô hình tương tác đa quy mô giữa xoáy thuận nhiệt đới với địa hình chắn và front lạnh cực đới. Những đóng góp lý thuyết cốt lõi được cụ thể hóa qua hệ thống các mệnh đề khoa học:

  1. Mệnh đề 1 (P1 - Về sự phá vỡ cấu trúc đối xứng do địa hình): Lực cản ma sát và hiệu ứng nâng cơ học của địa hình duyên hải phá vỡ tính đối xứng phương vị của trường gió tiếp tuyến và trường ẩm. Tại sườn đón gió, gia tốc hội tụ Ekman cưỡng bức làm tăng cực đại tỉ số xáo trộn ngưng kết và độ phản hồi vô tuyến radar ($\ge 45\text{ dBZ}$) ở góc phần tư phía trước của bão, trong khi sườn khuất gió chịu hiệu ứng fohn làm tiêu tán các dải mây xoắn phía sau.
  2. Mệnh đề 2 (P2 - Về cơ chế phân hóa tác động của không khí lạnh): Bản chất tác động của KKL lên cấu trúc bão phụ thuộc nghiêm ngặt vào chỉ số tương phản nhiệt độ bề mặt $\Delta T$ và khoảng cách tương tác. Như Trần Đình Bá (1979) đã chỉ ra: "Nếu tương phản nhiệt front vùng gần XTNĐ lớn (> 6°C) và nhiệt độ không khí hút vào XTNĐ thấp (<27°C) thì 24h sau đó XTNĐ sẽ yếu hẳn đi. Nếu tương phản nhiệt front nhỏ ($\le 6^\circ\text{C}$) nhiệt độ không khí hút vào XTNĐ cao ($\ge 27^\circ\text{C}$) thì XTNĐ tiếp tục phát triển hoặc ít thay đổi cường độ". Luận án bổ sung cơ chế động lực: KKL gián tiếp làm tăng gradient áp lực ngoại vi dẫn đến mở rộng bán kính gió mạnh cấp 6-8, trong khi KKL trực tiếp làm mất cân bằng gió nhiệt và triệt tiêu đối lưu sâu ở thành mắt bão.
  3. Mệnh đề 3 (P3 - Về sự lệch pha quỹ đạo do hiệu ứng hội tụ cưỡng bức): Sự tương tác giữa hoàn lưu xoáy bão và sườn núi hướng Tây Bắc - Đông Nam của dải Trường Sơn tạo ra một vùng dị thường áp thấp cục bộ và tập trung mômen động lượng ở phía nam tâm bão, giải thích hiện tượng bão có xu thế uốn cong quỹ đạo lệch về phía nam trước thời điểm đổ bộ.

Khung phân tích độc đáo

Khung phân tích của nghiên cứu là sự tích hợp chặt chẽ giữa ba trường phái lý thuyết: (1) Lý thuyết động lực học xoáy lớp biên Ekman, (2) Lý thuyết nhiệt động lực học đối lưu mây ẩm (Moist Convective Thermodynamics), và (3) Lý thuyết front học quy mô vừa (Mesoscale Frontogenesis).

Phương pháp tiếp cận mới nằm ở việc thiết kế ma trận độ nhạy địa hình lũy tiến (Topographic Sensitivity Matrix) kết hợp khảo sát biến trình thẳng đứng của trường dị thường nhiệt độ ($\Delta T$), véc-tơ gió tiếp tuyến, và tỉ số xáo trộn pha ngưng kết ($q_c, q_r, q_i$). Điều kiện biên được kiểm soát chặt chẽ trên lưu vực Biển Đông, giới hạn cho các cơn bão có cường độ từ cấp bão nhiệt đới (TS) đến bão rất mạnh (TY/STY) theo phân loại tại Quyết định số 44/2014/QĐ-TTg của Thủ tướng Chính phủ.

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

Thiết kế nghiên cứu

Luận án tuân thủ thế giới quan thực chứng (Positivism) và chủ nghĩa hiện thực phản biện (Critical Realism), dựa trên hệ phương trình vi phân đạo hàm riêng phi thủy tĩnh mô tả hoàn lưu khí quyển. Nghiên cứu sử dụng thiết kế đa tầng bậc (multi-level design) kết hợp giữa phương pháp mô phỏng số trị quy mô vừa, phương pháp phân tích synốp chuẩn mực và phương pháp thống kê kiểm định sai số.

Mô hình số trị được lựa chọn là mô hình quy mô vừa thế hệ mới WRF-ARW (Advanced Research WRF), tích hợp mô-đun đồng hóa dữ liệu WRFDA và kỹ thuật ban đầu hóa xoáy động lực học. Miền tính được thiết kế lồng ghép hai chiều (nested domains) bao phủ toàn bộ Biển Đông và bán đảo Đông Dương với độ phân giải lưới tinh bảo đảm mô phỏng chi tiết các quá trình động lực học thành mắt bão.

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu bao gồm các bước kiểm chuẩn nghiêm ngặt:

  1. Thu thập và tiền xử lý dữ liệu:
    • Số liệu tái phân tích toàn cầu GFS của NCEP với độ phân giải $0.5^\circ \times 0.5^\circ$ hoặc $1.0^\circ \times 1.0^\circ$ tại các thời điểm 00:00 UTC và 12:00 UTC làm điều kiện ban đầu và điều kiện biên.
    • Số liệu quan trắc quỹ đạo và cường độ bão thực tế khách quan (best-track) từ Trung tâm Cảnh báo Bão Liên hợp Hải quân Hoa Kỳ (JTWC), Cơ quan Khí tượng Nhật Bản (JMA) và Trung tâm Dự báo Khí tượng Thủy văn Quốc gia Việt Nam (NCHMF).
    • Số liệu viễn thám: Ảnh mây vệ tinh đa phổ (VIS, IR), kênh vi ba 89H, và dữ liệu phản xạ radar Doppler NEXRAD / trạm radar thời tiết ven biển Việt Nam.
  2. Kỹ thuật Ban đầu hóa xoáy (Vortex Initialization):
    • Tách lọc hoàn lưu xoáy quy mô lớn khỏi trường nền quy mô synốp.
    • Hiệu chỉnh vị trí tâm xoáy (vortex relocation) và đưa vào xoáy nhân tạo cân bằng động lực học (bogus vortex) dựa trên quan trắc khí áp trung tâm và bán kính gió cực đại ($R_{\max}$).
    • Đồng hóa số liệu quan trắc gió thám sát trên cao và gió chuyển động mây (cloud-drift winds) từ Đại học Wisconsin (Mỹ) thông qua thuật toán 3D-VAR.
  3. Thiết kế ma trận thí nghiệm độ cao địa hình:
    • CTL: Thí nghiệm đối chứng giữ nguyên độ cao địa hình thực tế của bán đảo Đông Dương.
    • TER2m: Giảm toàn bộ độ cao địa hình đất liền xuống 2m (loại bỏ hoàn toàn ảnh hưởng cơ học của núi).
    • TER50: Giảm 50% độ cao địa hình so với thực tế.
    • TER75: Giảm 75% độ cao địa hình.
    • TER150: Gia tăng 150% độ cao địa hình thực tế nhằm khuếch đại phản ứng động lực học của hoàn lưu bão.

Data và phân tích

Tập dữ liệu nghiên cứu bao gồm 18 cơn bão hoạt động trên Biển Đông trong giai đoạn 5 năm (2014-2018), chia thành hai phân nhóm rõ rệt: nhóm bão đơn thuần đổ bộ và nhóm bão chịu tác động kết hợp của KKL mùa đông (tiêu biểu như bão Kalmaegi 2014, bão Mujigae 2015, bão Wutip 2013, bão Damrey 2017, bão Usagi 2018).

Tổng cộng 228 thí nghiệm mô phỏng số trị đã được thực hiện và phân tích thống kê. Các chỉ số đánh giá độ chính xác bao gồm sai số trung bình (Mean Error - ME), sai số trung bình tuyệt đối (Mean Absolute Error - MAE), và độ lệch chuẩn đối với các biến:

  • Sai số vị trí tâm bão (km) tại các hạn dự báo 06h, 12h, 18h, 24h, 36h, 48h, 72h.
  • Sai số khí áp cực tiểu mực biển (Sea Level Pressure - SLP, đơn vị hPa).
  • Sai số vận tốc gió cực đại mực 10m ($V_{\max}$, đơn vị m/s).

Phát hiện đột phá và implications

Những phát hiện then chốt

Nghiên cứu mang lại 4 phát hiện đột phá về mặt cấu trúc và động lực học khí tượng:

  1. Hiệu ứng tái cấu trúc bất đối xứng trường mưa và gió do sườn núi đón gió: Khi xoáy bão tiến vào vùng biển ven bờ cách bờ biển dưới 150 km, lực ma sát bề mặt gia tăng đột ngột kết hợp với địa hình dâng cao của dãy Trường Sơn kích hoạt quá trình nâng cơ học cưỡng bức. Kết quả mô phỏng mặt cắt Đông - Tây qua tâm bão Damrey (2017) cho thấy tỉ số xáo trộn nước mưa và ngưng kết ($q_r, q_c$) tăng vọt tại sườn đón gió phía Tây và Tây Bắc của tâm bão, tạo ra vùng phản hồi vô tuyến cực đại vượt ngưỡng 45-50 dBZ. Sau khi đổ bộ 3 giờ, hoàn lưu xoáy bão bị cắt đứt chân tiếp xúc với nguồn nhiệt ẩm đại dương, cấu trúc thành mắt bão tan rã nhanh chóng trong thí nghiệm CTL và TER150, nhưng vẫn duy trì cấu trúc xoáy đối xứng nguyên vẹn trong thí nghiệm phẳng TER2m.
  2. Quy luật lệch hướng nam của quỹ đạo bão gần bờ: So sánh quỹ đạo bão trong các kịch bản địa hình chứng minh rằng độ cao địa hình càng lớn thì bão càng di chuyển chậm lại và quỹ đạo bão càng bị uốn cong lệch về phía nam trước thời điểm đổ bộ. Nguyên nhân động lực học được xác lập là do sự tích tụ mômen xoáy và hội tụ gió mạnh mực thấp ở rìa phía bắc giữa hoàn lưu bão và sườn núi, tạo ra một gradient khí áp bất đối xứng đẩy hoàn lưu tâm bão dịch chuyển dần xuống các vĩ độ thấp hơn.
  3. Phân hóa cơ chế tác động của không khí lạnh (Trực tiếp vs Gián tiếp):
    • Tác động trực tiếp: Khi KKL khô lạnh tràn sâu vào trung tâm bão (trường hợp bão Damrey 2017 và Mujigae 2015), một "nêm lạnh" hình thành ở tầng 1000-850 hPa làm sụt giảm nhiệt độ thế vị tương đương $\theta_e$, triệt tiêu dòng thăng thẳng đứng ($w$) trong thành mắt bão và làm dị thường nhiệt độ lõi nóng ở tầng 300-200 hPa giảm từ $+8\text{ K}$ xuống dưới $+2\text{ K}$, dẫn đến quá trình bão đầy lên và suy yếu đột ngột.
    • Tác động gián tiếp: Khi KKL duy trì ở khoảng cách 1000-1500 km phía bắc, tương tác sóng làm gia tăng đột biến gradient khí áp mực biển ($\nabla P$) giữa Áp cao lục địa và tâm bão, mở rộng phạm vi có gió mạnh cấp 6-8 lên tới hàng trăm kilômét ở rìa phía Bắc bão. Đúng như Lê Văn Thảo (1996) đã tổng kết: "Trong sáu trường hợp bão đổi hướng xuống phía Nam có năm trường hợp liên quan đến sự xâm nhập của KKL và một trường hợp nhiễu động sóng đông trên cao ở phía Đông của bão".
  4. Hiệu năng vượt trội của sơ đồ ban đầu hóa xoáy: Việc áp dụng kỹ thuật ban đầu hóa xoáy động lực học giúp loại bỏ hiện tượng xoáy bão quá yếu và lệch vị trí trong dữ liệu GFS nguyên thủy. Phân tích mặt cắt vĩ hướng bão Mujigae (2015) và Damrey (2017) tại thời điểm 00Z cho thấy mô hình có ban đầu hóa xoáy tái lập hoàn hảo cấu trúc phễu sâu của khí áp mực biển và phân bố gió 10m đối xứng quanh tâm, giúp giảm sai số khí áp cực tiểu tại tâm bão từ 15-20 hPa xuống dưới 4 hPa ngay tại hạn dự báo ban đầu.

Implications đa chiều

  • Về lý thuyết khí tượng học: Làm sáng tỏ cơ chế tương tác nhiệt - động lực học giữa xoáy thuận nhiệt đới vùng nhiệt đới với khối không khí cực đới biến tính và địa hình chắn duyên hải, giải quyết dứt điểm các tranh luận học thuật kéo dài về tác động hai mặt của KKL.
  • Về phương pháp luận mô hình số: Xác lập bộ tham số hóa vật lý vi mô, đối lưu (Kain-Fritsch/Betts-Miller-Janjic) và lớp biên hành tinh (YSU) tối ưu kết hợp với sơ đồ ban đầu hóa xoáy cho mô hình WRF tại khu vực Biển Đông Việt Nam.
  • Về ứng dụng thực tiễn và nghiệp vụ: Cung cấp cơ sở khoa học định lượng giúp các dự báo viên tại Trung tâm Dự báo Khí tượng Thủy văn Quốc gia và các Đài Khí tượng Thủy văn khu vực nhận diện chính xác các ổ mưa đặc biệt lớn và vùng gió giật cực đoan ở sườn đón gió duyên hải Trung Bộ.
  • Về chính sách và phòng chống thiên tai: Hỗ trợ Ban Chỉ đạo Quốc gia về Phòng, chống thiên tai hoạch định kịch bản sơ tán dân cư, nâng cấp đê biển và quy hoạch không gian xây dựng ven bờ dựa trên bản đồ phân vùng rủi ro bão tương tác KKL và địa hình.

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, nghiên cứu còn tồn tại một số hạn chế:

  1. Quy mô mẫu thời gian: Dữ liệu khảo sát tập trung chủ yếu trong giai đoạn 5 năm (2014-2018) với 18 cơn bão; mặc dù có tính đại diện cao nhưng chưa bao phủ hết các chu kỳ dao động khí hậu liên năm (ENSO, IOD, MJO).
  2. Cơ chế tương tác biển - khí quyển chưa liên kết hoàn toàn: Mô hình WRF-ARW được chạy dưới dạng mô hình khí quyển đơn lẻ với nhiệt độ bề mặt biển (SST) tĩnh từ dữ liệu phân tích, chưa tích hợp mô hình đại dương ba chiều (như ROMS hoặc POM) để tính toán phản hồi làm lạnh nước biển do bão khuấy đảo (sea surface cooling).
  3. Độ phân giải mô hình đối với vi địa hình: Mặc dù lưới tính đã được tinh chỉnh, các đèo núi hẹp và thung lũng cục bộ của dãy Trường Sơn vẫn bị làm mịn nhất định trong tập số liệu địa hình địa lý số.

Hướng nghiên cứu tiếp theo bao gồm:

  • Mở rộng xây dựng hệ thống dự báo tổ hợp đa mô hình kết hợp đồng hóa dữ liệu radar Doppler ven biển và thám không vô tuyến thời gian thực.
  • Phát triển mô hình liên kết hoàn toàn Khí quyển - Đại dương - Sóng (WRF - ROMS - SWAN) để mô phỏng tương tác ba chiều khi bão tương tác với KKL trên thềm lục địa Biển Đông.
  • Ứng dụng trí tuệ nhân tạo (Deep Learning) trong việc hiệu chỉnh sai số sau mô hình số trị (post-processing) cho các trường gió cực trị và mưa cực đoan do địa hình.

Tác động và ảnh hưởng

Luận án là công trình nghiên cứu hệ thống, chuyên sâu hàng đầu tại Việt Nam về động lực học bão tương tác địa hình và KKL. Về mặt học thuật, công trình đóng góp nguồn dữ liệu thực nghiệm mô phỏng số trị có giá trị cao, tạo tiền đề cho hàng loạt công bố quốc tế trên các tạp chí khí tượng chuyên ngành thuộc danh mục ISI/Scopus.

Về mặt chuyển giao công nghệ, khung mô hình WRF-ARW tích hợp sơ đồ ban đầu hóa xoáy đã được định hướng chuyển giao cho hệ thống dự báo thời tiết nghiệp vụ tại Việt Nam, nâng cao độ tin cậy của bản tin cảnh báo sớm trước 24-48 giờ. Đối với xã hội, việc xác định chính xác trọng tâm mưa lũ cục bộ tại miền Trung giúp giảm thiểu thiệt hại hàng trăm tỷ đồng mỗi năm cho cơ sở hạ tầng giao thông, thủy điện và sản xuất nông nghiệp ven biển.

Đối tượng hưởng lợi

  • Nghiên cứu sinh và học viên cao học ngành Khí tượng, Thủy văn, Biến đổi khí hậu: Tiếp cận khung phương pháp luận mô hình hóa số trị tiên tiến, kỹ thuật ban đầu hóa xoáy và phương pháp thiết kế thí nghiệm độ nhạy địa hình.
  • Các nhà khoa học và chuyên gia động lực học khí quyển: Kế thừa các phát hiện lý thuyết về tương tác sóng - xoáy và cơ chế nhiệt động học nêm lạnh.
  • Cơ quan dự báo Khí tượng Thủy văn Quốc gia (NCHMF) và khu vực: Sở hữu bộ công cụ tham số hóa tối ưu phục vụ trực tiếp công tác dự báo nghiệp vụ bão đổ bộ.
  • Các cơ quan quản lý nhà nước (Bộ Tài nguyên và Môi trường, Bộ Nông nghiệp và Phát triển nông thôn): Có căn cứ khoa học chuẩn xác để xây dựng quy chuẩn kỹ thuật công trình ven biển và chiến lược phòng chống thiên tai quốc gia.

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì? Trả lời: Luận án đã định lượng hóa và mở rộng lý thuyết nhiệt động lực học bão của Kleinschmidt (1951) và Miller (1958) bằng cách giải mã bản chất hai mặt của tương tác giữa KKL và bão nhiệt đới trên Biển Đông: phân định ranh giới vật lý giữa hiệu ứng tăng cường gradient khí áp ngoại vi (tác động gián tiếp từ xa) và hiệu ứng nêm lạnh triệt tiêu đối lưu sâu thành mắt bão (xâm nhập trực tiếp), đồng thời chứng minh bằng mô hình số quy luật địa hình Trường Sơn cưỡng bức lệch quỹ đạo bão về phía nam.

2. Điểm mới về mặt phương pháp nghiên cứu so với các công trình trước đây là gì? Trả lời: So với các nghiên cứu trước đây (như Yeh và Elsberry 1993 sử dụng mô hình lưới thô 45 km, hay Bùi Hoàng Hải và Phan Văn Tân 2002 dùng mô hình WBAR chính áp), luận án đã thiết lập hệ thống thí nghiệm ma trận độ nhạy địa hình 5 cấp độ (CTL, TER2m, TER50, TER75, TER150) trên mô hình phi thủy tĩnh độ phân giải cao WRF-ARW, tích hợp thuật toán ban đầu hóa xoáy động lực học giúp hiệu chỉnh đồng thời cấu trúc 3D của cả trường gió, áp và ẩm.

3. Phát hiện bất ngờ nhất từ dữ liệu thực nghiệm là gì? Trả lời: Phát hiện bất ngờ nhất là trong thí nghiệm triệt tiêu hoàn toàn địa hình (TER2m), cơn bão khi tiến sát bờ biển không hề bị giảm tốc độ di chuyển và cấu trúc trường gió cực đại duy trì đối xứng hoàn hảo quanh tâm mắt bão. Ngược lại, chỉ cần địa hình đạt 50% độ cao thực tế (TER50), sự bất đối xứng trong phân bố ẩm và trường mưa ngưng kết đã biểu hiện rõ rệt, chứng minh độ nhạy phi tuyến tính cực cao của hoàn lưu bão trước độ gồ ghề của địa hình ven biển.

4. Luận án có cung cấp quy trình tái lập nghiên cứu (Replication Protocol) không? Trả lời: Hoàn toàn có. Toàn bộ cấu hình namelist của mô hình WRF-ARW, tọa độ miền tính lồng ghép, các sơ đồ tham số hóa vật lý (Kain-Fritsch, WSM6, YSU), thuật toán ban đầu hóa xoáy và danh mục 228 ca thí nghiệm cho 18 cơn bão giai đoạn 2014-2018 đều được quy chuẩn hóa chi tiết trong luận án, cho phép các nhà nghiên cứu độc lập tái lập và kiểm chứng kết quả mô phỏng.

5. Chương trình nghiên cứu 10 năm tiếp theo được định hình như thế nào? Trả lời: Tác giả định hình lộ trình nghiên cứu mở rộng sang việc xây dựng hệ thống đồng hóa dữ liệu tổ hợp đa biến kết hợp trí tuệ nhân tạo (AI-NWP Hybrid), tích hợp mô hình hoàn lưu đại dương ba chiều có tính đến triều cường - nước dâng do bão và đánh giá tác động của biến đổi khí hậu toàn cầu đến tần suất xuất hiện các cơn bão dị thường tương tác KKL mạnh đổ bộ vào miền Trung Việt Nam.

Kết luận

  1. Luận án đã giải quyết trọn vẹn và tường minh bài toán động lực học về ảnh hưởng của địa hình duyên hải Việt Nam và không khí lạnh đến sự biến đổi cấu trúc ba chiều của bão trên Biển Đông.
  2. Chứng minh khoa học cơ chế sườn núi đón gió gây bất đối xứng nghiêm trọng trong trường gió mực 10m và tập trung cực đại lượng mưa ngưng kết ở góc phần tư phía trước bão, đồng thời làm sáng tỏ nguyên nhân quỹ đạo bão uốn cong lệch nam khi tiếp cận bờ biển miền Trung.
  3. Xác lập hệ thống tiêu chí nhiệt động lực học phân định rõ ràng giữa hai hình thế tác động trực tiếp (suy yếu nhanh do nêm lạnh cắt nguồn ẩm) và gián tiếp (tăng cường gió mạnh rìa bắc do chênh lệch áp suất) của không khí lạnh lên xoáy bão.
  4. Khẳng định hiệu năng vượt trội của phương pháp ban đầu hóa xoáy động lực học trong việc cải thiện trường phân tích ban đầu, giảm thiểu sai số khí áp tâm và dự báo quỹ đạo bão trong mô hình số trị WRF-ARW.
  5. Mở ra hướng tiếp cận mới trong mô phỏng khí tượng thực nghiệm thông qua thiết kế ma trận biến đổi độ cao địa hình (CTL, TER2m, TER50, TER75, TER150).
  6. Cung cấp cơ sở khoa học chuẩn mực và luận cứ thực tiễn vững chắc phục vụ đắc lực cho công tác dự báo nghiệp vụ cảnh báo thiên tai bão, mưa lớn lũ quét tại Việt Nam, đóng góp thiết thực vào kho tàng tri thức khí tượng học nhiệt đới khu vực và quốc tế.