Emergent asperity localization enables physics-based prediction of surface rupture displacement
تقدم هذه الدراسة إطاراً قائماً على الفيزياء يثبت أن النتوءات الموضعية ذات الانزلاق العالي تنشأ طبيعياً من حفظ العزم الزلزالي وهندسة الصدع المحدودة بدلاً من كونها افتراضات مفروضة، مما يتيح التنبؤ الدقيق بأنماط إزاحة التصدع السطحي عبر مختلف الزلازل التاريخية والأنظمة الصدعية.
تخيل قشرة الأرض كأنها أحجية (جيكسو) ضخمة تتحرك ببطء شديد. القطع، التي تُسمى الصفائح التكتونية، تحتك ببعضها البعض باستمرار، لكنها لا تنزلق بسلاسة؛ بل تتعثر وتتوقف، مما يؤدي إلى تراكم الضغط مثل شريط مطاطي يتم شده أكثر فأكثر. وفي النهاية، ينقطع الشريط المطاطي، وتتحرر الطاقة المخزنة في هزة مفاجئة وعنيفة نسميها زلزالاً. وبينما يعلم العلماء منذ زمن طويل مدى كبر حجم الزلزال (قدره)، فإن التنبؤ بدقة بمدى انفتاح الأرض عند السطح هو لغز أكثر تعقيداً بكثير. وذلك لأن الأرض لا تتحرك بشكل متساوٍ؛ فبعض البقع تتحرك قليلاً، بينما تقفز مناطق أخرى لأمتار بعيدة. هذه "البقع الساخنة" من الحركة الهائلة تُسمى النتوءات (asperities) (تخيلها كأنها "أهداف الثيران" أو مراكز التركيز في الزلزال). لعقود من الزمن، افترض الباحثون أن هذه النتوءات كانت ناتجة عن عيوب خفية وغير منتظمة في الصخور — مثل ندبة أو بقعة ضعيفة في إطار سيارة. ولكن ماذا لو لم تكن هذه النتوءات عيوباً على الإطلاق؟ ماذا لو كانت نتيجة طبيعية وحتمية لفيزياء عملية الانكسار نفسها؟ هذا السؤال يقع في قلب دراسة جديدة لعلماء الجيوفيزياء "جيا تشينج" و"شي وي شو"، الذين يحاولون معرفة ما إذا كانت الأجزاء الأكثر خطورة في الزلزال مكتوبة في قوانين الفيزياء أم أنها مجرد حوادث جيولوجية عشوائية.
في ورقتهم البحثية الجديدة، يقترح "تشينج" و"شو" فكرة مفاجئة: النتوءات لا تحتاج إلى تخطيط مسبق. إذ يشيران إلى أن مناطق الانزلاق العالي هذه تظهر بشكل طبيعي، مثل نمط يتشكل في كوب من القهوة الدوامة، ببساطة لأن الأرض يجب أن توازن ميزانية طاقتها. ولاختبار ذلك، بنى الباحثون نموذجاً حاسوبياً قائماً على الفيزياء يعمل كأنه "محاسب كوني". لم يخبروا الحاسوب بمكان وجود "النتوءات"، ولم يبرمجوا فيه أي بقع ضعيفة أو غير منتظمة في الصخور. بدلاً من ذلك، أعطوا النموذج قاعدتين بسيطتين: يجب أن يتطابق إجمالي الطاقة المنبعثة مع الطاقة المتراكمة على مدى آلاف السنين (حفظ العزم الزلزالي)، وأن يتناسب حجم الزلزال مع الطول الفيزيائي لخط الصدع.
وعندما قاموا بتشغيل المحاكاة، حدث شيء ساحر. على الرغم من أنهم بدأوا بتوزيع ضغط منتظم ومتساوٍ، إلا أن النظام "انكسر" طبيعياً إلى حالة تركزت فيها الحركة في مناطق محددة وموضعية. لقد "اختار" الحاسوب إنشاء النتوءات لأنها كانت الطريقة الأكثر كفاءة من حيث الطاقة لتفريغ الضغط المتراكم. وجد المؤلفون أن هذا السلوك الناشئ نجح في إعادة إنتاج أنماط الانزلاق المعقدة وغير المنتظمة الملحوظة في 23 زلزالاً تاريخياً من نوع الصدع المضرب، بدءاً من زلزال كونلون عام 2001 في الصين وصولاً إلى زلازل تركيا عام 2023. لم يكن النموذج مجرد تخمين؛ بل استطاع رصد المواقع الدقيقة لأكبر القفزات والتلاشي التدريجي للحركة عند نهايات التصدع، وكل ذلك دون أن يُطلب منه البحث عن تلك المواقع.
يجادل الباحثون بأن هذا يغير نظرتنا لخطوط الصدع. فبدلاً من رؤية النتوءات كـ "ندبات" غامضة وموروثة يجب علينا البحث عنها، يقترحون أنها خصائص ناشئة — نتيجة طبيعية لكيفية تفاعل خطوط الصدع المحدودة مع قوانين الفيزياء. ولإثبات أن هذا لم يكن مجرد صدفة للأحداث الماضية، طبقوا نموذجهم على صدعين نشطين لم يشهدوا تمزقاً هائلاً في الذاكرة القريبة: صدع سان أندرياس المركزي في كاليفورنيا، وصدع آنينغ هي في الصين. وباستخدام الهندسة المعروفة لهذه الصدوع ومعدلات انزلاقها طويلة الأمد فقط، توقع النموذج أين من المرجح أن تتشكل "النتوءات" المستقبلية للإزاحة. وأظهرت النتائج أنه حتى في الصدوع ذات الإجهاد الموحد، فإن هندسة الصدع نفسه تجبر الحركة على التكتل في نقاط محددة.
هذا لا يعني أن النموذج يتنبأ بدقة بـ متى سيحدث الزلزال أو التفاصيل الدقيقة لكل صدع. ويؤكد المؤلفون بحذر أن عملهم هو إطار عمل قائم على الفيزياء يشرح تنظيم التصدع، وليس بديلاً عن عمليات المحاكاة الديناميكية المعقدة في الوقت الفعلي. ومع ذلك، فإن نتائجهم تشير إلى أننا قد لا نحتاج لمعرفة كل عيب خفي في الصخر للتنبؤ بأكثر الأماكن التي ستتحرك فيها الأرض. ومن خلال فهم أن النتوءات هي نتيجة طبيعية لحفظ الطاقة، يمكننا بناء أدوات أفضل لتقييم المخاطر الزلزالية، مما يساعد المهندسين على تصميم جسور وأنابيب أكثر أماناً يمكنها تحمل التمزقات غير المنتظمة لسطح الأرض.
بيان المشكلة تُعد إزاحة التصدع السطحي معلماً حرجاً لتقييم الأضرار الناجمة عن الزلازل والمخاطر الزلزالية الخاصة بالصدوع. ومع ذلك، لا يزال التباين على طول مسار الصدع أمراً يصعب التنبؤ به. تعتمد النهج الحالية عادةً على علاقات قياس تجريبية أو تفرض ظروفاً إجهادية واحتكاكية غير متجانسة (مثل النتوءات المحددة مسبقاً) لإعادة إنتاج أنماط الانزلاق المرصودة. وثمة سؤال جوهري لم يُحل بعد في فيزياء الزلازل، وهو ما إذا كانت نطاقات الانزلاق العالي الموضعية (النتوءات) هي نتائج متأصلة في فيزياء الزلازل ناتجة عن القيود العالمية، أم أنها مجرد آثار ناتجة عن افتراضات النموذج المفروضة. علاوة على ذلك، تتطلب المحاكاة الديناميكية الواقعية معرفة تفصيلية بالإجهادات المتغيرة مكانياً والمعاملات الاحتكاكية، وهي معلومات نادراً ما تتوفر للصدوع الطبيعية، مما يحد من فائدتها التنبؤية للأنظمة الصدعية النشطة.
المنهجية طوّر المؤلفون إطار عمل قائماً على الفيزياء حيث ينبثق التمركز الموضعي للانزلاق الشبيه بالنتوءات بشكل طبيعي من قيود حفظ العزم الزلزالي وهندسة الصدع المحدودة، بدلاً من كونه مفروضاً مسبقاً. تدمج المنهجية خمسة مكونات متتالية:
هندسة الصدع والتقطيع: يتم نمذجة الصدوع كأسطح مستطيلة محدودة يتم تقسيمها إلى أجزاء فرعية، مع تحديد التقسيم بناءً على الخرائط الجيولوجية.
حدوث زلزال متوازن العزم: يتم ضبط علاقة "غوتنبرغ-ريشتر" للتردد والمقدار المقطوعة بواسطة معدل العزم التكتوني طويل الأمد. يضمن هذا أن مجموع العزم الزلزالي المتحرر خلال دورة زلزالية يوازن التحميل التكتوني (M˙tectonic=μLWVs).
قياس تمزق الصدع المعتمد على المقدار: يتم تقدير أبعاد التمزق (الطول والعرض) للزلازل المحاكات باستخدام قوانين القياس التجريبية (مثل Leonard, 2014)، والتي تربط المقدار بمساحة التمزق.
محاكاة التمزق العشوائية: يتم محاكاة تمزقات الزلازل عبر أخذ عينات "مونت كارلو". يتم أخذ مواقع النواة عشوائياً ضمن القيود الهندسية، وتُخصص أبعاد التمزق بناءً على المقدار. لا يتم تحديد أي مواقع محددة للنتوءات أو رقع احتكاكية غير متجانسة.
تراكم الإزاحة التراكمية: يقوم النموذج بتجميع مساهمات الإزاحة المتزامنة مع الزلزال من جميع التمزقات المحاكات عبر الدورة الزلزالية. يتم تقدير الإزاحة السطحية عن طريق إسقاط الانزلاق التراكمي للصدع إلى السطح باستخدام دالة تخميد أسية لمراعاة التشوه الضحل.
ولإثبات الآلية الفيزيائية، تم تقديم نموذج متغير مختزل الرتبة. يعمل هذا النموذج على تقليل دالة الطاقة المقيدة (الطاقة المرنة + طاقة الاحتكاك) مع الخضوع لقيد ثابت للعزم الزلزالي. ويظهر أن توزيع الانزلاق المتجانس غير مستقر ميكانيكياً؛ حيث تتطور الاضطرابات الصغيرة نحو حالات انزلاق موضعية مفضلة طاقياً دون قوة خارجية.
النتائج الرئيسية
بروز النتوءات: يؤكد النموذج المتغير أن النظام يتطور طبيعياً نحو نطاقات انزلاق موضعية (نتوءات) لتقليل طاقة النظام، وذلك تحت ميزانية عزم ثابتة، مما يثبت أن الانزلاق المتجانس غير مستقر.
التحقق من البيانات التاريخية: تمت عملية التحقق من الإطار مقابل 23 زلزالاً تاريخياً من نوع الإزاحة الجانبية شملت بيئات تكتونية متنوعة، بما في ذلك أحداث كونلون 2001، ومادوو 2021، وتركيا 2023.
نجح النموذج في إعادة إنتاج التباين من الدرجة الأولى على طول مسار إزاحة السطح، بما في ذلك مواقع ذروة الإزاحة القصوى وانخفاض الانزلاق نحو نهايات التمزق.
طابقت أغلفة الإزاحة المتوقعة الملاحظات الجيولوجية والجيوديسية دون الحاجة إلى معايرة خاصة بكل حدث أو تحديد مواقع النتوءات مسبقاً.
استطاع النموذج التقاط كل من أنماط التمزق المجزأة (مثل صدع آنينغهي) وأنماط التمزق المستمرة (مثل صدع سان أندرياس) بناءً على هندسة الصدع وقيود العزم فقط.
التنبؤات المستقبلية: تم تطبيق الإطار على وسط صدع سان أندرياس وصدع آنينغهي للتنبؤ بتوزيعات الإزاحة المستقبلية.
ولدت عمليات المحاكاة نطاقات عالية الإزاحة موضعية محكومة بتداخل التمزقات وتقسيم الصدع.
أظهر وسط صدع سان أندرياس نطاقات إزاحة أوسع وأكثر سلاسة بسبب الاستمرارية الهندسية، بينما أظهر صدع آنينغهي تمركزاً أقوى مرتبطاً بالقطع الفردية بسبب الانقطاعات الهندسية.
تعتمد هذه التنبؤات فقط على هندسة الصدع المخططة، ومعدلات الانزلاق، وإحصائيات التكرار، مما يثبت قدرة النموذج على تقييم المخاطر الاستباقية.
الأهمية والادعاءات يدعي البحث تقديم أساس فيزيائي موحد للتنبؤ بإزاحة التصدع السطحي، مما ينقل تفسير النتوءات من كونها سمات هيكلية مفروضة إلى خصائص ناشئة عن الأنظمة ذات الصدوع المحدودة.
إطار عمل موحد: يربط البحث بين التحميل التكتوني طويل الأمد، وإحصائيات الزلازل، وهندسة التمزق، والإزاحة السسطحية ضمن صيغة واحدة. ويوضح أن أنماط الإزاحة الموضعية يمكن تفسيرها من خلال التفاعل بين هندسة الصدع المحدودة وحفظ العزم الزلزالي، بدلاً من الحاجة إلى عدم تجانس خاص بكل حدث.
القدرة التنبؤية: يوفر الإطار أداة عملية لتقييم التباين المكاني لتصدع السطح المستقبلي، خاصة في الصدوع التي تفتقر إلى سجلات تاريخية كبيرة. كما يتيح تقييماً للمخاطر الزلزالية القائمة على الإزاحة باستخدام الكميات الجيولوجية القابلة للقياس فقط.
التكامل: يضع المؤلفون هذا النهج كعنصر مكمل لمحاكاة التمزق الديناميكي. فبينما تُعد النماذج الديناميكية ضرورية لدراسة سرعة التمزق وحركة الأرض في المجال القريب، فإن هذا الإطار المختزل القائم على الحفظ يعالج التنظيم واسع النطاق لإزاحة التمزق دون الحاجة إلى معاملات احتكاكية تفصيلية وغالباً ما تكون مجهولة.
تخلص الدراسة إلى أن تمركز النتوءات ليس مجرد سمة موروثة من عدم تجانس الصدع، بل هو خاصية ناشئة عن الأنظمة ذات الصدوع المحدودة، مما يوفر فهماً ميكانيكياً لتعقيد التمزق الزلزالي ويحسن تقييم المخاطر الزلزالية الخاصة بالصدوع.