Tropospheric Correction for InSAR Using Interpolated ECMWF Data and GPS Zenith Total Delay From the Southern California Integrated GPS Network

A tropospheric correction method for Interferometric Synthetic Aperture Radar (InSAR) was developed using profiles from the European Centre for Medium-Range Weather Forecasts (ECMWF) and Zenith Total Delay (ZTD) from the Global Positioning System (GPS). The ECMWF data were interpolated into a finer grid with the Stretched Boundary Layer Model (SBLM) using a Digital Elevation Model (DEM) with a horizontal resolution of 1 arcsecond. The output were converted into ZTD and combined with the GPS ZTD in order to achieve tropospheric correction maps utilizing both the high spatial resolution of the SBLM and the high accuracy of the GPS. These maps were evaluated for three InSAR images, with short temporal baselines (implying no surface deformation), from Envisat during 2006 on an area stretching northeast from the Los Angeles basin towards Death Valley. The RMS in the InSAR images was greatly reduced, up to 32%, when using the tropospheric corrections. Two of the residuals showed a constant gradient over the area, suggesting a remaining orbit error. This error was reduced by reprocessing the troposphere corrected InSAR images with the result of an overall RMS reduction of 15 − 68%

Damage Proxy Map from Interferometric Synthetic Aperture Radar Coherence

A method, apparatus, and article of manufacture provide the ability to generate a damage proxy map. A master coherence map and a slave coherence map, for an area prior and subsequent to (including) a damage event are obtained. The slave coherence map is registered to the master coherence map. Pixel values of the slave coherence map are modified using histogram matching to provide a first histogram of the master coherence map that exactly matches a second histogram of the slave coherence map. A coherence difference between the slave coherence map and the master coherence map is computed to produce a damage proxy map. The damage proxy map is displayed with the coherence difference displayed in a visually distinguishable manner

Kinematic fault slip evolution source models of the 2008 M7.9 Wenchuan earthquake in China from SAR interferometry, GPS and teleseismic analysis and implications for Longmen Shan tectonics

The M_w 7.9 2008 Wenchuan earthquake ruptured about 280 km of faults in the Longmen Shan of Sichuan province, China, at the eastern edge of the Tibetan Plateau. We use teleseismic waveforms with geodetic data from Global Positioning System, synthetic aperture radar interferometry and image amplitude correlation to produce a source model of this earthquake. The model describes evolution of fault slip during the earthquake. The geodetic data constrains the spatial distribution of fault slip and the seismic waveforms constrain mostly the time evolution of slip. We find that the earthquake started with largely thrust motion on an imbricate system of faults beneath the central Longmen Shan, including the Beichuan Fault and Pengguan Fault, with fault slip at depth extending up to 50 km northwest of the mountain front. The fault ruptures continued northeast along the Beichuan Fault with more oblique slip (right-lateral and thrust) and the proportion of lateral motion increasing in the northern Longmen Shan. The northernmost fault segment has a much steeper dip, consistent with nearly pure strike-slip motion. The kinematic source model shows that the rupture propagated to the northeast at about 2.5–3.0 km s^(−1), producing a cascade of subevents with a total duration of about 110 s. The complex fault ruptures caused shortening and uplift of the extremely steep central Longmen Shan, which supports models where the steep edge of the plateau is formed by thrusting over the strong crust of the Sichuan Basin

Widespread initiation, reactivation, and acceleration of landslides in the northern California Coast Ranges due to extreme rainfall

Episodically to continuously active slow‐moving landslides are driven by precipitation. Climate change, which is altering both the frequency and magnitude of precipitation worldwide, is therefore predicted to have a major impact on landslides. Here we examine the behavior of hundreds of slow‐moving landslides in northern California in response to large changes in annual precipitation that occurred between 2016 and 2018. We quantify the landslide displacement using repeat‐pass radar interferometry and pixel offset tracking techniques on a novel dataset from the airborne NASA/JPL Uninhabited Aerial Vehicle Synthetic Aperture Radar. We found that 312 landslides were moving due to extreme rainfall during 2017, compared to 119 during 2016, which was the final year of a historic multi‐year drought. However, with a return to below‐average rainfall in 2018, only 146 landslides remained in motion. The increased number of landslides during 2017 was primarily accommodated by landslides that were smaller than the landslides that remained active between 2016 and 2018. Furthermore, by examining a subset of 51 landslides, we found that 49 had increased velocities during 2017 when compared to 2016. Our results show that slow‐moving landslides are sensitive to large changes in annual precipitation, particularly the smaller and thinner landslides that likely experience larger basal pore‐water pressure changes. Based on climate model predictions for the next century in California, which include increases in average annual precipitation and increases in the frequency of dry‐to‐wet extremes, we hypothesize that there will be an overall increase in landslide activity

High-Resolution Soil-Moisture Maps Over Landslide Regions in Northern California Grassland Derived From SAR Backscattering Coefficients

Slow-moving landslides are destabilized by accumulated precipitation and consequent soil moisture. Yet, the continuous high-resolution soil-moisture measurements needed to aid the understanding of landslide processes are generally absent in steep terrain. Here, we produce soil-moisture time-series maps for a seasonally active grassland landslide in the northern California coast ranges, USA, using backscattering coefficients from NASA's uninhabited aerial vehicle synthetic aperture radar at 6-m resolution. A physically based radar scattering model is used to retrieve the near-surface (5-cm depth) soil moisture for the landslide. Both forward modeling (backscattering estimation) and the retrieval (soil-moisture validation) show good agreement. The root-mean-square errors (RMSE) for vertical transmit vertical receive (VV) and horizontal transmit horizontal receive (HH) polarizations in forward model comparison are 1.93 dB and 1.88 dB, respectively. The soil-moisture retrieval shows unbiased RMSE of 0.054 m³/m³. Our successful retrieval benefits from the surface and double-bounce scattering, which is common in grasslands. The retrieved maps show saturated wetness conditions within the active landslide boundaries. We also performed sensitivity tests for incidence angle and found that the retrieval is weakly dependent on the angle, especially while using copolarized HH and VV together. Using the two copolarized inputs, the retrieval is also not sensitive to the change of orientation angles of grass cylinders. The physical model inversion presented here can be generally applied for soil-moisture retrieval in areas with the same vegetation cover types in California

Delayed solidification of soft glasses: New experiments, and a theoretical challenge

When subjected to large amplitude oscillatory shear stress, aqueous Laponite suspensions show an abrupt solidification transition after a long delay time tc. We measure the dependence of tc on stress amplitude, frequency, and on the age-dependent initial loss modulus. At first sight our observations appear quantitatively consistent with a simple soft-glassy rheology (SGR)-type model, in which barrier crossings by mesoscopic elements are purely strain-induced. For a given strain amplitude {\gamma}0 each element can be classified as fluid or solid according to whether its local yield strain exceeds {\gamma}0. Each cycle, the barrier heights E of yielded elements are reassigned according to a fixed prior distribution {\rho}(E): this fixes the per-cycle probability R({\gamma}0) of a fluid elements becoming solid. As the fraction of solid elements builds up, {\gamma}0 falls (at constant stress amplitude), so R({\gamma}0) increases. This positive feedback accounts for the sudden solidification after a long delay. The model thus appears to directly link macroscopic rheology with mesoscopic barrier height statistics: within its precepts, our data point towards a power law for {\rho}(E) rather than the exponential form usually assumed in SGR. However, despite this apparent success, closer investigation shows that the assumptions of the model cannot be reconciled with the extremely large strain amplitudes arising in our experiments. The quantitative explanation of delayed solidification in Laponite therefore remains an open theoretical challenge.Comment: 16 pages, 6 figures, to appear in Faraday Discussion

Structural health monitoring of engineered structures using a space-borne synthetic aperture radar multi-temporal approach: from cultural heritage sites to war zones

Structural health monitoring (SHM) of engineered structures consists of an automated or semi-automated survey system that seeks to assess the structural condition of an anthropogenic structure. The aim of an SHM system is to provide insights into possible induced damage or any inherent signals of deformation affecting the structure in terms of detection, localization, assessment, and prediction. During the last decade there has been a growing interest in using several remote sensing techniques, such as synthetic aperture radar (SAR), for SHM. Constellations of SAR satellites with short repeat time acquisitions permit detailed surveys temporal resolution and millimetric sensitivity to deformation that are at the scales relevant to monitoring large structures. The all-weather multi-temporal characteristics of SAR make its products suitable for SHM systems, especially in areas where in situ measurements are not feasible or not cost effective. To illustrate this capability, we present results from COSMO-SkyMed (CSK) and TerraSAR-X SAR observations applied to the remote sensing of engineered structures. We show how by using multiple-geometry SAR-based products which exploit both phase and amplitude of the SAR signal we can address the main objectives of an SHM system including detection and localization. We highlight that, when external data such as rain or temperature records are available or simple elastic models can be assumed, the SAR-based SHM capability can also provide an interpretation in terms of assessment and prediction. We highlight examples of the potential for such imaging capabilities to enable advances in SHM from space, focusing on dams and cultural heritage areas. © (2016) COPYRIGHT Society of Photo-Optical Instrumentation Engineers (SPIE). Downloading of the abstract is permitted for personal use only

Rapid Imaging of Earthquake Ruptures with Combined Geodetic and Seismic Analysis

Rapid determination of the location and extent of earthquake ruptures is helpful for disaster response, as it allows prediction of the likely area of major damage from the earthquake and can help with rescue and recovery planning. With the increasing availability of near real-time data from the Global Positioning System (GPS) and other global navigation satellite system receivers in active tectonic regions, and with the shorter repeat times of many recent and newly launched satellites, geodetic data can be obtained quickly after earthquakes or other disasters. We have been building a data system that can ingest, catalog, and process geodetic data and combine it with seismic analysis to estimate the fault rupture locations and slip distributions for large earthquakes