남온타리오 지역 강우 강도‑기간‑빈도 변화 평가

본 논문은 남온타리오 지역의 강우 강도‑기간‑빈도(IDF) 곡선 변화를 정밀하게 평가하기 위해 NA‑CORDEX 프로젝트에 포함된 3개의 지역기후모델(RCM) 앙상블을 활용하였다. 연구는 크게 네 단계로 진행되었다. 첫 번째 단계는 관측 데이터와 모델 출력의 품질 평가이다. 캐나다 남온타리오 지역에 위치한 8개 강우 관측소에서 1970‑2010년 기간의 시간별 강우 데이터를 수집하고, 동일 기간에 대한 RCM 시뮬레이션 결과와 비교하였다. 초기 비교에서 모델은 강우 강도와 빈도에서 체계적인 편향을 보였으며, 이는 고해상도 지역 기후 모델이 대규모 대기 순환을 정확히 재현하더라도 강우 강도의 미세구조를 과소/과대 평가하는 전형적인 문제였다. 두 번째 단계는 편향보정이다. 저자들은 분위수 기반 파라메트릭 보정(Quantile Mapping)과 비파라메트릭 보정(Distribution‑Free) 두 가지 방법을 적용하였다. 각 보정 방법에 대해 평균제곱오차(MSE), 피어슨 상관계수(R), 콜모고로프‑스미르노프(K‑S) 검정 등을 수행해 성능을 평가한 결과, 분위수 매핑이 대부분의 관측소와 지속시간에서 가장 낮은 MSE와 높은 R 값을 기록하였다. 따라서 최종 분석에서는 분위수 매핑을 기본 보정 방법으로 채택하였다. 세 번째 단계는 시간 해상도 향상 및 극값 추출이다. RCM은 일일 강우 데이터를 제공하므로, 1‑24 시간 지속시간별 강우를 재현하기 위해 멀티플리케이티브 캐스케이드(Multiplicative Cascade) 기법을 적용하였다. 이 기법은 대규모 강우 패턴을 작은 시간 스케일로 분해하면서 공간·시간 일관성을 유지한다. 보정·다운스케일링된 시계열에서 연간 최대 강우량(Annual Maxima, AM)을 추출하고, 이를 일반화 극값분포(GEV)로 모델링하였다. GEV 파라미터 추정에는 두 가지 접근법을 사용하였다. 첫 번째는 전통적인 정상(stationary) 가정으로, 위치(μ), 규모(σ), 형태(ξ) 파라미터를 시간에 대해 고정하였다. 두 번째는 비정상(non‑stationary) 가정으로, 위치 파라미터 μ를 선형 시간 트렌드(μ(t)=μ0+β·t)와 연결시켜 미래 기후 변화가 극값 평균에 미치는 영향을 직접 반영하였다. 두 모델 모두 최대우도법(Maximum Likelihood Estimation)으로 파라미터를 추정했으며, AIC와 BIC를 통해 모델 적합성을 비교하였다. 네 번째 단계는 IDF 곡선 생성 및 비교이다. GEV 파라미터와 반환주기(T) 관계식 I(T)=μ+σ·