Comparing geospatial interpolation techniques for estimating ambient PM2.5 mass concentrations from sugarcane burning in South Florida
باستخدام شبكة كثيفة من المستشعرات منخفضة التكلفة في منطقة إيفرجليدز الزراعية بجنوب فلوريدا، قيمت هذه الدراسة تقنيات الاستيفاء الجيومكاني لتقدير الجسيمات الدقيقة (PM2.5) الناتجة عن حرق قصب السكر، ووجدت أنه بينما حققت طريقة الوزن بالمسافة العكسية أفضل أداء بشكل فردي، فإن النهج الديناميكي الذي يختار النموذج الأمثل يومياً وفر تقديرات التعرض الأكثر قوة لأبحاث الصحة المستقبلية.
المؤلفون الأصليون:Anne Mielnik, Olivia Sablan, Jeffrey R. Pierce, Lilah M. Besser, Juyoung Park, Diana Mitsova, Diane Cook, Sherry WeMott, Christine L. Williams, Lisa A.K. Wiese, Sheryl Magzamen
المؤلفون الأصليون: Anne Mielnik, Olivia Sablan, Jeffrey R. Pierce, Lilah M. Besser, Juyoung Park, Diana Mitsova, Diane Cook, Sherry WeMott, Christine L. Williams, Lisa A.K. Wiese, Sheryl Magzamen
في الزوايا الريفية لجنوب فلوريدا، حيث تهيمن حقول قصب السكر الشاسعة على المشهد الطبيعي، يشكل طقس موسمي طبيعة الهواء الذي يتنفسه الناس. ففي كل عام، يقوم المزارعون بحرق أكثر من 400,000 فدان من هذه الحقول لإزالة الأوراق القديمة وتجهيز المحصول للحصاد. هذه الممارسة، رغم كفاءتها في الزراعة، تطلق سحباً كثيفة من الدخان تحتوي على جزيئات دقيقة. هذه الجزيئات متناهية الصغر، بحيث يمكنها التسلل إلى أعماق الرئتين وحتى دخول مجرى الدم، مما يشكل مخاطر على صحة الإنسان. ولعقود من الزمن، كافح العلماء لقياس كمية هذا الدخان التي تصل بالضبط إلى المنازل والأحياء الفردية. تكمن المشكلة في أن أجهزة مراقبة جودة الهواء الرسمية التي تستخدمها الحكومة قليلة ومتباعدة في هذه المناطق الريفية، مما يترك فجوات كبيرة في البيانات. ومن دون شبكة كثيفة من أجهزة الاستشعار، يصعب معرفة ما إذا كان الشخص الذي يعيش على بعد أميال قليلة من حقل محترق يتنفس هواءً أنقى من شخص آخر يعيش بجوار الحقل تماماً. ولفهم التأثير الصحي الحقيقي لهذا الدخان، يحتاج الباحثون إلى وسيلة لملء تلك القطع المفقودة من اللغز، لإنشاء صورة كاملة لجودة الهواء عبر المنطقة بأكملها.
لقد شرع فريق من الباحثين في حل هذه المشكلة من خلال نشر شبكة من أجهزة استشعار الهواء منخفضة التكلفة عبر منطقة إيفرجليدز الزراعية. ومن يونيو 2024 إلى ديسمبر 2025، وضعوا 44 جهازاً صغيراً في منازل المتطوعين وفي مواقع ثابتة، حيث سجلت هذه الأجهزة قراءات هوائية كل ساعة. كما جمعوا السجلات الرسمية لمواعكن وأماكن حرق المزارعين لحقولهم، بالإضافة إلى صور الأقمار الصناعية التي تتبعت حركة سحب الدخان. لم يكن الهدف مجرد قياس جودة الهواء حيث توجد أجهزة الاستشعار، بل استخدام هذه القياسات لتقدير حالة جودة الهواء في كل مكان آخر في المنطقة، بما في ذلك أفنية منازل الأشخاص الذين لا يملكون أجهزة استشعار. ولتحقيق ذلك، اختبر الفريق عدة طرق رياضية مختلفة، تُعرف بالاستيفاء المكاني (geospatial interpolation)، والتي تعمل بمثابة جسر يربط بين نقاط البيانات المعروفة والنق points غير المعروفة. وقارنوا بين طريقة تعتمد ببساطة على إعطاء وزن كبير للقياسات القريبة، وبين نماذج إحصائية أكثر تعقيداً تحاول التنبؤ بكيفية انتشار التلوث عبر المناظر الطبيعية.
وجد الباحثون أن جودة الهواء في هذه المنطقة تعتمد بشكل كبير على وجود الدخان. ففي الأيام التي حدث فيها حرق لقصب السكر، ارتفع تركيز الجزيئات الدقيقة في الهواء بشكل ملحوف، بزيادة متوسطها 1.65 ميكروغرام لكل متر مكعب مقارنة بالأيام التي لم يحدث فيها حرق. ولم يكن هذا الارتفاع محسوساً بالتساوي في كل مكان؛ فقد كان الدخان أكثر تركيزاً في المناطق الداخلية، بالقرب من الحقول حيث يتم الحرق وحول بحيرة أوكيتشوبي. وفي المقابل، ظل الهواء على طول الساحل أكثر نقاءً. كشف هذا النمط عن تفاوت صارخ: فالمجتمعات التي تعيش بالقرب من حقول الحرق، والتي تضم غالباً سكانًا من ذوي الدخل المنخفض وغير البيض، كانت معرضة لأعلى مستويات التلوث، بينما تنفس أولئك الذين يعيشون بالقرب من المحيط هواءً أنقى بكثير. وقد أكدت الدراسة أن حرق قصب السكر يعد مصدراً رئيسياً وموضعياً للتلوث، مما يخلق لوحة فسيفسائية من جودة الهواء عبر المنطقة.
وعندما تعلق الأمر بتحديد أفضل طريقة لتقدير مستويات التلوث في المناطق غير المراقبة، اكتشف الباحثون أن البساطة غالباً ما كانت تعمل بشكل أفضل من التعقيد. فالطريقة التي حققت أداءً أكثر اتساقاً كانت تقنية تسمى "الترجيح بالمسافة العكسية" (inverse distance weighting). يفترض هذا النهج أن المواقع القريبة من جهاز الاستشعار تكون أكثر تشابهاً مع قراءة ذلك الجهاز من المواقع البعيدة عنه. وبينما أنتجت النماذج الأكثر تقدماً نتائج أفضل أحياناً في أيام معينة، إلا أنها أدت أيضاً في كثير من الأحيان إلى أرقام مستحيلة، مثل مستويات تلوث سالبة أو تركيزات عالية جداً تخالف الواقع، خاصة عندما تكون البيانات شحيحة. ومع ذلك، نادراً ما ارتكبت الطريقة الأبسط هذه الأخطاء، ووفرت تقديرات موثوقة في حوالي 97 بالمائة من الأيام التي تمت دراستها. وخلص الفريق إلى أن الطريقة الأكثر متانة لتقدير التعرض لدراسات الصحة المستقبلية هي استخدام نظام مرن وشبه آلي يختار الطريقة الأفضل لكل يوم محدد بناءً على كمية البيانات المتاحة ومدى ملاءمة النماذج للظروف.
يوفر هذا النهج طريقة جديدة لدراسة الآثار الصحية للحرق الزراعي دون الحاجة إلى أجهزة مراقبة باهظة الثمن من الدرجة الحكومية في كل حي. ومن خلال الجمع بين أجهزة الاستشعار منخفضة التكلفة وبيانات الأقمار الصناعية والخيارات الإحصائية الذكية، أنشأ الباحثون إطاراً يمكن استخدامه في مناطق ريفية أخرى حيث تكون مصادر التلوث عرضية ومراقبة البيانات فيها شحيحة. لم تثبت الدراسة أن دخان قصب السكر يسبب أمراضاً معينة لدى هذا السكان، لكنها نجحت في توفير بيانات التعرض التفصيلية والمحددة مكانياً التي يحتاجها علماء الصحة لطرح تلك الأسئلة في المستقبل. وتسلط النتائج الضوء على أنه بينما يعد الدخان جزءاً طبيعياً من الدورة الزراعية، فإن تأثيره ليس موزعاً بالتساوي، مما يترك بعض المجتمعات تتحمل عبئاً أثقل من التلوث مقارنة بغيرها.
ملخص تقني: مقارنة تقنيات الاستيفاء المكاني لتقدير الجسيمات الدقيقة (PM2.5) المحيطة الناتجة عن حرق قصب السكر في جنوب فلوريدا
بيان المشكلة يُعد الحرق الزراعي، وتحديداً حرق قصب السكر قبل الحصاد في منطقة إيفرجليد الزراعية (EAA) بجنوب فلوريدا، مصدراً هاماً للمواد الجسيمية الدقيقة (PM2.5). وبينما ربطت دراسات سابقة بين هذا التلوث والنتائج الصحية الضارة (مثل الربو، وأمراض الجهاز التنفسي، والوفيات) على مستوى المقاطعة أو المنطقة الإحصائية، إلا أن هناك نقصاً حاداً في البيانات المتعلقة بالتعرض المحدد للموقع (سكني أو مهني). وينبع هذا الفجوة من ندرة محطات الرصد التنظيمية لجودة الهواء في المناطق الريفية، ومن الطبيعة العارضة وقصيرة الأمد لسحب الدخان، والتي تفشل نماذج التعرض طويلة الأمد القياسية في رصدها. وبناءً على ذلك، فإن الدراسات الصحية على مستوى الأفراد التي تبحث في الآثار قصيرة المدى للحرق الزراعي تعوقها عدم القدرة على تقدير تركيزات الجسيمات الدقيقة (PM2.5) اليومية المحددة للموقع بدقة.
المنهجية لمعالجة هذه الفجوات، نشر الباحثون شبكة كثيفة مكونة من 44 مستشعراً منخفض التكلفة من نوع (PurpleAir) في مقاطعات غليدز، وهيندري، وأوكيتشوبي، وبالم بيتش، مدعومة بمراقبين اثنين من نظام جودة الهواء (AQS) التابع لوكالة حماية البيئة الأمريكية (EPA) باستخدام الطريقة المرجعية الفيدرالية (FRM). استمرت عملية جمع البيانات من 8 يونيو 2024 إلى 1 ديسمبر 2025، لتغطي مواسم الحرق وغير الحرق.
استخدمت الدراسة إطار عمل لمعالجة البيانات والنمذجة متعدد الخطوات:
معالجة البيانات: تم تصحيح بيانات المستشعرات الخام (بفاصل 10 دقائق) باستخدام عامل معايرة "باركجون" (Barkjohn) لتتوافق مع معايير وكالة حماية البيئة للجسيمات الدقيقة (PM2.5). تم تجميع البيانات في متوسطات ساعية و24 ساعة، مع إزالة القيم المتطرفة باستخدام طريقة المدى الربيعي (IQR).
تحديد تأثير الدخان: تم تحديد الأيام "المتأثرة بالدخان" باستخدام أربعة معايير: (1) PM2.5 ≥ 15 ميكروجرام/م³ لمدة ≥ ساعة واحدة؛ (2) PM2.5 ≥ 20 ميكروجرام/م³ لمدة ≥ ساعة واحدة؛ (3) PM2.5 ≥ 30 ميكروجرام/م³ لمدة ≥ ساعة واحدة؛ و(4) نهج هجين يجمع بين عتبة 20 ميكروجرام/م³ وبيانات سحب الدخان من نظام رسم خرائط المخاطر (HMS) التابع لـ NOAA. كما استُخدم معيار ثانٍ أكثر شمولاً (C2) يعتمد على سجلات تصاريح الحرق المفتوح (OBA) من وزارة الزراعة والمستهلكات في فلوريدا.
الاستيفاء المكاني: قام الباحثون بملاءمة 11 نموذجاً مكانياً مختلفاً لبيانات متوسط الـ 24 ساعة المحولة لوغاريتمياً على شبكة بمساحة 5 كم. شملت النماذج:
طريقة الوزن العكسي للمسافة (IDW).
الكريجينج العادي (Ordinary Kriging - OK) مع خمسة أنواع من ملاءمة التباين (Bessel، Exponential، Gaussian، Linear، Nugget).
الكريجينج العام (Universal Kriging - UK) مع خمسة أنواع من ملاءمة التباين.
تقييم النموذج: تم تقييم الأداء باستخدام التحقق المتبادل بترك عنصر واحد خارج المجموعة (LOOCV)، وحساب معامل التحديد (R2)، وجذر متوسط مربع الخطأ (RMSE)، ومتوسط الخطأ المطلق (MAE).
نهج التجميع (Ensemble): تم اختبار نهج "مجمع" شبه مؤتمت، يقوم باختيار النموذج "الأفضل" لكل يوم بناءً على عدد القياسات، وملاءمة التباين، ومقاييس التحقق المتبادل (LOOCV).
التحليل الإحصائي: استُخدمت اختبارات "ويلكوكسون" للرتب الموقعة (Wilcoxon signed-rank tests) لمقارنة تركيزات PM2.5 المقدرة بين الأيام المتأثرة بالدخان والأيام غير المتأثرة عبر منطقة الدراسة.
النتائج الرئيسية
جمع البيانات: أنتجت الدراسة 12,814 قياساً لمتوسط 24 ساعة من الجسيمات الدقيقة (PM2.5). كان وسيط عدد المستشعرات النشطة يومياً هو 24 (المدى: 2–34).
تأثير الدخان: كانت تركيزات PM2.5 المقدرة أعلى بكثير في الأيام المتأثرة بالدخان مقارنة بالأيام غير المتأثرة في جميع النماذج (p<0.001). في المتوسط، أدى وجود دخان حرق قصب السكر إلى زيادة تركيزات PM2.5 لمدة 24 ساعة بمقدار 1.65 ميكروجرام/م³ (الخطأ المعياري = 0.13 ميكروجرام/م³).
التوزيع المكاني: كانت تركيزات PM2.5 هي الأعلى في المناطق الداخلية، بالقرب من منطقة إيفرجليد الزراعية (EAA) وبحيرة أوكيتشوبي، وهي الأدنى على طول الساحل. يسلط هذا الضوء على التفاوتات في جودة الهواء، حيث تتركز الفئات ذات الدخل المنخفض وغير البيضاء بشكل أكبر في مناطق الحرق الداخلية.
أداء النماذج:
IDW مقابل الكريجينج: أدت نماذج IDW أداءً أفضل بشكل عام، حيث حافظت على أعلى قيم R2 وأقل مقاييس خطأ في معظم الأيام. نجح IDW في تقدير التركيزات في 527 يوماً من أصل 541 يوماً (97.4%). في المقابل، أنتجت نماذج OK وUK مع ملاءمات Bessel وGaussian غالباً قيم غير واقعية (تركيزات سالبة أو قيم متطرفة تتجاوز 1000 ميكروجرام/م³) وفشلت في تقدير البيانات في الأيام التي بها قياسات أقل.
التباين: اختلف أداء النموذج حسب الموسم ووجود الدخان. تراوحت قيم R2 بين 0 و0.99. أظهرت نماذج OK مع ملاءمة "Nugget" قيمة R2=1، ولكن اعتُبر ذلك مضللاً لأن النموذج توقع قيماً متطابقة لجميع المواقع، مما فشل في رصد التباين المكاني.
النهج المجمع: حقق النهج شبه المؤتمت لاختيار أفضل نموذج يومياً أداءً فائقاً، حيث حققت 82.9% من الأيام قيمة R2>0.90.
تقديرات التعرض: عند تطبيقها على 50 موقعاً سكنياً، أظهرت نماذج IDW أكبر تمايز بين أيام الدخان وغير الدخان (متوسط 5.52 ميكروجرام/م³ مقابل 3.31 ميكروجرام/م³)، بينما أظهرت بعض نماذج الكريجينج تمايزاً ضئيلاً، مما يشير إلى أنها قد تعمل على "تنعيم" (smoothing) الذروات العارضة.
الأهمية والادعاءات يزعم المؤلفون أن هذه هي الدراسة الأولى التي تقدر التعرضات المحددة للموقع لـ PM2.5 الناتجة عن الحرق الزراعي في المواقع غير المراقبة في جنوب فلوريدا. تكمن الأهمية الأساسية للدراسة في تطوير إطار عمل لدعم الدراسات الوبائية المستقبلية على مستوى الأفراد، وتحديداً تلك التي تبحث في النتائج الصحية قصيرة المدى مثل الوظائف الإدراكية في مناطق الحرق المتقطع.
المساهمات الرئيسية تشمل:
التحقق المنهجي: إثبات أن IDW هي طريقة قوية وفعالة حاسوبياً للمصادر الملوثة العارضة في المناطق الريفية حيث تكون كثافة البيانات متغيرة، وغالباً ما تتفوق على نماذج الكريجينج الأكثر تعقيداً في هذا السياق المحدد.
تكامل البيانات الهجين: النجاح في دمج بيانات المستشعرات منخفضة التكلفة المعايرة مع اكتشافات سحب الدخان المستمدة من الأقمار الصناعية وتصاريح الحرق التنظيمية للتغلب على قيود الرصد الشحيح في المناطق الزراعية الريفية.
رؤية العدالة الصحية: تقديم دليل على أن التعرض لـ PM2.5 الناتج عن حرق قصب السكر ليس مجرد قضية إقليمية فحسب، بل إنه يؤثر بشكل غير متناسب على المجتمعات ذات الدخل المنخفض وغير البيضاء التي تقطن في المناطق الداخلية بالقرب من EAA، وذلك على عكس مستويات التعرض المنخفضة على طول الساحل.
يظل المؤلفون متواضعين بشأن ادعاءاتهم، حيث أقروا بأن الأخطاء في نماذج التعرض قد تنتقل إلى تقديرات التأثير في الدراسات الصحية المستقبلية. كما أشاروا إلى أن تعيين متوسط واحد لمدة 24 ساعة قد يؤدي إلى التقليل من تقدير الآثار الصحية بسبب عدم رصد التغيرات الزمنية بالساعة. علاوة على ذلك، حذروا من أنه على الرغم من إمكانية تعميم نهجهم على مناطق قصب سكر أخرى (لويزيانا، تكساس) ومناطق الحرق الزراعي، إلا أنه قد يتطلب تعديلاً للمناطق ذات التضاريس المعقدة أو التشتت القوي المدفوع بالرياح.