A lightweight automatic workflow for expanding volcano-tectonic seismic catalogs: application to Villarrica Volcano, Chile
تقدم هذه الورقة سير عمل خفيف الوزن وآلياً بالكامل يوسع بشكل كبير سجلات الزلازل البركانية التكتونية من خلال اكتشاف وتحديد مواقع الأحداث التي فات رصدها سابقاً، كما تم إثباته من خلال تطبيقه الناجح على بركان فيلاريكا في تشيلي، والذي ضاعف السجل الحالي تقريباً مع الحفاظ على دقة تضاهي التحليل اليدوي الذي يقوم به الخبراء.
المؤلفون الأصليون:Sergio Morales, Sarah Jaye Oliva, Andrés I. Ávila
تخيل أن الأرض تشبه منزلاً ضخماً عابساً يئن ويصدر أصوات صرير ويتحرك في أساساته بين الحين والآخر. أحياناً، تكون هذه الأصوات مجرد استقرار للمنزل، ولكن في أحيان أخرى، تكون علامات تحذيرية على أن شيئاً كبيراً على وشك الحدوث، مثل استيقاظ بركان. العلماء الذين يدرسون "أصوات المنزل" هذه يُطلق عليهم علماء الزلازل، وهم يستمعون لأنواع محددة من الارتجاجات تسمى الزلازل "البركانية التكتونية". فكر في هذه الزلازل كأنها أصوات طقطقة حادة ناتجة عن تكسر الصخور في أعماق جدران البركان بينما يشق الماغما (الصخور المنصهرة شديدة الحرارة) طريقه للأعلى. إن رصد هذه الطقطقات أمر بالغ الأهمية لأنها تخبرنا ما إذا كان البركان يتعرض لضغط وقد ينفجر قريباً.
ومع ذلك، فإن الاستماع إلى بركان يشبه محاولة سماع صوت سقوط دبوس واحد في ملعب مزدحم وصاخب. فالإشارات غالباً ما تكون ضئيلة، ويمكن لـ "الضجيج" الناتج عن الرياح أو حركة المرور البعيدة أن يحجبها. تقليدياً، اعتمد العلماء على فرق من الخبراء البشريين للجلوس أمام شاشات الكمبيوتر، والاستماع لساعات من التسجيلات الصوتية للعثور على تلك الطقطقات الضئيلة. لكن البشر يتعبون، ولا يمكنهم الاستماع إلى كل شيء في وقت واحد، وخلال الأوقات المزدحمة عندما يبدأ البركان في "ثرثرته" فعلياً، فإنهم يفتقدون حتماً أهدأ وأصغر الطقطقات. تسأل هذه الورقة سؤالاً بسيطاً: هل يمكننا بناء روبوت يستمع بشكل أفضل من إنسان متعب، دون الحاجة إلى كمبيوتر خارق بحجم مدينة؟
قرر الباحثون وراء هذه الدراسة بناء ذلك بالضبط: سير عمل خفيف الوزن وتلقائي يعمل كـ "محقق" دؤوب وشديد الانتباه لبركان فيلاريكا في تشيلي. لم يكتفوا ببناء روبوت للاستماع فحسب؛ بل بنوا فريقاً كاملاً من المتخصصين الرقميين الذين يعملون معاً. أولاً، استخدموا نموذج "تعلم عميق" (نوع من الذكاء الاصطناعي المدرب على آلاف الأصوات الزلزالية) ليعمل كـ "آذان"، حيث يقوم بمسح البيانات لرصد أضعف موجات P وموجات S (الارتجاج الأول والثاني للزلزال). بعد ذلك، استخدموا خوارزمية تجميع ذكية تسمى MeanShift لتعمل كـ "عقل"، حيث تقوم بتجميع تلك الارتجاجات المتناثرة لمعرفة: "مهلاً، هذه الأصوات الثلاثة حدثت في نفس الوقت تقريباً؛ لا بد أنها من نفس الحدث!". أخيراً، استخدموا أداة تحديد موقع احتمالية لتعمل كـ "رسام خرائط"، حيث تحسب بدقة أين حدثت تلك الطقطقة تحت الأرض.
نتائج هذا العمل الاستقصائي الرقمي مثيرة للإعجاب. فعندما اختبروا نظامهم على بيانات مستمرة لمدة عامين من بركان فيلاريكا (من يناير 2023 إلى ديسمبر 2024)، قارنوا قائمة الروبوت بالقائمة الرسمية التي أعدها الخبراء البشريون. نجح الروبوت في العث de 72.5% من الأحداث التي وجدها البشر بالفعل، مما يثبت أنه لا يقل كفاءة في رصد الأحداث الكبيرة والواضحة. لكن السحر الحقيقي حدث مع الأحداث التي فاتها البشر. فقد اكتشف الروبوت 695 زلزالاً إضافياً لم تكن موجودة في سجل البشر، وهذا يمثل توسعاً بنسبة 91.4% في النشاط الزلزالي المعروف!
لم تكن هذه الزلازل "الإضافية" ضجيجاً عشوائياً؛ بل كانت في معظمها أحداثاً صغيرة وضحلة بالقرب من سطح البركان، بقوة تتراوح بين 1 و2 درجة. هذه هي الطقطقات الصغيرة والمراوغة التي غالباً ما يفتقدها البشر عندما يغرقون في حجم كبير من النشاط. لم يكتفِ الروبوت بالعثور على المزيد فحسب، بل وجد تلك التي كانت تختبئ في وضح النهار. كما كان توقيت هذه الأحداث دقيقاً للغاية، حيث كان الفارق في توقيت الروبوت حوالي 0.44 ثانية فقط عن تقديرات الخبراء البشر — وهو فرق ضئيل جداً لا يغير من حقيقة ما يفعله البركان.
توضح الورقة أن هذا ليس عصا سحرية تحل كل المشكلات. فالنظام مصمم خصيصاً للزلازل "البركانية التكتونية" (نوع تكسر الصخور الحاد) وقد يفتقد أنواعاً أخرى من الأصوات البركانية التي لا تمتلك ارتداد موجة S ثانية واضحة. كما أنه يعتمد على نفس خرائط سرعة الصخور تحت الأرض التي يستخدمها البشر، لذا إذا كانت تلك الخرائط غير مثالية، فستكون تخمينات الروبوت للموقع غير مثالية أيضاً. ومع ذلك، تثبت الدراسة أن هذا النظام الخفيف والتلقائي يمكنه العمل على خوادم الكمبيوتر القياسية الموجودة في معظم المراصد دون الحاجة إلى بطاقات رسوميات باهظة الثمن أو مراقبة بشرية مستمرة. وهذا يشير إلى أنه من خلال ترك المهمة الشاقة للروبوت في عملية الاستماع، يمكن لعلماء البراكين توسيع سجلاتهم بشكل كبير، ورصد الهمسات الهادئة للبركان التي قد يفتقدها البشر، وكل ذلك باستخدام الأجهزة التي يمتلكونها بالفعل.
بيان المشكلة تعد النشاطات الزلزالية البركانية التكتونية (VT) مؤشراً رئيسياً لتغيرات الإجهاد وهجرة السوائل داخل البناء البركاني، مما يجعل الفهرسة الدقيقة أمراً ضرورياً لتقييم المخاطر. ومع ذلك، تواجه المراقبة التشغيلية في براكين مثل "فيلاريكا" (تشيلي) قيوداً مزدوجة متراكمة: ندرة شبكات المراقبة مقترنة بالضجيج الخلفي المرتفع، والطبيعة الكثيفة للعمل اليدوي. غالباً ما تفشل سير العمل اليدوية في اكتشاف أصغر الأحداث الزلزالية (التي تتراوح عادة بين ML 1–2) خلال فترات النشاط الزلزالي المرتفع بسبب تداخل الموجات وضيق وقت المحللين. يؤدي هذا إلى سجلات غير مكتملة قد تحجب التسلسلات التحذيرية وتؤدي إلى تحيز في تقييمات المخاطر. وبينما أظهر التعلم العميق (DL) إمكانات واعدة، فإن النماذج الحالية غالباً ما تتطلب تدريباً خاصاً بالموقع، وتفتقر إلى القابلية للنقل عبر البيئات الجيولوجية المختلفة، ولم يتم التحقق من صحتها بصرامة مقابل سجلات الخبراء المستقلة في بيئات تشغيلية مستمرة.
المنهجية يقترح المؤلفون سير عمل آلياً، نموذجياً، وخفيف الوزن مصمماً للعمل على البنية التحتية القياسية للمرصد دون تدخل يدوي. يدمج النظام ثلاثة مكونات مفتوحة المصدر:
كشف الطور (Phase Detection): يستخدم سير العمل نموذج PhaseNet، وهو نموذج تعلم عميق مدرب لتحديد وصول موجات P و S. ولتعظيم استعادة الأحداث ذات القدر المنخفض، تم ضبط عتبة الكشف عند 0.4 (أقل من القيمة الافتراضية)، مع إسناد الرقابة الصارمة على الجودة إلى المراحل اللاحقة.
ربط الأطوار (Phase Association): بدلاً من القواعد المعقدة على مستوى الشبكة، يستخدم النظام خوارزمية MeanShift clustering (خوارزمية تجميع تعتمد على الوقت) لتجميع وصولات طور P من المحطات القريبة في أحداث مرشحة. وتعتمد عرض نطاق التجميع (clustering bandwidth) على معلمة واحدة قابلة للتفسير الفيزيائي مشتقة من فتحة الشبكة وسرعة موجة P نموذجية. يسمح هذا النهج للنظام بالتكيف مع أي هندسة للشبكة دون الحاجة لإعادة تدريب النموذج.
التحديد الموقعي الاحتمالي (Probabilistic Location): يتم تحديد موقع الحدث باستخدام خوارزمية NonLinLoc، وهي خوارزمية تحديد مواقع زلزالية احتمالية غير خطية. يتبع سير العمل إجراءً متعدد المراحل:
الموقع الأولي: يتم حسابه باستخدام أطوار P من الشبكة القريبة.
البحث عن طور S: بحث مستهدف عن أطوار S في المحطات القريبة بناءً على المركز البؤري الأولي. المحطات البعيدة: يتم دمج أطوار P و S من المحطات البعيدة فقط إذا استوفت معايير اتساق السعة (لضمان انتمائها لنفس الحدث بدلاً من مصدر أقرب).
التحسين التكراري: لا تُقبل الحلول إلا إذا أدت إلى تقليل عدم اليقين في الموقع.
يميز سير العمل بين شبكة "قريبة" (المحطات ضمن 15 كم من القمة) وشبكة "بعيدة" لموازنة قابلية اكتشاف الطور مع دقة العمق. طُبق الخط بأكمله على بيانات مستمرة لمدة عامين (يناير 2023 – ديسمبر 2024) من بركان فيلاريكا، باستخدام شبكة مكونة من 16 محطة.
النتائج الرئيسية تم التحقق من صحة سير العمل الآلي مقابل السجل الذي راجعه الخبراء من مرصد "أنديز ديل سور البركاني" (OVDAS)، مع قصر البحث على الأحداث التي يقل خطأ تحديد موقعها عن 5 كم.
توسيع السجل: استعاد سير العمل 72.5% (551 حدثاً) من 760 حدثاً يدوياً تم تحديد موقعها بدقة. والأهم من ذلك، فقد حدد 695 حدثاً إضافياً لم تكن موجودة في السجل اليدوي، مما يمثل توسعاً بنسبة 91.4% في سجل الأحداث البركانية التكتونية المحددة بدقة.
خصائص الأحداث: تتركز الأحداث المكتشفة حديثاً بالقرب من البناء البركاني عند أعماق ضحلة وضمن نطاق قدر ضيق يتراوح بين ML 1–2. هذه الفئة عادة ما تكون ممثلة بنقص في عمليات المعالجة اليدوية الروتينية.
دقة زمن الأصل: أظهرت أزمنة الأصل الآلية انحيازاً سالباً ثابتاً صغيراً قدره −0.44 ثانية بالنسبة للنماذج اليدوية. ويُعزى ذلك إلى وضع PhaseNet للوصولات بعد بدايات النماذج اليدوية بقليل، وهو اتجاه لا يؤثر على تفسير معدلات النشاط الزلزالي أو الأنماط المكانية.
جودة الموقع: بالنسبة للأحداث المتطابقة، اتفقت المواقع المركزية (epicentral locations) بشكل وثيق بين الطريقتين. أظهرت تقديرات العمق تشتتاً أكبر للأحداث الضحلة (أقل من 5 كم)، وهو ما يتماشى مع القيود المعروفة لدقة العمق في نماذج السرعة القريبة من السطح.
الكفاءة الحسابية: عالج سير العمل بيانات عامين على جهاز افتراضي قياسي (4 نوى معالجة، 16 جيجابايت رام، بدون وحدة معالجة رسومية) بمعدل إنتاجية يبلغ حوالي ساعة واحدة من الحساب لكل يوم من البيانات المستمرة.
الأهمية والادعاءات تضع الورقة البحثية سير العمل هذا كبديل عملي وقابل للتوسع للمراصد البركانية التي تسعى لتوسيع سجلاتها الزلزالية دون زيادة متناسبة في عبء عمل المحللين أو البنية التحتية الحسابية.
القابلية للتشغيل: تم تصميم النظام ليعمل على أجهزة قياسية دون تدخل يدوي، مما يجعله مناسباً للنشر في الوقت الفعلي تقريباً.
القابلية للنقل: لا يتطلب سير العمل سوى مجموعة صغيرة من المعلمات القابلة للتفسير الفيزيائي (مثل قوائم المحطات، عرض نطاق التجميع، وعتبات الجودة) للتكيف مع الشبكات البركانية المختلفة، مما يتجنب الحاجة إلى إعادة تدريب النماذج الخاصة بكل موقع.
الاستعادة المستهدفة: يستهدف النظام تحديداً الفئة "المفقودة" من الأحداث البركانية التكتونية الصغيرة (ML 1–2) التي غالباً ما تُفقد خلال فترات النشاط العالي، مما يحسن دقة رصد تغيرات الإجهاد في القشرة الضحلة.
القيود والنطاق: يشير المؤلفون صراحةً إلى أن سير العمل منحاز ضمنياً لمصادر الفشل الهش ذات بدايات موجة S النبضية (الأحداث البركانية التكتونية) ويقلل من تمثيل الأحداث طويلة الفترة أو الهجينة. علاوة على ذلك، لم يتم تحديد مقدار الاكتمال (magnitude of completeness) رسمياً، ويُقدم التوسع بنسبة 91.4% كحد أدنى للتحسن في حساسية الكشف. كما تقر الدراسة بأن نموذج السرعة أحادي البعد المستخدم هو قيد مشترك مع المعالجة اليدوية، لا سيما بالنسبة للمحطات البعيدة.
في الختام، تثبت الدراسة أن خط معالجة خفيف الوزن ومفتوح المصدر يمكن أن يعزز بشكل كبير اكتشاف وتحديد مواقع النشاط الزلزالي البركاني التكتوني، مما يوفر أداة قوية للمراصد لمراقبة الاضطرابات البركانية بشكل أكثر شمولاً.