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スコアを持つ3つの独立変数、すなわち相対湿度 ($X_1$)、降水量 ($X_2$)、および平均土壌温度 ($X_3$) が保持されました。
  4. モデル構築: 選択された気候因子に基づいて土壌VWCを予測する重回帰（MLR）モデルを開発しました。
  5. モデル検証: 残差プロット（適合値 vs 残差）および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)$。
  • 有意性: 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年までに降水量が10〜20%減少し、気温が2〜3℃上昇するという気候予測が進む中で、水資源管理と農業計画を支援するものであると断言しています。本研究は、モデルが変動の半分近くを説明していること（環境研究において良好な値とされる）に触れつつ、降水の直接的な影響と温度の逆の影響を効果的に浮き彫りにしていると結論付けています。著者らは、これらの主要な変数を将来の乾燥地域における評価に使用することを推奨し、予測精度をさらに向上させるために、より高度なモデル（例：時系列モデルや非線形モデル）の使用を提案しています。