OPTIMIZING SEISMIC B-VALUES IN THE JAVA REGION THROUGH VORONOI-BASED OK1993 MODELLING
DOI:
https://doi.org/10.23960/jge.v11i2.489Keywords:
b-value Mapping, Seismicity, Voronoi Modeling, Subduction Zones, OK1993 AlgorithmAbstract
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
How to Cite
Issue
Section
Citation Check
License
Copyright (c) 2025 JGE (Jurnal Geofisika Eksplorasi)

This work is licensed under a Creative Commons Attribution-NonCommercial 4.0 International License.
Authors/Readers/Third Parties can read, print and download, redistribute or republish the article (e.g. display in a repository), translate the article, download for text and data mining purposes, reuse portions or extracts from the article in other works, sell or re-use for commercial purposes, remix, transform, or build upon the material, they must distribute their contributions under the same license as the original Creative Commons Attribution-NonComercial (CC BY-NC).