Satellite InSAR Surface Deformation and Coulomb Stress Transfer Modeling across Active Tectonic Fault Systems

By TechIDaily Earth Observation & Geophysical Computing Group · Published 2026-10-20

Active continental deformation and earthquake hazards require sub-centimeter observation of ground surface kinematics over continental scales. Historical geodetic networks relied on sparse permanent GNSS stations, leaving vast fault systems poorly sampled between sparse benchmarks.

Spaceborne Interferometric Synthetic Aperture Radar (InSAR) has revolutionized crustal kinematics by converting phase differences between repeat radar passes into continuous line-of-sight (LOS) displacement maps with millimeter-level precision.

1. Interferometric Phase Decomposition

The raw interferometric phase $Delta phi$ across two SAR acquisitions separated by baseline $B_perp$ is given by:

$$Delta phi = phi_{ ext{def}} + phi_{ ext{topo}} + phi_{ ext{atm}} + phi_{ ext{orbit}} + phi_{ ext{noise}}$$

To isolate crustal deformation $phi_{ ext{def}}$, atmospheric phase screens ($phi_{ ext{atm}}$) must be estimated using numerical weather models (ERA5) and spatial-temporal filtering:

$$d_{ ext{LOS}} = rac{lambda}{4pi} Delta phi_{ ext{unwrapped}}$$

Where $lambda$ is the radar wavelength (e.g., $lambda = 5.6, ext{cm}$ for C-band Sentinel-1; $lambda = 24, ext{cm}$ for L-band NISAR).

2. Coulomb Failure Stress ($Delta ext{CFS}$) Modeling

Static stress transfer from seismic rupture or interseismic creep alters failure conditions on neighboring fault segments:

$$Delta ext{CFS} = Delta au + mu (Delta sigma_n + Delta p)$$

Where $Delta au$ is shear stress change along the slip direction, $Delta sigma_n$ is normal stress change (positive for unclamping), $mu$ is the apparent friction coefficient, and $Delta p$ is pore fluid pressure change.

Positive $Delta ext{CFS}$ values exceeding $0.01, ext{MPa}$ (0.1 bar) are widely established to trigger subsequent rupture events on critically stressed adjacent faults.

3. Geophysical Computing Verdict

The integration of Sentinel-1 constellation archives, the NASA-ISRO SAR (NISAR) radar platform, and cloud-native interferogram processing pipelines (ISCE3) enables automated near-real-time slip inversion and Coulomb stress forecasting, transforming earthquake hazard assessment into a continuous geodetic observation science.