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

ملخص تقني: النمذجة القائمة على الوكلاء لتبني الأسمدة منخفضة الانبعاثات في إطار إزالة الكربون من مزارع الألبان

بيان المشكلة
تعد الأنظمة الزراعية، ولا سيما مزارع الألبان الأيرلندية، مساهمًا رئيسيًا في انبعاثات غازات الاحتباس الحراري (GHG)، مدفوعة بشكل أساسي باستخدام الأسمدة النيتروجينية، وإدارة الماشية، وممارسات استخدام الأراضي. وبينما توجد تدابير عملية منخفضة الانبعاثات — مثل اعتماد اليوريا المحمية (اليوريا المعالجة بمثبط اليورياز NBPT) — فإن تبنيها عبر المزارع غير متكافئ. وينبع هذا التباين من الاختلافات في القدرة الإدارية، والقيود الهيكلية، والتفاعلات الاجتماعية. غالبًا ما تفشل نماذج التقييم التقليدية من الأعلى إلى الأسفل في استيعاب هذا التعقيد، حيث تعتمد على افتراضات السلوك الموحد والقبول الفوري مع تجاهل تأثيرات الأقران، والاعتماد على المسار، واتخاذ القرار العشوائي. وبناءً على ذلك، تبرز الحاجة إلى إطار عمل قادر على محاكاة التطور الاجتماعي التقني لإزالة الكربون، ومعاملة التبني ليس مجرد عملية حسابية مالية، بل كعملية ديناميكية تتأثر بالعدوى الاجتماعية والخصائص الهيكلية للمزارع.

المنهجية
يطور البحث إطار عمل للنمذجة القائمة على الوكلاء (ABM) مستند إلى أسس تجريبية لمحاكاة انتشار التقنيات منخفضة الانبعاثات عبر 295 مزرعة ألبان أيرلندية على مدى أفق زمني يمتد لـ 15 عامًا.

• الأساس البياني: يستخدم النموذج مجموعة بيانات شاملة على مستوى المزرعة مقدمة من (Teagasc) (الهيئة الأيرلندية لتطوير الزراعة والغذاء)، مجمعة من الاتحاد الأيرلندي لتربية الماشية (ICBF) و(Bord Bia). وتتضمن هذه البيانات أداء الإنتاج (إنتاج الحليب، المواد الصلبة، هيكل القطيع)، والاستخدام المدخلات (أنواع الأسمدة، الجير، الأعلاف)، والمتغيرات الاقتصادية.
• هيكل الوكيل: يتم نمذجة كل مزرعة من المزارع الـ 295 كوكيل مستقل بخصائص محددة: مساحة الأرض، الوحدات الحيوانية، إنتاج الحليب، ومدخلات النيتروجين. في البداية، تستخدم جميع الوكلاء نترات الأمونيوم والكالسيوم (CAN) التقليدية.
• بنية الشبكة الاجتماعية: لمحاكاة ديناميكيات النقاش الزراعي، يستخدم النموذج شبكة "واتس-ستروغاتس" (Watts-Strogatz) ذات العالم الصغير. تتميز هذه البنية بتكتل محلي عالٍ (يمثل دوائر المشورة المستقرة بين الأقران) وأطوال مسارات متوسطة قصيرة (تمثل انتقال المعلومات عبر المجموعات عبر المستشارين). تم تثبيت بارامترات الشبكة ($k=4$ اتصالات، واحتمالية إعادة التوصيل $p=0.1$) لتمثيل بنية أقران متفرقة ومتكتلة بدلاً من ضبطها لتطابق منحنيات تبني محددة.
• آلية التبني: يتم نمذجة تبني اليوريا المحمية كعملية قرار احتمالية وعشوائية. يكون احتمال التبني ($P_{i,t}$) لمزرعة غير متبنية $i$ عند الوقت $t$ دالة لـ:
  1. الميل الأساسي: نزعة ثابتة للابتكار ($\beta_0$).
  2. تأثير الأقران: إشارة متأخرة ($S_{i,t-1}$) تمثل متوسط معدل التبني للمزارع المجاورة، موازنة بمعامل تأثير اجتماعي مُعاير ($\omega = 0.85$).
  3. العوامل الهيكلية: حجم المزرعة وإنتاج الحليب بعد المعايرة.
  4. صدمات السياسة: تعديلات خارجية لضرائب الكربون أو الإعانات.
• سيناريوهات السياسة: يحاكي النموذج ثلاثة سيناريوهات: السيناريو الأساسي (لا توجد سياسة)، ضريبة الكربون (71 يورو/طن من انبعاثات الأسمدة المضمنة)، والإعانة (200 يورو/طن على اليوريا المحمية).
• الانبعاثات والاقتصاد: يحسب النموذج انبعاثات أكسيد النيتروز المباشرة بناءً على عوامل انبعاث خاصة بالأسمدة (CAN: 0.0149؛ اليوريا المحمية: 0.0040) ويحولها إلى مكافئات ثاني أكسيد الكربون. كما يحسب تكاليف التخفيف الهامشية (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 بمقدار خمسة أضعاف (من 0.0274 إلى 0.1476)، مما يؤكد أن تأثير الأقران هو آلية مهيمنة في تبني الأسمدة.
• ديناميكيات التبني: اتبع التبني منحنى "S" لوجستي. وصل السيناريو الأساسي إلى 50% من التبني في حوالي 7.2 سنة. وقد سرّعت ضريبة الكربون هذا المعدل إلى ~3.6 سنة، بينما كان سيناريو الإعانة هو الأسرع، حيث وصل إلى 50% في ~3.0 سنوات.
• التخفيف التراكمي: خلال الأفق الزمني لـ 15 عاماً، حقق سيناريو الإعانة تجنباً قدره 121.35 جيجا جرام من مكافئ ثاني أكسيد الكربون، وهو أعلى بنسبة 24% تقريباً من سيناريو ضريبة الكربون (97.73 جيجا جرام. ويعزى هذا الفرق إلى وصول ذروة التبني مبكراً في سيناريو الإعانة، مما يطيل فترة توفير الانبعاثات.
• التحولات التوزيعية: أكد اختبار "كولموغوروف-سميرنوف" أن التدخلات السياساتية تغير بشكل جوهري توزيع كثافة الكربون في القطاع ($D = 0.2407, p < 0.001$)، مما يؤدي فعلياً إلى كبح ذيل الكثافة العالية وسحب المزارع "المتأخرة" نحو المنوال المرجعي.
• الجدوى الاقتصادية: تم حساب تكلفة التخفيف الهامشية الخاصة عند 0.80 يورو لكل طن من مكافئ ثاني أكسيد الكربون، وهو أقل بكثير من معايير أسعار الكربون في الاتحاد الأوروبي. ومع ذلك، عند تضمين نفقات الإعانة الحكومية، ارتفعت تكلفة التخفيف الاجتماعية إلى 4.54 يورو لكل طن من مكافئ ثاني أكسيد الكربون.
• التباين: كشف التحليل الهيكلي أن كثافة الكربون مدفوعة بكفاءة استخدام المغذيات وتوازن التخزين أكثر من حجم المزرعة وحده. أظهرت المزارع "المتأخرة" (أعلى 20% من الملوثين) معدلات تخزين عالية (~4.5 وحدة حيوانية/هكتار) دون مكاسب مقابلة في إنتاج الحليب، مما يشير إلى تناقص العوائد وفقدان النيتروجين.

الأهمية والادعاءات
يزعم البحث أن أهميته الأساسية تكمن في توفير إطار عمل مُثبت وقابل للتعميم لاستقصاء انتشار التكنولوجيا المتغاير والمدفوع بالشبكات في الأنظمة الزراعية. ومن خلال دمج البيانات التجريبية مع النمذجة القائمة على الوكلاء، يوضح البحث أن:

• العدوى الاجتماعية أمر حاسم: تأثير الأقران ليس عاملاً ثانوياً بل هو محرك رئيسي لسرعة التبني. فالسياسات التي تفشل في مراعاة ديناميكيات الشبكة قد تقلل من تقدير سرعة الانتشار أو تأثير المتبنين الأوائل.
• الإعانات مقابل الضرائب: بينما يمكن لكل من ضرائب الكربون والإعانات دفع التبني، فإن التدخلات القائمة على الإعانات تخلق انتشاراً أسرع وتخفيفاً تراكمياً أكبر بسبب قدرتها على تحفيز نقاط التحول الاجتماعي ذاتية التعزيز في وقت مبكر من عملية الانتقال.
• الاعتماد على المسار: يوضح النموذج أن توقيت السياسة يتفاعل مع العدوى الشبكية لخلق مسارات انتقال مستدامة. يمكن للحوافز المبكرة أن "تثبّت" السلوك من خلال تقليل المخاطر المدركة عبر التحقق من صحة الأقران، وهي ديناميكية قد تجد الآليات السعرية الثابتة (مثل الضرائب) صعوبة في محاكاتها ما لم تكن قوية بما يكفي لتجاوز نفس العتبة الاجتماعية.
• التدخلات المستهدفة: يشير الإطار إلى أن استهداف المزارع ذات التأثير الهيكلي العالية وعالية الانبعاثات يحقق أكبر أثر هامشي، مما يدعم تصميم تدخلات متميزة قائمة على الأدلة بدلاً من النهج الموحد للجميع.

يضع المؤلفون النموذج كأداة لـ "مختبرات السياسة" التجريبية، مما يسمح لأصحاب المصلحة باستباق وتحسين استراتيجيات إزالة الكربون قبل تنفيذها، مع الاعتراف بالقيود المتعلقة بالتعقيد السلوكي المحدود والمتغيرات الاقتصادية الثابتة.