Parcel-Level Deforestation Hazard Rates in the Paraguayan Chaco: A Multi-Scale Empirical Analysis of Annual Clearing Dynamics on Large Cattle Ranches

기술 요약: 파라과이 차코(Chaco) 지역의 필지 단위 산림 벌채 위험률

문제 제기
본 논문은 검증된 탄소 표준(VCS) 모듈 VMD0006(BL-PL)에 따라 도출된 필지 단위 산림 벌채 위험률과 파라과이 차코 지역의 지역적 산림 벌채 통계 사이에서 인지되는 불일치를 다룬다. 비판론자들은 프록시(proxy) 필지를 사용하여 구축된 REDD+ 베이스라인이 과도한 산림 손실을 나타낼 수 있다고(베이스라인 과다 크레딧 발행) 주장한다. 이는 활성 필지(active parcels)에서 관찰되는 연간 위험률(예: 연간 약 15%)이 수십 년간의 지역 평균(예: 연간 약 2~4%)보다 현저히 높게 나타나기 때문이다. 핵심 문제는 분모의 불일치에 있다. 필지 단위의 비율은 개발이 활발히 진행 중인 단계의 총 필지 면적을 기준으로 계산되는 반면, 지역 단위의 비율은 활성 및 비활성 자산 전체에 걸쳐 산림 벌채를 희석하여 계산하기 때문이다. 본 연구는 이러한 비율들이 구조적으로 불일치하는 것인지, 아니면 명시적인 스케일링(scaling)을 통해 해결 가능한 것인지를 규명하고자 한다.

방법론
저자들은 원격 탐사, 지적 데이터, 그리고 통계적 스케일링을 결합한 다중 규모 실증 분석을 수행한다:

1. 데이터 소스: 주요 데이터셋은 2012년부터 2020년까지 파라과이 차코 서부 지역(북부 유네스코 생물권 보전지역 제외)의 9개 목장 부지에 걸친 22개의 프록시 필지로 구성된다. 검증을 위해 24개의 추가 지적 필지를 포함한 46개의 확장 샘플이 사용된다.
2. 비율 도출: 연간화된 필지 위험률($D_{parcel}$)은 VMD0006 BL-PL 공식을 사용하여 계산된다: $D_{pn} = D\%_{pn} / Yrs_{pn}$. 여기서 $D\%$는 총 필지 면적에 대한 누적 산림 벌채 백분율이며, $Yrs$는 벌채가 0보다 엄격히 큰 해의 수이다.
3. 스케일링 항등식: 본 논문은 필지와 지역 규모를 조정하기 위한 결정론적 회계 항등식을 개발한다:
   $$r_{region,conv} = (f / g) \times D_{parcel}$$
   여기서:
   • $r_{region,conv}$는 법적 전환 가능 산림 대비 지역 연간 산림 벌채율이다.
   • $g$는 법적 전환 가능 산림의 비율이다. Chaco Vivo 참조 사례의 경우 $g = 0.575$이며, 지역 총합 스케일링 시나리오를 위해 지역 산림 총량에 0.70의 계수를 적용하여 중앙 전환 면적($A_{conv}$)을 도출한다.
   • $f$는 활성 필지 비율(특정 해에 실제로 벌채가 진행 중인 필지의 비중)이다.
   • $D_{parcel}$은 필지 단위 위험률이다.
4. 통계 분석: 이산 시간 위험 분석(discrete-time hazard analysis)을 사용하여 시간적 벌채 패턴을 모델링하고, $g$, $f$, $D_{parcel}$의 불확실성을 전파하기 위해 몬테카를로 시뮬레이션(10,000회 추출)을 수행하며, 평균 위험률의 신뢰 구간을 생성하기 위해 부트스트랩 재표본 추출을 활용한다.
5. 검증: 스케일링 항등식으로부터 도출된 암묵적 활동 비중($f$)을 확장된 46개 필지 유니버스에서 직접 측정한 9년 평균 활동 비율과 비교한다.

주요 결과

• 필지 단위 비율: 22개 프록시 필지의 평균 연간화된 필지 위험률은 연간 14.99% (95% 부트스트랩 CI: 9.86% ~ 20.68%)이다. 분포는 오른쪽으로 치우쳐 있으며 상당한 이질성(SD = 13.15)을 보인다.
• 시간적 역학: 벌채는 초기에 집중되는 경향이 있다. 이산 시간 위험 분석 결과, 평균 위험은 1년 차의 0.28에서 3년 차의 0.05로 감소한다(81% 감소, $p=0.006$). 중앙값 필지는 관찰된 벌채의 89%를 시작 후 2년 이내에 완료한다.
• 스케일링 조정:
  • Chaco Vivo 참조 값인 $g = 0.575$(또는 지역 총 산림량에 적용된 0.70 계수)와 필지 비율 14.99%를 스케일링 항등식에 적용했을 때, 역사적 지역 산림 벌채량(160,000~300,000 ha/yr)과 일치하기 위해 필요한 암묵적 활동 비중($f$)은 0.08에서 0.18 사이로 나타난다.
  • 몬테카를로 역식별(inverse identification)을 통해 도출된 중앙값 암묵적 $f$는 0.12 (95% 구간: 0.07–0.19)이다.
  • 이는 확장된 46개 필지 세트에서 (연간 벌채 임계값 5% 기준) 직접 측정한 9년 평균 활동 비율인 0.121과 일치한다.
• 확정적 전환(Committed Conversion): 확장된 46개 필지 세트를 "확정적 전환"(누적 벌채 $\ge$ 20%) 조건으로 필터링했을 때, 결합 평균 비율은 15.78%/yr로, 기존의 14.99% 추정치를 포함한다.
• 지역적 투영: 전방향 불확실성 실험에 따르면 평균 암묵적 지역 산림 벌채량은 약 238,000 ha/yr이며, 시뮬레이션 추출값의 58.4%가 역사적 범위인 160,000~300,000 ha/yr 내에 존재한다.

주요 기여

1. 공식적 조정: 본 논문은 필지 규모의 위험 강도와 지역적 벌채 규모를 수학적으로 연결하는 결정론적 스케일링 항등식을 제공함으로써, 겉으로 보이는 불일치가 위험의 부풀리기가 아니라 분모의 정의(총 면적 대 전환 가능 산림)와 활동 혼합(활성 필지 대 비활성 필지)의 함수임을 입증한다.
2. $f$의 실증적 검증: 본 연구는 활동 필지 비율($f$)을 이론적 잔차에서 실증적으로 측정 가능한 양으로 격상시켰으며, 매년 약 12%의 전환 가능 자산이 활성 벌채 단계에 있다는 것을 보여줌으로써 높은 필지 비율과 지역 총량을 일치시킨다.
3. 시간적 특성 규명: 본 연구는 차코 지역 목장 개발의 "초기 집중적(front-loaded)" 특성을 정량화하여, 높은 연간화된 비율이 지속적인 장기 벌채가 아닌 짧고 강렬한 벌채 펄스(1~2년)에 의해 발생함을 보여준다.

의의 및 주장
저자들은 본 연구가 필지 단위 위험률이 반드시 과장된 지역 산림 벌채 위험을 의미하지 않는다는 점을 보여줌으로써 "베이스라인 과다 크레딧 발행" 논쟁을 해결한다고 주장한다. 논문의 핵심 주장은 다음과 같다:

• 14.99%/yr의 비율은 활성 개발 단계에 특화된 조건부 강도이며, 스케일링 없이는 지역 평균과 직접 비교할 수 없다.
• 법적 전환 가능 비율($g$)과 활성 필지 비중($f$)을 명시적으로 고려할 때, 필지 단위 비율은 역사적 지역 벌채 기록과 구조적으로 일치한다.
• 본 프레임워크는 공개된 원격 탐사, 지적 경계, 승인된 토지 이용 계획을 사용하여 재현 가능하다.
• 본 연구는 국가 전체를 대표하는 센서스라고 주장하지 않는다. 샘플 크기는 작고($N=22$) 활성 전선에 집중되어 있다. 그러나 프록시로부터 도출된 비율과 지역 데이터 사이의 구조적 일관성은 VMD0006 하에서의 베이스라인 구축을 위한 방법론이 타당함을 시사한다.

결론적으로, 본 논문은 두 수치 사이의 차이가 적절한 스케일링 항등식이 적용될 경우 단순한 집계상의 회계적 아티팩트(artifact)이지, REDD+ 베이스라인의 방법론적 실패 증거가 아니라고 밝히고 있다.