Tổng quan nghiên cứu

Khu vực Biển Đông hàng năm đón nhận trung bình từ 10 đến 12 cơn bão và áp thấp nhiệt đới, gây thiệt hại nghiêm trọng về người cùng tổn thất kinh tế ước tính chiếm từ 1% đến 1,5% GDP của Việt Nam. Công tác dự báo quỹ đạo và cường độ bão đóng vai trò then chốt trong việc giảm thiểu rủi ro thiên tai. Tuy nhiên, phần lớn các mô hình số trị nghiệp vụ truyền thống tại Việt Nam như MM5, ETA, RAMS hay WRF đơn thuần chỉ giải quyết các quá trình động lực khí quyển mà bỏ qua sự phản hồi hai chiều từ bề mặt đại dương.

Mục tiêu cốt lõi của nghiên cứu là đánh giá định lượng tác động của tương tác đại dương - khí quyển đến quá trình biến đổi cường độ và quỹ đạo bão tại khu vực Biển Đông bằng cách thiết lập mô hình ghép nối phân giải cao. Nghiên cứu thực hiện chuyển đổi và tùy biến hệ thống mô hình HWRF (vốn sử dụng mô hình POM tại khu vực Đại Tây Dương) sang sự kết hợp giữa mô hình khí quyển WRF và mô hình đại dương ROMS thông qua khung liên kết COAWST.

Phạm vi không gian của nghiên cứu bao phủ toàn bộ vùng Biển Đông và bán đảo Đông Dương trong tọa độ từ 3°N đến 30°N và 90°E đến 130°E. Phạm vi thời gian tập trung mô phỏng chi tiết 3 cơn bão nhiệt đới tiêu biểu gồm Mindulle (năm 2010), Nock-ten (năm 2011) và Nalgae (năm 2011). Kết quả nghiên cứu có ý nghĩa khoa học và thực tiễn sâu sắc, cung cấp cơ sở định lượng giúp giảm sai số dự báo áp suất tâm bão từ 14 đến 24 mb ở các hạn dự báo từ 48 đến 72 giờ, nâng cao độ tin cậy cảnh báo thiên tai cho các vùng ven biển.

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 nền tảng lý thuyết tương tác động lực và nhiệt dịch học giữa lớp biên khí quyển và lớp hòa lẫn trên mặt biển. Khi bão di chuyển trên biển, trường gió mạnh truyền động lượng xuống bề mặt nước, gây ra hiện tượng xáo trộn rối thẳng đứng và hiện tượng nước trồi (upwelling), làm đưa khối nước lạnh từ tầng sâu lên làm mát bề mặt biển (SST). Sự suy giảm nhiệt độ bề mặt biển này lập tức tạo ra phản hồi âm (negative feedback), làm cắt giảm thông lượng nhiệt hiện (HFX) và thông lượng ẩm/ẩn nhiệt (QFX, LH) bốc hơi từ biển vào hoàn lưu bão, trực tiếp kìm hãm cường độ bão phát triển.

Khung mô hình ghép nối vận dụng hai hệ phương trình thủy động lực học quy mô vừa:

  • Mô hình khí quyển WRF (phiên bản ARW): Giải hệ phương trình phi thủy tĩnh nén được đầy đủ trên lưới so le Arakawa C theo phương ngang và hệ tọa độ khối theo địa hình với 27 mực thẳng đứng trên cao.
  • Mô hình đại dương khu vực ROMS: Giải hệ phương trình thủy động lực thủy tĩnh và mặt nước tự do trên hệ tọa độ cong trực giao theo phương ngang và hệ tọa độ s thích ứng địa hình theo phương thẳng đứng với 20 tầng lớp nước.

Các khái niệm then chốt trong nghiên cứu bao gồm: Nhiệt độ bề mặt nước biển (SST), Khí áp mực biển cực tiểu (SLP min), Thông lượng nhiệt ẩn bề mặt (LH), Sai số dọc quỹ đạo (ATE) và Sai số ngang quỹ đạo (CTE).

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

Nghiên cứu lựa chọn cỡ mẫu gồm 3 cơn bão nhiệt đới đặc trưng hoạt động trên Biển Đông trong giai đoạn 2010 - 2011 với 18 mốc thời gian đánh giá từ 00h đến 72h. Phương pháp chọn mẫu tập trung vào các cơn bão có diễn biến phức tạp, cường độ biến đổi nhanh từ cấp 10 đến cấp 15 và đổ bộ trực tiếp vào khu vực Bắc Bộ và Trung Bộ Việt Nam.

Nguồn dữ liệu đầu vào bao gồm:

  • Trường khí tượng phân tích và dự báo toàn cầu GFS từ Trung tâm Dự báo Môi trường Quốc gia Hoa Kỳ (NCEP) với độ phân giải không gian phục vụ điều kiện ban đầu và điều kiện biên cho WRF.
  • Trường nhiệt độ, độ mặn và mực nước biển trung bình tháng từ bộ số liệu khí hậu COADS kết hợp trường SST cập nhật tuần của NCEP cho mô hình ROMS.
  • Địa hình đáy biển toàn cầu ETOPO-2 với độ phân giải 2 phút góc (xấp xỉ 4 km) dùng để nội suy địa hình đáy trong ROMS.

Phương pháp phân tích dựa trên kỹ thuật tích phân số trị ghép nối trực tiếp hai chiều (online two-way coupling) qua công cụ MCT trong hệ thống COAWST. Bước thời gian trao đổi thông tin giữa WRF và ROMS được thiết lập là 10 phút (600 giây), bước tích phân của WRF là 120 giây và ROMS là 150 giây, vận hành song song trên cụm máy tính bó 8 nodes mạng MPI. Lý do lựa chọn phương pháp ghép nối này là nhằm tái lập chính xác phản ứng nhiệt động học tức thời của mặt biển dưới tác động của bão mà các mô hình chạy đơn lẻ không thể thực hiện được.

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

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

Quá trình phân tích thực nghiệm so sánh giữa phương án ghép nối WRF-ROMS và phương án khí quyển đơn thuần WRF mang lại những phát hiện quan trọng:

Thứ nhất, mô hình ghép nối tái hiện chân thực hiệu ứng làm lạnh bề mặt biển (SST cooling) dưới tâm bão. Tại thời điểm ban đầu (00h), SST trong mô hình ROMS thấp hơn từ 2°C đến 3°C so với trường SST tĩnh của WRF. Khi bão phát triển và di chuyển, độ chênh lệch nhiệt độ bề mặt biển tăng nhanh lên mức 3°C đến 4°C ở hạn 12h và 24h, đạt mức chênh lệch 4°C đến 5°C (cá biệt trên 5°C tại khu vực xung quanh đảo Hải Nam trong bão Nock-ten) ở hạn dự báo 48h và 72h.

Thứ hai, việc suy giảm SST tác động trực tiếp đến việc điều chỉnh cường độ khí áp bão. Trong phương án đơn thuần WRF, do SST được giữ cố định theo thời gian nên nguồn năng lượng nhiệt ẩm liên tục bơm vào hoàn lưu bão, khiến mô hình dự báo bão sâu hơn và mạnh hơn rất nhiều so với thực tế. Ngược lại, phương án WRF-ROMS kìm hãm cường độ bão: trị số khí áp cực tiểu tại tâm (SLP min) của WRF-ROMS cao hơn phương án WRF từ 2 đến 3 mb ở hạn 12h, tăng lên 6 đến 8 mb ở hạn 24h, và chênh lệch lớn từ 14 đến 16 mb ở hạn 48h, đạt cực đại 22 đến 24 mb ở hạn dự báo 72h.

Thứ ba, tương tác biển - khí quyển góp phần cải thiện độ chính xác dự báo quỹ đạo bão ở các hạn dự báo dài. Sai số khoảng cách vị trí tâm bão trung bình (MPE) của phương án WRF-ROMS giảm khoảng 15% đến 25% ở các hạn 48h và 72h so với mô hình WRF đơn lẻ. Sai số ngang (CTE) và sai số dọc (ATE) cho thấy xu thế di chuyển của bão trong phương án ghép nối có độ ổn định cao hơn, hạn chế hiện tượng bão di chuyển quá nhanh do hoàn lưu xoáy quá mạnh gây ra.

Thảo luận kết quả

Cơ chế vật lý đằng sau sự khác biệt giữa hai phương án bắt nguồn từ chu trình trao đổi thông lượng tại lớp biên. Dữ liệu mô phỏng cho thấy khi nhiệt độ mặt biển giảm từ 2°C đến 4°C, thông lượng nhiệt hiện (HFX) và thông lượng ẩn nhiệt (LH) trung bình trên một ô lưới tâm bão trong WRF-ROMS giảm từ 30% đến 45% so với mô hình WRF đơn lẻ. Profile thẳng đứng của nhiệt độ và gió nhiễu động tại tâm bão cho thấy cấu trúc lõi ấm (warm core) của WRF-ROMS phát triển ở mức độ vừa phải, phù hợp với số liệu thám sát thực tế của bão Mindulle, Nock-ten và Nalgae.

Kết quả này hoàn toàn tương đồng với các công trình nghiên cứu quốc tế của Chun-Chieh Wu và cộng sự khi khẳng định thành phần làm lạnh SST tạo ra phản hồi âm kìm hãm bão phát triển cực đoan. Đồng thời, nghiên cứu kế thừa và chứng minh bằng định lượng kết luận của các nhà khoa học khí tượng Việt Nam về hiện tượng "vết nước lạnh" sau bão và tác động tiêu tán năng lượng của các vùng nước trồi ven bờ miền Trung.

Dữ liệu mô phỏng trong nghiên cứu được thể hiện trực quan qua bản đồ phân bố không gian hiệu nhiệt độ bề mặt biển (SST Difference Map), biểu đồ biến trình thời gian của khí áp tâm bão SLP min và biểu đồ phân rã sai số quỹ đạo ATE/CTE theo từng hạn dự báo từ 12h đến 72h, giúp minh chứng rõ nét sự vượt trội của hệ thống mô hình ghép nối.

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

Dựa trên kết quả nghiên cứu thực nghiệm, các giải pháp và khuyến nghị kỹ thuật được đề xuất như sau:

  • Tích hợp mô hình ghép nối WRF-ROMS vào hệ thống nghiệp vụ dự báo: Trung tâm Dự báo Khí tượng Thủy văn Quốc gia cần chủ trì triển khai cài đặt mô hình ghép nối hai chiều trên hệ thống siêu máy tính nghiệp vụ, đặt mục tiêu khống chế sai số khí áp tâm bão dưới 5 mb ở hạn dự báo 24 giờ và dưới 10 mb ở hạn 48 giờ trong lộ trình từ năm 2024 đến 2026.
  • Tăng độ phân giải không gian miền tính đại dương và khí quyển: Viện Khoa học Khí tượng Thủy văn và Biến đổi khí hậu cần thực hiện tinh chỉnh bước lưới tính từ 27,75 km xuống mức 3 km đến 9 km cho khu vực trọng điểm ven bờ Biển Đông trước năm 2026 nhằm mô phỏng chính xác tương tác địa hình và các hoàn lưu xoáy quy mô nhỏ.
  • Nâng cấp mạng lưới đồng hóa dữ liệu biển thời gian thực: Cục Khí tượng Thủy văn cần bổ sung từ 15 đến 20 trạm phao đại dương đo sâu tự động và tích hợp dữ liệu nhiệt - mặn vệ tinh vào mô-đun ROMSTOOL và WRFDA, mục tiêu nâng độ chính xác vị trí tâm bão thêm 10% đến 15% trong giai đoạn 2025 - 2027.
  • Phát triển hệ thống dự báo tích hợp đa thành phần (COAWST đầy đủ): Các nhóm nghiên cứu chuyên sâu tại Đại học Quốc gia Hà Nội cần mở rộng ghép nối mô hình sóng SWAN và mô hình vận chuyển bùn cát CSTMS với WRF-ROMS, nhằm đưa ra dự báo đồng thời về bão, nước dâng và sóng lớn với thời gian xuất kết quả dưới 30 phút sau khi nhận dữ liệu biên toàn cầu.

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

  • Chuyên viên dự báo khí tượng thủy văn tại các đài khu vực: Tiếp cận quy trình vận hành hệ thống mô hình tương tác biển - khí quyển để ứng dụng trực tiếp vào công tác phân tích, nhận định quỹ đạo và cường độ bão hạn 24 - 72 giờ trên Biển Đông.
  • Học viên cao học, nghiên cứu sinh ngành Khí tượng và Hải dương học: Sử dụng tài liệu như một tài liệu tham khảo chuẩn mực về cấu hình mô hình WRF, mô hình ROMS, kỹ thuật ghép nối công cụ MCT trong COAWST và phương pháp thiết lập lưới s 20 tầng thích ứng địa hình đáy biển.
  • Cơ quan quản lý và Ban Chỉ đạo Quốc gia về Phòng chống thiên tai: Khai thác các đánh giá khoa học để xây dựng kịch bản ứng phó bão sát thực tế, tránh các sai lệch về cường độ cực đoan gây lãng phí nguồn lực sơ tán dân cư.
  • Kỹ sư công trình biển, dầu khí và các đơn vị vận tải hàng hải: Tham khảo các phân tích về tương tác trường gió, dòng chảy mặt và biến động khí áp để tối ưu hóa thiết kế giàn khoan và lập lộ trình hàng hải an toàn trong mùa bão.

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

Tại sao mô hình khí quyển WRF đơn lẻ thường dự báo cường độ bão mạnh hơn thực tế?
Mô hình WRF đơn lẻ sử dụng trường SST tĩnh cố định trong suốt quá trình dự báo, bỏ qua quá trình nước trồi lạnh do gió bão xới trộn mặt biển. Do đó, nguồn ẩn nhiệt từ mặt biển liên tục cấp năng lượng cho bão, khiến khí áp tâm bão bị tính toán quá sâu (sai lệch từ 14 đến 24 mb ở hạn 48h - 72h).

Điểm ưu việt khi thay thế mô hình đại dương POM bằng ROMS trong hệ thống HWRF là gì?
Mô hình ROMS sử dụng thuật toán số trị bậc cao, hệ tọa độ s thích ứng địa hình đáy biển phức tạp và cơ chế tham số hóa xáo trộn rối Mellor-Yamada 2.5 cải tiến. Cấu trúc mã nguồn của ROMS cho phép ghép nối tối ưu hai chiều với WRF qua bộ công cụ MCT, phù hợp hoàn hảo với địa hình độ sâu biến đổi từ 10m đến 4500m của Biển Đông.

Tương tác đại dương - khí quyển ảnh hưởng thế nào đến quỹ đạo bão?
Sự phân bố bất đối xứng của trường làm lạnh SST làm thay đổi cấu trúc hoàn lưu đối lưu xung quanh tâm bão. Khi cường độ bão được kìm hãm hợp lý trong mô hình WRF-ROMS, lực dẫn môi trường tầng sâu tác động chuẩn xác hơn, giúp giảm sai số khoảng cách tâm bão (MPE) từ 15% đến 25% ở hạn dự báo 48h - 72h.

Miền tính và độ phân giải nào được thiết lập trong nghiên cứu này?
Nghiên cứu thiết lập miền tính kích thước 110x160 điểm lưới, độ phân giải 1/4 độ kinh vĩ (khoảng 27,75 km), bao phủ từ 3°N đến 30°N và 90°E đến 130°E. WRF cấu hình 27 mực thẳng đứng, còn ROMS chia thành 20 lớp s theo phương thẳng đứng kết hợp địa hình đáy ETOPO-2 độ phân giải 2 phút góc.

Bước thời gian ghép nối 10 phút giữa WRF và ROMS có vai trò gì?
Khoảng thời gian 600 giây (10 phút) là chu kỳ trao đổi dữ liệu tối ưu giữa hai mô hình: WRF gửi ứng suất gió, khí áp, bức xạ, mưa xuống ROMS; ROMS tính toán và phản hồi trường SST tức thời lên WRF. Bước thời gian này vừa đảm bảo tính ổn định số trị theo chuẩn CFL, vừa phản ánh kịp thời phản hồi nhiệt động lực học của mặt biển.

Kết luận

  • Luận văn đã thiết lập thành công hệ thống mô hình ghép nối hai chiều WRF-ROMS thông qua khung COAWST cho khu vực Biển Đông và bán đảo Đông Dương.
  • Khẳng định tương tác đại dương - khí quyển tạo ra hiệu ứng phản hồi âm làm giảm nhiệt độ bề mặt biển từ 2°C đến 5°C, trực tiếp kìm hãm cường độ bão không bị phát triển quá mức so với thực tế.
  • Khí áp tâm bão dự báo bằng WRF-ROMS sát với số liệu thực tế hơn, khắc phục sai số khí áp quá sâu từ 14 đến 24 mb tại các hạn dự báo từ 48 đến 72 giờ của mô hình WRF đơn lẻ.
  • Cải thiện độ chính xác dự báo quỹ đạo bão Mindulle, Nock-ten và Nalgae, giảm sai số khoảng cách vị trí tâm từ 15% đến 25% ở các hạn dự báo xa.
  • Đặt nền móng khoa học vững chắc cho việc chuyển giao công nghệ mô hình ghép nối số trị vào nghiệp vụ dự báo bão hạn 3 ngày tại Việt Nam.

Giai đoạn từ năm 2025 đến 2030 đòi hỏi việc tiếp tục nâng cấp hệ thống mô hình tích hợp sóng và nước dâng. Hãy ứng dụng ngay các kết quả nghiên cứu của luận văn để nâng cao năng lực phân tích số trị và tối ưu hóa hệ thống cảnh báo thiên tai khí tượng thủy văn.