산불 발생 시 장기 PM2.5 노출이 사망 위험을 크게 높인다

초록

본 연구는 2007‑2016년 미국 6천만 명 이상의 65세 이상 Medicare 수혜자를 대상으로, 연간 전체 소스 PM2.5와 산불 연관 PM2.5 일수의 상호작용이 전체 사망률에 미치는 영향을 분석하였다. 연간 평균 PM2.5가 현재 NAAQS 기준(9 µg/m³)보다 높을 때, 연간 비산불 PM2.5 일수가 35일 이상인 지역에서는 사망 위험이 현저히 증가함을 확인했다. 또한 빈곤율이 높은 ZIP 코드와 지역별 차이도 존재한다는 점을 제시한다.

상세 분석

이 연구는 대규모 관찰코호트를 활용해 장기 PM2.5 노출과 산불 연관 미세먼지( wildfire PM2.5 ) 일수 간의 효과 수정(effect modification)을 정량화하였다. 주요 노출 변수는 Di et al.의 일일 예측 모델을 ZIP 코드 수준으로 평균한 연간 전체 소스 PM2.5이며, 산불 일수는 Childs et al.의 10 km² 격자 데이터를 ZIP 코드에 면적 가중 평균한 후 0 µg/m³ 초과 일수를 연간 합산해 0‑20, 21‑35, 36 일 이상 세 구간으로 구분하였다.

통계적 접근은 각 구간별로 Poisson 회귀모델을 적용해 사망률 위험비(HR)를 추정했으며, 자연 스플라인을 이용해 비선형 노출‑반응 관계를 허용하였다. 모델은 개인 수준(연령, 성별, 인종·민족, Medicaid 가입 여부)과 지역·시간 고정효과, 그리고 10개의 ZIP/카운티 수준 사회경제·인구학적 변수와 4개의 기상 변수(여름·겨울 평균 최고 기온·습도)를 동시에 조정하였다. 이러한 다층 조정은 혼란 변수에 의한 편향을 최소화하고, 실제 PM2.5와 산불 연관 PM2.5의 독립적·상호작용 효과를 분리하는 데 기여한다.

결과적으로, 전체 코호트에서 연간 평균 PM2.5가 12 µg/m³에 도달했을 때 HR은 약 1.04(95 % CI 1.02‑1.06)였으며, 산불 일수가 0‑20일 구간에서는 비슷한 수준을 보였다. 그러나 21‑35일 구간에서는 HR이 약 1.06으로 상승했고, 36일 이상 구간에서는 HR이 1.12 이상으로 급격히 증가하였다. 이는 산불 연관 미세먼지 노출이 장기 PM2.5와의 독립적 위험을 넘어, PM2.5에 대한 민감도를 증폭시킨다는 강력한 증거이다.

빈곤율에 따른 서브그룹 분석에서는, 저소득 ZIP 코드(빈곤율 <15 %)에서 낮은 PM2.5 농도 구간에서도 위험곡선이 급격히 상승하는 반면, 고소득 지역에서는 위험이 비교적 완만하게 증가하였다. 이는 사회경제적 취약성이 PM2.5와 산불 연관 위험을 복합적으로 악화시킬 수 있음을 시사한다.

지역별 분석에서는 동북부, 중서부, 남부, 서부 네 개의 미국 인구통계학적 구역마다 노출‑반응 형태가 다르게 나타났으며, 특히 서부와 남부에서 산불 일수가 많은 지역에서 위험이 가장 크게 증폭되었다. 이는 기후·지형·산불 빈도 차이가 지역별 위험에 영향을 미친다는 점을 뒷받침한다.

연구의 강점은 1억 명 규모의 인구 데이터를 ZIP 코드 수준으로 정밀 매핑하고, 최신 고해상도 대기 모델과 산불 PM2.5 추정치를 결합한 점이다. 또한, 장기 추적 설계와 다중 혼란 변수 조정으로 인과 추론의 신뢰성을 높였다. 한계로는 관측 연구 특성상 잔여 교란(예: 개인 행동, 실내 공기질)과 노출 측정 오류가 존재할 수 있다. 또한, 산불 PM2.5를 0 µg/m³ 초과 일수만으로 정의함으로써 농도 강도 차이를 충분히 반영하지 못했을 가능성이 있다.

정책적 함의는 명확하다. 현재 NAAQS 기준 이하라 하더라도 산불 빈도가 높은 지역에서는 추가적인 위험 관리가 필요하다. 특히 사회경제적 취약계층을 대상으로 한 맞춤형 공공보건 전략과, 산불 예방·진압 정책이 장기 대기오염 관리와 연계되어야 한다.