Tổng quan nghiên cứu

Khu vực Biển Đông nằm trong khoảng tọa độ 5°N - 25°N và 105°E - 120°E, là một trong những cái nôi hoạt động bão phức tạp nhất trên thế giới. Theo thống kê khí hậu nhiều năm, trung bình mỗi năm Biển Đông đón nhận khoảng 11,9 đến 12,0 cơn xoáy thuận nhiệt đới (XTNĐ), trong đó hơn 68% tổng số cơn bão và áp thấp nhiệt đới tập trung dồn dập vào các tháng từ tháng 6 đến tháng 11. Các cơn bão kèm theo gió giật mạnh và lượng mưa lớn luôn tiềm ẩn nguy cơ gây thiệt hại nghiêm trọng về tính mạng con người, cơ sở hạ tầng ven biển và các hoạt động kinh tế biển của Việt Nam.

Vấn đề nghiên cứu trọng tâm của luận văn là giải quyết bài toán dự báo hạn mùa (dự báo trước từ 1 đến 3 tháng) về tổng số lượng XTNĐ hoạt động trên Biển Đông bằng các công cụ thống kê khí hậu hiện đại. Mục tiêu cụ thể là xây dựng các phương trình dự báo thống kê định lượng có độ tin cậy cao, xác định mối liên hệ phi tuyến giữa các hoàn lưu khí quyển - đại dương quy mô lớn với tần suất bão. Phạm vi nghiên cứu sử dụng chuỗi quan trắc 55 năm (1961 - 2015) để phân tích quy luật khí hậu và tập trung xây dựng mô hình trên chuỗi số liệu 30 năm (1981 - 2010), sau đó kiểm nghiệm độc lập trên chuỗi số liệu 5 năm (2011 - 2015).

Nghiên cứu có ý nghĩa thực tiễn to lớn trong việc nâng cao chỉ số kỹ năng dự báo mùa, giảm thiểu sai số tuyệt đối trung bình (MAE) và sai số quân phương (RMSE) trong công tác nghiệp vụ. Kết quả nghiên cứu cung cấp căn cứ khoa học định lượng giúp các cơ quan chỉ đạo phòng chống thiên tai tại 28 tỉnh, thành phố ven biển chủ động xây dựng kế hoạch ứng phó, bảo vệ an toàn cho hàng vạn tàu thuyền và hàng triệu cư dân ven bờ.

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 thiết lập dựa trên 2 nền tảng lý thuyết cốt lõi trong khoa học khí quyển: lý thuyết tương tác biển - khí quyển (Atmosphere-Ocean Interaction) và lý thuyết hoàn lưu khí quyển quy mô hành tinh. Các hệ thống thời tiết nhiệt đới chịu sự điều khiển trực tiếp bởi các dao động quy mô lớn, trong đó năng lượng nhiệt ẩm từ bề mặt đại dương là nguồn cấp động lực chính cho sự hình thành và duy trì xoáy thuận nhiệt đới.

Các khái niệm chuyên ngành then chốt bao gồm:

  1. Xoáy thuận nhiệt đới (XTNĐ): Hệ thống áp thấp quy mô lớn hình thành trên vùng biển nhiệt đới với hoàn lưu xoáy khép kín, được phân cấp theo tiêu chuẩn của Tổ chức Khí tượng Thế giới (WMO) thành áp thấp nhiệt đới (vận tốc gió cực đại 39 - 62 km/h), bão thường (62 - 88 km/h), bão mạnh (89 - 117 km/h) và bão rất mạnh (từ 118 km/h trở lên).
  2. Hiện tượng ENSO (El Niño - Southern Oscillation): Dao động đại dương - khí quyển nhiệt đới Thái Bình Dương với 3 pha: El Niño (pha nóng), La Niña (pha lạnh) và Trung tính (Neutral), được xác định dựa trên ngưỡng chuẩn sai nhiệt độ mặt nước biển (SSTA) vùng Nino 3.4 lệch ±0,5°C trong 6 tháng mùa bão.
  3. Độ đứt gió thẳng đứng (Vertical Wind Shear - WS500-850): Độ chênh lệch vectơ vận tốc gió giữa mực 500 hPa và 850 hPa, yếu tố quyết định đến sự tích tụ năng lượng xoáy đối lưu.
  4. Mô hình hồi quy độ lệch tuyệt đối nhỏ nhất (Least Absolute Deviation - LAD): Phương pháp ước lượng tham số tối ưu hóa tổng giá trị tuyệt đối của sai số, giảm thiểu tối đa ảnh hưởng tiêu cực của các giá trị ngoại lai (outliers).

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

Nghiên cứu sử dụng 3 nguồn dữ liệu lớn: bộ số liệu quỹ đạo bão Unisys (giai đoạn 1961 - 2015), bộ số liệu tái phân tích khí quyển - đại dương CFSR từ NCEP/NOAA với lưới phân giải 0,5° x 0,5° (giai đoạn 1981 - 2015), và 12 chỉ số dao động khí hậu toàn cầu (NINO3.4, SOI, ONI, MEI, PDO, QBO, PNA, WP...). Cỡ mẫu nghiên cứu gồm 30 năm số liệu phụ thuộc (1981 - 2010) dùng để huấn luyện mô hình và 5 năm số liệu độc lập (2011 - 2015) dùng để kiểm chứng.

Quy trình phân tích gồm các bước chặt chẽ:

  1. Phương pháp phân tích tương quan không gian: Lập bản đồ tương quan giữa số lượng XTNĐ mùa bão (tháng 6 - 11) với các trường khí tượng tháng 1, 2, 3 (nhiệt độ mặt biển SST, gió vĩ hướng U500, U850, áp suất mực biển PRMSL, độ cao địa thế vị HGT500, độ đứt gió WS500-850) để khoanh vùng các khu vực có tín hiệu vật lý rõ nét.
  2. Phân tích thành phần chính (PCA): Giảm số chiều dữ liệu không gian nhiều chiều về các vector trực giao (PCs), giữ lại trên 75% tổng phương sai giải thích mà không làm mất thông tin khí hậu nền.
  3. Hồi quy từng bước (Stepwise Regression): Tự động sàng lọc các nhân tố dự báo tối ưu, loại bỏ hiện tượng cộng tuyến đa biến.
  4. Xây dựng và kiểm định mô hình: Thiết lập phương trình hồi quy tuyến tính (HQTT) và hồi quy LAD cho từng hạn dự báo (1, 2 và 3 tháng). Đánh giá sai số thông qua 4 chỉ số thống kê tiêu chuẩn: sai số trung bình (ME), sai số tuyệt đối trung bình (MAE), sai số quân phương (RMSE) và hệ số tương quan (CORR).

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 đã lượng hóa chi tiết tác động của chu kỳ ENSO đến hoạt động của XTNĐ trên Biển Đông trong 55 năm (1961 - 2015). Kết quả thống kê chỉ ra có 14 năm El Niño (chiếm 25,5%), 13 năm La Niña (chiếm 23,6%) và 28 năm Trung tính (chiếm 50,9%). Tần suất XTNĐ biến động mạnh theo từng pha: những năm La Niña ghi nhận trung bình 13,8 cơn/năm (chiếm tỷ trọng 45% - 48% tổng số bão khu vực Tây Bắc Thái Bình Dương), cao hơn rõ rệt so với những năm Trung tính (11,86 cơn/năm) và thấp nhất trong các năm El Niño với chỉ khoảng 10,4 cơn/năm (chênh lệch tới hơn 32%).

Thứ hai, cấu trúc cường độ bão trên Biển Đông cho thấy tỷ lệ bão rất mạnh chiếm ưu thế lớn nhất trong toàn bộ chuỗi quan trắc với tỷ lệ trung bình đạt 47,3%, tiếp theo là bão thường và bão mạnh, trong khi áp thấp nhiệt đới chỉ chiếm tỷ lệ thấp nhất là 14,3%. Đáng chú ý, trong các năm El Niño, tỷ lệ bão rất mạnh đạt mức cực đại lên tới 48,9%, cao hơn năm La Niña (47,1%) và năm Trung tính (44,4%).

Thứ ba, việc trích xuất nhân tố dự báo bằng PCA và Stepwise Regression đã xác lập được tập nhân tố tiền đề có ý nghĩa vật lý cao. Đối với hạn dự báo 1 tháng, các trường khí quyển tháng 3 như áp suất mực biển khu vực Tây Bắc Thái Bình Dương (20°N - 40°N, 150°E - 180°E) và gió vĩ hướng U500 vùng vĩ độ thấp (5°N - 10°N, 160°E - 175°E) đóng góp tới hơn 70% phương sai giải thích của phương trình dự báo.

Thứ tư, khi thử nghiệm trên tập số liệu kiểm tra độc lập 2011 - 2015, phương pháp hồi quy LAD thể hiện sự vượt trội về độ ổn định so với hồi quy tuyến tính cổ điển. Mô hình LAD duy trì hệ số tương quan CORR đạt trên 0,72 đối với hạn 1 tháng và giữ mức sai số MAE dao động trong ngưỡng 1,1 đến 1,4 cơn/mùa, hạn chế tối đa các sai số đột biến.

Thảo luận kết quả

Cơ chế vật lý giải thích cho sự gia tăng bão trong pha La Niña là do sự dịch chuyển của hoàn lưu Walker, làm tăng cường dòng xoáy đối lưu nhiệt đới và giảm độ đứt gió thẳng đứng ở khu vực Tây Thái Bình Dương và Biển Đông, tạo điều kiện lý tưởng cho phôi bão hình thành. Ngược lại, trong các năm El Niño, dòng xiết gió tây tầng trên hoạt động mạnh làm tăng đứt gió thẳng đứng (WS500-850), ức chế quá trình hội tụ ẩm tầng thấp khiến số lượng bão giảm đi, tuy nhiên nền nhiệt đại dương tích tụ lại dễ tạo ra các cơn bão có cường độ cực mạnh.

Khi so sánh với các nghiên cứu tiền tiêu của các tác giả quốc tế (Camargo et al., 2007; tác giả nghiên cứu vùng biển Đông Á năm 2010), việc kết hợp lọc nhân tố PCA với mô hình LAD giúp nâng cao kỹ năng dự báo lên khoảng 15% đến 20% so với phương pháp dự báo dựa trên chuỗi khí hậu tham chiếu.

Dữ liệu kết quả khi biểu diễn trên đồ thị phân tán giữa số lượng bão thực tế và số lượng bão dự báo cho thấy đường hồi quy LAD bám rất sát đường chuẩn 45 độ, các điểm biểu diễn số liệu độc lập giai đoạn 2011 - 2015 phân bố tập trung với biên độ dao động nhỏ, minh chứng cho tính ứng dụng thực tiễn cao của mô hình trong công tác khí tượng nghiệp vụ.

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

  1. Tích hợp mô hình dự báo LAD vào hệ thống tác nghiệp chính thức: Tổng cục Khí tượng Thủy văn cần triển khai cài đặt phương trình dự báo hạn mùa bằng phương pháp LAD vào quy trình nghiệp vụ thường quy trước ngày 15 tháng 5 hàng năm. Đặt mục tiêu nâng cao độ chính xác dự báo số lượng bão mùa sớm thêm 10% - 15% trong lộ trình giai đoạn 2026 - 2028.
  2. Mở rộng và cập nhật liên tục bộ chỉ số tiền đề đại dương: Trung tâm Dự báo Khí tượng Thủy văn Quốc gia cần định kỳ 6 tháng một lần bổ sung các chỉ số khí hậu thế hệ mới (như ENSO Modoki, IOD, dao động MJO) kết hợp dữ liệu tái phân tích ERA5 và CFSv2 độ phân giải cao 0,25°, nhằm kiểm soát sai số quân phương RMSE ở mức dưới 1,2 cơn/mùa.
  3. Xây dựng phần mềm tự động hóa chuỗi xử lý PCA và Stepwise: Đơn vị nghiên cứu công nghệ thông tin ngành tài nguyên môi trường cần phát triển công cụ tự động trích xuất dữ liệu số trị dạng lưới, thực hiện phân tích thành phần chính và chạy mô hình dự báo tự động, giúp rút ngắn thời gian tính toán từ 48 giờ xuống dưới 2 giờ cho các chuyên viên dự báo.
  4. Tăng cường cơ chế chia sẻ thông tin dự báo sớm cho công tác phòng chống thiên tai: Ban Chỉ đạo Quốc gia về Phòng, chống thiên tai và Ủy ban chỉ huy cấp tỉnh tại 28 địa phương ven biển cần sử dụng bản tin dự báo hạn mùa trước 3 tháng (tháng 3 hàng năm) để xây dựng phương án di dời dân cư vùng trũng thấp, điều phối lịch thời vụ nuôi trồng thủy sản và bố trí lực lượng cứu hộ cứu nạn thường trực trước ngày 1 tháng 6 hàng năm.

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

  1. Học viên cao học, nghiên cứu sinh và sinh viên ngành Khí tượng, Khí hậu học, Thủy văn: Luận văn là tài liệu tham khảo chuẩn mực về phương pháp luận xử lý số liệu lưới khí tượng nhiều chiều, kỹ thuật phân tích thành phần chính PCA và cách thiết lập mô hình hồi quy độ lệch tuyệt đối nhỏ nhất (LAD).
  2. Dự báo viên và chuyên viên phân tích tại các Đài Khí tượng Thủy văn khu vực: Nghiên cứu cung cấp tài liệu hướng dẫn quy trình thực tế từ việc bóc tách số liệu tái phân tích CFSR đến bước đánh giá sai số ME, MAE, RMSE phục vụ công tác xây dựng bản tin dự báo khí hậu mùa.
  3. Cán bộ quản lý tại các cơ quan phòng chống thiên tai và tìm kiếm cứu nạn: Nắm bắt được quy luật phân bố tần suất bão, xu thế cực đoan theo từng pha El Niño/La Niña để xây dựng chiến lược ứng phó rủi ro thiên tai trung và dài hạn cho từng vùng kinh tế trọng điểm.
  4. Doanh nghiệp vận tải hàng hải, khai thác cảng biển, dầu khí và nuôi trồng thủy sản xa bờ: Ứng dụng kết quả dự báo bão trước 1 đến 3 tháng để lập kế hoạch điều động tàu thuyền, gia cố công trình biển và bảo hiểm rủi ro sản xuất an toàn trong các tháng cao điểm từ tháng 6 đến tháng 11.

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

1. Phương pháp hồi quy LAD có ưu thế vượt trội nào so với phương pháp hồi quy tuyến tính bình phương tối thiểu (HQTT)?
Hồi quy bình phương tối thiểu (HQTT) tối thiểu hóa tổng bình phương sai số nên rất nhạy cảm với các năm có số lượng bão bất thường (giá trị ngoại lai), dễ dẫn đến sai số mô hình bị lệch lớn. Phương pháp LAD tối thiểu hóa tổng độ lệch tuyệt đối, giúp đường hồi quy ổn định hơn, giảm hiện tượng quá khớp (overfitting) và duy trì sai số tuyệt đối MAE thấp trên tập kiểm tra độc lập.

2. Tại sao nghiên cứu lại tập trung dự báo tổng lượng bão trong khoảng thời gian từ tháng 6 đến tháng 11?
Số liệu khí hậu giai đoạn 1961 - 2015 chỉ ra rằng hơn 68% tổng số XTNĐ hoạt động trên Biển Đông tập trung vào các tháng từ tháng 6 đến tháng 11, với đỉnh điểm rơi vào các tháng 7, 8, 9 và 10. Đây là giai đoạn gió mùa mùa hè và dải hội tụ nhiệt đới hoạt động mạnh nhất, gây ra trên 85% tổng thiệt hại do thiên tai bão lũ tại Việt Nam.

3. Hiện tượng ENSO tác động cụ thể ra sao đến số lượng và cường độ bão trên Biển Đông?
Trong pha La Niña, số lượng XTNĐ trên Biển Đông tăng cao rõ rệt, đạt trung bình 13,8 cơn/năm do điều kiện nhiệt ẩm và đối lưu thuận lợi. Ngược lại, pha El Niño làm giảm số lượng bão xuống trung bình 10,4 cơn/năm nhưng lại làm tăng tỷ lệ bão rất mạnh (vận tốc gió trên 118 km/h) lên mức cao nhất, chiếm tới 48,9% tổng số bão.

4. Vai trò then chốt của phương pháp phân tích thành phần chính (PCA) trong luận văn là gì?
Các trường khí quyển dạng lưới như SST hay gió U500 có hàng nghìn điểm nút không gian, gây hiện tượng dư thừa dữ liệu và đa cộng tuyến. PCA biến đổi không gian dữ liệu nhiều chiều thành một số ít các thành phần chính (PCs) trực giao không tương quan với nhau, giữ lại hơn 75% thông tin quan trọng nhất để đưa vào phương trình dự báo một cách tối ưu.

5. Ý nghĩa thực tiễn của hạn dự báo từ 1 đến 3 tháng trong công tác phòng chống thiên tai là gì?
Dự báo hạn 1 đến 3 tháng (thực hiện từ tháng 3, 4, 5 cho mùa bão tháng 6 - 11) cho phép các cấp chính quyền và doanh nghiệp có đủ thời gian chuẩn bị nguồn lực, điều chỉnh kế hoạch sản xuất kinh doanh, kiểm tra an toàn hệ thống đê điều hồ đập và triển khai phương án sơ tán dân cư trước khi bão đổ bộ.

Kết luận

  • Luận văn đã phân tích toàn diện quy luật hoạt động của xoáy thuận nhiệt đới trên Biển Đông giai đoạn 1961 - 2015, xác định trung bình hàng năm có 11,9 - 12,0 cơn XTNĐ với tỷ lệ bão rất mạnh chiếm 47,3%.
  • Xác lập rõ rệt mối quan hệ giữa các pha ENSO và hoạt động của bão: năm La Niña ghi nhận tần suất bão cao nhất (13,8 cơn/năm), trong khi năm El Niño ghi nhận số lượng thấp nhất (10,4 cơn/năm) nhưng tỷ lệ bão rất mạnh đạt đỉnh (48,9%).
  • Ứng dụng thành công kỹ thuật lọc nhân tố kết hợp giữa phân tích thành phần chính (PCA) và hồi quy từng bước, trích xuất chính xác các tín hiệu tiền đề từ trường SST, gió U500, U850 và áp suất PRMSL.
  • Chứng minh phương pháp hồi quy độ lệch tuyệt đối nhỏ nhất (LAD) có độ tin cậy và tính ổn định cao hơn hồi quy tuyến tính cổ điển trên tập mẫu độc lập, mở ra hướng ứng dụng hiệu quả cho bài toán dự báo mùa.
  • Các đơn vị khí tượng nghiệp vụ và cơ quan quản lý thiên tai cần sớm đưa mô hình LAD vào thử nghiệm tác nghiệp định kỳ trước mùa bão năm 2026, đồng thời tiếp tục nghiên cứu tích hợp thêm các công cụ học máy (Machine Learning) để tối ưu hóa độ chính xác dự báo trong tương lai.