A Field-Calibrated Composite Indicator of Urban Vegetation Vulnerability to Compound Heat–Drought Stress: Spatially Blocked Validation and an Exploratory Ecological Application to Cardiovascular Hospitalizations in Isfahan, Iran
تقدم هذه الدراسة مؤشر DUVRI-Eco، وهو مؤشر مركب مُعاير ميدانياً يدمج الحرارة، والجفاف، وحالة الغطاء النباتي، وبنية المناظر الطبيعية لتقييم هشاشة الغطاء النباتي الحضري في أصفهان بإيران، والذي تم التحقق من صحته مقابل نفوق الأشجار ووُجد أنه مرتبط إيجابياً بمعدلات الاستشفاء من أمراض القلب والأوعية الدموية في تحليل بيئي استكشافي.
غالبًا ما تُوصف المدن بأنها غابات خرسانية، ولكن بالنسبة للأشجار التي تصطف في شوارعها، فإن الواقع هو صراع مستمر ضد تهديد مزدوج: الحرارة الحارقة ونقص المياه. وفي المناطق الجافة وشبه القاحلة، حيث يندر هطول الأمطار وترتفع درجات الحرارة، تواجه النباتات الحضرية ضغطًا فريدًا من نوعه. فعندما يتزامن جفاف ممتد مع موجة حر شديدة، لا يكون الإجهاد الواقع على الشجرة مجرد مجموع الاثنين، بل هو قوة تراكمية يمكن أن تسبب فشلًا هيدروليكيًا، مما يؤدي أساسًا إلى تجفيف الشجرة من الداخل. وهذا التراجع في المساحات الخضراء يهم صحة الإنسان بعمق؛ فالأشجار تبرد الهواء، وتفلتر الملوثات، وتوفر حاجزًا ضد الحرارة ذاتها التي تزيد من سوء حالات القلب. ومع ذلك، لا يزال من غير الواضح في العديد من المدن كيف تجتمع هذه الضغوط البيئية معًا لتقتل النباتات، أو كيف قد يترجم فقدان المساحات الخضاء إلى المزيد من الزيارات للمستشفيات بسبب مشاكل القلب.
في مدينة أصفهان بإيران، سعى فريق من الباحثين لفك رموز هذه الشبكة المعقدة. أرادوا تجاوز النظر إلى الحرارة والجفاف كقضايا منفصلة، أو الاعتماد على صور الأقمار الصناعية التي تظهر فقط مدى اخضرار مكان ما في لحظة واحدة. وبدلاً من ذلك، سعوا لبناء أداة عملية واحدة يمكنها قياس مدى تعرض الأشجار الحضرية لهذا الإجهاد المشترك. كان هدفهم إنشاء مؤشر يمكن التحقق منه مقابل أدلة من أرض الواقع: حالات موت الأشجار الفعلية المرصودة على الأرض. ومن خلال القيام بذلك، أملوا أن يروا ما إذا كانت المناطق التي تعاني فيها الأشجار كانت هي أيضًا المناطق التي يعاني فيها الناس أكثر من أزمات القلب والأوعية الدموية.
قام الباحثون بتطوير نظام قياس جديد أطلقوا عليه اسم DUVRI-Eco. ولبناء هذا النظام، جمعوا أربعة معلومات متميزة لكل من مناطق أصفهان الخمس عشرة التابعة للبلدية. أولاً، نظروا في تاريخ الحرارة الشديدة، وتتبعوا أعلى درجات الحرارة المسجلة كل عام. ثانيًا، قاسوا استمرار الجفاف من خلال عد أقصى عدد من الأيام الجافة المتتالية. ثالثًا، قيموا الحالة الصحية الراهنة للنباتات باستخدام صور الأقمار الصناعية لمعرفة مدى اخضرار النباتات. وأخيرًا، فحصوا بنية المناظر الطبيعية، وتحديدًا كيفية ترتيب المساحات الخضراء — سواء كانت موجودة كمساحات خضراء كبيرة ومتصلة أو كبقع صغيرة ومتفرقة من الأشجار. ودمجوا هذه العوامل الأربعة في درجة واحدة، لكنهم لم يخمنوا فحسب مقدار وزن كل عامل؛ بل قاموا بمعايرة النظام باستخدام بيانات من أربعة وثلاثين موقعًا محددًا في جميع أنحاء المدينة حيث قاموا فعليًا بعد الأشجار التي ماتت.
كانت عملية المعايرة هذه هي جوهر عملهم. فمن خلال مقارنة عواملهم البيئية الأربعة مع العدد الفعلي للأشجار الميتة في هذه المواقع، تمكنوا من تحديد أي الضغوطات كانت الأكثر ارتباطًا بالنفوق. وكشف التحليل أن طول فترة الجفاف كان العامل الأكثر تأثيرًا، يليه عن كثب شدة الحرارة، وحالة النباتات، وترتيب البقع الخضراء. ووجد الباحثون أن مؤشرهم الجديد كان جيدًا جدًا في التنبؤ بالمناطق الأكثر عرضة للخطر. وعندما اختبروه عبر استبعاد مناطق كاملة لمعرفة ما إذا كان لا يزال بإمكانه تحديد المناطق الضعيفة، أظهرت الأداة أداءً جيدًا، حيث ميزت بشكل صحيح بين المناطق عالية المخاطر والمنخفضة المخاطر في معظم الحالات. وقد أشار ذلك إلى أن الجمع بين الحرارة والجفاف وبنية المناظر الطبيعية قدم صورة واضحة للأماكن التي يرجح أن تفشل فيها الأشجار.
وبعد أن أسسوا طريقة موثوقة لقياس إجهاد النباتات، اتخذ الفريق خطوة إضافية لمعرفة ما إذا كان لهذا الضعف البيئي تكلفة بشرية. فقد نظروا في السجلات المستشفيات الخاصة بحالات القلب والأوعية الدموية، مثل النوبات القلبية والسكتات الدماغية، عبر نفس المناطق الخمس عشرة على مدار فترة أربع سنوات. وحددوا ارتباطًا بيئيًا وصفيًا: فالمناطق ذات الدرجات الأعلى في مؤشر ضعف النباتات تميل إلى امتلاك معدلات أعلى من حالات دخول المستشفى المسجلة لمشاكل القلب. وتحديدًا، مقابل كل زيادة طفيفة في درجة الضعف، ارتفع معدل دخول المستشفى بنحو خمسة بالمائة. وقد ظل هذا الارتباط قائمًا حتى بعد أن أخذ الباحثون في الاعسبار مستويات تلوث الهواء، المعروفة بتأثيرها على صحة القلب. ومع ذلك، يؤكد الباحثون أن هذه النتيجة هي استكشافية وتوليدية للفرضيات؛ فهي تصف نمطًا عبر الأحياء، لكنها لا يمكن أن تدعم تقييم المخاطر على المستوى الفردي أو تثبت أن فقدان المساحات الخضراء تسبب مباشرة في زيادة زيارات المستشفى.
وعلى الرغم من هذه القيود، يقدم العمل طريقة جديدة لمخططي المدن ومسؤولي الصحة العامة للتفكير في المرونة الحضرية. فهو يتجاوز مجرد الخرائط البسيطة للمساحات الخضراء إلى فهم أكثر دقة لكيفية صمود تلك المساحات تحت وطأة الإجهاد المناخي. ومن خلال تحديد المناطق الأكثر عرضة للآثار المشتركة للحرارة والجفاف، يمكن للمدن تحديد أولويات الاستثمار في الري، أو استبدال الأشجار، أو استراتيجيات التبريد. كما تسلط الدراسة الضوء على رابط محتمل بين صحة الغابة الحضرية وصحة سكانها، مما يشير إلى أن حماية الأشجار من الإجهاد المناخي ليست مجرد هدف بيئي، بل هي ضرورة للصحة العامة. ففي أصفهان، وربما في مدن أخرى شبه قاحلة، قد يكون الحفاظ على حياة الأشجار واحدًا من أكثر الطرق فعالية للحفاظ على نبض المدينة قويًا.
ملخص تقني: مؤشر مركب مُعاير ميدانياً لهشاشة الغطاء النباتي الحضري تجاه الإجهاد المركب الناتج عن الحرارة والجفاف
بيان المشكلة تواجه النباتات الحضرية في المدن القاحلة وشبه القاحلة إجهاداً مركباً من الحرارة الشديدة والجفاف المطول، مما يؤدي إلى تدهور الغطاء النباتي الذي يؤثر على خدمات النظم البيئية وربما على صحة الإنسان. إن المؤشرات الحالية للاستشعار عن بُعد (مثل VCI، وVHI، وSPEI) عادة ما ترصد بعداً واحداً فقط من هذا الإجهاد، أو تفشل في دمج التعرض الهيدرو-مناخي، وحالة النباتات، وهيكل المناظر الطبيعية، والنفوق المرصود ميدانياً ضمن إطار عمل واحد مُعاير. علاوة على ذلك، غالباً ما تعاني طرق التحقق القياسية من تسرب البيانات المكانية، حيث تظهر الملاحظات المتقاربة جغرافياً في كل من مجموعتي التدريب والاختبار، مما يؤدي إلى تضخيم تقديرات الأداء. وهناك فجوة محددة في فهم التأثير المشترك لهذه العوامل على صحة القلب والأوعية الدموية في البيئات الحضرية شبه القاحلة مثل أصفهان، إيران.
المنهجية طورت الدراسة مؤشر DUVRI-Eco، وهو مؤشر مركب قائم على المفهوم (construct-first) مصمم لقياس هشاشة النباتات الحضرية. مرت المنهجية بأربع مراحل:
تكامل البيانات والمعالجة المسبقة:
البيانات الهيدرو-مناخية: استُخدمت سجلات مناخية متجانسة لمدة 72 عاماً (1951–2022) من المحطة الرئيسية في أصفهان والمحطات المجاورة. شملت المتغيرات الرئيسية درجة الحرارة القصوى السنوية (TXx) والعدد الأقصى للأيام الجافة المتتالية (CDD). تم حساب مؤشر شذوذ الرطوبة المركب (CMAI) ولكن تم استبعاده من المؤشر النهائي بسبب التكرار العالي مع CDD (بمعامل ارتباط ρ=0.69).
بيانات النباتات والمناظر الطبيعية: استُخدمت صور أقمار صناعية متعددة العقود (1991–2022). كانت ملاحظات Landsat المباشرة متاحة لسنوات محددة، بينما تمت إعادة بناء السنوات الأخرى باستخدام نموذج لوحة (panel model) يتضمن درجة الحرارة، وهطول الأمطار، وNDVI المتأخر. تم قياس حالة النباتات عبر مؤشر الاختلاف الطبيعي للغطاء النباتي (NDVI). كما تم قياس هيكل المناظر الطبيعية باستخدام مؤشر أكبر رقعة (LPI)، والذي أُطلق عليه اسم تغطية الرقعة المهيمنة (DPC) في هذه الدراسة، وتم حسابه عبر برنامج FRAGSTATS v4.2 على قناع نباتي (NDVI>0.20).
البيانات الصحية: جُمعت سجلات الاستشفاء لأمراض القلب والأوعية الدموية (2019–2022) من المستشفيات التابعة للجامعات في أصفهان على مستوى المناطق (n=15 منطقة).
بناء المؤشر ومعايرته:
يدمج المؤشر أربعة مكونات: TXx، وCDD، وNDVI، وDPC.
تم تطبيع (normalization) المكونات إلى مقيء [0,1] باستخدام تطبيع (min-max) قوي (المئويات من 1 إلى 99) وعكس ترميز المتغيرات الواقية (NDVI، وDPC) بحيث تشير القيم الأعلى إلى هشاشة أكبر.
استُخدم نموذج أقصى احتمالية مقيد لمعايرة أوزان المكونات مقابل نفوق الأشجار المرصود ميدانياً في 34 موقعاً متداخلة ضمن 15 منطقة بلدية. تم تحديد النموذج على مستوى المنطقة، حيث لُخص النفوق كنسبة مجمعة.
كانت الأوزان مقيدة بأن تكون غير سالبة ومجموعها يساوي الواحد الصحيح. عاملت المعايرة الأوزان كمعلمات تنظيم نطاقية (domain-specific regularization parameters) بدلاً من كونها تأثيرات سببية قابلة للتقدير على مستوى السكان.
استراتيجية التحقق:
الحظر المكاني (Spatial Blocking): لتجنب تسرب البيانات، استخدمت الدراسة التحقق المتقاطع بالحظر على مستوى المناطق.
التحقق على مستوى الموقع: ترك موقع واحد خارج الاختبار (34 طية)، مع الإقرار بأن المواقع داخل المنطقة الواحدة تشترك في المتنبئات على مستوى المنطقة.
التحقق على مستوى المنطقة: ترك منطقة واحدة خارج الاختبار (15 طية)، حيث يتم تدريب النموذج على 14 منطقة واختباره على المنطقة المستبعدة. هذا هو المقياس الأساسي للتمييز الجغرافي الداخلي.
تشخيصات قابلية التحديد (Identifiability Diagnostics): استُخدمت أسطح الاحتمالية القصوى (profile-likelihood surfaces) و10,000 متجه وزن تم أخذ عينات منها من المنطقة المدعومة بالاحتمالية لتقييم استقرار درجة المؤشر مقابل دقة الأوزان الفردية.
التطبيق الصحي الاستكشافي:
ربللت دراسة لوحية إيكولوجية (201 9–2022) بين درجات DUVRI-Eco على مستوى المنطقة ومعدلات الاستشفاء لأمراض القلب والأوعية الدموية.
استخدم النموذج توزيع بواسون (Poisson) مع رابط لوغاريتمي، متضمناً تأثيرات ثابتة للشهر والسنة، وتقاطعاً عشوائياً للمنطقة، وPM2.5 كمتغير مربك محتمل تم قياسه.
أقر التحليل صراحةً بأن العدد الفعال لوحدات التعرض المستقلة هو 15 (المناطق)، وليس 720 ملاحظة (منطقة-شهر).
النتائج الرئيسية
أوزان المعايرة: خصص حل أقصى احتمالية المقيد أعلى وزن لاستمرار الجفاف (CDD: 0.31)، يليه الحرارة الشديدة (TXx: 0.28)، ثم حالة النباتات (NDVI: 0.26)، وتغطية الرقعة المهيمنة (DPC: 0.15).
قابلية تحديد الوزن: كشفت تشخيصات الاحتمالية القصوى عن ضعف قابلية التحديد للأوزان الفردية، خاصة لـ CDD وTXx (أسطح احتمالية مسطحة). ومع ذلك، كانت درجة المؤشر المركب والتصنيف المكاني قويين: بلغ متوسط ارتباط الرتب عبر 10,000 متجه وزن محتمل 0.94، وحافظت 13 منطقة من أصل 15 على تصنيف مستقر بالنسبة لنقطة التحول (τ≈0.65).
الأداء التنبئي:
AUC لمستوى الموقع (ترك موقع واحد خارج الاختبار): 0.81 (فاصل ثقة 95%: 0.73–0.88).
AUC لمستوى المنطقة (حظر المنطقة - ترك منطقة واحدة خارج الاختبار): 0.79 (فاصل ثقة 95%: 0.73–0.85). هذا هو التقدير الأساسي للتمييز الجغرافي الداخلي.
أظهر تحليل الاستئصال (Ablation analysis) أن إزالة DPC أدت إلى تغيير طفيف في AUC (ΔAUC=−0.02)، مما يؤكد دوره كمكون معرف للمفهوم وليس كمتنبئ مستقل أساسي.
الارتباط الصحي: ارتبطت زيادة قدرها 0.10 وحدة في DUVRI-Eco بزيادة بنسبة 5.2% في معدل الاستشفاء المسجل لأمراض القلب والأوعية الدموية (نسبة معدل الحدوث [IRR] = 1.052؛ فاصل ثقة 95%: 1.018–1.087). ظل هذا الارتباط متسقاً اتجاهياً عبر تحليلات الحساسية، بما في ذلك التعديل لـ PM2.5 وتحريك الأوزان.
الأنماط المكانية: تركزت الهشاشة العالية في المناطق الوسطى الشرقية (7، 10، 11، 12، 15)، بينما أظهرت المناطق الغربية هشاشة أقل. لوحظ وجود ارتباط مكاني إيجابي بين DUVRI-Eco ومعدلات الاستشفاء (Global Moran's I=0.31).
الأهمية والادعاءات تضع الورقة DUVRI-Eco كـ أداة فحص مُعايرة ميدانياً ومحددة بنطاق الدراسة بدلاً من كونها نموذجاً تنبؤياً عالمياً. وتتمثل مساهماتها الرئيسية في:
بنية قائمة على المفهوم: على عكس النماذج التي تعطي الأولوية للتنبؤ والتي تسعى لتحسين الأداء التدريجي، تم بناء DUVRI-Eco لتمثيل أربعة أبعاد نظرية، تمت معايرتها مقابل نفوق النباتات لترسيخ المفهوم.
التحقق المكاني الصارم: من خلال استخدام التحقق المتقاطع بالحظر على مستوى المناطق، تقدم الدراسة تقييماً أكثر واقعية للتعميم الجغرافي الداخلي، متجنبة المبالغة في مقاييس الأداء الشائعة في التقسيمات غير المحظورة مكانياً.
التمييز بين استقرار الدرجة واستقرار الوزن: توضح الدراسة أنه بينما كانت الأوزان الفردية للمكونات ضعيفة التحديد في عينة صغيرة (15 منطقة)، فإن الدرجة المركبة والتصنيفات المكانية الناتجة كانت مستقرة. وهذا يسمح للمؤشر بالعمل كإطار فحص حتى عندما يكون تحديد الأوزان الدقيقة غير مؤكد.
الارتباط الصحي المولد للفرضيات: يُقدم الارتباط بين هشاشة النباتات وحالات الاستشفاء لأمراض القلب والأوعية الدموية كإيجاد إيكولوجي استكشافي. وتذكر المؤلفات صراحة أن الدراسة لا يمكنها دعم الاستدلال السببي أو الفردي بسبب محدودية وحدات التعرض المستقلة (15 منطقة) والطبيعة الثابتة زمنياً للتعرض في نافذة التحليل الصحي.
القيود والحدود يؤكد المؤلفون أن DUVRI-Eco ليس إطاراً شاملاً لتقييم الهشاشة (حيث يفتقر لمقاييس القدرة التكيفية) وأن أوزانه ليست ثوابت قابلة للنقل. ويعد التحقق الخارجي مع إعادة المعايرة المحلية شرطاً مسبقاً لتطبيق المؤشر أو عتباته على مدن أخرى. إن النتائج الصحية هي حصراً نتائج مولدة للفرضيات، مقيدة بالارتباك الإيكولوجي وعدم القدرة على فصل خصائص المنطقة عن تأثير التعرض. كما أن نقطة التحول المرشحة (τ≈0.65) هي سمة وصفية لعينة معايرة أصفهان، وليست عتبة بيولوجية مثبتة أو حدّاً للإدارة.