Search NASA⌕ Search

SEARCH · Search NASA

Results for “zone of influence”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Integrating State Data Assimilation and Innovative Model Parameterization Reduces Simulated Carbon Uptake in the Arctic and Boreal Region

Model representation of carbon uptake and storage is essential for accurate projection of the response of the arctic-boreal zone to a rapidly changing climate. Land model estimates of LAI and aboveground biomass that can have a marked influence on model projections of carbon uptake and storage vary substantially in the arctic and boreal zone, making it challenging to correctly evaluate model estimates of Gross Primary Productivity (GPP). To understand and correct bias of LAI and aboveground biomass in the Community Land Model (CLM), we assimilated the 8-day Moderate Resolution Imaging Spectroradiometer (MODIS) LAI observation and a machine learning product of annual aboveground biomass into CLM using an Ensemble Adjustment Kalman Filter (EAKF) in an experimental region including Alaska and Western Canada. Assimilating LAI and aboveground biomass reduced these model estimates by 58% and 72%, respectively. The change of aboveground biomass was consistent with independent estimates of canopy top height at both regional and site levels. The International Land Model Benchmarking system assessment showed that data assimilation significantly improved CLM's performance in simulating the carbon and hydrological cycles, as well as in representing the functional relationships between LAI and other variables. Here, to further reduce the remaining bias in GPP after LAI bias correction, we re-parameterized CLM to account for low temperature suppression of photosynthesis. The LAI bias corrected model that included the new parameterization showed the best agreement with model benchmarks. Combining data assimilation with model parameterization provides a useful framework to assess photosynthetic processes in LSMs.

58 GEOSCIENCES↗
DOE OSTI2023

Influence of Codeposition on Methylammonium Iodide Degradation during Vapor Transport Deposition

Vapor transport deposition (VTD) is one manufacturing-compatible technique for large-scale deposition of perovskite films. However, to date, little understanding has been given to degradation products formed by perovskites grown under vacuum. Here, we use methylammonium iodide (MAI) as an example organic precursor to identify degradation patterns through both the sublimation process and transport through a zone with temperatures much higher than MAI can thermodynamically survive. We show that degradation due to sublimation effects MAI powders and subsequent integration into methylammonium lead iodide (MAPbI3) films differently than degradation solely due to material transport through a very high-temperature zone. Additionally, we categorize common post-deposition treatments by their influence on the degradation products seen in MAPbI3 films. Furthermore, with these results we give recommendations for VTD reactor design and the feasibility of co-deposition as a film growth technique.

14 SOLAR ENERGY↗
DOE OSTI2020

Microbial Ecology of the Plant Rhizosphere (PlantMicrobe)

The growth of bioenergy crops has the potential to help fulfill the nation’s demand for energy. However, growing these crops on arable lands may lead to socioeconomic concerns. One proposed solution is to grow biofuel crops in nutrient poor soils that aren’t well suited for food-crop growth. Therefore, in order to understand the impact of nutrient poor environments on plants, there is a need to investigate model plants and their rhizospheres in these environments using high throughput molecular analyses. Here the rhizosphere can be defined as the zone around underground portion of the plant that is influenced by the plant’s roots, and can include chemical, biological, and physical influences.

09 BIOMASS FUELS↗
DOE OSTI2025

Impact of a dynamic grid mix and climate on operational carbon emissions modeling for different building typologies and climate zones

Calculating operational carbon emissions through a building’s lifecycle is complex due to the dynamic nature of influencing factors such as climate and energy grid mix. This paper introduces a novel methodology for modeling 30-year operational carbon impacts of buildings and applies this method to mid-rise office and residential typologies across various US climate zones. The method accounts for these temporal variabilities using new and scarcely cited data sources. Key findings indicate that future changes in the climate, while impactful, play a relatively modest role in operational carbon emissions compared to significant reductions with modeling scenarios using the projected decarbonization of the electricity grid. Here, the study also finds that using annual, month-hourly, or hourly grid emission factors have a minimal impact on carbon accounting, except in certain climates and program types where emission patterns do not align with a building’s energy consumption. Warmer climates like Miami, Florida and Tucson, Arizona, which rely heavily on cooling, demonstrate larger variations in carbon emissions when using higher temporal resolution emission factors. Ultimately, this study underscores the critical role of grid decarbonization in reducing long-term emissions and the importance of incorporating this variable in life cycle assessment (LCA) modeling.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2026

Capturing the fractocohesive length scale in elastomers through a statistical mechanics-based gradient enhanced damage model

Here, this study aims to examine modeling flaw sensitivity in elastomers. The direct incorporation of polymer chain statistical mechanics considerations into a continuum stretch-based gradient-enhanced damage formulation, in turn, allows a representation of diffuse chain damage and fracture events that align with known micromechanical mechanisms. Through a series of numerical experiments, we simulate crack propagation and extract the fracture energy as an output of the model, while keeping track of the micromechanical signatures of diffuse chain damage that accommodate fracture propagation and eventually influence flaw sensitivity. Finally, by combining the fracture toughness and the work to rupture, we identify a fractocohesive length of the material, corresponding to the full width of the damage process zone. As the damage-to-fracture cascade in the proposed GED model is influenced by the introduction of a length scale associated with network imperfection and long-range load transfer, the emerging relationship of the two length scales is discussed, providing a potential link between microscopic damage mechanisms and the observed macroscopic fracture response.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Impact of Pressure-Dependent Interfacial Tension and Contact Angle on Capillary Heterogeneity Trapping of CO2 in Storage Aquifers

Summary Carbon dioxide (CO2) capillary trapping increases the total amount of CO2 that can be effectively immobilized in storage aquifers. This trapping, manifesting itself as accumulated CO2 columns at a continuum scale, is because of capillary threshold effects that occur below low-permeability barriers. Considering that capillary pressure is dictated by heterogeneous pore throat size, the trapped CO2 column height and associated CO2 saturation will vary spatially within a storage aquifer. This variation will be influenced by two pressure-dependent interfacial parameters—CO2/brine interfacial tension (IFT) and CO2/brine/rock contact angle. Our objective is to understand how the pressure dependence of these two parameters affects the heterogeneity of capillary trapped CO2 at a continuum scale. Our conceptual model is a 1D two-zone system with the upper zone being a flow barrier (low permeability) and the lower zone being a flow path (high permeability). The inputs to this model include microfacies-dependent capillary pressure vs. saturation curves and permeability values. The input capillary pressure curves were collected in the literature that represents carbonate microfacies (e.g., dolograinstone) in a prevalent formation in the Permian Basin. We then used the Leverett j-function to scale the capillary pressure curve for the two zones that are assigned with the same or different microfacies. During scaling, we considered the influence of pressure on both the IFT and contact angle of CO2/brine/dolomite systems. We varied the zone permeability contrast ratio from 2 to 50. We then assumed capillary gravity equilibriums and calculated the CO2 saturation buildup corresponding to various trapped CO2 column heights. The CO2 saturation buildup is defined as the CO2 saturation in the lower layer minus that in the upper one. We found that the saturation buildup can be doubled when varying pressure in a storage aquifer, after considering pressure-dependent IFT and contact angles. Thus, assuming these two parameters to be constant across such aquifers would cause large errors in the quantification of capillary trapping of CO2. The whole study demonstrates the importance of considering pressure-dependent interfacial properties in predicting the vertical distribution of capillary trapped CO2. It has important implications in developing a better understanding of leakage risks and consequent storage safety.

Engineering↗
DOE OSTI2021

Measuring Fault Zone and Host Rock Hydraulic Properties Using Tidal Responses

Pore pressure diffusion along faults influences induced seismicity and rupture mechanics. In situ hydraulic diffusivity measurements along faults are rare and generally lower than inferred from seismicity migration. Here we use the tidal response of deep geothermal boreholes to measure fault diffusivity and permeability. Initial interpretations of the observation with a homogeneous confined aquifer model result in diffusivities of 10 -3 –10 -1 m 2 /s. However, this model mixes signals from both the conduit and the host rock. We develop a model for tidal response with a fault passing through the aquifer based on the fault-guided fracture network and solve for hydraulic properties in both the fault and the host rock. The resulting fault permeability is 2 × 10 -14 –7 × 10 -14 m 2 (90% confidence interval[CI]) and fault diffusivity is 0.08–0.33 m 2 /s (90% CI), which is 2 orders of magnitude higher than the host rock diffusivity in some wells, thus highlighting the role of faults as fluid conduits.

58 GEOSCIENCES↗
DOE OSTI2025

Linked nitrogen and carbon dynamics reveal distinct pools and patterns in a deep, weathered bedrock rhizosphere

Nitrogen is one of the most limiting nutrients to forest productivity worldwide. Recently, it has been established that diverse ecosystems source a substantial fraction of their water from weathered bedrock, leading to questions about whether root-driven nitrogen cycling extends into weathered bedrock as well. In this study, we specifically examined nitrogen dynamics using specialized instrumentation distributed across a 16 m weathered bedrock vadose zone (WBVZ) underlying an old growth forest in northern California where the rhizosphere—composed of plant roots and their associated microbiome—extends meters into rock. We documented total dissolved nitrogen (TDN), dissolved organic carbon (DOC), inorganic nitrogen (ammonium and nitrate), and CO 2 and O 2 gases every 1.5 m to 16 m depth for 2 y. We found that TDN concentrations increased with depth, were an order of magnitude greater at 15 m than in the upper 30 cm, and that the majority of TDN throughout the weathered bedrock vadose zone was organic. We also found that TDN concentrations are influenced by depth, season, and interannual precipitation patterns. Carbon isotope composition of the DOC suggests that dissolved organic matter in the WBVZ is primarily derived from plant sources, and not the nitrogen-rich bedrock. We conclude that nitrogen dynamics in the WBVZ may be driven, in part, by an active rhizosphere, meters below the base of soil, and we argue that weathered bedrock horizons may play a key role in C-N cycling in ecosystems with deep-rooted plants.

Science & Technology - Other Topics↗
DOE OSTI2022

A novel flow cell for optical particle analyzers—application to measurements of Malvern Insitec under high pressure and temperature

A flow cell is a necessary measurement interface for some important optical analyzers. In our application of the flow cell, we utilize a state-of-the-art optical analyzer (Malvern Insitec) to measure particle size distribution and concentration in-situ in a sampled flow from a novel pressurized oxygen-fuel combustion process. However, since this sampling flow is a flow of moist flue gas under high temperature and pressure, and the flow contains particles and corrosive acid gases, it is an extreme challenge to obtain a flow cell with a high optical quality that does not perturb the measurement. To address this challenge, we propose a new design for an optical flow cell. By using a unique flow field in the proposed flow cell, the measurement zone can be well defined by the sampling flow, minimizing the influence of purge flow. To demonstrate this flow cell, we have built a test system, and conduct measurements utilizing polydisperse-particle standards (10-100 µm and 1-10 µm). The results reveal that the optical windows are well protected by the purge flow field, without risk of deposition from the sampling flow, and the Malvern Insitec can measure the particle size distribution by using this flow cell, without generating sample bias.

Cheng, Mao↗
DOE OSTI2021

Coupled Biotic-Abiotic Processes Control Biogeochemical Cycling of Dissolved Organic Matter in the Columbia River Hyporheic Zone

A critical component of assessing the impacts of climate change on watershed ecosystems involves understanding the role that dissolved organic matter (DOM) plays in driving whole ecosystem metabolism. The hyporheic zone—a biogeochemical control point where ground water and river water mix—is characterized by high DOM turnover and microbial activity and is responsible for a large fraction of lotic respiration. Yet, the dynamic nature of this ecotone provides a challenging but important environment to parse out different DOM influences on watershed function and net carbon and nutrient fluxes. We used high-resolution Fourier-transform ion cyclotron resonance mass spectrometry to provide a detailed molecular characterization of DOM and its transformation pathways in the Columbia river watershed. Samples were collected from ground water (adjacent unconfined aquifer underlying the Hanford 300 Area), Columbia river water, and its hyporheic zone. The hyporheic zone was sampled at five locations to capture spatial heterogeneity within the hyporheic zone. Our results revealed that abiotic transformation pathways (e.g., carboxylation), potentially driven by abiotic factors such as sunlight, in both the ground water and river water are likely influencing DOM availability to the hyporheic zone, which could then be coupled with biotic processes for enhanced microbial activity. The ground water profile revealed high rates of N and S transformations via abiotic reactions. The river profile showed enhanced abiotic photodegradation of lignin-like molecules that subsequently entered the hyporheic zone as low molecular weight, more degraded compounds. While the compounds in river water were in part bio-unavailable, some were further shown to increase rates of microbial respiration. Together, river water and ground water enhance microbial activity within the hyporheic zone, regardless of river stage, as shown by elevated putative amino-acid transformations and the abundance of amino-sugar and protein-like compounds. This enhanced microbial activity is further dependent on the composition of ground water and river water inputs. Our results further suggest that abiotic controls on DOM should be incorporated into predictive modeling for understanding watershed dynamics, especially as climate variability and land use could affect light exposure and changes to ground water essential elements, both shown to impact the Columbia river hyporheic zone.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Finite-rate entrainment effects on nitrogen oxide (NO x ) emissions in staged combustors

Axially staged combustors offer the potential for achieving low NO x emissions from gas turbines at the elevated temperatures required to achieve at least 65% combined cycle efficiencies. While there is evidence that premixing between the main burner products and secondary stream reduces NO x , it is important to characterize the specific requirements needed to achieve low emissions at gas turbine conditions. This study examines the sensitivity of NO x to finite-rate, large-scale entrainment of the main and secondary streams under the simplifying assumption of infinitely-fast small-scale mixing. Essentially, this allows us to isolate the effects of limited large-scale entrainment rates by using a homogeneous/uniform condition for the entrained and reacting gases. We use a reduced-order reactor network model to examine a generic staged-combustor whose inputs are physical time scales that embody the finite-rate entrainment characteristics. With simulations conducted over a large parameter space at typical operating conditions (25 atm, 650 K air), the results show that, even when small-scale mixing is infinitely fast and the reaction zone is uniform, entrainment significantly affects NO x emissions due to its influence on the equivalence ratio — and thus the temperature and time — at which the entrained mixture ignites. Fuel-air staging is shown to be vital to NO x reduction in most practical cases where entrainment times exceed roughly 1 ms and the secondary stream finishes entraining before the main burner products. A constrained NO x minimization indicates that using the leanest possible main burner and re-routing as much air as possible to the secondary stage is key to low NO x with finite-rate entrainment; for example, less than 10 ppm NO x can be achieved if the secondary stage contains 20% of the total mass flow.

03 NATURAL GAS↗
DOE OSTI2025

Exploring the impact of hydro-mechanical coupling on compaction localization

In deep subsurface environments, where porous rocks are saturated with fluids, hydromechanical coupling plays a critical role in the formation and propagation of localized deformation zones, such as compaction bands. In this study, we numerically investigate the influence of fluid-rock interactions on the development of compaction localization in fluid-saturated rocks. For this purpose, it uses a numerical framework that adapts the heat equation to model the diffusion of pore pressure caused by deformation. This is achieved with specialized subroutines for material behavior and heat transfer formulated within the ABAQUS platform and by allowing changes in rock permeability during strain localization. The results highlight that hydromechanical coupling significantly alters the compaction behavior compared to dry conditions. Under high confining stress, excess pore pressure inhibits the progression of compaction fronts, while in time-dependent processes like rock creep, pore pressure transients accelerate deformation by creating feedback loops between mechanical forces and fluid flow. These findings underscore the importance of saturation in controlling the dynamics of compaction localization and provide insights into the interplay between pore pressure, deformation, and permeability. The implications of this study extend to applications in geological storage, resource extraction, and subsurface geomechanics, offering a deeper understanding of coupled hydromechanical processes in porous rocks.

58 GEOSCIENCES↗
DOE OSTI2023

Modeling Gray Zone Tactics Used Against Green Energy Supply Chains

This report documents a project demonstrating capabilities within Lawrence Livermore National Laboratory’s Cyber and Infrastructure Resilience program to place newly developed numeric simulations into the current, larger geopolitical context. The emergence of state-backed gray zone threats introduces an unprecedented challenge because they fall short of armed conflict and involve ambiguous tactics, requiring a fresh analytical approach to understanding their potential impact and countering them. In the context of our project, we analyze how China might use gray zone tactics to degrade the power of the United States, including our influence over global efforts to combat climate change.

25 ENERGY STORAGE↗
DOE OSTI2023

Joint Inversion of Surface Electrical Resistivity Tomography and Seismic Refraction Data between the 200 Areas

Geologic stratigraphy on the Hanford Site influences groundwater and contaminant migration through the aquifer system and the vadose zone. The current geologic framework model (GFM) relies heavily on a sparse distribution of borehole data in some locations to map geologic contacts and hydrologic properties in the subsurface. Non-invasive geophysical methods such as electrical resistivity tomography (ERT), transient electromagnetic surveying, and seismic imaging are being used at Hanford to map subsurface structure in areas with limited well observations. This is to develop and mature the capability of geophysical methods to aid in GFM refinement, to identify regions of subsurface complexity, and for optimal well siting. A joint inversion of co-located seismic refraction and ERT data was carried out for data collected on a ~2.3-km profile between the 200 Areas on the Hanford Site. While ERT and seismic refraction images have sensitivity to overlapping physical properties (porosity, moisture content, lithology), the resolution and physics used to acquire each of these datasets are different and therefore information can be different or mutually complementary. Performing a joint inversion provides a reasonable option for a coherent, coupled interpretation for mutually complementary datasets. Between the 200 Areas, there are few boreholes to interpret the geologic framework model, and these data sets were obtained to provide a first line of evidence toward identifying stratigraphic structure. The seismic refraction and ERT data were independently inverted during fiscal year 2022 and broadly showed a two-layer structure with a trough-like feature that is ~1 km wide and upwards of 150 m deep. The depth of the trough feature was greater in the ERT image compared to the seismic image, which indicated a maximum depth of approximately 110 m. The objective of the joint inversion described in this report was to invert the seismic refraction and ERT data together while constraining the ERT image to be structurally similar to the seismic refraction image. The approach was applied using the geophysical inverse modeling program E4D, which has the capability to invert first-arrival times from seismic refraction data and ERT resistances using a “cross-gradient” constraint. The application of cross-gradient constraints with different weights produces ERT models that show a high degree of similarity within the upper 100 m (above ~120 m elevation). None of the ERT models show an improved structural similarity to the seismic result; therefore, it is recommended that further attempts to jointly interpret these models focus on petrophysics and image resolution. Petrophysical measurements of core samples would improve knowledge of what drives the ERT response in this region and, along with downhole geophysical measurements, could be used to “ground truth” the surface-based geophysical results. Image resolution studies would provide insight into which regions of the inverted images are reliable and which regions are poorly constrained.

58 GEOSCIENCES↗
DOE OSTI2023

Interannual variability of spring and summer monsoon growing season carbon exchange at a semiarid savanna over nearly two decades

Eddy covariance measurements of land-atmosphere energy, carbon, and water exchange now span multiple decades at some sites, supporting an improved understanding of flux interannual variability (IAV) and its ecophysiological and physical controls. Most eddy covariance IAV studies have focused on temperate forest ecosystems, where carbon fluxes are large and flux records are longest – but also where IAV is much lower than in dryland regions, which have been identified as an essential driver of the trend and variability in the global terrestrial carbon sink. In this study, we leveraged 19 years of continuous micrometeorological measurements at the AmeriFlux US-SRM mesquite savanna site in southern Arizona, USA to quantify the IAV, trends, and drivers of carbon fluxes during the distinct spring and summer growing seasons. We also assessed the ability of modern satellite and land surface models to capture the IAV of seasonal water and carbon fluxes. Annual net ecosystem production (NEP) was small and highly variable (23 +/- 64 gC m –2 yr –1 ). Precipitation and associated measures of water availability determined most of the variability in NEP, largely through their influence on annual and seasonal gross ecosystem productivity (GEP) as opposed to ecosystem respiration (ER). Root-zone soil moisture captured between 73% (spring) and 85% (summer) of GEP variability and between 73% (spring) and 58% (summer) of ER variability. Throughout the study period, soil moisture and greenness increased with associated increases in GEP, ER and NEP. These trends were strongly influenced by very productive and wet summer growing seasons during the last two years, which were characterized by abundant understory grass cover. Typically, less than half of the variability in growing season GEP and evapotranspiration was captured by satellite-based estimates and land surface model simulations with local site forcing and calibration, highlighting the ongoing utility of long-term datasets to support careful model testing and improvement.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Modulation of Tropical Convection-Circulation Interaction by Aerosol Indirect Effects in Convective Self-Aggregation Simulations of a Gray Zone Global Model

Disentangling the response of tropical convective updrafts to enhanced aerosol concentrations has been challenging. Leading theories for explaining the influence of aerosol concentrations on tropical convection are based on the dynamical response of convection to changes in cloud microphysics, neglecting possible changes in the environment. In recent years, global convection-permitting models (GCPM) have been developed to circumvent problems arising from imposing artificial scale separation on physical processes associated with deep convection. Here, we use a global model in the convective gray zone that partially simulates deep convection to investigate how enhanced concentrations of aerosols that act as cloud condensate nuclei (CCN) impact tropical convection features by modulating the convection-circulation interaction. Results from a pair of idealized non-rotating radiative-convective equilibrium simulations show that the enhanced CCN concentration leads to weaker large-scale circulation, the closeness of deep convective systems to the moist cluster edges, and more mid-level cloud water at an equilibrium state in which convective self-aggregation occurred. Correspondingly, the enhanced CCN concentration modulates how the physical processes that support or oppose convective aggregation maintain the aggregated state at equilibrium. Overall, the enhanced CCN concentration facilitates the development of deep convection in a drier environment but reduces mean precipitation. Our results emphasize the importance of allowing atmospheric phenomena to evolve continuously across spatial and temporal scales in simulations when investigating the response of tropical convection to changes in cloud microphysics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Detailed chemistry modeling of rotating detonations with dilute n -heptane sprays and preheated air

Utilization of liquid fuels is crucial to enabling commercialization of rotating detonation engines (RDE) in the near future. In this study, Eulerian-Lagrangian simulations are conducted for rotating detonative combustion with dilute n-heptane sprays and preheated air. Two-dimensional flattened configuration is used and a skeletal chemical mechanism with 44 species and 112 elementary reactions for n-heptane combustion is adopted. The flow structure, droplet distribution, and thermochemical parameters in the refill zone are first analyzed. It is shown that the mixture in the refill zone is heterogeneous, including evaporating droplets, vapor, and air. When the total temperature is below 950 K, the average equivalence ratio increases with the total temperature. When it is higher than 950 K, the average equivalence ratio is almost constant. Subsequently, the chemical explosive mode analysis is applied to identify the controlling reactions and dominant combustion modes in the fuel refill zone and reaction fronts. Results demonstrate that the initiation reaction (R104: n-C 7 H 16 + O 2 2-C 7 H 15 + HO 2 ) and low-temperature reaction (R107: RO 2 R'O 2 H) are dominant in the upstream and downstream of the refill zone, respectively. The intermediate species from low-temperature chemistry, R'O 2 H, is found to be important for the chemical explosive mode in the undetonated mixture. The influence of species diffusion and dispersed droplets is further analyzed. Results show that vapor autoignition facilitated by droplet evaporation occurs in the refill zone. Finally, the effects of the air total temperature on the detonation propagation speed and RDE propulsion performance are investigated. Further, it is found that the detonation propagation speed and specific impulse increase with air total temperature. The total pressure ratio first increases and then decreases as the air total temperature increases. Moreover, when the total temperature of the preheated air is above 1,300 K, the effects of low-temperature chemistry are negligible.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2025

Coupled geomechanical investigation of depletion-induced fault reactivation

Fault reactivation during subsurface fluid production pose significant challenges to safe and sustainable resource extraction. Here, this study presents a three-dimensional coupled geomechanical framework to investigate the processes driving fault reactivation, capturing the interactions between reservoir dynamics and geomechanical responses. Verification against theoretical estimations based on linear poroelasticity confirms the model's capacity in representing reservoir background stress responses. However, the study reveals that relying solely on background stress states can underestimate or overestimate fault reactivation potential, emphasizing the importance of including localized stress perturbations such as differential compaction and stress redistribution. Applied to a fault (M1) inspired by the geological characteristics of the Groningen field, the model shows slip initiation at 2965 m depth with 16.0 MPa depletion, aligning with field observations where seismicity occurred at approximately 3 km depth after 15.8 MPa depletion. Parametric studies reveal: (1) inelastic reservoir compaction delays fault reactivation and mitigates fault slip by reducing stress concentration, (2) higher intermediate in-situ stress magnitudes decrease the Coulomb Failure Stress (CFS) increase rate and reduce fault slip, (3) larger fault offsets amplify shear stress near the offset zone, promoting earlier reactivation and longer rupture propagation, and (4) fault permeability significantly influences pressure diffusion, with low-permeability faults leading to sharper stress changes and earlier fault destabilization. These insights highlight the critical role of geological and mechanical parameters in fault reactivation and provide a predictive framework for mitigating induced seismicity risks.

Coupled geomechanical modeling↗