Constraints on the Geometry and Frictional Properties of the Main Himalayan Thrust Using Coseismic, Postseismic, and Interseismic Deformation in Nepal

The geometry and frictional properties of a fault system are key parameters required to understand its seismic behavior. The Main Himalayan Thrust in Nepal is the type example of a continental megathrust and forms part of a fault system which accommodates a significant fraction of India‐Eurasia convergence. Despite extensive study of this zone of shortening, the geometry of the fault system remains controversial. Here, we use interseismic, coseismic, and postseismic geodetic data in Nepal to investigate the proposed downdip geometries. We use interseismic and coseismic data from previous studies, acquired before and during the 2015 urn:x-wiley:jgrb:media:jgrb53987:jgrb53987-math-0001 7.8 Gorkha earthquake. We then supplement these by processing our own postseismic deformation data, acquired following the Gorkha earthquake. We find that kinematic modeling of geodetic data alone cannot easily distinguish between the previously proposed geometries. We therefore develop a mechanical joint coseismic‐postseismic slip inversion which simultaneously solves for the distribution of coseismic slip and rate‐strengthening friction parameters. We run this inversion using the proposed geometries and find that they are all capable of explaining the majority of geodetic data. We find values for the rate parameter, urn:x-wiley:jgrb:media:jgrb53987:jgrb53987-math-0002, from the rate‐and‐state friction law that are between 0.8 and urn:x-wiley:jgrb:media:jgrb53987:jgrb53987-math-0003, depending on the geometry used. These values are in agreement with results from laboratory studies and those inferred from other earthquakes. We suggest that the limitations of earthquake cycle geodesy partly explain the continued controversy over the geometry and role of various faults in the Nepal Himalaya

A Network Inversion Filter combining GNSS and InSAR for tectonic slip modeling

Studies of the earthquake cycle benefit from long-term time-dependent slip modeling, as it can be a powerful means to improve our understanding on the interaction of earthquake cycle processes such as interseismic, coseismic, postseismic, and aseismic slip. Observations from Interferometric Synthetic Aperture Radar (InSAR) allow us to model slip at depth with a higher spatial resolution than when using GNSS alone. While the temporal resolution of InSAR has typically been limited, the recent fleet of SAR satellites including Sentinel-1, COSMO-SkyMED, and RADARSAT-2 permits the use of InSAR for time-dependent slip modeling, at intervals of a few days when combined. With the vast amount of SAR data available, simultaneous data inversion of all epochs becomes challenging. Here, we expanded the original Network Inversion Filter to include InSAR observations of surface displacements in addition to GNSS. In the NIF framework, geodetic observations are limited to those of a given epoch, with a stochastic model describing slip evolution over time. The combination of the Kalman forward filtering and backward smoothing allows all geodetic observations to constrain the complete observation period. Combining GNSS and InSAR allows modeling of time-dependent slip at unprecedented spatial resolution. We validate the approach with a simulation of the 2006 Guerrero slow slip event. We highlight the importance of including InSAR covariance information, and demonstrate that InSAR provides an additional constraint on the spatial extent of the slow slip

Statistical comparison of InSAR tropospheric correction techniques

Correcting for tropospheric delays is one of the largest challenges facing the interferometric synthetic aperture radar (InSAR) community. Spatial and temporal variations in temperature, pressure, and relative humidity create tropospheric signals in InSAR data, masking smaller surface displacements due to tectonic or volcanic deformation. Correction methods using weather model data, GNSS and/or spectrometer data have been applied in the past, but are often limited by the spatial and temporal resolution of the auxiliary data. Alternatively a correction can be estimated from the interferometric phase by assuming a linear or a power-law relationship between the phase and topography. Typically the challenge lies in separating deformation from tropospheric phase signals. In this study we performed a statistical comparison of the state-of-the-art tropospheric corrections estimated from the MERIS and MODIS spectrometers, a low and high spatial-resolution weather model (ERA-I and WRF), and both the conventional linear and new power-law empirical methods. Our test-regions include Southern Mexico, Italy, and El Hierro. We find spectrometers give the largest reduction in tropospheric signal, but are limited to cloud-free and daylight acquisitions. We find a ~ 10–20% RMSE increase with increasing cloud cover consistent across methods. None of the other tropospheric correction methods consistently reduced tropospheric signals over different regions and times. We have released a new software package called TRAIN (Toolbox for Reducing Atmospheric InSAR Noise), which includes all these state-of-the-art correction methods. We recommend future developments should aim towards combining the different correction methods in an optimal manner

The Sentinel-1 constellation for InSAR applications: Experiences from the InSARAP project

The two-satellite Copernicus Sentinel-1 (S1) constellation became operational in Sep 2016, with the successful in-orbit commissioning of the S1B unit. During, the commissioning phase and early operational phase it has been confirmed that the interferometric performance of the constellation is excellent, with no observed phase anomalies. In this work, we show an analysis of selected performance parameters for the S1 constellation, as well as initial results based on the available data from the first months of operations

LiCSBAS: An Open-Source InSAR Time Series Analysis Package Integrated with the LiCSAR Automated Sentinel-1 InSAR Processor

For the past five years, the 2-satellite Sentinel-1 constellation has provided abundant and useful Synthetic Aperture Radar (SAR) data, which have the potential to reveal global ground surface deformation at high spatial and temporal resolutions. However, for most users, fully exploiting the large amount of associated data is challenging, especially over wide areas. To help address this challenge, we have developed LiCSBAS, an open-source SAR interferometry (InSAR) time series analysis package that integrates with the automated Sentinel-1 InSAR processor (LiCSAR). LiCSBAS utilizes freely available LiCSAR products, and users can save processing time and disk space while obtaining the results of InSAR time series analysis. In the LiCSBAS processing scheme, interferograms with many unwrapping errors are automatically identified by loop closure and removed. Reliable time series and velocities are derived with the aid of masking using several noise indices. The easy implementation of atmospheric corrections to reduce noise is achieved with the Generic Atmospheric Correction Online Service for InSAR (GACOS). Using case studies in southern Tohoku and the Echigo Plain, Japan, we demonstrate that LiCSBAS applied to LiCSAR products can detect both large-scale (>100 km) and localized (~km) relative displacements with an accuracy of <1 cm/epoch and ~2 mm/yr. We detect displacements with different temporal characteristics, including linear, periodic, and episodic, in Niigata, Ojiya, and Sanjo City, respectively. LiCSBAS and LiCSAR products facilitate greater exploitation of globally available and abundant SAR datasets and enhance their applications for scientific research and societal benefit

Rift Focusing and Magmatism During Late-Stage Rifting in Afar

Processes that facilitate the transition between continental rifting and sea-floor spreading remain unclear. Variations in the spatial distribution of extension and magmatism through Afar and into the Red Sea are indicative of temporal evolution of the rifting process. We develop a time series of Sentinel-1 interferometric synthetic aperture radar (InSAR) observations of ground deformation covering the whole Afar rift zone from 2014 to 2019, to study the distribution of extension. By incorporating Global Navigation Satellite System observations, we resolve 3D average velocities in the vertical, rift-perpendicular, and rift-parallel directions. Results show the spatial distribution of long-wavelength deformation over the rift zone, as well as deformation at individual volcanic centers, including Dallol, Nabro, and Erta ’Ale. We find that in northern and central Afar, the majority of extension is accommodated close to the rift axis (±15–30 km). In southern Afar, near the Nubia-Arabia-Somalia triple junction, amagmatic extension is distributed over 80–160 km, which may indicate an increase in rift focusing with rift maturity. We also observe rapid surface uplift and rift-perpendicular extension at the Dabbahu-Manda-Hararo rift segment with velocities of 33 ± 4 mm/yr and 37 ± 4 mm/yr respectively. These are higher than the background extension rate of 18–20 mm/yr, but have decreased by 55%–70% since 2006–2010. The data suggests that this is due to an ongoing long-lived response to the 2005–2010 rifting episode, with potential continued processes below the rift segment including a lower-crustal viscous response and magma movement. Continued observations of surface deformation provide key constraints on tectono-magmatic processes involved in rift development

Improving the Resolving Power of InSAR for Earthquakes Using Time Series: A Case Study in Iran

Interferometric Synthetic Aperture Radar (InSAR) is an established method to measure earthquake surface displacements. However, due to decorrelation and atmospheric noise, only a certain fraction of earthquakes is readily observable with single interferograms. To enhance the potential of retrieving InSAR earthquake observations, we apply InSAR time series analysis and use several recent earthquakes (Mw 5.6–6.3, 2018–2019) in Iran as case studies. We find that the coseismic displacement signals of these earthquakes, which might not be discernible within single interferograms, are better resolved using our approach. We reconstruct the coseismic deformation fields by fitting surface displacements using a time series approach. We find that the reconstructed coseismic deformation fields yield more robust and seismologically consistent earthquake modeling results when compared to single coseismic interferograms. Our work suggests that a time series approach is an effective way to improve the resolving power of InSAR for earthquake studies

Evaluation of the Multilook Size in Polarimetric Optimization of Differential SAR Interferograms

The interferometric coherence is a measure of the correlation between two SAR images and constitutes a commonly used estimator of the phase quality. Its estimation requires a spatial average within a 2-D window, usually named as multilook. The multilook processing allows reducing noise at the expenses of a resolution loss. In this letter, we analyze the influence of the multilook size while applying a polarimetric optimization of the coherence. The same optimization algorithm has been carried out with different multilook sizes and also with the nonlocal SAR filter filter, which has the advantage of preserving the original resolution of the interferogram. Our experiments have been carried out with a single pair of quad-polarimetric RADARSAT-2 images mapping the Mount Etna's volcanic eruption of May 2008. Results obtained with this particular data set show that the coherence is increased notably with respect to conventional channels when small multilook sizes are employed, especially over low-vegetated areas. Conversely, very decorrelated areas benefit from larger multilook sizes but do not exhibit an additional improvement with the polarimetric optimization

Topological Modes in One Dimensional Solids and Photonic Crystals

This is the final version of the article. Available from American Physical Society via the DOI in this record.It is shown theoretically that a one-dimensional crystal with time reversal symmetry is characterized by a Z_{2} topological invariant that predicts the existence or otherwise of edge states. This is confirmed experimentally through the construction and simulation of a photonic crystal analogue in the microwave regime.We thank an anonymous referee for helpful comments. The authors acknowledge financial support from the EPSRC through the QUEST program grant (Grant No. EP/I034548/1), C.A.M.B. was supported by QinetiQ and the EPSRC through the Industrial CASE scheme (Grant No. 08000346). H.M. is supported by a grant from the U. S. Department of Energy to the Particle Astrophysics Theory group at CWRU. T.J.A. is funded by the Research Corporation for Science Advancement through a Cottrell Award

Novel corner-reflector array application in essential infrastructure monitoring

High-precision monitoring of infrastructure using artificial reflectors is possible with freely available Sentinel-1 data, but large reflectors are needed. We find that a triangular trihedral corner reflector should typically have at least 1-m inner leg length. As such large reflectors are often not feasible for use in urban areas for essential infrastructure monitoring, we designed a multiple corner-reflector array to replace a single corner reflector with an inner leg length of 1 m. In this case, we use four reflectors where each of them is a truncated triangular trihedral with an inner leg length of 0.33 m. We measured interferometric synthetic aperture radar (InSAR) amplitude, phase, and coherence of this reflector array with various configurations of alignments of the array. We find that as long as great care is taken in the relative positioning of the four corner reflectors, so that they constructively interfere, each horizontal or vertical configuration provides the expected amplitude, coherence, and phase stability. Applications of multiple small corner reflectors in urban areas range from essential infrastructure monitoring (e.g., bridges, overpasses, and tunnel constructions), through assessment of structural health of buildings, to monitoring highway and railway embankments. We show that the multiple corner array works when placed in a single InSAR resolution cell, but depending on the application, the number and projection of corner reflectors can be varied, as long as sufficient signal-to-clutter ratio is achieved in the area of interest