🔬 physics

Exceedance Probabilities for Large Earthquakes From DIY Local Earthquake Ensemble Nowcasting and Forecasting

تقدم هذه الورقة طريقة "تحويل التنبؤ الآني" لتعديل إحصائيات غوتنبرغ-ريشتر ضمن مجموعات الزلازل الإقليمية لتحسين التنبؤ الآني والتنبؤ المستقبلي لاحتمالات تجاوز الزلازل الكبيرة، مما يثبت اتساقها وإمكانية تطبيقها على منطقة لوس أنجلوس في أعقاب زلزال نورث ريدج عام 1994.

المؤلفون الأصليون: John B Rundle, Ian Baughman, Andrea Donnellan, Lisa Grant Ludwig, Geoffrey Fox, Kazuyoshi Nanjo

نُشر 2026-04-07

📖 4 دقيقة قراءة☕ قراءة في استراحة قهوة

CC BY 4.0

المؤلفون الأصليون: John B Rundle, Ian Baughman, Andrea Donnellan, Lisa Grant Ludwig, Geoffrey Fox, Kazuyoshi Nanjo

البحث الأصلي مرخَّص بموجب CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ هذا شرح مولَّده بالذكاء الاصطناعي للبحث أدناه. لم يكتبه المؤلفون ولم يصادقوا عليه. وللتحقق من الدقة التقنية، يرجى الرجوع إلى البحث الأصلي. اقرأ إخلاء المسؤولية الكامل

تخيل أنك تحاول التنبؤ بموعد هبوب العاصفة الكبرى القادمة على بلدتك. لا يمكنك رؤية العاصفة وهي قادمة، ولكن يمكنك عدّ عدد قطرات المطر الصغيرة التي تسقط كل يوم. هذه الورقة البحثية تتحدث عن طريقة جديدة لاستخدام تلك "قطرات المطر الصغيرة" (الزلازل الصغيرة) لتخمين متى قد تضرب "العاصفة الكبرى" (الزلزال الكبير)، وتحديداً لمنطقة لوس أنجلوس.

إليك تفصيل طريقتهم باستخدام تشبيهات بسيطة:

١. الفكرة الجوهرية: عدّ القطرات

يعتقد العلماء أنه قبل حدوث زلزال هائل، لا تبقى الأرض ساكنة فحسب، بل تصبح "متوترة". ومع ازدياد توترها، تطلق هزات صغيرة (زلازل صغيرة).

التشبيه: فكر في شريط مطاطي يتم شده. قبل أن ينقطع (الزلزال الكبير)، يصدر أصوات "تزييق" و"طقطقة" صغيرة (زلازل صغيرة). وكلما زادت أصوات "التزييق" التي تسمعها، اقتربت المطاط من الانقطاع.
الهدف: بدلاً من التخمين بناءً على التقويم (مثلاً: "مرت ١٠ سنوات، لذا هو موعده الآن")، هم يقومون بعدّ عدد الزلازل الصغيرة منذ آخر زلزال كبير. ويسمون هذا "الزمن الطبيعي".

٢. المشكلة: عدم كفاية البيانات

لتقديم تخمين جيد، تحتاج إلى الكثير من التاريخ. ولكن بالنسبة للوس أنجلوس، لم نشهد زلزالاً ضخماً منذ عام ١٩٩٤ (زلزال نورثريج). وهذا ليس كافياً لبناء نموذج موثوق لتلك المدينة وحدها.

التشبيه: تخيل أنك تحاول تعلم كيفية خبز كعكة مثالية، لكنك لم تخبز سوى كعكة واحدة فقط في حياتك. لا يمكنك أن تكون واثقاً جداً في وصفتك.

٣. الحل: "المجموعة" (دراسة الحي)

لحل مشكلة نقص البيانات، نظر العلماء إلى لوس أنجلوس ثم رسموا مستطيلات أكبر فأكبر حولها، مثل التموجات في بركة ماء.

التشبيه: نظروا إلى "حي لوس أنجلوس"، ثم "منطقة لوس أنجلوس الكبرى"، ثم "جنوب كاليفورنيا"، وهكذا. وافترضوا أن الطريقة التي تتصرف بها الزلازل في هذه المناطق الأكبر هي مشابهة إحصائياً لكيفية تصرفها في الدائرة الصغيرة حول لوس أنجلوس.
الطريقة: عاملوا كل هذه المناطق ذات الأحجام المختلفة كـ "فصل تدريبي". فقد درسوا كيف تؤدي الزلازل الصغيرة إلى زلازل كبيرة في كل هذه المناطق لبناء "وصفة رئيسية" للتنبؤ بالزلزال الكبير القادم في لوس أنجلوس.

٤. "التنبؤ الآني" (تقرير الطقس للوقت الحالي)

لقد طوروا أداة تسمى "التنبؤ الآني" (Nowcast). هذا ليس توقعاً طويل الأمد (مثل توقع لمدة ١٠ سنوات)؛ بل هو "تقرير حالة حالي".

التشبيه: فكر في تطبيق للطقس يخبرك: "من المحتمل بنسبة ٤٠٪ أن تمطر خلال الساعة القادمة".
كيف يعمل: هم يعدّون الزلازل الصغيرة منذ عام ١٩٩٤. ومع زيادة هذا "العد"، يرتفع تقرير "التنبؤ الآني"، مما يخبرنا أن النظام يقترب من لحظة "الانقطاع".

٥. "المرآة السحرية" (تحويل التنبؤ الآني)

هنا يكمن الذكاء. أدرك العلماء أن المناطق الكبيرة (التموجات) قد يكون لها قواعد مختلفة قليلاً عن دائرة لوس أنجلوس الصغيرة. وللتأكد من مطابقة بيانات التدريب للهدف، اخترعوا "تحويل التنبؤ الآني".

التشبيه: تخيل أنك تحاول تعلم قيادة سيارة فيراري، لكنك تملك سيارة تويوتا فقط في مدرسة القيادة الخاصة بك. التويوتا تتعامل مع الطريق بشكل مختلف. "التحويل" يشبه جهاز محاكك سحري يعدل مقود وسرعة التويتوتا بحيث تشعر وكأنك تقود الفيراري تماماً.
النتيجة: قاموا بتعديل البيانات من المناطق الكبيرة لكي تطابق القواعد الإحصائية لدائرة لوس أنجلوس تماماً. ووجدوا أنه حتى مع هذا التعديل، ظلت توقعاتهم ثابتة ومتسقة.

٦. الساعتان: الزمن التقويمي مقابل الزمن الطبيعي

تقارن الورقة بين طريقتين لقياس الزمن:
١. الزمن التقويمي: كم عدد السنوات التي مرت منذ عام ١٩٩٤؟
٢. الزمن الطبيعي: كم عدد الزلازل الصغيرة التي حدثت منذ عام ١٩٩٤؟

النتيجة المفاجئة:

الزمن التقويمي: مع مرور السنوات، تبدو فرصة حدوث زلزال كبير وكأنها تنخفض في البداية، ثم ترتفع ببطء. الأمر يشبه انتظار حافلة تأخرت كثيراً؛ تبدأ في الاعتقاد أنها قد لا تأتي أبداً، لكن في النهاية، أنت تعلم أنها قادمة.
الزمن الطبيعي: مع زيادة عد الزلازل الصغيرة، ترتفع فرصة حدوث زلزال كبير بشكل حاد.
لماذا الاختلاف؟ الأرض تصبح "أكثر صلابة" بمرور الوقت. فمع تعافي القشرة الأرضية وتصلبها بعد زلزال كبير، تتوقف عن إحداث شقوق صغيرة (زلازل صغيرة) لفترة من الوقت. لذا، حتى لو مرت سنوات عديدة، إذا كان "العد" للزلازل الصغيرة منخفضاً، فإن النظام ليس "مشحوناً" بالتوتر بعد. "الزمن الطبيعي" هو مؤشر أفضل لتوتر الأرض من ساعة الحائط.

ملخص النتائج للوس أنجلوس

الوضع الحالي: منذ زلزال نورثريج عام ١٩٩٤، وقع حوالي ٤٤٨ زلزالاً صغيراً في دائرة قطرها ١٢٥ كم حول لوس أنجلوس.
التوقع: بناءً على هذا العد، هناك احتمال يتراوح بين ٢٩٪ إلى ٤٣٪ (اعتماداً على الطريقة المحددة المستخدمة) لحدوث زلزال بقوة ٦.٠ أو أكثر في العام القادم.
الخلاصة: الأرض لا تنتظر مجرد تاريخ على التقويم؛ إنها تنتظر عدداً معيناً من "التزييقات" (الزلازل الصغيرة) قبل أن تكون جاهزة لـ "الانقطاع".

باختصار: تقدم هذه الورقة طريقة أفضل للاستماع إلى "تزييق" الأرض لفهم مدى قربنا من "الانقطاع" الكبير القادم، وذلك باستخدام بيانات من المنطقة بأكملها لتقديم تخمين أكثر ذكاءً للمدينة.

إليك ملخص تقني مفصل للورقة البحثية بعنوان: "احتمالات التجاوز للزلازل الكبيرة: من الآنية (Nowcasting) والإنذار المبكر لمجموعات الزلازل المحلية ذاتية التنفيذ: عبر الحجم، والزمن الطبيعي، والزمن التقويمي."

1. بيان المشكلة

تتناول الورقة البحثية تحدي توقع حدوث الزلازل الكبيرة محلياً ( $M_T \geq 6$ ) ضمن منطقة جغرافية محددة (مُعرفة كدائرة بنصف قطر عشوائي يحيط بنقطة اهتمام، مثل مدينة ما). تكمن المشكلة الجوهرية في تحديد كمية احتمالية وقوع هذه الأحداث ضمن زمن طبيعي محدد (عدد الزلازل الصغيرة) وزمن تقويمي (وقت الانتظار) أعقب آخر حدث كبير.

يهدف المؤلفون إلى التغلب على ندرة البيانات في المناطق المحلية من خلال استخدام نهج المجموعات (Ensemble Approach)، حيث تُستخدم البيانات الإحصائية من مناطق مستطيلة أكبر محيطة لتدريب وتوقع الاحتمالات للدائرة الصغيرة محل الاهتمام. وهناك فرضية حرجة يتم اختبارها وهي ما إذا كانت إحصائيات تكرار الحجم (Gutenberg-Richter) لهذه المناطق الكبيرة "متشابهة في القياس" مع الدائرة المستهدفة.

2. المنهجية

تعتمد المنهجية على سلسلة مكونة من ثلاثة أجزاء، تدمج بين الزمن الطبيعي (عد الزلازل)، وتعلم المجموعات (Ensemble Learning)، وتحليل منحنى خصائص تشغيل المستقبِل (ROC).

أ. الفرضيات والتعريفات الأساسية

قانون غوتنبرغ-ريشتر (GR): الفرضية الأساسية هي أن علاقة GR ( $N = 10^{a - bM}$ ) تسري على جميع المناطق. ويُفترض أن قيمة $b$ (الميل) متشابهة عبر المجموعة والدائرة المستهدفة، وإن كانت قد تختلف قليلاً ( $\pm 0.05$ ).
الزمن الطبيعي: يُعرف بأنه عدد الزلازل الصغيرة ( $M_S$ ) التي وقعت منذ آخر زلزال كبير ( $M_T$ ).
مجموعة المناطق: تعمل سلسلة من المناطق المستطيلة المتوسعة والمتمركزة حول نقطة الاهتمام (لوس أنجلوس) كمجموعة تدريب. ويُفترض أن هذه المناطق هي "نسخ" إحصائية للدائرة المستهدبة عند قياسها بشكل مناسب.

ب. دالة الآنية (Nowcast Function)

يستخدم المؤلفون دالة آنية $\Phi(n)$ تعتمد على إحصاءات بواسون (Poisson statistics) لنمذجة احتمال وقوع زلزال كبير بعد $n$ من الأحداث الصغيرة:
$\Phi(n) = 1 - e^{-(n/N)}$
حيث $N$ هو العدد المتوقع من الأحداث الصغيرة بين الزلازل الكبيرة، والمستمد من قانون GR.

ج. سير عمل المجموعة (Ensemble Workflow)

بناء الدورات: يتم تقسيم البيانات إلى "دورات زلزالية"، حيث تبدأ كل دورة وتنتهي بحدث ذي حجم مستهدف.
القياس (Scaling): يتم قياس أوقات التوقع ( $T_F$ ) لأعضاء المجموعة بالنسبة للدائرة المستهدفة بناءً على نسبة معدلات الأحداث الصغيرة ( $R_i / R_C$ ).
بناء منحنى ROC: لكل عضو في المجموعة، يتم بناء منحنيات ROC شرطية عن طريق مسح قيم العتبة مقابل دالة الآنية:
- الإيجابي الحقيقي (TP): الآنية > العتبة وَ وقوع الزلزال الكبير التالي ضمن نافذة التوقع.
- الإيجابي الكاذب (FP): الآنية > العتبة وَ عدم وقوع زلزال.
- السالب الحقيقي/الكاذب: يتم تعريفهما عكسياً.
متوسط المجموعة: يتم حساب متوسط منحنى ROC والانحراف المعياري عبر جميع أعضاء المجموعة لتحديد القيمة التنبؤية الإيجابية (PPV)، والتي تمثل احتمالية التوقع.

د. "تحويل الآنية" (الابتكار الرئيسي)

لاختبار فرضية أن إحصائيات المجموعة تطابق الدائرة المستهدبة بصرامة، قدم المؤلفون تحويل الآنية (Nowcast Transform):

الغرض: إجبار إحصائيات GR (تحديداً قيمة $b$ ) لأعضاء المجموعة على مطابقة تلك الخاصة بالمنطقة الدائرية تماماً.
الآلية:
1. حساب قيمة الآنية $\Phi(n)$ لفاصل زمني في المجموعة باستخدام قيمة $b$ المحلية الخاصة بها.
2. عكس الدالة باستخدام قيمة $b$ الخاصة بالدائرة المستهدبة ( $b_C$ ) لإيجاد فاصل زمني محول $n'$ .
3. تعديل الوقت: إذا كان $n' \neq n$ $n^{'} \neq = n$ ، تقوم الخوارزمية بتعديل الطوابع الزمنية للأحداث الصغيرة داخل الفاصل الزمني:
  - إذا كان $n' < n$ : إزالة الأحداث عشوائياً.
  - إذا كان $n' > n$ : استكمال (Interpolation) أوقات أحداث جديدة بين الجيران الحاليين.
- هذا يحافظ على هيكل الكثافة الزمنية (مثل تجمعات الهزات الارتدادية) مع فرض التجانس الإحصائي.

هـ. حساب احتمالية التجاوز

تحسب الورقة منحنيات التجاوز الشرطية (توزيعات الناجي) لـ:

الحجم: احتمال أن يتجاوز الحدث التالي حجماً معيناً $M$ .
الزمن الطبيعي: وقت الانتظار من حيث عدد الأحداث الصغيرة.
الزمن التقويمي: وقت الانتظار بالسنوات، والمستمد من قياس فترات الزمن الطبيعي عكسياً باستخدام عامل قياس المعدل.

3. المساهمات الرئيسية

التوقع المحلي عبر المجموعات: إثبات إمكانية حساب احتمالات الزلازل المحلية بقوة باستخدام فقط أعداد الزلازل الصغيرة من مجموعة محيطة من المناطق، مما يتجاوز الحاجة إلى سجلات تاريخية محلية طويلة.
تحويل الآنية: تطوير خوارزمية مبتكرة لتعديل إحصائيات المجموعة لتتطابق مع المنطقة المستهدفة، مما يؤكد صحة فرضية "التشابه في القياس". وأظهرت النتائج أن التوقعات باستخدام البيانات المحولة متسقة عموماً مع البيانات غير المحولة، مما يعزز موثوقية الطريقة.
تحليل المقياس الزمني المزدوج: توفير إطار شامل لحساب احتمالات التجاوز في كل من الزمن الطبيعي (عد الأحداث) والزمن التقويمي، مع تسليط الضوء على العلاقة غير الخطية بينهما.
ملفات مخاطر الحجم والزمن: توليد منحنيات تجاوز محددة للحجم (مستويات 25%، 50%، 75%) وأوقات الانتظار، مما يسمح بتقييم المخاطر الديناميكي مع وقوع زلازل صغيرة جديدة.

4. النتائج

طُبقت المنهجيات على دائرة بنصف قطر 125 كم حول لوس أنجلوس، كاليفورنيا، في أعقاب زلزال نورث ريدج ( $M6.7$ ) عام 1994.

اتساق المنهجيات: وُجد أن منحنيات القيمة التنبؤية الإيجابية (PPV) للمجموعات الأصلية، والمفلترة (إزالة القيم المتطرفة)، والمحولة متسقة بشكل عام.
- قيم PPV النهائية كانت 29% (الأصلية)، و35% (المفلترة)، و43% (المحولة).
- يشير هذا إلى أن منهج المجموعة قوي حتى لو لم تكن الإحصائيات الأساسية موحدة تماماً.
تجاوز الحجم: مع زيادة عدد الزلازل الصغيرة (بالابتعاد عن آخر زلزال كبير)، تزداد احتمالية وقوع حدث عالي الحجم بشكل حاد، لا سيما لمستويات التجاوز المتوسطة (50%) والمنخفضة (25%).
التباعد بين الزمن التقويمي والزمن الطبيعي:
- الزمن الطبيعي: مع زيادة عدد الأحداث الصغيرة، ينخفض الزمن الطبيعي المتبقي المتوقع للزلزال الكبير التالي (بديهياً، النظام أصبح "أقرب" إلى الفشل).
- الزمن التقويمي: يزداد وقت الانتظار المتوقع بالسنوات التقويمية مع مرور الوقت منذ آخر حدث.
- التفسير: يعزو المؤلفون ذلك إلى تصلب القشرة الأرضية (Crustal Stiffening). مع مرور الوقت، تغلق الشقوق الصغيرة، مما يزيد من صلابة القشرة. هذا يغير آلية التشوه من الانزلاق غير المستقر (stick-slip) إلى الانزلاق المستقر، مما يسبب "تباطؤاً شاذًا" في الزمن التقويمي مقارنة بعملية بواسون البسيطة.

5. الأهمية

الآنية التشغيلية: توفر الورقة أداة عملية قائمة على البيانات لـ "الآنية" (تقييم الحالة الراهنة) والتنبؤ قصير المدى بالمخاطر الزلزالية في مراكز حضرية محددة، حتى في الحالات التي تفتقر فيها المنطقة إلى سجلات تاريخية طويلة للأحداث الكبيرة.
إدارة المخاطر: من خلال تحديد احتمالات التجاوز لكل من الحجم والزمن، توفر الطريقة ملف تعريف دقيق للمخاطر يمكن أن يوجه الاستعداد للطوارئ والتخطيط للبنية التحتية.
الرؤية الفيزيائية: يقدم التباين بين توقعات الزمن الطبيعي والزمن التقويمي تفسيراً فيزيائياً لتطور منطقة الصدع (تقوية القشرة) وصولاً إلى الزلازل الكبيرة، مما يربل بين علم الزلازل الإحصائي والآليات الجيوفيزيائية.
قابلية التكرار: يوفر المؤلفون كود بايثون مفتوح المصدر ويستخدمون بيانات متاحة مجاناً من هيئة المسح الجيولوجي الأمريكية (USGS)، مما يجعل المنهجية سهلة التكرار والتوسع في مناطق أخرى.

في الختام، تصادق هذه الورقة على نهج "مجموعة المناطق" للتنبؤ بالزلازل المحلية، وتقدم تحويلاً إحصائياً صارماً لضمان الاتساق، مما يوفر إطاراً قوياً لتقدير احتمالية الزلازل المستقبلية من حيث الحجم، وعدد الأحداث، والزمن التقويمي.