Tectonics of the central Andes

Acquisition of nearly complete coverage of Thematic Mapper data for the central Andes between about 15 to 34 degrees S has stimulated a comprehensive and unprecedented study of the interaction of tectonics and climate in a young and actively developing major continental mountain belt. The current state of the synoptic mapping of key physiographic, tectonic, and climatic indicators of the dynamics of the mountain/climate system are briefly reviewed

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

Rapid Damage Mapping for the 2015 M_w 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

Evolution of dike opening during the March 2011 Kamoamoa fissure eruption, Kīlauea Volcano, Hawai\u27i

The 5–9 March 2011 Kamoamoa fissure eruption along the east rift zone of Kīlauea Volcano, Hawai`i, followed months of pronounced inflation at Kīlauea summit. We examine dike opening during and after the eruption using a comprehensive interferometric synthetic aperture radar (InSAR) data set in combination with continuous GPS data. We solve for distributed dike displacements using a whole Kīlauea model with dilating rift zones and possibly a deep décollement. Modeled surface dike opening increased from nearly 1.5 m to over 2.8 m from the first day to the end of the eruption, in agreement with field observations of surface fracturing. Surface dike opening ceased following the eruption, but subsurface opening in the dike continued into May 2011. Dike volumes increased from 15, to 16, to 21 million cubic meters (MCM) after the first day, eruption end, and 2 months following, respectively. Dike shape is distinctive, with a main limb plunging from the surface to 2–3 km depth in the up-rift direction toward Kīlauea’s summit, and a lesser projection extending in the down-rift direction toward Pu`u `Ō`ō at 2 km depth. Volume losses beneath Kīlauea summit (1.7 MCM) and Pu`u `Ō`ō (5.6 MCM) crater, relative to dike plus erupted volume (18.3 MCM), yield a dike to source volume ratio of 2.5 that is in the range expected for compressible magma without requiring additional sources. Inflation of Kīlauea’s summit in the months before the March 2011 eruption suggests that the Kamoamoa eruption resulted from overpressure of the volcano’s magmatic system

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