Tổng quan về luận án
Biến đổi khí hậu toàn cầu đang tác động sâu sắc đến các hệ thống hoàn lưu khí quyển quy mô hành tinh, trong đó hệ thống gió mùa mùa hè (GMMH) châu Á đóng vai trò là một trong những cơ chế tuần hoàn nhiệt lực - động lực quan trọng nhất. Việt Nam nằm trọn trong vùng giao thoa phức tạp giữa ba phân hệ gió mùa chính: Gió mùa mùa hè Nam Á (Indian Summer Monsoon - ISM), Gió mùa mùa hè Đông Á (East Asian Summer Monsoon - EASM) và Gió mùa mùa hè Tây Bắc Thái Bình Dương (Western North Pacific Summer Monsoon - WNPSM). Trong hệ thống phân vùng khí hậu Việt Nam, Tây Nguyên và Nam Bộ là hai khu vực chịu ảnh hưởng trực tiếp, điển hình và mạnh mẽ nhất của khối khí gió mùa Tây Nam xuất phát từ áp cao cận nhiệt đới Nam Ấn Độ Dương (áp cao Mascarene) và vịnh Bengal.
Mặc dù có vai trò thiết yếu đối với an ninh nguồn nước, hệ sinh thái nông nghiệp và kinh tế - xã hội, các nghiên cứu trước đây tại Việt Nam chủ yếu tiếp cận GMMH thông qua các chuỗi quan trắc mưa cục bộ hoặc khảo sát hoàn lưu quy mô lớn bằng số liệu tái phân tích độ phân giải thấp. Tồn tại một khoảng trống nghiên cứu (research gap) cốt lõi: sự thiếu hụt các nghiên cứu định lượng phân định rành mạch giữa khái niệm "mùa mưa" (rainy season) mang tính địa phương và "mùa gió mùa mùa hè" (summer monsoon season) mang tính động lực hoàn lưu quy mô synop; đồng thời thiếu vắng các mô phỏng chi tiết hóa động lực phân giải cao (high-resolution dynamical downscaling) để dự tính xu thế biến đổi của GMMH và các chỉ số mưa cực đoan trong các giai đoạn tương lai của thế kỷ 21 theo các kịch bản nồng độ khí nhà kính mới của Ủy ban Liên chính phủ về Biến đổi khí hậu (IPCC).
Luận án tiến sĩ chuyên ngành Khí tượng và Khí hậu học (Mã số: 62 44 02 22) của nghiên cứu sinh Ngô Thị Thanh Hương với đề tài "Biến đổi một số đặc trưng gió mùa mùa hè ở Tây Nguyên và Nam Bộ", dưới sự hướng dẫn khoa học của PGS. Vũ Thanh Hằng và PGS. Nguyễn Hướng Điền tại Trường Đại học Khoa học Tự nhiên – Đại học Quốc gia Hà Nội, đã giải quyết trọn vẹn khoảng trống học thuật này thông qua hệ thống câu hỏi nghiên cứu và giả thuyết khoa học:
- Câu hỏi nghiên cứu 1 (RQ1): Ngày bắt đầu (BĐGM), ngày kết thúc (KTGM) và độ dài mùa GMMH tại Tây Nguyên và Nam Bộ có quy luật phân hóa không gian - thời gian như thế nào trong thời kỳ lịch sử 1981–2014, và sự khác biệt bản chất giữa mùa mưa và mùa GMMH là gì?
- Câu hỏi nghiên cứu 2 (RQ2): Các trường hoàn lưu nhiệt động lực quy mô lớn (thông lượng ẩm tích phân thẳng đứng VMF, gradient nhiệt độ tầng đối lưu trên 300 hPa, nhiệt độ mặt nước biển SST) chi phối cơ chế bùng phát và thoái lui của GMMH khu vực như thế nào?
- Câu hỏi nghiên cứu 3 (RQ3): Dưới tác động của kịch bản phát thải trung bình thấp RCP4.5, các đặc trưng động lực GMMH và chế độ mưa cực đoan tại hai khu vực này sẽ biến đổi ra sao trong ba giai đoạn: đầu thế kỷ (2016–2035), giữa thế kỷ (2046–2065) và cuối thế kỷ 21 (2080–2099)?
Hệ thống giả thuyết khoa học tương ứng được kiểm định:
- Giả thuyết 1 (H1): Ngày bắt đầu mùa mưa và ngày bắt đầu GMMH không đồng nhất về mặt thời gian; sự bắt đầu của GMMH đòi hỏi sự kết hợp đồng thời giữa nguồn ẩm đối lưu cục bộ và sự thiết lập ổn định của đới gió Tây Nam tầng thấp (850 hPa).
- Giả thuyết 2 (H2): Sự chuyển pha của GMMH chịu sự điều hòa trực tiếp bởi sự đảo chiều của gradient nhiệt độ kinh hướng tại tầng đối lưu trên (300 hPa) giữa cao nguyên Tây Tạng/lục địa Á-Âu và bắc Ấn Độ Dương, kết hợp với các pha ENSO trên Thái Bình Dương.
- Giả thuyết 3 (H3): Mô hình khí hậu khu vực RegCM4.3 sau khi được hiệu chỉnh sai số bằng hàm phân bố xác suất Gamma có khả năng tái hiện chân thực cơ chế GMMH và đưa ra dự tính đáng tin cậy về sự gia tăng của các cực trị giáng thủy trong tương lai.
Phạm vi không gian nghiên cứu bao trùm 13 trạm quan trắc bề mặt đại diện tại Tây Nguyên (Đăk Tô, Kon Tum, Pleiku, An Khê, Ayunpa, Buôn Hồ, Buôn Ma Thuột, M'Đrắk, Đăk Nông, Đà Lạt, Bảo Lộc, Liên Khương, Cát Tiên) và 4 trạm trọng điểm tại Nam Bộ (Tây Ninh, Đồng Phú, Tân Sơn Hòa, Cần Thơ) trong giai đoạn lịch sử 34 năm (1981–2014), cùng miền tính động lực học khí hậu khu vực ($15^\circ\text{S} - 27^\circ\text{N}, 90^\circ\text{E} - 145^\circ\text{E}$) độ phân giải 25 km cho thời kỳ cơ sở (1986–2005) và tương lai (2016–2099).
Literature Review và Positioning
Khung lý thuyết về gió mùa đã trải qua quá trình tiến hóa học thuật từ các quan điểm nhiệt lực cổ điển đến các mô hình động lực học đối lưu hiện đại. Khromov (1967) đặt nền móng với định nghĩa kinh điển: "Gió mùa là chế độ dòng khí của hoàn lưu khí quyển trên một phạm vi rộng lớn của bề mặt trái đất, trong đó, ở mọi nơi trong khu vực gió mùa, gió thịnh hành chuyển ngược hướng hay gần như ngược hướng từ mùa đông sang mùa hè và từ mùa hè sang mùa đông", với góc chuyển hướng gió thịnh hành $\ge 120^\circ$. Ramage (1971) đã chuẩn hóa 4 tiêu chí định lượng về góc đổi hướng gió ($120^\circ - 180^\circ$), tần suất gió thịnh hành tháng 1 và tháng 7 ($>40%$), sự thay thế xoáy thuận - xoáy nghịch mặt đất và tốc độ gió hợp thành ($>3\text{ m/s}$).
Cấu trúc GMMH Nam Á và Đông Á đã được phân tích sâu sắc bởi Krishnamurti et al. (1976) và Domroes & Peng (1988), nhận diện các thành phần hoàn lưu then chốt: áp cao Mascarene, dòng xiết vượt xích đạo Đông Phi (Somali Jet), rãnh gió mùa Bắc Ấn Độ, áp cao Tây Tạng trên tầng đối lưu 200 hPa, dòng xiết gió Đông nhiệt đới và áp cao cận nhiệt đới Tây Bắc Thái Bình Dương (WPSH). Murakami & Matsumoto (1994), Wang & Ho (2002) và Wang et al. (2003) phân chia hệ thống gió mùa châu Á - Thái Bình Dương thành ba tiểu hệ thống tương tác: ISM ($5^\circ\text{N}-27^\circ\text{N}, 65^\circ\text{E}-105^\circ\text{E}$), WNPSM ($5^\circ\text{N}-22.5^\circ\text{N}, 105^\circ\text{E}-150^\circ\text{E}$) và EASM ($22.5^\circ\text{N}-45^\circ\text{N}, 105^\circ\text{E}-140^\circ\text{E}$), trong đó bán đảo Đông Dương giữ vị trí ranh giới bản lề điều tiết sự tương tác giữa các hệ thống này.
+---------------------------------------------------+
| HỆ THỐNG HOÀN LƯU GIÓ MÙA CHÂU Á |
+---------------------------------------------------+
|
+---------------------------------+---------------------------------+
| | |
v v v
+-------------------+ +-------------------+ +-------------------+
| GMMH NAM Á | | BÁN ĐẢO ĐÔNG DƯƠNG| | GMMH ĐÔNG Á |
| (ISM) | | (Tây Nguyên/NamBộ)| | (EASM) |
+-------------------+ +-------------------+ +-------------------+
| • Áp cao Mascarene| | • Vùng chuyển tiếp| | • Áp cao cận nhiệt|
| • Dòng xiết Somali| | nhiệt động lực | | Tây Bắc TBD |
| • Áp cao Tây Tạng | | • Tương tác gió | | • Front Meiyu/Baiu|
| • Rãnh gió mùa Ấn | | tây nam & đệm | | • Dòng xiết vượt |
+-------------------+ +-------------------+ +-------------------+
Trong y văn tồn tại hai trường phái tranh luận học thuật sâu sắc:
- Tranh luận về sự đồng pha giữa bùng phát gió mùa và bắt đầu mùa mưa: Các nghiên cứu tại Ấn Độ và Đông Á quy mô lớn (Wang & Ho, 2002; Zhang et al., 2002) cho rằng ngày bắt đầu mùa mưa đồng nhất với ngày bùng phát GMMH. Ngược lại, Nguyễn Hướng Điền et al. (2000) và Nguyễn Thị Hiền Thuận (2007) chỉ ra rằng tại bán đảo Đông Dương, gió mùa Tây Nam có thể đã thiết lập ở tầng thấp nhưng mùa mưa chưa thể bắt đầu nếu thiếu các nhiễu động synop (xoáy thuận, rãnh thấp, hội tụ nhiệt đới ITCZ hoặc hiệu ứng nâng cưỡng bức địa hình) đẩy khối khí xích đạo nóng ẩm lên cao đạt mức ngưng kết; hoặc khi một lưỡi áp cao cận nhiệt đới Tây Thái Bình Dương vẫn ngự trị đè nén đối lưu.
- Tranh luận về tác động của biến đổi khí hậu đến thời điểm bắt đầu GMMH: Inoue & Ueda (2011) khi phân tích các mô hình AOGCMs thuộc dự án CMIP3 theo kịch bản SRES A1B kết luận ngày bắt đầu GMMH trên vịnh Bengal và Biển Đông vào cuối thế kỷ 21 sẽ đến muộn hơn 5–10 ngày do sự chậm trễ đảo chiều gradient nhiệt độ kinh hướng tầng đối lưu trên. Trái lại, Zuo & Zhou (2015) và Ruiqing et al. (2016) dựa trên tổ hợp đa mô hình CMIP5 theo kịch bản RCP8.5 lại chỉ ra xu thế GMMH bùng phát sớm hơn đi kèm gia tăng cường độ hoàn lưu ẩm.
Về mặt định vị nghiên cứu quốc tế, công trình của Ngô Thị Thanh Hương đặt cạnh nghiên cứu của Fasullo & Webster (2003) tại khu vực Kerala (Ấn Độ) và Kajikawa & Wang (2012) trên khu vực Biển Đông. Trong khi Fasullo & Webster (2003) sử dụng thông lượng ẩm tích phân VMF thuần túy để xác định ngày bắt đầu trung bình ngày 4/6 ($\pm 7.4$ ngày) và kết thúc ngày 7/9 ($\pm 11$ ngày), và Kajikawa & Wang (2012) áp dụng chỉ số gió vĩ hướng tầng thấp đơn lẻ $U_{850}$ ($5^\circ\text{N}-15^\circ\text{N}, 110^\circ\text{E}-120^\circ\text{E}$ - USCS), luận án này đã khắc phục hạn chế của các chỉ số đơn lẻ bằng việc thiết lập một bộ chỉ tiêu kết hợp chặt chẽ giữa trường động lực gió vĩ hướng tầng 850 hPa và trường giáng thủy mặt đất, phản ánh chính xác tính chất vi khí hậu nhiệt đới gió mùa gió Tây Nam tại miền Nam 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 mở rộng lý thuyết hoàn lưu gió mùa nhiệt đới của Ramage (1971) và Khromov (1967) trong điều kiện mặt đệm phức tạp của bán đảo Đông Dương. Đóng góp lý luận quan trọng của luận án là việc chứng minh sự phân tách (decoupling) và tái kết hợp (coupling) giữa trường nhiệt động lực học quy mô lớn và hoàn lưu ẩm quy mô địa phương.
+------------------------------------------+
| SỰ CHÊNH LỆCH NHIỆT ĐỘ ĐỐI LƯU (300 hPa) |
| Lục địa Á-Âu vs Bắc Ấn Độ Dương |
+------------------------------------------+
|
v
+----------------------------------------------------------------------------------+
| CHỈ TIÊU KẾT HỢP XÁC ĐỊNH GMMH (LUẬN ÁN) |
| |
| +------------------------------------+ +--------------------------------+ |
| | ĐIỀU KIỆN ĐỘNG LỰC | | ĐIỀU KIỆN NHIỆT ẨM | |
| | Gió vĩ hướng $u_{850} > 0.5$ m/s| & | Lượng mưa trượt 5 ngày > 5 mm | |
| | Duy trì $\ge 5$ ngày liên tiếp | | Duy trì $\ge 5$ ngày liên tiếp| |
| +------------------------------------+ +--------------------------------+ |
+----------------------------------------------------------------------------------+
|
v
+------------------------------------------+
| BÙNG PHÁT VÀ TIẾN TRIỂN GMMH CHUẨN |
| TẠI TÂY NGUYÊN VÀ NAM BỘ |
+------------------------------------------+
Mô hình lý thuyết được phát biểu qua ba mệnh đề khoa học:
- Mệnh đề 1 (P1 - Về cơ chế chuyển dịch hoàn lưu): Sự bùng phát GMMH tại Tây Nguyên và Nam Bộ xảy ra đồng thời với quá trình suy yếu và rút lui về phía đông của sống áp cao cận nhiệt đới Tây Bắc Thái Bình Dương, tạo điều kiện cho rãnh thấp gió mùa vịnh Bengal mở rộng sang phía đông và dòng xiết gió tây kinh hướng mực thấp tràn qua kinh độ $100^\circ\text{E}$.
- Mệnh đề 2 (P2 - Về động lực nhiệt đối lưu trên cao): Sự đảo chiều của gradient nhiệt độ kinh hướng tại mực 300 hPa giữa miền lục địa Á-Âu ($20^\circ\text{N}-40^\circ\text{N}, 50^\circ\text{E}-100^\circ\text{E}$) và miền xích đạo ($0^\circ-20^\circ\text{N}, 50^\circ\text{E}-100^\circ\text{E}$) là cơ chế nhiệt lực vĩ mô chi phối thời điểm sớm/muộn của GMMH khu vực.
- Mệnh đề 3 (P3 - Về sự dịch chuyển thập kỷ): Tồn tại bước nhảy thập kỷ (decadal shift) rõ nét trong cơ chế bùng phát GMMH giữa hai thời kỳ 1981–1997 và 1998–2014, gắn liền với sự tăng cường dị thường ấm của nhiệt độ bề mặt nước biển khu vực xích đạo Tây Thái Bình Dương và sự chuyển dịch hoàn lưu Walker.
Khung phân tích độc đáo
Khung phân tích của luận án tích hợp ba trụ cột lý thuyết: 1) Hoàn lưu khí quyển quy mô hành tinh; 2) Động lực học mây đối lưu nhiệt đới; và 3) Thống kê khí hậu phi tham số kết hợp mô hình hóa khí hậu số trị phân giải cao.
Trọng tâm của khung phân tích là Bộ chỉ tiêu kết hợp (Combined Criterion) do tác giả đề xuất:
- Ngày bắt đầu GMMH (BĐGM): Là ngày đầu tiên trong chuỗi 5 ngày liên tiếp đồng thời thỏa mãn: (i) Lượng mưa trung bình trượt 5 ngày đạt $\bar{R}5 > 5.0\text{ mm/ngày}$; và (ii) Thành phần gió vĩ hướng trung bình ngày tại mực 850 hPa đạt $u{850} > 0.5\text{ m/s}$.
- Ngày kết thúc GMMH (KTGM): Được xác định bằng quy trình quét ngược từ cuối năm (31/12) trở về trước thỏa mãn đồng thời các điều kiện ngưỡng trên.
Điều kiện biên lý thuyết (boundary conditions) được xác lập nghiêm ngặt: Miền tính gió vĩ hướng đại diện cho Tây Nguyên là $[12.0^\circ\text{N}-14.5^\circ\text{N}; 107.5^\circ\text{E}-109.0^\circ\text{E}]$ và Nam Bộ là $[9.0^\circ\text{N}-11.5^\circ\text{N}; 104.5^\circ\text{E}-107.0^\circ\text{E}]$. Ngưỡng $u_{850} > 0.5\text{ m/s}$ đóng vai trò là điều kiện đủ để loại bỏ các đợt mưa dông nhiệt cục bộ tiền gió mùa (pre-monsoon convective showers) vốn không phản ánh sự biến đổi của hoàn lưu khí quyển quy mô lớn.
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ới phương pháp luận định lượng quy chuẩn trong khí tượng học động lực. Thiết kế nghiên cứu tích hợp đa cấp độ không gian và thời gian: từ quy mô hoàn lưu toàn cầu (Global reanalysis, GCMs) đến quy mô khu vực (RCM downscaling) và đối chuẩn tại các điểm trạm quan trắc bề mặt.
Mẫu số liệu nghiên cứu được thu thập đồng bộ và kiểm soát chất lượng nghiêm ngặt:
- Số liệu quan trắc bề mặt: Chuỗi số liệu mưa ngày tại 13 trạm Tây Nguyên và 4 trạm Nam Bộ trong 34 năm liên tục (1981–2014) từ Trung tâm Khí tượng Thủy văn Quốc gia.
- Số liệu tái phân tích hoàn lưu: Số liệu NCEP/NCAR phân giải $2.5^\circ \times 2.5^\circ$ (1981–2014) và số liệu ERA-Interim của ECMWF phân giải $1.5^\circ \times 1.5^\circ$ (1986–2005).
- Số liệu mô hình khí hậu: Mô hình khí hậu toàn cầu CNRM-CM5 (Pháp) thuộc dự án CMIP5 cung cấp điều kiện biên cho mô hình khu vực RegCM4.3 dưới kịch bản phát thải RCP4.5.
Quy trình nghiên cứu rigorous
Quy trình mô phỏng số trị chi tiết hóa động lực kế thừa các cấu hình tối ưu từ dự án CORDEX-SEA (Juneng et al., 2016; Ngo Duc Thanh et al., 2016):
- Mô hình khí hậu khu vực: RegCM phiên bản 4.3 (Regional Climate Model version 4.3).
- Cấu trúc thẳng đứng: 18 mực sigma ($\sigma$) với đỉnh mô hình tại 5 hPa (mb).
- Miền tính tích phân: Tọa độ $15^\circ\text{S} - 27^\circ\text{N}, 90^\circ\text{E} - 145^\circ\text{E}$, bước lưới ngang 25 km.
- Cấu hình vật lý tối ưu (Physics Parameterizations):
- Sơ đồ bức xạ sóng ngắn và sóng dài CCM3.
- Sơ đồ tương tác bề mặt đất - sinh quyển BATS1e (Biosphere-Atmosphere Transfer Scheme 1e).
- Sơ đồ tham số hóa vận chuyển ẩm quy mô lớn SUBEX (Subgrid Explicit Moisture Scheme).
- Sơ đồ đối lưu và trao đổi bề mặt đại dương MIT-Emanuel (được chứng minh có kỹ năng mô phỏng mưa đối lưu vượt trội trên vùng biển nhiệt đới Đông Nam Á so với các sơ đồ Grell, Kuo hay Fritsch-Chappell).
+-----------------------------------------------------------------------------------+
| QUY TRÌNH MÔ HÌNH HÓA VÀ HIỆU CHỈNH SAI SỐ KHÍ HẬU |
+-----------------------------------------------------------------------------------+
| Điều kiện biên: ERA-Interim (1.5°) & CNRM-CM5 (GCM) |
| | |
| v |
| Mô hình khu vực: RegCM4.3 (Miền tính: 15°S-27°N, 90°-145°E; 18 mực Sigma; 25 km) |
| Tham số hóa: MIT-Emanuel (Đối lưu) + BATS1e (Mặt đất) + SUBEX (Ẩm) + CCM3 (Bức xạ)|
| | |
| v |
| Sản phẩm thô: NetCDF mưa ngày & gió vĩ hướng 850 hPa (1986-2005 & 2016-2099) |
| | |
| v |
| Hiệu chỉnh sai số hệ thống: Phương pháp Gamma Quantile Mapping (Piani et al. 2010)|
| | |
| v |
| Phân tích thống kê: Kiểm nghiệm Mann-Kendall + Sen's slope + Student t-test (95%)|
+-----------------------------------------------------------------------------------+
Tính hợp thức và độ tin cậy được bảo đảm thông qua tam giác đạc phương pháp (methodological triangulation): kết quả của bộ chỉ tiêu kết hợp được đối chiếu song song với chỉ tiêu thuần mưa của Zhang et al. (2002), chỉ tiêu gió của Cheang et al. (1988) và chỉ số USCS của Wang et al. (2004).
Data và phân tích
Nhằm triệt tiêu sai số hệ thống (systematic bias) của mô hình số trị, luận án ứng dụng phương pháp hiệu chỉnh hàm phân bố xác suất lý thuyết Gamma (Gamma Quantile Mapping) theo Piani et al. (2010) và Argüeso et al. (2013). Hàm mật độ xác suất Gamma hai tham số ($k$: hình dạng, $\theta$: tỷ lệ) của lượng mưa ngày mô phỏng ($M$) được ánh xạ chuẩn hóa về hàm phân bố thực nghiệm của chuỗi quan trắc ($O$):
$$f(x; k, \theta) = \frac{x^{k-1} e^{-\frac{x}{\theta}}}{\theta^k \Gamma(k)} \quad \text{với } x > 0$$
$$P_{corr} = F^{-1}{obs}\left( F{mod}(P_{raw}; k_{mod}, \theta_{mod}); k_{obs}, \theta_{obs} \right)$$
Hệ thống 8 chỉ số mưa cực đoan theo chuẩn của Tổ chức Khí tượng Thế giới (WMO/ETCCDI) và Peterson et al. (2001) được tính toán chi tiết:
| Ký hiệu chỉ số |
Tên chỉ số khí hậu |
Định nghĩa toán học / Ngưỡng khí hậu |
Đơn vị |
| WD |
Wet Days |
Số ngày trong năm/mùa có lượng mưa ngày $R \ge 1.0\text{ mm}$ |
ngày |
| PRCPTOT |
Total Precipitation |
Tổng lượng mưa tích lũy từ các ngày mưa ẩm ($R \ge 1.0\text{ mm}$) |
mm |
| SDII |
Simple Daily Intensity Index |
Cường độ mưa trung bình ngày mưa: $\text{SDII} = \text{PRCPTOT} / \text{WD}$ |
mm/ngày |
| RX1day |
Max 1-day Precipitation |
Lượng mưa cực đại tích lũy trong 1 ngày của năm hoặc mùa |
mm |
| RX5day |
Max 5-day Precipitation |
Lượng mưa cực đại tích lũy trong 5 ngày liên tiếp |
mm |
| R50 |
Heavy Precipitation Days |
Số ngày có lượng mưa rất lớn $R \ge 50.0\text{ mm}$ |
ngày |
| CDD |
Consecutive Dry Days |
Độ dài chuỗi ngày khô hạn liên tiếp dài nhất ($R < 1.0\text{ mm}$) |
ngày |
| CWD |
Consecutive Wet Days |
Độ dài chuỗi ngày mưa ẩm liên tiếp dài nhất ($R \ge 1.0\text{ mm}$) |
ngày |
Xu thế biến đổi dài hạn được lượng hóa bằng kiểm nghiệm phi tham số Mann-Kendall ($S, \tau, \text{Var}(S)$) và độ dốc xu thế trung vị Sen (Sen’s Slope Estimator $Q$):
$$S = \sum_{k=1}^{n-1} \sum_{j=k+1}^n \text{sign}(x_j - x_k)$$
$$Q = \text{median}\left( \frac{x_j - x_k}{j - k} \right), \quad \forall j > k$$
Mức độ khác biệt có ý nghĩa thống kê giữa các thời kỳ (1998–2014 so với 1981–1997; và các giai đoạn tương lai so với thời kỳ cơ sở 1986–2005) được đánh giá bằng kiểm nghiệm thống kê Student hai mẫu ($t$-test) ở mức ý nghĩa $\alpha = 0.05$ (độ tin cậy 95%).
Phát hiện đột phá và implications
Những phát hiện then chốt
- Phát hiện sự lệch pha giữa mùa mưa và mùa GMMH: Phân tích thực nghiệm trên 34 năm chuỗi số liệu khẳng định ngày bắt đầu mùa mưa (BĐMM) theo chỉ tiêu Zhang et al. (2002) tại Tây Nguyên và Nam Bộ thường diễn ra vào cuối tháng 4 đến đầu tháng 5, sớm hơn ngày bùng phát gió mùa chính thức (BĐGM) từ 1 đến 2 pentad (khoảng 5–10 ngày). Ngày bắt đầu GMMH trung bình chuẩn hóa theo chỉ tiêu kết hợp diễn ra vào giữa tháng 5 (pentad 27–28, tương đương ngày 13–18/5), và kết thúc vào giữa tháng 10 (pentad 57–58, ngày 10–15/10), với thời gian kéo dài mùa gió mùa trung bình từ 150 đến 160 ngày.
- Xác lập bước nhảy thập kỷ về thời điểm bắt đầu GMMH: So sánh giữa hai giai đoạn 1981–1997 và 1998–2014 phát hiện xu thế bùng phát GMMH có dấu hiệu dịch chuyển sớm hơn trong giai đoạn sau. Nguyên nhân nhiệt động lực học cốt lõi là sự gia tăng chênh lệch nhiệt độ tầng đối lưu trên (300 hPa) giữa miền lục địa Á-Âu ($20^\circ\text{N}-40^\circ\text{N}, 50^\circ\text{E}-100^\circ\text{E}$) và vùng nhiệt đới ($0^\circ-20^\circ\text{N}, 50^\circ\text{E}-100^\circ\text{E}$) diễn ra sớm hơn trong mùa xuân, thúc đẩy sự hình thành sớm của dòng xiết gió Tây Nam và thông lượng ẩm VMF từ vịnh Bengal sang Biển Đông.
- Tác động bất đối xứng của hiện tượng ENSO: Tương quan giữa nhiệt độ mặt nước biển (SST) các tháng mùa đông (DJF) và mùa xuân (MAM) tại vùng Nino 3.4 với chuỗi thời gian BĐGM/KTGM cho thấy: trong các năm El Niño (ví dụ năm 1987, 1998), áp cao cận nhiệt đới Tây Thái Bình Dương duy trì cường độ mạnh và lấn sâu về phía tây, làm chậm ngày bắt đầu GMMH, mùa gió mùa ngắn hơn và tổng lượng mưa mùa GMMH (SMrain) suy giảm nghiêm trọng. Ngược lại, trong các năm La Niña (như năm 2000), xoáy kép xích đạo (vortices twin) phát triển mạnh trên vùng biển Ấn Độ Dương, đẩy gió mùa bùng phát sớm và kéo dài với lượng mưa dồi dào.
- Dự tính sự gia tăng vượt trội của các cực trị giáng thủy tương lai (Kịch bản RCP4.5): Mô phỏng RegCM4.3 chỉ ra rằng đến giữa thế kỷ (2046–2065) và cuối thế kỷ 21 (2080–2099), mặc dù thời điểm bắt đầu và kết thúc GMMH có mức biến động trung bình không quá lớn ($\pm 3$ đến 7 ngày), nhưng cấu trúc phân bố mưa mùa gió mùa thay đổi sâu sắc theo hướng cực đoan hóa:
- Chỉ số mưa 1 ngày lớn nhất (RX1day) và 5 ngày lớn nhất (RX5day) tăng phổ biến từ $10% - 25%$ trên toàn khu vực Tây Nguyên (đặc biệt tại các sườn đón gió như Bảo Lộc, Đăk Nông) và Nam Bộ.
- Số ngày mưa rất lớn ($R \ge 50\text{ mm}$ - R50) gia tăng có ý nghĩa thống kê ($p < 0.05$).
- Độ dài đợt khô hạn liên tiếp lớn nhất (CDD) trong mùa GMMH có xu hướng giảm nhẹ hoặc không đổi, nhưng cường độ mưa ngày (SDII) tăng rõ rệt, chứng minh hiện tượng tập trung giáng thủy cường suất cao trong các đợt mưa ngắn.
+---------------------------------------------------------------------------------------+
| BẢNG TỔNG HỢP BIẾN ĐỔI CÁC ĐẶC TRƯNG GMMH VÀ MƯA CỰC ĐOAN |
+--------------------------+------------------------------+-----------------------------+
| Đặc trưng khí hậu | Thời kỳ Lịch sử (1981-2014) | Dự tính Tương lai (RCP4.5) |
+--------------------------+------------------------------+-----------------------------+
| Ngày bắt đầu (BĐGM) | Giữa tháng 5 (Pentad 27-28) | Biến động sớm hơn 3-5 ngày |
| Ngày kết thúc (KTGM) | Giữa tháng 10 (Pentad 57-58) | Có xu hướng rút lui muộn hơn|
| Độ dài mùa gió mùa | 150 - 160 ngày | Kéo dài thêm 4 - 8 ngày |
| Lượng mưa cực đại RX1day | Tăng nhẹ (Kiểm nghiệm Sen) | Tăng mạnh 12% - 25% (p<0.05)|
| Lượng mưa cực đại RX5day | Xu thế tăng tại Bảo Lộc, M'Đrắk| Tăng phổ biến 10% - 20% |
| Cường độ mưa SDII | Tăng tại phía Nam | Tăng 0.5 - 1.8 mm/ngày/thập kỷ|
+--------------------------+------------------------------+-----------------------------+
Implications đa chiều
- Về mặt lý luận khoa học: Cung cấp bằng chứng định lượng khẳng định bán đảo Đông Dương đóng vai trò là "cầu nối động lực" trung gian trong việc giải phóng ẩn nhiệt đối lưu, kích hoạt sự chuyển dịch mùa của hệ thống gió mùa mùa hè châu Á.
- Về mặt phương pháp luận: Chuẩn hóa quy trình chi tiết hóa động lực học kết hợp hiệu chỉnh sai số Gamma hai tham số, thiết lập tiền lệ phương pháp luận tin cậy cho các nghiên cứu dự tính khí hậu thủy văn độ phân giải cao tại các khu vực nhiệt đới có địa hình phức tạp.
- Về mặt thực tiễn và chính sách:
- Nông nghiệp: Cung cấp cơ sở khoa học chính xác để tái cấu trúc lịch thời vụ gieo trồng lúa vụ Hè Thu tại Nam Bộ và chăm sóc cây công nghiệp dài ngày (cà phê, hồ tiêu, cao su) tại Tây Nguyên, tránh các đợt hạn bà chằng đầu mùa hoặc úng lụt cuối mùa.
- Quy hoạch tài nguyên nước & Năng lượng tái tạo: Cảnh báo mức gia tăng của các cực trị RX1day và RX5day hỗ trợ việc tính toán lại lũ thiết kế và quy trình vận hành liên hồ chứa thủy điện - thủy lợi trên các lưu vực sông lớn (sông Sê San, Srepok, sông Ba, hệ thống sông Đồng Nai).
Limitations và Future Research
Luận án thẳng thắn chỉ ra các giới hạn nghiên cứu mang tính khách quan:
- Hạn chế về tổ hợp mô hình biên: Nghiên cứu mới chỉ sử dụng điều kiện biên từ một mô hình khí hậu toàn cầu duy nhất (CNRM-CM5) cho các dự tính tương lai do hạn chế về tài nguyên tính toán tại thời điểm thực hiện, chưa triển khai được phương pháp tổ hợp đa mô hình (multi-model ensemble) để lượng hóa đầy đủ biên độ bất định (uncertainty range).
- Mạng lưới trạm quan trắc bề mặt tại Nam Bộ: Số lượng trạm quan trắc khí tượng bề mặt dài hạn tại Nam Bộ được đưa vào tính toán chi tiết còn tương đối thưa (4 trạm đại diện), chưa bao phủ hết tính biến động vi khí hậu ven biển Tây Nam và bán đảo Cà Mau.
- Độ phân giải mô hình khu vực: Lưới tính 25 km tuy đã cải thiện đáng kể so với GCMs ($150-200\text{ km}$) nhưng vẫn chưa phân giải trực tiếp (convection-permitting) được hoàn toàn các đối lưu nhiệt địa phương quy mô đám mây đối lưu riêng lẻ trong địa hình đồi núi hiểm trở dạng hẻm vực.
Chương trình nghiên cứu tương lai (Future Research Agenda) cần tập trung vào:
- Mở rộng tích hợp các kịch bản phát thải mới nhất theo khung CMIP6/CORDEX-SEA Phase 2 (các kịch bản SSP1-2.6, SSP2-4.5, SSP5-8.5).
- Ứng dụng mô hình khí hậu khu vực phân giải siêu cao (< 5 km - Convection Permitting Modeling) nhằm mô phỏng chính xác tương tác giữa địa hình sườn đón gió dãy Trường Sơn Nam và cơ chế đối lưu giông sét nhiệt đới.
- Phát triển hệ thống dự báo hạn mùa nghiệp vụ về ngày bùng phát GMMH dựa trên các chỉ số động lực học phân tích thời gian thực.
Tác động và ảnh hưởng
- Tác động học thuật (Academic Impact): Luận án là một trong những công trình tiên phong tại Việt Nam tích hợp hoàn chỉnh từ lý thuyết hoàn lưu nhiệt động lực quy mô hành tinh đến phương pháp mô phỏng số trị phân giải cao về GMMH. Các bài báo công bố từ luận án trên các tạp chí chuyên ngành uy tín trong nước và quốc tế tạo nền tảng trích dẫn vững chắc cho các nghiên cứu sinh và học viên cao học ngành Khí tượng, Thủy văn và Biến đổi khí hậu.
- Chuyển đổi ngành công nghiệp và năng lượng: Cung cấp thông số đầu vào quan trọng cho các kỹ sư thủy văn và các nhà phát triển năng lượng tái tạo (thủy điện, điện gió, điện mặt trời) trong việc tối ưu hóa mô hình dự báo dòng chảy vào hồ chứa và đánh giá biến động trường gió tầng thấp tại Tây Nguyên và Nam Bộ.
- Tác động chính sách cấp quốc gia và địa phương: Đóng góp trực tiếp cơ sở dữ liệu khoa học cho Cục Biến đổi khí hậu (Bộ Tài nguyên và Môi trường) trong việc cập nhật Kịch bản Biến đổi khí hậu và Nước biển dâng cho Việt Nam; hỗ trợ Ủy ban Nhân dân các tỉnh Tây Nguyên và Tây Nam Bộ xây dựng Kế hoạch hành động ứng phó với biến đổi khí hậu giai đoạn 2021–2030, tầm nhìn 2050.
Đối tượng hưởng lợi
+-------------------------------------------------------------------------------------+
| ĐỐI TƯỢNG HƯỞNG LỢI TRỰC TIẾP |
+-------------------------------------------------------------------------------------+
| |
| +------------------------------------+ +------------------------------------+ |
| | GIỚI VIỆN - TRƯỜNG | | CƠ QUAN HOẠCH ĐỊNH | |
| | • Nghiên cứu sinh & Giảng viên | | • Bộ Tài nguyên và Môi trường | |
| | • Khung lý thuyết & chỉ tiêu kết | | • Cục Quản lý Tài nguyên nước | |
| | hợp gió - mưa chuẩn hóa | | • Dữ liệu kịch bản BĐKH cấp tỉnh | |
| +------------------------------------+ +------------------------------------+ |
| |
| +------------------------------------+ +------------------------------------+ |
| | KHỐI DOANH NGHIỆP | | CỘNG ĐỒNG NÔNG NGHIỆP | |
| | • Tập đoàn Điện lực Việt Nam (EVN) | | • Ban Chỉ đạo Phòng chống thiên tai| |
| | • Chủ hồ chứa thủy điện/thủy lợi | | • Hợp tác xã trồng lúa & cây công | |
| | • Thiết kế dung tích phòng lũ hồ | | nghiệp: Tránh hạn - né mặn | |
| +------------------------------------+ +------------------------------------+ |
+-------------------------------------------------------------------------------------+
- Nghiên cứu sinh và Giới nghiên cứu Khí hậu: Tiếp cận bộ công cụ phân tích hoàn lưu synop, phương pháp hiệu chỉnh sai số Gamma cho mô hình khí hậu khu vực, và bộ chỉ số mưa cực đoan chuẩn hóa WMO.
- Cơ quan Dự báo Khí tượng Thủy văn Quốc gia: Ứng dụng bộ chỉ tiêu kết hợp giữa gió vĩ hướng tầng thấp và lượng mưa để nâng cao độ tin cậy trong các bản tin nhận định thời điểm bắt đầu/kết thúc mùa mưa và mùa GMMH nghiệp vụ hàng năm.
- Ngành Nông nghiệp và Phát triển Nông thôn: Hàng triệu hộ nông dân canh tác lúa tại Đồng bằng sông Cửu Long và cây công nghiệp tại Tây Nguyên hưởng lợi gián tiếp thông qua việc điều chỉnh khung lịch thời vụ canh tác chính xác theo xu thế biến đổi của mùa gió mùa.
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ì và đã mở rộng lý thuyết nào?
Đóng góp độc đáo nhất là việc làm sáng tỏ cơ chế nhiệt động lực học phân tách giữa "mùa mưa" và "mùa gió mùa mùa hè", từ đó mở rộng lý thuyết gió mùa nhiệt đới cổ điển của Khromov (1967) và Ramage (1971). Luận án chứng minh rằng trên bán đảo Đông Dương, gió mùa bùng phát không đơn thuần là sự xuất hiện của mưa mà là sự thiết lập đồng bộ của đới gió Tây Nam ổn định tầng thấp ($u_{850} > 0.5\text{ m/s}$) dưới sự kích hoạt của việc đảo chiều gradient nhiệt độ kinh hướng tầng đối lưu 300 hPa giữa cao nguyên Tây Tạng và bắc Ấn Độ Dương.
2. Đột phá phương pháp luận của luận án khi so sánh với các nghiên cứu quốc tế trước đây?
So với phương pháp thuần mưa của Zhang et al. (2002) hay thuần gió vĩ hướng USCS của Wang et al. (2004) và Fasullo & Webster (2003), luận án đã thiết lập thành công Bộ chỉ tiêu kết hợp đồng thời ngưỡng động lực ($u_{850} > 0.5\text{ m/s}$) và ngưỡng ẩm ($\bar{R}_5 > 5\text{ mm/ngày}$). Đồng thời, luận án tiên phong ứng dụng kỹ thuật hiệu chỉnh sai số phân bố xác suất Gamma cho sản phẩm mô hình RegCM4.3 phân giải cao 25 km tại Việt Nam, khắc phục triệt để sai số hệ thống do địa hình sườn dốc gây ra.
3. Phát hiện nào gây bất ngờ nhất từ chuỗi số liệu thực nghiệm và mô hình?
Phát hiện bất ngờ nhất là sự tồn tại của bước nhảy thập kỷ giữa hai giai đoạn 1981–1997 và 1998–2014: GMMH có xu hướng bùng phát sớm hơn trong giai đoạn sau nhưng độ dài mùa gió mùa không kéo dài tương ứng mà có sự phân hóa mạnh mẽ, đi kèm với sự suy giảm số ngày mưa ẩm (WD) nhưng lại tăng mạnh cường độ mưa ngày cực đoan (SDII, RX1day) và số ngày mưa rất lớn (R50).
4. Giao thức tái lập nghiên cứu (Replication Protocol) có được cung cấp đầy đủ không?
Luận án cung cấp tường minh toàn bộ giao thức thực nghiệm: từ tọa độ kinh vĩ độ của 17 trạm quan trắc bề mặt, miền tính kinh vĩ độ mô hình RegCM4.3 ($15^\circ\text{S}-27^\circ\text{N}, 90^\circ\text{E}-145^\circ\text{E}$), 18 mực sigma, cấu hình các sơ đồ tham số hóa vật lý (MIT-Emanuel, BATS1e, SUBEX, CCM3), các phương trình toán học hiệu chỉnh Gamma và kiểm nghiệm Mann-Kendall / Student $t$-test. Bất kỳ nhà nghiên cứu nào cũng có thể tái lập hoàn toàn quy trình này trên hệ thống máy tính hiệu năng cao.
5. Luận án đã phác thảo chương trình nghiên cứu 10 năm tiếp theo như thế nào?
Chương trình 10 năm được định hướng xoay quanh: 1) Tích hợp mô hình khí hậu phân giải cực cao convection-permitting (< 3 km) để mô phỏng tương tác vi khí hậu - địa hình; 2) Sử dụng tổ hợp đa mô hình CMIP6 để giảm thiểu độ bất định của các kịch bản cực đoan; và 3) Xây dựng hệ thống đồng hóa số liệu và học máy (machine learning) hỗ trợ dự báo sớm sự bùng phát GMMH phục vụ vận hành thông minh hồ chứa thủy điện.
Kết luận
- Chuẩn hóa thành công Bộ chỉ tiêu kết hợp: Luận án đã đề xuất và kiểm chứng thành công chỉ tiêu kết hợp giữa gió vĩ hướng tầng thấp ($u_{850} > 0.5\text{ m/s}$) và lượng mưa trượt 5 ngày ($\bar{R}_5 > 5.0\text{ mm/ngày}$), phân định rành mạch giữa mùa mưa đối lưu cục bộ và mùa GMMH quy mô synop tại Tây Nguyên và Nam Bộ.
- Lượng hóa toàn diện quy luật khí hậu lịch sử 1981–2014: Xác định ngày bắt đầu GMMH trung bình rơi vào giữa tháng 5 (pentad 27–28), kết thúc vào giữa tháng 10 (pentad 57–58) với độ dài mùa trung bình 150–160 ngày; đồng thời phát hiện bước nhảy thập kỷ bùng phát sớm hơn sau năm 1997.
- Làm sáng tỏ cơ chế tương tác hoàn lưu đa quy mô: Chứng minh mối liên kết động lực chặt chẽ giữa sự đảo chiều gradient nhiệt độ 300 hPa trên cao, sự suy thoái của áp cao cận nhiệt đới Tây Bắc Thái Bình Dương, thông lượng ẩm VMF và tác động chi phối của hiện tượng ENSO đến tính biến động của mùa gió mùa.
- Làm chủ công nghệ mô phỏng khí hậu khu vực phân giải cao: Ứng dụng xuất sắc mô hình RegCM4.3 tích hợp sơ đồ đối lưu MIT-Emanuel và phương pháp hiệu chỉnh sai số Gamma để tái hiện chân thực trường hoàn lưu và giáng thủy khu vực.
- Cung cấp bức tranh dự tính tin cậy theo kịch bản RCP4.5: Cảnh báo xu thế gia tăng có ý nghĩa thống kê của các chỉ số giáng thủy cực đoan (RX1day, RX5day, R50, SDII) đến cuối thế kỷ 21, cung cấp cơ sở khoa học nền tảng cho chiến lược thích ứng biến đổi khí hậu và an ninh nguồn nước quốc gia.
