Parcel-Level Deforestation Hazard Rates in the Paraguayan Chaco: A Multi-Scale Empirical Analysis of Annual Clearing Dynamics on Large Cattle Ranches
تحلل هذه الورقة ديناميكيات إزالة الغابات على مستوى القطع الأرضية في منطقة تشاكو الباراغوانية لتوضيح كيف تؤدي دورات الإزالة المضغوطة في مزارع الماشية النشطة إلى توليد معدلات خطر سنوية عالية، والتي عندما يتم قياسها بنسبة القطع الأرضية التي تتحول فعلياً، تتوافق مع الإحصاءات الإقليمية الأوسع وتساهم في حسابات الخط المرجعي لمشاريع الحد من انبعاثات إزالة الغابات وتدهورها (REDD+).
المؤلفون الأصليون:William C. Schuman, Ionuț-Alexandru Ciocîrlie
في الغابات الجافة الشاسعة في تشاكو باراغواي، نما توتر هادئ بين طريقتين لحساب الأشجار المفقودة. فمن ناحية، تُظهر الخرائط الإقليمية تقلص الغابة بوتيرة ثابتة ومقدور عليها، حيث تختفي آلاف الهكتارات كل عام مع توسع مزارع الماشية. ومن ناحية أخرى، تروي صور الأقمار الصناعية للممتلكات الزراعية الفردية قصة مختلفة: فعندما يقرر صاحب مزرعة إزالة الأرض، فإنه يفعل ذلك بسرعة مذهلة، حيث يزيل في كثير من الأحيان جزءًا ضخمًا من الغابة في عام واحد فقط. وقد أدى هذا التباين إلى جدل جاد حول مشاريع ائتمان الكربون المصممة لحماية هذه الغابات. ويجادل النقاد بأنه إذا نظرت إلى مدى سرعة قيام كل مزرعة بإزالة أرضها، فقد تستنتج أن المنطقة بأكملها تتعرض للدمار بسرعة أكبر بكثير مما هي عليه في الواقع، مما يؤدي إلى المبالغة في تقدير التهديد وربما إنتاج ائتمانات حفظ مشوبة بالخلل. والسؤال الجوهري هو ما إذا كانت هذه الموجات السريعة والكثيفة من الإزالة في المزارع الفردية هي مجرد حالة استثنائية أم أنها علامة على أزمة أكبر وأكثر عدوانية.
ولحل هذه المعضلة، وجه الباحثان ويليام شومان وإيونوتس-ألكساندرو تشوكيرلي اهتمامهما إلى الآليات المحددة لكيفية عمل أصحاب المزارع في هذا المشهد. فقد ركزا على اثنتين وعشرين مزرعة ماشية كبيرة في الجزء الغربي من منطقة تشاكو، وهي منطقة يسمح فيها القانون بتحويل الغابة إلى مراعي، بشرط الاحتفاظ بنسبة مئوية معينة من الأشجار كمحمية قانونية. وباستخدام بيانات الأقمار الصناعية من عام 2012 إلى 2020، تتبعا بدقة مقدار الغابة التي تمت إزالتها من كل ملكية والمدة التي استغرقها ذلك الإزالة. ووجدا أنه عندما تكون المزرعة في خضم مرحلة تطوير نشطة، فإن معدل الإزالة يكون مكثفًا بالفعل. ففي المتوسط، قامت هذه المزارع النشطة بإزالة ما يقرب من خمسة عشر بالمائة من إجمالي مساحة أرضها في كل عام. وهذا الرقم مرتفع، وهو تمامًا حيث توقع النقاد أن يكون. ومع ذلك، اكتشف الباحثان أن هذا الرقم المرتفع ينطبق فقط على الشريحة الصغيرة من المزارع التي تعمل بالفعل على إزالة الأراضي في تلك اللحظة المحددة.
إن المفتاح لفهم الصورة الكاملة يكمن في توقيت ونطاق العملية. فقد كشفت الدراسة أن الإزالة ليست قطرات بطيئة وثابتة عبر المنطقة بأكملها؛ بل هي حدث يتركز في البداية. فعندما يبدأ صاحب المزرعة في تحويل الغابة إلى مرعى، فإنه عادة ما يقوم بالعمل الشاق في العام أو العامين الأولين، حيث يزيل مساحة ضخمة من الأرض قبل الانتقال إلى غيرها أو التوقف مؤقتًا. وبحلول العام الثالث، تتباطأ الوتيرة بشكل كبير. وهذا يعني أنه بينما قد تكون مزرعة معينة بمثابة "آلة إزالة" لفترة قصيرة، فإن معظم المزارع الأخرى في المنطقة تظل خاملة، تنتظر دورها، أو أنها انتهت بالفعل. وقد حسب الباحثون أنه في أي عام معين، تكون حوالي اثني عشر بالمائ^ة فقط من المزارع المؤهلة في مرحلة الإزالة النشطة وعالية السرعة هذه. أما الثمانية والثمانون بالمائ}^ة الأخرى فلا تقوم بإزالة أي شيء مهم في ذلك العام.
ومن خلال الجمع بين هاتين الحقيقتين — السرعة العالية للمزارع النشطة والعدد القليل للمزارع النشطة في أي وقت مضى — تمكن الباحثون من سد الفجوة بين المستوى المحلي والإقليمي. فقد أظهروا أنه إذا أخذت معدل الإزالة السريع للمزارع النشطة وضربته في الكسر الصغير من المزارع النشطة فعليًا، فإن النتيجة تطابق تمامًا المعدل الأبطأ والثابت الذي تظهره الخرائط الإقليمية. إن الأرقام المرتفعة للمزارع الفردية حقيقية، لكنها تخف بفضل حقيقة أن معظم المزارع لا تقوم بالإزالة في أي عام معين. ويشير هذا الاكتشاف إلى أن القلق بشأن المبالغة في تقدير التهديد بناءً على بيانات قطع الأراضي هو أمر لا أساس له من الصحة، بشرط أن تأخذ الرياضيات في الاعالح حساب عدد المزارع التي تعمل في آن واحد.
كما نظرت الدراسة إلى مجموعة أكبر مكونة من ستة وأربعين مزرعة للتأكد من أن نتائجهم الأولية لم تكن مجرد صدفة ناتجة عن عينة صغيرة. وحتى مع هذه المجموعة الموسعة، ظل النمط ثابتًا: كانت الإزالة الأكثر عدوانية تحدث في موجات قصيرة، وظل متوسط معدل الإزالة عبر جميع المزارع، عند تعديله وفقًا لعدد المزارع النشطة، متسقًا مع السجلات الإقليمية. وأشار الباحثون إلى أن رقم خمسة عشر بالمائ}^ة لكل عام هو متوسط للممتلكات الأكثر نشاطًا، وأن المعدل الفعلي يختلف بشكل كبير اعتمادًا على حجم المزرعة ورأس المال المتاح للمالك. فالمزارع الأصغر والأكثر تمويلًا تميل إلى الإزالة بشكل أسرع، بينما المزارع الأكبر أو الأقل تمويلًا تسير بشكل أبطأ.
في نهاية المطاف، يقدم هذا العمل تفسيرًا واضحًا لسبب اختلاف الأرقام اعتمادًا على كيفية حسابها. إن إزالة الغابات السريعة التي تُرى في المزارع الفردية هي ظاهرة حقيقية مدفوعة بالحاجة الاقتصادية لإنشاء المراعي بسرعة، ولكن هذا لا يعني أن الغابة بأكملها تتلاشى بنفس هذه السرعة. إن فقدان الغابة هو سلسلة من الأحداث المكثفة وقصيرة العمر المنتشرة عبر مشهد واسع، وليس موجة مستمرة ومنتظمة. ومن خلال التوفيق بين الكثافة المحلية والمتوسط الإقليمي، تقدم الدراسة طريقة لبناء خطوط أساس حفظ أكثر دقة تعكس الديناميكيات الحقيقية لتشاكو دون المبالغة في المخاطر أو التقليل من سرعة التغيير حيث يحدث فعليًا.
ملخص تقني: معدلات مخاطر إزالة الغابات على مستوى القطع الأرضية في منطقة تشاكو الباراغوايانية
بيان المشكلة تتناول الورقة البحثية تناقضاً ملحوظاً بين معدلات مخاطر إزالة الغابات على مستوى القطعة الأرضية (parcel-level) المستمدة بموجب نموذج المعيار المعتمد للتحقق (VCS) - الوحدة VMD0006 (BL-PL)، وبين الإحصائيات الإقليمية لإزالة الغابات في منطقة تشاكو الباراغوايانية. يجادل النقاد بأن خطوط الأساس لبرامج REDD+ التي تُبنى باستخدام قطع أرضية نموذجية قد تبالغ في تقدير خسارة الغابات (الإفراط في منح الائتمان للأساس) لأن معدلات المخاطر السنوية المرصودة في القطع النشطة (على سبيل المثال، ~15% سنوياً) تبدو أعلى بكثير من المتوسطات الإقليمية الممتدة لعقود (على سبيل المثال، ~2-4% سنوياً). وتكمن المشكلة الجوهرية في عدم تطابق المقام: حيث تُحسب معدلات القطعة الأرضية على إجمالي مساحة القطعة خلال مراحل التطوير النشطة، بينما تعمل المعدلات الإقليمية على تخفيف عملية الإزالة عبر كامل المساحة الشاسعة للعقارات النشطة وغير النشطة. تسعى الدراسة إلى تحديد ما إذا كانت هذه المعدلات غير متسقة هيكلياً أم أنه يمكن حل هذا التفاوت من خلال عملية قياس (scaling) صريحة.
المنهجية يستخدم المؤلفون تحليلاً تجريبياً متعدد المقاييس يجمع بين الاستشعار عن بُعد، والبيانات العقارية (الكاتاستر)، والقياس الإحصائي:
مصدر البيانات: تتكون مجموعة البيانات الأولية من 22 قطعة أرض نموذجية عبر تسعة مواقع للمزارع في المنطقة الغربية من تشاكو الباراغوايانية (باستثناء محمية اليونسكو للمحيط الحيوي الشمالية)، والتي تمت مراقبتها من عام 2012 إلى 2020. وتُستخدم عينة ثانوية موسعة مكونة من 46 قطعة أرض (تتضمن 24 قطعة أرض عقارية إضافية) لأغراض التحقق.
اشتقاق المعدلات: تُحسب معدلات مخاطر القطعة الأرضية السنوية (Dparcel) باستخدام صيغة VMD0006 BL-PL: Dpn=D%pn/Yrspn، حيث أن D% هي إجمالي إزالة الغابات كنسبة مئوية من إجمالي مساحة القطعة، و $Yrs$ هو عدد السنوات التي شهدت إزالة فعلية أكبر من الصفر.
هوية القياس (Scaling Identity): تطور الورقة هوية محاسبية حتمية للتوفيق بين مقاييس القطعة الأرضية والمقياس الإقليمي: rregion,conv=(f/g)×Dparcel حيث أن:
rregion,conv هو معدل إزالة الغابات السنوي الإقليمي بالنسبة للغابات القابلة للتحويل قانوناً.
g هو كسر الغابات القابلة للتحويل قانوناً. بالنسبة لحالة مرجع "تشاكو فيفو" (Chaco Vivo)، فإن g=0.575؛ ولتطبيق سيناريوهات القياس الإجمالي الإقليمي، يتم تطبيق عامل قدره 0.70 على إجمالي الغابات الإقليمية لاشتقاق المساحة المركزية القابلة للتحويل (Aconv).
f هو كسر القطعة النشطة (حصة القطع التي تقوم بعمليات إزالة نشطة في سنة معينة).
Dparcel هو معدل مخاطر القطعة الأرضية.
التحليل الإحصائي: تستخدم الدراسة تحليل المخاطر في الوقت المنفصل لنمذجة أنماط الإزالة الزمنية، ومحاكاة مونت كارلو (10,000 سحب) لنشر عدم اليقين في g و f و Dparcel، وإعادة أخذ العينات بطريقة "بوتستراب" (bootstrap) لتوليد فترات الثقة لمتوسط معدل المخاطر.
التحقق: تتم مقارنة حصة النشاط الضمنية (f) المستمدة من هوية القياس مقابل متوسط حصة النشطة لمدة تسع سنوات من مجموعة الـ 46 قطعة أرض الموسعة.
النتائج الرئيسية
معدلات مستوى القطعة: بلغ متوسط معدل مخاطر القطعة الأرضية السنوي للقطع الـ 22 النموذجية 14.99% سنوياً (بفاصل ثقة بوتستراب 95%: من 9.86% إلى 20.68%). التوزيع ملتوي جهة اليمين، مع تباين كبير (الانحراف المعياري = 13.15).
الديناميكيات الزمنية: عملية الإزالة تتركز بقوة في المراحل الأولى (front-loaded). يُظهر تحليل المخاطر في الوقت المنفصل أن متوسط المخاطر ينخفض من 0.28 في السنة الأولى إلى 0.05 بحلول السنة الثالثة (انخفاض بنسبة 81%، p=0.006). وتكمل القطعة الأرضية الوسيطة 89% من إزالتها المرصودة خلال سنتين من بدء العملية.
التوفيق عبر القياس:
باستخدام هوية القياس مع معامل g الخاص بمرجع "تشاكو فيفو" البالغ 0.575 (أو عامل التجميع الإقليمي 0.70 المطبق على مخزون الغابات الإقليمي لاشتقاق المساحة المركزية القابلة للتحويل)، فإن حصة النشاط الضمنية (f) المطلوبة لمطابقة معدلات إزالة الغابات التاريخية (160,000–300,000 هكتار/سنة) تقع بين 0.08 و 0.18.
ينتج عن تحديد الهوية العكسي عبر محاكاة مونت كارلو وسيط حصة f الضمنية البالغ 0.12 (بفترة 95%: 0.07–0.19).
يتوافق هذا مع حصة النشاط المتوسطة المقاسة مباشرة لمدة تسع سنوات والبالغة 0.121 (عند حد إزالة سنوي > 5%) من مجموعة الـ 46 قطعة أرض الموسعة.
التحويل الملتزم به: عند تصفية مجموعة الـ 46 قطعة أرض الموسعة بحثاً عن "التحويل الملتزم به" (إزالة تراكمية ≥ 20%)، فإن المتوسط المشترك يبلغ 15.78%/سنة، وهو ما يشمل تقدير الـ 14.99%/سنة الأصلي.
الإسقاط الإقليمي: تشير تجربة عدم اليقين الاستشرافي إلى أن متوسط إزالة الغابات الضمني يبلغ حوالي 238,000 هكتار/سنة، حيث وقعت 58.4% من عمليات المحاكاة ضمن النطاق التاريخي (160,000–300,000 هكتار/سنة).
المساهمات الرئيسية
التوفيق الرسمي: تقدم الورقة هوية قياس حتمية تربط رياضياً بين شدة مخاطر مستوى القطعة الأرضية وحجم الإزالة الإقليمي، مما يثبت أن التناقض الظاهري هو وظيفة لتعريف المقام (المساحة الكلية مقابل الغابة القابلة للتحويل) وخلط النشاط (القطع النشطة مقابل غير النشطة)، وليس تضخيماً للمخاطر.
التحقق التجريبي من f: تنقل الدراسة "حصة القطعة النشطة" (f) من كونها متبقياً نظرياً إلى كمية مقاسة تجريبياً، حيث تُظهر أن حوالي 12% من المساحة القابلة للتحويل تكون في مرحلة إزالة نشطة في أي سنة معينة، وهو ما يكفي للتوفيق بين معدلات القطعة العالية والإجماليات الإقليمية.
التوصيف الزمني: توضح الدراسة الطبيعة "المتركزة في المقدمة" لتطوير مزارع تشاكو، حيث تُظهر أن معدلات الإزالة السنوية العالية مدفوعة بنبضات إزالة قصيرة ومكثفة (1-2 سنة) بدلاً من الإزالة المستمرة طويلة الأمد.
الأهمية والادعاءات يزعم المؤلفون أن هذه الدراسة تحل الجدل حول "الإفراط في منح الائتمان للأساس" من خلال إظهار أن معدلات مخاطر مستوى القطعة الأرضية المستمدة من منطق BL-PL لا تعني بالضرورة تضخيماً لمخاطر إزالة الغابات الإقليمية. وتؤكد الورقة ما يلي:
أن معدل 14.99%/سنة هو شدة مشروطة خاصة بمراحل التطوير النشطة، ولا يمكن مقارنته مباشرة بالمعدلات الإقليمية دون إجراء عملية القياس (scaling).
عندما يتم حساب كسر الغابات القابلة للتحويل (g) وحصة القطعة النشطة (f) بشكل صريح، فإن معدلات مستوى القطعة الأرضية تكون متسقة هيكلياً مع سجلات إزالة الغابات الإقليمية الممتدة لعقود.
الإطار العملي قابل للتكرار باستخدام بيانات الاستشعار عن بُعد المتاحة للجمهور، والحدود العقارية، وخطط استخدام الأراضي المعتمدة.
لا تدعي الدراسة أنها تعداد وطني ممثل؛ فعينتها صغيرة (N=22) وغير عشوائية، وتركز على الحدود النشطة. ومع ذلك، فإن الاتساق الهيكلي بين المعدلات المشتقة من النماذج والبيانات الإقليمية يشير إلى أن المنهجية سليمة لبناء خطوط الأساس تحت معيار VMD0006.
تخلص الورقة إلى أن التفاوت بين معدلات القطعة الأرضية والمعدلات الإقليمية هو مجرد أثر ناتج عن تجميع البيانات (bookkeeping artifact) وليس دليلاً على فشل منهجي في خطوط أساس REDD+، شريطة تطبيق هوية القياس بشكل صحيح.