OPTIMIZING SEISMIC B-VALUES IN THE JAVA REGION THROUGH VORONOI-BASED OK1993 MODELLING

Authors

  • Yudha Styawan Geophysical Engineering, Faculty of Industrial Technology, Institut Teknologi Sumatera, Indonesia

DOI:

https://doi.org/10.23960/jge.v11i2.489

Keywords:

b-value Mapping, Seismicity, Voronoi Modeling, Subduction Zones, OK1993 Algorithm

Abstract

The spatial variation of b-values in seismically active regions provides critical insight into the stress state and rupture potential of fault systems. This study focuses on the Java region and surrounding subduction zones, where detailed mapping of b-values remains uncertain despite high seismic risk. A Voronoi-based ensemble modelling framework is implemented, incorporating the Ogata-Katsura 1993 (OK1993) formulation and spatial sampling via Sobol sequences to ensure uniform partitioning. Earthquake data from 1995 onward were compiled and harmonized into moment magnitude (Mw) using conversion equations from the Indonesian Earthquake Source and Hazard Map 2017. The OK1993 model enables estimation of b-values optimized via trust-constr and initialized with maximum likelihood estimates. The results reveal that high b-values (b > 1.2) dominate offshore southwest Lampung and south of Bali, whereas low b-values (b < 0.8) appear parts of the Sumatra fault near the Sunda Strait, faults across Java, and thrusts north of Bali and Lombok. Moderate b-values (0.8–1.0) extend along the southern Java trench and may represent partially coupled megathrust segments. Interestingly, the low b-value zones may indicate locked asperities and potential seismic gap segments, especially along southern Java, where large ruptures have not occurred in recent decades. This study demonstrates the utility of spatially adaptive, data-driven approaches in capturing complex tectonic segmentation and supports their integration into future seismic hazard assessments in Indonesia, particularly in Java and its surrounding regions.

References

Aki, K. (1965). Maximum Likelihood Estimate of b in The Formula log N = a-bM and Its Confidence Limits. Bulletin of Earthquake Research, 43, 237–239.

Arimuko, A., Rohadi, S., & Rahman, A. S. (2023). Seismotectonic Studies to Determine The Recurrence of Earthquakes Mw > 7 Using A Statistical Approach and Plate Motion in The Megathrust Western Part of Java. Geotechnical and Geological Engineering, 41(2), 1397–1406. https://doi.org/10.1007/s10706-022-02342-z

Arubi, D., Zulfakriza, Rosalia, S., Sahara, D. P., & Puspito, N. T. (2022). Estimation of b-value Variation As Earthquake Precursor in Java Region With Maximum Likelihood Method. IOP Conference Series: Earth and Environmental Science, 1047(1), 012027. https://doi.org/10.1088/1755-1315/1047/1/012027

Conn, A. R., Gould, N. I. M., & Toint, P. L. (2000). Trust Region Methods. Society for Industrial and Applied Mathematics. https://doi.org/10.1137/1.9780898719857

Irsyam, M., Cummins, P. R., Asrurifak, M., Faizal, L., Natawidjaja, D. H., Widiyantoro, S., Meilano, I., Triyoso, W., Rudiyanto, A., Hidayati, S., & Ridwan, M. (2020). Development of The 2017 National Seismic Hazard Maps of Indonesia. Earthquake Spectra, 36, 112–136.

Jiang, C., Han, L., Long, F., Lai, G., Yin, F., Bi, J., & Si, Z. (2021). Spatiotemporal Heterogeneity of b Values Revealed By A Data-Driven Approach for The 17 June 2019 MS 6.0 Changning Earthquake Sequence, Sichuan, China. Natural Hazards and Earth System Sciences, 21(7), 2233–2244. https://doi.org/10.5194/nhess-21-2233-2021

Kagan, Y. Y. (1999). Universality of The Seismic Moment–Frequency Relation. Pure and Applied Geophysics, 155(2), 537–573.

Kamer, Y. & Hiemer, S. (2013). Comment on “Analysis of the b-values before and after the 23 October 2011 Mw 7.2 Van–Erciş, Turkey, earthquake”. Tectonophysics, 608, 1448–1451. https://doi.org/10.1016/j.tecto.2013.07.040

Kamer, Y. & Hiemer, S. (2015). Data-driven Spatial b Value Mapping. Geophysical Journal International, 203(1), 1043–1054.

Kelleher, J., Sykes, L. R., & Oliver, J. (1973). Possible Criteria for Predicting Earthquake Locations and Their Application to Major Plate Boundaries of the Pacific and the Caribbean. Journal of Geophysical Research, 78(14), 2547–2585.

McCaffrey, R. (2009). The Tectonic Framework of The Sumatran Subduction Zone. Annual Review of Earth and Planetary Sciences, 37, 345–366. https://doi.org/10.1146/annurev.earth.031208.100212

McCann, W. R., Nishenko, S. P., Sykes, L. R., & Krause, J. (1979). Seismic Gaps and Plate Tectonics: Seismic Potential for Major Boundaries. Pure and Applied Geophysics, 117(6), 1082–1147.

Ogata, Y. & Katsura, K. (1993). Analysis of Temporal and Spatial Heterogeneity of Magnitude Frequency Distribution Inferred From Earthquake Catalogues. Geophysical Journal International, 113(3), 727–738.

Rahayu, S. R. & Madrinovella, I. (2024). Spatial and Temporal b-value Analysis of The Yogyakarta Region Using Earthquake Data 1960–2024. JGE (Jurnal Geofisika Eksplorasi), 10(3), 191–203. https://doi.org/10.23960/jge.v10i3.468

Schorlemmer, D. & Wiemer, S. (2005). Microseismicity Data Forecast Rupture Area. Nature, 434, 1086–1086.

Schorlemmer, D., Wiemer, S., & Wyss, M. (2005). Variations in Earthquake-Size Distribution Across Different Stress Regimes. Nature, 437(7058), 539–542. https://doi.org/10.1038/nature04094

Scholz, C. H. (2015). On The Stress Dependence of The Earthquake b Value. Geophysical Research Letters, 42(5), 1399–1402.

Simons, W. J. F., Socquet, A., Vigny, C., Ambrosius, B. A. C., Haji Abu, S., Promthong, C., & Spakman, W. (2007). A Decade of GPS in Southeast Asia: Resolving Sundaland Motion and Boundaries. Journal of Geophysical Research: Solid Earth, 112(B6), 2005JB003868. https://doi.org/10.1029/2005JB003868

Sobol’, I. M. (1967). On The Distribution of Points in A Cube And The Approximate Evaluation of Integrals. USSR Computational Mathematics and Mathematical Physics, 7(4), 86–112. https://doi.org/10.1016/0041-5553(67)90144-9

Thapa, N., Dresen, G., & Goebel, T. H. (2025). Does b ‐value Increase With Pore‐Pressure? : Insights From Laboratory Experiments and Induced Seismicity. Geophysical Research Letters, 52(11), e2025GL115740. https://doi.org/10.1029/2025GL115740

Tormann, T., Wiemer, S., & Mignan, A. (2014). Systematic Survey of High‐Resolution b Value Imaging Along Californian Faults: Inference on Asperities. Journal of Geophysical Research: Solid Earth, 119(3), 2029–2054. https://doi.org/10.1002/2013JB010867

Virtanen, P., Gommers, R., Oliphant, T.E., & Haberland, M. (2020). SciPy 1.0: Fundamental Algorithms For Scientific Computing in Python. Nature Methods, 17, 261–272.

Wiemer, S. & Wyss, M. (2002). Mapping Spatial Variability of The Frequency-Magnitude Distribution of Earthquakes. Advances in Geophysics, 45, 259–302.

Downloads

Published

2025-08-22

How to Cite

Styawan, Y. (2025). OPTIMIZING SEISMIC B-VALUES IN THE JAVA REGION THROUGH VORONOI-BASED OK1993 MODELLING. JGE (Jurnal Geofisika Eksplorasi), 11(2), 109–121. https://doi.org/10.23960/jge.v11i2.489

Issue

Section

Articles

Citation Check

Similar Articles

1 2 3 4 > >> 

You may also start an advanced similarity search for this article.

Most read articles by the same author(s)