Agent-Based Modeling of Low-Emission Fertilizer Adoption for Dairy Farm Decarbonisation using Empirical Farm Data

技術概要：牛乳生産の脱炭素化に向けた低排出肥料の導入に関するエージェントベースモデル

問題定義
農業システム、特にアイルランドの酪農は、主に窒素肥料の使用、家畜管理、土地利用慣行に起因する温室効果ガス（GHG）排出の主要な寄与者である。保護尿素（尿素分解酵素阻害剤 NBPT で処理された尿素）の導入など、実用的な低排出対策は存在するが、農場間での導入にはばらつきがある。この不均質性は、管理能力の違い、構造的制約、社会的相互作用に起因する。従来のトップダウン型の評価モデルは、均質な行動と即時的な導入を前提とした仮定に依存し、同僚効果、経路依存性、確率的な意思決定を無視しているため、この複雑さを捉えきれない。したがって、導入を単なる財務計算ではなく、社会的感染と構造的な農場特性に影響される動的プロセスとして扱う、脱炭素化の社会技術的進化をシミュレーション可能な枠組みが必要である。

方法論
本研究は、15 年間の期間にわたり 295 件のアイルランド酪農農場における低排出技術の普及をシミュレーションする、経験的に裏付けられたエージェントベースモデル（ABM）枠組みを開発した。

• データ基盤: モデルは、アイルランド家畜育種連合（ICBF）および Bord Bia から集約された、アイルランドの農業・食料開発庁（Teagasc）が提供する包括的な農場レベルのデータセットを利用する。このデータセットには、生産性能（乳量、固形分、畜群構造）、投入資材の使用（肥料の種類、石灰、飼料）、および経済変数が含まれる。
• エージェント構造: 295 件の農場のそれぞれは、土地面積、家畜単位、乳生産量、窒素投入量という特定の属性を持つ自律的なエージェントとしてモデル化される。初期状態では、すべてのエージェントが従来の硝酸アンモニウムカルシウム（CAN）を使用する。
• 社会ネットワーク構造: 農業討論グループの社会的ダイナミクスを再現するため、モデルはワッツ＝ストロガッツの小世界ネットワークを採用する。この構造は、高い局所的クラスタリング（安定した同僚助言サークルを表す）と短い平均経路長（アドバイザーを介したグループ間での情報伝達を表す）を特徴とする。ネットワークパラメータ（$k=4$ の接続、再結合確率 $p=0.1$）は、特定の導入曲線に適合させるのではなく、疎でクラスタリングされた同僚構造を表すために固定された。
• 導入メカニズム: 保護尿素の導入は、確率的な意思決定プロセスとしてモデル化される。時間 $t$ における非導入農場 $i$ の導入確率（$P_{i,t}$）は、以下の関数である：
  1. ベースラインの傾き: 革新への固定的な傾向（$\beta_0$）。
  2. 同僚の影響: 隣接農場の平均導入率を表す遅延信号（$S_{i,t-1}$）を、較正された社会的影響パラメータ（$\omega = 0.85$）で重み付けしたもの。
  3. 構造的要因: 正規化された農場規模と乳生産量。
  4. 政策ショック: 炭素税や補助金に対する外生的な調整。
• 政策シナリオ: モデルは 3 つのシナリオをシミュレーションする。ベースライン（政策なし）、炭素税（埋め込まれた肥料排出に適用される 1 トンあたり 71 ユーロの CO2 税）、および補助金（保護尿素に対して 1 トンあたり 200 ユーロ）である。
• 排出量と経済性: モデルは、肥料固有の排出係数（CAN: 0.0149、保護尿素: 0.0040）に基づいて直接的な亜酸化窒素排出量を計算し、CO2 換算値に変換する。また、私的投入コストの変化と総削減量を比較して限界削減費用（MAC）を計算し、私的コストと政府支出を含む社会的コストを区別する。
• 検証と不確実性: モデルは、歴史的な導入データ（2019–2024 年）を用いて較正され、2023–2024 年のサンプル外データに対して、二乗平均平方根誤差（RMSE）、平均絶対誤差（MAE）、および決定係数（$R^2$）指標を用いて検証された。堅牢性は、250 回のモンテカルロ反復と、社会的影響の重みおよび補助金レベルに関する感度分析を通じて評価された。

主な貢献

1. 経験的較正: 本研究は、現実世界の農場レベルデータを較正された ABM に接続し、高い予測精度（$R^2 = 0.979$、RMSE = 0.0274）で肥料導入の現実的な普及経路をモデル化することに成功した。
2. 生物物理的・経済的モデルの統合: 経済的コスト計算と生物物理的排出量モデルを組み合わせ、シミュレーション環境内で限界削減費用を推定し、私的費用と社会的費用を区別する。
3. 不確実性の定量化: モンテカルロ法を利用することで、導入時期と累積排出量における不確実性を定量化し、決定論的な単一点推定ではなく確率的な洞察を提供する。
4. 政策実験室枠組み: このモデルは、気候戦略の堅牢性と普及速度を事前にテストする「政策実験室」として機能し、脱炭素化を純粋な金銭的計算ではなく、社会技術的進化として概念化する。

結果

• 予測精度: モデルは高い予測精度を示し、決定係数（$R^2$）は 0.9793、RMSE は 0.0274 であった。社会ネットワーク効果の取り込みは決定的であった。これを取り除くと RMSE は 5 倍に増加（0.0274 から 0.1476）し、同僚間の影響が肥料導入における支配的なメカニズムであることを確認した。
• 導入ダイナミクス: 導入はロジスティック S 字曲線に従った。ベースラインシナリオでは、約 7.2 年で 50% の導入に達した。炭素税はこれを約 3.6 年に加速し、補助金シナリオは最も迅速で、約 3.0 年で 50% に達した。
• 累積削減: 15 年間の期間において、補助金シナリオは 121.35 Gg CO2 換算の回避を達成し、炭素税シナリオ（97.73 Gg CO2 換算）より約 24% 高かった。この差は、補助金による早期の導入ピークに起因し、これが排出削減の期間を延長するためである。
• 分布のシフト: コルモゴロフ＝スミルノフ検定により、政策介入が部門の炭素強度分布を根本的に変化させることが確認された（$D = 0.2407, p < 0.001$）。これにより高強度の尾部が効果的に抑制され、「遅れ」農場が基準モードへと引きずり込まれた。
• 経済的妥当性: 私的限界削減費用は 1 トン CO2 換算あたり 0.80 ユーロと計算され、現在の EU 炭素価格基準より著しく低かった。しかし、政府の補助金支出を含めると、社会的削減費用は 1 トン CO2 換算あたり 4.54 ユーロに上昇した。
• 不均質性: 構造分析により、炭素強度は農場規模だけでなく、栄養利用効率と放牧バランスによって駆動されることが明らかになった。「遅れ」農場（排出量の上位 20%）は、乳量の増加に対応しない高い放牧率（約 4.5 頭/ヘクタール）を示し、逓減収益と窒素損失を示唆した。

意義と主張
本論文は、その主な意義が、農業システムにおける不均質でネットワーク媒介型の技術普及を調査するための、検証済みかつ一般化可能な枠組みの提供にあると主張する。経験的データとエージェントベースモデルを統合することで、本研究は以下のことを実証する。

• 社会的感染の重要性: 同僚の影響は二次的要因ではなく、導入速度の主要な駆動力である。ネットワークダイナミクスを考慮しない政策は、普及速度や初期導入者の影響を過小評価する可能性がある。
• 補助金対炭素税: 炭素税も補助金も導入を促進し得るが、補助金ベースの介入は、移行の初期段階で自己強化型の社会的臨界点をより早く引き起こす能力により、より迅速な普及とより大きな累積削減を生み出す。
• 経路依存性: モデルは、政策のタイミングがネットワーク感染と相互作用して持続的な移行経路を創出することを示す。初期のインセンティブは、同僚による検証を通じて知覚されるリスクを軽減することで行動を「固定化」する。このダイナミクスは、静的な価格メカニズム（例えば税）が、同じ社会的閾値を越えるのに十分なほど強力でない限り、複製するのが困難である。
• 標的化された介入: この枠組みは、構造的に影響力が高く排出量の多い農場を標的とすることが最大の限界効果をもたらすことを示唆しており、画一的なアプローチではなく、差別化された証拠に基づく介入の設計を支持する。

著者らは、このモデルを「政策実験室」の実験ツールとして位置づけ、ステークホルダーが実施前に脱炭素戦略を予測し最適化することを可能にする一方、限定された行動の複雑さと静的な経済変数に関する限界を認識している。