Tổng quan về luận án

Biển Đông và vùng duyên hải Việt Nam nằm trọn trong tâm điểm hoạt động của một trong những ổ bão nhiệt đới phức tạp nhất hành tinh. Thống kê khí hậu học khẳng định: "Việt Nam nằm ở khu vực của ổ bão ở Tây Bắc Thái Bình Dương (30 cơn/năm, chiếm 38% số bão trên toàn cầu)". Sự tương tác đa quy mô giữa các hoàn lưu quy mô lớn (áp cao cận nhiệt đới Tây Thái Bình Dương, gió mùa Tây Nam, rãnh gió mùa) với địa hình phức tạp tạo nên độ bất định phi tuyến tính sâu sắc trong chuyển động và cường độ của các xoáy thuận nhiệt đới (XTND). Trong bối cảnh biến đổi khí hậu toàn cầu, việc mở rộng hạn dự báo quỹ đạo bão từ ngắn hạn (1–3 ngày) lên hạn vừa (5 ngày / 120 giờ) là đòi hỏi cấp bách nhằm giảm thiểu thiệt hại nhân mạng và kinh tế.

                  ┌────────────────────────────────────────┐
                  │ TÍCH HỢP ĐỘNG LỰC & ĐA CƠ CHẾ VẬT LÝ  │
                  │   - Mô hình khu vực RAMS               │
                  │   - 3 Sơ đồ đối lưu: Kuo, KF, KFCT     │
                  └──────────────────┬─────────────────────┘
                                     │
            ┌────────────────────────┴────────────────────────┐
            ▼                                                 ▼
┌───────────────────────┐                         ┌───────────────────────┐
│   NHIỄU MÔI TRƯỜNG   │                         │      NHIỄU XOÁY       │
│ - Chuẩn hóa Gauss     │                         │ - Cài xoáy giả 3D     │
│ - Trực giao hóa       │                         │ - Bất đối xứng động lực│
└───────────┬───────────┘                         └───────────┬───────────┘
            │                                                 │
            └────────────────────────┬────────────────────────┘
                                     ▼
                  ┌─────────────────────────────────────┐
                  │    CHU TRÌNH NUÔI NHIỄU PHÁT TRIỂN  │
                  │       (BGM/BV Cycle - 24 giờ)       │
                  └──────────────────┬──────────────────┘
                                     ▼
                  ┌─────────────────────────────────────┐
                  │   TỔ HỢP 39 THÀNH PHẦN BAN ĐẦU      │
                  └──────────────────┬──────────────────┘
                                     ▼
                  ┌─────────────────────────────────────┐
                  │    SIÊU TỔ HỢP (SUPERENSEMBLE)      │
                  │   - Tuyển chọn 25 thành phần tối ưu │
                  │   - Hiệu chỉnh hồi quy đa biến      │
                  │   - Đánh giá độ tin cậy A, B, C     │
                  └─────────────────────────────────────┘

Khoảng trống nghiên cứu then chốt mà luận án giải quyết xuất phát từ giới hạn căn bản của các trường phân tích toàn cầu và các mô hình tất định:

  • Sai số vị trí và cường độ bão trong các trường phân tích ban đầu còn lớn do thiếu hụt số liệu quan trắc bề mặt biển.
  • Kỹ thuật cài xoáy giả (bogusing) truyền thống dạng đối xứng 3 chiều chỉ định vị lại tâm bão mà hoàn toàn bỏ qua cấu trúc bất đối xứng động lực vốn quyết định hướng dịch chuyển thực tế.
  • Các trung tâm quốc tế hàng đầu như NCEP (Hoa Kỳ), ECMWF (Châu Âu), JMA (Nhật Bản) đã vận hành hệ thống dự báo tổ hợp (EPS) 5 ngày, nhưng tại Việt Nam, các mô hình số trị khu vực nghiệp vụ trước đây chưa từng thiết lập quy trình dự báo tổ hợp quỹ đạo bão hạn 5 ngày có khả năng phân tách và đồng hóa đồng thời nhiễu môi trường và nhiễu xoáy.

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

  • Câu hỏi nghiên cứu 1 (RQ1): Làm thế nào để xây dựng và thích ứng kỹ thuật nuôi nhiễu (Breeding of Growing Modes - BGM) cho trường ban đầu chứa xoáy thuận nhiệt đới trên mô hình khí quyển khu vực RAMS?
  • Câu hỏi nghiên cứu 2 (RQ2): Sự tương tác phi tuyến giữa nhiễu xoáy (vortex perturbation) và nhiễu môi trường (environmental perturbation) chi phối như thế nào đến độ tán tổ hợp (ensemble spread) và sai số quỹ đạo hạn 120 giờ?
  • Câu hỏi nghiên cứu 3 (RQ3): Phương pháp siêu tổ hợp (superensemble) kết hợp hiệu chỉnh thống kê đa biến có khả năng triệt tiêu sai số hệ thống để nâng cao chất lượng dự báo bão hạn 5 ngày tại Biển Đông hay không?
  • Giả thuyết 1 (H1): Chuỗi nhiễu ban đầu kết hợp đồng thời nhiễu môi trường được trực giao hóa và nhiễu xoáy bất đối xứng sẽ duy trì tốc độ phát triển động lực tối ưu, bao trùm trạng thái khí quyển thực tế tốt hơn phương pháp nhiễu ngẫu nhiên hoặc nuôi nhiễu đơn lẻ.
  • Giả thuyết 2 (H2): Việc tích hợp đa cơ chế vật lý (3 sơ đồ tham số hóa đối lưu: Kuo, Kain-Fritsch [KF], và Kain-Fritsch cải tiến [KFCT]) vào tập thành phần nuôi nhiễu sẽ làm tăng tính phân kỳ động lực, mở rộng không gian nghiệm có ích của tổ hợp.
  • Giả thuyết 3 (H3): Tuyển chọn 25 thành phần ưu việt từ không gian 39 thành phần thông qua phương trình hồi quy siêu tổ hợp sẽ giảm thiểu sai số khoảng cách vị trí tâm bão ở hạn 120 giờ xuống mức tương đương hoặc vượt trội so với các sản phẩm dự báo tổ hợp đơn lẻ từ các trung tâm khí tượng quốc tế.

Luận án kế thừa và phát triển nền tảng lý thuyết hỗn loạn khí quyển của Lorenz (1963, 1969), khung lý thuyết vector nuôi (Bred Vectors - BV) của Toth và Kalnay (1993, 1997), cùng nguyên lý dự báo tổ hợp hạn ngắn/hạn vừa của Du (2007). Nghiên cứu khảo sát tập mẫu gồm 199 trường hợp thử nghiệm thực tế đối với các cơn bão hoạt động trên Tây Bắc Thái Bình Dương và Biển Đông trong giai đoạn 5 năm (2009–2013), cung cấp luận cứ thực nghiệm vững chắc cho mô hình số trị khí tượng nhiệt đới.

Literature Review và Positioning

Lý thuyết dự báo thời tiết bằng mô hình số trị (NWP) chịu sự chi phối tuyệt đối của "hiệu ứng cánh bướm" (Butterfly Effect). Lorenz (1963, 1969) đã chứng minh tính hữu hạn của khả năng dự báo khí quyển bắt nguồn từ hai nguồn sai số bất khả kháng: sai số phân tích trường ban đầu (sự chênh lệch giữa trạng thái thực của khí quyển và trường phân tích) và sai số mô hình (do sự xấp xỉ các quá trình vật lý vi mô). Nhận định kinh điển của Toth và Kalnay (1993) khẳng định: "Nếu dự báo cho độ chính xác thấp tại một khu vực nào đó là dấu hiệu của sự bất định của dòng cơ bản. Các sai số nhỏ của trường ban đầu (hoặc không hoàn hảo trong mô hình) sẽ phát triển với biên độ lớn trong thời gian hữu hạn".

                       TIẾN TRÌNH TIỄN HÓA CÁC KỸ THUẬT TẠO TỔ HỢP BAN ĐẦU
                       
   [Phương pháp Cổ điển]                [Trường phái Tối ưu]                [Trường phái Lọc thích nghi]
   - Monte Carlo ngẫu nhiên              - Véctơ kỳ dị (SV)                 - Ensemble Kalman Filter (EnKF)
   - Trễ thời gian (LAF/SLAF)            (ECMWF - Buizza 1994)              (Evensen 1994; Houtekamer 1998)
   (Hoffman & Kalnay 1983)                        │                                       │
            │                                     │                                       │
            └─────────────────────────┐           │           ┌───────────────────────────┘
                                      ▼           ▼           ▼
                                   [Trường phái Động lực học (BGM/BV)]
                                   - Toth & Kalnay (1993, 1997); Kalnay (2003)
                                   - Biến đổi tổ hợp ET (Wei et al. 2007)
                                                  │
                                                  ▼
                               [ĐÓNG GÓP CỦA LUẬN ÁN - CÔNG THANH (2014)]
                               - Tách & ghép: Nhiễu Môi trường + Nhiễu Xoáy
                               - Tích hợp Đa sơ đồ đối lưu (Kuo, KF, KFCT)
                               - Siêu tổ hợp (Superensemble) trên mô hình RAMS

Tổng quan các nhánh tiếp cận tạo trường nhiễu ban đầu trong dự báo tổ hợp toàn cầu và khu vực:

  1. Phương pháp trễ thời gian (LAF & SLAF): Hoffman và Kalnay (1983) đề xuất Lagged-Average Forecasting (LAF) tận dụng các dự báo quá khứ làm nhiễu ban đầu, phản ánh được "sai số ngày" nhưng nhược điểm là chứa các dự báo quá cũ làm giảm độ chính xác. Ebisuzaki và Kalnay (1991) cải tiến thành Scaled Lagged-Average Forecasting (SLAF) bằng cách chuẩn hóa biên độ sai số, song quy mô thành phần vẫn bị giới hạn nghiêm ngặt.
  2. Phương pháp Véctơ kỳ dị (Singular Vectors - SV): Buizza (1994), Palmer và cộng sự (1998), Errico và Baumhefner (1998) ứng dụng SV trên mô hình tiếp tuyến tuyến tính (TLM) và mô hình liên hợp (adjoint) tại ECMWF. SV tối ưu hóa toán học tốc độ tăng trưởng năng lượng tổng trong tương lai (0–48 giờ). Tuy nhiên, phương pháp này đòi hỏi chi phí tính toán khổng lồ (gấp 3 lần số lượng vector kỳ dị cần tìm) và cấu trúc nhiễu cực kỳ nhạy cảm với việc lựa chọn chuẩn tối ưu.
  3. Phương pháp Véctơ nuôi (Breeding of Growing Modes - BGM/BV): Toth và Kalnay (1993, 1997), Kalnay (2003) xây dựng chu trình tích phân phi tuyến phi đối xứng. Véctơ nuôi phản ánh chân thực sai số phân tích động lực học của quá khứ đi vào chu trình đồng hóa số liệu. Wang và Bishop (2003) cùng Martin và cộng sự (2007) phản biện rằng các thành phần BV truyền thống có xu hướng suy giảm tính trực giao, dẫn đến tương quan cao giữa các thành phần và làm giảm dung lượng thông tin tổ hợp.
  4. Phương pháp Lọc Kalman tổ hợp (EnKF & ETKF): Evensen (1994), Houtekamer và Mitchell (1998), Bishop (2001), Whitaker và Hamill (2002), Wang và Bishop (2003) phát triển Ensemble Transform Kalman Filter (ETKF) và Local Ensemble Kalman Filter (LEKF). Hamill và cộng sự (2000b) cùng Bowler (2006) và Descamps và Talagrand (2007) so sánh trên các mô hình Lorenz và khẳng định EnKF/ETKF cung cấp thống kê sai số nền tối ưu nhưng đòi hỏi hạ tầng siêu máy tính phức tạp.

Đối với bài toán xoáy nhiệt đới, các trung tâm quốc tế đã phát triển các hệ thống chuyên biệt:

  • JMA (Nhật Bản): Vận hành hệ thống Typhoon Ensemble Prediction System (TEPS) và One-Week EPS (WEPS), ứng dụng kỹ thuật đồng hóa xoáy giả 3 chiều kết hợp SV (Yamaguchi, 2012).
  • CMA / Viện Nghiên cứu Khí tượng Thượng Hải (Trung Quốc): Xây dựng hệ thống GRAPES-TCM sử dụng BGM và phương pháp cực đại hóa nhiễu phi tuyến (CNOP) (Chen và nnk, 2009; Huang và nnk, 2007).
  • NCEP (Hoa Kỳ): Triển khai Short-Range Ensemble Forecast (SREF) kết hợp biến đổi tổ hợp ET và đa mô hình động lực (Du, 2007, 2011).

Luận án định vị trọng tâm vào khoảng trống: Thay vì chỉ áp dụng BGM cổ điển lên trường quy mô lớn như NCEP hay cấy xoáy nhân tạo tĩnh như JMA, nghiên cứu đề xuất giải pháp tách biệt không gian trạng thái thành nhiễu môi trường và nhiễu xoáy động lực, sau đó đồng thời chuẩn hóa, trực giao hóa và tích hợp vào mô hình RAMS với 3 sơ đồ tham số hóa đối lưu độc lập cho khu vực Biển Đông.

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

┌────────────────────────────────────────────────────────────────────────────────────────┐
│                   KHUNG PHÂN TÍCH TÍCH HỢP TỔ HỢP ĐA YẾU TỐ VÀ ĐỘNG LỰC                 │
├──────────────────────────────────────────┬─────────────────────────────────────────────┤
│ TRỤ CỘT ĐỘNG LỰC KHÍ QUYỂN               │ TRỤ CỘT THỐNG KÊ VÀ HẬU XỬ LÝ               │
├──────────────────────────────────────────┼─────────────────────────────────────────────┤
│ 1. Lý thuyết BGM/BV phi tuyến            │ 1. Trực giao hóa Gram-Schmidt / Ma trận ET │
│    - Toth & Kalnay (1993, 1997)          │    - Triệt tiêu tương quan thành phần       │
│ 2. Lý thuyết Nhiễu động lực học xoáy     │ 2. Hiệu chuẩn sai số phân tích thực         │
│    - Phân tách trường sóng và hoàn lưu  │    - Phân bố chuẩn Gauss: 3 m/s & 1°C       │
│ 3. Đa cơ chế vật lý mây đối lưu          │ 3. Siêu tổ hợp thống kê đa biến (Superens)  │
│    - Sơ đồ Kuo, KF, KFCT (Nguyễn M. Trung│    - Tuyển lựa 25 thành phần từ 39 nghiệm   │
└──────────────────────────────────────────┴─────────────────────────────────────────────┘

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

Luận án mở rộng khung lý thuyết dự báo tổ hợp của Du (2007) đối với các hệ thống thời tiết nguy hiểm quy mô vừa trong điều kiện nhiệt đới gió mùa:

"Dự báo tổ hợp là một phương pháp động lực. Dòng cơ bản phụ thuộc vào định lượng của độ bất định dự báo (sai số ngày) và cung cấp thông tin cơ sở để liên kết giữa độ bất định dự báo và độ tin cậy dự báo cho người dùng để sau đó có thể được chuẩn bị tốt nhất... Ba đặc tính cơ bản cần phải tuân theo khi nghiên cứu những dao động của nhiễu phát triển nhanh trong trường ban đầu: Tính sát thực, tính phân kỳ và tính trực giao. Trong đó tính sát thực là độ lớn của nhiễu cần nằm