기온 측정 모델의 선택이 폭염 사망 위험 예측에 미치는 영향 분석

초록

이 연구는 이탈리아 라치오 지역을 대상으로 세 가지 서로 다른 기온 재구성 모델(베이지안 양자 회귀, 베이블리안 가우시안 회귀, ERA5-Land)이 폭염 관련 사망 위험 추정에 미치는 영향을 비교 분석했습니다. 연구 결과, 모델 선택에 따라 최소 위험 온도와 상대 위험도가 크게 달라지며, 특히 베이지안 모델이 고온 극단값을 더 정밀하게 포착하여 보건 정책 수립 및 조기 경보 시스템 구축에 있어 정밀한 온도 모델링의 중요성을 강조합니다.

상세 분석

본 논문은 기후 변화에 따른 보건 위기 대응을 위해 온도 모델링의 정밀도가 역학적 추정치에 미치는 영향을 통계학적 관점에서 심도 있게 다루고 있습니다. 연구의 핵심 기술적 차별점은 1단계 온도 재구성 과정에서 적용된 ‘베이지안 공간 양자회귀(QR-GP)’ 모델에 있습니다. 이 모델은 비대칭 라플라스(Asymmetric Laplace) likelihood를 사용하여 온도 분포의 특정 분위수(0.05~0.95)를 추정함으로써, 단순 평균값이 놓치기 쉬운 기온의 극단적 변동성을 포착합니다. 특히 공간적 Gaussian Process(GP)를 통해 고도와 계절성 변수를 포함한 복잡한 공간적 자기상관 구조를 모델링함으로써, ERA5-Land와 같은 격자형 재분석 자료가 가진 지형적 평활화(smoothing) 한계를 극복했습니다.

2단계 역학 분석에서는 케이스-크로스오버(Case-crossover) 설계를 채택하여 시간불변 교란 요인을 통제하고, 베이지안 조건부 포아송 모델을 통해 온도와 사망률 간의 비선형적 관계를 스플라인(Spline) 함수로 정밀하게 추정했습니다. 분석 결과, 베이지안 모델 기반의 추정치는 ERA5-Land에 비해 최소 위험 온도($T_{min}$)를 약 4°C 가량 낮게 산출(24-26°C vs 28-30°C)했는데, 이는 기존의 재분석 자료를 기반으로 한 보건 정책이 실제 위험을 과소평가할 위험이 있음을 시사합니다. 또한, 열파(Heatwave)의 정의(임계값 및 지속 기간)에 따른 위험도 변동을 정량화함으로써, 단순 고온 노출뿐만 아니라 지속적인 열 스트레스가 미치는 영향을 통계적으로 입증했습니다. 이는 향후 기후 적응 전략 수립 시, 단순 기온 상승뿐만 아니라 국지적 기온 변동성을 반영한 정밀한 모델링이 필수적임을 기술적으로 뒷받침합니다.

본 연구는 이탈리아 라치오(Lazio) 지역의 2008년부터 2022년까지의 데이터를 바탕으로, 온도 재구성 모델의 차이가 폭염 관련 사망 위험 추정에 미치는 영향을 두 단계의 분석 프레임워크를 통해 실증적으로 분석하였습니다.

첫 번째 단계에서는 일일 최고기온을 시·군 단위로 재구성하기 위해 세 가지 모델을 비교하였습니다. 첫째, 베이지안 공간 양자회귀(QR-GP) 모델은 비대칭 라플라스 분포와 Gaussian Process를 결합하여 기온의 분위수별 분포와 공간적 상관성을 정밀하게 추정하였습니다. 둘째, 베이지안 가우시안 회귀 모델은 관측소 데이터를 기반으로 공간적 상관구조를 모델링하였습니다. 셋째, ERA5-Land 재분석 자료를 사용하였습니다. 분석 결과, 베이지안 모델들은 ERA5-Land에 비해 평균 최고기온이 0.8~1.5°C 높게 나타났으며, 특히 복잡한 지형을 반영하여 공간적 변동성이 훨씬 크게 나타났습니다. 이는 기존의 격자형 재분석 자료가 국지적인 고온 현상을 과소평가할 가능성이 있음을 보여줍니다.

두 번째 단계에서는 개인별 사망 기록(심혈관 및 호흡기 질환)을 활용하여 온도와 사망률 간의 연관성을 분석하였습니다. 케이스-크로스오버 설계를 통해 요일, 월, 연도 등 시간불변 교란 요인을 제거하였으며, 3일 평균 최고기온을 노출 변수로 설정하였습니다. 또한, 다양한 임계값과 지속 기간을 조합한 열파 정의를 적용하여 분석을 수행했습니다.

주요 연구 결과는 다음과 같습니다. 첫째, 모든 모델에서 고온 시 상대위험(RR)이 급격히 상승하는 경향은 일치했으나, 사망 위험이 가장 낮은 온도인 ‘최소 위험 온도($T_{min}$)‘는 모델 간에 큰 차이를 보였습니다. 베이지안 모델에서는 2426°C로 나타난 반면, ERA5-Land에서는 2830°C로 나타나 모델 선택이 위험 임계값 설정에 결정적인 영향을 미침을 확인했습니다. 둘째, 취약 계층 분석 결과, 여성과 80세 이상의 고령층에서 상대위험(RR)의 상승 폭이 1.3~1.6배로 나타나 성별 및 연령에 따른 뚜렷한 취약성을 확인했습니다. 셋째, 열파의 영향력은 정의에 따라 달랐으나, 90th percentile 이상의 기온이 3일 이상 지속될 때 RR이 1.45(95% CI 1.20-1.73)로 가장 높게 나타나, 지속적인 열 노출의 위험성을 입증했습니다.

결론적으로, 본 연구는 온도 노출 모델의 선택이 역학적 결과와 공중보건 권고 사항을 완전히 바꿀 수 있음을 보여줍니다. 특히 베이지안 모델이 고온의 극단값을 더 잘 포착한다는 점을 고려할 때, 향후 기후 변화에 대응하는 조기 경보 시스템과 보건 적응 전략은 국지적 기온 변동성을 정밀하게 반영할 수 있는 고도화된 모델링을 기반으로 설계되어야 합니다.