An official website of the United States government
Abstract We investigate the relative importance of injection and production on the spatial‐temporal distribution of induced seismicity at the Raft River geothermal field. We use time‐series of InSAR measurements to document surface deformation and calibrate a hydro‐mechanical model to estimate effective stress changes imparted by injection and production. Seismicity, located predominantly in the basement, is induced primarily by poroelastic stresses from cold water reinjection into a shallower reservoir. The poroelastic effect of production from a deeper reservoir is minimal and inconsistent with observed seismicity, as is pore‐pressure‐diffusion in the basement and along reactivated faults. We estimate an initial strength excess of ∼20 kPa in the basement and sedimentary cover, but the seismicity rate in the sedimentary cover is four times lower, reflecting lower density of seed‐points for earthquake nucleation. Our modeling workflow could be used to assess the impact of fluid extraction or injection on seismicity and help design or guide operations.
more »
« less
InSAR-observed surface deformation in New Mexico’s Permian Basin shows threats and opportunities presented by leaky injection wells
Abstract Knowledge of aquifer dynamics, including groundwater storage changes, is key to effective groundwater resource and reservoir management. Resolving and accurate modeling of these processes requires knowledge of subsurface poroelastic properties and lateral heterogeneity within units of interest. Computationally demanding methods for determining lateral heterogeneity in poroelastic properties exist but remain difficult to practically employ. The InSAR-based detection of uplift over a New Mexico well with a casing breach provides an opportunity to determine poroelastic properties using a tractable 2D analytical plane strain solution for surface uplift created by a pressurized reservoir with overburden. Using a Bayesian inversion framework, we calculate poroelastic properties under deep (depth of well-screen) and shallow (depth of well-breach) conditions. We find that shallow injection is necessary to produce the observed deformation. However, pressure-varying forward solutions for uplift are required to reproduce the temporal evolution of deformation. For this we use realistic shallow poroelastic properties and well dynamics, which reflect the evolving injection conditions at the well breach as the casing further erodes. Analysis of individual interferograms or InSAR time series may provide insights into shallow subsurface heterogeneity or anomalous injection conditions at operating wells more rapidly than scheduled field inspections.
more »
« less
- PAR ID:
- 10518545
- Publisher / Repository:
- Scientific Reports
- Date Published:
- Journal Name:
- Scientific Reports
- Volume:
- 13
- Issue:
- 1
- ISSN:
- 2045-2322
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Groundwater mixing dynamics play a crucial role in the biogeochemical cycling of shallow wetlands. In this paper, we conducted groundwater simulations to investigate the combined effects of evaporation and local heterogeneity on mixing dynamics in shallow wetland sediments. The results show that evaporation causes groundwater and solutes to upwell from deep sediments to the surface. As the solute reaches the surface, evaporation enhances the accumulation of the solute near the surface, resulting in a higher solute concentration than in deep sediments. Mapping of flow topology reveals that local heterogeneity generates spatially varied mixing patterns mainly along preferential flow pathways. The upwelling of groundwater induced by surface evaporation through heterogeneous sediments is likely to create distinct mixing hotspots that differ spatially from those generated by lateral preferential flows driven by large‐scale hydraulic gradients, which enhances the overall mixing in the subsurface. These findings have strong implications for biogeochemical processing in wetlands.more » « less
-
Abstract Uturuncu volcano in southern Bolivia last erupted around 250 ka but is exhibiting signs of recent activity, including over 50 yr of surface uplift, elevated seismic activity, and fumarolic activity. We studied the spatial and temporal scales of surface deformation from 1992 to 2021 to better understand subsurface activity. We tracked Uturuncu’s recent deformation using interferometric synthetic aperture radar (InSAR) data and the global navigation satellite system (GNSS) station UTUR, located near Uturuncu’s summit. We observed a spatially coherent signal of uplift from 2014 to 2021 from Sentinel-1 A/B satellites that indicates the Altiplano-Puna magma body, located 19–24 km below ground level, and previously noted as the source of the large region of deformation, is still active. The ground is now uplifting at a rate of ~3 mm/yr compared to prior rates of ~10 mm/yr. We corroborated this waning uplift with in situ data from station UTUR. We combined the Sentinel-1 data with TerraSAR-X interferograms to constrain an ~25 km2 region of subsidence located 11 km SSW of Uturuncu, with a source depth of 2.1 km below ground level to an active period of ~2.5 yr with ~5 mm/yr subsidence. We developed a conceptual model that relates these varying depths and time scales of activity in a transcrustal magmatic system. We associate the surface uplift with pressurization from ascending gases and brines from magmatic reservoirs in the midcrust. We infer the existence of brine lenses in the shallow hydrothermal system based on low subsurface resistivity correlated with surface subsidence.more » « less
-
Hydrologic behavior and soil properties across forested landscapes with complex topography exhibit high variability. The interaction of groundwater with spatially distinct soils produces and transports solutes across catchments, however, the spatiotemporal relationships between groundwater dynamics and soil solute fluxes are difficult to directly evaluate. While whole-catchment export of solutes by shallow subsurface flow represents an integration of soil environments and conditions but many studies compartmentalize soil solute fluxes as hillslope vs. riparian, deep vs. shallow, or as individual soil horizon contributions. This potentially obscures and underestimates the hillslope variation and magnitude of solute fluxes and soil development across the landscape. This study determined the spatial variation and of shallow soil base cation fluxes associated with weathering reactions (Ca, Mg, and Na), soil elemental depletion, and soil saturation dynamics in upland soils within a small, forested watershed at the Hubbard Brook Experimental Forest, NH. Base cation fluxes were calculated using a combination of ion-exchange resins placed in shallow groundwater wells (0.3 – 1 m depth) located across hillslope transects (ridges to lower backslopes) and measurements of groundwater levels. Groundwater levels were also used to create metrics of annual soil saturation. Base cation fluxes were positively correlated with soil saturation frequency and were greatest in soil profiles where primary minerals were most depleted of base cations (i.e., highly weathered). Spatial differences in soil saturation across the catchment were strongly related to topographic properties of the upslope drainage area and are interpreted to result from spatial variations in transient groundwater dynamics. Results from this work suggest that the structure of a catchment defines the spatial architecture of base cation fluxes, likely reflecting the mediation of subsurface stormflow dynamics on soil development. Furthermore, this work highlights the importance of further compartmentalizing solute fluxes along hillslopes, where certain areas may disproportionately contribute solutes to the whole catchment. Refining catchment controls on base cation generation and transport could be an important tool for opening the black box of catchment elemental cycling.more » « less
-
Abstract The shallow portion of a megathrust represents the zone of first contact between two colliding plates, and its rheological properties control the seismic and tsunami hazards generated by the fault. The high cost of underwater geodetic data collection results in sparse observations, leading to limited constraints on the interseismic behavior of megathrusts. The Rakhine‐Bangladesh megathrust offers a unique opportunity to probe the behavior of the shallow megathrust as it is the only ocean‐continent subduction zone where the near‐trench region is fully accessible on land. Here, we use observations from ALOS‐2 wide‐swath imagery spanning 2015 to 2022 to conduct an InSAR timeseries analysis of the overriding plate within Bangladesh and the Indo‐Myanmar Ranges. We identify a narrow pattern of alternating uplift and subsidence associated with mapped anticlines but show that it cannot be explained by slip on the megathrust or other fault structures. Instead, we argue that the deformation is likely caused by active aseismic folding within the wedge above a shallow decollement. We show that estimates of the decollement depth derived from a viscous folding model and the observed anticline spacing are in agreement with previous seismic observations of the decollement depth across the fold belt. We suggest that the role of ductile deformation in the overriding plate in subduction zones may be more important than previously recognized.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1038/s41598-023-42696-9
