Pietro Milillo

Geodetic Constraints on the 2014 M 6.0 South Napa Earthquake

On 24 August 2014, the M 6.0 South Napa earthquake shook much of the San Francisco Bay area, leading to significant damage in the Napa Valley. The earthquake occurred in the vicinity of the West Napa fault (122.313° W, 38.22° N, 11.3 km), a mapped structure located between the Rodger’s Creek and Green Valley faults, with nearly pure right‐lateral strike‐slip motion (strike 157°, dip 77°, rake –169°; http://comcat.cr.usgs.gov/earthquakes/eventpage/nc72282711#summary, last accessed December 2014) (Fig. 1). The West Napa fault previously experienced an M 5 strike‐slip event in 2000 but otherwise exhibited no previous definitive evidence of historic earthquake rupture (Rodgers et al., 2008; Wesling and Hanson, 2008). Evans et al. (2012) found slip rates of ∼9.5  mm/yr along the West Napa fault, with most slip rate models for the Bay area placing higher slip rates and greater earthquake potential on the Rodger’s Creek and Green Valley faults, respectively (e.g., Savage et al., 1999; d’Alessio et al., 2005; Funning et al., 2007).

Coseismic deformation and triggered landslides of the 2016 Mw 6.2 Amatrice earthquake in Italy

The Central Apennines in Italy have had multiple moderate-size but damaging shallow earthquakes. In this study, we optimize the fault geometry and invert for fault slip based on coseismic GPS and interferometric synthetic aperture radar (InSAR) for the 2016 Mw 6.2 Amatrice earthquake in Italy. Our results show that nearly all the fault slip occurred between 3 and 6 km depth but extends 20 km along strike. There was less than 4 cm static surface displacement at the town Amatrice where the most devastating damage occurred. Landslides triggered by earthquake ground shaking are not uncommon, but triggered landslides with submeter movement are challenging to be observed in the field. We find evidence of coseismically triggered deep-seated landslides northwest and northeast of the epicenter where coseismic peak ground acceleration was estimated >0.5 g. By combining ascending and descending InSAR data, we are able to estimate the landslide thickness as at least 100 and 80 m near Monte Vettore and west of Castelluccio, respectively. The landslide near Monte Vettore terminates on the preexisting fault Monte Vettore Fault (MVEF) scarp. Our results imply that the long-term fault slip rate of MVEF estimated based on paleoseismic studies could potentially have errors due to triggered landslides from nearby earthquake events.

COSMO-SkyMed Spotlight Interferometry Over Rural Areas: The Slumgullion Landslide in Colorado, USA

In the last 7 years, spaceborne synthetic aperture radar (SAR) data with resolution of better than a meter acquired by satellites in spotlight mode offered an unprecedented improvement in SAR interferometry (InSAR). Most attention has been focused on monitoring urban areas and man-made infrastructure exploiting geometric accuracy, stability, and phase fidelity of the spotlight mode. In this paper, we explore the potential application of the COSMO-SkyMed® Spotlight mode to rural areas where decorrelation is substantial and rapidly increases with time. We focus on the rapid repeat times of as short as one day possible with the COSMO-SkyMed® constellation. We further present a qualitative analysis of spotlight interferometry over the Slumgullion landslide in southwest Colorado, which moves at rates of more than 1 cm/day.

An aseismic slip transient on the North Anatolian Fault

Constellations of Synthetic Aperture Radar (SAR) satellites with short repeat time acquisitions allow exploration of active faults behavior with unprecedented temporal resolution. Along the North Anatolian Fault (NAF) in Turkey, an 80 km long section has been creeping at least since the 1944, M_w 7.3 earthquake near Ismetpasa, with a current Interferometric Synthetic Aperture Radar (InSAR)-derived average creep rate of 8 ± 3 mm/yr (i.e., a third of the NAF long-term slip rate). We use a dense set of SAR images acquired by the COSMO-SkyMed constellation to quantify the spatial distribution and temporal evolution of creep over 1 year. We identify a major burst of aseismic slip spanning 31 days with a maximum slip of 2 cm, between the surface and 4 km depth. This result shows that fault creep along this section of the NAF does not occur at a steady rate as previously thought, highlighting a need to revise our understanding of the underlying fault mechanics.

On the Synergistic Use of SAR Constellations’ Data Exploitation for Earth Science and Natural Hazard Response

Several current and expected future SAR satellites missions (e.g., TanDEM-X (TDX)/PAZ, COSMO-SkyMed (CSK), and Sentinel-1A/B) are designed as constellations of SAR sensors. Relative to single satellite systems, such constellations can provide greater spatial coverage and temporal sampling, thereby enabling better control on interferometric decorrelation and lower latency data access. These improvements lead to more effective near real-time disaster monitoring, assessment and response, and a greater ability to constrain dynamically changing physical processes. Using observations from the CSK system, we highlight examples of the potential for such imaging capabilities to enable advances in Earth science and natural hazards response.

Rapid Damage Mapping for the 2015 Mw 7.8 Gorkha Earthquake Using Synthetic Aperture Radar Data from COSMO–SkyMed and ALOS-2 Satellites

The 25 April 2015 M_w 7.8 Gorkha earthquake caused more than 8000 fatalities and widespread building damage in central Nepal. The Italian Space Agency’s COSMO–SkyMed Synthetic Aperture Radar (SAR) satellite acquired data over Kathmandu area four days after the earthquake and the Japan Aerospace Exploration Agency’s Advanced Land Observing Satellite-2 SAR satellite for larger area nine days after the mainshock. We used these radar observations and rapidly produced damage proxy maps (DPMs) derived from temporal changes in Interferometric SAR coherence. Our DPMs were qualitatively validated through comparison with independent damage analyses by the National Geospatial-Intelligence Agency and the United Nations Institute for Training and Research’s United Nations Operational Satellite Applications Programme, and based on our own visual inspection of DigitalGlobe’s WorldView optical pre- versus postevent imagery. Our maps were quickly released to responding agencies and the public, and used for damage assessment, determining inspection/imaging priorities, and reconnaissance fieldwork.

Monitoring dam structural health from space: Insights from novel InSAR techniques and multi-parametric modeling applied to the Pertusillo dam Basilicata, Italy

Abstract The availability of new constellations of synthetic aperture radar (SAR) sensors is leading to important advances in infrastructure monitoring. These constellations offer the advantage of reduced revisit times, providing low-latency data that enable analysis that can identify infrastructure instability and dynamic deformation processes. In this paper we use COSMO-SkyMed (CSK) and TerraSAR-X (TSX) data to monitor seasonal induced deformation at the Pertusillo dam (Basilicata, Italy) using multi-temporal SAR data analysis. We analyzed 198 images spanning 2010–2015 using a coherent and incoherent PS approach to merge COSMO-SkyMed adjacent tracks and TerraSAR-X acquisitions, respectively. We used hydrostatic-seasonal-temporal (HST) and hydrostatic-temperature-temporal (HTT) models to interpret the non-linear deformation at the dam wall using ground measurements together with SAR time-series analysis. Different look geometries allowed us to characterize the horizontal deformation field typically observed at dams. Within the limits of our models and the SAR acquisition sampling we found that most of the deformation at the Pertusillo dam can be explained by taking into account only thermal seasonal dilation and hydrostatic pressure. The different models show slightly different results when interpreting the aging term at the dam wall. The results highlight how short-revisit SAR satellites in combination with models widely used in the literature for interpreting pendulum and GPS data can be used for supporting structural health monitoring and provide valuable information to ground users directly involved in field measurements.

Space geodetic monitoring of engineered structures: The ongoing destabilization of the Mosul dam, Iraq

We present a detailed survey of the ongoing destabilization process of the Mosul dam. The dam is located on the Tigris river and is the biggest hydraulic structure in Iraq. From a geological point of view the dam foundation is poor due to a site geology formed by alternating strata of highly soluble materials including gypsum, anhydrite, marl and limestone. Here we present the first multi-sensor cumulative deformation map for the dam generated from space-based interferometric synthetic aperture radar measurements from the Italian constellation COSMO-SkyMed and the European sensor Sentinel-1a over the period 2014–2016 that we compare to an older dataset spanning 2004–2010 acquired with the European Envisat satellite. We found that deformation was rapid during 2004–2010, slowed in 2012–2014 and increased since August 2014 when grouting operations stopped due to the temporary capture of the dam by the self proclaimed Islamic State. We model the inferred deformation using a Markov chain Monte Carlo approach to solve for change in volume for simple tensile dislocations. Results from recent and historical geodetic datasets suggests that the volume dissolution rate remains constant when the equivalent volume of total concrete injected during re-grouting operations is included in the calculations.

Tidally induced variations in vertical and horizontal motion on Rutford Ice Stream, West Antarctica, inferred from remotely sensed observations

To better understand the influence of stress changes over floating ice shelves on grounded ice streams, we develop a Bayesian method for inferring time-dependent 3-D surface velocity fields from synthetic aperture radar (SAR) and optical remote sensing data. Our specific goal is to observe ocean tide-induced variability in vertical ice shelf position and horizontal ice stream flow. Thus, we consider the special case where observed surface displacement at a given location can be defined by a 3-D secular velocity vector, a family of 3-D sinusoidal functions, and a correction to the digital elevation model used to process the SAR data. Using nearly 9 months of SAR data collected from multiple satellite viewing geometries with the COSMO-SkyMed 4-satellite constellation, we infer the spatiotemporal response of Rutford Ice Stream, West Antarctica, to ocean tidal forcing. Consistent with expected tidal uplift, inferred vertical motion over the ice shelf is dominated by semidiurnal and diurnal tidal constituents. Horizontal ice flow variability, on the other hand, occurs primarily at the fortnightly spring-neap tidal period (M_(sf)). We propose that periodic grounding of the ice shelf is the primary mechanism for translating vertical tidal motion into horizontal flow variability, causing ice flow to accelerate first and most strongly over the ice shelf. Flow variations then propagate through the grounded ice stream at a mean rate of ∼29 km/d and decay quasi-linearly with distance over ∼85 km upstream of the grounding zone.

Remote Sensing of Ground Deformation for Monitoring Groundwater Management Practices: Application to the Santa Clara Valley During the 2012–2015 California Drought

Groundwater management typically relies on water level data and spatially limited deformation measurements. While InSAR has been used to study hydrological deformation, its limited temporal sampling can lead to biases in rapidly changing systems. Here we use 2011-2017 COSMO-SkyMed data with revisit intervals as short as one day to study the response of the Santa Clara Valley (SCV) aquifer in California to the unprecedented 2012-2015 drought. Cross correlation and independent component analyses of deformation time series enable tracking water though the aquifer system. The aquifer properties are derived prior to and during the drought to assess the success of water-resource management practices. Subsidence due to groundwater withdrawal dominates during 2011-2017, limited to the confined aquifer and west of the Silver Creek Fault, similar to pre-drought summer periods. Minimum water levels and elevations were reached in mid-2014, but thanks to intensive groundwater management efforts the basin started to rebound in summer 2015, during the deepening drought. By 2017, water levels were back to their pre-drought levels, while elevations had not yet fully rebounded due to the delayed poroelastic response of aquitards and their large elastic compressibility. As water levels did not reach a new low stand, the drought led to only elastic and recoverable changes in the SCV. The SCV lost 0.09 km3 during the drought while seasonal variations amount to 0.02 km3. Analysis of surface loads due to water mass changes in the aquifer system suggests that groundwater drawdowns could influence the stress on nearby faults.