ALBERT

Description: Recent observations suggesting the influence of creep on earthquakes nucleation and arrest are strong incentives to investigate the physical mechanisms controlling how active faults slip. We focus here on deriving generic characteristics of shallow creep along the Haiyuan fault, a major strike-slip fault in China, by investigating the relationship between fault slip and geometry. We use optical images and time series of Synthetic Aperture Radar data to map the surface fault trace and the spatiotemporal distribution of surface slip along the creeping section of the Haiyuan fault. The fault trace roughness shows a power-law behavior similar to that of the aseismic slip distribution, with a 0.8 roughness exponent, typical of a self-affine regime. One possible interpretation is that fault geometry controls to some extent the distribution of aseismic slip, as it has been shown previously for coseismic slip along active faults. Creep is characterized by local fluctuations in rates that we define as creep bursts. The potency of creep bursts follows a power-law behavior similar to the Gutenberg–Richter earthquake distribution, whereas the distribution of bursts velocity is non-Gaussian, suggesting an avalanche-like behavior of these slip events. Such similarities with earthquakes and lab experiments lead us to interpret the rich dynamics of creep bursts observed along the Haiyuan fault as resulting from long-range elastic interactions within the heterogeneous Earth’s crust.

Description: Oblique convergence across Tibet leads to slip partitioning with the coexistence of strike-slip, normal and thrust motion on major fault systems. A key point is to understand and model how faults interact and accumulate strain at depth. Here, we extract ground deformation across the Haiyuan Fault restraining bend, at the northeastern boundary of the Tibetan plateau, from Envisat radar data spanning the 2001–2011 period. We show that the complexity of the surface displacement field can be explained by the partitioning of a uniform deep-seated convergence. Mountains and sand dunes in the study area make the radar data processing challenging and require the latest developments in processing procedures for Synthetic Aperture Radar interferometry. The processing strategy is based on a small baseline approach. Before unwrapping, we correct for atmospheric phase delays from global atmospheric models and digital elevation model errors. A series of filtering steps is applied to improve the signal-to-noise ratio across high ranges of the Tibetan plateau and the phase unwrapping capability across the fault, required for reliable estimate of fault movement. We then jointly invert our InSAR time-series together with published GPS displacements to test a proposed long-term slip-partitioning model between the Haiyuan and Gulang left-lateral Faults and the Qilian Shan thrusts. We explore the geometry of the fault system at depth and associated slip rates using a Bayesian approach and test the consistency of present-day geodetic surface displacements with a long-term tectonic model. We determine a uniform convergence rate of 10 [8.6–11.5] mm yr –1 with an N89 [81–97]°E across the whole fault system, with a variable partitioning west and east of a major extensional fault-jog (the Tianzhu pull-apart basin). Our 2-D model of two profiles perpendicular to the fault system gives a quantitative understanding of how crustal deformation is accommodated by the various branches of this thrust/strike-slip fault system and demonstrates how the geometry of the Haiyuan fault system controls the partitioning of the deep secular motion.

Description: Great earthquakes rarely occur within active accretionary prisms, despite the intense long-term deformation associated with the formation of these geologic structures. This paucity of earthquakes is often attributed to partitioning of deformation across multiple structures as well as aseismic deformation within and at the base of the prism ( Davis et al. , 1983 ). We use teleseismic data and satellite optical and radar imaging of the 2013 M w  7.7 earthquake that occurred on the southeastern edge of the Makran plate boundary zone to study this unexpected earthquake. We first compute a multiple point-source solution from W -phase waveforms to estimate fault geometry and rupture duration and timing. We then derive the distribution of subsurface fault slip from geodetic coseismic offsets. We sample for the slip posterior probability density function using a Bayesian approach, including a full description of the data covariance and accounting for errors in the elastic structure of the crust. The rupture nucleated on a subvertical segment, branching out of the Chaman fault system, and grew into a major earthquake along a 50° north-dipping thrust fault with significant along-strike curvature. Fault slip propagated at an average speed of 3.0 km/s for about 180 km and is concentrated in the top 10 km with no displacement on the underlying décollement. This earthquake does not exhibit significant slip deficit near the surface, nor is there significant segmentation of the rupture. We propose that complex interaction between the subduction accommodating the Arabia–Eurasia convergence to the south and the Ornach Nal fault plate boundary between India and Eurasia resulted in the significant strain gradient observed prior to this earthquake. Convergence in this region is accommodated both along the subduction megathrust and as internal deformation of the accretionary wedge. Online Material: Figures showing waveform fits, focal mechanism, root mean square misfit, variation of the shear modulus, covariance functions, and model predictions.

Description: Despite remarkable successes achieved by Differential InSAR, estimations of low tectonic strain rates remain challenging in areas where deformation and topography are correlated, mainly because of the topography-related atmospheric phase screen (APS). In areas of high relief, empirical removal of the stratified component of the APS may lead to biased estimations of tectonic deformation rates. Here we describe a method to correct interferograms from the effects of the spatial and temporal variations in tropospheric stratification by computing tropospheric delay maps coincident with SAR acquisitions using the ERA-Interim global meteorological model. The modeled phase delay is integrated along vertical profiles at the ERA-I grid nodes and interpolated at the spatial sampling of the interferograms above the elevation of each image pixel. This approach is validated on unwrapped interferograms. We show that the removal of the atmospheric signal before phase unwrapping reduces the risk of unwrapping errors in areas of rough topography.

Description: The rise and support of the ~5000 m topographic scarp at the front of Indian-Eurasian collision in the Himalaya involves long-term uplift above a mid-crustal ramp within the Main Himalayan Thrust (MHT) system. Locking of the shallower portion of the flat-ramp-flat during the interseismic period also produces transient uplift above the transition zone. However, spatial and temporal relationships between permanent and transient vertical deformation in the Himalaya are poorly constrained, leading to an unresolved causal relationship between the two. Here, we use interferometric synthetic aperture radar (InSAR) to measure interseismic uplift on a transect crossing the whole Himalaya in central Nepal. The uplift velocity of 7 mm/yr at the front of the Annapurna mountain range is explained by an 18–21 mm/yr slip rate on the deep shallow-dipping portion of the MHT, with full locking of the mid-crustal ramp underlying the High Himalaya. The transient uplift peak observed by InSAR matches spatially with the long-term uplift peak deduced from the study of trans-Himalayan river incision, although models of the seismic cycle involving thrusting over a ramp of fixed geometry predict an ~20 km separation between the two peaks. We argue that this coincidence indicates that today’s mid-crustal ramp in central Nepal is located southward with respect to its average long-term location, suggesting that mountain growth proceeds by frontward migration of the ramp driven by underplating of material from the Indian plate under the Himalaya.

Description: Interferometric synthetic aperture radar time series methods aim to reconstruct time-dependent ground displacements over large areas from sets of interferograms in order to detect transient, periodic, or small-amplitude deformation. Because of computational limitations, most existing methods consider each pixel independently, ignoring important spatial covariances between observations. We describe a framework to reconstruct time series of ground deformation while considering all pixels simultaneously, allowing us to account for spatial covariances, imprecise orbits, and residual atmospheric perturbations. We describe spatial covariances by an exponential decay function dependent of pixel-to-pixel distance. We approximate the impact of imprecise orbit information and residual long-wavelength atmosphere as a low-order polynomial function. Tests on synthetic data illustrate the importance of incorporating full covariances between pixels in order to avoid biased parameter reconstruction. An example of application to the northern Chilean subduction zone highlights the potential of this method. ©2018. American Geophysical Union. All Rights Reserved.

Description: Large earthquakes in stable continental regions remain puzzling as; unlike at plate boundaries, they do not result from the local buildup of strain driven by plate tectonics. The 2017 Mw6.5, Bostwana normal faulting earthquake occurred in a region devoid from recent tectonic activity and where present-day deformation is negligible. The depth of the event (29 ± 4 km), in a felsic lower crust where ductile deformation is expected, likely requires a transient pulse of fluids from a deep source to activate brittle faulting. The mainshock was preceded by two foreshock swarm-like sequences that may be further evidence for fluid movement in a critically loaded fault network. Contrary to plate boundary events, the Mw6.5 Botswana earthquake did not require prior localized stress or strain accumulation. We propose that the crust in stable continental regions, even long after the last tectonic episode, constitutes a reservoir of elastic stress that can be released episodically, for instance, as a result of deep fluid migration. ©2018. American Geophysical Union. All Rights Reserved.

Description: On 12 November 2017, a MW=7.3 earthquake struck near the Iranian town of Ezgeleh, at the Iran-Iraq border. This event was located within the Zagros fold and thrust belt which delimits the continental collision between the Arabian and Eurasian Plates. Despite a high seismic risk, the seismogenic behavior of the complex network of active faults is not well documented in this area due to the long recurrence interval of large earthquakes. In this study, we jointly invert interferometric synthetic aperture radar and near-field strong motions to infer a kinematic slip model of the rupture. The incorporation of these near-field observations enables a fine resolution of the kinematic rupture process. It reveals an impulsive seismic source with a strong southward rupture directivity, consistent with significant damage south of the epicenter. We also show that the slip direction does not match plate convergence, implying that some of the accumulated strain must be partitioned onto other faults. ©2019. The Authors.