Higher road-traffic mortality on warm days in Thailand: a nationwide time-stratified case-crossover study, 2011–2025

기술 요약: 태국의 온난한 날씨에서의 높은 도로 교통 사망률

문제 제기 및 동기
도로 교통 외상은 전 세계적으로, 특히 WHO 서태평양 지역에서 주요 사망 원인이지만, 주변 온도와 도로 교통 사망 사이의 구체적인 연관성은 비외인성(non-accidental) 원인에 비해 그동안 충분히 탐구되지 않았다. 고온 효과에 대한 메타 분석을 포함한 기존 문헌들은 주로 비외인성 이환율 및 사망률에 집중해 왔거나, 국지적 온도 분포상 어디에서 위험이 축적되는지를 특성화하지 않은 채 총체적인 글로벌 부담 추정치만을 제공해 왔다. 또한, 동남아시아의 증거는 부족한 실정이다. 이전의 지역 연구들은 전인적 사망(all-cause), 심혈관 또는 호흡기 사망에 초점을 맞추었으며, 태국 및 더 넓은 ASEAN 지역의 열-건강 관련 문헌에서 도로 교통 사망은 눈에 띄게 부재했다.

현재의 열-건강 경보 시스템(대상이 되는 극심한 폭염)이 타겟으로 하는 극심한 더위의 날에 부담이 집중되는지, 아니면 완만하게 따뜻한 날들에 걸쳐 축적되는지에 대한 이해에는 결정적인 공백이 존재한다. 또한, 연구마다 참조 온도의 가변성이 보고된 상대 위험도(relative risks)의 비교를 저해한다. 본 연구는 태국(높은 도로 사망률과 오토바이가 주를 이루는 차량 보유 대수를 가진 국가)을 대상으로 온도-사망 연관성의 형태와 국지적 온도 범위에 따른 초과 사망의 분포를 조사함으로써 이러한 공백을 해결하고자 한다.

연구 방법론
본 연구는 2011년 1월 1일부터 2025년 8월 24일까지의 데이터를 사용하여 전국 단위의 시계열 층화 케이스-크로스오버(time-stratified case-crossover) 설계를 채택하였다.

• 데이터 소스:
  • 결과(Outcome): 보건부 감시 체계, 태국 왕립 경찰 기록, e-claim 데이터를 통합한 RTDDI(Road Traffic Death Data Integration) 파일을 통해 도출된 태국 77개 주의 283,847건의 도로 교통 사망 사례.
  • 노출(Exposure): 1,347개의 조화된 기상 관측소와 일일 최고 기온($T_{max}$)의 0.1° 격자 재구성 데이터를 포함하는 커스텀 노출 데이터셋(총 450만 개의 일일 관측치).
• 통계 모델:
  • 시계열 층화 케이스-크로스오버 설계에 조건부 포아송(conditional Poisson) 모델을 적용하였다.
  • 층화(Stratification): 장기적 추세, 계절성 및 요일 효과를 통제하기 위해 '주(province) × 연(year) × 월(month) × 요일(day-of-week)'로 층을 정의하였다.
  • 노출 지표: 국지적 기후 변동(예: 우본 라차타니에서는 36°C가 일상적이지만 푸켓에서는 극단적인 경우)을 고려하기 위해 주별 백분위 밴드(<P5, P5–25, P25–40, P40–60 [참조], P60–75, P75–90, P90–95, P95–99, ≥P99)를 사용하였다.
  • 교란 변수(Confounders): 사망 피크가 발생하는 것으로 알려진 강수량과 특정 공휴일(송크란 및 신년)을 조정하였다.
  • 추정(Estimation): 밀집된 층 고정 효과(dense stratum fixed effects)를 가진 표준 포아송 GLM과 검증을 마친 커스텀 조건부 포아송 추정치를 구현하였으며, 두 모델 간의 일치도는 $5.8 \times 10^{-15}$였다. 분산은 주 단위의 클러스터 강건(cluster-robust) 방법을 사용하여 추정하였다.
• 민감도 분석:
  • 공간 해상도: 사고 좌표를 사용하여 0.1° 격자 노출을 할당하여 107,313건의 지리적 위치가 파악된 사망 사례를 재분석하였다.
  • 시간적 할당: 사망 날짜 대신 사고 날짜를 기준으로 노출을 할당하여 94,437건의 사망 사례를 재분석하였다.
  • 대안 지표: 복합 습열 지수(열지수, 체감온도, 습구온도) 및 15가지 다른 폭염 정의를 테스트하였다.
  • 공휴일 제외: 송크란과 신년 기간을 완전히 제거하여 교란 효과를 테스트하였다.

주요 결과

1. 노출-반응 형태: 도로 교통 사망 위험은 분포의 낮은 온도 쪽에서 시작하여 약 75번째 국지 백분위수까지 상승한다. 이 지점 이후 위험은 정체(plateau)된다.
   • 참조군(P40–60) 대비 상대 위험(RR) P75–90: 1.045 (95% CI 1.023–1.068).
   • RR P95–99: 1.050 (1.018–1.083).
   • RR ≥P99: 1.041 (0.979–1.108).
   • 이진 대조(예: ≥P95 vs. <P95)는 이미 상승된 위험을 포함하는 그룹(P60–90)과의 비교이기 때문에 유의미하지 않거나 약한 결과(RR 1.015)를 보였다.
2. 공간적 및 시간적 강건성: 관측소 데이터와 격자 데이터 간의 평균 절대 차이가 2.95°C임에도 불구하고, 사고 좌표를 통한 0.1° 미세 공간 해상도 할당 시에도 정체되는 형태가 유지되었다. 마찬가지로, 사망 날짜가 아닌 사고 날짜를 기준으로 노출을 할당해도 연관성의 형태는 변하지 않았다.
3. 습도 및 폭염: 복합 습열 지수는 건구 온도(dry-bulb temperature)보다 모델 적합도를 개선하지 못했다. 다양한 임계값과 지속 시간을 가진 폭염 정의들은 일일 온도 외에 유의미한 총 효과나 지속성에 따른 추가 효과를 보여주지 않았다.
4. 기여 부담:
   • P40–60 참조 밴드 이상의 온도와 관련된 초과 사망은 연구 기간 동안 4,509명(전체 도로 교통 사망의 1.59%)이었다.
   • 부담의 분포: 이 초과 사망의 약 72%는 60번째와 90번째 백분위수 사이(완만하게 따뜻한 날들)에서 발생했다. 극심한 더위(≥P95)에서 발생한 초과 사망은 14.8%에 불 old했다.
5. 검증: 모델은 송크란(RR 1.311) 및 신년(RR 1.194) 기간의 알려진 사망률 상승을 성공적으로 회복하였다. 송크란 기간을 제외하더라도 온도-반응 곡선은 변하지 않았다.

의의 및 주장
본 논문은 태도의 도로 교통 사망이 중간 온도인 날보다 따뜻한 날에 유의미하게 높으며, 완만하게 따뜻한 온도에서 증가하다가 75번째 백분위수 위에서 평탄해진다고 주장한다.

주된 의의는 부담의 분포에 있다. 극심한 더위의 날(P90 이상)보다 완만하게 따뜻한 날(P60–90)이 훨씬 더 빈번하기 때문에, 이들이 온도 관련 초과 사망의 대다수(72%)를 차지한다. 결과적으로, 본 논문은 극심한 폭염만을 감지하도록 설계된 현재의 열-건강 경보 시스템이 문제와 일치하지 않을 수 있다고 주장한다. 이러한 결과는 전체 온도 관련 부담 중 아주 작은 부분만을 차지하는 극단적인 날들에 국한되기보다는, 따뜻한 계절 전체에 걸친 지속적인 도로 안전 개입이 필요함을 시사한다.

또한 본 연구는 이전의 글로벌 추정치에서 동남아시아 데이터의 부재가 중요한 공백이었음을 강조하며, 특정 연관성의 형태(P75 이후의 정체)가 이 결과에 대한 이진적 "극심한 폭염" 대조의 효용성에 도전한다는 점을 밝힌다.