Estimasi Magnitude Of Completeness (Mc), A-Value, dan B-value untuk menilai dinamika seismotektonik gempa dangkal di wilayah Nusa Tenggara Barat menggunakan metode Entire Magnitude Range (EMR)

Indonesia

Wilayah Nusa Tenggara Barat (NTB) merupakan kawasan dengan seismisitas tinggi akibat interaksi zona subduksi Indo-Australia–Eurasia dan sistem Flores Back Arc Thrust. Penelitian ini bertujuan mengestimasi nilai Magnitude of Completeness (Mc), a-value, dan b-value gempa dangkal (kedalaman ≤ 70 km) di wilayah NTB menggunakan metode Entire Magnitude Range (EMR), serta menganalisis variasi spasial dan temporalnya terhadap dinamika seismotektonik wilayah tersebut. Data yang digunakan adalah data sekunder katalog gempa bumi dari United States Geological Survey (USGS) periode 2010–2025 sebanyak 354 kejadian, direduksi menjadi 233 kejadian melalui declustering metode Reasenberg menggunakan perangkat lunak ZMAP pada MATLAB. Hasil analisis Frequency Magnitude Distribution dengan metode EMR menghasilkan nilai Mc sebesar 4,40, b-value sebesar 1,2591 ± 0,1072, dan a-value sebesar 7,6799, dengan misfit sebesar 0,040940 yang menunjukkan kecocokan model baik. Nilai b-value lebih besar dari satu mengindikasikan dominasi gempa bermagnitudo kecil hingga menengah di wilayah NTB. Analisis spasial menunjukkan bagian barat NTB memiliki nilai a-value dan b-value relatif lebih rendah dibandingkan bagian tengah-timur, mengindikasikan akumulasi tegangan batuan lebih tinggi di wilayah barat. Analisis temporal pada tiga periode (2010–2014, 2015–2019, dan 2020–2025) menunjukkan penurunan b-value pada 2015–2019 berkaitan dengan akumulasi tegangan menjelang rangkaian Gempa Lombok 2018, diikuti peningkatan b-value dan a-value pada 2020–2025 yang mengindikasikan pelepasan tegangan serta peningkatan aktivitas kegempaan. Hasil ini menunjukkan karakteristik seismotektonik gempa dangkal di wilayah NTB bersifat dinamis secara spasial dan temporal, sehingga dapat menjadi bahan pertimbangan dalam mitigasi bencana gempa bumi di wilayah tersebut.

English

West Nusa Tenggara (NTB) region has high seismicity resulting from interaction between the Indo-Australian–Eurasian subduction zone and the Flores Back Arc Thrust system. This study aims to estimate the Magnitude of Completeness (Mc), a-value, and b-value of shallow earthquakes (depth ≤ 70 km) in NTB using the Entire Magnitude Range (EMR) method, and to analyze their spatial and temporal variations in relation to seismotectonic dynamics. The data used are secondary earthquake catalog data from the United States Geological Survey (USGS) for 2010–2025, comprising 354 events, reduced to 233 events through declustering using the Reasenberg method in ZMAP on MATLAB. The Frequency Magnitude Distribution analysis using EMR yielded an Mc of 4.40, a b-value of 1.2591 ± 0.1072, and an a-value of 7.6799, with a misfit of 0.040940 indicating a good model fit. A b-value greater than one indicates dominance of small- to moderate-magnitude earthquakes in NTB. Spatial analysis revealed that western NTB has relatively lower a-value and b-value than the central–eastern part, suggesting higher rock-stress accumulation in the west. Temporal analysis across three periods (2010–2014, 2015–2019, and 2020–2025) showed a decrease in b-value during 2015–2019, associated with stress accumulation prior to the 2018 Lombok earthquake sequence, followed by an increase in both b-value and a-value during 2020–2025, indicating stress release and increased seismic activity. These findings show that shallow-earthquake seismotectonic characteristics in NTB are dynamic both spatially and temporally, providing a basis for earthquake disaster mitigation in the region.

Arabic

تعتبر منطقة نوسا تينجارا الغربية (NTB) منطقة ذات نشاط زلزالي مرتفع نتيجة لتفاعل منطقة اندماج الصفائح الهندية-الأسترالية مع الأوراسية ونظام فلوريس باك أرك ترست. تهدف هذه الدراسة إلى تقدير قيمة حجم الاكتمال (Mc)، وقيمة a، وقيمة b للزلازل الضحلة (العمق ≤ 70 كم) في منطقة NTB باستخدام طريقة النطاق الكامل للحجم (EMR)، وكذلك تحليل التغيرات المكانية والزمنية لها بالنسبة للديناميكيات الزلزالية التكتونية للمنطقة. البيانات المستخدمة هي بيانات ثانوية من كتالوج الزلازل لهيئة المسح الجيولوجي الأمريكية (USGS) للفترة من 2010 إلى 2025 وعددها 354 حدثًا، تم تقليلها إلى 233 حدثًا من خلال إزالة التكتلات باستخدام طريقة Reasenberg عبر برنامج ZMAP علىMATLAB. نتائج تحليل توزيع تكرار القوى باستخدام طريقة EMR أسفرت عن قيمة Mc تساوي 4.40، وقيمة b تساوي 1.2591 ± 0.1072، وقيمة a تساوي 7.6799، مع خطأ قدره 0.040940 مما يشير إلى توافق جيد للنموذج. إن قيمة b الأكبر من واحد تشير إلى سيطرة الزلازل ذات القوي الصغيرة إلى المتوسطة في منطقة NTB. يشير التحليل المكاني إلى أن الجزء الغربي من NTB يحتوي على قيم a وb أقل نسبيًا مقارنة بالجزء الأوسط-الشرقي، مما يدل على تراكم أكبر للضغوط الصخرية في المنطقة الغربية. أما التحليل الزمني لثلاث فترات (2010–2014، 2015–2019، و2020–2025) فقد أظهر انخفاض قيمة b في 2015–2019 مرتبطًا بتراكم الضغوط قبيل سلسلة زلازل لومبوك 2018، تلاه ارتفاع في قيمتي b وa في 2020–2025 مما يشير إلى تحرير الضغوط وزيادة النشاط الزلزالي. تظهر هذه النتائج أن الخصائص الزلزالية التكتونية للزلازل الضحلة في منطقة NTB ديناميكية من حيث المكان. تُظهر هذه النتائج أن خصائص الزلازل الضحلة في منطقة NTB ديناميكية من الناحيتين المكانية والزمنية، لذلك يمكن أن تُعتبر مادة للنظر فيها عند التخفيف من كوارث الزلازل في تلك المنطقة