Tổng quan nghiên cứu

Việt Nam nằm trọn trong khu vực nhiệt đới gió mùa châu Á, sở hữu đường bờ biển dài hơn 3.260 km và địa hình phức tạp, hàng năm phải gánh chịu từ 10 đến 15 đợt mưa lớn diện rộng cùng nhiều hình thái thời tiết cực đoan. Các đợt mưa lớn kéo dài thường gây ra ngập lụt nghiêm trọng tại các đô thị đồng bằng, lũ quét và sạt lở đất tại các vùng núi dốc, làm thiệt hại kinh tế ước tính hàng nghìn tỷ đồng mỗi năm. Vấn đề nghiên cứu cốt lõi được đặt ra là làm thế nào để nâng cao độ chính xác trong dự báo hạn vừa và hạn mở rộng đối với các đợt mưa lớn diện rộng, khi các mô hình thời tiết quy mô synop truyền thống thường gặp giới hạn về thời gian dự báo trước từ 3 đến 5 ngày.

Mục tiêu cụ thể của đề tài là khảo sát toàn diện quy luật hoạt động của dao động Madden – Julian, xác định các đặc trưng động lực học khí quyển và làm sáng tỏ mối liên hệ định lượng giữa 8 pha hoạt động của dao động này với tần suất, cường độ và phân bố không gian của các đợt mưa lớn diện rộng trên toàn lãnh thổ Việt Nam. Phạm vi nghiên cứu về mặt không gian bao phủ dải vĩ độ từ 15°N đến 15°S trên quy mô toàn cầu đối với trường hoàn lưu khí quyển và 6 khu vực khí hậu trọng điểm của Việt Nam.

Về mặt thời gian, luận văn sử dụng chuỗi số liệu khí hậu toàn cầu kéo dài 42 năm từ năm 1974 đến năm 2016, cùng chuỗi quan trắc các đợt mưa lớn diện rộng chi tiết trong 14 năm từ năm 2000 đến năm 2013. Ý nghĩa thực tiễn của công trình thể hiện qua việc cung cấp các căn cứ khoa học giúp nâng cao năng lực cảnh báo sớm mưa lớn trước từ 10 đến 20 ngày, hỗ trợ giảm thiểu tối đa rủi ro thiên tai và tối ưu hóa quy trình điều tiết hồ chứa thủy điện, thủy lợi.

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

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng dựa trên lý thuyết dao động nội mùa nhiệt đới do hai nhà khoa học Roland Madden và Paul Julian phát hiện vào các năm 1971 và 1972. Dao động Madden – Julian là hiện tượng ghép đôi động lực học giữa hoàn lưu khí quyển quy mô lớn và hoạt động đối lưu nhiệt đới sâu, đặc trưng bởi chu kỳ biến động từ 30 đến 90 ngày, phổ biến nhất trong khoảng 40 đến 50 ngày. Hệ thống này có quy mô không gian theo phương ngang rộng từ 12.000 km đến 20.000 km, lan truyền có hệ thống từ tây sang đông dọc theo xích đạo với vận tốc trung bình khoảng 5 m/s, nhanh hơn tốc độ chuyển động của các khối mây riêng lẻ nhưng chậm hơn đáng kể so với sóng Kelvin tự do đạt từ 15 đến 17 m/s.

Khung mô hình nghiên cứu thứ hai được áp dụng là mô hình không gian pha hai chiều của Wheeler và Hendon công bố năm 2004. Mô hình này phân chia chu trình tiến triển của dao động thành 8 pha không gian tương ứng với vị trí tâm đối lưu tăng cường: pha 1 (Bán cầu Tây và châu Phi), pha 2 và 3 (Ấn Độ Dương), pha 4 và 5 (khu vực Lục địa biển), pha 6 và 7 (Tây Thái Bình Dương), pha 8 (tiến vào Tây Bán Cầu). Các khái niệm then chốt gồm: Phát xạ sóng dài đỉnh khí quyển OLR đo lường mức độ phát triển của mây đối lưu, trường gió vĩ hướng tầng đối lưu thấp tại mực 850 hPa, trường gió vĩ hướng tầng đối lưu trên tại mực 200 hPa, và chỉ số cường độ dao động tổng hợp RMM với ngưỡng hoạt động mạnh được xác lập khi giá trị vượt qua 1,0.

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

Nguồn dữ liệu đầu vào bao gồm: bộ số liệu tái phân tích hoàn lưu khí quyển toàn cầu từ Trung tâm Dự báo Môi trường Quốc gia Mỹ và Trung tâm Nghiên cứu Khí quyển Quốc gia Mỹ NCEP/NCAR trên lưới độ phân giải 2,5° x 2,5°; dữ liệu quan trắc phát xạ sóng dài OLR từ Cơ quan Quản lý Khí quyển và Đại dương Quốc gia Mỹ NOAA từ năm 1981 đến năm 2013; chuỗi chỉ số RMM thực tế từ Cục Khí tượng Úc BoM gồm 9.299 ngày có dao động hoạt động mạnh giai đoạn 1974 – 2016; và tài liệu tổng kết các đợt mưa lớn diện rộng từ Trung tâm Dự báo Khí tượng Thủy văn Quốc gia giai đoạn 2000 – 2013. Cỡ mẫu nghiên cứu bao gồm hàng trăm đợt mưa lớn diện rộng được chọn mẫu toàn bộ theo tiêu chí chuẩn: mưa xảy ra trên 50% diện tích của một vùng địa lý và trùng khớp với những ngày có chỉ số dao động lớn hơn 1,0.

Phương pháp phân tích chủ đạo là phân tích hàm trực giao tự nhiên EOF kết hợp kỹ thuật lọc dải tần số Lanczos chu kỳ 20 đến 100 ngày với 201 trọng số, được lập trình tự động hóa thông qua ngôn ngữ NCL để tái tạo bộ chỉ số ReCal. Lý do lựa chọn phương pháp phân tích trực giao trên dữ liệu lọc dải tần là khả năng triệt tiêu triệt để các sóng nhiễu tần số cao dưới 10 ngày và khử dao động liên năm của chu kỳ ENSO, giúp hai vector riêng đầu tiên EOF1 và EOF2 nắm bắt trọn vẹn 43,2% tổng phương sai biến động của hoàn lưu nhiệt đới.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Thứ nhất, nghiên cứu đã xác định quy luật phân bố số ngày hoạt động mạnh của dao động nội mùa. Thống kê 9.299 ngày trong chuỗi 1974 – 2016 của BoM cho thấy số ngày xuất hiện nhiều nhất ở pha 5 với tỷ lệ 13,53% và ít nhất ở pha 8 với tỷ lệ 11,53%. Khi so sánh với chuỗi chỉ số ReCal giai đoạn 1981 – 2013, hệ số tương quan giữa hai thành phần trực giao PC1-RMM1 và PC2-RMM2 đạt mức rất cao, dao động từ 0,717 đến 0,945 ở hầu hết các năm, khẳng định tính ổn định vượt trội của quy trình lọc dữ liệu.

Thứ hai, tần suất xuất hiện mưa lớn diện rộng tại Việt Nam có sự phụ thuộc mật thiết vào vị trí các pha. Trong giai đoạn 2000 – 2013, các đợt mưa lớn diện rộng tập trung với mật độ cao áp đảo khi dao động hoạt động tại pha 4, pha 5, pha 6 và pha 7, chiếm trên 68% tổng số đợt mưa lớn có dao động mạnh. Ngược lại, tại pha 1 và pha 8, tần suất xảy ra mưa lớn diện rộng suy giảm rõ rệt, chỉ chiếm khoảng 12% đến 15%.

Thứ ba, sự phân bố mưa lớn có tính phân hóa không gian và mùa vụ sâu sắc. Mưa lớn diện rộng có liên quan đến dao động xuất hiện chủ yếu từ tháng 5 đến tháng 12 hàng năm. Tại khu vực Trung Trung Bộ và Nam Trung Bộ, mưa lớn cực đoan đạt đỉnh vào các tháng 10 và 11 khi dao động ở pha 6 và 7. Điển hình như đợt mưa lịch sử từ ngày 12 đến ngày 19 tháng 10 năm 2010 tại Nghệ An và Quảng Bình ghi nhận lượng mưa tích lũy từ 500 mm đến 900 mm, diễn ra trùng khớp với thời điểm dao động ở pha 6 có cường độ RMM đạt mức 1,84.

Thảo luận kết quả

Cơ chế vật lý dẫn đến sự gia tăng đột biến của lượng mưa tại Việt Nam bắt nguồn từ sự dịch chuyển của vùng hội tụ ẩm và chuyển động thẳng đứng. Khi tâm đối lưu của dao động di chuyển qua khu vực Lục địa biển ở pha 4 và 5, đới gió tây nam tầng thấp tại mực 850 hPa được tăng cường mạnh mẽ, vận chuyển khối lượng hơi nước khổng lồ từ Ấn Độ Dương vào Nam Bộ và Tây Nguyên. Khi dao động tiến sang Tây Thái Bình Dương ở pha 6 và 7, trường chuẩn sai khí áp bề mặt chuyển sang giá trị âm, kích hoạt rãnh áp thấp và tạo điều kiện thuận lợi cho các xoáy thuận nhiệt đới kết hợp với không khí lạnh tràn xuống từ phía bắc, gây mưa cực đoan kéo dài tại miền Trung.

Kết quả này hoàn toàn nhất quán với các công bố quốc tế của Wheeler và Hendon năm 2009 tại khu vực châu Úc cũng như Zhang năm 2009 tại khu vực Đông Nam Trung Quốc, đồng thời bổ sung bằng chứng thực nghiệm quan trọng về ảnh hưởng của hiệu ứng địa hình chắn gió dãy Trường Sơn tại Việt Nam. Về mặt trình bày dữ liệu, các kết quả nghiên cứu có thể được trực quan hóa tối ưu thông qua biểu đồ không gian trạng thái hai chiều giữa PC1 và PC2 mô tả quỹ đạo ngược chiều kim đồng hồ của dao động, biểu đồ cột so sánh tần suất mưa theo 8 pha, và bản đồ trường chuẩn sai khí áp mực biển với các đường đẳng trị cách nhau 20 hPa kết hợp vector gió tầng 850 hPa để hiển thị rõ các tâm xoáy quy mô hành tinh.

Đề xuất và khuyến nghị

Thứ nhất, tích hợp module giám sát và dự báo chỉ số pha RMM và ReCal thời gian thực vào quy trình nghiệp vụ của Trung tâm Dự báo Khí tượng Thủy văn Quốc gia. Hành động cụ thể là xây dựng công cụ giải đoán tự động tích hợp trong các mô hình dự báo số trị toàn cầu, hướng đến mục tiêu nâng cao hạn dự báo các đợt mưa lớn diện rộng từ 3 ngày lên từ 10 đến 15 ngày, thực hiện lộ trình áp dụng thử nghiệm và chuẩn hóa trong giai đoạn 2026 – 2027.

Thứ hai, cập nhật quy trình vận hành liên hồ chứa trên các lưu vực sông trọng điểm miền Trung như sông Mã, sông Vu Gia – Thu Bồn và sông Ba. Cơ quan chủ trì là Bộ Tài nguyên và Môi trường phối hợp cùng Bộ Nông nghiệp và Phát triển Nông thôn cùng Tập đoàn Điện lực Việt Nam, áp dụng chỉ số cảnh báo pha 6 và 7 để chủ động hạ dung tích đón lũ trước từ 7 đến 10 ngày, hướng đến mục tiêu giảm thiểu 25% nguy cơ ngập lụt vùng hạ du và nâng cao hiệu quả tích nước phát điện trong giai đoạn 2026 – 2028.

Thứ ba, xây dựng cẩm nang chỉ dẫn ứng phó thiên tai chuyên biệt dựa trên chu kỳ nội mùa cho Ban Chỉ đạo Quốc gia về Phòng, chống Thiên tai. Giải pháp này giúp các địa phương từ Thanh Hóa đến Bình Định chủ động rà soát các điểm có nguy cơ lũ quét, sạt lở đất ngay khi nhận bản tin dao động bước vào pha kích hoạt đối lưu, đặt mục tiêu giảm thiểu ít nhất 20% thiệt hại về người và tài sản trong mùa mưa bão, hoàn thành ban hành tài liệu trong năm 2026.

Thứ tư, thiết lập cơ sở dữ liệu mở về 8 dạng hình thế synop kết hợp dao động nội mùa phục vụ nghiên cứu và đào tạo đại học. Trường Đại học Khoa học Tự nhiên phối hợp cùng Viện Khoa học Khí tượng Thủy văn và Biến đổi khí hậu chủ trì triển khai nhằm nâng cao năng lực nghiên cứu ứng dụng cho cán bộ khí tượng trẻ, thời gian thực hiện liên tục từ năm 2026 đến năm 2029.

Đối tượng nên tham khảo luận văn

Nhóm dự báo viên và chuyên viên khí tượng thủy vă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: Luận văn cung cấp công cụ định lượng và tiêu chí phân tích pha dao động để nhận diện sớm các đợt mưa diện rộng trước từ 2 đến 3 tuần, ứng dụng trực tiếp vào việc lập các bản tin dự báo khí hậu và thời tiết hạn mùa.

Nhóm nhà nghiên cứu, giảng viên và học viên cao học chuyên ngành Khí tượng, Khí hậu học và Hải dương học: Công trình là tài liệu tham khảo giá trị về phương pháp trích xuất tín hiệu nội mùa thông qua phân tích trực giao EOF và kỹ thuật xử lý số liệu lưới toàn cầu bằng ngôn ngữ lập trình khoa học NCL.

Nhóm cán bộ quản lý thiên tai và kỹ sư điều hành hồ đập thủy điện, thủy lợi: Luận văn giúp nắm vững cơ chế xuất hiện các đợt mưa cực đoan kết hợp giữa dao động nội mùa và các hình thế synop nguy hiểm, làm cơ sở xây dựng kịch bản điều tiết dòng chảy và phương án sơ tán dân cư an toàn.

Nhóm sinh viên chuyên ngành Khoa học Trái đất và Tài nguyên Nước: Tài liệu mang lại bức tranh tổng quan toàn diện và dễ hiểu về động lực học khí quyển nhiệt đới, làm tài liệu học tập và nghiên cứu đồ án tốt nghiệp chuyên ngành khí hậu học ứng dụng.

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

Dao động Madden – Julian là gì và có chu kỳ biến động như thế nào? Đây là hiện tượng dao động nội mùa quy mô lớn trong vùng nhiệt đới, kết hợp giữa hoàn lưu khí quyển và mây đối lưu sâu. Dao động di chuyển từ tây sang đông với vận tốc khoảng 5 m/s, sở hữu chu kỳ biến động bất quy tắc từ 30 đến 90 ngày, phổ biến nhất trong khoảng 40 đến 50 ngày.

Làm thế nào để xác định dao động đang ở trạng thái hoạt động mạnh? Cường độ dao động được tính toán thông qua chỉ số biên độ tổng hợp RMM từ hai thành phần chính PC1 và PC2. Khi chỉ số RMM lớn hơn 1,0, dao động được xác định là đang hoạt động mạnh; ngược lại, khi giá trị này nhỏ hơn 1,0, dao động được xem là suy yếu hoặc triệt tiêu.

Tại Việt Nam, dao động ở những pha nào dễ gây ra mưa lớn diện rộng nhất? Mưa lớn diện rộng tại Việt Nam xuất hiện với tần suất cao nhất khi dao động ở pha 4, 5, 6 và 7. Trong đó, pha 4 và 5 kích hoạt mưa mạnh ở Nam Bộ và Tây Nguyên, còn pha 6 và 7 kích hoạt các đợt mưa lũ đặc biệt lớn tại khu vực Trung Bộ.

Sự khác biệt kỹ thuật giữa chỉ số RMM của BoM và bộ chỉ số ReCal là gì? Chỉ số RMM của BoM được tính toán trực tiếp từ chuỗi số liệu thời gian thực sau khi trừ trung bình 120 ngày trước đó, trong khi chỉ số ReCal sử dụng phép lọc dải tần số 20 đến 100 ngày giúp loại bỏ hoàn toàn các dao động ngắn hạn, tạo quỹ đạo dịch chuyển mượt mà hơn.

Luận văn đã phân loại bao nhiêu nhóm hình thế synop gây mưa lớn tại Việt Nam? Nghiên cứu đã hệ thống hóa và phân loại 8 nhóm hình thế synop chính ở tầng đối lưu thấp, bao gồm: rãnh áp thấp mặt đất, dải hội tụ nhiệt đới, không khí lạnh, xoáy thuận nhiệt đới, nhiễu động gió đông, gió mùa tây nam mạnh, rãnh xích đạo, và đới gió đông nam hoạt động mạnh.

Kết luận

Nghiên cứu đã thực hiện thành công việc lượng hóa mối liên hệ giữa dao động Madden – Julian và các đợt mưa lớn diện rộng tại Việt Nam với 5 đóng góp học thuật nổi bật:

  • Xác lập hoàn chỉnh bộ chỉ số ReCal với độ tương quan cao từ 0,717 đến 0,945 so với chuẩn quốc tế của Cục Khí tượng Úc, chứng minh hiệu quả vượt trội của phương pháp lọc dải chu kỳ 20 đến 100 ngày.
  • Phát hiện và chứng minh quy luật phân bố dị thường của mưa lớn tập trung ở các pha 4, 5, 6 và 7, giải thích rõ cơ chế tương tác giữa sóng nội mùa quy mô hành tinh với địa hình khu vực.
  • Thống kê chi tiết diễn biến mưa lớn diện rộng trong 14 năm từ năm 2000 đến năm 2013, chỉ ra các đợt mưa cực đoan có MJO mạnh thường xuất hiện trong giai đoạn từ tháng 5 đến tháng 12.
  • Phân loại chuẩn xác 8 nhóm hình thế thời tiết synop tầng đối lưu thấp, làm sáng tỏ cơ chế kết hợp giữa xoáy thuận nhiệt đới, không khí lạnh và dao động nội mùa gây mưa lớn tại miền Trung.
  • Đặt nền móng khoa học vững chắc cho việc ứng dụng chỉ số dao động nội mùa vào quy trình nghiệp vụ dự báo khí tượng hạn vừa và mở rộng tại Việt Nam.

Hướng nghiên cứu tiếp theo trong giai đoạn 2026 – 2030 là mở rộng đánh giá tương tác phi tuyến giữa dao động Madden – Julian và các pha nóng lạnh của hiện tượng ENSO đến cấu trúc mưa cực đoan cục bộ. Các cơ quan quản lý và dự báo viên khí tượng thủy văn được khuyến khích áp dụng ngay kết quả nghiên cứu này vào thực tiễn tác nghiệp để nâng cao hiệu quả phòng chống thiên tai quốc gia.