Bối cảnh và vấn đề nghiên cứu

Ô nhiễm không khí do bụi lơ lửng, đặc biệt là các hạt bụi có đường kính khí động học nhỏ hơn hoặc bằng 10 micromet ($\text{PM}{10}$), đang là một trong những thách thức môi trường lớn tại các đô thị công nghiệp hóa nhanh. Do kích thước nhỏ, $\text{PM}{10}$ có khả năng thâm nhập sâu vào đường hô hấp và hệ tuần hoàn của con người, làm gia tăng nguy cơ mắc các bệnh tim mạch, hô hấp cũng như tỷ lệ nhập viện và tử vong trong cộng đồng.

Sự biến thiên của nồng độ $\text{PM}_{10}$ trong khí quyển chịu sự chi phối phức tạp bởi các nguồn phát thải nhân sinh và các yếu tố khí tượng thủy văn (nhiệt độ, độ ẩm tương đối, tốc độ gió) cũng như cấu trúc động lực học của tầng khí quyển sát mặt đất, đại diện bởi chiều cao lớp biên hành tinh (Planetary Boundary Layer Height - PBLH). Lớp biên khí quyển là nơi diễn ra các quá trình phát tán, xáo trộn rối và tích tụ chất ô nhiễm.

Thành phố Đào Viên (Taoyuan) là một đô thị trực thuộc trung ương ở phía tây bắc Đài Loan với diện tích $1.221\text{ km}^2$ và dân số đạt $2.116.988$ người (năm 2016). Đào Viên đã chuyển mình từ một khu vực nông nghiệp thành trung tâm kinh tế, công nghiệp và thương mại trọng điểm phía bắc Đài Loan. Sự gia tăng nhanh chóng về dân số, lưu lượng giao thông đường bộ và các hoạt động vận tải hàng không nhộn nhịp tại Sân bay Quốc tế Đào Viên đã làm tăng tiêu thụ năng lượng và phát thải hạt vào bầu khí quyển.

Trước thực tế đó, khóa luận tốt nghiệp "The correlations between particulate matter concentrations, planetary boundary layer height and meteorological parameters" (Mối tương quan giữa nồng độ bụi mịn, chiều cao lớp biên và các thông số khí tượng) được tác giả Đỗ Minh Hồng thực hiện nhằm giải quyết 2 câu hỏi/nhiệm vụ nghiên cứu chính:

  1. Xác định hàm lượng và đặc điểm biến thiên (theo giờ, theo tháng và theo mùa) của nồng độ bụi $\text{PM}_{10}$ tại khu vực đô thị thành phố Đào Viên trong năm 2015.
  2. Phân tích định lượng mức độ ảnh hưởng và mối tương quan giữa chiều cao lớp biên hành tinh (PBLH) cùng các yếu tố khí tượng (nhiệt độ, độ ẩm, tốc độ gió) đối với sự tích tụ và phân tán nồng độ $\text{PM}_{10}$.

Đối tượng và phạm vi nghiên cứu:

  • Đối tượng nghiên cứu: Nồng độ bụi mịn $\text{PM}_{10}$, chiều cao lớp biên hành tinh (PBLH), các thông số khí tượng bề mặt (nhiệt độ không khí, độ ẩm tương đối, tốc độ gió).
  • Phạm vi không gian: Trạm quan trắc chất lượng không khí Đào Viên (kinh độ $121.99^\circ\text{E}$), thuộc Mạng lưới Quan trắc Chất lượng Không khí Đài Loan (TAQMN).
  • Phạm vi thời gian: Dữ liệu chuỗi thời gian liên tục 1 năm, từ ngày 01/01/2015 đến ngày 31/12/2015.

Cơ sở lý thuyết và phương pháp

Khung lý thuyết và khái niệm khoa học

  • Bụi mịn $\text{PM}_{10}$: Hỗn hợp các hạt rắn và giọt chất lỏng lơ lửng trong không khí có đường kính khí động học từ $2,5\text{ }\mu\text{m}$ đến $10\text{ }\mu\text{m}$ (hạt thô hít thở được) và $\le 2,5\text{ }\mu\text{m}$ (hạt mịn) theo định nghĩa của Cơ quan Bảo vệ Môi trường Hoa Kỳ (US EPA). Nguồn gốc bao gồm bụi khoáng, khí thải nhà máy điện, phương tiện giao thông, sol khí biển và muối biển. Chu trình ô nhiễm gồm ba giai đoạn: phát thải tại nguồn, vận chuyển - khuếch tán trong khí quyển, và tiếp nhận bởi sinh vật.
  • Tiêu chuẩn chất lượng không khí: Quy chuẩn của Cơ quan Bảo vệ Môi trường Đài Loan (Taiwan EPA) quy định nồng độ $\text{PM}_{10}$ trung bình 24 giờ tối đa là $125\text{ }\mu\text{g/m}^3$ và nồng độ trung bình năm là $65\text{ }\mu\text{g/m}^3$.
  • Cơ chế tác động của vi khí hậu:
    • Gió: Tốc độ gió điều tiết quá trình pha loãng và xáo trộn rối. Gió mạnh làm tăng xáo trộn không khí giúp phân tán bụi, trong khi gió lặng làm bụi tích tụ cục bộ.
    • Nhiệt độ: Tác động đến tốc độ phản ứng hóa học hình thành sol khí thứ cấp và áp suất hơi bão hòa của các hợp chất bán bay hơi (như amoni nitrat).
    • Độ ẩm: Chi phối quá trình ngưng kết, tăng kích thước hạt và thúc đẩy quá trình lắng đọng ướt (wet deposition) qua giáng thủy.
  • Lớp biên hành tinh (Planetary Boundary Layer - PBL): Tầng thấp nhất của khí quyển đối lưu chịu tác động trực tiếp của ma sát bề mặt Trái Đất. Độ dày PBL thay đổi theo ngày đêm (thấp vào ban đêm do giảm dòng đối lưu nhiệt từ bề mặt; cao vào ban ngày do đối lưu phát triển mạnh) và theo mùa (thấp vào mùa lạnh, cao vào mùa nóng). Hầu hết các chất ô nhiễm phát thải từ bề mặt bị giam giữ bên trong lớp biên này.

Nguồn dữ liệu và trang thiết bị thu thập

  1. Dữ liệu chất lượng không khí và khí tượng bề mặt:
    • Thu thập từ trạm quan trắc Đào Viên thuộc Mạng lưới Quan trắc Chất lượng Không khí Đài Loan (TAQMN) - hệ thống gồm 76 trạm quan trắc tự động trên toàn lãnh thổ Đài Loan (60 trạm tổng hợp, 5 trạm công nghiệp, 2 trạm công viên quốc gia, 4 trạm nền, 6 trạm giao thông và 2 trạm khác).
    • Thiết bị quan trắc $\text{PM}_{10}$: Máy đo nồng độ bụi tự động VEREWA F701 (DURAG Group). Thiết bị hoạt động dựa trên nguyên lý suy giảm hấp thụ tia beta (electron) phát ra từ nguồn phóng xạ đi qua lớp bụi thu gom trên băng lọc; độ khả dụng đạt $>95%$, giới hạn phát hiện $<1\text{ }\mu\text{g/m}^3$. Thiết bị được lắp đặt trên mái công trình ở độ cao $14,5\text{ m}$ so với mặt đất và tự động hiệu chuẩn điểm 0 (zero check) cùng khoảng đo (span check) hàng ngày.
    • Các thông số khí tượng thu thập đồng thời tại trạm: Nhiệt độ không khí ($^\circ\text{C}$), Độ ẩm tương đối ($%$), Tốc độ gió ($\text{m/s}$).
  2. Dữ liệu chiều cao lớp biên khí quyển (PBLH):
    • Trích xuất từ Hệ thống Mô hình hóa và Đồng hóa Dữ liệu Toàn cầu (GMAO) của NASA thông qua dòng sản phẩm Goddard Earth Observing System-Forward Processing (GEOS-FP phiên bản 5.1).
    • Dữ liệu có độ phân giải ngang $0,3125^\circ$ kinh độ và độ phân giải thời gian 1 giờ dạng lưới lồng phủ khắp Đài Loan.

Phương pháp xử lý và phân tích dữ liệu

  • Tổng hợp và sàng lọc được $8.489$ giá trị đo theo giờ hợp lệ của $\text{PM}_{10}$, PBLH và các biến khí tượng trong suốt năm 2015.
  • Sử dụng phần mềm MATLAB để giải mã, xử lý các tệp dữ liệu viễn thám định dạng HDF4 từ GMAO và tính toán thống kê.
  • Áp dụng phân tích tương quan tuyến tính Pearson ($r$) và hệ số xác định ($R^2$) để đánh giá mức độ liên kết giữa nồng độ $\text{PM}_{10}$ và từng gradient của các thông số khí tượng cũng như chiều cao lớp biên:

$$\text{Pearson } r = \frac{\sum (X_i - \bar{X})(Y_i - \bar{Y})}{\sqrt{\sum (X_i - \bar{X})^2 \sum (Y_i - \bar{Y})^2}}$$


Nội dung chính theo từng chương

Phần I: Đặt vấn đề và Mục tiêu nghiên cứu (Research Rationale and Objectives)

  • Trình bày tổng quan về tác hại sức khỏe của bụi mịn $\text{PM}_{10}$, sự cần thiết của việc hiểu rõ hành vi động lực học của bụi trong mối tương quan với các biến số khí tượng và lớp biên khí quyển tại khu vực đô thị công nghiệp Đào Viên.
  • Thiết lập 2 mục tiêu nghiên cứu cụ thể và nêu rõ giới hạn nghiên cứu khi chỉ sử dụng chuỗi số liệu từ một trạm quan trắc đại diện.

Phần II: Tổng quan tài liệu (Literature Review)

  • Hệ thống hóa các kiến thức nền tảng về bụi $\text{PM}_{10}$, cơ chế vận chuyển khí quyển và các quy chuẩn môi trường của US EPA và Taiwan EPA.
  • Phân tích chi tiết các cơ chế vật lý - hóa học khí quyển: vai trò pha loãng của gió, ảnh hưởng của nhiệt độ tới nồng độ amoni nitrat (ở nhiệt độ ban đầu $>17^\circ\text{C}$ hoặc $<17^\circ\text{C}$), vai trò của độ ẩm trong quá trình lắng đọng ướt.
  • Trình bày cấu trúc lớp biên hành tinh (PBL), cơ chế nâng cao/hạ thấp tầng xáo trộn và tác động đến sự giam giữ sol khí bề mặt.
  • Giới thiệu mạng lưới TAQMN Đài Loan, cơ cấu các sản phẩm dữ liệu mô hình GMAO/GEOS-FP của NASA (GEOS-5 FP, GEOS-5 FP-IT, MERRA-2, 7km-G5NR, SMAP L4) và cơ sở toán học của phép phân tích tương quan Pearson cùng mô hình hồi quy bình phương tối thiểu trên nền tảng MATLAB.

Phần III: Vật liệu và Phương pháp nghiên cứu (Materials and Methods)

  • Cung cấp bức tranh tổng thể về địa lý, khí hậu cận nhiệt đới ẩm của Đào Viên (nằm cách Đài Bắc $40\text{ km}$ về phía tây nam, diện tích $1.221\text{ km}^2$, khí hậu đặc trưng mùa đông ấm/ôn hòa, mùa hè nóng ẩm).
  • Mô tả chi tiết thông số kỹ thuật của trạm quan trắc Đào Viên, thiết bị đo bụi VEREWA F701 đặt ở độ cao $14,5\text{ m}$.
  • Thiết lập sơ đồ quy trình nghiên cứu (Flowchart): Nhập dữ liệu 1 năm (2015) từ TAQMN và GMAO $\rightarrow$ Xử lý định dạng trên MATLAB $\rightarrow$ Phân tích chuỗi thời gian, biến thiên ngày đêm và biến thiên mùa $\rightarrow$ Phân tích hệ số tương quan Pearson $\rightarrow$ Đánh giá tương quan giữa $\text{PM}_{10}$ và các biến số.

Phần IV: Kết quả và Thảo luận (Results and Discussion)

  • Đặc trưng khí tượng và lớp biên: Nhiệt độ trung bình năm đạt $23^\circ\text{C}$ (tháng 7 nóng nhất, cực đại $38^\circ\text{C}$); Độ ẩm tương đối trung bình $78%$ (dao động tháng từ $74%$ đến $82%$); Tốc độ gió trung bình thấp, với hơn $50%$ số lần quan trắc có vận tốc $<1\text{ m/s}$.
  • Biến thiên của PBLH: PBLH trung bình ban đêm duy trì khoảng $700\text{ m}$, bắt đầu dâng cao từ 06:00 sáng, đạt đỉnh cực đại trên $900\text{ m}$ (trung bình mùa đông đạt đỉnh cao nhất $1.255\text{ m}$ do gió mùa đông bắc mạnh thúc đẩy xáo trộn cơ học) vào lúc 12:00 trưa, sau đó hạ dần về mức nền ban đêm vào khoảng 21:00.
  • Biến thiên nồng độ $\text{PM}_{10}$: Nồng độ giờ dao động từ $1$ đến $148\text{ }\mu\text{g/m}^3$, giá trị trung bình là $40\text{ }\mu\text{g/m}^3$; nồng độ phân vị thứ 90 nhỏ hơn $71\text{ }\mu\text{g/m}^3$. Nồng độ trung bình ngày nằm trong khoảng $10 - 104\text{ }\mu\text{g/m}^3$, đạt tiêu chuẩn chất lượng không khí của Đài Loan ($<125\text{ }\mu\text{g/m}^3$).
  • Quy luật biến thiên thời gian của $\text{PM}_{10}$:
    • Theo mùa: Nồng độ đạt đỉnh vào mùa đông (tháng 1 cao nhất: $53\text{ }\mu\text{g/m}^3$) và thấp nhất vào mùa hè (tháng 8 thấp nhất: $28\text{ }\mu\text{g/m}^3$). Mùa hè có mưa lớn kéo dài (từ tháng 5 đến tháng 8) giúp rửa trôi bụi hiệu quả qua lắng đọng ướt; mùa đông chịu ảnh hưởng của gió mùa đông bắc vận chuyển bụi từ vùng trung tâm châu Á tới.
    • Theo ngày đêm: Nồng độ thấp nhất lúc 06:00 - 07:00 sáng (giao thông thưa thớt), xuất hiện đỉnh chính vào giờ cao điểm buổi sáng (09:00) và đỉnh phụ vào giờ cao điểm buổi chiều (17:00).

Kết quả và đóng góp

Kết quả định lượng chi tiết

Nghiên cứu đã xác lập các chỉ số tương quan thống kê cụ thể giữa nồng độ $\text{PM}_{10}$ với chiều cao lớp biên và các thông số khí tượng tại thành phố Đào Viên trong năm 2015:

Yếu tố phân tích Điều kiện / Phân tầng Hệ số tương quan Pearson ($r$) Hệ số xác định ($R^2$)
Chiều cao lớp biên (PBLH) Toàn bộ quan trắc $-0,48$ $0,23$
Tầng sát mặt đất ($\text{PBLH} < 900\text{ m}$) $-0,84$ $0,71$
Tầng trên cao ($\text{PBLH} > 900\text{ m}$) $-0,32$ $0,10$
Nhiệt độ không khí Toàn bộ chuỗi dữ liệu $-0,75$ $0,56$
Độ ẩm tương đối Toàn bộ chuỗi dữ liệu $-0,92$ $0,85$
Tốc độ gió Gió nhẹ ($< 2\text{ m/s}$) Tương quan nghịch chặt chẽ Tích tụ ô nhiễm cục bộ
Gió mạnh ($> 4\text{ m/s}$) $-0,68$ Pha loãng / Vận chuyển xa

Đóng góp của nghiên cứu

  • Đóng góp học thuật: Làm rõ cơ chế phi tuyến tính của chiều cao lớp biên hành tinh đối với sự phân tán bụi $\text{PM}{10}$. Tác giả chứng minh rằng mối tương quan giữa PBLH và $\text{PM}{10}$ không đồng nhất trên toàn bộ cột khí quyển: ở độ cao dưới $900\text{ m}$ gần mặt đất, tương quan nghịch diễn ra rất mạnh ($r = -0,84$), phản ánh rõ sự nén ép nồng độ ô nhiễm khi tầng xáo trộn bị thu hẹp; trong khi ở độ cao trên $900\text{ m}$, tương quan trở nên yếu ($r = -0,32$).
  • Đóng góp thực tiễn: Khẳng định chất lượng không khí Đào Viên năm 2015 nằm trong ngưỡng quy chuẩn quốc gia của Đài Loan ($10 - 104\text{ }\mu\text{g/m}^3$ so với giới hạn $125\text{ }\mu\text{g/m}^3$), đồng thời chỉ ra các giai đoạn rủi ro ô nhiễm cao nhất (mùa đông và các khung giờ cao điểm giao thông 09:00, 17:00) để cơ quan quản lý có phương án cảnh báo.

Hạn chế và hướng nghiên cứu tiếp

  • Hạn chế đã được tác giả chỉ rõ: Toàn bộ dữ liệu thực nghiệm nồng độ bụi $\text{PM}_{10}$ và khí tượng bề mặt chỉ được thu thập từ một trạm quan trắc đơn lẻ (trạm Đào Viên). Điều này làm hạn chế tính đại diện không gian của kết quả nghiên cứu đối với toàn bộ vùng đô thị Đào Viên vốn có địa hình trải dài từ vùng núi Tuyết Sơn ra eo biển Đài Loan.
  • Hướng nghiên cứu tiếp: Mở rộng tích hợp số liệu từ nhiều trạm quan trắc khác nhau trong mạng lưới TAQMN (trạm công nghiệp, trạm giao thông, trạm nền); bổ sung các thông số ô nhiễm quan trọng khác như $\text{PM}_{2.5}$ và các hợp chất dạng khí; kết hợp chuỗi dữ liệu nhiều năm để phân tích xu thế dài hạn dưới tác động của biến đổi khí hậu.

Giá trị tham khảo

  • Đối tượng thụ hưởng: Khóa luận là tài liệu tham khảo cho sinh viên, học viên cao học và nhà nghiên cứu thuộc các chuyên ngành Khoa học Môi trường, Quản lý Môi trường, Khí tượng - Khí hậu học, và Kỹ thuật Môi trường.
  • Nội dung giá trị:
    • Phương pháp trích xuất, xử lý dữ liệu đồng hóa khí quyển viễn thám toàn cầu (GMAO/GEOS-FP của NASA) kết hợp dữ liệu quan trắc trạm bề mặt thông qua MATLAB.
    • Phương pháp phân đoạn dữ liệu lớp biên hành tinh ($<900\text{ m}$ và $>900\text{ m}$) trong phân tích tương quan thống kê môi trường.
    • Cơ sở dữ liệu thực nghiệm chi tiết về khí hậu và chất lượng không khí tại khu vực đô thị công nghiệp Đài Loan.

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

1. Nồng độ $\text{PM}_{10}$ trung bình và phạm vi biến thiên tại trạm Đào Viên trong năm 2015 được ghi nhận như thế nào?

Theo kết quả phân tích $8.489$ giá trị đo theo giờ trong năm 2015, nồng độ $\text{PM}_{10}$ trung bình năm là $40\text{ }\mu\text{g/m}^3$, giá trị theo giờ dao động từ $1$ đến $148\text{ }\mu\text{g/m}^3$, với $90%$ số liệu ghi nhận dưới mức $71\text{ }\mu\text{g/m}^3$. Nồng độ trung bình ngày nằm trong khoảng từ $10$ đến $104\text{ }\mu\text{g/m}^3$, không vượt quá quy chuẩn chất lượng không khí 24 giờ của Đài Loan ($125\text{ }\mu\text{g/m}^3$).

2. Nguyên nhân nào khiến nồng độ $\text{PM}_{10}$ tại Đào Viên đạt mức thấp nhất vào mùa hè và cao nhất vào mùa đông?

Vào mùa hè, lượng mưa giáng thủy tăng đột biến và kéo dài từ tháng 5 đến tháng 8 do gió mùa tây nam, giúp loại bỏ các hạt bụi lơ lửng thông qua cơ chế lắng đọng ướt (wet deposition), làm nồng độ bụi giảm xuống mức thấp nhất vào tháng 8 ($28\text{ }\mu\text{g/m}^3$). Ngược lại, vào mùa đông, gió mùa đông bắc mang theo không khí lạnh cùng lượng bụi và chất ô nhiễm vận chuyển tầm xa từ lục địa Trung Á đến Đài Loan, đồng thời các điều kiện nhiệt độ thấp làm hạn chế phát tán, dẫn đến nồng độ bụi đạt đỉnh cao nhất vào tháng 1 ($53\text{ }\mu\text{g/m}^3$).

3. Mối tương quan giữa $\text{PM}_{10}$ và chiều cao lớp biên hành tinh (PBLH) có đặc điểm gì đặc biệt?

Tương quan tuyến tính Pearson trên toàn bộ tập dữ liệu cho thấy mối liên hệ nghịch ở mức yếu ($r = -0,48; R^2 = 0,23$). Tuy nhiên, khi phân tầng phân tích theo độ cao, trong dải lớp biên sát mặt đất ($\text{PBLH} < 900\text{ m}$), $\text{PM}_{10}$ có tương quan nghịch rất chặt chẽ ($r = -0,84; R^2 = 0,71$). Khi $\text{PBLH} > 900\text{ m}$, mức độ tương quan giảm rõ rệt ($r = -0,32; R^2 = 0,10$).

4. Thiết bị nào được sử dụng để đo nồng độ $\text{PM}_{10}$ tại trạm Đào Viên và hoạt động theo nguyên lý nào?

Trạm quan trắc sử dụng máy đo nồng độ bụi tự động VEREWA F701 do DURAG Group chế tạo, được đặt ở độ cao $14,5\text{ m}$ trên mái tòa nhà. Thiết bị hoạt động dựa trên nguyên lý đo độ suy giảm hấp thụ tia beta (electron) phát ra từ nguồn phóng xạ khi đi qua lớp bụi đã được gom lại trên màng lọc trong một chu kỳ xác định, với giới hạn phát hiện dưới $1\text{ }\mu\text{g/m}^3$ và độ khả dụng đạt trên $95%$.

5. Dữ liệu chiều cao lớp biên hành tinh (PBLH) được trích xuất từ nguồn nào?

Dữ liệu PBLH được lấy từ hệ thống mô hình đồng hóa dữ liệu Goddard Earth Observing System-Forward Processing (GEOS-FP phiên bản 5.1 DAS) do Văn phòng Mô hình hóa và Đồng hóa Toàn cầu (GMAO) thuộc NASA cung cấp trực tuyến, với độ phân giải không gian $0,3125^\circ$ kinh độ và độ phân giải thời gian 1 giờ bao phủ khu vực Đài Loan.


Kết luận

Khóa luận tốt nghiệp của tác giả Đỗ Minh Hồng đã phân tích toàn diện chuỗi số liệu thực nghiệm năm 2015 để làm rõ quy luật biến thiên của nồng độ bụi $\text{PM}{10}$ tại thành phố Đào Viên, Đài Loan. Kết quả nghiên cứu khẳng định sự suy giảm rõ rệt của nồng độ bụi mịn dưới tác động của độ ẩm cao ($r = -0,92$) và nhiệt độ tăng ($r = -0,75$), đồng thời phát hiện mối tương quan nghịch rất mạnh giữa $\text{PM}{10}$ và chiều cao lớp biên hành tinh ở tầng sát mặt đất dưới $900\text{ m}$ ($r = -0,84$). Công trình cung cấp cơ sở khoa học và phương pháp luận hữu ích cho việc kết hợp dữ liệu mô hình toàn cầu với quan trắc bề mặt trong công tác quản lý và đánh giá chất lượng môi trường không khí đô thị.