Tổng quan nghiên cứu

Ô nhiễm không khí xuyên biên giới là một trong những thách thức môi trường cấp bách nhất hiện nay, đe dọa trực tiếp đến an ninh sinh thái và sức khỏe cộng đồng tại khu vực Đông Nam Á. Khu vực miền Bắc Việt Nam, với tọa độ địa lý trải dài từ 19°53’ Bắc đến 23°23’ Bắc, sở hữu địa hình lòng chảo phức tạp và là nơi sinh sống của hơn 19,9 triệu người tại vùng đồng bằng sông Hồng với mật độ dân số cao lên tới 1.225 người/km2 (theo số liệu thống kê dân số năm 2009). Đặc điểm tự nhiên này khiến khu vực chịu tác động kép từ các hoạt động phát thải nội tại của hơn 14.800 km2 diện tích châu thổ và sự dịch chuyển của các khối khí ô nhiễm công nghiệp từ các quốc gia láng giềng.

Vấn đề nghiên cứu trọng tâm được đặt ra là sự xuất hiện bất thường của các chất ô nhiễm như SO2, NOx và hiện tượng mưa axit tại những vùng núi cao như Lạng Sơn, Lào Cai – nơi hầu như không có các nguồn thải công nghiệp quy mô lớn. Mục tiêu cốt lõi của đề tài là ứng dụng hệ thống mô hình chất lượng không khí đa quy mô CMAQ (Community Multi-scale Air Quality) để mô phỏng cơ chế vận chuyển động lực học, tính toán định lượng mức độ và xác định tần suất ảnh hưởng của các chất ô nhiễm không khí xuyên biên giới đến miền Bắc Việt Nam.

Phạm vi nghiên cứu tập trung vào chuỗi dữ liệu phát thải năm 2010 với 4 tháng đại diện cho 4 mùa đặc trưng (tháng 1, 4, 7, 10), đồng thời đối chiếu với chuỗi số liệu quan trắc giai đoạn 2004 - 2008 tại các trạm môi trường quốc gia. Kết quả nghiên cứu cung cấp hệ thống số liệu khoa học nền tảng, làm rõ tỷ lệ đóng góp của nguồn phát thải ngoại biên vượt trên 50% đối với lắng đọng SO2 trong mùa đông, tạo tiền đề vững chắc cho việc hoạch định chính sách kiểm soát chất lượng không khí và đàm phán môi trường quốc tế.

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 2 nền tảng lý thuyết thủy khí động lực học và hóa học khí quyển chuyên sâu:

  • Thuyết khuếch tán rối khí quyển và phương trình liên tục Navier-Stokes: Mô hình hóa sự biến đổi nồng độ chất ô nhiễm theo không gian 3 chiều và thời gian, tích hợp các quá trình bình lưu hướng ngang, bình lưu hướng dọc, khuếch tán rối và lắng đọng trọng trường. Mô hình vận hành trên 6 giả thiết thủy động lực học chuẩn của Srivastava nhằm đơn giản hóa tính toán mà vẫn đảm bảo độ chính xác cao.
  • Cơ chế hóa học khí quyển CB-IV (Carbon Bond IV) và RADM2 (Regional Acid Deposition Model version 2): Tính toán chi tiết động học phản ứng pha khí, pha lỏng trong mây và sự hình thành sol khí thứ cấp (sulfate, nitrate, amoni) từ các tiền chất phát thải.

Các khái niệm then chốt được chuẩn hóa trong nghiên cứu bao gồm: bình lưu quy mô lớn (advection), khuếch tán rối lớp biên (turbulent diffusion), tốc độ quang phân (photolysis rate qua bộ xử lý JPROC), và phân loại 11 dạng bề mặt sử dụng đất (mô đun LUPROC) phục vụ tính toán lắng đọng khô.

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

Nghiên cứu ứng dụng tổ hợp các phương pháp định lượng và mô hình hóa hiện đại với quy trình xử lý dữ liệu nghiêm ngặt:

  • Nguồn dữ liệu đầu vào: Sử dụng cơ sở dữ liệu kiểm kê phát thải châu Á do Trung tâm Nghiên cứu Biến đổi Toàn cầu (FRCGC) và Cục Khoa học Kỹ thuật Biển - Địa cầu Nhật Bản (JAMSTEC) cung cấp; trường dữ liệu khí tượng tái phân tích từ Trung tâm Nghiên cứu Khí quyển Quốc gia Hoa Kỳ (NCAR) kết hợp mô hình WRF/MM5. Số liệu thực nghiệm mặt đất được thu thập liên tục từ trạm quan trắc tự động Nguyễn Văn Cừ (Hà Nội) và dữ liệu vệ tinh quan sát Trái Đất Aura-OMI của NASA.
  • Quy mô mẫu và lưới tính toán: Thiết lập miền lưới tính đa cấp bao phủ toàn bộ Đông Á và miền Bắc Việt Nam với độ phân giải ngang từ 36 km đến 12 km. Tập hợp dữ liệu bao gồm 8.760 bước thời gian giờ cho 4 tháng mô phỏng trọng điểm năm 2010, kết hợp phân tích chuỗi dữ liệu quan trắc 5 năm (2004 - 2008) tại 5 đô thị trọng điểm vùng I (Hà Nội, Hải Phòng, Hạ Long, Vinh, Huế).
  • Lý do lựa chọn phương pháp: Mô hình CMAQ kết hợp bộ giao diện MCIP (Meteorology-Chemistry Interface Processor) cho phép giải đồng thời hệ phương trình bảo toàn khối lượng và các phản ứng hóa học phi tuyến tính. Phương pháp kịch bản triệt tiêu nguồn thải (Zero-Emission Scenario) được lựa chọn để bóc tách chính xác mối quan hệ nguồn thải - nơi tiếp nhận (Source-Receptor relationship) giữa nguồn phát thải nội địa và nguồn phát thải quốc tế.

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 trên mô hình CMAQ đã mang lại 4 phát hiện mang tính bước ngoặt về bức tranh ô nhiễm không khí tại miền Bắc Việt Nam:

  • Đóng góp vượt trội của SO2 xuyên biên giới: Trong các tháng mùa đông (đặc biệt là tháng 1), lượng SO2 vận chuyển từ các khu công nghiệp phía Nam Trung Quốc sang miền Bắc Việt Nam chiếm tỷ trọng từ 50% đến 65% tổng lượng lắng đọng toàn vùng. Lượng SO2 lắng đọng đo đạc xấp xỉ ngưỡng tiêu chuẩn an toàn của Châu Âu.
  • Sự phân hóa rõ rệt theo mùa khí quyển: Vào tháng 1, dưới ảnh hưởng của hoàn lưu gió mùa Đông Bắc, nồng độ CO và NO2 tại các khu vực biên giới phía Bắc tăng đột biến, lan sâu vào khu vực đồng bằng sông Hồng. Ngược lại, vào tháng 7 khi gió mùa Tây Nam và Đông Nam chiếm ưu thế, tỷ lệ đóng góp ô nhiễm xuyên biên giới từ phương Bắc giảm xuống dưới 15%, lúc này nguồn phát thải nội sinh chiếm ưu thế trên 85%.
  • Hiện tượng ô nhiễm nền tại vùng đồi núi: Tại các khu vực trạm quan trắc vùng núi như Lạng Sơn, Lào Cai và Tam Đảo (Vĩnh Phúc), nồng độ trung bình của SO2 và NO2 duy trì ở mức cao xấp xỉ 0,030 - 0,050 mg/m3 vào mùa khô, dù mật độ phát thải cục bộ tại chỗ rất thấp.
  • Độ tương thích giữa mô phỏng và thực tế: Kết quả kiểm định mô hình tại trạm Nguyễn Văn Cừ (Hà Nội) cho thấy hệ số tương quan (R) giữa nồng độ mô phỏng và nồng độ thực đo đạt mức từ 0,50 đến 0,65 đối với các chất CO, SO2 và NO2, phản ánh đúng chu kỳ biến thiên ngày đêm của các chất ô nhiễm.

Thảo luận kết quả

Nguyên nhân chính dẫn đến sự gia tăng nồng độ các chất ô nhiễm trong mùa đông bắt nguồn từ sự kết hợp giữa các đợt gió mùa Đông Bắc tầng thấp với hoàn lưu áp cao lục địa Châu Á, vận chuyển các khối khí chứa nồng độ hạt vật chất và khí axit cao qua quãng đường hàng trăm kilomet.

Kết quả này hoàn toàn tương đồng với các nghiên cứu quốc tế trong khu vực: nghiên cứu của Shimadera tại Nhật Bản ghi nhận 33% SO2 và 69% - 82% sol khí có nguồn gốc xuyên lục địa; nghiên cứu của Azkar tại Bangladesh ghi nhận 21% - 53% CO và 35% - 67% PM2.5 vận chuyển từ phía Tây Bengal và Bắc Ấn Độ; hay nghiên cứu của Seungmin Lee tại Seoul (Hàn Quốc) khẳng định dòng bụi PM10 ở độ cao 3 km bắt nguồn từ các sa mạc và khu công nghiệp Trung Quốc.

Để biểu diễn trực quan các trường nồng độ, kết quả nghiên cứu đã tích hợp dữ liệu vào các bản đồ phân bố không gian 2D/3D dạng ô lưới (raster grid) và biểu đồ chuỗi thời gian đối chiếu đa kênh. Biểu đồ nồng độ giờ tại trạm Nguyễn Văn Cừ cho thấy sự trùng khớp về pha cực đại giữa dữ liệu mô hình CMAQ và số liệu đo trạm tự động. Ma trận nguồn thải - nơi tiếp nhận (Source-Receptor Matrix) thể hiện rõ khi kịch bản phát thải Việt Nam bằng 0 (Zero VN Emission) được kích hoạt, vệt nồng độ SO2 và NO2 cao vẫn bao phủ hơn 70% diện tích các tỉnh biên giới phía Bắc trong tháng 1.

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

Nhằm chủ động kiểm soát và giảm thiểu tác động tiêu cực của ô nhiễm không khí xuyên biên giới, luận văn đề xuất 4 nhóm giải pháp chiến lược:

  • Thiết lập mạng lưới 8 trạm quan trắc tự động chuyên sâu tại khu vực biên giới phía Bắc: Bộ Tài nguyên và Môi trường cần chủ trì lắp đặt các trạm quan trắc liên tục đo đa chỉ tiêu (SO2, NOx, CO, O3, PM2.5, PM10) tại các tỉnh Lạng Sơn, Lào Cai, Quảng Ninh, Cao Bằng và Hà Giang, hoàn thành trong giai đoạn 2025 - 2027 với mục tiêu đảm bảo tỷ lệ thu nhận dữ liệu đạt trên 95%.
  • Tích hợp hệ thống WRF-CMAQ vào nghiệp vụ dự báo ô nhiễm quốc gia: Trung tâm Dự báo Khí tượng Thủy văn Quốc gia cùng các Viện nghiên cứu chuyên ngành cần nâng cấp quy trình dự báo nghiệp vụ hàng ngày, cung cấp bản tin cảnh báo ô nhiễm không khí trước 48 đến 72 giờ với độ chính xác đạt trên 80% từ năm 2026.
  • Thúc đẩy đàm phán cơ chế hợp tác môi trường đa phương khu vực Đông Bắc Á - ASEAN: Bộ Ngoại giao phối hợp cùng Bộ Tài nguyên và Môi trường tham gia xây dựng thỏa thuận kiểm soát ô nhiễm không khí xuyên biên giới tương tự mô hình EMEP của Châu Âu hoặc mạng lưới EANET mở rộng, hướng tới mục tiêu giảm 20% lượng phát thải các tiền chất axit hóa liên vùng vào năm 2030.
  • Hoàn thiện hệ thống kiểm kê phát thải quốc gia chuẩn hóa theo mô hình SMOKE: Các Ủy ban Nhân dân tỉnh, thành phố trực thuộc đồng bằng Bắc Bộ cần tiến hành kiểm kê phát thải chi tiết (giao thông, công nghiệp, nông nghiệp, dân sinh) với độ phân giải lưới 1x1 km trước năm 2028 nhằm kiểm soát triệt để các nguồn thải nội sinh.

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

Nội dung luận văn mang lại giá trị thực tiễn và học thuật cao cho 4 nhóm đối tượng chính:

  • Cơ quan quản lý nhà nước và nhà hoạch định chính sách: Cán bộ thuộc Cục Kiểm soát ô nhiễm môi trường, Vụ Môi trường (Bộ TN&MT) và các Sở TN&MT địa phương có thể sử dụng số liệu làm căn cứ khoa học xây dựng Quy hoạch bảo vệ môi trường quốc gia và đàm phán các hiệp định môi trường quốc tế.
  • Các nhà nghiên cứu và giảng viên chuyên ngành Môi trường/Khí tượng: Giảng viên, nghiên cứu sinh và học viên cao học có thể tham khảo khung phương pháp ghép nối mô hình động lực học khí tượng (WRF/MM5) với mô hình lan truyền hóa học (CMAQ/CCTM) làm tài liệu giảng dạy và nghiên cứu mở rộng.
  • Chuyên gia quan trắc và kỹ sư mô hình hóa chất lượng không khí: Kỹ sư môi trường tại các trung tâm quan trắc có thể học hỏi quy trình xử lý dữ liệu phát thải diện/điểm qua mô đun ECIP, hiệu chỉnh điều kiện biên BCON/ICON và xử lý ảnh vệ tinh Aura-OMI.
  • Các tổ chức phi chính phủ và hiệp hội bảo vệ sức khỏe cộng đồng: Các tổ chức hoạt động vì không khí sạch có thể khai thác các kịch bản lan truyền ô nhiễm theo mùa để triển khai các chương trình truyền thông nâng cao nhận thức và khuyến cáo bảo vệ sức khỏe người dân vào các tháng cao điểm mùa đông.

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

Mô hình CMAQ là gì và tại sao lại là lựa chọn tối ưu cho nghiên cứu ô nhiễm không khí? CMAQ là hệ thống mô hình chất lượng không khí đa quy mô thế hệ thứ 3 do Cơ quan Bảo vệ Môi trường Hoa Kỳ (EPA) và NOAA phát triển. Mô hình có khả năng tích hợp đồng thời các quá trình vật lý, động lực học khí quyển và hơn 100 phản ứng hóa học phức tạp, cho phép mô phỏng lan truyền đa chất trên quy mô không gian rộng lớn từ đô thị đến toàn cầu.

Ô nhiễm không khí xuyên biên giới tác động mạnh nhất đến miền Bắc vào thời gian nào? Tác động lớn nhất diễn ra vào các tháng mùa đông (từ tháng 10 đến tháng 4 năm sau, đỉnh điểm là tháng 1). Trong giai đoạn này, các đợt gió mùa Đông Bắc hoạt động mạnh mẽ, mang theo lượng lớn khí thải SO2, NOx và bụi mịn từ lục địa phía Bắc tràn xuống miền Bắc Việt Nam.

Tỷ lệ đóng góp của nguồn thải xuyên biên giới đối với lượng SO2 tại miền Bắc là bao nhiêu? Theo kết quả mô phỏng từ mô hình CMAQ kết hợp dữ liệu kiểm kê phát thải quốc tế, lượng SO2 lắng đọng có nguồn gốc xuyên biên giới từ các nước láng giềng đóng góp trên 50% tổng lượng lắng đọng tại khu vực Bắc Bộ trong mùa đông, xấp xỉ ngưỡng tiêu chuẩn giới hạn của Châu Âu.

Làm thế nào để phân biệt giữa ô nhiễm nội sinh và ô nhiễm từ quốc gia láng giềng? Nghiên cứu áp dụng phương pháp chạy kịch bản đối chứng trên mô hình CMAQ. Bằng cách thiết lập kịch bản phát thải nội địa bằng 0 (Zero VN Emission) và kịch bản phát thải nước ngoài bằng 0 (Zero China Emission), các nhà khoa học có thể bóc tách và định lượng chính xác nồng độ các chất ô nhiễm do từng nguồn đóng góp.

Mức độ tương thích giữa kết quả tính toán của mô hình CMAQ và số liệu thực tế ra sao? Kết quả kiểm định giữa dữ liệu mô phỏng từ mô hình và số liệu quan trắc thực nghiệm tại trạm tự động Nguyễn Văn Cừ (Hà Nội) cho thấy hệ số tương quan R dao động trong khoảng 0,50 đến 0,65 đối với các khí CO, SO2 và NO2, chứng minh mô hình phản ánh sát thực tế xu thế biến động chất lượng không khí.

Kết luận

  • Ứng dụng thành công hệ thống mô hình CMAQ kết hợp mô hình khí tượng để mô phỏng toàn diện cơ chế lan truyền các chất ô nhiễm không khí (SO2, NO2, CO) đến khu vực miền Bắc Việt Nam.
  • Khẳng định nguồn phát thải xuyên biên giới đóng góp trên 50% lượng lắng đọng SO2 tại Bắc Bộ trong mùa đông, làm sáng tỏ nguyên nhân gây mưa axit tại các tỉnh miền núi phía Bắc.
  • Xác định rõ tính quy luật mùa của các chất ô nhiễm với sự chi phối áp đảo của gió mùa Đông Bắc vào tháng 1 và sự suy giảm ảnh hưởng ngoại biên vào tháng 7.
  • Cung cấp phương pháp luận chuẩn xác về việc thiết lập các kịch bản triệt tiêu nguồn thải để đánh giá ma trận nguồn thải - nơi tiếp nhận trong quản lý môi trường.
  • Đặt nền tảng dữ liệu khoa học tin cậy phục vụ công tác dự báo chất lượng không khí nghiệp vụ và xây dựng chiến lược ngoại giao môi trường của Việt Nam giai đoạn 2025 - 2030.

Luận văn đã đóng góp một công trình khoa học toàn diện, giải quyết trọn vẹn bài toán định lượng ô nhiễm không khí xuyên biên giới tại Việt Nam bằng công cụ mô hình hóa hiện đại. Để tiếp tục hoàn thiện và phát huy các kết quả nghiên cứu, việc cập nhật cơ sở dữ liệu phát thải mới nhất và mở rộng mạng lưới trạm quan trắc tự động liên tục tại các tỉnh biên giới là bước đi tất yếu trong giai đoạn tiếp theo. Các nhà quản lý, nhà nghiên cứu và bạn đọc quan tâm có thể tìm đọc chi tiết toàn văn luận văn tại Thư viện Đại học Quốc gia Hà Nội để cùng chung tay thúc đẩy các giải pháp bảo vệ bầu không khí trong lành cho cộng đồng.