skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


  • NSF Public Access
  • Search Results
  • Stress‐Based and Convolutional Forecasting of Injection‐Induced Seismicity: Application to the Otaniemi Geothermal Reservoir Stimulation
Title: 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
Author(s) / Creator(s):
;
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
  1. ; (, Geophysical Research Letters)
    Abstract Oscillatory stresses are ubiquitous on Earth and other solid‐surface bodies. Tides and seasonal signals perpetually stress faults in the crust. Relating seismicity to these stresses offers fundamental insight into earthquake triggering. We present a simple model that describes seismicity rate due to perpetual oscillatory stresses. The model applies to large‐amplitude, nonharmonic, and quasiperiodic stressing. However, it is not valid for periods similar to the characteristic timeta. We show that seismicity rate from short‐period stressing scales with the stress amplitude, but for long periods with the stressing rate. Further, that background seismicity rateris equal to the average seismicity rate during short‐period stressing. We suggest thatAσ0may be underestimated if stresses are approximated by a single harmonic function. We revisit Manga et al. (2019,https://doi.org/10.1029/2019GL082892), which analyzed the tidal triggering of marsquakes and provide a rescaling of their seismicity rate response that offers a self‐consistent comparison of different hydraulic conditions. 
    more » « less
  2. ; (, Politics & Policy)
    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
  3. (, Geophysical Research Letters)
    Abstract This letter compares the predictions of two expressions proposed for the porosity evolution in the context of rate and state friction. One (Segall & Rice, 1995,https://doi.org/10.1029/95jb02403) depends only on the sliding velocity; the other (Sleep, 1995,https://doi.org/10.1029/94jb03340) depends only on the state variable. Simulations of both are similar for velocity stepping and slide‐hold‐slide experiments. They differ significantly for normal effective stress jumps at constant sliding velocity. Segall and Rice (1995,https://doi.org/10.1029/95jb02403) predicts no change in the porosity; Sleep (1995,https://doi.org/10.1029/94jb03340) does. Simulation with a spring‐block model indicates that the magnitude of rapid slip events is essentially the same for the two formulations. Variations of porosity and induced pore pressure near rapid slip events are similar and consistent with experimental observations. Predicted porosity variations during slow slip intervals and the time at which rapid slip events occur are significantly different. The simulation indicates that changes in friction stress due to pore pressure changes exceed those due to rate and state effects. 
    more » « less
  4. ; (, Journal of Geophysical Research: Biogeosciences)
    Abstract Most earth system models fail to capture the seasonality of carbon fluxes in radiation‐limited tropical evergreen forests (TEF) in the Amazon. Kim et al. (2012,https://doi.org/10.1111/j.1365-2486.2011.02629.x) first statistically incorporated a light‐controlled phenology module into an ecosystem model to improve carbon flux simulations at one TEF site. However, it is not clear how their approach can be extended to other TEF sites with different climatic conditions. Here we evaluated temporal variability in plant functional traits at three different TEF sites using a data‐conditioned stochastic parameterization method. We showed that previously studied links—between seasonal photosynthetically active radiation (PAR) and the traitsVcmax25and leaf longevity—occur across sites. We further determined that seasonal PAR could similarly drive variations in the stomatal conductance slope parameter. Differences found in temporal trait estimates among sites indicate that dynamic trait parameters cannot be applied uniformly over space, but it may be possible to extrapolate them based on climatic factors. Motivated by recent observations that physiological capacity develops as leaves mature, we built new regression models for predicting traits that not only include PAR but also an autoregressive lag term to capture observed physiological delays behind PAR‐driven phenology shifts. With our stochastic parameterization, we predicted the three sites to be carbon neutral or carbon sinks under the RCP 8.5 future climate scenario. In contrast, projections using standard static trait parameters show most of the Amazonian TEF region becoming a carbon source. We further approximated that variable traits may allow at least a third of the radiation‐limited TEF region in the Amazon to serve as a future net carbon sink. 
    more » « less
  5. ; (, Journal of Geophysical Research: Solid Earth)
    Abstract Synchronization behavior of large earthquakes (rupture of nearby faults close in time for many cycles) has been reported in many fault systems. The general idea is that the faults in the system have similar repeating intervals and are positively coupled through stress interaction. However, many details of such synchronization remain unknown. Here, we built a numerical model in the framework of rate‐and‐state friction to simulate earthquake cycles on the west Gofar fault, East Pacific Rise. Our model consists of two seismic patches separated by a barrier patch, which are constrained by seismic observations. We varied the parameters in the barrier to understand its role on earthquake synchronization. First, we found that when the barrier is relatively weak, synchronization can be achieved by afterslip or post‐seismic creep in the barrier patch. Second, static stress transfer can lead to synchronization, opposite to the suggestion by Scholz (2010,https://doi.org/10.1785/0120090309), which was based on results from a spring‐slider model using rate‐and‐state friction. Third, the width of the barrier is more important than its strength. When the barrier is narrow enough (no more than half the width of the seismic patch in our model), the system can achieve synchronization even with a very strong barrier. Fourth, for certain simulations, the interaction between the two seismic patches promotes partial rupture in the seismic patches and leads to complex behavior: the system switches from synchronized to unsynchronized over 10–20 cycles. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email