An official website of the United States government
Abstract Dynamic earthquake triggering is commonly identified through the temporal correlation between increased seismicity rates and global earthquakes that are possible triggering events. However, correlation does not imply causation. False positives may occur when unrelated seismicity rate changes coincidently occur at around the time of candidate triggers. We investigate the expected false positive rate in Southern California with globalM ≥ 6 earthquakes as candidate triggers. We compute the false positive rate by applying the statistical tests used by DeSalvio and Fan (2023),https://doi.org/10.1029/2023jb026487to synthetic earthquake catalogs with no real dynamic triggering. We find a false positive rate of ∼3.5%–8.5% when realistic earthquake clustering is present, consistent with the 95% confidence typically used in seismology. However, when this false positive rate is applied to the tens of thousands of spatial‐temporal windows in Southern California tested in DeSalvio and Fan (2023),https://doi.org/10.1029/2023jb026487, thousands of false positives are expected. The expected false positive occurrence is large enough to explain the observed apparent triggering following 70% of large global earthquakes (DeSalvio & Fan, 2023,https://doi.org/10.1029/2023jb026487), without requiring any true dynamic triggering. Aside from the known triggering from the nearby El Mayor‐Cucapah, Mexico, earthquake, the spatial and temporal characteristics of the reported triggering are indistinguishable from random false positives. This implies that best practice for dynamic triggering studies that depend on temporal correlation is to estimate the false positive rate and investigate whether the observed apparent triggering is distinguishable from the correlations that may occur by chance.
more »
« less
Stress‐Based and Convolutional Forecasting of Injection‐Induced Seismicity: Application to the Otaniemi Geothermal Reservoir Stimulation
Abstract Induced seismicity observed during Enhanced Geothermal Stimulation at Otaniemi, Finland is modeled using both statistical and physical approaches. The physical model produces simulations closest to the observations when assuming rate‐and‐state friction for shear failure with diffusivity matching the pressure build‐up at the well‐head at onset of injections. Rate‐and‐state friction implies a time‐dependent earthquake nucleation process which is found to be essential in reproducing the spatial pattern of seismicity. This implies that permeability inferred from the expansion of the seismicity triggering front (Shapiro et al., 1997,https://doi.org/10.1111/j.1365-246x.1997.tb01215.x) can be biased. We suggest a heuristic method to account for this bias that is independent of the earthquake magnitude detection threshold. Our modeling suggests that the Omori law decay during injection shut‐ins results mainly from stress relaxation by pore pressure diffusion. During successive stimulations, seismicity should only be induced where the previous maximum of Coulomb stress changes is exceeded. This effect, commonly referred to as the Kaiser effect, is not clearly visible in the data from Otaniemi. The different injection locations at the various stimulation stages may have resulted in sufficiently different effective stress distributions that the effect was muted. We describe a statistical model whereby seismicity rate is estimated from convolution of the injection history with a kernel which approximates earthquake triggering by fluid diffusion. The statistical method has superior computational efficiency to the physical model and fits the observations as well as the physical model. This approach is applicable provided the Kaiser effect is not strong, as was the case in Otaniemi.
more »
« less
- Award ID(s):
- 1822214
- PAR ID:
- 10536162
- Publisher / Repository:
- American Geophysical Union
- Date Published:
- Journal Name:
- Journal of Geophysical Research: Solid Earth
- Volume:
- 128
- Issue:
- 4
- ISSN:
- 2169-9313
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract In 2020, Arizonans approved Proposition 207, the Smart and Safe Arizona Act, which legalized recreational marijuana sales. Previous research has typically used non‐spatial survey data to understand marijuana legalization voting patterns. However, voting behavior can, in part, be shaped by geographic context, or place, which is unaccounted for in aspatial survey data. We use multiscale geographically weighted regression to analyze how place shaped Proposition 207 voting behavior, independently of demographic variations across space. We find significant spatial variability in the sensitivity of voting for Proposition 207 to changes in several of the predictor variables of opposition and support for recreational marijuana legalization. We argue that local statistical modeling approaches provide a more in‐depth understanding of ballot measure voting behavior than the current use of global models. Related ArticlesBranton, Regina, and Ronald J. McGauvran. 2018. “Mary Jane Rocks the Vote: The Impact of Climate Context on Support for Cannabis Initiatives.”Politics & Policy46(2): 209–32.https://doi.org/10.1111/polp.12248.Brekken, Katheryn C., and Vanessa M. Fenley. 2020. “Part of the Narrative: Generic News Frames in the U.S. Recreational Marijuana Policy Subsystem.”Politics & Policy49(1): 6–32.https://doi.org/10.1111/polp.12388.Fisk, Jonathan M., Joseph A. Vonasek, and Elvis Davis. 2018. “‘Pot'reneurial Politics: The Budgetary Highs and Lows of Recreational Marijuana Policy Innovation.”Politics & Policy46(2): 189–208.https://doi.org/10.1111/polp.12246.more » « less
-
Abstract Research Highlight:Dri, G. F., Bogdziewicz, M., Hunter, M., Witham, J., & Mortelliti, A. (2025). Coupled effects of forest growth and climate change on small mammal abundance and body weight: Results of a 39‐year field study.Journal of Animal Ecology.https://doi.org/10.1111/1365‐2656.70114. Biodiversity is declining due to global environmental change, yet it remains challenging to assess how specific drivers, such as climate change, affect the dynamics and trends of individual species. While many studies correlate climate variables with species abundance or occurrence, few explicitly link environmental drivers to demographic processes to uncover the mechanisms behind population trends. Such insight requires long‐term data capable of revealing slow‐moving, nonlinear trends and disentangling natural variability from directional change. In a 39‐year study, Dri et al. (2025) demonstrate the power of sustained observation and mechanistic approaches by linking climate warming and forest maturation to increased acorn production, which enhanced body condition and survival in white‐footed mice, ultimately driving population increases. Their findings underscore the importance of long‐term data for identifying meaningful ecological trends and tracing the causal pathways by which biodiversity changes. Effective conservation under global change depends on two key shifts: greater investment in long‐term biodiversity monitoring and broader adoption of frameworks that explicitly connect environmental drivers to demographic responses. Together, these approaches provide the foundation for adaptive, evidence‐based conservation strategies in a rapidly changing world.more » « less
-
Abstract We provide a new algorithm for mass‐balance calculations in petrology and geochemistry based on the log‐ratio approach championed initially by John Aitchison (e.g., Aitchison, 1982,https://doi.org/10.1111/j.2517-6161.1982.tb01195.x; Aitchison, 1984,https://doi.org/10.1007/bf01029316) along with the underlying principles, mathematical frameworks, and data requirements. Log‐ratio Inversion of Mixed End‐members (LIME) is written in MATLAB and calculates phase proportions in an experiment or rock given a bulk composition, the composition of each phase, and the associated compositional uncertainties. An important advantage of LIME is that performing the mass‐balance calculation in inverse log‐ratio space constrains phase proportions to be between 0 and 100 wt.%. Further, the resulting LIME phase proportions provide realistic estimates of uncertainty regardless of data distribution. These two characteristics of LIME improve upon standard multiple linear regression techniques, which may yield negative values for phase proportions if non‐constrained or report oversimplified symmetric errors. Primary applications of LIME include estimating phase abundances, calculating melting and metamorphic reaction stoichiometries, and checking for open system behavior in phase equilibria experiments. The technique presented here covers whole‐rock analysis, mineralogy, and phase abundance, but could be extended to isotopic tracers, trace element modeling, and regolith component un‐mixing. We highlight the importance of uncertainty estimations for phase abundances to the fields of petrology and geochemistry by comparing our results from LIME to previously published mass‐balance calculations. Furthermore, we present case studies that demonstrate the role of mass‐balance calculations in determining magma crystallinity and defining melting reactions.more » « less
-
Induced seismicity in Oklahoma and South Kansas has been widely attributed to wastewater disposal into the deep Arbuckle formation. However, the relative contributions of pore-pressure diffusion and poroelastic stress changes to earthquake triggering remain debated. In this study, we apply the Coulomb threshold rate-and-state seismicity forecasting model of Heimisson et al. (2022) to induced seismicity in the region from 2000 to 2024. Our model is informed by poroelastic stress changes resulting from wastewater injection between 1995 and 2024 and is benchmarked against existing seismicity forecast models. Despite its simplicity, our model accurately reproduces the onset, peak, and decline of seismicity, demonstrating strong agreement with the observed earthquake activity in space and time. It provides robust constraints on permeability, yielding a range consistent with previously reported values. Based on the fit to the data, the model informed by poroelastic stress changes performs better. However, regardless of the assumed mechanism, both models yield similarly reliable seismicity forecasts, indicating that the choice of mechanism has a limited impact on forecasting performance. Finally, we estimate the probability of an >= 5 event occurring between 2021 and 2024 to range from 7% to 18% and conclude that seismic risk will remain elevated if wastewater injection volumes into the Arbuckle persist at similar levels in the coming years.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript
- Journal Article:
- https://doi.org/10.1029/2022jb024960
