강우와 식생 이질성이 토양수분·증발산에 미치는 규모 의존적 효과

초록

본 연구는 식생 유형과 강우의 공간적 이질성이 토양수분과 증발산의 평균값 사이에 비선형·히스테리시스 관계를 만들며, 이러한 관계는 평균 영역의 크기와 이질성 규모에 따라 달라진다는 것을 실험적 시뮬레이션을 통해 입증한다. 기존의 평균 파라미터 방식보다 제시된 업스케일링 함수가 텍사스 건조 관목림에서 토양수분과 증발산을 더 정확히 예측한다는 결과를 제시한다.

상세 분석

이 논문은 식생-토양-강우 시스템을 식물 규모에서 지역 규모로 확대할 때 발생하는 비선형성 및 히스테리시스 현상을 정량적으로 분석한다. 먼저, 식물 수준에서 토양수분(θ)과 증발산(ET) 사이의 함수는 단일값 관계로 가정되지만, 실제 현장에서는 강우 강도와 지속시간, 식생 종류(관목, 초목 등)의 공간적 분포가 서로 다르기 때문에 평균화 과정에서 다중값(히스테리시스) 관계가 나타난다. 저자들은 2차원 격자 모델에 강우 이벤트를 무작위로 배치하고, 각 격자 셀에 서로 다른 식생 유형 파라미터(근권 깊이, 최대 전도성 등)를 할당함으로써 ‘명시적 이질성 모델’을 구축하였다.

시뮬레이션 결과, 평균 토양수분이 일정 범위 내에서 변동할 때 평균 증발산은 동일한 θ값에 대해 두 개 이상의 ET 값을 가질 수 있었으며, 이는 강우 강도와 식생 비율에 따라 순환 경로가 달라지기 때문이다. 특히, 평균 영역이 작을수록(예: 1 km²) 히스테리시스 루프가 크게 나타났고, 영역이 커질수록(예: 100 km²) 루프가 축소되었지만 완전히 사라지지는 않았다. 이는 공간적 이질성이 ‘스케일 의존적’이라는 핵심 결론을 뒷받침한다.

업스케일링 함수는 기존의 ‘effective‑parameter’ 접근법(즉, 토양·식생 파라미터를 평균값으로 대체)과 달리, 평균 θ와 평균 ET 사이에 히스테리시스 영역을 명시적으로 포함한다. 함수 형태는 θ에 대한 다항식 기반의 두 개의 분기(증가 단계와 감소 단계)로 구성되며, 각 분기의 계수는 강우 이질성 지표(표준편차)와 식생 비율에 따라 조정된다. 검증 결과, 텍사스의 건조 관목림 사례에서 평균 토양수분 RMSE가 0.032 m³ m⁻³, 평균 증발산 RMSE가 0.18 mm day⁻¹로, 기존 방법 대비 각각 약 30 %와 25 % 개선되었다.

또한, ‘leakage’(지표면을 통한 물 손실) 예측에 있어 단일값 함수는 큰 오차를 보였으며, 히스테리시스 구조를 포함한 함수가야만 평균 leakage를 정확히 재현할 수 있었다. 마지막으로, 통계‑동적 접근법과 Montaldo & Albertson 모델을 적용했을 때도 동일한 히스테리시스 현상이 나타났으며, 이는 이 현상이 모델 구조에 독립적인 보편적 특성임을 시사한다.