Tổng quan nghiên cứu

Lượng phát thải NOx tại khu vực châu Á sau năm 1990 đã vượt ngưỡng 28 triệu tấn/năm, gia tăng nhanh chóng do quá trình công nghiệp hóa và đốt sinh khối. Theo ước tính từ các báo cáo khí quyển quốc tế, xon khí làm suy giảm bức xạ bề mặt trung bình 3,4 W/m2 trên phạm vi toàn cầu và dao động từ 10 đến 20 W/m2 tại khu vực gió mùa châu Á, trực tiếp làm thay đổi gradient khí áp và hoàn lưu khí quyển. Tại Việt Nam, các mô hình dự báo thời tiết truyền thống thường bỏ qua tương tác hai chiều giữa các quá trình động lực và hóa học khí quyển. Tuy nhiên, các nghiên cứu thực nghiệm đã chứng minh xon khí cacbon đen có thể làm lượng mưa tăng từ 1 đến 4 mm/ngày ở miền Bắc và giảm từ 1 đến 4 mm/ngày ở miền Nam.

Vấn đề cốt lõi đặt ra là việc thiếu hụt một công cụ mô hình hóa tích hợp có khả năng mô phỏng đồng thời trường khí tượng và quá trình phát tán, chuyển hóa của các chất ô nhiễm. Luận văn hướng đến mục tiêu thử nghiệm và đánh giá khả năng ứng dụng mô hình số trị WRF/Chem phiên bản 3.4 tại khu vực Việt Nam, tập trung kiểm chứng khả năng tái tạo trường nhiệt độ, lượng mưa, hoàn lưu gió và nồng độ bụi 5 cấp kích thước. Phạm vi nghiên cứu bao phủ không gian từ 1,5°S đến 32°N và 88°E đến 124°E, với thời gian tích phân thử nghiệm liên tục trong giai đoạn từ 00h ngày 01/01/2006 đến 00h ngày 05/01/2006. Kết quả nghiên cứu có ý nghĩa thực tiễn quan trọng, đặt nền móng cho việc xây dựng hệ thống dự báo chất lượng không khí đồng thời, kiểm soát sai số nhiệt độ dưới 1°C và nâng cao độ chính xác cảnh báo ô nhiễm môi trường.

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 tương tác xon khí - bức xạ - mây và hệ phương trình động lực học khí quyển bất thủy tĩnh. Khung mô hình trung tâm là WRF/Chem phiên bản 3.4, cho phép tích hợp đồng thời các quá trình động lực khí tượng với các phản ứng hóa học khí quyển trong cùng một bước tích phân.

Ba khái niệm học thuật then chốt được vận hành bao gồm:

  • Mô hình ghép nối đồng thời (Online Coupling): Khác với mô hình không đồng thời (Offline), phương pháp này tính toán tương tác hóa học và khí tượng trên cùng một lưới tính tại từng bước thời gian, phản ánh chính xác tác động hồi tiếp của xon khí lên hệ thống khí hậu.
  • Độ dày quang học của khí quyển (AOD): Đại lượng đo mức độ suy giảm bức xạ mặt trời khi truyền qua lớp xon khí trong khí quyển.
  • Hồi tiếp bức xạ xon khí: Cơ chế xon khí trực tiếp hấp thụ, tán xạ ánh sáng và gián tiếp làm thay đổi vi vật lý mây, làm thay đổi phân bố nhiệt độ và lượng mưa.

Về cơ chế hóa học, luận văn áp dụng mô hình GOCART mô phỏng 5 dải bán kính hạt bụi từ 0,1 đến 10 µm, cơ chế hóa học pha khí RADM2 xử lý 24 phản ứng quang hóa và mô hình xon khí MADE/SORGAM mô phỏng các sol khí hữu cơ thứ cấp.

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

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

  • Dữ liệu khí tượng: Phân tích toàn cầu NCEP FNL độ phân giải 1,0 độ cập nhật 6 giờ/lần với 26 mực áp suất từ 1000 hPa đến 10 hPa làm điều kiện biên và ban đầu.
  • Dữ liệu kiểm định: Số liệu tái phân tích NCEP NNRP độ phân giải 2,5 độ và chuỗi số liệu quan trắc lưới APHRODITE của Nhật Bản độ phân giải 0,25 độ.
  • Dữ liệu phát thải: Cơ sở dữ liệu tái phân tích đối lưu RETRO độ phân giải 0,5 độ và dữ liệu phát thải toàn cầu EDGAR độ phân giải 1x1 độ, được xử lý qua công cụ Prep_chem_sources.

Về quy mô mẫu và lưới tính, miền mô phỏng được thiết kế gồm 27 mực thẳng đứng và độ phân giải ngang 30 km trên lưới so le Arakawa C. Miền tính bao trùm tọa độ 1,5°S đến 32°N và 88°E đến 124°E, trong khi miền phân tích tập trung từ 0°N đến 30°N và 90°E đến 120°E nhằm loại trừ nhiễu động biên. Phương pháp chọn miền tính mở rộng sang phía đông bắc và tây nam nhằm nắm bắt đầy đủ các hệ thống hoàn lưu gió mùa chi phối vận chuyển ô nhiễm. Lý do lựa chọn mô hình đồng thời là nhằm loại bỏ hoàn toàn sai số nội suy thời gian và không gian, vốn gây ra độ lệch từ 30% đến 40% trong các mô hình ngoại tuyến truyền thống. Thời gian thực hiện mô phỏng kéo dài liên tục 96 giờ.

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

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

Thứ nhất, mô hình tái hiện xuất sắc trường nhiệt độ không khí mực 2m tại khu vực Việt Nam với sai số chỉ dao động trong khoảng ±1°C so với số liệu quan trắc APHRODITE. Trường gió tại các mực áp suất 850 hPa, 500 hPa và 200 hPa được mô phỏng chuẩn xác cả về hướng và vận tốc với độ lệch vận tốc dưới 2 m/s so với dữ liệu tái phân tích NCEP NNRP.

Thứ hai, nồng độ bụi loại 1 (DUST_1) thể hiện sự phân tầng rõ rệt theo chiều thẳng đứng. Ở tầng cao 500 hPa và 200 hPa, nồng độ bụi đạt mức rất cao từ 10 đến 12 x 10^12 µg/kg, cao gấp 6 đến 10 lần so với tầng sát bề mặt 1000 hPa. Khu vực Nam Bộ và Nam Trung Bộ chịu tác động xâm nhập bụi mạnh nhất từ dòng vận chuyển phía đông, trong khi miền Bắc duy trì nồng độ thấp chỉ từ 1 đến 2 x 10^12 µg/kg do nằm trong tâm áp cao.

Thứ ba, trường lượng mưa mô phỏng có độ sai lệch trung bình ±5 mm trên phần lớn lãnh thổ Việt Nam. Tuy nhiên, tại khu vực Trung Trung Bộ và Tây Nam Bộ trong 48 giờ đầu tiên, mô hình ghi nhận sai số thiên dương cục bộ lên tới ±10 mm.

Thứ tư, khi kích hoạt tương tác bức xạ và hóa học đầy đủ (tùy chọn chem_opt = 11 và chem_opt = 300), nhiệt độ mô phỏng tại miền Bắc Việt Nam giảm khoảng 3°C so với trường hợp không tích hợp hóa học (WRF_NOCHEM).

Thảo luận kết quả

Mức giảm nhiệt độ 3°C tại miền Bắc khi tích hợp mô-đun hóa học phản ánh chính xác hiệu ứng làm mát bề mặt của xon khí thông qua việc hấp thụ và tán xạ bức xạ mặt trời xuống đất. Phát hiện này hoàn toàn phù hợp với kết luận của Menon và cộng sự khi chỉ ra xon khí làm giảm bức xạ mặt đất từ 10 đến 20 W/m2 tại châu Á, cũng như nghiên cứu của Paolo Tuccella tại châu Âu với mức sai số nhiệt độ tương đương 3°C.

Về động lực học, dòng xiết gió tây tại mực 200 hPa với vận tốc lên tới 50 m/s đóng vai trò là hành lang vận chuyển bụi chính từ phía tây sang đông, kết hợp với dòng gió bắc đưa bụi lan tỏa xuống vĩ độ thấp từ 0°N đến 18°N. Toàn bộ các phát hiện này có thể được biểu diễn trực quan và khoa học thông qua bảng ma trận sai số thống kê (chỉ số BIAS và RMSE giữa các tùy chọn hóa học 300, 301, 11 so với APHRODITE), kết hợp biểu đồ mặt cắt kinh hướng (trung bình từ 102°E đến 110°E) và mặt cắt vĩ hướng (trung bình từ 8°N đến 24°N) thể hiện thông lượng vận chuyển bụi theo trục thẳng đứng.

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

  • Chuẩn hóa cơ sở dữ liệu kiểm kê phát thải địa phương có độ phân giải cao 0,1x0,1 độ thay thế nguồn dữ liệu toàn cầu EDGAR 1x1 độ, nhằm giảm sai số phát thải đầu vào xuống dưới 15% trong vòng 12 tháng, do Cục Kiểm soát Ô nhiễm Môi trường chủ trì thực hiện.
  • Tích hợp mô-đun đồng hóa số liệu quang học xon khí (AOD) từ vệ tinh MODIS và mạng lưới quan trắc mặt đất AERONET vào hệ thống WRFDA-Chem, nâng cao hệ số tương quan dự báo bụi mịn PM2.5 lên trên 0,85 theo lộ trình 18 tháng, do Trung tâm Dự báo Khí tượng Thủy văn Quốc gia phụ trách.
  • Hiệu chỉnh các tham số vi vật lý mây và sơ đồ rửa trôi ẩm trong mô hình WRF/Chem, nhằm giảm sai số thiên dương lượng mưa tại dải ven biển miền Trung từ ±10 mm xuống dưới ±3 mm trong thời hạn 2 năm, thực hiện bởi nhóm nghiên cứu mô hình hóa tại Đại học Quốc gia Hà Nội.
  • Thiết lập hệ thống dự báo chất lượng không khí nghiệp vụ trực tuyến với lưới phân giải 3 đến 9 km cho các đô thị lớn như Hà Nội và Thành phố Hồ Chí Minh, đảm bảo thời gian dự báo trước 72 giờ và độ trễ tính toán dưới 4 giờ, do Sở Tài nguyên và Môi trường các địa phương vận hành trong vòng 24 tháng.

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

  • Học viên cao học và nghiên cứu sinh ngành Khí tượng học, Khí hậu học và Khoa học Môi trường: Nắm vững quy trình biên dịch mô hình WRF/Chem phiên bản 3.4, kỹ thuật cấu hình namelist 27 tầng khí quyển và quy trình sử dụng công cụ Prep_chem_sources để xử lý tệp phát thải NetCDF.
  • Dự báo viê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 khu vực: Khai thác bộ tham số hóa tối ưu cho khu vực Đông Nam Á để nâng cao độ chính xác dự báo nhiệt độ bề mặt với sai số dưới 1°C và theo dõi các đợt lan truyền bụi xuyên biên giới.
  • Cán bộ quản lý môi trường tại Bộ Tài nguyên và Môi trường: Tham khảo bức tranh định lượng về tác động suy giảm bức xạ 3,4 W/m2 và phân bố nồng độ bụi để xây dựng chính sách kiểm soát khí thải công nghiệp và kế hoạch hành động quốc gia về quản lý chất lượng không khí.
  • Giảng viên và kỹ sư phát triển phần mềm mô hình hóa: Sử dụng các phân tích chi tiết về cơ chế RADM2, MADE/SORGAM và GOCART để xây dựng bài giảng chuyên ngành và thiết kế hệ thống tự động hóa xử lý dữ liệu khí tượng - môi trường.

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

Mô hình WRF/Chem có ưu điểm vượt trội gì so với mô hình khí tượng WRF truyền thống? Mô hình WRF/Chem tích hợp đồng thời các quá trình động lực khí tượng và hóa học khí quyển trong cùng một bước tích phân. Nhờ đó, mô hình loại bỏ hoàn toàn sai số nội suy dữ liệu thời gian - không gian (vốn chiếm 30% đến 40% sai số ở mô hình offline), đồng thời tính toán được tác động hồi tiếp của xon khí lên bức xạ và vi vật lý mây.

Tại sao nhiệt độ mô phỏng tại miền Bắc Việt Nam giảm 3°C khi tích hợp mô-đun hóa học? Khi kích hoạt các cơ chế hóa học như RADM2 và MADE/SORGAM, mô hình tính toán sự hấp thụ và tán xạ bức xạ mặt trời của các hạt xon khí lơ lửng. Quá trình này ngăn cản năng lượng mặt trời chiếu xuống mặt đất, làm nhiệt độ lớp sát bề mặt giảm khoảng 3°C so với mô hình không có hóa học.

Dữ liệu phát thải đầu vào cho WRF/Chem tại Việt Nam được xây dựng từ những nguồn nào? Nghiên cứu khai thác hai cơ sở dữ liệu toàn cầu gồm RETRO có độ phân giải 0,5 độ với chuỗi số liệu 40 năm và EDGAR có độ phân giải 1x1 độ. Toàn bộ dữ liệu thô được chuyển đổi sang định dạng chuẩn NetCDF tương thích với lưới tính 30 km thông qua công cụ chuyển đổi Prep_chem_sources.

Trường gió và khí áp do mô hình WRF/Chem mô phỏng tại Việt Nam đạt độ chính xác ra sao? Mô hình tái tạo rất chuẩn xác trường gió tại các mực 850, 500 và 200 hPa cả về hướng và tốc độ so với dữ liệu NCEP NNRP. Trường độ cao địa thế vị nhận diện đúng vị trí các tâm áp, dù có sai số biên độ khoảng 10 hPa tại mực 850 hPa.

Phân bố nồng độ bụi tại Việt Nam trong giai đoạn thử nghiệm có đặc điểm gì nổi bật? Nồng độ bụi tập trung cao nhất ở các tầng đối lưu trên cao (500 và 200 hPa) với giá trị đạt từ 10 đến 12 x 10^12 µg/kg. Khu vực Nam Bộ và Nam Trung Bộ chịu tác động lớn nhất từ dòng vận chuyển bụi phía đông, trong khi miền Bắc duy trì nồng độ thấp chỉ từ 1 đến 2 x 10^12 µg/kg.

Kết luận

  • Đã thử nghiệm thành công mô hình trực tuyến WRF/Chem phiên bản 3.4 với độ phân giải ngang 30 km và cấu hình 27 tầng thẳng đứng cho khu vực Việt Nam.
  • Khẳng định độ tin cậy cao của mô hình trong tái hiện trường nhiệt độ bề mặt (sai số dao động ±1°C) và trường gió các tầng 850 hPa, 500 hPa, 200 hPa so với dữ liệu APHRODITE và NNRP.
  • Chứng minh rõ nét hiệu ứng hồi tiếp khí hậu - hóa học khi nồng độ xon khí làm giảm nhiệt độ bề mặt miền Bắc tới 3°C so với mô hình không hóa học.
  • Lượng hóa chi tiết sự phân tầng nồng độ 5 cấp hạt bụi với mức cực đại từ 10 đến 12 x 10^12 µg/kg tại tầng cao và xác định Nam Bộ là tâm điểm tiếp nhận dòng bụi bình lưu.
  • Đề xuất lộ trình 12 đến 24 tháng nhằm hoàn thiện cơ sở dữ liệu phát thải nội địa chi tiết, đồng thời khuyến nghị các cơ quan khí tượng và môi trường đưa ngay khung mô hình WRF/Chem vào hệ thống giám sát chất lượng không khí quốc gia.