Search NASA⌕ Search

SEARCH · Search NASA

Results for “linear estimation”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
NASA NTRS1983

Polarized-interferometer feasibility study

The feasibility of using a polarized-interferometer system as a rendezvous and docking sensor for two cooperating spacecraft was studied. The polarized interferometer is a radio frequency system for long range, real time determination of relative position and attitude. Range is determined by round trip signal timing. Direction is determined by radio interferometry. Relative roll is determined from signal polarization. Each spacecraft is equipped with a transponder and an antenna array. The antenna arrays consist of four crossed dipoles that can transmit or receive either circularly or linearly polarized signals. The active spacecraft is equipped with a sophisticated transponder and makes all measurements. The transponder on the passive spacecraft is a relatively simple repeater. An initialization algorithm is developed to estimate position and attitude without any a priori information. A tracking algorithm based upon minimum variance linear estimators is also developed. Techniques to simplify the transponder on the passive spacecraft are investigated and a suitable configuration is determined. A multiple carrier CW signal format is selected. The dependence of range accuracy and ambiguity resolution error probability are derived and used to design a candidate system. The validity of the design and the feasibility of the polarized interferometer concept are verified by simulation.

Raab, F. H.↗
NASA NTRS1974

Identification of aircraft stability and control derivatives in the presence of turbulence

A maximum likelihood estimator for a linear system with state and observation noise is developed to determine stability and control derivatives from flight data obtained in the presence of turbulence. The formulation for the longitudinal short-period mode is presented briefly, including a special case that greatly simplifies the problem if the measurement noise on one signal is negligible. The effectiveness and accuracy of the technique are assessed by applying it first to simulated flight data, in which the true parameter values and state noise are known, then to actual flight data obtained in turbulence. The results are compared with data obtained in smooth air and with wind-tunnel data. The complete maximum likelihood estimator, which accounts for both state and observation noise, is shown to give the most accurate estimate of the stability and control derivatives from flight data obtained in turbulence. It is superior to the techniques that ignores state noise and to the simplified method that neglects the measurement noise on the angle-of-attack signal.

Iliff, K. W.↗
NASA NTRS1991

Characterization of frequency standard instability by estimation of their covariance matrix

The popular 3-cornered hat method used for evaluating the noise contributions of individual frequency standards is revisited. This method is used in several cases, but sometimes the results are not consistent because one or more estimated clock variances turn out to be negative. Different causes of this unacceptable result have been conjectured: among them one regards the hypothesis of uncorrelated clocks, essential in this method. Since recently realistic cases of correlation between clocks, mainly due to the environmental conditions, have been observed, this paper proposes an entirely revisited version of the 3-cornered hat method which permits to evaluate the individual variances and also the possible covariances between clocks, by relaxing the hypothesis of uncorrelation. The uncertainty and the lack of contemporaneity of the measurement series are assumed to be negligible. The lack of the uncorrelation hypothesis calls for a more general mathematical model leading to an underdetermined linear system. The estimates of the (co)variances of the measurement series us well us those of the individual clocks are introduced by means of the scalar product of the related time series and arranged in the respective covariance matrices S and R. Since covariance matrix is positive definite by definition, the problem consists in estimating the unknown R, subject to the constraint of positive definiteness, from the known S. Unfortunately, this constraint is not sufficient to estimate R. Therefore a suitable optimization criterion is proposed, which assures the positive definiteness of R and, at the same time, minimizes the global correlation among clocks. Examples of frequency instability measurements processed by the "classical" 3-cornered hat method and the here-revisited method are presented showing that the solutions are identical only when the uncorrelation hypothesis doesn't violate the positive definiteness of R.

PatriziaTravella↗
DOE OSTI2018

Acute HBV infection in humanized chimeric mice has multiphasic viral kinetics

Background: Chimeric uPA/SCID mice reconstituted with humanized livers are useful for studying HBV infection in the absence of an adaptive immune response. However, the detailed characterization of HBV infection kinetics necessary to enable in-depth mechanistic studies in this novel in vivo HBV infection model is lacking. Methods: To characterize HBV kinetics post-inoculation (p.i.) to steady state, 42 mice were inoculated with HBV. Serum HBV DNA was frequently measured from 1 minute to 63 days p.i. Total intrahepatic HBV DNA, HBV cccDNA, and HBV RNA was measured in a subset of mice at 2, 4, 6, 10, and 13 weeks p.i. HBV half-life (t 1/2 ) was estimated using a linear mixed-effects model. Results: During the first 6 h p.i. serum HBV declined in repopulated uPA/SCID mice with a t 1/2 =62 min [95%CI=59-67min]. Thereafter, viral decline slowed followed by a 2 day lower plateau. Subsequent viral amplification was multiphasic with an initial mean doubling time of t 2 = 8±3 h followed by an interim plateau before prolonged amplification (t 2 =2±0.5 days) to a final HBV steady state of 9.3 ± 0.3 log copies/ml. Serum HBV and intrahepatic HBV DNA were positively correlated (R2=0.98). Conclusions: HBV infection in uPA/SCID chimeric mice is highly dynamic despite the absence of an adaptive immune response. The serum HBV t 1/2 in humanized uPA/SCID mice was estimated to be ~1 h regardless of inoculum size. Finally, the HBV acute infection kinetics presented here is an important step in characterizing this experimental model system so that it can be effectively used to elucidate the dynamics of the HBV lifecycle and thus possibly reveal effective antiviral drug targets.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1992

Dynamically stratified Monte Carlo forecasting

A new method for performing Monte Carlo forecasts is introduced. The method, called dynamic stratification, selects initial perturbations based on a stratification of the error distribution. A simple implementation is presented in which the error distribution used for the stratification is estimated from a linear model derived from a large ensemble of 12-h forecasts with the full dynamic model. The stratification thus obtained is used to choose a small subsample of initial states with which to perform the dynamical Monte Carlo forecasts. Several test cases are studied using a simple two-level general circulation model with uncertain initial conditions. It is found that the method provides substantial reductions in the sampling error of the forecast mean and variance when compared to the more traditional approach of choosing the initial perturbations at random. The degree of improvement, however, is sensitive to the nature of the initial error distribution and to the base state. In practice the method may be viable only if the computational burden involved in obtaining an adequate estimate of the error distribution is shared with the data-assimilation procedure.

Schubert, Siegfried↗
DOE OSTI2025

Simultaneous Observation of Ion-scale Wave Packets with Opposite Polarizations and Their Implications on the Generation Region in the Inner Heliosphere

This paper reports a dispersion analysis of two wave packets simultaneously observed near the local proton gyrofrequency by the Parker Solar Probe. The observed wave event exhibits clear two-banded wave packets both propagating along the magnetic field, characterized by left-handed (L-mode) and right-handed (R-mode) polarizations simultaneously. By incorporating the Doppler shift effect into a linear dispersion analysis, we find two possible scenarios that explain these simultaneous opposite polarizations: (1) Two inherently L-mode waves in the plasma frame, propagate parallel and antiparallel to the solar wind velocity, with similar wave frequencies and wave numbers. The polarization of the antiparallel propagating wave reverses as it moves sunward in the plasma frame while still comoving with the solar wind in the stationary frame. This reversal manifests the polarization of the wave as an R-mode in the spacecraft frame. (2) Simultaneous L-mode and R-mode waves propagate parallel to the solar wind velocity, with different wave frequencies and wave numbers. Concurrent proton observations during the wave event reveal a dominant anisotropic ($T$⟂/$T$ ∥ > 1) core distribution with a drifting beam population. Estimation of the linear growth rate for both L-mode and R-mode waves suggests that both scenarios are plausible, indicating that the observation is near the wave-generation region. We explore the potential impact of these simultaneous waves on solar wind heating and scattering effects, hypothesizing that such waves might enhance efficiency compared to waves with a single wave packet, contingent upon the statistical significance of such waves.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Monitoring of Temperature Measurements for Different Flow Regimes in Water and Galinstan with Long Short-Term Memory Networks and Transfer Learning of Sensors

Temperature sensing is one of the most common measurements of a nuclear reactor monitoring system. The coolant fluid flow in a reactor core depends on the reactor power state. We investigated the monitoring and estimation of the thermocouple time series using machine learning for a range of flow regimes. Measurement data were obtained, in two separate experiments, in a flow loop filled with water and with liquid metal Galinstan. We developed long short-term memory (LSTM) recurrent neural networks (RNNs) for sensor predictions by training on the sensor’s own prior history, and transfer learning LSTM (TL-LSTM) by training on a correlated sensor’s prior history. Sensor cross-correlations were identified by calculating the Pearson correlation coefficient of the time series. The accuracy of LSTM and TL-LSTM predictions of temperature was studied as a function of Reynolds number (Re). The root-mean-square error (RMSE) for the test segment of time series of each sensor was shown to linearly increase with Re for both water and Galinstan fluids. Using linear correlations, we estimated the range of values of Re for which RMSE is smaller than the thermocouple measurement uncertainty. For both water and Galinstan fluids, we showed that both LSTM and TL-LSTM provide reliable estimations of temperature for typical flow regimes in a nuclear reactor. The LSTM runtime was shown to be substantially smaller than the data acquisition rate, which allows for performing estimation and validation of sensor measurements in real time.

Pantopoulou, Stella↗
DOE OSTI2024

Point Sensor Networks Struggle to Detect and Quantify Short Controlled Releases at Oil and Gas Sites

This study evaluated multiple commercially available continuous monitoring (CM) point sensor network (PSN) solutions under single-blind controlled release testing conducted at operational upstream and midstream oil and natural gas (O&G) sites. During releases, PSNs reported site-level emission rate estimates of 0 kg/h between 38 and 86% of the time. When non-zero site-level emission rate estimates were provided, no linear correlation between the release rate and the reported emission rate estimate was observed. The average, aggregated across all PSN solutions during releases, shows 5% of the mixing ratio readings at downwind sensors were greater than the site’s baseline plus two standard deviations. Four of seven total PSN solutions tested during this field campaign provided site-level emission rate estimates with the site average relative error ranging from −100% to 24% for solution D, −100% to −43% for solution E, −25% for solution F (solution F was only at one site), and −99% to 430% for solution G, with an overall average of −29% across all sites and solutions. Of all the individual site-level emission rate estimates, only 11% were within ±2.5 kg/h of the study team’s best estimate of site-level emissions at the time of the releases.

42 ENGINEERING↗
NASA NTRS1974

The structure and circulation of the deep Venus atmosphere

A simple model for the structure of a nonrotating Hadley regime in an atmosphere with large thermal inertia is developed. The radiative fluxes are estimated by using a linearization about the radiative equilibrium state, and the dynamical fluxes are estimated by using scaling analysis. The requirement that differential heating by these fluxes be in balance in both the meridional and vertical directions leads to two equations for the mean static stability and meridional temperature contrast. The solution depends on two parameters: the strength of the radiative heating, as measured by the static stability of the radiative equilibrium state; and the ratio of the time it takes an external gravity wave to traverse the atmosphere to the time it would take the atmosphere to cool off radiatively. It is shown that it is not necessary to invoke convection to explain the approximate adiabatic lapse rate in the Venus atmosphere, but a greenhouse effect is necessary to explain the high surface temperatures.

Stone, P. H.↗
NASA NTRS1992

Tidal capture of stars by a massive black hole

The processes leading to tidal capture of stars by a massive black hole and the consequences of these processes in a dense stellar cluster are discussed in detail. When the amplitude of a tide and the subsequent oscillations are sufficiently large, the energy deposited in a star after periastron passage and formation of a bound orbit cannot be estimated directly using the linear theory of oscillations of a spherical star, but rather numerical estimates must be used. The evolution of a star after tidal capture is discussed. The maximum ratio R of the cross-section for tidal capture to that for tidal disruption is about 3 for real systems. For the case of a stellar system with an empty capture loss cone, even in the case when the impact parameter for tidal capture only slightly exceeds the impact parameter for direct tidal disruption, tidal capture would be much more important than tidal disruption.

Novikov, I. D.↗
NASA NTRS2016

Seasonal Monitoring and Estimation of Regional Aerosol Distribution over Po Valley, Northern Italy, Using a High-Resolution MAIAC Product

In this work, the new 1-km-resolved Multi-Angle Implementation of Atmospheric Correction (MAIAC) algorithm is employed to characterize seasonal AOD-PM10 correlations over northern Italy. The accuracy of the new dataset is assessed versus the widely used Moderate Resolution Imaging Spectroradiometer (MODIS) Collection 5.1 Aerosol Optical Depth (AOD) data, retrieved at 0.55 microns with spatial resolution of 10 km (MYD04). We focused on evaluating the ability of these two products to characterize both temporal and spatial distributions of aerosols within urban and suburban areas. Ground PM10 measurements were obtained from 73 of the Italian Regional Agency for Environmental Protection (ARPA) monitoring stations, spread across northern Italy, for a three-year period from 2010 to 2012. The Po Valley area (northern Italy) was chosen as the study domain because of severe urban air pollution, resulting from the highest population and industrial manufacturing density in the country, being located in a valley where two surrounding mountain chains favor the stagnation of pollutants. We found that the global correlations between PM10 and AOD are R(sup 2) = 0.83 and R(sup 2) = 0.44 for MYD04_L2 and for MAIAC, respectively, suggesting for a greater sensitiveness of the high-resolution product to small-scale deviations. However, the introduction of Relative Humidity (RH) and Planetary Boundary Layer (PBL) depth corrections gave a significant improvement to the PM AOD correlation, which led to similar performance: R(sup 2) = 0.96 for MODIS and R(sup 2) = 0.95 for MAIAC. Furthermore, the introduction of the PBL information in the corrected AOD values was found to be crucial in order to capture the clear seasonal cycle shown by measured PM10 values. The study allowed us to define four seasonal linear correlations that estimate PM10 concentrations satisfactorily from the remotely sensed MAIAC AOD retrieval. Overall, the results show that the high resolution provided by MAIAC retrieval data is much more relevant than 10km MODIS data to characterize PM10 in this region of Italy which has a pretty limited geographical domain, but a broad variety of land usages and consequent particulate concentrations.

MODIS↗
NASA NTRS1997

Global Geodetic Signatures of the Antarctic Ice Sheet

Four scenarios of present day Antarctic ice sheet mass change are developed from comprehensive reviews of the available glaciological and oceanographic evidence. The gridded scenarios predict widely varying contributions to secular sea level change xi ranging from -1.1 to 0.45 mm/yr, and predict polar motion m and time-varying low-degree gravitational coefficients J(sub l) that differ significantly from earlier estimates. A reasonably linear relationship between the rate of sea level change from Antarctica xi(sub A) and the predicted Antarctic J(sub l) is found for the four scenarios. This linearity permits a series of forward models to be constructed that incorporate the effects of ice mass changes in Antarctica, Greenland, and distributed smaller glaciers, as well as postglacial rebound (assuming the ICE-3G deglaciation history), with the goal of obtaining optimum reconciliation between observed constraints on J(sub l) and sea level rise xi. Numerous viable combinations of lower mantle viscosity and hydrologic sources are found that safely "observed" in the range of 1 to 2-2.5 mm/yr and observed J(sub l) for degrees 2, 3, and 4. In contrast, rates of global sea level rise above 2.5 mm/yr are inconsistent with available J(sub l) observations. The successful composite models feature a pair of lower mantle viscosity solutions arising from the sensitivity of J(sub l) to glacial rebound. The paired values are well separated at xi = 1 mm/yr, but move closer together as xi is increased, and, in fact, merge around xi = 2 - 2.5 mm/yr, revealing an intimate relation between xi and preferred lower mantle viscosity. This general pattern is quite robust and persists for different J(sub l) solutions, for variations in source assumptions, and for different styles of lower mantle viscosity stratification. Tighter J(sub l) constraints for l greater than 2 may allow some viscosity stratification schemes and source assumptions to be excluded in the future. For a given total observed xi the sea level rise from Antarctica xi(sub A) is tightly constrained and ranges from 0 to + 1 mr/yr (corresponding to an ablating ice sheet) as estimates of xi are raised from 1 to 2.5 mm/yr. However, when the degree 3 zonal harmonic constraint is removed, the solutions show little sensitivity to Antarctic mass balance, emphasizing the need for a well determined odd-degree secular zonal harmonic for determining polar ice mass balance.

James, Thomas S.↗
NASA NTRS2015

Cloud-to-Ground Lightning Estimates Derived from SSMI Microwave Remote Sensing and NLDN

Lightning observations are collected using ground-based and satellite-based sensors. The National Lightning Detection Network (NLDN) in the United States uses multiple ground sensors to triangulate the electromagnetic signals created when lightning strikes the Earth's surface. Satellite-based lightning observations have been made from 1998 to present using the Lightning Imaging Sensor (LIS) on the NASA Tropical Rainfall Measuring Mission (TRMM) satellite, and from 1995 to 2000 using the Optical Transient Detector (OTD) on the Microlab-1 satellite. Both LIS and OTD are staring imagers that detect lightning as momentary changes in an optical scene. Passive microwave remote sensing (85 and 37 GHz brightness temperatures) from the TRMM Microwave Imager (TMI) has also been used to quantify characteristics of thunderstorms related to lightning. Each lightning detection system has fundamental limitations. TRMM satellite coverage is limited to the tropics and subtropics between 38 deg N and 38 deg S, so lightning at the higher latitudes of the northern and southern hemispheres is not observed. The detection efficiency of NLDN sensors exceeds 95%, but the sensors are only located in the USA. Even if data from other ground-based lightning sensors (World Wide Lightning Location Network, the European Cooperation for Lightning Detection, and Canadian Lightning Detection Network) were combined with TRMM and NLDN, there would be enormous spatial gaps in present-day coverage of lightning. In addition, a globally-complete time history of observed lightning activity is currently not available either, with network coverage and detection efficiencies varying through the years. Previous research using the TRMM LIS and Microwave Imager (TMI) showed that there is a statistically significant correlation between lightning flash rates and passive microwave brightness temperatures. The physical basis for this correlation emerges because lightning in a thunderstorm occurs where ice is first present in the cloud and electric charge separation occurs. These ice particles efficiently scatter the microwave radiation at the 85 and 37 GHz frequencies, thus leading to large brightness temperature depressions. Lightning flash rate is related to the total amount of ice passing through the convective updraft regions of thunderstorms. Confirmation of this relationship using TRMM LIS and TMI data, however, remains constrained to TRMM observational limits of the tropics and subtropics. Satellites from the Defense Meteorology Satellite Program (DMSP) have global coverage and are equipped with passive microwave imagers that, like TMI, observe brightness temperatures at 85 and 37 GHz. Unlike the TRMM satellite, however, DMSP satellites do not have a lightning sensor, and the DMSP microwave data has never been used to derive global lightning. In this presentation, a relationship between DMSP Special Sensor Microwave Imager (SSMI) data and ground-based cloud-to-ground (CG) lightning data from NLDN is investigated to derive a spatially complete time history of CG lightning for the USA study area. This relationship is analogous to the established using TRMM LIS and TMI data. NLDN has the most spatially and temporally complete CG lightning data for the USA, and therefore provides the best opportunity to find geospatially coincident observations with SSMI sensors. The strongest thunderstorms generally have minimum 85 GHz Polarized Corrected brightness Temperatures (PCT) less than 150 K. Archived radar data was used to resolve the spatial extent of the individual storms. NLDN data for that storm spatial extent defined by radar data was used to calculate the CG flash rate for the storm. Similar to results using TRMM sensors, a linear model best explained the relationship between storm-specific CG flash rates and minimum 85 GHz PCT. However, the results in this study apply only to CG lightning. To extend the results to weaker storms, the probability of CG lightning (instead of the flash rate) was calculated for storms having 85 GHz PCT greater than 150 K. NLDN data was used to determine if a CG strike occurred for a storm. This probability of CG lightning was plotted as a function of minimum 85 GHz PCT and minimum 37 GHz PCT. These probabilities were used in conjunction with the linear model to estimate the CG flash rate for weaker storms with minimum 85 GHz PCTs greater than 150 K. Results from the investigation of CG lightning and passive microwave radiation signals agree with the previous research investigating total lightning and brightness temperature. Future work will take the established relationships and apply them to the decades of available DMSP data for the USA to derive a map of CG lightning flash rates. Validation of this method and uncertainty analysis will be done by comparing the derived maps of CG lightning flash rates against existing NLDN maps of CG lightning flash rates.

Winesett, Thomas↗
NASA NTRS1985

Wavelength error analysis in a multiple-beam Fizeau laser wavemeter having a linear diode array readout

An estimate of the wavelength accuracy of a laser wavemeter is performed for a system consisting of a multiple-beam Fizeau interferometer and a linear photosensor array readout. The analysis consists of determining the fringe position errors which result when various noise sources are included in the fringe forming and detection process. Two methods of estimating the fringe centers are considered: (1) maximum pixel current location, and (2) average pixel location for two detectors with nearly equal output currents. Wavelength error results for these two methods are compared for some typical wavemeter parameters.

Robinson, D. M.↗
DOE OSTI2023

Granger Causality for prediction in Dynamic Mode Decomposition: Application to power systems

Here, the dynamic mode decomposition (DMD) technique extracts the dominant modes characterizing the innate dynamical behavior of the system within the measurement data. For appropriate identification of dominant modes from the measurement data, the DMD algorithm necessitates ensuring the quality of the input measurement data sequences. On that account, for validating the usability of the dataset for the DMD algorithm, the paper proposed two conditions: Persistence of excitation (PE) and the Granger Causality Test (GCT). The virtual data sequences are designed with the hankel matrix representation such that the dimensions of the subspace spanning the essential system modes are increased with the addition of new state variables. The PE condition provides the lower bound for the trajectory length, and the GCT provides the order of the model. Satisfying the PE condition enables estimating an approximate linear model, but the predictability with the identified model is only assured with the temporal causation among data searched with GCT. The proposed methodology is validated with the application for coherency identification (CI) in a multi-machine power system (MMPS), an essential phenomenon in transient stability analysis. The significance of PE condition and GCT is demonstrated through various case studies implemented on 22 bus six generator system.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Generalized fiducial inference on differentiable manifolds

We introduce a novel approach to inference on parameters that take values in a Riemannian manifold embedded in a Euclidean space. Parameter spaces of this form are ubiquitous across many fields, including chemistry, physics, computer graphics, and geology. Here, this new approach uses generalized fiducial inference (GFI) to obtain a posterior-like distribution on the manifold, without needing to know local parameterizations that map to the constrained space from an unconstrained Euclidean space. Using mathematical tools from Riemannian geometry, we construct a constrained generalized fiducial distribution (CGFD). A Bernstein-von Mises-type result for the CGFD, which provides intuition for how the desirable asymptotic qualities of the unconstrained generalized fiducial distribution are inherited by the CGFD, is provided. To illustrate the practical use of the CGFD, we provide a proof-of-concept example in the context of a linear logspline density estimation problem, and demonstrate that CGFD-based confidence sets exhibit desirable coverage properties via simulation. As an application, we fit a CGFD to COVID-19 case count data from North Carolina, USA.

97 MATHEMATICS AND COMPUTING↗
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↗