Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Geothermal”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Optimizing district energy systems by integrating Borehole Thermal Energy Storage Using a Mixed-Integer Linear Programming g-function framework with a Multi-Timescale Rolling Horizon method

Shallow geothermal has gained increasing attention in recent years; however, a reliable framework for its accurate incorporation into large-scale energy system optimization remains lacking. This study proposes a Mixed-Integer Linear Programming (MILP) framework combined with the g-function approach to integrate Borehole Thermal Energy Storage (BTES) technology into energy system optimization. Validation against a Modelica-based reservoir network simulation demonstrates that the proposed framework effectively captures the ground thermal response under varying energy loads and accurately estimates the borefield energy supply. To enhance scalability, a Rolling Horizon with Multi-Timescale (RH-MTS) method is further introduced, reducing computational time by 73 % for the 1-year optimization model with only minor loss of optimality. The framework is demonstrated through the case study of the UC Berkeley campus. Results indicate that BTES is a cost-effective and low-carbon solution: two borefields comprising 382 boreholes can meet 8.0 % and 6.6 % of the total campus heating and cooling demand, respectively, at an average energy rate of 0.70–0.77 USD/kWh and carbon intensity of 0.54 kg-CO2/kWh. Short-term analysis reveals a 35%–65% decline in BTES energy flow after 3–6 months of continuous heating/cooling operation, while long-term simulation shows that annual energy production of BTES can vary by up to 12.0 % after four years before stabilizing. Overall, this study develops a novel optimization framework that couples physics-based g-function method with MILP optimization framework, thereby advancing methodological development for shallow-geothermal integration and providing actionable guidance for BTES deployment in district-energy systems.

Yang, Jiahui↗
DOE OSTI2025

Enhancing building resilience in cold climates: Integrating heat pump technologies with renewable energy

As electrification advances and Cold Climate Heat Pump technology progresses, ensuring grid stability becomes increasingly critical for effective heating in cold climates. However, natural disasters, especially during winter, pose significant threats to grid stability, impacting the reliability of air-source heat pumps. Despite these challenges, the integration of renewable energy sources and storage solutions in heating systems has not been extensively studied within the context of resilience. Here, this paper delves into the literature on renewable-powered heat pumps to assess their potential in enhancing building resilience in U.S. cold climate zones, which are particularly susceptible to extreme weather and grid disruptions. By leveraging renewable sources—solar, geothermal, and water—in conjunction with heat pump technology and supported by thermal or battery storage, this approach aims to provide a dependable solution for maintaining indoor heating during grid failures. Our analysis begins with a review of various renewable energy sources suitable for heat pumps, followed by an exploration of their application in cold climate regions across the U.S., and discussions on potential integration strategies with heat pump systems. This study highlights the advantages and suitability of solar irradiance and geothermal resources, emphasizing the importance of tailored, site-specific assessments to maximize energy efficiency and resilience. Additionally, it outlines the economic and environmental considerations necessary for implementing such systems and identifies potential challenges and areas for future research to facilitate the broader integration of renewable energy in heating solutions for enhanced resilience.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

Effect of initial water saturation on the performance of fracturing fluids with and without polyallylamine under simulated EGS conditions

Objectives/Scope: StimuFrac (US Patents 9,873,828 B2 and 9,447,315 B2), a CO 2 - reactive polymer aqueous solution [polyallylamine (PAA) 1wt% in water] combined with CO 2 , can be used as a potentially less water-intensive fracturing fluid for enhanced geothermal systems (EGS). Our previous results show that in hot dry rock (HDR), PAA/CO 2 fracturing fluids outperformed other fluids such as water, CO 2 , and CO 2 /water in generating large fractures with less fluid consumed. The objective of this work is to investigate the effect of initial water saturation on the performance of StimuFrac by conducting hydraulic fracturing tests with ½ foot cubic rock samples held under representative EGS stress/temperature conditions and by using cyclic injection strategies (under constant injection rate). The resulting fracture hydraulic conductivities, breakdown pressures, and volumes of fluids required are compared. Methods/Procedures/Process: To simulate geothermal reservoir conditions, in all tests, the rock sample was held under triaxial confinement and at 200 °C, and different volumes of water were initially injected into the rock sample before any fracturing processes were initiated. For the single-cycle PAA (or water) alternating CO 2 (PAG or WAG) injection fracturing experiments, one complete cycle consisted of two steps: (1) injecting a PAA slug (or water slug) followed by (2) injecting CO 2 to initiate and propagate the fracture. For experiments involving multiple injection cycles, the CO 2 injection pressure is increased until it peaks and begins to decline (indicating fracture initiation at this moment), and then continued being injected for another 30 seconds to propagate the fracture. Then, these two-step cycles [injection of PAA (or water) followed by CO 2 injection (up to 2-4 mL/min)] are repeated. Applications/Significance/Novelty: The results of this study suggest that water saturation significantly affects the fracturing fluid transmission into the rock pore space, thus affecting the fracture initiation and propagation. In this study, fracturing tests via a single injection cycle or multiple injection cycles were performed. Splitting the rock samples in half after testing reveals that fracture propagation is significantly limited under high water saturation conditions (three-day initial water injection) compared to stimulation experiments performed in hot dry rock. The fractures propagate less than 1/3 of the distance from the wellbore to the outer rock surface, and in some cases, no fracture is generated. This may be caused by leak-off dominating the fracturing process and the fluid injection rate is insufficient to overcome leak-off, even under high injection rate conditions. Additionally, CO 2 could be leaking off into the wellbore annulus and this may be making it more difficult to generate sufficiently high-pressure gradients away from the near-wellbore region. Under low water saturation conditions (dry rock or after 1-day initial water injection), PAA/CO 2 consistently generated significantly larger fractures compared with the other fluids. CO 2 generated large fractures only in the hot dry rock and only when using high injection rates, though data variability is high.

58 GEOSCIENCES↗
DOE OSTI2021

4D Proxy Imaging of Fracture Dilation and Stress Shadowing Using Electrical Resistivity Tomography During High Pressure Injections into a Dense Rock Formation

Fluid flow through fractured rock systems is governed in large part by the distribution, interconnectivity, and size of fracture apertures. In-situ stress is one of the primary factors controlling fracture aperture, and one that is altered significantly during high-pressure fluid injections or extractions. Interactions between stress, pore pressure, aperture, and fluid flow can result in complex and evolving poroelastic behavior with significant implications regarding the predictability and risk of developing and managing deep subsurface reservoirs (geothermal, fossil energy, and geologic carbon sequestration). In saturated crystalline rocks, bulk electrical conductivity is highly sensitive to porosity, and therefore to fracture aperture size and distribution. We demonstrate the use of time-lapse 3D electrical resistivity tomography for remotely monitoring stress induced changes in aperture distribution during high pressure injections into a fractured rock system at a scale of tens of meters. Results reveal a complex and continuously evolving stress field involving aperture dilations in the natural fracture system and aperture contractions in adjacent zones of shadow stress. Results provide information about the spatiotemporal changes in the system behavior and point to the potential of electrical imaging for autonomously and remotely monitoring evolving stress conditions by proxy through changes in bulk electrical conductivity.

electrical resisitivity tomography, stress, 4D Tim↗
DOE OSTI2024

Thermal Equations of State of Magnesite: Implication for the Complex Mid‐Lower Mantle Seismic Scatterers

Abstract Magnesite (MgCO 3 ) entering the lower mantle together with the subducted oceanic crust is an important carbon carrier. The reaction between magnesite and mantle minerals has been documented, but its influence on the density and velocity profiles of lower mantle remains unexplored. To decipher the deep carbon transportation and its associated effect, here we determined the thermal equations of state of magnesite up to 120 GPa and 2600 K using X‐ray diffraction in laser‐heated diamond anvil cells. The obtained thermal elastic parameters of magnesite facilitated a comprehensive understanding on the influence of magnesite‐SiO 2 reaction, variation of carbon and SiO 2 content, and temperature on the origin of lower‐mantle scatterers at 1,000–1,800 km depth. Our modeling revealed that the depth of the lower‐mantleV S scatterers is mainly controlled by the Al 2 O 3 content in SiO 2 , while its magnitude depends on the SiO 2 content. Along normal geotherm, the magnesite‐SiO 2 reaction would occur before the post‐stishovite transition, consuming substantial SiO 2 in the subducted oceanic crust. Depending on the amount of residual SiO 2 , the post‐stishovite transition can produce a 2.5–5.2 (2)%V S reduction, compatible with the observed seismic scatterers in Izu‐Bonin and Mariana subduction zones. Along slab geotherm, this reaction occurs after the post‐stishovite transition, generating a greaterV S reduction of 4.4–6.4 (4)%. We thus propose that the reaction between sinking MgCO 3 and SiO 2 in the slab is one of the potential factors influencing the magnitude of the lower‐Vs scatterers at 1,000–1,900 km depth. Our results provide new insights into the deep‐mantle carbonate transportation influencing regional geophysics.

Geochemistry & Geophysics↗
DOE OSTI2026

Sequential Fracture Activation and Stress Evolution During EGS Stimulation at Utah FORGE Revealed by Waveform Cross‐Correlation

Mapping fracture networks in Enhanced Geothermal Systems (EGS) is essential for optimizing reservoir performance, yet complex fracture evolution during stimulation remains difficult to resolve. This study examines the evolution of microseismicity and fracture networks during stage 3 of the 2022 EGS stimulation at the Utah Frontier Observatory for Research in Geothermal Energy site. We map the fracture network represented by 20 clusters of seismic events identified by waveform similarities with cross-correlation. We characterize their geometric properties such as strike, dip, length, and width, and analyze the time evolution of activated fractures. The results reveal a systematic fracture evolution: early activation of pre-existing natural fractures, complex network development during peak injection, and continued activation of less favorably oriented fractures post-injection. Magnitude calibration using the Principal Component Analysis of cross-correlated waveforms improves relative amplitude measurements, refining estimations of the Gutenberg-Richter b-values with spatial variations in b-values suggesting stress re-distribution across the stimulated area. Analysis of the stress state of selected fractures further shows that fractures requiring higher excess pore pressure primarily activate at the end of injection and post-injection, highlighting stress transfer due to pore pressure as a dominant triggering mechanism. These findings provide insights into fracture propagation, stress evolution, and seismic hazard assessment in EGS reservoirs.

Asirifi, Richard [Texas A & M Univ., College Stati↗
DOE OSTI2021

The stability of subducted glaucophane with the Earth’s secular cooling

The blueschist to eclogite transition is one of the major geochemical–metamorphic processes typifying the subduction zone, which releases fluids triggering earthquakes and arc volcanism. Although glaucophane is an index hydrous mineral for the blueschist facies, its stability at mantle depths in diverse subduction regimes of contemporary and early Earth has not been experimentally determined. Here, we show that the maximum depth of glaucophane stability increases with decreasing thermal gradients of the subduction system. Along cold subduction geotherm, glaucophane remains stable down ca. 240 km depth, whereas it dehydrates and breaks down at as shallow as ca. 40 km depth under warm subduction geotherm or the Proterozoic tectonic setting. Our results imply that secular cooling of the Earth has extended the stability of glaucophane and consequently enabled the transportation of water into deeper interior of the Earth, suppressing arc magmatism, volcanism, and seismic activities along subduction zones.

58 GEOSCIENCES↗
DOE OSTI2021

Brockarchaeota, a novel archaeal phylum with unique and versatile carbon cycling pathways

Geothermal environments, such as hot springs and hydrothermal vents, are hotspots for carbon cycling and contain many poorly described microbial taxa. Here, we reconstructed 15 archaeal metagenome-assembled genomes (MAGs) from terrestrial hot spring sediments in China and deep-sea hydrothermal vent sediments in Guaymas Basin, Gulf of California. Phylogenetic analyses of these MAGs indicate that they form a distinct group within the TACK superphylum, and thus we propose their classification as a new phylum, ‘Brockarchaeota’, named after Thomas Brock for his seminal research in hot springs. Based on the MAG sequence information, we infer that some Brockarchaeota are uniquely capable of mediating non-methanogenic anaerobic methylotrophy, via the tetrahydrofolate methyl branch of the Wood-Ljungdahl pathway and reductive glycine pathway. The hydrothermal vent genotypes appear to be obligate fermenters of plant-derived polysaccharides that rely mostly on substrate-level phosphorylation, as they seem to lack most respiratory complexes. In contrast, hot spring lineages have alternate pathways to increase their ATP yield, including anaerobic methylotrophy of methanol and trimethylamine, and potentially use geothermally derived mercury, arsenic, or hydrogen. Their broad distribution and their apparent anaerobic metabolic versatility indicate that Brockarchaeota may occupy previously overlooked roles in anaerobic carbon cycling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Hot springs viruses at Yellowstone National Park have ancient origins and are adapted to thermophilic hosts

Geothermal springs house unicellular red algae in the class Cyanidiophyceae that dominate the microbial biomass at these sites. Little is known about host-virus interactions in these environments. We analyzed the virus community associated with red algal mats in three neighboring habitats (creek, endolithic, soil) at Lemonade Creek, Yellowstone National Park (YNP), USA. We find that despite proximity, each habitat houses a unique collection of viruses, with the giant viruses, Megaviricetes, dominant in all three. The early branching phylogenetic position of genes encoded on metagenome assembled virus genomes (vMAGs) suggests that the YNP lineages are of ancient origin and not due to multiple invasions from mesophilic habitats. The existence of genomic footprints of adaptation to thermophily in the vMAGs is consistent with this idea. The Cyanidiophyceae at geothermal sites originated ca. 1.5 Bya and are therefore relevant to understanding biotic interactions on the early Earth.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Identifying geologic characteristics and operational decisions to meet global carbon sequestration goals

Geologic carbon sequestration is the process of injecting and storing CO 2 in subsurface reservoirs and is an essential technology for global environmental security (e.g., climate change mitigation) and economic security (e.g., CO 2 tax credits). To meet energy, economic, and environmental goals, society will have to identify vast volumes of high-capacity, low-cost, and viable storage reservoirs for sequestering CO 2 . In turn, this requires understanding how major geologic characteristics (such as reservoir depth, thickness, permeability, porosity, and temperature) and design and operational decisions (such as injection well spacing) impact CO 2 injection rates, storage capacity, and economics. Although many numerical simulation tools exist, they cannot repeat the required thousands or millions of simulations to identify ideal reservoir properties and the sensitivity and interaction between geologic parameters and operational decisions. Here, we use SCO 2 T—a fast-running, reduced-order modeling framework—to explore the sensitivity of major geologic parameters and operational decisions to engineering (CO 2 injection rates, plume dimensions, and storage capacities and effectiveness) and costs. Our results show, for the first time, benefits and impacts such as allowing CO 2 plumes to overlap, how different well spacing patterns affect CO 2 sequestration, the effects on costs of including brine treatment and disposal, and the effect of restricting injection rates to 1 MtCO 2 per y based on well limitations. We reveal multiple novel and unintuitive findings including: (i) deeper reservoirs have reduced carbon sequestration costs until injection rates reach 1 MtCO 2 per y, at which point deeper reservoirs become more expensive, (ii) thicker formations allow for increased injection rates and storage capacity, but thickness barely impacts plume areas, (iii) higher geothermal gradients result in reduced sequestration costs, unless brine treatment/disposal costs are included, at which point reservoirs having lower geothermal gradients are more economical because they produce less brine for each unit of injected CO 2 , and (iv) allowing plumes to overlap has a significantly positive impact of increasing storage capacities but has only a small influence on reducing sequestration costs. Altogether, our results illustrate new scientific conclusions to help identify suitable sites to inject and store CO 2 , to help understand the complex interaction between geology and resulting costs, and to help support the pursuit of meeting global sequestration targets.

58 GEOSCIENCES↗
DOE OSTI2024

Novel order-level lineage of ammonia-oxidizing archaea widespread in marine and terrestrial environments

Ammonia-oxidizing archaea (AOA) are among the most ubiquitous and abundant archaea on Earth, widely distributed in marine, terrestrial, and geothermal ecosystems. However, the genomic diversity, biogeography, and evolutionary process of AOA populations in subsurface environments are vastly understudied compared to those in marine and soil systems. Here, we report a novel AOA order Candidatus (Ca.) Nitrosomirales which forms a sister lineage to the thermophilic Ca. Nitrosocaldales. Metagenomic and 16S rRNA gene-read mapping demonstrates the abundant presence of Nitrosomirales AOA in various groundwater environments and their widespread distribution across a range of geothermal, terrestrial, and marine habitats. Terrestrial Nitrosomirales AOA show the genetic capacity of using formate as a source of reductant and using nitrate as an alternative electron acceptor. Nitrosomirales AOA appear to have acquired key metabolic genes and operons from other mesophilic populations via horizontal gene transfer, including genes encoding urease, nitrite reductase, and V-type ATPase. The additional metabolic versatility conferred by acquired functions may have facilitated their radiation into a variety of subsurface, marine, and soil environments. We also provide evidence that each of the four AOA orders spans both marine and terrestrial habitats, which suggests a more complex evolutionary history for major AOA lineages than previously proposed. Together, these findings establish a robust phylogenomic framework of AOA and provide new insights into the ecology and adaptation of this globally abundant functional guild.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Distributed Thermal Response Multi-Source Modeling to Evaluate Heterogeneous Subsurface Properties

A thorough assessment of thermal properties in heterogeneous subsurface is necessary in design of low-temperature borehole heat exchangers (BHEs). For this study, a distributed thermal response test (DTRT), which combines distributed temperature sensing (DTS) with a conventional thermal response test (TRT), was conducted in a U-bend geothermal loop installed in an open borehole at the University of Illinois at Urbana-Champaign to estimate thermal properties by analyzing the thermal response of different geologic materials while applying a constant heat input rate. Fiber-optic cables in the DTRT were deployed both inside the U-bend geothermal loop and in the center of the borehole to improve the accuracy of calculated heat-loss rates and borehole temperature profile measurements. To assess the subsurface thermal conductivity during the heating phase of the DTRT, a single-source model and a multi-source model, both based on the infinite line source method, were developed using the borehole temperature data and temperatures inside and along the outside of the loop, separately. The two models returned similar thermal conductivity values. The multi-source modeling has the advantage of predicting the thermal conductivity of heterogeneous geologic materials from borehole temperature profiles during the DTRT heating phase. Additionally, based on the distributed thermal conductivity measured in the borehole, estimates were made for both radial thermal impacts and the rate of heat loss in the BHE.

58 GEOSCIENCES↗
DOE OSTI2021

Comparison of Rotary and Linear Cutting Methodology in Determining Specific Cutting Energy of Granite

Abstract Single polycrystalline diamond compact (PDC) cutting is a practical technique to understand the rock–tool interactions in drag-bit-type geothermal drilling operations. This paper introduces a rotary cutting method to determine specific cutting energy (SCE) and compares it with the conventional linear cutting method. In this work, granite is selected to represent hard rock formations in geothermal drilling. Cutting tests are conducted on a Computer Numerical Control (CNC) machine with a realistic cutting speed of 12.7 m/min and several chip loads ranging from 0.08 to 0.25 mm. The cutting force is measured using a dynamometer and then converted to specific cutting energy (SCE). The results show that the rotary method produces an inverse relationship between SCE and chip load, whereas the linear method shows the opposite. As a result, the produced SCE by the rotary method tends to be lower than that of the linear method at a higher chip load at and over 0.16 mm. The difference may be attributed to the cutting configuration and associated force components.

Engineering↗
DOE OSTI2020

Fat Crayon Toolkit

The Fat Crayon Toolkit (FCT) provides python modules for assisting in the design and visualization of boreholes and tunnels and their relationship with natural features such as fractures. Distances of closest approach among the boreholes can be calculated, along with other key design metrics. Additionally, the shear propensity of natural fractures (either deterministic or stochastic) can be visualized in either 3D (as VTK files) or in 2D (stereonets). The example provided demonstrates a test bed developed for the Enhanced Geothermal Systems (EGS) Collab project. This material is based upon work supported by the Geothermal Office of the Department of Energy and the Sanford Underground Research Facility and was prepared by LLNL under Contract DE-AC52-07NA27344.

Morris, JosephP.↗
DOE OSTI2023

T2Well-ECO2NV2 v1.0 beta

T2Well-ECO2NV2 combines the T2Well code for simulating coupled well-reservoir fluid and heat flow with the ECO2N Version 2.0 equation of state module for calculating the fluid properties of water-CO2-salt mixtures at high temperatures. The combination of these two capabilities in T2Well-ECO2NV2 creates the capability of simulating injection and withdrawal of fluids in high-Temperature geothermal or geologic carbon sequestration systems. While the previously released T2Well-ECO2N Version 1.0 was useful for geologic carbon sequestration, the new T2Well-ECO2NV2 has capabilities for higher temperature and thereby allows the simulation of CO2-rich geothermal systems.

Oldenburg, CurtisM↗
DOE OSTI2024

Advancing subsurface analysis: Integrating computer vision and deep learning for the near real-time interpretation of borehole image logs in the Illinois Basin-Decatur Project

The accurate quantification and mapping of subsurface natural fracture systems using borehole imaging logs are critical for the success of CO 2 sequestration in geologic formations, optimization of engineered geothermal systems, and hydrocarbon production enhancement. However, traditional interpretation processes suffer from time-consuming procedures and human bias. To address these challenges and expedite fracture analysis, we investigated the application of integrated computer vision and DL workflows to automate image log analysis. Specifically, the design of our workflow was crafted to swiftly detect fractures and baffles by using actual electrical resistivity of borehole wall from microresistivity imaging device alongside their binary representation. This novel approach significantly reduces computational time while providing invaluable insights. By incorporating conventional logging and microseismic data, we present a regional subsurface natural fracture mapping technique. Through the minimization of human bias in image log analysis, our automated workflow achieves reduced fracture interpretation time and costs while ensuring robust and reproducible results. We demonstrated the efficacy of our approach by applying the workflow to the Illinois Basin-Decatur Project site. The automated workflow successfully identified major fractured zones, multiple baffles, and an interbedded layer with a high resolution of 0.01 ft or 0.12 in. (0.3 cm) and can be upscaled to any desired resolution. Validation through microseismic and image log interpretations allows for accurate and near-real-time mapping of fractures and baffles, significantly enhancing CO 2 pressure forecasting and postinjection site care. Our approach stands out due to its robustness, consistency, and reduced computational cost compared with alternative feature extraction technologies. It presents exciting possibilities for advancing CO 2 sequestration and engineered geothermal efforts by offering comprehensive and efficient fracture mapping solutions. This technology can contribute significantly to the optimization of CO 2 sequestration projects, facilitating sustainable environmental practices, and combating climate change.

Geochemistry & Geophysics↗
DOE OSTI2026

Limited Dynamic Earthquake Triggering in Nevada

Dynamic triggering occurs when seismic waves from distant large earthquakes temporarily alter stress conditions along faults, potentially triggering new earthquakes hundreds to thousands of kilometers away from the source. Previous studies have linked triggered seismicity to anthropogenic activities such as geothermal, oil, and gas production. Although these activities are present in Nevada, little work has been conducted to explore dynamically triggered seismicity in Nevada. Here, we analyze a newly published, high-resolution earthquake catalog for Nevada to identify local seismicity dynamically triggered by teleseismic events (Mw≥7) from 2008 to 2023. We identify 94 dynamically triggered earthquakes concentrated in four distinct regions, which qualitatively show a modest positive correlation with geothermal well locations. Triggered seismicity in Nevada is predominantly delayed, with some instantaneously triggered by Rayleigh waves. The prevalence of delayed triggering indicates that pore fluid interactions may play a critical role in controlling dynamic triggering susceptibility in Nevada. Our results demonstrate that dynamic triggering can provide valuable insight to help identify critically stressed regions.

58 GEOSCIENCES↗
DOE OSTI2020

More Than 40 yr of Potentially Induced Seismicity Close to the San Andreas Fault in San Ardo, Central California

Evidence for fluid–injection–induced seismicity is rare in California hydrocarbon basins, despite widespread injection close to seismically active faults. We investigate a potential case of injection–induced earthquakes associated with San Ardo oilfield operations that began in the early 1950s. The largest potentially induced events occurred in 1955 (M L 5.2) and 1985 (M w 4.5) within ~6 km from the oilfield. We analyze Synthetic Aperture Radar interferometric images acquired by Sentinel–1A/B satellites between 2016 and 2020 and find surface deformation of up to 1.5 cm/yr, indicating pressure–imbalance in parts of the oilfield. Fluid injection in San Ardo is concentrated within highly permeable rocks directly above the granitic basement at a depth of ~800 m. Seismicity predominantly occurs along basement faults at 6–13 km depths. Seismicity and wastewater disposal wells are spatially correlated to the north of the oilfield. Temporal correlations are observed over more than 40 yr with correlation coefficients of up to 0.71 for seismicity within a 24 km distance from the oilfield. Such large distances have not previously been observed in California but are similar to the large spatial footprint of injection in Oklahoma. Furthermore, the San Ardo seismicity shows anomalous clustering with earthquakes consistently occurring at close spatial proximity but long interevent times. Similar clustering has previously been reported in California geothermal fields and may be indicative of seismicity driven by long–term, spatially persistent external forcing.

58 GEOSCIENCES↗