論文要約：多変量格子変形モデルによる作物輪作の土壌栄養動態への影響評価

タイトル: Multivariate Lattice Deformation: A Spatially Explicit Framework for Assessing Crop Rotation Impacts on Soil Nutrient Dynamics
著者: Marco Mandap
発表日: 2026 年 3 月

1. 背景と課題 (Problem)

持続可能な農業において、作物輪作は土壌肥沃度を維持するための重要な戦略ですが、その効果を評価する既存の方法には以下の重大な限界があります。

• 空間的均質性の仮定: 従来の評価は圃場平均値や単一の栄養素（N, P, K のいずれか）に焦点を当てており、土壌の空間的異質性（テクスチャ、鉱物組成による緩衝能の違い）を無視しています。
• 多変量相互作用の欠如: 作物は N, P, K を同時に吸収・利用しますが、多くのモデルはこれらを独立して扱っており、多変量な相互作用や、特定の栄養素の枯渇がシステム全体に与える複合的な影響を見逃しています。
• 空間的パターンと管理の乖離: 土壌栄養素はメーターからキロメートルスケールで複雑な空間パターンを示しますが、均質な試験区で得られた結果が不均質な農地全体に適用されることで、過剰施肥や不足施肥が発生しています。

これらの課題を解決し、空間的・多変量的な視点から作物輪作の持続可能性を予測・評価するための新たな枠組みが必要です。

2. 手法と枠組み (Methodology)

本研究は、物理的な格子（ラティス）の变形にアナロジーを用いた**「多変量格子変形モデル（Multivariate Lattice Deformation Model）」**を提案しました。このモデルは、土壌を 4 次元テンソルとして表現し、以下の要素を統合しています。

2.1 モデルの基本概念

• 4 次元テンソル表現: 土壌状態を $S(x, y, t, c)$ で定義します。
  • $x, y$: 空間座標（20×20 のグリッド）
  • $t$: 時間（作物輪作サイクル）
  • $c$: 栄養素チャネル（N, P, K）
• 作物の力ベクトル: 各作物は、N, P, K に対する同時的な需要（または供給、例えばマメ科作物の窒素固定）を表す「力ベクトル」として作用します。
• 空間的変異（剛性）: 土壌の緩衝能（栄養変化に対する抵抗性）を「剛性（Stiffness）」として表現します。砂質土（剛性低）は変化しやすく、粘土質土（剛性高）は変化しにくいと設定し、これが空間的に変化するパラメータとして組み込まれます。
• 空間平滑化（カーネル畳み込み）: 拡散、土壌移動、耕起による混合を表現するため、ガウスカーネルによる平滑化を適用し、現実的な空間的自己相関を導入します。

2.2 評価指標と統計的検定

• 多変量ストレス指標（Deformation Energy）: 初期状態からの現在の N-P-K 空間におけるユークリッド距離を計算し、累積的な栄養状態の変化を単一のスカラー値（$D$）として定量化します。
• 統計的有意性: クラム＝フォン・ミーゼス（Cramér–von Mises）検定を用いて、初期分布と最終分布の差異がランダムな変動ではなく、体系的な劣化であることを非パラメトリックに検証します。
• 空間的自己相関: モランの I（Moran's I）を用いて、生成された空間パターンの現実性を評価します。

2.3 シミュレーション設定

• 対象: 3 年輪作（トウモロコシ→ダイズ→コムギ）を 1 サイクル実施。
• 比較: 連続トウモロコシ栽培など、異なる管理シナリオとの比較。
• 感度分析: 平滑化スケール、剛性コントラスト、作物の吸収量（力）の大きさを変化させてモデルの挙動を検証。

3. 主要な結果 (Key Results)

3.1 空間的異質性とストレス分布

• 平均ストレス: 1 サイクル後の平均多変量ストレスは 0.63（正規化スケール）でした。
• 剛性の影響: 低緩衝能（砂質）の領域では平均ストレスが 0.84 と高く、高緩衝能（粘土質）の領域では 0.42 と低く、土壌の物理的特性が栄養枯渇の空間パターンを強く制御することが確認されました。
• ホットスポット: 砂質土の領域でストレスが集中し、最大値 0.91 に達しました。

3.2 多変量動態の発見

• リン（P）の枯渇: ダイズの窒素固定（N 増加）にもかかわらず、リンの平均枯渇率（17.9%）は窒素（10.8%）を上回りました。
• 支配的ストレス要因: 全体の 81% のグリッドで窒素が主要なストレス要因でしたが、19% のセルではリンが支配的でした。これは単一栄養素分析では見逃される重要なリスクです。
• 相関の空間的変動: 低緩衝能地域では N, P, K の枯渇が強く相関していましたが、高緩衝能地域では相関が弱く、栄養素ごとの緩衝メカニズムの違いが空間的に現れました。

3.3 輪作多様性の効果

• 連続栽培の影響: 連続トウモロコシ栽培は、多様な 3 年輪作と比較して平均ストレスを 41%、最大ストレスを 27% 増加させました。
• 空間的均質化: 輪作が単純化されると、空間的なストレスのばらつき（変動係数）が減少し、圃場全体が均一的に劣化する傾向が見られました。

3.4 感度分析と統計的検証

• 剛性コントラスト: 土壌の緩衝能の差が小さい（均一な砂質土など）場合、ストレスは空間的に均質化し、平均ストレスレベルが連続栽培に近い高値に達しました。
• 統計的有意性: クラム＝フォン・ミーゼス検定により、すべての栄養素で初期状態からの有意な逸脱（$p \le 0.002$）が検出されました。また、モランの I（0.29–0.30, $p < 0.001$）により、現実的な空間的自己相関が再現されていることが確認されました。

4. 主な貢献と意義 (Contributions and Significance)

4.1 概念的・方法的貢献

• 新たな枠組みの提案: 土壌栄養動態を「多変量な力による格子の变形」として捉えるアナロジーは、空間的・多変量的な相互作用を直感的かつ定量的に理解する新しい視点を提供しました。
• 地理統計とメカニズムモデルの統合: 土壌の空間的変動を記述する地理統計学的アプローチと、作物の影響を記述するメカニズム的アプローチを統合し、両者のギャップを埋めました。
• 定量的な持続可能性評価: クラム＝フォン・ミーゼス検定を用いることで、輪作が「統計的に安定」か「システム劣化を招くか」を仮説検定する形式を提供しました。

4.2 実務的・政策的意義

• 空間的リスクゾーンの特定: 多変量ストレスが臨界値に近づく領域を特定でき、精密農業における変量施肥（Variable Rate Application）や、特定の栄養素に焦点を当てた管理介入の優先順位付けを可能にします。
• 「貧乏がさらに貧しくなる」動態の解明: 初期の土壌特性（緩衝能）が低い地域では、同じ管理でも劣化が加速することを示し、均一な施肥推奨がなぜ失敗するかを説明しました。
• 政策への示唆: 作物多様化の恩恵を空間的・多変量的な指標で定量化でき、持続可能な農業へのインセンティブ設計や認証基準の科学的根拠となり得ます。

5. 結論

本研究で提案された多変量格子モデルは、作物輪作が土壌栄養に与える影響が本質的に「空間的」かつ「多変量的」であることを実証しました。従来の圃場平均ベースの評価では見逃されていた、特定の栄養素の局所的枯渇や、土壌緩衝能によるリスクの空間的偏りを明らかにしました。この枠組みは、計算コストを抑えつつ複雑な土壌 - 作物相互作用を捉えるための実用的なツールとして、持続可能な農業管理の設計と、サイト固有の介入策の立案に重要な指針を提供します。