Tổng quan nghiên cứu

Mưa đá là một trong những loại hình thiên tai cực đoan có sức tàn phá nghiêm trọng nhưng vô cùng khó dự báo tại Việt Nam. Theo thống kê của Ban Chỉ đạo Quốc gia về Phòng, chống Thiên tai, chỉ riêng trong 7 tháng đầu năm 2020, tổ hợp thiên tai dông lốc, sét và mưa đá đã khiến 44 người thiệt mạng, 127 người bị thương, làm hư hại hơn 50.000 ngôi nhà và phá hủy gần 40.000 hecta hoa màu, gây thiệt hại kinh tế ước tính vượt mức 500 tỷ đồng. Trên phạm vi toàn cầu, mưa đá gây tổn thất trung bình hơn 10 tỷ USD mỗi năm tại Hoa Kỳ và hơn 4,5 tỷ USD vào năm 2015 tại Trung Quốc.

Vấn đề cốt lõi của bài toán dự báo mưa đá nằm ở quy mô không gian nhỏ (dưới 10 km), vòng đời tồn tại ngắn và diễn biến đối lưu diễn ra cực kỳ nhanh chóng. Khu vực miền Bắc Việt Nam, đặc biệt là vùng núi Tây Bắc, sở hữu địa hình chia cắt phức tạp kết hợp hoàn lưu gió mùa tạo nên một trung tâm mưa đá có tần suất xuất hiện cao trong thời kỳ chuyển mùa.

Mục tiêu cụ thể của luận văn là ứng dụng mô hình số trị độ phân giải cao WRF kết hợp mô-đun Hailcast và chỉ số độ xoáy dòng thăng (UH) nhằm thử nghiệm, đánh giá khả năng dự báo vị trí và kích thước hạt mưa đá. Nghiên cứu tập trung vào hai đợt thiên tai điển hình xảy ra vào ngày 17/03/2020 và ngày 24/03/2020 tại các tỉnh Lai Châu, Lào Cai, Sơn La, Yên Bái và Điện Biên. Ý nghĩa của công trình thể hiện ở việc ứng dụng chỉ số kỹ năng không gian FSS cùng thuật toán radar phân cực đôi HSDA để lượng hóa sai số dự báo, mở ra giải pháp cảnh báo sớm hiệu quả cho ngành Khí tượng Thủy văn Việt Nam.

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 nhiệt động lực học khí quyển và vi vật lý mây đối lưu sâu. Quá trình hình thành mưa đá trong mây vũ tích đòi hỏi 4 điều kiện môi trường then chốt: dòng thăng mạnh trên nền khí quyển bất ổn định, lượng ẩm dồi dào ở tầng thấp, độ đứt gió thẳng đứng lớn và cơ chế kích phát động lực từ địa hình. Theo phân loại khí hậu mưa đá toàn cầu, mưa đá tại miền Bắc Việt Nam chủ yếu thuộc loại III, đặc trưng bởi vĩ độ trung bình - cận nhiệt đới trên địa hình phức tạp gắn liền với độ đứt gió trung bình.

Khung mô hình nghiên cứu tích hợp các khái niệm chính bao gồm:

  • Phổ phân bố kích thước hạt ngưng kết (PSD) dạng hàm gamma mô tả sự chuyển pha giữa nước mây, nước mưa, băng ngưng kết và mưa đá.
  • Chỉ số độ xoáy dòng thăng (Updraft Helicity - UH), đại diện cho tích phân độ xoáy thẳng đứng nhân với vận tốc dòng thăng trong lớp 2 - 5 km, biểu thị cường độ xoay của ổ dông đối lưu.
  • Mô hình tăng trưởng kích thước hạt mưa đá một chiều HAILCAST nhúng trong mô hình WRF, tính toán sự gia tăng kích thước hạt dựa trên cân bằng nhiệt và khối lượng của 5 kích thước nhân ngưng kết ban đầu.
  • Năng lượng đối lưu tiềm tàng khả dụng (CAPE) và hội tụ thông lượng ẩm tích lũy thẳng đứng (VIMFC) làm thước đo điều kiện quy mô synop thuận lợi.

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

Nguồn dữ liệu đầu vào phục vụ mô phỏng bao gồm số liệu dự báo toàn cầu GFS độ phân giải 0,25 độ với tần suất 3 giờ/lần của NCEP và số liệu tái phân tích ERA5 của ECMWF. Dữ liệu thực nghiệm kiểm chứng sử dụng số liệu radar thời tiết Doppler phân cực đôi băng C đặt tại Pha Đin (Điện Biên) ở độ cao 1470 mét, quét chu kỳ 10 phút/lần trong bán kính hữu dụng 100 - 120 km. Dữ liệu này được xử lý bằng thư viện Py-ART để tính toán thuật toán phân loại kích thước hạt mưa đá (HSDA) gồm 3 cấp: mưa đá nhỏ (5 - 25 mm), mưa đá lớn (25 - 50 mm) và mưa đá rất lớn (trên 50 mm). Ngoài ra, báo cáo thiệt hại và quan trắc mặt đất tại hơn 20 điểm dân cư được thu thập để đối chiếu.

Cấu hình mô hình WRF phiên bản 4.3 được thiết lập với 3 miền lưới lồng có độ phân giải ngang lần lượt là 13,5 km, 4,5 km và 1,5 km, trải trên 51 mực thẳng đứng. Miền tính 1,5 km bao phủ toàn bộ khu vực Bắc Bộ nhằm giải quyết trực tiếp bài toán đối lưu hiển mà không cần sơ đồ tham số hóa đối lưu. Bốn thí nghiệm vi vật lý bậc cao được lựa chọn để so sánh bao gồm: Thompson (MP008), Milbrandt-Yau (MP009), Morrison-graupel (MP101 với mật độ graupel 500 kg/m³) và Morrison-hail (MP102 với mật độ 900 kg/m³).

Phương pháp phân tích sử dụng điểm kỹ năng không gian Fractions Skill Score (FSS) với ngưỡng hữu ích L = 0,5 để đánh giá trường nhị phân theo các bán kính lân cận n. Việc lựa chọn FSS thay vì các chỉ số đánh giá điểm truyền thống giúp loại bỏ lỗi phạt kép (double penalty), cho phép đánh giá chính xác kỹ năng dự báo hình thái không gian của các hiện tượng cực đoan quy mô nhỏ dưới 10 km.

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

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

Kết quả phân tích định lượng từ các thử nghiệm mô phỏng chỉ ra những phát hiện quan trọng:

Thứ nhất, trong đợt mưa đá ngày 17/03/2020, mô hình WRF thể hiện kỹ năng dự báo không gian vượt trội. Chỉ số UH đạt quy mô kỹ năng tối thiểu (scalemin) ở khoảng cách 25 - 30 điểm lưới, tương đương bán kính không gian 30 - 40 km. Biến kích thước hạt cực đại HC từ mô hình Hailcast đạt scalemin ở quy mô 30 - 40 điểm lưới (khoảng 40 - 50 km). Dải mưa đá chính kéo dài theo trục đông - tây qua Lai Châu và Lào Cai được mô phỏng sát thực tế quan trắc radar HSDA.

Thứ hai, độ nhạy của các sơ đồ vi vật lý tạo ra sự phân hóa rõ rệt về diện tích và cường độ dự báo. Sơ đồ Morrison-graupel (MP101) và Milbrandt-Yau (MP009) mô phỏng dòng thăng phát triển rất mạnh, xuyên phá lên tới mực áp suất 200 mb với bề rộng dòng thăng đạt khoảng 8 km và vùng mây băng mở rộng trên 25 km. Ngược lại, sơ đồ Thompson (MP008) mô phỏng dòng thăng yếu hơn, chỉ đạt tới mực 300 mb với bề rộng 5 km, đồng thời tạo ra vùng dự báo khống rất lớn tại Sơn La chiếm hơn 66% tổng diện tích cảnh báo, khiến điểm FSS của biến HC không vượt qua được ngưỡng kỹ năng 0,5.

Thứ ba, đợt mưa đá ngày 24/03/2020 có tính chất đa ổ đối lưu phân tán, khiến điểm kỹ năng FSS giảm đáng kể. Các sơ đồ MP009 và MP102 cần quy mô không gian lân cận lên tới 70 - 80 điểm lưới (khoảng 100 - 120 km) mới đạt ngưỡng kỹ năng FSS trên 0,5. Điều này phản ánh thách thức lớn khi mô hình đối lưu hiển mô phỏng các ổ dông độc lập có kích thước dưới 10 km.

Thứ tư, môi trường nhiệt động quy mô lớn kích phát mưa đá miền Bắc được định lượng rõ nét: năng lượng MCAPE ở mức trung bình từ 1000 - 1500 J/kg, độ đứt gió thẳng đứng lớp 0 - 6 km đạt giá trị cao từ 20 - 25 m/s, kết hợp với hội tụ ẩm mực thấp VIMFC vượt ngưỡng 200 × 10⁻⁶ kg m⁻² s⁻¹.

Thảo luận kết quả

Khác với các đợt mưa đá tại vùng Đại bình nguyên Hoa Kỳ vốn đòi hỏi chỉ số CAPE bùng nổ trên 3000 J/kg, mưa đá tại miền Bắc Việt Nam hình thành chủ yếu nhờ sự tương tác giữa độ đứt gió thẳng đứng mạnh (20 - 25 m/s) và quá trình nâng cưỡng bức địa hình của dãy Hoàng Liên Sơn. Điều này giải thích tại sao giản đồ thiên khí Skew-T tại tâm ổ dông ghi nhận mức năng lượng MCAPE giải phóng chỉ quanh mức 500 J/kg nhưng dòng thăng vẫn đạt vận tốc thẳng đứng cực đại. Cơ chế này hoàn toàn tương đồng với các nghiên cứu mưa đá tại khu vực Vân Nam (Trung Quốc) của các tác giả quốc tế.

Khi biểu diễn dữ liệu qua các mặt cắt thẳng đứng - vĩ hướng, cấu trúc đối lưu dạng dông siêu ổ (supercell) hình chữ T được tái hiện rõ nét ở các sơ đồ MP009 và MP101. Vùng lõi có tỉ lệ xáo trộn pha graupel và nước đá dày đặc đã đẩy đường đẳng nhiệt 0°C lên cao hơn môi trường xung quanh, phản ánh quá trình giải phóng nhiệt ẩn ngưng kết mãnh liệt. Việc chỉ số UH có xu hướng dự báo cường độ dông quá mức trong khi Hailcast (HC) đưa ra kích thước hạt thực tế hơn (phổ biến 20 - 35 mm) cho thấy vai trò thiết yếu của việc tích hợp các phương trình động lực hạt đá thay vì chỉ dựa vào các chỉ số thay thế thuần túy.

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

  1. Tích hợp mô-đun Hailcast vào hệ thống dự báo nghiệp vụ số trị WRF quốc gia: Trung tâm Dự báo Khí tượng Thủy văn Quốc gia cần chuẩn hóa cấu hình mô hình WRF phân giải 1,5 km sử dụng sơ đồ vi vật lý Morrison-graupel (MP101), phấn đấu nâng tỷ lệ dự báo sớm dông sét và mưa đá trước 18 - 24 giờ đạt độ chính xác không gian trên 70% trong giai đoạn 2024 - 2026.
  2. Nâng cấp và hiệu chuẩn sản phẩm HSDA từ mạng lưới radar phân cực đôi: Cục Khí tượng Thủy văn cần triển khai thuật toán HSDA tự động trên toàn bộ 05 trạm radar Doppler băng C phân cực đôi (Pha Đin, Việt Trì, Pleiku, Quy Nhơn, Nha Trang), hoàn thiện quy trình lọc nhiễu địa hình vùng núi nhằm giảm tỷ lệ báo động giả xuống dưới 15% trước năm 2026.
  3. Thiết lập mạng lưới quan trắc bề mặt chuyên dụng cho mưa đá: Ban Chỉ đạo Quốc gia về Phòng, chống Thiên tai phối hợp với các tỉnh miền núi phía Bắc lắp đặt hệ thống tấm đo mưa đá tiêu chuẩn (hailpads) và phát triển ứng dụng di động cho phép cộng đồng gửi dữ liệu hình ảnh tọa độ thời gian thực, mục tiêu thu thập tối thiểu 50 điểm kiểm chứng độc lập cho mỗi sự kiện thiên tai từ năm 2025.
  4. Triển khai các biện pháp công trình bảo vệ sản xuất nông nghiệp: Sở Nông nghiệp và Phát triển Nông thôn tại các tỉnh Lai Châu, Lào Cai, Sơn La cần hỗ trợ nông dân lắp đặt hệ thống lưới che chắn mưa đá cho trên 65% diện tích vườn cây ăn quả giá trị cao và vùng trồng hoa màu chuyên canh trước mùa dông bão tháng 3 - tháng 5 hằng năm.

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

  • Dự báo viên khí tượng tại các đài khu vực và trung ương: Nắm vững phương pháp khai thác chỉ số độ xoáy dòng thăng UH và sản phẩm kích thước hạt HC từ mô hình WRF để đưa ra các bản tin cảnh báo thiên tai cực đoan sớm trước 12 - 36 giờ.
  • Nhà nghiên cứu và học viên cao học chuyên ngành Khí tượng, Khí hậu học: Sử dụng làm tài liệu tham khảo chuyên sâu về thiết lập mô hình đối lưu hiển, lựa chọn sơ đồ tham số hóa vi vật lý bậc cao và phương pháp chấm điểm không gian FSS cho hiện tượng hiếm.
  • Cán bộ quản lý rủi ro thiên tai và tìm kiếm cứu nạn cấp tỉnh/huyện: Hiểu rõ các quy luật phân bố không gian, ngưỡng kích thước hạt và thời gian xuất hiện mưa đá tại khu vực Tây Bắc để xây dựng phương án phòng ngừa, sơ tán và cứu trợ kịp thời.
  • Chủ trang trại nông nghiệp công nghệ cao và doanh nghiệp bảo hiểm nông nghiệp: Ứng dụng các phân tích về tần suất, quy mô và mức độ thiệt hại của mưa đá để xác định vùng an toàn đầu tư, thiết kế kết cấu nhà màng và xây dựng biểu phí bảo hiểm cây trồng phù hợp.

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

Mô hình WRF độ phân giải 1,5 km có ưu điểm gì vượt trội so với mô hình dự báo truyền thống?
Độ phân giải siêu cao 1,5 km cho phép mô hình mô phỏng trực tiếp các quá trình đối lưu quy mô mây mà không cần dùng sơ đồ tham số hóa thô sơ. Điều này giúp nắm bắt chính xác dòng thăng mạnh trong mây vũ tích, nâng cao khả năng cảnh báo vị trí dông lốc trước 24 giờ.

Sự khác biệt cốt lõi giữa chỉ số UH và biến HC trong dự báo mưa đá là gì?
Chỉ số UH đánh giá gián tiếp cường độ xoáy của dòng thăng trong ổ dông (ngưỡng trên 50 m²s⁻² cho mưa đá), trong khi biến HC từ mô hình Hailcast mô phỏng trực tiếp sự lớn lên và tan chảy của hạt đá, cung cấp đường kính hạt mưa đá cụ thể (mm) tại mặt đất.

Tại sao sơ đồ vi vật lý Morrison (MP101) lại cho kết quả mô phỏng tốt hơn sơ đồ Thompson (MP008)?
Sơ đồ Morrison dự báo bậc hai cho cả tỉ lệ xáo trộn và tổng số hạt ngưng kết, giúp tái hiện chân thực cấu trúc dòng thăng lên tới mực 200 mb và vùng mây băng rộng trên 25 km, trong khi Thompson tạo diện dự báo khống lớn tại Sơn La vượt quá 66%.

Thuật toán HSDA trên radar Pha Đin đóng vai trò gì trong nghiên cứu?
Thuật toán HSDA tận dụng các biến phân cực đôi (Z, Zdr, CC) để phân loại kích thước hạt mưa đá ở góc quét thấp 0,5 độ. Dữ liệu này tạo ra bức tranh quan trắc không gian khách quan, làm cơ sở đối chiếu và chấm điểm kỹ năng FSS cho mô hình số.

Những dấu hiệu thiên khí nào cho thấy nguy cơ cao xuất hiện mưa đá ở miền Bắc?
Sự kết hợp giữa không khí lạnh nén rãnh áp thấp ở tầng thấp với rãnh gió tây sâu trên mực 500 hPa, đi kèm hội tụ ẩm VIMFC trên 200 × 10⁻⁶ kg m⁻² s⁻¹ và độ đứt gió 0 - 6 km đạt 20 - 25 m/s là điều kiện kinh điển gây mưa đá diện rộng.

Kết luận

  • Luận văn đã thiết lập và thử nghiệm thành công mô hình WRF 1,5 km tích hợp mô-đun Hailcast và chỉ số UH để dự báo mưa đá tại miền Bắc Việt Nam với độ chính xác cao ở hạn dự báo 18 - 30 giờ.
  • Đánh giá không gian bằng chỉ số FSS khẳng định kỹ năng dự báo vị trí mưa đá đạt độ tin cậy trong phạm vi bán kính 30 - 40 km đối với các đợt dông đối lưu quy mô lớn.
  • Khẳng định sơ đồ vi vật lý Morrison-graupel (MP101) và Milbrandt-Yau (MP009) có khả năng mô phỏng cấu trúc vi vật lý mây dông siêu ổ và phân bố kích thước hạt vượt trội so với các sơ đồ truyền thống.
  • Làm sáng tỏ vai trò quyết định của độ đứt gió tầng trung bình (20 - 25 m/s) và dòng hội tụ ẩm mực thấp kết hợp địa hình Hoàng Liên Sơn trong cơ chế phát sinh mưa đá nhiệt đới.
  • Định hướng nghiên cứu giai đoạn 2025 - 2030 tập trung vào việc đồng hóa số liệu radar phân cực đôi đa trạm và ứng dụng thuật toán học máy nhằm hoàn thiện hệ thống cảnh báo sớm thiên tai thời gian thực phục vụ cộng đồng.