Modeling Soil Volumetric Water Content Based on Climatic Factors in the Tripoli Region (Libya) Using Multiple Linear Regression

기술 요약: 트리폴리(리비아)의 토양 체적 함수율 모델링

문제 제기
토양 수분은 지구의 물과 에너지 순환의 핵심 요소로, 농업 생산성, 수자원 관리, 그리고 가뭄 및 홍수와 같은 극한 기상 현상의 빈도에 영향을 미칩니다. 반건조 내지 건조 기후를 특징으로 하는 리비아에서 토양 수분 역학을 이해하는 것은 식량 안보와 재난 위험 완화에 필수적입니다. 그러나 현장 측정을 통한 장기적인 전 지구적 토양 수분 모니터링은 제한적입니다. 원격 탐사가 대안을 제공할 수는 있지만, 통계적 방법을 사용하여 국지적 기후 요인과 토양 체적 함수율(VWC) 사이의 관계를 정량적으로 모델링할 필요가 있습니다. 본 연구는 북아프리카의 강수량 감소 및 기온 상승을 나타내는 기후 변화 전망의 맥락에서, 특정 기후 변수가 트리폴리 지역의 토양 수분에 어떻게 영향을 미치는지에 대한 정량적 분석의 공백을 다룹니다.

연구 방법론
본 연구는 ERA5 데이터셋(유럽 중기 예보 센터)에서 추출한 2001년부터 2024년까지의 297개 월간 데이터를 활용하였습니다. 분석은 R 통계 소프트웨어를 사용하여 수행되었습니다.

• 변수:
  • 종속 변수 ($Y$): 4개 토양 깊이에서의 평균 체적 토양 함수율.
  • 초기 독립 변수: 상대 습도, 강수량, 평균 토양 온도(4개 깊이 전체), 기온, 풍속, 잠재 증발량, 지하수 깊이.
• 통계적 접근 방식:
  1. 기술 통계 분석: 모든 변수의 분포와 중심 경향성을 조사하였습니다.
  2. 다중공선성 평가: 상관 행렬과 분산 팽창 요인(VIF) 분석을 수행하였습니다. 초기 분석 결과, 여러 변수(예: 기온, 토양 온도, 잠재 증발량) 사이에서 심각한 다중공선성이 발견되었습니다(VIF 값이 허용 임계치인 5를 훨씬 초과하는 12 이상으로 나타남).
  3. 변수 선택: 모델의 안정성을 확보하기 위해 높은 VIF 값을 가진 변수들을 제외하였습니다. 최종 모델에는 허용 가능한 VIF 점수를 가진 세 가지 독립 변수인 상대 습도($X_1$), 강수량($X_2$), 평균 토양 온도($X_3$)가 포함되었습니다.
  4. 모델 구축: 선택된 기후 요인을 기반으로 토양 VWC를 예측하기 위한 다중 선형 회귀(MLR) 모델을 개발하였습니다.
  5. 모델 검증: 잔차 도표(잔차 대 적합값)와 Q-Q 도표를 사용하여 잔차의 등분산성과 정규성을 확인하여 가설을 검증하였습니다.

주요 결과

• 기술 통계: 평균 체적 토양 함수율은 0.02로 표준 편차가 낮았으며, 이는 일관성을 나타냅니다. 강수량은 변동성이 높았고(평균 0.92 mm), 상대 습도는 평균 44.10%였습니다.
• 상관 분석:
  • 토양 VWC와 상대 습도 사이에 강한 양의 상관관계가 발견되었습니다 ($r = 0.638$).
  • 토양 VWC와 강수량 사이에 중간 정도의 양의 상관관계가 존재했습니다 ($r = 0.528$).
  • 토양 VWC와 평균 토양 온도 사이에 중간 정도의 음의 상관관계가 관찰되었습니다 ($r = -0.559$).
• 회귀 모델 성능:
  • 최종 모델 방정식은 다음과 같습니다: $\hat{Y} = 0.01378 + (7.0 \times 10^{-5} \times X_1) + (4.22 \times 10^{-4} \times X_2) - (9.8 \times 10^{-5} \times X_3)$.
  • 유의성: 세 가지 예측 변수 모두 통계적으로 유의했습니다 ($p < 0.001$).
    • 강수량: 가장 영향력 있는 요인으로 식별되었습니다 ($\beta = 4.22 \times 10^{-4}$). 강수량이 1 mm 증가하면 평균 토양 수분은 0.000422 직접 증가합니다.
    • 상대 습도: 증발 감소에 기인하여 양(+)의 효과를 보입니다 ($\beta = 7.0 \times 10^{-5}$).
    • 토양 온도: 증발 및 침투 증가에 기인하여 유의미한 음(-)의 효과를 보입니다 ($\beta = -9.8 \times 10^{-5}$).
  • 설명력: 조정된 $R^2$ 값은 0.477로, 이는 모델이 토양 체적 함수율 변동의 약 48%를 설명함을 의미합니다.
  • 전반적인 적합도: F-검정은 모델의 전반적인 통계적 유의성을 확인했습니다 ($F \approx 91$, $p < 0.001$). 진단 도표를 통해 잔차가 무작위로 분포하고 정규 분포를 따름을 확인하여 모델의 가설을 검증하였습니다.

의의 및 주장
본 논문은 본 연구가 선택된 기후 요인과 트리폴리 지역의 토양 체적 함수율 사이의 통계적으로 유의미한 관계를 성공적으로 구축했음을 주장합니다. 강수량, 상대 습도, 토양 온도를 주요 동인으로 분리함으로써, 본 연구는 반건조 환경에서 토양 수분 역학을 이해하기 위한 예측 도구를 제공합니다.

저자들은 이 모델링 접근 방식이 2050년까지 강수량 1020% 감소 및 기온 23°C 상승을 예측하는 북아프라카의 기후 전망 상황에서 수자원 관리 및 농업 계획을 지원한다고 단언합니다. 본 연구는 모델이 변동성의 거의 절반을 설명하며(환경 연구에서 양호한 수치로 간주됨), 강수의 직접적인 영향과 온도의 역(逆) 영향을 효과적으로 강조한다는 결론을 내립니다. 저자들은 건조 지역의 향후 평가를 위해 이러한 핵심 변수들을 사용할 것을 권장하며, 예측 정확도를 더욱 향상시키기 위해 더 발전된 모델(예: 시계열 또는 비선형 모델)을 사용할 것을 제안합니다.