Tổng quan nghiên cứu

Biến đổi khí hậu toàn cầu và sự gia tăng cường độ của các hiện tượng cực đoan như ENSO đang đặt ra những thách thức nghiêm trọng đối với công tác dự báo khí tượng thủy văn tại khu vực Đông Nam Á. Số liệu thống kê trong giai đoạn 1961 đến 2007 cho thấy tần số bão hoạt động trên khu vực Biển Đông có sự biến động phức tạp từ mức tối thiểu 6 cơn vào năm 1997 lên tới khoảng 15 cơn mỗi năm vào các năm 1996, 1999 và 2001. Tại Việt Nam, công tác dự báo hạn mùa trước đây chủ yếu dựa trên các phương pháp thống kê kinh nghiệm, vốn bộc lộ nhiều hạn chế do không phản ánh đầy đủ các cơ chế vật lý phi tuyến và gặp khó khăn khi xuất hiện các pha chuyển tiếp khí hậu đột biến.

Nhằm giải quyết bài toán cấp thiết này, nghiên cứu đặt mục tiêu đánh giá toàn diện khả năng ứng dụng mô hình khí hậu khu vực RegCM phiên bản 3 (RegCM3) kết hợp với sản phẩm đầu ra từ hệ thống mô hình toàn cầu đại dương - khí quyển CAM-SOM phục vụ bài toán dự báo hạn mùa tại Việt Nam. Phạm vi nghiên cứu bao trùm toàn bộ lãnh thổ Việt Nam và Biển Đông với miền tính tâm đặt tại 11.5°N, 108°E, sử dụng chuỗi số liệu thực nghiệm liên tục 10 năm từ 1996 đến 2005 cùng các thí nghiệm dự báo thử nghiệm giai đoạn 2000-2005. Ý nghĩa khoa học và thực tiễn của công trình thể hiện ở việc thiết lập quy trình tích hợp mô hình động lực hiện đại, giúp kiểm soát sai số nhiệt độ bề mặt trung bình trong ngưỡng 1°C đến 2°C và nâng cao độ chính xác dự báo lượng mưa mùa hè, tạo cơ sở khoa học vững chắc giúp các ngành nông nghiệp và quản lý tài nguyên nước chủ động giảm thiểu từ 15% đến 20% thiệt hại do thiên tai khí hậu gây ra.

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 trên nền tảng lý thuyết hoàn lưu chung khí quyển và lý thuyết tương tác nhiệt động lực học giữa khí quyển - đại dương - mặt đất. Khí hậu hạn mùa trên quy mô khu vực chịu sự chi phối mạnh mẽ bởi các dao động quy mô lớn như hiện tượng ENSO, dao động nhiệt độ mặt nước biển (SST) và biến thiên độ ẩm đất.

Mô hình nghiên cứu cốt lõi bao gồm mô hình khí hậu khu vực RegCM3 kết hợp hệ thống mô hình toàn cầu CAM-SOM (Community Atmosphere Model phiên bản 3.0 tích hợp mô hình đại dương lớp mỏng Slab Ocean Model).

Bốn khái niệm khoa học then chốt được vận dụng xuyên suốt bao gồm:

  • Dự báo hạn mùa (Seasonal forecast): Dự báo trạng thái dị thường của khí quyển với thời trễ (leadtime) từ 0 đến 3 tháng.
  • Hệ tọa độ thẳng đứng thủy tĩnh sigma: Hệ tọa độ bám sát địa hình bề mặt với 18 mực thẳng đứng, đỉnh mô hình đạt 50 mb.
  • Tham số hóa đối lưu mây tích: Cơ chế tính toán giáng thủy quy mô dưới lưới thông qua các sơ đồ Kuo sửa đổi, Grell (khép kín Arakawa - Schubert) và MIT Emanuel.
  • Mô hình đại dương lớp xáo trộn (SOM): Biểu diễn tầng mặt đại dương với độ sâu ấn định từ 70 m đến 100 m nhằm mô phỏng thông lượng nhiệt bề mặt Q và trao đổi năng lượng khí quyển - biển.

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

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

  • Bộ số liệu tái phân tích toàn cầu NCEP/DOE Reanalysis 2 (NNRP2) với độ phân giải lưới 2.5° x 2.5°.
  • Nhiệt độ mặt nước biển phân tích tối ưu OISST độ phân giải 1.0°.
  • Dữ liệu quan trắc nhiệt độ và lượng mưa phân tích CRU độ phân giải 0.5°.
  • Số liệu quan trắc thực tế tại 48 trạm khí tượng điển hình trải đều trên toàn lãnh thổ Việt Nam với tần suất 4 obs quan trắc tiêu chuẩn mỗi ngày trong giai đoạn 1996-2005.

Cỡ mẫu nghiên cứu gồm 48 trạm quan trắc đại diện cho 7 vùng khí hậu đặc trưng của Việt Nam: Tây Bắc (B1), Đông Bắc (B2), Đồng bằng Bắc Bộ (B3), Bắc Trung Bộ (B4), Nam Trung Bộ (N1), Tây Nguyên (N2) và Nam Bộ (N3). Phương pháp chọn mẫu là phân tầng địa lý phi ngẫu nhiên, bảo đảm bao phủ đầy đủ các đặc trưng địa hình phức tạp từ đồng bằng duyên hải đến các đỉnh núi cao trên 1500 m.

Phương pháp phân tích chủ đạo là mô hình hóa số trị động lực học kết hợp thống kê đối sánh sai số (Mean Bias, sai số bình phương trung bình và biểu đồ phân tán). Lý do lựa chọn phương pháp mô hình động lực là khả năng giải quyết trực tiếp hệ phương trình vi phân đầy đủ của khí quyển, tái tạo chính xác các quá trình trao đổi nhiệt ẩm phi tuyến mà phương pháp thống kê truyền thống không thể mô tả.

Tiến trình nghiên cứu được thiết kế qua hai thí nghiệm chính:

  • Thí nghiệm 1 (TN1): Tích hợp RegCM3 trên lưới 144x105 điểm (độ phân giải 36 km) với 3 sơ đồ đối lưu khác nhau, chạy từ ngày 01 tháng 04 đến ngày 31 tháng 10 hàng năm trong chuỗi 10 năm (1996-2005) sử dụng điều kiện biên thực NNRP2.
  • Thí nghiệm 2 (TN2a & TN2b): Kết nối đầu ra CAM-SOM sang RegCM3 thông qua bộ thư viện NetCDF Fortran Interface (NFI) và các module hiệu chỉnh cam42.f, CCSM.f để thực hiện các dự báo mùa thử nghiệm giai đoạn 2000-2005 với leadtime từ 0 đến 3 tháng.

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

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

Quá trình mô phỏng và kiểm chứng thực nghiệm mang lại bốn phát hiện khoa học quan trọng:

Mô hình RegCM3 có khả năng tái tạo xuất sắc cấu trúc hoàn lưu quy mô lớn trên khu vực Đông Nam Á. Trường vectơ gió và độ cao địa thế vị tại các mực chuẩn 850 mb, 500 mb và 200 mb đều có sự tương đồng rất cao so với dữ liệu tái phân tích NNRP2, độ lệch geopotential height chỉ dao động trong khoảng 5 mb. Trường khí áp mực biển (MSLP) mô phỏng có dạng hình thế khớp hoàn toàn với thực tế, độ chênh lệch áp suất trung bình toàn miền duy trì ở mức khoảng 2 mb.

Sơ đồ tham số hóa đối lưu Grell cho hiệu quả mô phỏng trường mưa tối ưu nhất tại Việt Nam. Trong khi sơ đồ Kuo có xu hướng mô phỏng thiên khô rõ rệt (lượng mưa thấp hơn thực tế từ 50 mm đến 100 mm) và sơ đồ Emanuel mô phỏng vượt ngưỡng từ 100 mm đến 200 mm (sai số dương từ 30% đến 45%), sơ đồ Grell tái hiện trường mưa trung bình tháng 5 đến tháng 10 gần sát với dữ liệu lưới CRU, đặc biệt tại khu vực Bắc Trung Bộ (B4).

Trường nhiệt độ không khí mực 2m (T2m) thể hiện tính ổn định cao với sai số hệ thống dao động từ 1°C đến 2°C trên hầu hết các trạm quan trắc. Kết quả mô phỏng tại các vùng đồng bằng (B3, N3) đạt độ chính xác cao nhất, bám sát diễn biến nhiệt độ thực tế. Tại các vùng núi cao có địa hình chia cắt phức tạp như Sa Pa và Đà Lạt, mô hình có xu hướng mô phỏng nhiệt độ cao hơn quan trắc khoảng 2°C do hiệu ứng san bằng độ cao địa hình trong ô lưới 36 km.

Hệ thống ghép nối Reg_CAMSOM chứng minh tính khả thi cao trong bài toán dự báo nghiệp vụ hạn mùa với leadtime từ 0 đến 3 tháng. Kết quả dự báo các tháng 4, 5, 6 giai đoạn 2001-2005 phản ánh rõ nét bước chuyển mùa của hoàn lưu gió mùa mùa hè và giữ vững hệ số tương quan không gian trên 65% đối với trường nhiệt độ tại 48 trạm quan trắc.

Thảo luận kết quả

Nguyên nhân chính dẫn đến sai số nhiệt độ mực 2m thấp hơn quan trắc khoảng 1°C đến 2°C là do sự tương tác giữa bức xạ mặt trời và sơ đồ bề mặt đất BATS1E trong việc tính toán cân bằng nhiệt ẩm lớp đệm. Sự phân hóa sâu sắc giữa 3 sơ đồ đối lưu xuất phát từ bản chất giả thiết khép kín: sơ đồ Kuo dựa trên mức độ hội tụ ẩm quy mô lớn, sơ đồ Emanuel phụ thuộc vào thế năng đối lưu khả dụng (CAPE), trong khi sơ đồ Grell cân bằng giữa dòng thăng và dòng giáng mây tích, giúp thích ứng tốt nhất với khí hậu nhiệt đới gió mùa nhiều đồi núi tại Việt Nam.

So sánh với các nghiên cứu quốc tế như dự án DEMETER của Liên minh Châu Âu hay hệ thống dự báo khí hậu CFS của NCEP Hoa Kỳ, kết quả của luận văn khẳng định mô hình khí hậu khu vực có kỹ năng vượt trội trong việc chi tiết hóa các trường khí tượng địa phương so với các mô hình hoàn lưu chung toàn cầu GCM thuần túy.

Các tập dữ liệu nghiên cứu được biểu diễn trực quan qua hệ thống biểu đồ và bảng phân tích khoa học:

  • Đồ thị tụ điểm phân tán (scatter plot) 2 chiều so sánh giá trị quan trắc tại 48 trạm với giá trị mô phỏng của 3 sơ đồ đối lưu cho thấy mức độ tập trung sai số nhiệt độ quanh đường 1:1.
  • Bản đồ sai số không gian phân bố 2 chiều thể hiện mức chênh lệch nhiệt độ và giáng thủy trên từng vùng địa hình.
  • Biểu đồ mặt cắt thẳng đứng (vertical profile) của nhiệt độ, độ ẩm và biểu đồ Hovmoller trung bình trượt 5 ngày minh họa chi tiết sự dịch chuyển của các hệ thống thời tiết theo vĩ độ.

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

Dựa trên các kết quả nghiên cứu thực nghiệm, bốn nhóm giải pháp mang tính chiến lược được đề xuất nhằm nâng cao năng lực dự báo hạn mùa tại Việt Nam:

  • Tích hợp và vận hành chính thức mô hình RegCM3 trong quy trình nghiệp vụ dự báo hạn mùa. Ứng dụng sản phẩm dự báo toàn cầu từ hệ thống CAM-SOM và CFS làm điều kiện biên động lực. Mục tiêu giảm sai số dự báo nhiệt độ trung bình tháng xuống dưới 1.2°C và nâng tỷ lệ dự báo đúng xu thế lượng mưa đạt trên 75%. Lộ trình thực hiện trong vòng 12 tháng do Trung tâm Dự báo Khí tượng Thủy văn Quốc gia chủ trì.
  • Nâng cấp độ phân giải lưới tính toán từ 36 km xuống 10 km đến 12 km. Ứng dụng các phép chiếu bản đồ hiện đại và tinh chỉnh các module tiền xử lý dữ liệu Fortran nhằm mô phỏng chính xác các tâm mưa lớn cực đoan có địa hình đón gió phức tạp như Bắc Quang hay Trà My. Mục tiêu giảm sai số cục bộ giáng thủy vùng núi từ 25% đến 30%. Thời gian triển khai trong 18 đến 24 tháng bởi các nhóm nghiên cứu mô hình số trị chuyên sâu.
  • Phát triển hệ thống dự báo tổ hợp đa mô hình (Multi-Model Ensemble). Kết hợp kỹ thuật hiệu chỉnh sai số sau mô hình (Model Output Statistics - MOS) giữa các sơ đồ đối lưu Grell và Emanuel để loại bỏ sai số hệ thống. Đích đến là kéo dài thời trễ dự báo hiệu quả lên 3 đến 6 tháng với độ tin cậy tăng thêm 20%. Kế hoạch thực hiện trong giai đoạn 3 năm phối hợp cùng Tổ chức Khí tượng Thế giới (WMO).
  • Hiện đại hóa cơ sở dữ liệu vi khí hậu và số hóa tham số mặt đệm đất. Cập nhật bản đồ thảm phủ thực vật GLCC và độ cao địa hình USGS độ phân giải 30 giây cung góc cho toàn bộ 63 tỉnh thành. Mục tiêu bảo đảm 100% dữ liệu đầu vào thời gian thực, rút ngắn thời gian xử lý biên xuống dưới 30 phút mỗi chu kỳ dự báo. Trách nhiệm thực hiện thuộc Cục Khí tượng Thủy văn và các viện nghiên cứu chuyên ngành trong vòng 2 năm.

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

Luận văn là tài liệu tham khảo khoa học giá trị cao dành cho 4 nhóm đối tượng trọng tâm:

  • Nhà nghiên cứu và nghiên cứu sinh chuyên ngành Khí tượng, Khí hậu học: Tiếp cận toàn diện phương pháp luận ghép nối mô hình toàn cầu với mô hình khu vực, kỹ thuật can thiệp mã nguồn Fortran và xử lý cấu trúc dữ liệu NetCDF chuyên sâu phục vụ các đề tài nghiên cứu biến đổi khí hậu.
  • Dự báo viên tại các cơ quan khí tượng thủy văn: Nắm giữ bộ số liệu thực nghiệm kiểm chứng 10 năm tại 48 trạm quan trắc và quy tắc lựa chọn sơ đồ tham số hóa đối lưu thích hợp cho từng vùng khí hậu nhằm xây dựng các bản tin dự báo mùa nghiệp vụ chuẩn xác.
  • Cơ quan quản lý nông nghiệp và quy hoạch tài nguyên nước: Khai thác thông tin dự báo biến động nhiệt ẩm và xu thế mùa mưa trước từ 1 đến 3 tháng để lập kế hoạch điều tiết 100% dung tích hồ chứa thủy điện, bố trí cơ cấu mùa vụ và chủ động phòng ngừa hạn mặn.
  • Giảng viên và sinh viên các trường đại học khối Khoa học Trái đất: Sử dụng làm giáo trình tham khảo mẫu mực về mô hình hóa số trị, thiết kế thí nghiệm khí hậu và phương pháp xử lý số liệu quan trắc khí quyển.

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

Tại sao phương pháp mô hình động lực lại ưu việt hơn phương pháp thống kê trong dự báo hạn mùa?
Phương pháp mô hình động lực dựa trên việc giải các hệ phương trình vật lý đầy đủ mô tả tương tác phức tạp giữa đại dương và khí quyển. Trong khi phương pháp thống kê thất bại trước các hiện tượng chuyển pha đột ngột, mô hình động lực như RegCM3 có thể nắm bắt chính xác các đợt ENSO mạnh và tái tạo trường địa thế vị với sai số chỉ khoảng 5 mb.

Sơ đồ tham số hóa đối lưu nào phù hợp nhất cho điều kiện khí hậu Việt Nam trong mô hình RegCM3?
Sơ đồ Grell với giả thiết khép kín Arakawa - Schubert cho kết quả mô phỏng tốt nhất về cả hoàn lưu và phân bố mưa. Trong thực tế thí nghiệm 10 năm, sơ đồ Grell phản ánh lượng mưa trung bình tháng 5 đến tháng 10 khớp nhất với số liệu CRU, tránh được hiện tượng quá khô của sơ đồ Kuo (hụt 50-100 mm) và quá ẩm của sơ đồ Emanuel (thừa 100-200 mm).

Mô hình RegCM3 thể hiện độ chính xác nhiệt độ bề mặt như thế nào trên các vùng khí hậu khác nhau?
Mô hình kiểm soát sai số nhiệt độ mực 2m trung bình từ 1°C đến 2°C trên toàn bộ 48 trạm quan trắc. Độ chính xác đạt mức tối ưu tại các vùng đồng bằng Bắc Bộ (B3) và Nam Bộ (N3). Tuy nhiên, tại các trạm núi cao như Sa Pa hay Đà Lạt, mô hình cho giá trị cao hơn thực tế khoảng 2°C do độ phân giải 36 km làm mịn địa hình.

Việc tích hợp mô hình đại dương CAM-SOM giải quyết vấn đề gì trong bài toán dự báo mùa?
Mô hình CAM-SOM cung cấp trường nhiệt độ mặt nước biển SST tương tác liên tục thông qua tầng xáo trộn sâu 70 m đến 100 m. Điều này cho phép RegCM3 cập nhật động lực các thông lượng nhiệt ẩm bề mặt biển, tạo điều kiện mở rộng hạn dự báo từ 1 tháng lên 3 tháng với độ ổn định cao.

Cách thức xử lý kỹ thuật để kết nối dữ liệu từ mô hình toàn cầu CAM-SOM sang mô hình khu vực RegCM3?
Tác giả đã lập trình chỉnh sửa các module Fortran cam42.fCCSM.f, sử dụng thư viện NetCDF Fortran Interface (NFI) để giải mã trường phổ T42 (~2.81°). Hệ thống thực hiện phép nội suy thẳng đứng từ tọa độ lai sang các mực áp suất chuẩn và nội suy song tuyến tính sang lưới Mercator 36 km của RegCM3.

Kết luận

  • Luận văn đã xây dựng thành công quy trình công nghệ ứng dụng mô hình khí hậu khu vực RegCM3 kết nối với hệ thống mô hình toàn cầu CAM-SOM phục vụ dự báo hạn mùa tại Việt Nam.
  • Kết quả kiểm chứng 10 năm (1996-2005) khẳng định mô hình tái tạo chuẩn xác cấu trúc hoàn lưu quy mô lớn với sai số địa thế vị dưới 5 mb và sai số khí áp mực biển chỉ khoảng 2 mb.
  • Sơ đồ đối lưu Grell được chứng minh là lựa chọn tối ưu nhất cho khu vực gió mùa Đông Nam Á, kiểm soát sai số nhiệt độ mực 2m tại 48 trạm quan trắc trong ngưỡng ổn định 1°C đến 2°C.
  • Thử nghiệm dự báo hạn mùa giai đoạn 2001-2005 với leadtime từ 0 đến 3 tháng cho thấy khả năng nắm bắt rõ nét thời điểm bắt đầu và cường độ của gió mùa mùa hè.
  • Đóng góp then chốt của công trình là giải quyết triệt để bài toán tiền xử lý dữ liệu ghép nối mô hình, mở đường cho hướng tiếp cận dự báo số trị hạn dài hiện đại tại Việt Nam.

Trong vòng 24 tháng tới, việc tiếp tục nâng cấp độ phân giải lưới tính toán lên 10 km và tích hợp dự báo tổ hợp đa mô hình là bước đi tất yếu. Các đơn vị nghiên cứu và cơ quan khí tượng thủy văn cần đẩy mạnh ứng dụng ngay kết quả luận văn vào nghiệp vụ thực tế nhằm nâng cao chất lượng cảnh báo sớm thiên tai, phục vụ phát triển kinh tế - xã hội bền vững.