GIS-Integrated Analysis of Spatial and Vertical Variability in Airborne Pollen and Fungal Spore Concentrations Across an Urban Campus in Islamabad, Pakistan

기술 요약: 파키스탄 이슬라마바드 도시 캠퍼스 내 항공 화분 및 곰팡이 포자 농도의 공간적 및 수직적 변이성에 대한 GIS 통합 분석

문제 제기
항공 생물 에어로졸, 특히 화분과 곰팡이 포자는 도시 대기 질의 핵심 구성 요소이며 호흡기 및 알레르기 질환의 주요 원인이다. 파키스탄의 에어로바이올로지(aerobiology) 모니터링은 역사적으로 단일 스테이션 Burkard 트랩 시스템에 의존해 왔으나, 이러한 설정은 도시 환경에 내재된 미세한 공간적 불균질성을 포착하지 못한다. 또한, 기존 연구들은 종종 수직 농도 구배에 대한 정밀한 데이터와 국지적 토지 이용/토지 피복(LULC) 특성의 구체적인 영향력을 결여하고 있다. 특히 남아시아 도시들을 위한 현장 모니터링 데이터와 고해 resolution 모델링이 매우 부족한 실정이다. 본 연구는 단일 도시 캠퍼스 내 생물 에어로졸의 공간적 및 수직적 변이성에 대한 이해의 공백을 메우고자 하며, 단일 스테이션 모니터링이 충분한지, 그리고 국지적 환경 요인이 농도 차이를 어떻게 유발하는지 규명하는 것을 목표로 한다.

연구 방법론
본 연구는 2016년 4월 4일부터 20일까지 봄철 최성기 분산 기간 동안 파키스탄 이슬라마바드에 위치한 국립 과학기술대학교(NUST) 캠퍼스에서 수행되었다.

• 샘플링 설계: 다양한 토지 피복(입구 게이트, 개방된 공간, 주거 지역, 건설 현장, 건물 내부)을 대표하도록 3.28 km² 캠서스 내에 6개의 샘플링 지점을 전략적으로 선정하였다. 샘플링은 세 가지 서로 다른 수직 높이(지면 높이 0.91 m, 창문 높이 7.66 m, 옥상 높이 15.3 m)에서 수행되었다.
• 장비: 표준 장비(예: Burkard 트랩)의 비용 제약을 극복하기 위해, 저자는 저비용 휴대용 능동 충돌 샘플러(단가 약 미화 20달러)를 개발하여 배치하였다. 이 장치는 입자들을 관성 충돌을 통해 포집하기 위해 실리콘 겔이 코팅된 유리 막을 가진 회전 헤드를 사용하며, 10%의 가동률(duty cycle)로 작동한다. 총 144개의 샘-플이 수집되었다.
• 분석:
  • 현미경 관찰: 직접 현미경 관찰법을 통해 33종의 화분 유형과 14종의 곰팡이 속(genera)을 식별하였다.
  • 공간 분석: GIS를 사용하여 각 지점 주변 200 m 버퍼 내의 LULC 변수(자연 식생, 초지, 황무지, 건물 면적, 수목 수)를 정량화하였다. 누적 생물 에어로졸 농도와 LULC 변수 간의 스피어먼 서열 상관관계(Spearman rank correlations)를 계산하였다.
  • 기상 모델링: 기상 데이터(온도, 습도, 풍속, 풍향 등)는 캠퍼스 내 단일 기상 관측소로부터 얻었다. 풍향은 직교 성분(sin 및 cos)으로 분해되었다. Leave-One-Out 교차 검증(LOO-CV)을 사용하여 보통 최소 제곱(OLS) 및 랜덤 포레스트(RF) 회귀 모델을 적용하였다.
  • 해석 가능성: 비선형 패턴을 해석하기 위해, 연구진은 기상 동인의 방향성 효과를 결정하기 위해 순열 중요도(Permutation Importance)와 SHAP(SHapley Additive exPlanations)을 활용하였다.
  • 수직 구배: 농도 변화를 정량화하기 위해 멱함수 쇠퇴 모델($C = a \times h^{-b}$)을 적용하였으며, 윌콕슨 부호 순위 검정(Wilcoxon signed-rank tests)과 효과 크기(Cohen's d) 계산을 통해 보완하였다.

주요 결과

• 공간적 불균질성: 상당한 도시 내 변이성이 관찰되었다. 화분 농도는 약 1 km 이내에 위치한 지점들 사이에서 5배(377 ~ 1,985 grains m⁻³) 차이가 났으며, 곰팡이 농도는 3배(1,166 ~ 3,234 grains m⁻³) 차이를 보였다. 가장 높은 화분 수는 개방된 공간에서 나타났고, 가장 높은 곰팡이 부하는 수목 밀도가 높은 입구 게이트에서 발견되었다.
• 분류학적 우점도: Aspergillus와 Alternaria가 지배적인 곰팡이 분류군이었으며, 이들이 전체 곰팡이 검출량의 70~80%를 차지했다.
• LULC 동인: 스피어먼 상관관계 분석(n=6) 결과, 자연 식생 비율은 Aspergillus와 Alternaria 모두에 대해 유의미한 양의 예측 변수였다(ρ = 0.829, p = 0.042). 반대로, 황무지는 총 곰팡이 농도에 대해 유의미한 음의 예측 변수였다(ρ = −0.880, p = 0.021).
• 기상 동인: SHAP 분석은 서로 다른 곰팡이 분류군에 대한 뚜렷한 풍향 의존성을 확인하였다. Aspergillus 농도는 동풍 성분(cos(WD))과 음의 상관관계를 보였고, Alternaria는 남풍 성분(sin(WD))과 양의 상관관계를 보였다. 온도는 화분 농도와 음의 상관관계를 보였다.
• 모델 성능: 단일 기상 기록을 모든 지점에 적용하여 공간적 분산을 포착하지 못한 실험 설계로 인해, LOO-CV 하에서 OLS와 RF 모델 모두 음(-)의 $R^2$ 값을 나타냈다. 그러나 OLS가 RF보다 약간 더 나은 성능을 보였는데, 이는 해당 시간 척도에서의 선형적 관계를 시사한다.
• 수직 구배: 생물 에어로졸 농도는 높이에 따라 현저히 감소하였다. 평균 화분 농도는 지면에서 옥상 높이 사이에서 77% 감소하였고, 곰팡이 농도는 60% 감소하였다. 멱함수 쇠퇴 모델이 잘 적합되었다(예: Aspergillus $R^2 = 0.95$). 효과 크기(Cohen's d)는 매우 컸으며(1.13–1.96), 이는 옥상 모니터링이 인간의 호흡 구역 노출을 상당히 과소평가함을 나타낸다.

의의 및 주장
본 논문은 이 지역의 곰팡이 분류군에 대한 종 수준의 풍향 분산(divergence)을 세 가지 독립적인 중요도 측정치(MDI, 순열 중요도, SHAP)를 통해 확인한 최초의 정량화 사례임을 주장한다. 본 연구는 단일 스테이션 모니터링이 국지적 LULC와 미세 지형에 의해 상당한 농도 구배가 발생하는 도시 내 생물 에어로졸의 불균질성을 특징짓기에 불충분하다는 것을 입증한다.

주요 기여 중 하나는 표준 트랩 대비 98%의 비용 절감을 달एको 한 저비용 휴대용 샘플링 방법의 검증이며, 이는 자원이 제한된 남아시아 환경을 위한 복제 가능한 템플릿을 제공한다. 또한, 본 연구는 보안 및 방해 없는 노출을 위해 흔히 사용되는 옥상 모니터링 스테이션이 지면의 알레르기 유발 물질 노출을 60~77% 과소평가한다는 정량적 근거를 제공한다. 저자들은 인간의 노출 위험을 정확하게 평가하기 위해 모니터링 스테이션이 지면 높이(0.9–1.5 m)에 배치되어야 한다고 결론짓는다. 재현과 확장을 용이하게 하기 위해 모든 데이터, 코드 및 Jupyter 노트북은 공개적으로 제공된다.