6.4.4 — Earthquake Response — maturity: live

Aftershock Risk Modelling

Using satellite-derived surface deformation and stress-transfer maps to forecast the location, magnitude and timing of aftershocks following a major earthquake.

Auto-drafted reference page — content reviewed for shape, individual references not yet link-checked.

After a major earthquake, the question civil protection agencies cannot answer fast enough is: where will the next damaging shock hit, and how hard? Conventional seismic networks are ground-based, sparse in developing nations, and tell you what happened, not what is coming. Aftershock sequences are governed by stress redistribution across the fault system — information that is encoded in the centimetre-scale surface deformation field that InSAR satellites capture within hours of the main event.

A sovereign constellation pairing C-band or L-band SAR with a ground-truth seismic telemetry feed can generate Coulomb stress-transfer maps in near-real time, then feed them into operational aftershock forecasting models (ETAS, Coulomb3, OEF-based frameworks). The stack ingests deformation data, resolves the fault rupture geometry, calculates stress increments on optimally oriented receiver faults, and probabilistically forecasts M≥5 aftershock rates over 24-hour to 30-day windows. Revisit cadence is the decisive variable: a 12-satellite walker constellation at 520 km provides 6-hour revisit over seismically active corridors, letting the model update after every significant aftershock.

The operational payoff is measurable. Emergency managers get a probabilistic hazard map, updated every six hours, showing which districts remain under elevated risk. Search-and-rescue teams can be redeployed away from areas due a M6+ aftershock before it strikes. Engineers inspecting nominally standing buildings get a ranked list of sites where ground shaking is statistically most likely to recur. This is not academic seismology — it is the difference between a government that manages the disaster sequence and one that is perpetually surprised by it.

What matters

Aftershock sequences decay as a power law (Omori–Utsu), but stress-transfer anomalies produce off-fault clusters that pure statistical models miss without satellite deformation input.
A six-hour InSAR revisit over the rupture zone is sufficient to update Coulomb stress maps after each M≥4.5 aftershock and revise the 24-hour probabilistic hazard forecast.
Nations without sovereign SAR access waited 48–72 hours for commercial or allied imagery after the 2023 Kahramanmaraş sequence — enough delay to compromise the first critical search-and-rescue window.
Aftershock forecast maps carry legal weight: evacuation orders and building reoccupancy decisions made against them must be traceable to sovereign-controlled, auditable data chains.

Sovereignty score: 9/10 — Aftershock forecasts govern evacuation orders, reoccupancy permits and rescue redeployment decisions — authorities that cannot be delegated to a foreign satellite operator's tasking queue.

Tasking priority: during a major disaster, commercial and allied SAR constellations are overwhelmed with competing requests; sovereign satellites are legally and operationally directed by national civil protection authorities without negotiation or embargo risk.
Data latency is life-critical: a 6-hour InSAR revisit updates the probabilistic hazard map before the next significant aftershock; dependence on foreign tasking routinely introduces 48-72 hour gaps that void the forecast's operational value.
Legal accountability: evacuation and reoccupancy decisions derived from aftershock forecasts carry administrative and criminal liability under national disaster law, requiring an auditable, sovereign-controlled data chain that foreign commercial providers cannot guarantee.
Geopolitical exposure: seismically active nations in contested regions (Turkey, Iran, Pakistan, Central Asia) cannot assume uninterrupted access to US- or EU-controlled SAR assets during periods of diplomatic tension coinciding with seismic events.

Reference architecture

Payload: C-band SAR, 5 m stripmap resolution, 80 km swath, interferometric mode (IW/TOPS); optional L-band secondary payload for deeper soil penetration and improved coherence over vegetated terrain
Bus class: ESPA-class microsat, 160–200 kg wet mass, 900 W payload power, 3-axis stabilised to <0.05° pointing accuracy for repeat-pass interferometric coherence
Orbit: Sun-synchronous LEO at 510–530 km, 12-satellite walker constellation, ascending and descending geometry, 6-hour revisit over active seismic corridors between 60°N and 60°S
Ground segment: 3-station national network (X-band downlink, S-band TT&C) co-located with national seismological institute; SatNOGS UHF/VHF backup for housekeeping; direct downlink to forward-deployed mobile X-band terminal for disaster-zone operations
Software & data
Latency
Cost band
Lead time

References

Operational Earthquake Forecasting: State of Knowledge and Guidelines for Utilization — USGS — USGS Open-File Report 2013-1138 defines the scientific basis and operational protocols for short-term aftershock forecasting, including Coulomb stress transfer and ETAS models. It is the de facto reference framework for civil protection agencies implementing operational earthquake forecasting.
Sentinel-1 SAR for Earthquake Rapid Mapping — ESA Copernicus — ESA's Sentinel-1 mission routinely acquires interferometric SAR over earthquake zones within hours of major events, providing the surface deformation data that underpins fault-rupture inversion and Coulomb stress modelling. Access is free but controlled by ESA operational priorities, not national civil protection mandates.
Coulomb 3.4 Stress Transfer Software — USGS — The USGS Coulomb software calculates static stress changes on specified receiver faults following an earthquake rupture, widely used to identify fault segments brought closer to failure by the main shock. Integration with InSAR-derived slip models is the standard operational workflow.
ETAS Aftershock Forecasting — GFZ Helmholtz Centre Potsdam — GFZ operates real-time Epidemic-Type Aftershock Sequence (ETAS) models as part of European seismic hazard services, demonstrating that probabilistic aftershock rate forecasts are operationally viable when fed with rapid source parameters. Their work underpins the link between geodetic data ingestion and forecast update cycles.
The 2023 Kahramanmaraş Earthquake Sequence — USGS Earthquake Hazards Program — USGS documentation of the February 2023 M7.8–M7.7 doublet illustrates the complexity of multi-segment rupture sequences that generate spatially distributed aftershock fields extending hundreds of kilometres. The event underscored the operational need for rapid, high-revisit satellite geodetic coverage to constrain aftershock hazard across an unexpectedly wide region.