Educational inequalities in heat-related mortality across Belgian districts

기술 요약: 벨기에 지역별 열 관련 사망률의 교육적 불평등

문제 제기
기후 변화는 극한 폭염의 빈도와 강도를 심화시키며, 취약 계층에게 불균형적인 건강 위험을 초래하고 있다. 사회경제적 지위가 노출 및 적응 능력에 영향을 미친다는 사실은 잘 알려져 있으나, 특히 개인적 차원에서 교육 수준이 열 관련 사망률을 어떻게 조절하는지에 대한 증거는 제한적이다. 기존 연구들은 주로 지역 단위의 대리 지표에 의존해 왔으며, 이는 개인적 상황과 지리적 맥락이 각각 어떻게 취약성을 유발하는지에 대한 이해의 공백을 남겼다. 본 연구는 사회경제적 및 공간적 건강 불평등이 뚜렷하게 나타나는 지역인 벨기에의 42개 NUTS-3 구역을 대상으로, 고령자(60세 이상)의 열 취약성에서 나타나는 교육 수준의 경사도를 정량화하고자 한다.

연구 방법론
본 연구는 2001년 인구 조사에서 도출된 개인별 교육 수준과 연결된 종합적인 벨기에 사망 기록 데이터셋(2002~2019년)을 활용한다. 분석은 일관된 노출 분류를 보장하기 위해 폐쇄 코호트(closed cohort)로 제한된 60세 이상 성인을 대상으로 한다.

• 데이터 통합:

  • 사망: 교육 수준이 확인된 1,392,134건의 사망 사례를 분석하였으며, 이를 세 그룹(저학력: 무학/초등/중등 하위, 중등 교육: 상급 중등/비고등 후기 교육, 고등 교육: 학사 학위 이상)으로 분류하였다.
  • 온도: 벨기에 왕립 기상 연구소의 일일 최고 기온(1992~2022년) 데이터를 NUTS-3 구역 경계에 공간적으로 결합하였다.
  • 맥락 변수: 벨기에 다차원 결핍 지수(Belgian Index of Multiple Deprivation)에서 추출한 구역 단위 결핍 점수를 메타 예측 변수로 구성하였다.

• 분석 프레임워크:
  저자들은 2단계 분포 지연 비선형 모델링(two-stage distributed lag non-linear modelling, DLNM) 프레임워크를 적용하였다:

  1. 1단계: 126개 구역-교육 조합(42개 구역 × 3개 교육 수준) 각각에 대해 최대 7일의 지연 시간을 갖는 DLNM을 사용하여 비선형 및 지연 온도-사망 연관성을 추정하였다.
  2. 2단계: 다변량 메타 회귀 분석을 통해 이 추정치들을 통합하여 교육 그룹별 인구 평균 노출-반응 함수를 산출하였다. 이 프레임워크는 최적 선형 불편 예측값(BLUPs)을 활용하여 지역별 세부 수치를 조정하였으며, 1단계의 불확실성을 명시적으로 전파하였다.

• 정의 및 지표:

  • 극한 폭염: 국가 상위 97.5 백분위수(28.3°C)를 초과하는 날로 정의하였다.
  • 최소 사망 온도(MMT): 교육 그룹 간의 비교 가능성을 확보하기 위해, 그룹별 MMT(15.8°C에서 21.4°C 사이)의 차이에도 불구하고 공통 MMT인 21.0°C(국가 상위 77 백분위수)를 공유 중심점으로 설정하였다.
  • 결과: 2011년 벨기에 인구 조사 연령 분포에 따른 직접 표준화 방식을 사용하여 연령 표준화 열 관련 사망률(ASMR, 인구 10만 명당 연간 사망자 수)을 계산하였다. 불확실성은 몬테카를로 시뮬레이션(1,000회 반복)을 통해 정량화하였다.

주요 결과
분석 결과, 모든 벨기에 구역에서 극한 폭염에 대한 취약성의 명확하고 일관된 교육적 경사도가 발견되었다:

• 상대 위험도(RR): 97.5 백분위수(28.3°C)에서 열 관련 사망의 RR은 저학력 그룹이 1.12(95% CI: 1.10–1.15), 중등 교육 그룹이 1.11(1.08–1.14), 고등 교육 그룹이 1.10(1.05–1.16)이었다. 이 경사도는 더 극단적인 온도(99 백분위수, 31.3°C)에서 더욱 크게 벌어지는데, 이때 저학력 그룹의 RR은 1.28(1.24–1.31)까지 상승한 반면, 고등 교육 그룹은 1.16(1.09–1.22)에 머물렀다.
• 사망률: 극한 폭염 기간 동안의 연령 표준화 열 관련 사망률은 저학력 그룹이 26.5, 중등 교육 그룹이 19.4, 고등 교육 그룹이 13.0(인구 10만 명당 연간 사망자 수)으로 나타났다. 저학력 그룹의 신뢰 구간은 나머지 두 그룹과 겹치지 않으며, 이는 통계적으로 유의미한 격차가 존재함을 나타낸다.
• 공간적 일관성: 구역 간 이질성은 낮게 나타났으며($I^2 = 8.4\%$), 이는 교육적 경사도가 특정 예외 구역에 의해 주도되는 것이 아니라 국가적 현상임을 시사한다. 그러나 절대적 불평등은 서부(Tournai-Mouscron)와 남부(Virton, Arlon)에서 가장 두드러졌다.
• 기제: 스플라인 분해(Spline decomposition) 결과, 교육의 보호 효과는 가장 높은 온도에서만 통계적으로 유의미해졌으며(계수 $b_5$), 이는 고학력이 특히 극한 폭염 발생 시 위험을 완화함을 의미한다.

의의 및 주장
본 논문은 교육 수준이 벨기에 고령자 사이에서 극한 폭염 취약성을 조절하는 중요한 요인이며, 이는 지리적 맥락보다는 주로 개인적 차원의 경로(건강 문해력, 체온 조절 행동, 자원 접근성 등)를 통해 작용한다고 주장한다.

• 인구 통계적 특성과의 차별성: 연령 표준화 후에도 경사도가 지속된다는 점은 이러한 격차가 단순히 저학력 집단 내의 고령 구조에 의한 인구 통계적 결과물이 아니라, 실제 열 민감도의 차이를 반영한다는 것을 입증한다.
• 정책적 시사점: 기후 변화로 인해 현재의 극한 임계치를 초과하는 날이 증가함에 따라, 문서화된 불평등한 부담은 더욱 심화될 가능성이 높다. 저자들은 이것이 지리적 또는 지역 단위의 개입에만 의존하기보다 개인적 차원의 취약성을 다루는 표적화된 열 정책의 필요성을 강조한다고 주장한다.
• 한계점: 저자들은 본 연구가 교육의 효과를 쇠약성(frailty)이나 기저 질환 상태와 완전히 분리할 수 없음을 솔직하게 언급하였다. 이는 저학력층이 더 고령이며 기초 사망률이 높기 때문이다. 또한, 본 분석은 열의 짧은 지연 구조를 고려하여 열 관련 사망에 국한되었으므로, 한랭 관련 사망의 경사도를 신뢰성 있게 추정하는 데는 한계가 있다.

결론적으로, 본 연구는 온난화되는 기후 환경 속에서 낮은 교육 수준이 열 관련 사망 위험을 높이는 결정적인 증폭제로 작作하며, 벨기에의 고령 저학력 인구에게 누적된 불이익을 초래한다는 강력한 증거를 제공한다.