Performance comparison of the Penman-Monteith-based evapotranspiration models in a citrus orchard of Southwest China

기술 요약: 중국 남서부 감귤 과수원의 Penman-Monteith 기반 증발산 모델 성능 비교

문제 제ない (Problem Statement)
작물 증발산(ET)을 정확하게 추정하는 것은 관개 전략을 최적화하고 농업의 지속 가능한 물 사용을 보장하는 데 필수적이다. Penman-Monteith(P-M) 모델은 식생이 잘 발달한 생태계에서는 산업 표준이지만, 감귤 과수원과 같은 불균질한 경관에서는 고전적인 단일 소스 "빅 리프(big-leaf)" 프레임워크로 인해 어려움을 겪는다. 토양 증발(E)과 식물 증산(T)을 분리하지 못하고, 복잡한 기공 반응을 벌크 저항 항으로 단순화하는 모델의 한계는 특히 식생이 희박하거나 패치 형태를 띠는 캐노피에서 부정확성을 초래할 수 있다. 이러한 한계를 해결하기 위해 이중 소스 및 하이브리드 모델이 제안되었으나, 다양한 시간 해상도(30분 단위에서 일 단위까지)와 페놀로지(생육 단계)에 걸쳐 모델 성능을 평가하고 구조적 동인에 대한 정량적 민감도 분석을 수행한 아열대 과수 생태계에서의 엄격한 다규모 비교 연구는 부족한 실정이다.

연구 방법론 (Methodology)
본 연구는 2020년부터 2023년까지 중국 쓰촨성 칭라이시의 청귤(Citrus Tachibana) 과수원에서 수집된 4년간의 에디 공분산(EC) 데이터를 활용하였다. 연구에서는 다음과 같이 서로 다른 개념적 구조를 가진 5가지 P-M 기반 ET 모델을 평가하였다:

1. Penman-Monteith (P-M): 단일 소스 모델.
2. Shuttleworth-Wallace (S-W): 이중 소스 층상 모델.
3. Two-Patch (T-P): 독립적인 토양 및 식생 패치를 가정하는 이중 소스 패치 모델.
4. Topography and Vegetation-based surface energy partitioning algorithm (TVET): S-W 저항 네트워크와 T-P 집합 방식을 결합한 하이브리드 모델.
5. Hybrid-Dual Source (H-D): 에너지 분배를 위한 층상 체계와 플럭스 집합을 위한 패치 방법을 채택하고 대기 안정도 보정을 포함하는 모델.

모델 성능은 통계적 지표인 $R^2$, RMSE, MAE 및 일치도 지수($d_1$)를 사용하여 30분 단위, 일 단위 및 생육 단계별(개화/결실기, 유과기, 과실 비대기, 과실 성숙기) 규모에서 평가되었다. 또한, 전역 민감도 분석을 수행하여 주요 저항 매개변수(캐노피/엽 기공 저항, 공기역학적 저항, 토양 저항)와 환경 입력값(복사 에너지, 온도, VPD, 풍속, 토양 수분 함량, 엽면적 지수)이 시뮬레이션된 ET에 미치는 영향을 정량화하였다.

주요 결과 (Key Results)

• 시간적 규모 성능: 모델 성능은 시간 해상도에 따라 크게 달라졌다. 30분 단위 규모에서 H-D 모델은 가장 높은 정확도($R^2 = 0.84$, RMSE = 0.038 mm·0.5h⁻¹)를 달성하며 일변화 주기를 효과적으로 재현하였다. S-W 모델은 정오 시간대에 체계적인 과대 추정 경향을 보였다. 일 단위 규모에서는 패치 기반 모델(T-P, TVET, H-D)이 단일 소스(P-M) 및 층상(S-W) 모델보다 일관되게 우수한 성능을 보였다. H-D 모델은 역시 탁월한 정확도($R^2 = 0.91$, RMSE = 0.47 mm·d⁻¹)와 최소한의 평균 과소 추정(1.63%)을 나타냈다.
• 페놀로지 민감도: 모델의 효용성은 생육 단계에 따라 달랐다. H-D 모델은 개화/결실기 및 과실 성숙기에 뛰어난 성능을 보였다. 반면, 더 단순한 T-P 모델은 과실 비대기에 가장 좋은 성능($R^2 = 0.88$)을 보였는데, 이는 토양 증발이 억제되는 밀집된 캐노피 조건에서는 복잡한 저항 네트워크가 수확 체감이나 추가적인 오차를 유발할 수 있음을 시사한다.
• 민감도 분석: 시뮬레이션된 ET는 P-M 모델에서는 캐노피 저항($r_c$)에, 다른 모델들에서는 엽 기공 저항($r_{sc}$)에 가장 민감하게 반응하였다. 이러한 매개변수를 20% 감소시키면 ET가 3.71%~7.41% 증가하였다. 환경 입력값 중에서는 토양 수분 함량(SWC)이 가장 반응성이 높은 요인이었으며, SWC가 20% 증가할 때 모델 간 ET 변동폭은 10.43%에서 32.27% 사이였다. S-W 모델은 엽면적 지수(LAI)에 가장 높은 민감도를 보였는데, LAI가 20% 감소하면 ET가 11.38% 감소하였다. 반면 패치 기반 모델들은 증산과 증발 사이의 상쇄 관계 덕분에 상대적으로 낮은 민감도를 보였다.

핵심 기여 (Key Contributions)
본 연구는 단순한 성능 순위 산정을 넘어 모델 거동의 구조적 원인을 밝힘으로써, 아열대 과수원에서 P-M 모델군에 대한 체계적이고 다규모적인 평가를 제공한다. 본 연구는 대기 안정도 보정과 이중 소스 에너지 분할을 통합함으로써, 불균질한 과수원 환경에 가장 견고한 옵션으로 H-D 모델을 식별하였다. 또한, 모델 선택은 맥락 의존적이어야 함을 명확히 하였다. 즉, H-D 모델이 일반적으로 우수하지만, 밀집된 캐노피 기간에는 T-P와 같은 더 단순한 패치 기반 모델이 충분히 효율적일 수 있다. 민감도 분석은 특정 저항 매개변수와 환경 강제력에 대한 이 모델들의 의존성을 해체하여, SWC와 기공 저항이 주요 불확실성 동인임을 강조하였다.

의의 (Significance)
본 연구의 결과는 아열대 과수 생태계에서의 모델 선택 및 개선을 위한 기계론적 근거를 제공한다. 특정 시간적, 페놀로지적 맥락에 따른 최적의 모델을 식별함으로써, 본 연구는 감귤 생산을 위한 정밀한 물 관리 및 관개 전략을 지원한다. 또한, ET 추정 오차를 줄이기 위해 식물 생리적 반응과 토양 수리학적 과정을 정확하게 매개변수화하는 것이 필수적임을 강조하며, 궁극적으로 물에 민감한 농업 시스템의 물 이용 효율과 지역 수자원 관리에 기여한다.