Reconciling intrinsic and ecosystem water-use efficiency reveals physiological mechanisms underlying carbon-water coupling across biological scales

기술 요약: 내재적 및 생태계 수분 이용 효율의 조화

문제 제기
탄소-수분 결합이 개별 잎에서 전체 생태계로 어떻게 규모 확장(scaling)되는지를 이해하는 것은 생태수문학에서 여전히 중요한 과제로 남아 있다. 수분 이용 효율(WUE)은 식물의 생리적 기능에 대한 핵심 지표이지만, 규모 불변적(scale-invariant)이지 않다. 잎 수준에서 내재적 수분 이용 효율($iWUE$)은 순 광합성 동화량($A$)과 기공 전도도($g_s$)의 비로 정의되며, 종단 탄소 동위원소($\delta^{13}C$)로부터 재구성된다. 반면, 생태계 규모의 $WUE$($WUE_{eco}$)는 일반적으로 에디 공분산(EC) 측정값으로부터 총 일차 생산량(GPP)과 총 증발산량(ET)의 비를 통해 도출된다. 이 두 지표는 장기 연구에서 종종 약하거나 일관되지 않은 관계를 보인다. 이러한 불일치는 $WUE_{eco}$가 $iWUE$ 재구성에 존재하지 않는 비증산적 수분 플럭스(토양 증발, 캐노피 차단, 하층 식생 증산)와 구조적 복잡성을 통합하기 때문에 발생한다. 또한, 장기적인 $iWUE$ 추세가 기공 폐쇄에 의한 것인지 아니면 광합성 능력의 향상에 의한 것인지를 구별하는 것은, 여러 데이터 소스를 통합하는 통일된 프레임워크의 부재로 인해 복잡해진 과제이며, 생리적 신호와 수문학적 제어를 분리해내는 작업이다.

연구 방법론
본 연구는 스위스 레게른(Lägeren, CH-Lae) 지역의 유럽 너도밤나무(Fagus sylvatica)가 우점하는 성숙한 혼합 낙엽수림에서 20년(2004~2023년) 동안 세 가지 상호 보완적인 접근 방식을 통합하였다:

1. 나이테 분석: 15그루의 F. sylvatica에서 추출한 셀룰로오스의 안정 동위원소($\delta^{13}C$, $\delta^{18}O$, $\delta^{2}H$)를 측정하였다. $\delta^{13}C$는 관측된 내재적 수분 이용 효율($iWUE_{obs}$)을 계산하는 데 사용되었으며, 삼중 동위원소 접근법은 $iWUE$ 변동성을 결정하는 기공 조절과 광합성 능력 사이의 원인을 구분하는 데 활용되었다.
2. 에디 공분산 (EC): 연속적인 반시각적 플럭스 측정을 통해 생태계 수준의 GPP와 ET 데이터를 제공받았으며, 이를 통해 관측된 생태계 수분 이용 효율($WUE_{eco}^{obs}$)을 계산하였다. 이때 일일 비율 평균화와 관련된 편향을 피하기 위해 합의 비율(ratio-of-sums) 방식을 사용하였다.
3. 프로세스 기반 모델링: 3D-CMCC-FEM 모델을 사용하여 탄소 및 수분 순환을 시뮬레이션하였다. 모델은 성숙한 F. sylvatica에 맞춰 특화 구성되었으며, 이를 통해 모델링된 내재적 수분 이용 효율($iWUE_{mod}$)과 모델링된 생태계 수분 이용 효율($WUE_{eco}^{mod}$)을 생성하였다.
4. 통계적 통합: 연구는 이 데이터셋들을 일간, 월간, 연간 시간 척도에 따라 비교하였다. 또한, $\delta^{18}O$와 $\delta^{2}H$의 변동성이 생리적 조절보다는 수원 역학에 의해 제어되는지 확인하기 위해 이중 동위원소 프레임워크(Scheidegger et al., 2000)와 삼중 동위원소 상관관계를 활용하여 $iWUE$ 추세를 생리적 메커니즘(기공 조절 vs 광합성 능력)과 수문학적 제어(원수 역학)에 귀속시켰다.

주요 결과

• 잎 규모의 일관성: 동위원소 유래 $iWUE_{obs}$와 모델 시뮬레이션 $iWUE_{mod}$ 사이에 강한 양의 상관관계($r = 0.61, p < 0.05$)가 나타났다. 두 데이터셋 모두 연구 기간 동안 약 20%의 장기적 증가를 나타냈다. 그러나 절대값은 유의미하게 달랐는데, $iWUE_{mod}$의 평균($172.7 \mu\text{mol CO}_2 \text{ mol}^{-1} \text{ H}_2\text{O}$)이 $iWUE_{obs}$($98.9 \mu\text{mol CO}_2 \text{ mol}^{-1} \text{ H}_2\text{O}$)보다 높았다. 이는 모델이 연간 목질 형성 과정을 통합하는 나이테와 달리, 여름철 정점 조건을 포착하기 때문인 것으로 분석되었다.
• 생태계 규모의 괴리: 관측된 생태계 WUE와 모델링된 WUE 사이에서 뚜렷한 괴리가 관찰되었다. $WUE_{eco}^{obs}$는 유의미한 양의 추세(약 55% 증가)를 보인 반면, $WUE_{eco}^{mod}$는 안정적인 상태를 유지했다. 모델은 ET를 약 50% 과소 추정하여 $WUE_{eco}^{mod}$ 값을 부풀리는 결과를 초래했다. 이러한 과소 추정은 모델의 단일 종 구성에 기인하며, 공동 출현 종(예: 물푸레나무, 가문비나무), 어린 나무, 그리고 하층 식생(특히 Allium ursinum) 등 생태계 전체의 수분 손실에 크게 기여하는 요소들이 제외되었기 때문이다.
• 생리적 동인: 삼중 동위원소 분석 결과, $iWUE$의 장기적 증가는 지속적인 기공 폐쇄보다는 강화된 광합성 동화($A$)에 의해 주도되었다. 이 결론은 $\delta^{18}O$에서 양의 추세(기공 폐쇄에 의한 증발 농축을 나타냄)가 나타나지 않았다는 점과, $\delta^{18}O$ 및 $\delta^{2}H$의 변동성이 잎 수준의 생리적 조절보다는 수원 역학에 의해 주로 제어되었다는 사실에 의해 뒷받침되었다.
• 시간적 척도: $WUE_{eco}$에 대한 관측치와 모델 간의 일치도는 시간적 집계에 따라 개선되었다. 일간 규모에서는 상관관계가 약했으나 월간 규모에서는 강화되었는데, 이는 고빈도 기상 노이즈와 측정 불확실성이 단기적 결합을 가리는 반면, 계절적 패턴은 모델에 의해 더 신뢰성 있게 포착됨을 시사한다.

의의 및 주장
본 논문은 나이테 동위원소, 에디 공분산, 프로세스 기반 모델링을 통합함으로써 잎 수준과 생태계 규모의 WUE를 조화시키기 위한 최초의 메커니즘적 프레임워크를 제공한다고 주장한다. 저자들은 다음과 같이 단언한다:

1. 규모 의존성은 근본적이다: $iWUE$와 $WUE_{eco}$는 서로 교환 가능한 지표가 아니다. 이들은 서로 다른 과정, 공간 영역, 그리고 시간적 통합을 반영한다. 수분 분할과 군집 구성을 고려하지 않은 채 나이테 $iWUE$ 추세를 생태계 규모의 탄소-수분 교환으로 직접 외삽하는 것은 타당하지 않다.
2. 메커니즘적 귀속: 본 연구는 탄소-산소 동위원소 신호의 디커플링을 통해, 이 온대림에서의 장기적인 $iWUE$ 증가가 기공의 가뭄 회피가 아닌 $CO_2$ 시비 효과에 의한 광합성 강화에 의해 주도된다는 것을 입증하였다.
3. 모델의 한계: 생태계 WUE의 괴리는 프로세스 기반 모델이 전체 생태계 수분 플럭스를 정확하게 시뮬레이션하기 위해 다종 캐노피, 하층 식생, 명시적인 수분 분할을 반드시 포함해야 할 필요가 있음을 강조한다.
4. 방법론적 통합: 본 연구는 모델의 종 표현 및 플럭스 분할에 관한 한계를 명시적으로 인정하는 조건 하에, 프로세스 기반 모델을 관측 데이터를 해석하는 "메커니즘적 가교"로 사용하는 것을 검증한다.

저자들은 나이테 동위원소와 EC 관측치가 탄소-수분 순환의 서로 다른 측면을 포착하지만, 이들의 통합이 상승하는 대기 중 $CO_2$와 증가하는 대기 수분 요구에 대한 산림의 반응을 예측하는 데 있어 강력한 프레임워크를 제공한다고 결론짓는다.