Global non-cropland soil organic carbon loss risk under multi-level warming scenarios
تستخدم هذه الدراسة إطاراً مترابطاً بين العمليات والبيانات للكشف عن أن مخزونات الكربون العضوي في التربة غير الزراعية عالمياً تواجه خسائر صافية متزايدة في ظل سيناريوهات الاحترار متعددة المستويات، حيث تحدث أكبر الانخفاضات في النظم البيئية الحرجية والمناطق ذات خطوط العرض العليا مدفوعة بالاختلال بين التحلل المتسارع بفعل درجة الحرارة ومدخلات الكربون.
تحت أقدامنا يكمن مستودع شاسع وغير مرئي، يحمل من الكربون أكثر مما يحتويه الغلاف الجوي وجميع غابات العالم مجتمعة. هذا المستودع هو كربون التربة العضوي، وهو خليط معقد من المواد النباتية المتحللة، والجذور، والحياة المجهرية التي تراكمت على مدى آلاف السنين. إنها طريقة الأرض لتخزين الطاقة، واستقرارها أمر بالغ الأهمية للمناخ. فعندما يبقى هذا الكربون محبوساً في الأرض، فإنه يساعد في الحفاظ على برودة الكوكب. ولكن عندما يتحلل ويتسرب إلى الهواء على شكل ثاني أكسيد الكربون، فإنه يعمل كبطانية، تحبس الحرارة وتسرع من الاحتباس الحراري. لطالما عرف العلماء أن ارتفاع درجات الحرارة يمكن أن يسرع من تحلل هذه المادة العضوية، لكن الصورة الكاملة لكيفية تفاعل أجزاء مختلفة من العالم مع هذا الأمر ظلت ضبابية. فقد ركزت معظم الدراسات السابقة على الأراضي الزراعية، حيث تتدخل الأيدي البشرية باستمرار في حرث التربة وتسميدها، أو على تجارب صغيرة محكومة يصعب توسيع نطاقها لتصل إلى حجم قارة. أما الأراضي البرية الشاسعة — الغابات، والمراعي، والأحراش التي تغطي معظم الكوكب — فقد كانت أقل فهماً في هذا السياق.
وضع فريق من الباحثين نصب أعينهم رسم خريطة لمستقبل هذه التربة البرية في ظل عالم يزداد دفئاً. أرادوا معرفة مقدار الكربون الذي قد يُفقد بالضبط إذا ارتفعت حرارة الكوكب، وأين ستكون تلك الخسائر في أشد حالاتها. ولتحقيق ذلك، بنوا نموذجاً رقمياً يجمع بين طريقتين قويتين في التفكير حول الأرض. أولاً، استخدموا نموذجاً قائماً على العمليات، والذي يعمل مثل حاسبة بيولوجية، تحاكي كيفية استهلاك الميكروبات للمادة العضية وتحللها بناءً على درجة الحرارة والرطوبة. ثم دمجوا ذلك مع أداة تعلم آلي، وهي نوع من البرامج الحاسوبية التي تتعلم الأنماط من كميات هائلة من البيانات الواقعية. ومن خلال تغذية هذا النظام الهجين ببيانات من ما يقرب من 50,000 نقطة لعينات التربة عبر العالم، تمكنوا من محاكاة ما سيحدث لكربون التربة إذا ارتفعت درجات الحرارة بمقادير مختلفة، تتراوح من زيادة طفيفة إلى زيادة شديدة. وقد استبعدوا تحديداً الأراضي الزراعية للتركيز على المناظر الطبيعية غير الزراعية التي تهيمن على سطح الأرض.
ترسم النتائج صورة واضحة ومثيرة للقلق. ففي ظل كل سيناريو من سيناريوهات الاحترار الذي اختبروه، فقدت التربة في هذه النظم البيئية البرية كربونها. لم يكن هناك أي سيناريو اكتسبت فيه التربة كربوناً كافياً لتعويض الخسائر الناجمة عن الحرارة. ومع اشتداد الاحترار، نمت الخسائر في عمقها وسرعتها. وفي السيناريو الأكثر تطرفاً الذي صمّموه، حيث ترتفع درجات الحرارة بشكل كبير، انخفض المخزون العالمي للكربون في الطبقة العليا من التربة بنحو 18 بيتاغرام. ولوضع ذلك في المنظور، فإن هذا يعادل تقريباً إجمالي كمية الكربون التي أطلقت الناتجة عن جميع الأنشطة البشرية على مدى عدة سنوات، ولكن هذه المرة تأتي من الأرض نفسها. ولم تكن الخسائر موزعة بالتساوي؛ بل تركزت في غابات نصف الكرة الشمالي، لا سيما في المساحات الخشبية الشاسعة في شمال أمريكا الشمالية وأوراسيا، وكذلك في المناطق المرتفعة مثل هضبة التبت. وهذه المناطق، التي تخزن حالياً كميات هائلة من الكربون، تبدو الأكثر حساسية للحرارة المتزايدة.
كما كشفت الدراسة أن ليست كل المناظر الطبيعية تتفاعل بالطريقة ذاتها. فقد عانت الغابات من أكبر خسائر الكربون، تليها المراعي ثم الأحراش. وهذا يبدو منطقياً عند النظر في ميكانيكا التربة: فالغابات تمتلك طبقات غنية وعميقة من المواد العضوية التي تكون نشطة للغاية عند تسخينها، مما يؤدي إلى إطلاق سريع للكربون. وفي المقابل، تم تحديد القوى الدافعة وراء هذه التغييرات في كمية نمو النبات ومتوسط درجة الحرارة. فعندما تنمو النباتات أكثر، تضع مزيداً من الكربون في التربة، ولكن عندما يسخن الهواء، تعمل الميكروبات التي تحلل هذا الكربون بشكل أسرع، وتطلقه مرة أخرى في الهواء. وفي عمليات المحاكاة، تجاوزت سرعة التحلل قدرة النباتات على تعويض الكربون المفقود. ووجد الباحثون أن خصائص التربة، مثل قدرة معادن الطين على التمسك بالكربون، لعبت دوراً، لكنها لم تستطع وقف الاتجاه العام للخسارة في عالم دافئ.
وبينما أظهرت عمليات المحاكاة نمطاً ثابتاً من الفقدان، أشار الباحثون إلى أن عدم اليقين في توقعاتهم كان أعلى في الأجزاء الأكثر برودة وبعداً في العالم، حيث تندر البيانات. وهذا يشير إلى أنه بينما يبدو الاتجاه العام واضحاً، فإن الكمية الدقيقة للكربون المفقود في مناطق محددة عالية خطوط العرض قد تكون أكبر مما هو مقدر. لم تدّعِ الدراسة أنها حلت لغز كربون التربة، لكنها قدمت خريطة أكثر دقة للمخاطر. لقد أثبتت أن التوازن الطبيعي لتربة العالم البرية هش. ومع ارتفاع حرارة الكوكب، من المرجح أن تصبح هذه الخزانات الشاسعة مصادر للكربون بدلاً من أن تكون بالوعات له، مما يخلق حلقة مفرغة حيث يسبب الاحترار مزيداً من إطلاق الكربون، والذي بدوره يسبب مزيداً من الاحترار. وتعد هذه النتائج بمثابة تذكير صارخ بأن استقرار المناخ العالمي لا يعتمد فقط على ما نطلقه، بل على كيفية استجابة أنظمة الأرض نفسها للحرارة التي أطلقناها بالفعل.
ملخص تقني: مخاطر فقدان الكربون العضوي في التربة عالمياً في الأراضي غير الزراعية تحت سيناريوهات الاحترار متعددة المستويات
بيان المشكلة يشكل الكربون العضوي في التربة (SOC) أكبر مستودع للكربون البري، ومع ذلك، لا تزال استجابة هذا الكربون للاحترار العالمي غير مفهومة بشكل كافٍ بسبب التفاعلات المعقدة بين المناخ، والنظم البيئية، وخصائص التربة. وبينما استخدمت الدراسات السابقة التجارب الميدانية والنماذج القائمة على العمليات، فإن هذه النهج تواجه قيوداً في التعميم المكاني والقدرة على التوسع، لا سيية بالنسبة للنظم البيئية غير الزراعية (الغابات، والمراعي، والأحراش) التي تهيمن على السطح البري. غالباً ما تركز الأبحاث الحالية على رسم خرائط ثابتة للكربون العضوي في التربة أو على سياقات زراعية محددة، مما يترك فجوة في فهم كيفية استجابة مخزونات الكربون العضوي في الأراضي غير الزراعية العالمية لتغيرات شدة الاحترار عبر أنواع الغطاء الأرضي المختلفة. هناك حاجة إلى إطار عمل يدمج الفهم الميكانيكي للعمليات مع قدرات التنبؤ المكاني لتقييم ديناميكيات الكربون العضوي في التربة تحت سيناريوهات الاحترار متعددة التدرجات.
المنهجية طورت هذه الدراسة إطار عمل "دمج العمليات والبيانات" المقترن الذي يجمع بين نموذج RothC للعمليات وخوارزمية الغابة العشوائية (Random Forest) للتعلم الآلي. وقد مرت المنهجية بالخطوات التالية:
جمع البيانات ومعالجتها مسبقاً: استخدمت الدراسة 49,244 نقطة عينة تربة عالمية من خدمة معلومات التربة العالمية (WoSIS). تم توحيد البيانات لعمق يتراوح بين 0-20 سم باستخدام دالات "spline" متساوية المساحة. تم حصر نطاق الدراسة في المناطق العالمية غير الزراعية، مع استبعاد الأراضي الزراعية، والأراضي الرطبة، والأرض الجرداء، والمسطحات المائية صراحةً لعزل الاستجابات الناتجة عن درجة الحرارة عن المتغيرات المربكة المتعلقة بالإدارة أو الهيدرولوجيا.
تصميم السيناريوهات ومحاكاة العمليات: تم وضع سبعة سيناريوهات لدرجة الحرارة: سيناريو أساسي (0% احترار) وستة سيناريوهات احترار (T5–T30). تم حساب مقادير الاحترار كزيادات بنسبة 5% إلى 30% من الفرق بين درجة الحرارة الشهرية القصوى والمتوسطة. تم تشغيل نموذج RothC عند كل نقطة عينة لمدة عام واحد تحت هذه السيناريوهات لمحاكاة دوران الكربون العضوي في التربة. تم تهيئة النموذج للوصول إلى حالة الاستقرار على مدى 1,000 عام، ثم تم تشغيله للأمام باستخدام البيانات المناخية التاريخية (1981–2014) مع تطبيق اضطرابات درجات الحرارة. تم تعديل مدخلات الكربون بناءً على بيانات الإنتاجية الأولية الصافية (NPP).
النمذجة المكانية: تم حساب الفرق في الكربون العضوي في التربة بين كل سيناريو احترار والأساس (ΔSOC). ولمعالجة التوزيعات ذات الذيول الثقيلة والقيم المتطرفة، خضع المتغير المستهدف لتقليم المئينات (quantile truncation). بعد ذلك، تم تدريب نموذج الغابة العشوائية باستخدام 15 متغيرًا بيئيًا (مناخي، تربة، تضاريس، وعوامل نباتية) للتنبؤ بالتوزيع المكاني لـ ΔSOC عالمياً.
التحليل وعدم اليقين: تم تقييم أداء النموذج باستخدام R2، وRMSE، وMAE. كما تم استخدام تحليل SHAP (تفسيرات شابلي المضافة) لتفسير أهمية الميزات والتأثيرات الهامشية. وتم تقييم عدم اليقين باستخدام غابات الانحدار المئيني (QRF) لإنشاء فترات التنبؤ.
المساهمات الرئيسية
إطار دمج العمليات والبيانات: أرست الدراسة إطاراً جديداً يربط نموذج عملية RothC مع الغابة العشوائية، مما يوازن بين القيود الميكانيكية لدورة الكربون وقدرات التوسع المكاني للتعلم الآلي.
تقييم الاحترار متعدد التدرجات: على عكس الدراسات الثنائية أو أحادية السيناريو، غطت هذه الدراسة بشكل شامل استجابات الكربون العضوي عبر طيف مستمر من شدة الاحترار (5% إلى 30% من نطاق درجة الحرارة) للأراضي غير الزراعية العالمية.
التركيز المكاني المحدد على الأراضي غير الزراعية: استهدفت الدراسة تحديداً النظم البيئية البرية غير الزراعية التي غالباً ما يتم تجاهلها، مما وفر تقييماً دقيقاً لمخاطر فقدان الكربون العضوي في التربة في الغابات والمراعي والأحراش.
النتائج الرئيسية
صافي الخسارة العالمية: تحت جميع سيناريوهات الاحترار، أظهر الكربون العضفي في التربة في الأراضي العالمية غير الزراعية في طبقة 0–20 سم صافي خسارة. ومع زيادة شدة الاحترار من T5 إلى T30، انخفضت مخزونات الكربون العضفي في التربة بمقدار 7.58 إلى 17.60 بيتا جول (Pg)، حيث تراوحت قيم متوسط ΔSOC بين −40.23 و −92.73 جم كربون لكل متر مربع.
حساسية الغطاء الأرضي: شهدت النظم البيئية للغابات أكبر الخسائر المطلقة (تصل إلى 11.88 بيتا جول تحت T30)، تليها المراعي والأحراش. وأظهرت الغابات حساسية أعلى تجاه اضطرابات الاحترار مقارنة بأنواع الغطاء الأرضي الأخرى.
التباين المكاني: تركزت استجابات الكربون العضفي السالبة في خطوط العرض المتوسطة والعالية والمناطق المرتفعة، بما في ذلك شمال أمريكا الشمالية، وشمال أوراسيا، وهضبة تشينغهاي-التبت. وأظهرت المناطق ذات المناخ القاري أكبر الخسائر.
العوامل الدافعة: حدد تحليل SHAP أن الإنتاجية الأولية الصافية (NPP) ومتوسط درجة الحرارة السنوية (AMT) هما المحركان الرئيسيان. ويشير هذا إلى أن استجابات الكربون العضفي في التربة تنظمها التوازنات بين مدخلات النظام البيئي (NPP) وعمليات التحلل المدفوعة بدرجة الحرارة. كما لعبت خصائص التربة المتعلقة باستقرار الكربون (مثل السعة التبادلية الكاتيونية CEC) دوراً هاماً أيضاً.
أداء النموذج: أظهر نموذج الغابة العشوائية قدرة تنبؤية مستقرة عبر السيناريوهات، مع قيم R2 تتراوح بين 0.50 و 0.54. وأظهرت النتائج منحنى استجابة أكثر سلاسة واستمرارية لدرجة الحرارة مقارنة بالتباين الملاحظ في تجارب الاحترار المنعزلة.
الأهمية والادعاءات يدعي المؤلفون أن هذه الدراسة تثبت أن استجابات الكربون العضفي في التربة الناجمة عن الاحترار تنظمها التوازنات الديناميكية بين مدخلات الكربون، وخسائر التحلل، وعمليات استقرار كربون التربة. ويُقدم الإطار المقترن كتقدم منهجي يحسن دمج القيود الميكانيكية مع التنبؤ المكاني، مما يوفر رؤى جديدة حول استجابات الكربون العضفي في التربة للاحترار العالمي.
تخلص الدراسة إلى أن الإطار المقترح يدعم إدارة الكربون الإقليمية في ظل الظروف المناخية المتغيرة من خلال الكشف عن اتجاه متسارع في فقدان الكربون العضفي في التربة تحت سيناريوهات الاحترار القوي. وهي تسلط الضوء على أن النظم البيئية غير الزراعية، التي تشكل غالبية السطح البري، معرضة لفقدان الكربون مع اشتداد الاحترار، لا سيما في مناطق خطوط العرض العالية والمناطق القارية. ويضع المؤلفون عملهم كخطوة نحو الانتقال من رسم خرائط الكربون العضفي في التربة الساكنة إلى تقييمات الاستجابة القائمة على السيناريوهات، مما يوفر مرجعاً للدراسات المستقبلية حول التغذية المرتدة بين الكربون العضفي في التربة والمناخ.