Skip to main content

Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

  • Article
  • Published:

Stress-driven recurrence and precursory moment-rate surge in caldera collapse earthquakes

Abstract

Predicting the recurrence times of earthquakes and understanding the physical processes that immediately precede them are two outstanding problems in seismology. Although geodetic measurements record elastic strain accumulation, most faults have recurrence intervals longer than available measurements. Foreshocks provide the principal observations of processes before mainshocks, but variability between sequences limits generalizations of pre-failure behaviour. Here we analyse seismicity and deformation data for highly characteristic caldera collapse earthquakes from 2018 Kīlauea Volcano (Hawaii, USA), with a mean recurrence interval of 1.4 days. These events provide a unique test of stress-induced earthquake recurrence and document processes preceding mainshocks with magnitude greater than five. We show that recurrence intervals are well predicted by stress histories inferred from near-field deformation measurements and that cycle-averaged seismicity reveals a critical phase, minutes before mainshocks, where earthquakes grew larger and seismic moment rate surged dramatically. The average moment rate in the final 15 minutes (0.7% of the mean cycle duration) was 4.75 times the background, a highly significant change. We infer that as the average stress increased, ruptures were more likely to overcome geometric barriers and grow larger, leading to characteristic, whole-fault ruptures. These findings imply that stress heterogeneity influences both earthquake nucleation and growth, including on potentially hazardous tectonic faults.

This is a preview of subscription content, access via your institution

Access options

Buy this article

Prices may be subject to local taxes which are calculated during checkout

Fig. 1: Summary of observations.
Fig. 2: Test of inter-event recurrence times.
Fig. 3: Magnitude as a function of time to mainshock.
Fig. 4: Changes in earthquake behaviour before mainshocks.
Fig. 5: Schematic illustration of the effect of stress heterogeneity on earthquake nucleation and rupture extent.

Similar content being viewed by others

Data availability

GNSS data are available from the UNAVCO archive (https://www.unavco.org/data/data.html). Seismic data are from the USGS HVO (https://doi.org/10.7914/SN/HV). Data can be accessed at http://ds.iris.edu/mda/HV. Any use of trade, firm or product names is for descriptive purposes only and does not imply endorsement by the US Government.

References

  1. Reid, H. F. The elastic-rebound theory of earthquakes. Univ. Calif. Publ. Bull. Dept. Geol. 6, 413–444 (1911).

    Google Scholar 

  2. Scholz, C. H. The Mechanics Of Earthquakes and Faulting (Cambridge Univ, Press, 2019).

  3. Schwartz, D. P. & Coppersmith, K. J. Fault behavior and characteristic earthquakes: examples from the Wasatch and San Andreas fault zones. J. Geophys. Res. Solid Earth 89, 5681–5698 (1984).

    Article  Google Scholar 

  4. Weldon, R. et al. Wrightwood and the earthquake cycle: what a long recurrence record tells us about how faults work. GSA Today https://doi.org/10.1130/1052-5173 (2004).

  5. Sawai, Y. et al. Challenges of anticipating the 2011 Tohoku earthquake and tsunami using coastal geology. Geophys. Res. Lett. 39, L21309 (2012).

  6. Ogata, Y., Utsu, T. & Katsura, K. Statistical features of foreshocks in comparison with other earthquake clusters. Geophys. J. Int. 121, 233–254 (1995).

    Article  ADS  Google Scholar 

  7. Bouchon, M., Durand, V., Marsan, D., Karabulut, H. & Schmittbuhl, J. The long precursory phase of most large interplate earthquakes. Nat. Geosci. 6, 299–302 (2013).

    Article  ADS  CAS  Google Scholar 

  8. Gulia, L. & Wiemer, S. Real-time discrimination of earthquake foreshocks and aftershocks. Nature 574, 193–199 (2019).

    Article  ADS  CAS  PubMed  Google Scholar 

  9. van der Elst, N. J. B-positive: a robust estimator of aftershock magnitude distribution in transiently incomplete catalogs. J. Geophys. Res. Solid Earth 126, e2020JB021,027 (2021).

    Article  Google Scholar 

  10. Rivière, J., Lv, Z., Johnson, P. & Marone, C. Evolution of b-value during the seismic cycle: insights from laboratory experiments on simulated faults. Earth Planet. Sci. Lett. 482, 407–413 (2018).

    Article  ADS  Google Scholar 

  11. Fang, Z. & Dunham, E. M. Additional shear resistance from fault roughness and stress levels on geometrically complex faults. J. Geophys. Res. Solid Earth 118, 3642–3654 (2013).

    Article  ADS  Google Scholar 

  12. Schorlemmer, D., Wiemer, S. & Wyss, M. Variations in earthquake-size distribution across different stress regimes. Nature 437, 539–542 (2005).

    Article  ADS  CAS  PubMed  Google Scholar 

  13. Neal, C. et al. The 2018 rift eruption and summit collapse of Kīlauea Volcano. Science 363, 367–374 (2019).

  14. Lai, V. H., Zhan, Z., Brissaud, Q., Sandanbata, O. & Miller, M. S. Inflation and asymmetric collapse at Kīlauea summit during the 2018 eruption from seismic and infrasound analyses. J. Geophys. Res. Solid Earth 126, e2021JB022139 (2021).

    Article  ADS  Google Scholar 

  15. Murray, J. & Langbein, J. Slip on the San Andreas fault at Parkfield, California over two earthquake cycles and the implications for seismic hazard. Bull. Seismol. Soc. Am. 96, S283–S303 (2006).

    Article  Google Scholar 

  16. Anderson, K. & Johanson, I. Incremental caldera collapse at Kīlauea Volcano recorded in ground tilt and high-rate GNSS data, with implications for collapse dynamics and the magma system. Bull. Volcanol. 84, 1–26 (2022).

    Article  Google Scholar 

  17. Wang, T. et al. Physics‐based model reconciles caldera collapse induced static and dynamic ground motion: application to Kīlauea 2018. Geophys. Res. Lett. 49, e2021GL097440 (2022).

    Article  ADS  Google Scholar 

  18. Dietterich, H. R. et al. Lava effusion rate evolution and erupted volume during the 2018 Kīlauea lower east rift zone eruption. Bull. Volcanol. https://doi.org/10.1007/s00445-021-01443-6 (2021).

  19. Anderson, K. R. et al. The 2018 eruption of Kīlauea: insights, puzzles, and opportunities for volcano science. Annu. Rev. Earth. Planet. Sci. https://doi.org/10.1146/annurev-earth-031621-075925 (2023).

  20. Wang, T. A. et al. Ring fault creep drives volcano-tectonic seismicity during caldera collapse of Kīlauea in 2018. Earth Planet. Sci. Lett. 618, 118288 (2023).

    Article  CAS  Google Scholar 

  21. Fildes, R. A. et al. Interevent seismicity statistics associated with the 2018 quasiperiodic collapse events at Kīlauea, HI, USA. Earth Space Sci. 7, e2019EA000766 (2020).

    Article  ADS  Google Scholar 

  22. Tepp, G. Material failure and caldera collapse: insights from the 2018 Kilauea eruption. Earth Planet. Sci. Lett. 553, 116621 (2021).

    Article  CAS  Google Scholar 

  23. Roman, A. & Lundgren, P. Dynamics of large effusive eruptions driven by caldera collapse. Nature 592, 392–396 (2021).

    Article  ADS  CAS  PubMed  Google Scholar 

  24. Segall, P. & Anderson, K. Repeating caldera collapse events constrain fault friction at the kilometer scale. Proc. Natl Acad. Sci. USA 118, e2101469118 (2021).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Segall, P. et al. Caldera collapse geometry revealed by near‐field GPS displacements at Kīlauea Volcano in 2018. Geophys. Res. Lett. 47, e2020GL088867 (2020).

    Article  ADS  Google Scholar 

  26. Kumagai, H. Very-long-period seismic signals and caldera formation at Miyake Island, Japan. Science 293, 687–690 (2001).

    Article  ADS  CAS  PubMed  Google Scholar 

  27. Shelly, D. R. & Thelen, W. A. Anatomy of a caldera collapse: Kīlauea 2018 summit seismicity sequence in high resolution. Geophys. Res. Lett. 46, 14395–14403 (2019).

    Article  ADS  Google Scholar 

  28. Gibbons, S. J. & Ringdal, F. The detection of low magnitude seismic events using array-based waveform correlation. Geophys. J. Int. 165, 149–166 (2006).

    Article  ADS  Google Scholar 

  29. Butler, R. Volcanic earthquake foreshocks during the 2018 collapse of Kīlauea Caldera. Geophys. J. Int. 220, 71–78 (2020).

    Article  ADS  Google Scholar 

  30. Sagy, A., Brodsky, E. E. & Axen, G. J. Evolution of fault-surface roughness with slip. Geology 35, 283–286 (2007).

    Article  ADS  Google Scholar 

  31. Shi, Y. & Bolt, B. A. The standard error of the magnitude-frequency b value. Bull. Seismol. Soc. Am 72, 1677–1687 (1982).

    Article  Google Scholar 

  32. Cattania, C. & Segall, P. Precursory slow slip and foreshocks on rough faults. J. Geophys. Res. Solid 126, e2020JB020430 (2021).

    Article  ADS  Google Scholar 

  33. Gailler, L. et al. 3D electrical conductivity imaging of Halema'uma'u lava lake (Kīlauea volcano). J. Volcanol. Geotherm. Res. 381, 185–192 (2019).

    Article  ADS  CAS  Google Scholar 

  34. Langbein, J. Preliminary report on the 28 September 2004, M 6.0 Parkfield, California earthquake. Seismol. Res. Lett. 76, 10–26 (2005).

    Article  Google Scholar 

  35. Marzocchi, W. & Sandri, L. A review and new insights on the estimation of the b-value and its uncertainty. Ann. Geophys. 46, 1271–1282 (2003).

  36. Deichmann, N. Theoretical basis for the observed break in ML/Mw scaling between small and large earthquakes. Bull. Seismol. Soc. Am. 107, 505–520 (2017).

    Article  Google Scholar 

Download references

Acknowledgements

We thank B. Ellsworth for helpful discussions, A. Flinders and N. Van Der Elst of the USGS for internal reviews and W. Marzocchi for improving the manuscript. This study was supported by the National Science Foundation EAR-2040425 (P.S.).

Author information

Authors and Affiliations

Authors

Contributions

P.S.: conceptualization, methodology, writing—original draft, supervision. M.V.M.: data analysis, statistical computations, modelling and interpretation, contributions to writing and editing. D.R.S.: earthquake location and catalogue analysis, editing. T.A.W.: analyses of time/spectral-domain comparisons between VT and VLP events, editing. K.R.A.: analysis of geodetic data, editing.

Corresponding author

Correspondence to Paul Segall.

Ethics declarations

Competing interests

The authors declare no competing interests.

Peer review

Peer review information

Nature Geoscience thanks the anonymous reviewers for their contribution to the peer review of this work.

Additional information

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Extended data

Extended Data Fig. 1 Comparison of velocity waveforms from closely located VT and VLP earthquakes.

a, map view of Kīlauea summit, with VT and VLP epicenters, as well as three broadband stations. b,c,d, first 10 seconds of vertical velocity seismograms at stations STCD, MLOD, HLPD. VLP time shifted to match VT.

Extended Data Fig. 2 Comparison of closely located VT and VLP vertical component spectrograms at station MLOD.

a, VT b, VLP.

Extended Data Fig. 3 Cumulative seismic moment in the final 15 minute stacked catalog.

Moments are multiples of the moment of an M = 2.4 event, which is 4.47 × 1012 Nm. Red line shows exponential fit.

Extended Data Fig. 4 Mainshock epicentral regions are active in the 30 minutes prior to mainshocks.

Epicenters of earthquakes in 26 of the last 29 collapse cycles. (26 depicted VLP locations are from Shelly and Thelen (2019). Three VLPs locations, in Cycles 34, 41, and 46, are not available in the Shelley-Thelen catalog.) Epicenters colored by time to collapse. Light blue: background time to failure > 30 minutes, dark blue: last 30 minutes; red stars mainshocks epicenters. Note that nearly all mainshocks are preceded by nearby VT events in the final 30 minutes, supporting the hypothesis that mainshocks are runaway VT earthquakes.

Supplementary information

Supplementary Information

Supplementary Discussion, Figs. 1–5 and Table 1.

Rights and permissions

Springer Nature or its licensor (e.g. a society or other partner) holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Segall, P., Matthews, M.V., Shelly, D.R. et al. Stress-driven recurrence and precursory moment-rate surge in caldera collapse earthquakes. Nat. Geosci. 17, 264–269 (2024). https://doi.org/10.1038/s41561-023-01372-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1038/s41561-023-01372-3

Search

Quick links

Nature Briefing

Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.

Get the most important science stories of the day, free in your inbox. Sign up for Nature Briefing