Search NASA⌕ Search

SEARCH · Search NASA

Results for “Spatial correlation”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Atomic Structural Features of Stacking Faults and Domain Connections in the Li- and Mn-Rich Cathode

Li- and Mn-rich layered oxides (LMRs), a class of earth-abundant materials for rechargeable Li-ion battery cathodes, crystallize into layered structures of two different symmetries: C2/m represented by Li 2 MnO 3 and R$\overline{3}$m represented by LiMn 0.5 Ni 0.5 O 2 . Fundamental questions about how the C2/m and R$\overline{3}$m domains spatially correlate within the same oxide grain and how the C2/m stacking faults arrange themselves when this happens still remain. Here, by using integrated differential phase contrast imaging in scanning transmission electron microscopy (STEM-iDPC), we probe the structural and compositional details of a prototypical, cobalt-free LMR material, 0.3Li 2 MnO 3 •0.7LiMn 0.5 Ni 0.5 O 2 (Li 1.13 Mn 0.57 Ni 0.3 O 2 ). The connection between the C2/m and R$\overline{3}$m domains is found to be abrupt, facilitated by the small lattice mismatch between the two structures. Stacking faults in the C2/m domains feature atomic plane shifting that accommodates stacking sequence changes, which explains why the stacking faults form in a random manner. Furthermore, a local disordering mechanism was identified to correlate with the C2/m stacking faults. Chemically, it is found that Ni coexists with Mn at the transition metal sites within the nominal Li 2 MnO 3 domain. As a result, this study demonstrates that STEM-iDPC is a very useful tool for capturing all the elements in a single image, revealing atomic details on domain connections and stacking faults in the LMRs.

Li-rich and Mn-rich oxides↗
DOE OSTI2020

Predictive Inverse Model for Advective Heat Transfer in a Short–Circuited Fracture: Dimensional Analysis, Machine Learning, and Field Demonstration

Identifying fluid flow maldistribution in planar geometries is a well–established problem in subsurface science/engineering. Of particular importance to the thermal performance of enhanced (or “engineered”) geothermal systems is identifying the existence of nonuniform (i.e., heterogeneous) permeability and subsequently predicting advective heat transfer. Here, machine learning via a genetic algorithm (GA) identifies the spatial distribution of an unknown permeability field in a two–dimensional Hele–Shaw geometry (i.e., parallel plates). The inverse problem is solved by minimizing the L2 norm between simulated residence time distribution (RTD) and measurements of an inert tracer breakthrough curve (BTC) (C–Dot nanoparticle). Principal component analysis (PCA) of spatially correlated permeability fields enabled reduction of the parameter space by more than a factor of 10 and restricted the inverse search to reservoir–scale permeability variations. Thermal experiments and tracer tests conducted at the mesoscale Altona Field Laboratory (AFL) demonstrate that the method accurately predicts the effects of extreme flow channeling on heat transfer in a single bedding–plane rock fracture. However, this is only true when the permeability distributions provide adequate matches to both tracer RTD and frictional pressure loss. Without good agreement to frictional pressure loss, it is still possible to match a simulated RTD to measurements, but subsequent predictions of heat transfer are grossly inaccurate. Here, the results of this study suggest that it is possible to anticipate the thermal effects of flow maldistribution, but only if both simulated RTDs and frictional pressure loss between fluid inlets and outlets are in good agreement with measurements.

42 ENGINEERING↗
DOE OSTI2020

Scaling Analysis of Two–Phase Flow in Fractal Permeability Fields

Fluid mixing in permeable media is essential in many practical applications. The mixing process is a consequence of velocity fluctuations owing to geological heterogeneities and mobility contrast of fluids. Heterogeneities in natural rocks are often spatially correlated, and their properties, such as permeability, may be described using fractal distributions. This work models the fractal characteristics of such permeability fields in which the covariance function is expressed as a power-law function. A generalized scaling relation is derived relating various fractal permeability fields using the magnitude of their fluctuations. Here, this relation reveals the self-similar behavior of two-phase flow in such permeable media. To that end, a recently developed, high-resolution numerical simulator is employed to validate the analytically derived scaling relations. Two flow problems are considered in which flow is governed by 1) a linear, and 2) a nonlinear transport equation. Due to the probabilistic representation of the fractal permeability fields, a sensitivity study is conducted for each flow scenario to determine the number of realizations required for statistical convergence. Scaling analysis is performed using ensemble averages of simulated saturation profiles and their mixing lengths. Results support the validity of the developed scaling relation across the range of investigated flow conditions at intermediate times. The dynamics of linear flow in the asymptotic regime is affected by the correlation structure of heterogeneity. In nonlinear flow, scaling behavior appears to be dominated by the degree of nonlinearity.

58 GEOSCIENCES↗
DOE OSTI2022

Shallow Aseismic Slip in the Delaware Basin Determined by Sentinel-1 InSAR

The Delaware Basin, Texas is currently a hot-spot of induced seismicity and ground deformation due to fluid extraction and injection associated with horizontal drilling techniques; however, the driving mechanism behind the seismicity and deformation remains under debate. Here, using vertical and east-west horizontal surface deformation measurements derived from Sentinel-1 interferometric synthetic aperture radar (InSAR), we show that the subsurface responds differently to oil and gas activity in the northern and southeastern portions of the basin. In the north, where there is little seismicity, deformation patterns display long-wavelengths and equidimensional patterns. In contrast, the southeast region hosts most of the seismicity and displays spatial deformation patterns with narrow linear features that strike parallel to the maximum principal horizontal stress and to trends in seismicity, suggesting movement along normal faults. We model a linear deformation feature using edge dislocations and show that the InSAR observations can be reproduced by slip on normal faults contained within the Delaware Mountain Group (DMG), the formation that hosts local wastewater injection and the majority of earthquakes. Our model consists of three parallel, high-angle normal faults, with two dipping toward one another in a graben structure. Slip magnitudes reach up to 25 cm and are spatially correlated with injection wells. Measured seismicity can only explain ~2% of the fault motion predicted by our fault model, suggesting that slip leading to the deformation is predominantly aseismic. We conclude that seismic and aseismic fault motion in the southeastern Delaware Basin is likely driven by wastewater injection near critically-stressed normal faults within the DMG.

58 GEOSCIENCES↗
DOE OSTI2021

Impacts of Pore Network-Scale Wettability Heterogeneity on Immiscible Fluid Displacement: A Micromodel Study

The mixed-wet nature of reservoir formations imposes a wide range of rock wettability from strong resident-fluid wetting to strong invading-fluid wetting. The characteristics of two-phase flow in porous media composed of mixed-wetting surfaces remain poorly understood. In this study, we investigated the displacement of resident ethylene glycol (EG) by hexane in two mixed-wet micromodels of identical 2.5-D geometry heterogeneity, with uniformly or heterogeneously distributed patches strongly wetting to hexane. These patches are mixed among pores with unaltered EG-wetting surfaces. Along with control tests in the originally EG-wet micromodel, here we show the classic fingering and transitions in flow regimes at $logCa$ (capillary number) from -7.2 to -3.9. Moreover, pore-scale distributions of wettability and their spatial correlation influence displacement efficiency. In the two mixed-wet micromodels, we found (a) an increase of steady-state hexane saturation at the end of experiments by up to 0.12 in the capillary fingering regime and a decrease of at most by 0.06 in the viscous fingering regime, compared to the EG-wet micromodel, and (b) dispersed and fragmented hexane distribution after displacement. Brine drainage during supercritical CO 2 (scCO 2 ) injections in these micromodels occurs with lower wettability contrasts, and under similar viscosity ratios and interfacial tensions resulted in higher displacement efficiency relative to displacement of EG by hexane. While mixed-wettability can enhance displacement efficiency compared to uniform wettability, the dynamics of immiscible fluids in strong mixed-wet reservoirs are expected to be less pronounced in contributing to the efficiency of geological CO 2 sequestration, oil recovery, and remediation of hydrocarbon-contaminated aquifers.

2.5-D micromodel↗
DOE OSTI2022

Three‐Dimensional Permeability Inversion Using Convolutional Neural Networks and Positron Emission Tomography

Abstract Quantification of heterogeneous multiscale permeability in geologic porous media is key for understanding and predicting flow and transport processes in the subsurface. Recent utilization of in situ imaging, specifically positron emission tomography (PET), enables the measurement of three‐dimensional (3‐D) time‐lapse radiotracer solute transport in geologic media. However, accurate and computationally efficient characterization of the permeability distribution that controls the solute transport process remains challenging. Leveraging the relationship between local permeability variation and solute advection behavior, an encoder‐decoder based convolutional neural network (CNN) is implemented as a permeability inversion scheme using a single PET scan of a radiotracer pulse injection experiment as input. The CNN can accurately capture the 3‐D spatial correlation between the permeability and the radiotracer solute arrival time difference maps in geologic cores. We first test the inversion accuracy using synthetic test datasets and then test the accuracy on a suite of experimental PET imaging datasets acquired on four different geologic cores. The network‐predicted permeability maps from the geologic cores are used to parameterize forward numerical models that are directly compared with the experimental PET imaging data. The results indicate that a single trained network can generate robust 3‐D permeability inversion maps in seconds. Numerical models parameterized with these permeability maps closely capture the experimentally observed solute arrival time behavior. This work provides an unprecedented approach for efficiently characterizing multiscale permeability heterogeneity in complex geologic samples.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Magnetotelluric Investigations of the Kilauea Volcano, Hawaii

In 2002 and 2003 a collaborative effort was undertaken between Lawrence Berkeley National Laboratory, Sandia National Laboratories, the U.S. Geological Survey (USGS) Menlo Park, the USGS Hawaiian Volcano Observatory, and Electromagnetic Instruments Inc. to study the Kilauea volcano in Hawaii using the magnetotelluric (MT) technique. The work was motivated by a desire to improve understanding of the magma reservoirs and conduits within Kilauea and the East and Southwest Rift zones, which has implications for understanding Kilauea's plumbing system. An improved understanding of the rift zones has implications in understanding large-scale landslides that are generated in the Hilina Slump, which produce significant impacts on coastal communities. Up to eight stations operated simultaneously, with multiple remote reference sites, and data were processed using multi-station robust processing techniques. Data was acquired at 70 sites over the Southwest and East rift zones. Good to excellent quality data were obtained even in the harshest conditions, such as those encountered on the fresh lava flows of the East Rift Zone, where electrical contact resistances are on the order of 100 kΩ. A three-dimensional (3D) MT model study was done to guide interpretation of the observed MT measurements. Synthetic modeling demonstrates that conductive bodies in the upper 3 km can be spatially resolved where MT station sampling is good. Resistivity anomalies in the 3D inversions have a high degree of spatial correlation with previously published seismic velocity anomalies beneath Kilauea. Melt fractions between 0.096 and 0.117 are calculated for the Kilauea and Pu'u'ō'ō low resistivity anomalies, respectively.

58 GEOSCIENCES↗
DOE OSTI2023

Examining Infrasound Propagation at High Spatial Resolution Using a Nodal Seismic Array

Infrasound—acoustic waves in the atmosphere below 20 Hz—is a useful monitoring tool. Topography and atmospheric structure strongly control infrasound propagation, and at common source–receiver distances neither of these effects can be ignored when quantitative source constraints are sought. Detailed spatial measurements of the infrasound wavefield would inform propagation models and improve source estimates. However, the “large-N” deployment strategy now well-known in seismology has not yet been realized for infrasound studies. Here, we use the 900-node seismic array from the 2014 Imaging Magma Under St. Helens (iMUSH) experiment as a proxy for a large-N infrasound network, by leveraging acoustic–seismic coupled arrivals. The active-source component of iMUSH consisted of 23 shallowly buried explosions around Mount Saint Helens volcano; these explosions produced epicentral infrasound recorded on the nodes. We find that the bulk presence of ground-coupled infrasound on the nodes is controlled by wind noise and source–receiver distance, with observed arrivals for eight explosions. Explosions with the most extensive coupling produce complex spatial waveform patterns across the array. These patterns are related to both topographic and atmospheric propagation effects, as well as spatially variable site (coupling) effects. We compare our observations to simple topographic diffraction and high-resolution wind advection models, and full-wave numerical simulations. We find strong spatial correlations between (a) coupled arrival strength and modeled topographic obstruction and (b) coupled arrival time and along-path winds. Our seismoacoustic analyses and results are applicable to other existing and future nodal seismic data sets and can expand the utility of such deployments.

58 GEOSCIENCES↗
DOE OSTI2024

Trends of Sediment Resuspension and Budget in Southern Lake Michigan Under Changing Wave Climate and Hydrodynamic Environment

Sediment suspension and transport driven by waves and currents play a significant role in both the ecological and physical environments of large lakes. Lake Michigan has faced a rapidly increasing water level associated with intensified wind waves in the past decade. To investigate the spatiotemporal characteristics of suspended sediment concentration (SSC) and associated coastal sediment budgets in southern Lake Michigan, a 30-year (1991–2020) hindcast was performed using a coupled wave-current-sediment model (SWAN-FVCOM-CSTMS). We found that in southern Lake Michigan, the basin-wide mean SSC increased, and the coastal sediment loss accelerated dramatically, corresponding with intensified waves, currents and lake water level rises over the past decade. The basin-wide mean SSC, coastal sediment loss, wave height, wind speed, current speed, and water level in southern Lake Michigan are highly correlated. Spatially, the results reveal decreases in coastal SSC and sediment loss in the western portion of the southern basin, while the eastern sectors show an increase in both metrics. This reflects a clear shift in the wave climate and hydrodynamic environment. The alterations in long-term coastal sediment budgets imply that considerable shoreline transformations are being influenced by modifications in the wave climate. Understanding the spatiotemporal characteristics of SSC and coastal sediment budgets is crucial for strategic water resource management and coastal infrastructure planning.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Advancing Organized Convection Representation in the Unified Model: Implementing and Enhancing Multiscale Coherent Structure Parameterization

To address the effect of stratiform latent heating on meso- to large-scale circulations, an enhanced implementation of the Multiscale Coherent Structure Parameterization (MCSP) is developed for the Met Office Unified Model. MCSP represents the top-heavy stratiform latent heating from under-resolved organized convection in general circulation models. We couple the MCSP with a mass-flux convection scheme (CoMorph-A) to improve storm lifecycle continuity. The improved MCSP trigger is specifically designed for mixed-phase deep convective cloud, combined with a background vertical wind shear, both known to be crucial for stratiform development. We also test a cloud top temperature dependent convective-stratiform heating partitioning, in contrast to the earlier fixed partitioning. Assessments from ensemble weather forecasts and decadal simulations demonstrate that MCSP directly reduces cloud deepening and precipitation areas by moderating mesoscale circulations. Indirectly, it amends tropical precipitation biases, notably correcting dry and wet biases over India and the Indian Ocean, respectively. Remarkably, the scheme outperforms a climate model ensemble by improving seasonal precipitation cycle predictions in these regions. The scheme also improves Madden-Julian Oscillation (MJO) spectra, achieving better alignment with observational and reanalysis data by intensifying the simulated MJO over the Indian Ocean during phases 4 to 5. However, the scheme increases precipitation overestimation over the Western Pacific. Shifting from fixed to temperature-dependent convective-stratiform partitioning reduces the Pacific precipitation overestimation and further improves the seasonal cycle in India. Spatially correlated biases highlight the necessity for advances beyond deterministic approaches to align MCSP with environmental conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Development of the United States GReenhouse Gas and Air Pollutants Emissions System (GRA 2 PES)

In the U.S., emissions of greenhouse gases and air pollutants are often developed independently. Here, we describe the GReenhouse gas And Air Pollutants Emissions System (GRA 2 PES), which provides gridded emissions of fossil-fuel carbon dioxide (ffCO 2 ) and 93 air quality (AQ) species for 17 combustion and non-combustion sectors at 4 km × 4 km spatial resolution across the contiguous US. We find that the AQ emissions most spatially correlated with ffCO 2 are nitrogen oxides (NO x , ρ = 0.67), followed by sulfur dioxide (SO 2 , ρ = 0.51), carbon monoxide (CO, ρ = 0.44), and fine particulate matter (PM 2.5 , ρ = 0.38). We evaluate GRA 2 PES ffCO 2 emissions with an ensemble of publicly available regional and global inventories at national (Normalized Mean Bias (NMB) = +1.4%), state (NMB = +1.5%, R 2 = 0.98), and urban (NMB = +11.5%, R 2 = 0.97) scales. Nationally, the differences of publicly available inventories from the ensemble average range from −10.0% to +5.7%, and consistency diverges at state and urban scales. We simulate GRA 2 PES ffCO 2 in a particle dispersion model and compare to measurements of radiocarbon ( 14 C)-derived ffCO 2 collected in Los Angeles (August 2021), with results suggesting that GRA 2 PES ffCO 2 may be low by 19% for this city, but well within model-observation differences for other publicly available inventories (−43% to +94%). GRA 2 PES AQ/ffCO 2 ratios converted to concentration space generally agree with field observations (NMB = +4%, log R 2 = 0.90). Lastly, we present a method by which to utilize GRA 2 PES to derive AQ emission fluxes from ffCO 2 emissions.

Lyu, Congmeng [National Oceanic and Atmospheric Ad↗
DOE OSTI2021

A new hypothesis for the origin of Amazonian Dark Earths

Amazonian Dark Earths (ADEs) are unusually fertile soils characterised by elevated concentrations of microscopic charcoal particles, which confer their distinctive colouration. Frequent occurrences of pre-Columbian artefacts at ADE sites led to their ubiquitous classification as Anthrosols (soils of anthropic origin). However, it remains unclear how indigenous peoples created areas of high fertility in one of the most nutrient-impoverished environments on Earth. Here, we report new data from a well-studied ADE site in the Brazilian Amazon, which compel us to reconsider its anthropic origin. The amounts of phosphorus and calcium—two of the least abundant macronutrients in the region—are orders of magnitude higher in ADE profiles than in the surrounding soil. The elevated levels of phosphorus and calcium, which are often interpreted as evidence of human activity at other sites, correlate spatially with trace elements that indicate exogenous mineral sources rather than in situ deposition. Stable isotope ratios of neodymium, strontium, and radiocarbon activity of microcharcoal particles also indicate exogenous inputs from alluvial deposition of carbon and mineral elements to ADE profiles, beginning several thousands of years before the earliest evidence of soil management for plant cultivation in the region. Our data suggest that indigenous peoples harnessed natural processes of landscape formation, which led to the unique properties of ADEs, but were not responsible for their genesis. If corroborated elsewhere, this hypothesis would transform our understanding of human influence in Amazonia, opening new frontiers for the sustainable use of tropical landscapes going forward.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Grafting nanometer metal/oxide interface towards enhanced low-temperature acetylene semi-hydrogenation

Metal/oxide interface is of fundamental significance to heterogeneous catalysis because the seemingly “inert” oxide support can modulate the morphology, atomic and electronic structures of the metal catalyst through the interface. The interfacial effects are well studied over a bulk oxide support but remain elusive for nanometer-sized systems like clusters, arising from the challenges associated with chemical synthesis and structural elucidation of such hybrid clusters. We hereby demonstrate the essential catalytic roles of a nanometer metal/oxide interface constructed by a hybrid Pd/Bi 2 O 3 cluster ensemble, which is fabricated by a facile stepwise photochemical method. The Pd/Bi 2 O 3 cluster, of which the hybrid structure is elucidated by combined electron microscopy and microanalysis, features a small Pd-Pd coordination number and more importantly a Pd-Bi spatial correlation ascribed to the heterografting between Pd and Bi terminated Bi 2 O 3 clusters. The intra-cluster electron transfer towards Pd across the as-formed nanometer metal/oxide interface significantly weakens the ethylene adsorption without compromising the hydrogen activation. As a result, a 91% selectivity of ethylene and 90% conversion of acetylene can be achieved in a front-end hydrogenation process with a temperature as low as 44 °C.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Long-baseline quantum sensor network as dark matter haloscope

Ultralight dark photons constitute a well-motivated candidate for dark matter. A coherent electromagnetic wave is expected to be induced by dark photons when coupled with Standard-Model photons through kinetic mixing mechanism, and should be spatially correlated within the de Broglie wavelength of dark photons. Here we report the first search for correlated dark-photon signals using a long-baseline network of 15 atomic magnetometers, which are situated in two separated meter-scale shield rooms with a distance of about 1700 km. Both the network’s multiple sensors and the shields large size significantly enhance the expected dark-photon electromagnetic signals, and long-baseline measurements confidently reduce many local noise sources. Using this network, we constrain the kinetic mixing coefficient of dark photon dark matter over the mass range 4.1 feV-2.1 peV, which represents the most stringent constraints derived from any terrestrial experiments operating over the aforementioned mass range. Our prospect indicates that future data releases may go beyond the astrophysical constraints from the cosmic microwave background and the plasma heating.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Benchmarking quantum logic operations relative to thresholds for fault tolerance

Contemporary methods for benchmarking noisy quantum processors typically measure average error rates or process infidelities. However, thresholds for fault-tolerant quantum error correction are given in terms of worst-case error rates—defined via the diamond norm—which can differ from average error rates by orders of magnitude. One method for resolving this discrepancy is to randomize the physical implementation of quantum gates, using techniques like randomized compiling (RC). In this work, we use gate set tomography to perform precision characterization of a set of two-qubit logic gates to study RC on a superconducting quantum processor. We find that, under RC, gate errors are accurately described by a stochastic Pauli noise model without coherent errors, and that spatially correlated coherent errors and non-Markovian errors are strongly suppressed. We further show that the average and worst-case error rates are equal for randomly compiled gates, and measure a maximum worst-case error of 0.0197(3) for our gate set. Our results show that randomized benchmarks are a viable route to both verifying that a quantum processor’s error rates are below a fault-tolerance threshold, and to bounding the failure rates of near-term algorithms, if—and only if—gates are implemented via randomization methods which tailor noise.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Resolving length-scale-dependent transient disorder through an ultrafast phase transition

Abstract Material functionality can be strongly determined by structure extending only over nanoscale distances. The pair distribution function presents an opportunity for structural studies beyond idealized crystal models and to investigate structure over varying length scales. Applying this method with ultrafast time resolution has the potential to similarly disrupt the study of structural dynamics and phase transitions. Here we demonstrate such a measurement of CuIr 2 S 4 optically pumped from its low-temperature Ir-dimerized phase. Dimers are optically suppressed without spatial correlation, generating a structure whose level of disorder strongly depends on the length scale. The redevelopment of structural ordering over tens of picoseconds is directly tracked over both space and time as a transient state is approached. This measurement demonstrates the crucial role of local structure and disorder in non-equilibrium processes as well as the feasibility of accessing this information with state-of-the-art XFEL facilities.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Multi-modal Dataset of a Polycrystalline Metallic Material: 3D Microstructure and Deformation Fields

The development of high-fidelity mechanical property prediction models for the design of polycrystalline materials relies on large volumes of microstructural feature data. Concurrently, at these same scales, the deformation fields that develop during mechanical loading can be highly heterogeneous. Spatially correlated measurements of 3D microstructure and the ensuing deformation fields at the micro-scale would provide highly valuable insight into the relationship between microstructure and macroscopic mechanical response. They would also provide direct validation for numerical simulations that can guide and speed up the design of new materials and microstructures. However, to date, such data have been rare. Here, a one-of-a-kind, multi-modal dataset is presented that combines recent state-of-the-art experimental developments in 3D tomography and high-resolution deformation field measurements.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Permeability-controlled migration of induced seismicity to deeper depths near Venus in North Texas

Abstract Migration of seismic events to deeper depths along basement faults over time has been observed in the wastewater injection sites, which can be correlated spatially and temporally to the propagation or retardation of pressure fronts and corresponding poroelastic response to given operation history. The seismicity rate model has been suggested as a physical indicator for the potential of earthquake nucleation along faults by quantifying poroelastic response to multiple well operations. Our field-scale model indicates that migrating patterns of 2015–2018 seismicity observed near Venus, TX are likely attributed to spatio-temporal evolution of Coulomb stressing rate constrained by the fault permeability. Even after reducing injection volumes since 2015, pore pressure continues to diffuse and steady transfer of elastic energy to the deep fault zone increases stressing rate consistently that can induce more frequent earthquakes at large distance scales. Sensitivity tests with variation in fault permeability show that (1) slow diffusion along a low-permeability fault limits earthquake nucleation near the injection interval or (2) rapid relaxation of pressure buildup within a high-permeability fault, caused by reducing injection volumes, may mitigate the seismic potential promptly.

58 GEOSCIENCES↗