skip to main content

Attention:

The NSF Public Access Repository (NSF-PAR) system and access will be unavailable from 11:00 PM ET on Thursday, October 10 until 2:00 AM ET on Friday, October 11 due to maintenance. We apologize for the inconvenience.


Title: Hindcasting injection-induced aseismic slip and microseismicity at the Cooper Basin Enhanced Geothermal Systems Project
Abstract

There is a growing recognition that subsurface fluid injection can produce not only earthquakes, but also aseismic slip on faults. A major challenge in understanding interactions between injection-related aseismic and seismic slip on faults is identifying aseismic slip on the field scale, given that most monitored fields are only equipped with seismic arrays. We present a modeling workflow for evaluating the possibility of aseismic slip, given observational constraints on the spatial-temporal distribution of microseismicity, injection rate, and wellhead pressure. Our numerical model simultaneously simulates discrete off-fault microseismic events and aseismic slip on a main fault during fluid injection. We apply the workflow to the 2012 Enhanced Geothermal System injection episode at Cooper Basin, Australia, which aimed to stimulate a water-saturated granitic reservoir containing a highly permeable ($$k = 10^{-13} - 10^{-12}$$k=10-13-10-12$$\hbox {m}{^2}$$m2) fault zone. We find that aseismic slip likely contributed to half of the total moment release. In addition, fault weakening from pore pressure changes, not elastic stress transfer from aseismic slip, induces the majority of observed microseismic events, given the inferred stress state. We derive a theoretical model to better estimate the time-dependent spatial extent of seismicity triggered by increases in pore pressure. To our knowledge, this is the first time injection-induced aseismic slip in a granitic reservoir has been inferred, suggesting that aseismic slip could be widespread across a range of lithologies.

 
more » « less
Award ID(s):
1947448
NSF-PAR ID:
10380300
Author(s) / Creator(s):
;
Publisher / Repository:
Nature Publishing Group
Date Published:
Journal Name:
Scientific Reports
Volume:
12
Issue:
1
ISSN:
2045-2322
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    The elliptic flow$$(v_2)$$(v2)of$${\textrm{D}}^{0}$$D0mesons from beauty-hadron decays (non-prompt$${\textrm{D}}^{0})$$D0)was measured in midcentral (30–50%) Pb–Pb collisions at a centre-of-mass energy per nucleon pair$$\sqrt{s_{\textrm{NN}}} = 5.02$$sNN=5.02 TeV with the ALICE detector at the LHC. The$${\textrm{D}}^{0}$$D0mesons were reconstructed at midrapidity$$(|y|<0.8)$$(|y|<0.8)from their hadronic decay$$\mathrm {D^0 \rightarrow K^-\uppi ^+}$$D0K-π+, in the transverse momentum interval$$2< p_{\textrm{T}} < 12$$2<pT<12 GeV/c. The result indicates a positive$$v_2$$v2for non-prompt$${{\textrm{D}}^{0}}$$D0mesons with a significance of 2.7$$\sigma $$σ. The non-prompt$${{\textrm{D}}^{0}}$$D0-meson$$v_2$$v2is lower than that of prompt non-strange D mesons with 3.2$$\sigma $$σsignificance in$$2< p_\textrm{T} < 8~\textrm{GeV}/c$$2<pT<8GeV/c, and compatible with the$$v_2$$v2of beauty-decay electrons. Theoretical calculations of beauty-quark transport in a hydrodynamically expanding medium describe the measurement within uncertainties.

     
    more » « less
  2. Abstract

    A search for exotic decays of the Higgs boson ($$\text {H}$$H) with a mass of 125$$\,\text {Ge}\hspace{-.08em}\text {V}$$GeVto a pair of light pseudoscalars$$\text {a}_{1} $$a1is performed in final states where one pseudoscalar decays to two$${\textrm{b}}$$bquarks and the other to a pair of muons or$$\tau $$τleptons. A data sample of proton–proton collisions at$$\sqrt{s}=13\,\text {Te}\hspace{-.08em}\text {V} $$s=13TeVcorresponding to an integrated luminosity of 138$$\,\text {fb}^{-1}$$fb-1recorded with the CMS detector is analyzed. No statistically significant excess is observed over the standard model backgrounds. Upper limits are set at 95% confidence level ($$\text {CL}$$CL) on the Higgs boson branching fraction to$$\upmu \upmu \text{ b } \text{ b } $$μμbband to$$\uptau \uptau \text{ b } \text{ b },$$ττbb,via a pair of$$\text {a}_{1} $$a1s. The limits depend on the pseudoscalar mass$$m_{\text {a}_{1}}$$ma1and are observed to be in the range (0.17–3.3) $$\times 10^{-4}$$×10-4and (1.7–7.7) $$\times 10^{-2}$$×10-2in the$$\upmu \upmu \text{ b } \text{ b } $$μμbband$$\uptau \uptau \text{ b } \text{ b } $$ττbbfinal states, respectively. In the framework of models with two Higgs doublets and a complex scalar singlet (2HDM+S), the results of the two final states are combined to determine upper limits on the branching fraction$${\mathcal {B}}(\text {H} \rightarrow \text {a}_{1} \text {a}_{1} \rightarrow \ell \ell \text{ b } \text{ b})$$B(Ha1a1bb)at 95%$$\text {CL}$$CL, with$$\ell $$being a muon or a$$\uptau $$τlepton. For different types of 2HDM+S, upper bounds on the branching fraction$${\mathcal {B}}(\text {H} \rightarrow \text {a}_{1} \text {a}_{1} )$$B(Ha1a1)are extracted from the combination of the two channels. In most of the Type II 2HDM+S parameter space,$${\mathcal {B}}(\text {H} \rightarrow \text {a}_{1} \text {a}_{1} )$$B(Ha1a1)values above 0.23 are excluded at 95%$$\text {CL}$$CLfor$$m_{\text {a}_{1}}$$ma1values between 15 and 60$$\,\text {Ge}\hspace{-.08em}\text {V}$$GeV.

     
    more » « less
  3. Abstract

    Let Kbe a finite simplicial, cubical, delta or CW complex. The persistence map $$\textrm{PH}$$PHtakes a filter $$f:K\rightarrow \mathbb {R}$$f:KRas input and returns the barcodes of the sublevel set persistent homology of fin each dimension. We address the inverse problem: given target barcodes D, computing the fiber $$\textrm{PH}^{-1}(D)$$PH-1(D). For this, we use the fact that $$\textrm{PH}^{-1}(D)$$PH-1(D)decomposes as a polyhedral complex when Kis a simplicial complex, and we generalise this result to arbitrary based chain complexes. We then design and implement a depth-first search that recovers the polytopes forming the fiber $$\textrm{PH}^{-1}(D)$$PH-1(D). As an application, we solve a corpus of 120 sample problems, providing a first insight into the statistical structure of these fibers, for general CW complexes.

     
    more » « less
  4. Abstract

    We report on a series of detailed Breit-Pauli and Dirac B-spline R-matrix (DBSR) differential cross section (DCS) calculations for excitation of the$$5\,^2\textrm{S}_{1/2} \rightarrow 5\,^2\textrm{P}_{1/2}$$52S1/252P1/2and$$5\,^2\textrm{S}_{1/2}\rightarrow 5\,^2\textrm{P}_{3/2}$$52S1/252P3/2states in rubidium by 40 eV incident electrons. The early BP computations shown here were carried out with both 5 states and 12 states, while the DBSR models coupled 150 and 325 states, respectively. We also report corresponding results from a limited set of DCS measurements on the unresolved$$5\,^2\textrm{P}_{1/2,3/2}$$52P1/2,3/2states, with the experimental data being restricted to the scattered electron angular range 2–$$10^\circ $$10. Typically, good agreement is found between our calculated and measured DCS for excitation of the unresolved$$5\,^2\textrm{P}_{1/2,3/2}$$52P1/2,3/2states, with best accord being found between the DBSR predictions and the measured data. The present theoretical and experimental results are also compared with predictions from earlier 40 eV calculations using the nonrelativistic Distorted-Wave Born Approximation and a Relativistic Distorted-Wave model.

    Graphic abstract 
    more » « less
  5. Abstract

    We study the structure of the Liouville quantum gravity (LQG) surfaces that are cut out as one explores a conformal loop-ensemble$$\hbox {CLE}_{\kappa '}$$CLEκfor$$\kappa '$$κin (4, 8) that is drawn on an independent$$\gamma $$γ-LQG surface for$$\gamma ^2=16/\kappa '$$γ2=16/κ. The results are similar in flavor to the ones from our companion paper dealing with$$\hbox {CLE}_{\kappa }$$CLEκfor$$\kappa $$κin (8/3, 4), where the loops of the CLE are disjoint and simple. In particular, we encode the combined structure of the LQG surface and the$$\hbox {CLE}_{\kappa '}$$CLEκin terms of stable growth-fragmentation trees or their variants, which also appear in the asymptotic study of peeling processes on decorated planar maps. This has consequences for questions that do a priori not involve LQG surfaces: In our paper entitled “CLE Percolations” described the law of interfaces obtained when coloring the loops of a$$\hbox {CLE}_{\kappa '}$$CLEκindependently into two colors with respective probabilitiespand$$1-p$$1-p. This description was complete up to one missing parameter$$\rho $$ρ. The results of the present paper about CLE on LQG allow us to determine its value in terms ofpand$$\kappa '$$κ. It shows in particular that$$\hbox {CLE}_{\kappa '}$$CLEκand$$\hbox {CLE}_{16/\kappa '}$$CLE16/κare related via a continuum analog of the Edwards-Sokal coupling between$$\hbox {FK}_q$$FKqpercolation and theq-state Potts model (which makes sense even for non-integerqbetween 1 and 4) if and only if$$q=4\cos ^2(4\pi / \kappa ')$$q=4cos2(4π/κ). This provides further evidence for the long-standing belief that$$\hbox {CLE}_{\kappa '}$$CLEκand$$\hbox {CLE}_{16/\kappa '}$$CLE16/κrepresent the scaling limits of$$\hbox {FK}_q$$FKqpercolation and theq-Potts model whenqand$$\kappa '$$κare related in this way. Another consequence of the formula for$$\rho (p,\kappa ')$$ρ(p,κ)is the value of half-plane arm exponents for such divide-and-color models (a.k.a. fuzzy Potts models) that turn out to take a somewhat different form than the usual critical exponents for two-dimensional models.

     
    more » « less