A comprehensive, multi-model exploration of flow scenarios for supporting probabilistic hazard assessment of pyroclastic density currents: application at Fuego volcano (Guatemala)
تقدم هذه الدراسة إطار عمل شاملاً متعدد النماذج باستخدام TITAN2D وVolcflow لتقييم مخاطر تيارات الكثافة البركانية في بركان فويغو، كاشفةً أنه بينما تتحكم التضاريس بشكل أساسي في مسارات التدفق، فإن الافتراضات الريولوجية الخاصة بكل نموذج تؤثر بشكل كبير على بصمات الغمر والقدرة على الحركة، مما يؤسس نهجاً قابلاً للنقل لنشر عدم اليقين المصدري والهيكلي في تقييمات المخاطر الاحتمالية.
المؤلفون الأصليون:Elaine T. Spiller, Pablo Tierz, Sylvain Charbonnier, Eric C. P. Breard, Rüdiger Escobar-Wolf, Robert L. Wolpert, Roberto Mérida, Gustavo Chigna, Eliza S. Calder
المؤلفون الأصليون: Elaine T. Spiller, Pablo Tierz, Sylvain Charbonnier, Eric C. P. Breard, Rüdiger Escobar-Wolf, Robert L. Wolpert, Roberto Mérida, Gustavo Chigna, Eliza S. Calder
تشكل البراكين التي لا تنام حقاً تحدياً فريداً للناس الذين يعيشون في ظلها. فخلافاً للبركان الذي يثور مرة واحدة ثم ينام لقرون، تنتج هذه العمالقة المستمرة انفجارات صغيرة متكررة يمكن أن تتحول فجأة إلى أنهار سريعة الحركة من الغاز والرماد والصخور. يسمي العلماء هذه الأنهار "تيارات الكثافة البركانية الفتاتية". وهي مدمرة بكفاءة مرعبة لكل ما يعترض طريقها، وغالباً ما تسافر عشرات الكيلومترات بعيداً عن القمة. وبالنسبة للمجتمعات التي تعيش على منحدرات مثل هذه البراكين، فإن معرفة المسارات التي قد تسلكها هذه التدفقات هي مسألة حياة أو موت. ومع ذلك، فإن التنبؤ بمسارها أمر صعب للغاية لأن الأرض غير مستوية، وتتصرف التدفقات بطرق معقدة، كما أن الظروف التي تسببها يمكن أن تتفاوت بشكل هائل. ولاتخاذ قرارات آمنة بشأن عمليات الإخلاء أو الاستيطان طويل الأمد، يحتاج العلماء إلى ما هو أكثر من مجرد تخمين واحد؛ إنهم بحاجة إلى صورة واضحة لجميع الطرق الممكنة التي يمكن أن يسلكها التدفق، بما في ذلك مدى المسافة التي قد يقطعها وأين قد يفيض خارج قنواته الطبيعية.
لقد تصدى فريق من الباحثين مؤخراً لهذه المشكلة في بركان "فويغو" في غواتيمالا، وهو جبل نشط وشديد الانحدار تسبب في حدوث تدفقات مدمرة في الماضي، بما في ذلك حدث كارثي في عام 2018 أودى بحياة المئات. لم يعتمد الفريق على برنامج حاسوبي واحد للتنبؤ بالمستقبل، بل أجروا آلاف المحاكاة باستخدام نموذجين حاسوبيين مختلفين وموثوقين على نطاق واسع، يقومان بحساب كيفية تحرك هذه التدفقات الساخنة والثقيلة فوق التضاريس الحقيقية. أرادوا معرفة ما إذا كانت النماذج ستتفق على الصورة الكبيرة، حتى لو اختلفت في التفاصيل. ولتحقيق ذلك، أنشأوا مكتبة ضخمة من السيناريوهات؛ حيث قاموا بمحاكاة تدفقات تتراوح من الصغيرة جداً إلى الضخمة، وتستمر من ثوانٍ إلى عدة دقائق، وقاموا بتغيير "الاحتكاك" الداخلي لهذه التدفقات لمعرفة كيف يغير ذلك من حركتها. وقد أجروا هذه المحاكاة عبر ستة أودية رئيسية تقطع منحدرات البركان، منتجين إجمالي 3600 سيناريو تدفق مختلف.
كشفت النتائج عن اتساق مطمئن في الخطوط العريضة. فقد اتفقت كلا النموذجين الحاسوبيين على الحدود الخارجية للمناطق التي يمكن أن تصل إليها التدفقات. والأهم من ذلك، حدد كلاهما نفس النقاط بالضبط على طول الأودية حيث يتغير شكل الأرض بطريقة تجبر التدفق على الفيضان خارج قناته والانتشار عبر الأرض المحيطة. وهذا يشير إلى أن شكل التضاريس نفسه هو العامل الأكثر قوة في التحكم في مسار هذه التيارات الخطيرة. فإذا كان في الأرض انحناء حاد أو استواء مفاجئ، فمن المرجح أن يفيض التدفق، بغض النظر عن النموذج الحاسوبي المستخدم. وهذا يمنح مخططي الطوارئ الثقة في إمكانية تحديد المناطق الأكثر عرضة للخطر بناءً على التضاريس وحدها.
ومع ذلك، قدم النموذجان قصتين مختلفتين تماماً عما يحدث بمجرد فيضان التدفق. فقد أنتج أحد النموذجين تدفقات طويلة وضيقة، تظل متماسكة داخل وديان الأنهار لمعظم رحلتها. بينما أنتج النموذج الآخر تدفقات أقصر ولكنها أعرض بكثير، حيث تنتشر مثل ورقة رقيقة عبر الأرض بمجرد خروجها من المجرى الرئيسي. ويعود هذا الاختلاف إلى كيفية تعامل كل نموذج مع المقاومة التي تبذلها الأرض تجاه الصخور المتحركة. يفترض أحد النموذجين أنه كلما أصبح التدفق أقل سمكاً، فإنه يتباطأ بشكل أسرع، مما يبقيه محصوراً. بينما يفترض الآخر أن التدفق يستمر في الحركة طالما كانت الأرض شديدة الانحدار، حتى لو كانت طبقة الصخور رقيقة جداً. وهذا يعني أنه بينما يتفق النموذجان على أين قد يفيض التدفق، فإنهما يختلفان على مقدار الأرض التي سيغطيها بمجرد حدوث ذلك.
اكتشف الباحثون أيضاً أن حجم التدفق يلعب دوراً بطريقة محددة. فقد وجدوا نقطة تحول عند حجم قدره 10 ملايين متر مكعب تقريباً. فدون هذا الحجم، تميل التدفقات إلى البقاء داخل القنوات. أما فوق هذا الحجم، فقد بدأت التدفقات في أحد النموذجين بالانتشار بسرعة، متحولة إلى طبقات رقيقة يمكنها تغطية مساحات شاسعة من الأراضي بعيداً عن الوادي الرئيسي. يساعد هذا الاكتشاف في تفسير سبب تسبب بعض الثورات في أضرار محلية بينما تغطي أخرى مناطق بأكملة. علاوة على ذلك، استخدم الفريق أدوات إحصائية متقدمة لرسم خريطة توضح بالضبط ما هو مزيج حجم التدفق وسرعته ومدته اللازم للوصول إلى قرى محددة تقع بالقرب من ضفاف الأنهار. ووجدوا أن الوقت الذي يستغرقه التدفق لبدء الحركة يمكن أن يغير النتيجة، حتى بالنسبة للمواقع التي تبعد عشرة كيلومترات عن القمة.
في نهاية المطاف، لا يوفر هذا العمل تنبؤاً واحداً مثالياً بالكارثة القادمة، بل يوفر إطاراً قوياً لفهم نطاق الاحتمالات. فمن خلال تشغيل آلاف السيناريوهات والمقارنة بين طريقتين مختلفتين في فهم فيزياء الحركة، أظهر الباحثون أنه بينما تعتمد تفاصيل سلوك التدفق على الأدوات المستخدمة لقياسه، فإن المخاطر الأساسية تمليها تضاريس الأرض. ويسمح هذا النهج للعلماء بالانتقال من مجرد التخمينات البسيطة نحو فهم احتمالي أكثر اكتمالاً للمخاطر البركانية. وبالنسبة للمجتمعات التي تعيش على منحدرات "فويغو"، يعني هذا أن خرائط السلامة المستقبلية يمكن بناؤها على أساس يقر بوجود عدم اليقين بدلاً من تجاهله، مما يضمن أن القرارات المتعلقة بمكان العيش أو كيفية الإخلاء تستند إلى الطيف الكامل لما قد يفعله الجبل.
ملخص تقني: استكشاف شامل متعدد النماذج لسيناريوهات التدفق لتقييم المخاطر الاحتمالية لتيارات الكثافة البركانية في بركان فوغو
بيان المشكلة تشكل البراكين النشطة باستمرار والتي تولد تيارات الكثافة البركانية (PDCs) مخاطر عالية على السكان المحيطين، كما يتضح من ثوران بركان فوغو في غواتيمالا عام 2018، والذي أدى إلى خسائر كبيرة في الأرواح بسبب تجاوز التيارات لحواف الأودية (overbanking). يتطلب التخفيف الفعال للمخاطر تقيما علميا قويا يأخذ في الاعتبار أوجه عدم اليقين في ظروف المصدر، والتضاريس، والريولوجيا (علم التدفق)، وبنية النموذج. وبينما يعد نمذجة التدفق أداة قياسية، فإن النهج المعتمد على نموذج واحد غالباً ما يفشل في استيعاب النطاق الكامل لعدم اليقين المعرفي المرتبط ببنية النموذج والصياغة الريولوجية. هناك حاجة إلى إطار عمل منهجي لنشر أوجه عدم اليقين هذه في تقييمات المخاطر البركانية الاحتمالية (PVHA) لدعم القرارات المتعلقة بالإخلاء وإعادة التوطين.
المنهجية يقدم المؤلفون استكشافاً منهجياً متعدد النماذجي لسيناريوهات تيارات الكثافة البركانية (PDC) في بركان فوغو باستخدام نموذجين من نماذج التدفق متوسط العمق واسعي الاستخدام: TITAN2D و Volcflow. يستخدم هذان النموذجان صياغات ريولوجية وتنفيذاً عددياً متباينين.
تصميم المحاكاة: تم استخدام تصميم "لاتين-هايبركيوب" (Latin-hypercube) لتوليد إجمالي 3,600 محاكاة (300 لكل نموذج لكل "بارانكا" - أخدود) عبر ستة "بارانكاس" رئيسية تصب من البركان.
معلمات المدخلات: قامت عمليات المحاكاة بتغيير ثلاثة معلمات مدخلات رئيسية: حجم الـ PDC (يرمز له بـ V، ويتراوح من 0.3 إلى 300 مليون متر مكعب)، ومدة تدفق المصدر (d)، والمعلمات الريولوجية الخاصة بكل نموذج (احتكاك كولوم μ لنموذج TITAN2D؛ والإجهاد المثبط الثابت، CRS، لنموذج Volcflow). تم استخدام تصميم ملء الفراغ (space-filling design) على مقياس لوغاريتمي لاستكشاف فضاء السيناريوهات دون افتراض توزيعات احتمالية كثيفة محددة مسبقاً.
التضاريس والمصدر: استخدمت عمليات المحاكاة نموذج سطح رقمي (DSM) بدقة 8 أمتار مشتق من بيانات ALOS AW3D. تم نمذجة بدء الـ PDC كمصادر تدفق تقع عند رؤوس الأودية (barrancas)، على بعد حوالي 0.5 كم من حافة الفوهة، مع ربط الأحجام والمدد عبر علاقة خطية لوغاريتمية-لوغاريتمية لتعكس الارتباطات الملحوظة بين حجم التدفق ومدته.
المحاكاة الإحصائية (Emulation): لفك رموز عتبات الغمر بكفاءة دون الحاجة لتشغيل عشرات الآلاف من عمليات المحاكاة الإضافية، استخدم المؤلفون محاكيات العمليات الغاوسية ذات التصفير (zGP emulators). تم تدريب هذه المحاكيات على مخرجات المحاكاة لنمذجة العلاقة بين معلمات المدخلات والكمية المراد قياسها (أقصى عمق للتدفق) عند نقاط اهتمام محددة (A، B، و C) في "بارانكا ترينيداد". سمح ذلك بالتحديد السريع لـ "خطوط تساوي العمق الصفري" (zero-flow isolines) التي تفصل بين سيناريوهات الغمر وسيناريوهات عدم الغمر.
المساهمات الرئيسية
مقارنة منهجية للنماذج: يقدم هذا العمل المقارنة المنهجية الأكثر شمولاً بين TITAN2D و Volcflow حتى الآن، موضحاً كيف تؤثر صياغاتهما الريولوجية المتباينة (المعتمدة على كولوم مقابل الإجهاد المثبط الثابت) وتنفيذاتهما العددية على الجريان، وتشكيل القنوات، والانتشار الجانبي، وتجاوز الحواف.
إطار عمل قابل للنقل: يضع إطار عمل منهجياً وقابلاً للنقل لتطبيق نماذج التدفق العددية على البراكين النشطة والمتكررة إنتاجاً لتيارات الكثافة البركانية، والمصمم خصيصاً لدعم تقييمات المخاطر البركانية الاحتمالية (PVHA) المستقبلية.
التحكم التضاريسي مقابل التحكم بالنموذج: تقيم الدراسة مدى استقلالية مواقع تجاوز الحواف عن اختيار النموذج، مع التمييز بين الضوابط التضاريسية والسلوكيات المعتمدة على النموذج.
مقياس حركية جديد: قدم المؤلفون مقياس نسبة المساحة الجديد (A100cm:A10cm)، والذي يُعرف بأنه نسبة المساحة المغمورة بتدفقات بعمق ≥ 100 سم إلى المساحة بعمق ≥ 10 سم. يميز هذا المقياس الترقق والانتشار البعيد، معالجاً أوجه القصور في مقاييس الحركية القياسية التقليدية مثل H/L (ارتفاع السقوط/طول الجريان).
تطبيق المحاكي (Emulator): يستخدم العمل محاكيات العمليات الغاوسية لتحديد تركيبات محددة من الحجم والمدة والمعلمات الريولوجية التي تفصل بين سيناريوهات الغمر وعدمه عند مواقع مختارة، مما يسلط الضوء على تأثير مدة المصدر على الغمر لمسافات تصل إلى 10 كم من القمة.
النتائج الرئيسية
الأغلفة الخارجية: ينتج كلا النموذجين أغلفة غمر خارجية متقاربة بشكل عام، ويحددان باستمرار نفس المناطق التضاريسية (مثل انكسارات المنحدرات، وانحناءات الوديان) المعرضة لتجاوز الحواف. يشير هذا إلى أن هندسة الأودية (barranca geometry) تمثل التحكم الأساسي الأول على مسارات التدفق.
الاختلافات السلوكية المنهجية: داخل هذه الأغلفة، يظهر النموذجان سلوكيات متباينة:
TITAN2D: ينتج بصمات أوسع وأقصر عموماً مع انتشار أكثر اتساعاً خارج القنوات. مقاومة كولوم المهيمنة فيه مستقلة عن عمق التدفق، مما يسمح للطبقات الرقيقة بالاستمرار في الحركة إذا تجاوز المنحدر المحلي زاوية الاحتكاك.
Volcflow: ينتج تيارات أضيق وأطول وأكثر تشكلاً في القنوات. إن صياغة الإجهاد المثبط الثابت تجعل الحركية تعتمد على عمق التدفق؛ فمع ترقق التدفقات عند تجاوز الحواف، تزداد المقاومة، مما يؤدي إلى توقفها بسرعة.
حجم الانتقال: يحدد مقياس نسبة المساحة الجديد انتقالاً في سلوك TITAN2D حول 10 مليون متر مكعب (في الأودية المدروسة)، حيث تتوسع التيارات بسرعة كطبقات رقيقة متجاوزة للحواف. دون هذا الحجم، تظل التدفقات شديدة التشكل في القنوات؛ وفوقه، تنتشر بشكل كبير.
رؤى المحاكي (Emulator): تكشف نتائج zGP أن مدة المصدر تؤثر على الغمر لمسافات تصل إلى 10 كم من القمة. علاوة على ذلك، تغير بنية النموذج من مدى تعرض المواقع الواقعة داخل القناة أو خارجها. فعلى سبيل المثال، في TITAN2D، تتعرض المواقع خارج القناة للغمر بأحجام أصغر من المواقع داخل القناة، بينما يظهر Volcflow الترتيب المعاكس.
محدودية المقاييس: تؤكد الدراسة أن مقاييس الحركية القياسية (مثل H/L) توفر فقط وصفاً جزئياً لسلوك الـ PDC ولا يمكنها التمييز بين النقل داخل القنوات والانتشار الواسع للطبقات الرقيقة.
الأهمية والادعاءات يدعي البحث أنه يضع إطار عمل قابلاً للنقل لنشر عدم اليقين في المصدر والنموذج ضمن تقييمات مخاطر الـ PDC الاحتمالية. ومن خلال التعامل مع النموذجين كتمثيل لسلوكيات طرفية متباينة بدلاً من التنافس على نموذج واحد "أفضل"، يوفر المؤلفون طريقة لحصر عدم اليقين في الانتشار الجانبي والجريان.
تُقدم النتائج كخطوة تأسيسية نحو خريطة مخاطر احتمالية كاملة لبركان فوغو. ويؤكد المؤلفون أنه بينما يشكل العمل الحالي استكشافاً مصمماً لفضاء السيناريوهات بدلاً من كونه مجموعة احتمالية مرجحة، فإن كتالوج البصمات الناتجة والمحاكيات المدربة مصممة للسماح بنشر نماذج الاحتمالية اللاحقة لظروف المصدر والمعلمات عبر النظام. ويشدد العمل على ضرأ ضرورة حمل كل من ظروف المصدر (الحجم، المدة) وعدم اليقين الهيكلي للنموذج إلى حسابات المخاطر المستقبلية لتقييم المخاطر بدقة للسكان الذين يعيشون في المناطق البينية وعلى طول حواف التصريف. ويظل المؤلفون متواضعين بشأن القيود، مشيرين إلى أن النماذج لا تحل فيزياء الفشل الأولية، أو التآكل، أو اقتران التيارات المركزة مع السحب المخففة، وأن العمل المستقبلي يجب أن يدمج نتائج هذه المحاكاة مع بيانات السكان وتقدير عدم اليقين لإنتاج خرائط مخاطر تشغيلية.