Tổng quan về luận án

Sự biến động của hoàn lưu khí quyển và tương tác khí quyển - đại dương trong bối cảnh biến đổi khí hậu toàn cầu đang làm thay đổi sâu sắc chế độ mưa tại khu vực nhiệt đới gió mùa Đông Nam Á. Luận án tiến sĩ chuyên ngành Khí tượng học (Mã số: 9440222.01) của tác giả Phạm Thanh Hà, dưới sự hướng dẫn khoa học của GS. Phan Văn Tân tại Trường Đại học Khoa học Tự nhiên – Đại học Quốc gia Hà Nội (2022), mang tiêu đề: "Dự báo hạn mùa và nội mùa ngày bắt đầu mùa mưa ở Việt Nam trên cơ sở sản phẩm mô hình số". Công trình đã giải quyết bài toán cốt lõi trong khí tượng học nhiệt đới: xác lập cơ sở khoa học khách quan để định lượng ngày bắt đầu mùa mưa (NBĐMM) và phát triển công nghệ dự báo hạn nội mùa (7–40 ngày) và hạn mùa (dưới 6 tháng) dựa trên sự kết hợp giữa mô hình số trị và các phương pháp thống kê tiên tiến.

Trong nông nghiệp nhiệt đới, thời điểm bắt đầu mùa mưa chi phối trực tiếp đến mùa vụ canh tác lúa nước, an ninh lương thực và quy trình tích trữ, điều tiết nước của các công trình thủy lợi, thủy điện. Như nguồn tài liệu luận án đã trích dẫn: "Đối với lĩnh vực nông nghiệp, lượng mưa là một trong những nhân tố quan trọng ảnh hưởng trực tiếp đến năng suất và thu nhập của người nông dân [Lambert, 2014]. Tuy nhiên, nền nông nghiệp ở Việt Nam dường như dễ bị tổn thương bởi những biến đổi bất thường của thời tiết nói chung, mưa nói riêng". Mặc dù giữ vai trò sống còn, NBĐMM không phải là một đại lượng vật lý quan trắc trực tiếp được ghi nhận bằng khí cụ đơn lẻ mà phải được xác định thông qua các hệ thống chỉ tiêu gián tiếp.

Khoảng trống nghiên cứu (research gap) lớn nhất trong các công trình trước đây tại Việt Nam (Nguyễn Hướng Điền và cs, 2000; Nguyễn Thị Hiền Thuận, 2007; Lê Thị Xuân Lan và cs, 2017; Ngô Thị Thanh Hương, 2018) nằm ở tính cục bộ và chủ quan khi lựa chọn ngưỡng lượng mưa và số ngày mưa cố định cho từng trạm riêng lẻ hoặc chỉ tập trung vào khu vực Nam Bộ và Tây Nguyên mà chưa có một phương pháp luận toán học thống nhất cho toàn lãnh thổ. Đồng thời, việc ứng dụng các mô hình động lực toàn cầu thế hệ mới (như mô hình CFSv2 của NCEP và hệ thống S2S của ECMWF) vào bài toán dự báo hạn mùa và nội mùa NBĐMM tại Việt Nam chưa được khai thác một cách hệ thống.

Luận án đặt ra 3 câu hỏi nghiên cứu (Research Questions - RQ) và 3 giả thuyết khoa học (Hypotheses - H) tương ứng:

  • RQ1: Làm thế nào để xây dựng một bộ chỉ tiêu khách quan xác định NBĐMM vừa phản ánh được tính đồng nhất quy mô vùng khí hậu, vừa bảo toàn đặc trưng vi khí hậu cục bộ tại từng trạm?
    • H1: NBĐMM chịu sự chi phối đồng thời của hoàn lưu gió mùa quy mô lớn và địa hình cục bộ, do đó tồn tại một không gian tham số tối ưu hóa được mối tương quan giữa tín hiệu trạm và tín hiệu vùng.
  • RQ2: Mối quan hệ tương tác giữa dao động khí quyển - đại dương quy mô lớn (ENSO) với sự biến động liên năm của NBĐMM trên từng vùng khí hậu Việt Nam diễn ra như thế nào?
    • H2: Các pha nóng (El Niño) và pha lạnh (La Niña) tạo ra sự lệch pha có ý nghĩa thống kê đối với thời điểm bắt đầu mùa mưa trên các vùng khí hậu khác nhau.
  • RQ3: Khả năng dự báo hạn mùa và hạn nội mùa NBĐMM của các mô hình số toàn cầu (CFSv2, S2S-ECMWF) đạt mức độ chính xác và kỹ năng dự báo ($SS_{MAE}$) ra sao khi tiếp cận theo phương pháp trực tiếp và phương pháp thống kê - động lực?
    • H3: Phương pháp thống kê - động lực khai thác các trường hoàn lưu quy mô lớn của mô hình số sẽ cho kỹ năng dự báo vượt trội so với việc sử dụng trực tiếp trường mưa thô của mô hình.

Khung lý thuyết của luận án được xây dựng dựa trên sự tích hợp của Lý thuyết Động lực học Gió mùa Châu Á (Asian Summer Monsoon Dynamics), Lý thuyết Hoàn lưu Khí quyển - Đại dương (Coupled Atmosphere-Ocean Teleconnections) và Lý thuyết Dự báo Hạn nội mùa đến Mùa (S2S Prediction Framework). Về quy mô không gian và thời gian, nghiên cứu sử dụng chuỗi số liệu quan trắc 41 năm (1979–2019) tại 131 trạm khí tượng bề mặt phân bố trên 7 vùng khí hậu Việt Nam (R1: Tây Bắc, R2: Đông Bắc, R3: Đồng bằng Bắc Bộ, R4: Bắc Trung Bộ, R5: Nam Trung Bộ, R6: Tây Nguyên, R7: Nam Bộ), kết hợp với số liệu tái phân tích và dự báo lại từ CFSv2 (1983–2019, độ phân giải $1^\circ \times 1^\circ$) và S2S-ECMWF (2000–2019, độ phân giải $0.5^\circ \times 0.5^\circ$, 11 thành phần tổ hợp).

Literature Review và Positioning

Tổng quan y văn quốc tế và trong nước về xác định và dự báo NBĐMM phân hóa thành hai trường phái tiếp cận chính:

Trường phái thứ nhất dựa trên đặc trưng chuỗi số liệu mưa bề mặt tại quy mô địa phương (Stern et al., 1981; Sivakumar, 1988; Omotosho, 1992; Moron et al., 2009). Nhóm phương pháp này định nghĩa NBĐMM dựa trên tổng lượng mưa tích lũy $P$ trong chuỗi $N$ ngày liên tiếp và kiểm soát nguy cơ hạn giả bằng điều kiện không xuất hiện chuỗi ngày khô $D$ ngày tiếp theo. Nhằm khắc phục tính chủ quan của các ngưỡng cố định, Liebmann và Marengo (2001) cùng Bombardi et al. (2017) đã phát triển phương pháp giải tích đường chuẩn sai tích lũy lượng mưa ngày (anomalous accumulation). Tuy nhiên, phương pháp chuẩn sai tích lũy gặp hạn chế lớn trong các năm dị thường cực đoan (quá ẩm hoặc quá khô), dẫn đến việc nhận định ngày bắt đầu bị lệch pha nghiêm trọng so với thực tế quan trắc.

Trường phái thứ hai dựa trên sự chuyển biến của hoàn lưu khí quyển quy mô lớn và đặc trưng nhiệt - động lực học (Webster & Yang, 1992; Li & Yanai, 1996; Tong et al., 2009; Nguyen et al., 2014). Hướng tiếp cận này sử dụng các biến trường như gió vĩ hướng tầng 850 hPa (U850), gió kinh hướng tầng 200 hPa, bức xạ sóng dài đi ra (OLR) và gradient nhiệt độ tầng đối lưu để xác định Ngày bắt đầu gió mùa (NBĐGM).

Cuộc tranh luận học thuật cốt lõi nảy sinh từ sự không đồng nhất giữa NBĐMM và NBĐGM. Các nghiên cứu quốc tế tại Tây Phi (Omotosho et al., 2000) và nghiên cứu tại Việt Nam của Nguyễn Hướng Điền và cs (2000), Ngo-Thanh et al. (2018) đã chứng minh rằng NBĐMM và NBĐGM không hoàn toàn trùng khớp. Tại Tây Nguyên và Nam Bộ, NBĐMM quy mô địa phương thường đến sớm hơn NBĐGM mùa hè quy mô khu vực từ 2 đến 3 tuần do hoạt động đối lưu nhiệt địa phương hình thành trước khi hệ thống gió mùa mùa hè Tây Nam thiết lập ổn định.

So sánh với các nghiên cứu quốc tế điển hình:

  1. Khu vực Sahel / Tây Phi: Nghiên cứu của Vellinga et al. (2013) sử dụng mô hình GloSea4 của Met Office chỉ ra khả năng dự báo NBĐGM trước 2–3 tháng dựa trên trường nhiệt độ mặt biển (SST) nhiệt đới với diện tích dưới đường cong ROC đạt trên 0.7. Kumi et al. (2020) khẳng định các mô hình động lực có thể dự báo các pha mùa mưa Tây Phi ở hạn 20–60 ngày, nhưng chất lượng phụ thuộc sống còn vào chỉ tiêu định nghĩa.
  2. Khu vực Gió mùa Ấn Độ: Alessandri et al. (2015) khai thác mô hình CMCC và Chevuturi et al. (2018) sử dụng mô hình GloSea5-GC2 chứng minh rằng mô hình động lực mô phỏng các chỉ tiêu hoàn lưu quy mô lớn (gió U850, gradient nhiệt đối lưu) tốt hơn nhiều so với việc mô phỏng lượng mưa đối lưu điểm lưới địa phương.

Luận án của Phạm Thanh Hà định vị chính xác vào điểm giao thoa này: hóa giải khoảng cách giữa quy mô khu vực và quy mô địa phương bằng phương pháp phân tích thành phần chính nhiều biến (MPCA) kết hợp kiểm định thống kê Fisher, thiết lập khung dự báo hạn mùa và nội mùa chuyên biệt cho điều kiện địa hình phức tạp kéo dài trên 15 vĩ độ của Việt Nam.

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

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

Luận án đóng góp vào Lý thuyết Khí hậu học Động lực thông qua việc làm sáng tỏ cơ chế ghép nối giữa dao động hoàn lưu quy mô hành tinh với phản ứng thời tiết địa phương. Nghiên cứu đã hệ thống hóa và chứng minh bằng thực nghiệm rằng tác động của hiện tượng ENSO lên NBĐMM có sự phân hóa không gian sâu sắc:

  • Trong các năm El Niño (pha nóng), NBĐMM có xu thế xuất hiện muộn hơn đáng kể so với trung bình nhiều năm tại khu vực Tây Nguyên (R6) và Nam Bộ (R7).
  • Ngược lại, trong các năm La Niña (pha lạnh), gió mùa mùa hè và mùa mưa có xu thế bắt đầu sớm hơn từ 1 đến 2 tuần.
       +-----------------------------------------------------------+
       |   CƠ CHẾ ĐIỀU BIẾN QUY MÔ LỚN (ENSO / SST / PRMSL / U850) |
       +-----------------------------------------------------------+
                                     |
                                     v
       +-----------------------------------------------------------+
       |         TÍN HIỆU KHÍ HẬU VÙNG (MPCA -> EOF1 Mode)         |
       +-----------------------------------------------------------+
                                     |
                         [Tối ưu hóa tương quan]
                                     |
                                     v
       +-----------------------------------------------------------+
       |     BỘ CHỈ TIÊU ĐỊA PHƯƠNG TẠI TRẠM (P_opt, N_opt, D_opt) |
       |              (Kiểm định Fisher ANOVA f-test)              |
       +-----------------------------------------------------------+
                                     |
                                     v
       +-----------------------------------------------------------+
       |    NGÀY BẮT ĐẦU MÙA MƯA QUAN TRẮC CHUẨN XÁC (1979-2019)    |
       +-----------------------------------------------------------+

Khung lý thuyết thiết lập 4 mệnh đề khoa học (Propositions):

  • Mệnh đề 1 ($P_1$): Tín hiệu bắt đầu mùa mưa tại một điểm trạm là hàm phụ thuộc của hoàn lưu gió mùa quy mô vùng và hiệu ứng điều biến địa hình cục bộ.
  • Mệnh đề 2 ($P_2$): Tồn tại một ma trận 42 chỉ tiêu kết hợp ($P \times N \times D$) cho phép quét toàn diện không gian trạng thái mưa chuyển mùa, trong đó cực đại của hàm tương quan giữa chuỗi trạm và vector riêng EOF1 đại diện cho bộ chỉ tiêu tối ưu.
  • Mệnh đề 3 ($P_3$): Trường mưa dự báo thô của mô hình động lực bị suy giảm kỹ năng do sai số hệ thống trong tham số hóa đối lưu, trong khi các trường động lực quy mô lớn (khí áp mực biển PRMSL, gió vĩ hướng U850, nhiệt độ mặt biển SST) duy trì tín hiệu dự báo ổn định ở hạn mùa.
  • Mệnh đề 4 ($P_4$): Kỹ năng dự báo hạn nội mùa của sản phẩm mô hình số tỷ lệ nghịch với hạn dự báo và phụ thuộc trực tiếp vào việc lựa chọn chỉ tiêu xác định NBĐMM tương thích với độ phân giải mô hình.

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

Khung phân tích của luận án tích hợp 3 trụ cột phương pháp luận:

  1. Phương pháp Tối ưu hóa Chỉ tiêu Khí hậu: Kết hợp độ nhạy $SI$ (Sensitivity Index) với phân tích thành phần chính nhiều biến (MPCA) để trích xuất mode biến động vượt trội (EOF1 giải thích tỷ lệ phương sai lớn nhất) đại diện cho tín hiệu vùng.
  2. Kỹ thuật Ghép nối Thống kê - Động lực (Statistical-Dynamical Downscaling): Ứng dụng hồi quy tuyến tính từng bước (Stepwise Multiple Linear Regression) để sàng lọc các nhân tố dự báo tối ưu từ lưới dữ liệu CFSv2 toàn cầu ($40^\circ\text{E} - 100^\circ\text{W}, 40^\circ\text{S} - 40^\circ\text{N}$).
  3. Quy trình Kiểm chứng Nghiêm ngặt (Cross-Validation Protocol): Ứng dụng phương pháp kiểm chứng chéo loại một (One-Leave-One-Out - OLOO) nhằm triệt tiêu hiện tượng quá khớp (overfitting) và đảm bảo tính độc lập tuyệt đối của bộ dữ liệu thử nghiệm.

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

Thiết kế nghiên cứu

Nghiên cứu được thiết kế theo trường phái Thực chứng luận (Positivism), dựa trên phân tích định lượng chuỗi thời gian khí tượng thực nghiệm và mô phỏng số trị động lực học khí quyển. Quy trình nghiên cứu đa cấp độ (multi-level design) được triển khai qua 3 tầng phân tích: Cấp điểm trạm (131 trạm) $\rightarrow$ Cấp vùng khí hậu (7 vùng R1–R7) $\rightarrow$ Cấp hoàn lưu khu vực và toàn cầu.

Tiêu chí lựa chọn mẫu dữ liệu quan trắc đảm bảo tính đại diện nghiêm ngặt theo chuẩn WMO:

  • Tổng số ngày khuyết trong một năm không vượt quá 25 ngày ($\ge 340$ ngày quan trắc/năm).
  • Trong giai đoạn chuyển mùa trọng điểm (tháng 3 đến tháng 6), số ngày thiếu liên tiếp không quá 3 ngày. Những ngày thiếu ngẫu nhiên thỏa mãn điều kiện được bổ khuyết bằng giá trị trung bình nhiều năm của chính ngày đó.

Quy trình nghiên cứu rigorous

Quy trình xác định chỉ tiêu NBĐMM tối ưu được chuẩn hóa qua 3 bước toán học:

Bước 1: Quét không gian tham số và tính chỉ số độ nhạy ($SI$) Xây dựng 42 bộ chỉ tiêu kết hợp từ tập tham số thực nghiệm:

  • Ngưỡng lượng mưa 5 ngày tích lũy: $P \in {10, 15, 20, 25, 30, 40, 50}\text{ mm}$.
  • Số ngày có mưa trong 5 ngày: $N \in {2, 3, 4}\text{ ngày}$ (với mưa ngày $> 0.1\text{ mm}$).
  • Điều kiện kiểm soát hạn giả: Không xuất hiện chuỗi $\ge 7$ ngày khô liên tiếp trong $D \in {20, 30}\text{ ngày}$ sau NBĐMM.

Độ nhạy của NBĐMM đối với các chỉ tiêu tại từng trạm được lượng hóa qua chỉ số $SI$: $$SI = \frac{1}{L \times M \times (M-1)} \sum_{k=1}^{L} \sum_{j=1}^{M} \sum_{e=j+1}^{M} |X_{kj} - X_{ke}|$$ Trong đó $L = 41$ năm (1979–2019), $M = 42$ chỉ tiêu, $X_{kj}$ và $X_{ke}$ là NBĐMM xác định theo chỉ tiêu $j$ và $e$.

Bước 2: Trích xuất tín hiệu vùng bằng MPCA Ma trận dữ liệu đầu vào có kích thước $([42 \times N_{sta}] \times [42 \times N_{sta}])$ được phân tích thành phần chính để tìm các vector riêng EOFs và chuỗi hệ số thời gian PCs tương ứng. Mode EOF1 giải thích phần lớn tổng phương sai biến động liên năm của toàn vùng.

Bước 3: Tối ưu hóa tương quan và kiểm định ANOVA f-test Chỉ tiêu tối ưu được chọn cho từng trạm là bộ $(P, N, D)$ có hệ số tương quan cao nhất giữa chuỗi NBĐMM trạm và chuỗi PC1 của vùng. Đồng thời, phân tích phương sai một chiều (One-Way ANOVA, kiểm nghiệm Fisher) với mức ý nghĩa $\alpha = 0.05$ được thực hiện để đánh giá độ tin cậy thống kê của việc thay đổi các ngưỡng: $$F = \frac{MSB}{MSE} = \frac{\frac{1}{K-1}\sum_{i=1}^{K} n_i (\bar{Y}i - \bar{Y})^2}{\frac{1}{N-K}\sum{i=1}^{K}\sum_{j=1}^{n_i} (Y_{ij} - \bar{Y}_i)^2}$$

Xác định xu thế biến đổi dài hạn bằng ước lượng độ dốc phi tham số Sen (Sen's Slope Estimator) và kiểm định phi tham số Mann-Kendall ở mức tin cậy 95%.

Data và phân tích

+---------------------------------------------------------------------------------------+
| NGUỒN SỐ LIỆU QUAN TRẮC: 131 Trạm Khí tượng (1979-2019) / 7 Vùng Khí hậu (R1 - R7)    |
+---------------------------------------------------------------------------------------+
        |                                                       |
        v                                                       v
+------------------------------------+  +-----------------------------------------------+
| DỰ BÁO HẠN MÙA (< 6 Tháng)         |  | DỰ BÁO HẠN NỘI MÙA (7 - 40 Ngày)              |
| - CFSv2 Reforecast (1983-2019)     |  | - S2S-ECMWF Reforecast (2000-2019)            |
| - Biến: PRMSL, U850, SST, Rain     |  | - 11 Ensemble Members, CY46R1/CY47R1          |
| - Lưới: 1.0° x 1.0°                |  | - Lưới: 0.5° x 0.5°, Tco639/Tco319            |
+------------------------------------+  +-----------------------------------------------+
        |                                                       |
        +---------------------------+---------------------------+
                                    |
                                    v
+---------------------------------------------------------------------------------------+
| HỆ THỐNG ĐÁNH GIÁ CHẤT LƯỢNG NGHIỆM THU:                                              |
| - Kiểm chứng chéo One-Leave-One-Out (OLOO)                                            |
| - Bộ chỉ số kỹ năng: CORR (Pearson), ME (Sai số trung bình), MAE (Tuyệt đối TB)        |
| - Điểm kỹ năng chuẩn hóa: SS_MAE = 1 - (MAE_forecast / MAE_climatology)               |
+---------------------------------------------------------------------------------------+

Bộ chỉ số đánh giá kỹ năng dự báo bao gồm Hệ số tương quan Pearson ($CORR$), Sai số trung bình ($ME$), Sai số quân phương ($RMSE$), Sai số tuyệt đối trung bình ($MAE$) và Điểm kỹ năng sai số tuyệt đối ($SS_{MAE}$): $$SS_{MAE} = 1 - \frac{MAE_{mô hình}}{MAE_{khí hậu}}$$ Dự báo được xác nhận có kỹ năng thực chất khi $SS_{MAE} > 0$ (vượt trội hơn phương pháp dự báo dựa trên giá trị trung bình khí hậu nhiều năm).

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

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

  1. Quy luật phân bố không gian của NBĐMM và tính bất định cục bộ: NBĐMM xuất hiện sớm nhất tại khu vực Tây Bắc (R1) và Tây Nguyên (R6) vào khoảng cuối tháng 4 (trung bình từ 20/4 đến 5/5). Nam Bộ (R7) và Đồng bằng Bắc Bộ (R3) bắt đầu mùa mưa vào đầu tháng 5. Riêng khu vực Nam Trung Bộ (R5) và dải duyên hải Bắc Trung Bộ (R4), mùa mưa dịch chuyển đặc trưng sang mùa thu - đông (bắt đầu từ tháng 8 đến tháng 9). Chỉ số độ nhạy $SI$ tại các vùng địa hình chia cắt phức tạp (dãy Hoàng Liên Sơn ở R1, vùng núi Đông Bắc R2, sườn đông Trường Sơn R4) đạt giá trị rất cao ($SI > 15\text{ ngày}$), khẳng định việc áp đặt một chỉ tiêu duy nhất cho toàn quốc là sai lầm khoa học nghiêm trọng.

  2. Tác động điều biến có ý nghĩa thống kê của ENSO: Phân tích tương quan chuẩn sai giữa chỉ số $ONI$ trung bình các tháng 12–1–2 với NBĐMM trên 7 vùng khí hậu chứng minh: Trong các pha El Niño mạnh (như 1982–1983, 1997–1998, 2015–2016), xác suất NBĐMM đến muộn hơn trung bình nhiều năm tại R6 và R7 lên tới trên 75% với hệ số tương quan đạt mức tin cậy 95% ($p < 0.05$). Ngược lại, pha La Niña thúc đẩy đối lưu sớm, đưa NBĐMM đến sớm hơn từ 7 đến 15 ngày.

  3. Sự vượt trội của phương pháp thống kê - động lực trong dự báo hạn mùa (CFSv2): Khi sử dụng trực tiếp sản phẩm mưa thô từ CFSv2 ở các hạn từ 1 đến 6 tháng, điểm kỹ năng $SS_{MAE}$ tại hầu hết các vùng khí hậu đều nhận giá trị âm ($SS_{MAE} < 0$), hệ số tương quan thấp ($CORR < 0.2$). Tuy nhiên, khi chuyển sang phương pháp thống kê - động lực (sàng lọc các trường PRMSL, SST, U850 từ CFSv2 làm nhân tố đầu vào cho phương trình hồi quy từng bước), chất lượng dự báo tăng vọt: hệ số tương quan kiểm chứng độc lập (OLOO giai đoạn 1983–2011) đạt từ 0.45 đến 0.72, điểm kỹ năng $SS_{MAE} > 0$ đạt giá trị dương trên các vùng R3, R6, R7 ở các hạn 1, 2 và 3 tháng.

  4. Kỹ năng dự báo hạn nội mùa của hệ thống S2S-ECMWF: Hệ thống S2S-ECMWF (11 ensemble members) chứng minh khả năng dự báo hạn nội mùa (7–40 ngày) rất khả quan đối với sản phẩm mưa. Ở các hạn dự báo ngắn 7–14 ngày, mô hình đạt điểm kỹ năng cao ($SS_{MAE} > 0.35$, $MAE \approx 4\text{–}7\text{ ngày}$). Kỹ năng duy trì giá trị dương đến hạn 21–28 ngày đối với các vùng khí hậu có mùa mưa gió mùa điển hình (R1, R6, R7), mở ra triển vọng ứng dụng nghiệp vụ thực tế.

Implications đa chiều

  • Về mặt Lý thuyết: Cung cấp bằng chứng thực nghiệm giải thích cơ chế tương tác đa quy mô giữa dao động nội mùa nhiệt đới (MJO/BSISO), hệ thống hoàn lưu gió mùa Nam Á - Đông Á và địa hình chắn gió cục bộ tại bán đảo Đông Dương.
  • Về mặt Phương pháp luận: Thiết lập quy trình chuẩn hóa kết hợp MPCA - ANOVA - OLOO có khả năng chuyển giao áp dụng cho các bài toán khí hậu tương tự tại các quốc gia trong khu vực gió mùa Châu Á.
  • Về Ứng dụng Thực tiễn: Cung cấp công cụ tính toán tự động NBĐMM cho mạng lưới khí tượng bề mặt quốc gia; hỗ trợ Tập đoàn Điện lực Việt Nam (EVN) và các ban quản lý lưu vực sông (Sông Đà, Sê San, Srepok, Đồng Nai) tối ưu hóa kế hoạch trữ nước cuối mùa khô và xả lũ đầu mùa mưa.
  • Về Hoạch định Chính sách: Cung cấp cơ sở khoa học định lượng để Bộ Nông nghiệp và Phát triển Nông thôn điều chỉnh lịch xuống giống mùa vụ lúa Hè Thu tại Đồng bằng sông Cửu Long và vụ Mùa tại miền Bắc, chủ động phòng tránh rủi ro hạn đầu mùa và xâm nhập mặn.

Limitations và Future Research

Luận án thẳng thắn chỉ ra 4 giới hạn nghiên cứu (boundary conditions):

  1. Độ phân giải mô hình CFSv2 ($1^\circ \times 1^\circ$): Chưa mô phỏng chi tiết tương tác địa hình quy mô vi khí hậu tại các thung lũng hẹp và sườn đón gió phức tạp ở miền Trung.
  2. Chuỗi số liệu tái dự báo S2S-ECMWF: Thời đoạn 20 năm (2000–2019) tuy đáp ứng chuẩn tối thiểu nhưng còn tương đối ngắn để đánh giá đầy đủ các chu kỳ khí hậu đa thập kỷ như Dao động Thập kỷ Thái Bình Dương (PDO).
  3. Hiện tượng gián đoạn mưa (Dry Spells): Bộ chỉ tiêu kiểm soát hạn giả $D = 20\text{–}30\text{ ngày}$ đã loại bỏ hiệu quả hiện tượng "mùa mưa giả", nhưng chưa phân tích sâu đặc trưng cường độ và độ dài của các đợt hạn nội mùa xảy ra sau thời điểm bắt đầu.
  4. Hiệu chỉnh sai số mưa đối lưu: Chưa áp dụng các thuật toán hiệu chỉnh phi tuyến tính tiên tiến (Machine Learning downscaling) để khắc phục triệt để sai số hệ thống của trường mưa thô mô hình số.

Chương trình nghiên cứu tương lai (Future Research Agenda):

  • Phát triển mô hình tổ hợp đa mô hình (Multi-Model Ensemble - MME) kết hợp các hệ thống dự báo hạn nội mùa từ ECMWF, NCEP, JMA, UKMO trong khuôn khổ dự án quốc tế S2S.
  • Tích hợp mạng nơ-ron học sâu (Deep Learning / Physics-Informed Neural Networks - PINNs) vào bài toán hiệu chỉnh sai số và hạ quy mô không gian.
  • Mở rộng nghiên cứu đồng thời 3 đại lượng liên hoàn: Ngày bắt đầu mùa mưa (Onset), Ngày kết thúc mùa mưa (Cessation) và Độ dài mùa mưa (Length of Rainy Season - LRS).

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

  • Tác động Học thuật: Mở ra hướng nghiên cứu mới về dự báo hạn nội mùa (S2S) tại Việt Nam – lĩnh vực vốn tồn tại khoảng trống giữa dự báo thời tiết hạn vừa ($\le 10$ ngày) và dự báo khí hậu hạn mùa. Dự kiến đóng góp các công bố trên các tạp chí chuyên ngành quốc tế uy tín (Q1/Q2) thuộc hệ thống ISI/Scopus (như Journal of Climate, International Journal of Climatology, Meteorology and Atmospheric Physics).
  • Tác động Ngành Năng lượng và Nông nghiệp: Tối ưu hóa vận hành kinh tế của hệ thống thủy điện quốc gia, tiết kiệm hàng trăm tỷ đồng thông qua việc giảm thiểu xả tràn lãng phí hoặc thiếu hụt nước phát điện vào thời kỳ chuyển tiếp mùa khô - mưa.
  • Tác động Quản lý Rủi ro Thiên tai: Giúp Trung tâm Dự báo Khí tượng Thủy văn Quốc gia nâng cấp các bản tin dự báo khí hậu mùa vụ từ định tính sang định lượng xác suất, giảm nhẹ thiệt hại do hạn hán, thiếu nước đầu mùa tại Tây Nguyên và Nam Bộ.

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

  • Nghiên cứu sinh và Giới nghiên cứu Khí hậu: Tiếp cận phương pháp luận MPCA tối ưu hóa chỉ tiêu khí hậu và quy trình khai thác dữ liệu lớn S2S.
  • Giảng viên và Chuyên gia Khí tượng Thủy văn: Nguồn tài liệu tham khảo chuẩn mực trong đào tạo sau đại học về mô hình hóa khí hậu và dự báo thống kê - động lực.
  • Kỹ sư R&D Nông nghiệp và Năng lượng: Thuật toán dự báo thời điểm chuyển mùa phục vụ xây dựng phần mềm quản lý nông nghiệp thông minh và điều độ hồ chứa nước.
  • Cơ quan Quản lý Nhà nước (Bộ TN&MT, Bộ NN&PTNT): Luận cứ khoa học vững chắc để xây dựng kịch bản ứng phó biến đổi khí hậu và chính sách an ninh nguồn nước 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ì? Đóng góp độc đáo nhất là việc xây dựng thành công phương pháp luận khách quan xác định NBĐMM thông qua việc tối ưu hóa mối quan hệ giữa tín hiệu địa phương (từng trạm) và tín hiệu vùng (trích xuất qua mode EOF1 của phân tích MPCA từ ma trận 42 chỉ tiêu kết hợp), hóa giải mâu thuẫn tồn tại hàng thập kỷ giữa các chỉ tiêu kinh nghiệm chủ quan đơn lẻ.

2. Điểm đổi mới phương pháp luận so với các nghiên cứu quốc tế trước đây? So với Cherchi & Navarra (2003) hay Bombardi et al. (2017), luận án không áp đặt một chỉ tiêu cứng nhắc trên toàn lưới dữ liệu mà kết hợp kiểm định phương sai một chiều ANOVA f-test để nhận diện vai trò nhạy cảm của từng tham số ($P, N, D$) tại từng trạm, đồng thời áp dụng phương pháp thống kê - động lực đa trường (PRMSL, U850, SST) từ mô hình CFSv2 với quy trình kiểm chứng chéo OLOO đa chu kỳ.

3. Phát hiện bất ngờ nhất từ số liệu thực nghiệm là gì? Trường mưa dự báo thô của mô hình hiện đại CFSv2 gần như không có kỹ năng dự báo ($SS_{MAE} < 0$) ở hạn 1–6 tháng do sai số tham số hóa mây đối lưu nhiệt đới, nhưng các trường hoàn lưu động lực quy mô lớn (khí áp mực biển và gió U850) của cùng mô hình đó lại chứa đựng tín hiệu dự báo hạn mùa cực kỳ chuẩn xác khi được xử lý qua mô hình thống kê hồi quy từng bước.

4. Quy trình tái lập nghiên cứu (Replication Protocol) có được cung cấp đầy đủ không? Quy trình được cung cấp chi tiết tuyệt đối: bao gồm công thức toán học tính chỉ số $SI$, cấu trúc ma trận đầu vào MPCA $([42 \times N_{sta}] \times [42 \times N_{sta}])$, phương trình kiểm định Fisher, bộ tọa độ và chỉ số lựa chọn nhân tố dự báo cho từng vùng (R3, R6, R7), cùng hệ thống liên kết truy cập dữ liệu mở CFSv2 và S2S-ECMWF.

5. Lộ trình nghiên cứu 10 năm (10-Year Research Agenda) được định hình ra sao? Xây dựng hệ thống dự báo tổ hợp liền vết (Seamless Ensemble Prediction System) tích hợp Trí tuệ Nhân tạo (AI/Machine Learning) dự báo chuỗi hiện tượng: Bắt đầu mùa mưa $\rightarrow$ Hạn gián đoạn nội mùa $\rightarrow$ Đỉnh mưa cực đoan $\rightarrow$ Kết thúc mùa mưa trên toàn cõi Đông Dương.

Kết luận

  1. Xác lập thành công phương pháp luận khách quan xác định NBĐMM tại 131 trạm khí tượng trên 7 vùng khí hậu Việt Nam giai đoạn 1979–2019 dựa trên tích hợp tín hiệu trạm và tín hiệu vùng qua phân tích MPCA và kiểm định Fisher ANOVA.
  2. Định lượng hóa toàn diện ảnh hưởng của hiện tượng ENSO lên sự biến động liên năm của NBĐMM, chứng minh hiện tượng El Niño làm mùa mưa đến muộn và La Niña làm mùa mưa đến sớm có ý nghĩa thống kê ($p < 0.05$) tại Tây Nguyên và Nam Bộ.
  3. Chứng minh tính ưu việt của phương pháp thống kê - động lực khai thác các trường hoàn lưu (PRMSL, U850, SST) từ mô hình CFSv2 trong dự báo hạn mùa (dưới 6 tháng), đạt điểm kỹ năng $SS_{MAE} > 0$ và hệ số tương quan OLOO đạt tới 0.72.
  4. Tiên phong đánh giá và khẳng định tiềm năng dự báo hạn nội mùa (7–40 ngày) NBĐMM tại Việt Nam sử dụng sản phẩm mưa tổ hợp từ hệ thống mô hình S2S-ECMWF với kỹ năng dự báo vượt trội ở hạn 7–21 ngày.
  5. Mở ra 3 nhánh nghiên cứu học thuật mới: Dự báo hạn nội mùa (S2S) ứng dụng AI/Machine Learning; Dự báo liên hoàn Onset - Cessation - Dry Spells; và Khung liên kết dự báo khí hậu - điều tiết thủy điện - an ninh lương thực.
  6. Cung cấp sản phẩm khoa học có tính thực tiễn cao, sẵn sàng chuyển giao công nghệ cho công tác dự báo nghiệp vụ khí tượng thủy văn quốc gia và phục vụ phát triển kinh tế - xã hội bền vững của đất nước.