Heat, wind, and the calendar: why a warming climate makes its performance effects unidentifiable at the IRONMAN World Championship

기술 요약: IRONMAN 퍼포먼스 분석에서의 열, 바람, 그리고 달력

문제 제기
하와이 코나에서 열리는 IRONMAN 세계 선수권 대회는 환경 조건이 초장거리 지구력 퍼포먼스에 미치는 영향을 연구하기 위한 전형적인 데이터셋 역할을 한다. 지구 온도가 상승함에 따라 열과 바람이 경기 결과에 어떻게 영향을 미치는지 정량화하려는 관심이 높아지고 있다. 그러나 이 데이터셋에 대한 기존의 분석들은 순환성(circularity), 혼란 변수(confounding), 그리고 분산의 오귀속(misattribution of variance)을 유발할 수 있는 방법론적 선택에 의존해 왔다. 본 연구가 다루는 핵심 문제는 코나 레이스 아카이브에서 보고된 외견상의 환경 효과가 진정한 생리학적 신호인지, 아니면 달력 시간, 기후 트렌드, 그리고 퍼포먼스 향상 사이의 공선성(collinearity)에서 발생하는 통계적 인공물인지 여부이다.

방법론
본 연구는 2002년부터 2025년 사이(2020년 및 2021년 제외)의 IRONMAN 세계 선수권 대회 완주자 42,605명을 포함하는 현재까지 가장 포괄적인 데이터셋을 활용한다. 분석은 식별 가능성(identifiability)과 혼란 변수에 관한 네 가지 특정 연구 질문을 해결하기 위해 엄격한 통계적 프레임워크를 채택한다:

1. 데이터 구축 및 환경 재구성:

   • 코호트: 표준 연령 그룹(18~69세)의 일반인 완주자 및 프로 선수.
   • 환경 지표: 저자는 ERA5 재분석 데이터를 사용하여 180km 자전거 코스와 42.2km 러닝 코스의 지리적으로 위치한 웨이포인트를 따라 샘플링함으로써, 경로에 분포된 태양광 검증 습구흑구온도(WBGT)와 풍하중(wind loads)을 재구성하였다.
   • 지표 수정: 본 연구는 기존 문헌에서 사용된 "시간 가중 누적 열 부하(°C·h)"가 기계적으로 경기 지속 시간(duration)에 비례한다는 점을 지적하며 이를 명시적으로 거부한다. 대신 이를 지속 시간에 독립적인 강도 측정치(°C)와 외생적인 고정 창 참조값으로 대체한다.
   • 통계 모델링: 퍼포먼스는 얇은 판 스무스(thin-plate smooths)를 사용하는 일반화 가법 모델(GAMs)을 사용하여 모델링된다. 여기에는 세속적 트렌드(secular trends)를 흡수하기 위한 성별/연도별 스무스와 성별, 연령, 기술적 시대(특히 2017년 이후의 카본 플레이트 러닝화 및 2016년 이후의 사이클링 "최적화" 시대)에 대한 파라미터 제어가 포함된다.
   • 검증: 결과는 Leave-years-out 교차 검증과 누적 국소 효과(Accumulated Local Effects, ALE) 플롯을 사용하는 극단적 그래디언트 부스팅(XGBoost)을 통해 삼각 측량되어, 외삽 없이 비선형 관계를 특징짓는다.

2. 분석 전략:

   • 본 연구는 러닝 시간의 분산을 인구통계학적 특성, 세속적 트렌드, 사이클링 퍼포먼스 전이(carryover), 그리고 환경적 요인으로 구분한다.
   • 단일 기상 관측소가 분포된 샘플링에 비해 불충분하다는 가설을 테스트한다.
   • 달력 연도와 경기 당일 기후 변수 사이의 공선성을 조사하여 열 효과의 "식별 가능성"을 평가한다.

주요 기여
본 논문은 코나 데이터셋에 대한 방법론적 비판과 교정된 경험적 분석을 기여한다:

• 순환적 지표의 식별: 누적 열 부하 지표(°C·h)가 예측 변수에 결과 변수(지속 시간)를 내포함으로써 겉보기 상관관계를 부풀리는 인공물임을 입증한다.
• 달력에 의한 혼란: 단일 장소, 연간 개최 방식의 레이스 설계에서는 경기 당일의 열이 달력 연도의 세속적 트렌드와 통계적으로 구별 불가능함을 밝힌다. 기후가 연구 기간 동안 온난화되었기 때문에, "열"은 "시간"과 완벽하게 결합(aliased)되어 있으며, 피험자 내 변동(within-subject variation) 없이는 인과적인 열 효과를 격리하는 것이 불가능하다.
• 자전거-러닝 결합의 분해: 선수의 능력(퍼포먼스 전이)과 환경적 전이를 분리하여 자전거-러닝 관계를 정량화한다.
• 공간적 샘플링의 필요성: 분포된 샘플링이 열(공간적으로 일관됨)을 포착하는 데는 중복적이지만, 코스 전반에 걸쳐 크게 변하는 풍하중을 포착하는 데는 필수적임을 명확히 한다.

결과

1. 열 부하의 순환성: 시간 가중 누적 열 부하는 $r = 0.93$의 상관관계를 보였다. 이를 지속 시간 독립적 강도로 대체했을 때 상관관계는 $r = -0.24$로 떨어졌으며, 퍼포먼스와의 겉보기 연관성은 현저히 감소했다(부분 $r \approx 0.25$에서 $0.33$).
2. 식별 불가능한 열 효과: 경기 당일 WBGT는 10년당 $+1.28^\circ\text{C}$씩 상승했다($p < 0.001$). 결과적으로, 일일 참조 WBGT 분산의 99.5%는 달력 연도만으로 설명된다. 연도 트렌드를 통제하면 남은 열 변동이 없어 열 효과를 추정할 수 없다. 환경적 추론을 위한 유효 샘플 크기는 42,605명의 선수가 아니라 약 23번의 경기 당일 기후이다.
3. 자전거-러닝 결합: 사이클링 조건과 러닝 퍼포먼스 사이의 겉보기 결합은 압도적으로 선수의 능력에 의해 주도된다.
   • 러닝 모델에 사이클링 퍼포먼스(분할 시간 및 완주 시계)를 추가했을 때 설명된 분산($\Delta R^2$)은 0.379만큼 증가했다(이득의 99%).
   • 사이클링 환경 조건과 열-피로 상호작용을 추가했을 때 $\Delta R^2$는 0.003(1%)만 증가했다.
   • 환경적 결합에 대한 95% 신뢰 구간은 $[-0.006, +0.003]$으로, 0과 구별할 수 없다.
4. 식별 가능한 바람 효과: 열 부하와 달리, 사이클링 구간의 풍하중은 진정한 설명력을 유지했다. 바람 지수는 잔잔한 상태에서 바람이 부는 상태까지 약 50분의 부분 효과 범위를 설명했다. 이 효과는 기후가 단조로운 온난화 트렌드를 보이는 것과 달리, 연도 내에서 코스와 출발 웨이브에 따라 변하기 때문에 식별 가능했다.
5. 샘플링 해상도: 단일 기상 관측소는 열을 포착하는 데 있어 분포된 샘플링만큼 효과적이었으나($r = 0.81$에서 $0.87$), 바람의 경우 단일 관측소(예: 부두)는 코스 평균 바람을 포착하는 데 실패했다(상관관계$r = 0.20$). 즉, 식별된 바람 효과를 주도하는 특정 횡풍 체계를 포착하기 위해서는 분포된 샘플링이 필요했다.

의의 및 주장
본 논문은 IRONMAN 세계 선수 championship에서 기존에 열 효과로 귀속되었던 주요 환경적 퍼포먼스 효과들이 대부분 분석적 선택과 연구 설계의 한계로 인한 인공물이라고 주장한다.

• 방법론적 경고: 본 연구는 환경 퍼포먼스 연구를 위한 주의 깊은 템플릿 역할을 하며, 기후가 시간과 함께 변화하는 단일 장소 연간 설계에서는 샘플 크기에 관계없이 열 효과를 식별할 수 없음을 보여준다.
• 내러티브의 재구성: "자전거-러닝" 결합은 환경 조건이 아닌 일반적인 선수 능력의 함수라고 단언한다. 이 데이터셋에서 확인된 유일하고 견고하며 식별 가능한 환경 결정 요인은 사이클링 구간의 바람 부하이다.
• 생리학적 vs 통계적: 저자들은 열 효과에 대한 자신들의 귀무 결과(null result)가 통제된 연구에서 잘 입증된 열 손상의 생리학적 기전과 모순되는 것이 아니라, 특정 집계된 단일 장소 설계가 공선성 때문에 이러한 효과를 해결하기에 불충분함을 강조한다.
• 실무적 시사점: 실무자들에게 본 연구는 열 관리가 안전을 위해 여전히 중요하지만, 집계된 레이스 기록을 사용하여 특정 열에 따른 시간 페널티를 예측할 수는 없음을 시사한다. 반대로, 바람 부하는 사이클링 구간에서 정량화 가능하고 실행 가능한 환경적 요인이다.

본 논문은 열 효과에 대한 향와 연구가 기후와 달력 시간 사이의 공선성을 깨뜨리기 위해 피험자 내 반복 측정 설계 또는 다중 장소 비교를 필요로 한다고 결론짓는다.