Optical Offsets 7To determine the co-seismic horizontal displacement field due to the Gorkha earthquake, we use optical 8 image correlation to measure the displacement of pixels between pre-and post-earthquake satellite im-9 ages. We are able to resolve sub-pixel displacements of less than 1/15th of the Landsat8 pixel resolu-10 tion (i.e. < 1 m) using the COSI-Corr software package images, which helps to increase the signal-to-noise ratio, (4) the deformation field is resolved perpendicular to 16 the look angle (i.e. the horizontal plane for nadir images), thereby providing measurements complementary 17to InSAR (which is sensitive to vertical displacements), (5) the nadir look angle is insensitive to topographic 18 residuals produced during orthorectification of the satellite images (such residuals are produced when a lower 19 resolution digital elevation model, DEM, is used during the orthorectification process), and (5) Landsat8 20 images are freely available from the USGS as an orthorectified product -see 5 for additional details. 21Landsat8 images are typically acquired at 10am each morning. Consequently, the illumination charac-22 teristics (i.e. shadows) vary in every image acquired throughout the year according to the position of the 23 sun. Because shadows produce sharp edges in satellite images, they strongly influence the correlation. There-24 fore, to reduce the effect of differing shadows biasing the displacement field, we correlate Landsat8 images 25 acquired at a similar time of year, thereby yielding images with similar illumination characteristics (i.e. sun 26 azimuth and elevation). In addition to having similar illumination characteristics, we also require images with 27 minimal cloud cover. From the Landsat8 archive, we found two suitable images from the (pre-earthquake) 2813th May 2014 (sun azimuth: 109• , sun elevation: 68

6p.International audienceDespite 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

International audienceWe analyse the Mw 7.7 Balochistan earthquake of 09/24/2013 based on ground surface deformation measured from sub-pixel correlation of Landsat-8 images, combined with back-projection and finite source modeling of teleseismic waveforms. The earthquake nucleated south of the Chaman strike-slip fault and propagated southwestward along the Hoshab fault at the front of the Kech Band. The rupture was mostly unilateral, propagated at 3 km/s on average and produced a 200 km surface fault trace with purely strike-slip displacement peaking to 10 m and averaging around 6 m. The finite source model shows that slip was maximum near the surface. Although the Hoshab fault is dipping by 45° to the North, in accordance with its origin as a thrust fault within the Makran accretionary prism, slip was nearly purely strike-slip during that earthquake. Large seismic slip on such a non-optimally oriented fault was enhanced possibly due to the influence of the free surface on dynamic stresses or to particular properties of the fault zone allowing for strong dynamic weakening. Strike-slip faulting on thrust fault within the eastern Makran is interpreted as due to eastward extrusion of the accretionary prism as it bulges out over the Indian plate. Portions of the Makran megathrust, some thrust faults in the Kirthar range and strike-slip faults within the Chaman fault system have been brought closer to failure by this earthquake. Aftershocks cluster within the Chaman fault system north of the epicenter, opposite to the direction of rupture propagation. By contrast, few aftershocks were detected in the area of maximum moment release. In this example, aftershocks cannot be used to infer earthquake characteristics

[1] Interferometric synthetic aperture radar data are used to map the interseismic velocity field along the Haiyuan fault system (HFS), at the north-eastern boundary of the Tibetan plateau. Two M $ 8 earthquakes ruptured the HFS in 1920 and 1927, but its 260 km-long central section, known as the Tianzhu seismic gap, remains unbroken since $1000 years. The Envisat SAR data, spanning the 2003-2009 period, cover about 200 Â 300 km 2 along three descending and two ascending tracks. Interferograms are processed using an adapted version of ROI_PAC. The signal due to stratified atmospheric phase delay is empirically corrected together with orbital residuals. Mean line-of-sight velocity maps are computed using a constrained time series analysis after selection of interferograms with low atmospheric noise. These maps show a dominant left-lateral motion across the HFS, and reveal a narrow, 35 km-long zone of high velocity gradient across the fault in between the Tianzhu gap and the 1920 rupture. We model the observed velocity field using a discretized fault creeping at shallow depth and a least squares inversion. The inferred shallow slip rate distribution reveals aseismic slip in between two fully locked segments. The average creep rate is $5 mm yr À1 , comparable in magnitude with the estimated loading rate at depth, suggesting no strain accumulation on this segment. The modeled creep rate locally exceeds the long term rate, reaching 8 mm yr À1, suggesting transient creep episodes. The present study emphasizes the need for continuous monitoring of the surface velocity in the vicinity of major seismic gaps in terms of seismic hazard assessment.

Spatial and temporal variations of pressure, temperature, and water vapor content in the atmosphere introduce significant confounding delays in interferometric synthetic aperture radar (InSAR) observations of ground deformation and bias estimates of regional strain rates. Producing robust estimates of tropospheric delays remains one of the key challenges in increasing the accuracy of ground deformation measurements using InSAR. Recent studies revealed the efficiency of global atmospheric reanalysis to mitigate the impact of tropospheric delays, motivating further exploration of their potential. Here we explore the effectiveness of these models in several geographic and tectonic settings on both single interferograms and time series analysis products. Both hydrostatic and wet contributions to the phase delay are important to account for. We validate these path delay corrections by comparing with estimates of vertically integrated atmospheric water vapor content derived from the passive multispectral imager Medium-Resolution Imaging Spectrometer, onboard the Envisat satellite. Generally, the performance of the prediction depends on the vigor of atmospheric turbulence. We discuss (1) how separating atmospheric and orbital contributions allows one to better measure long-wavelength deformation and (2) how atmospheric delays affect measurements of surface deformation following earthquakes, and (3) how such a method allows us to reduce biases in multiyear strain rate estimates by reducing the influence of unevenly sampled seasonal oscillations of the tropospheric delay. IntroductionSynthetic aperture radar interferometry (InSAR) has been successfully used to measure ground deformations related to hydrologic, volcanic, and tectonic processes [e.g., Bawden et al., 2001;Beauducel et al., 2000;Massonnet et al., 1992]. Rapid, large-amplitude deformation signals such as coseismic displacement fields [e.g., Simons et al., 2002;Lasserre et al., 2005] or volcano-tectonic episodes [e.g., Pritchard and Simons, 2002;Wright et al., 2004;Doubre and Peltzer, 2007;Grandin et al., 2010] are now routinely measured by InSAR. Still, the detection of low-amplitude, long-wavelength deformation fields such as those due to interseismic strain accumulation or postseismic motion remains challenging because of interferometric decorrelation, inaccurate orbits, and atmospheric propagation delays [e.g., Peltzer et al., 2001;Wright et al., 2001;Ryder et al., 2007;Wen et al., 2012;Jolivet et al., 2012;Grandin et al., 2012;Béjar-Pizarro et al., 2013]. Here we focus on a specific method to mitigate the impact of atmospheric artifacts.Spatiotemporal variations of the refractivity of air can introduce a change in the measured interferometric phase, hereafter called the atmospheric phase screen (APS). This phase change, or phase delay, can be on the order of several centimeters and often overwhelms the deformation signal of interest [Hanssen, 2001]. These phase delays result from the combined effects of turbulent mixing in the atmosphere (hereaft...

International audienceWe use high-resolution Synthetic Aperture Radar- and GPS-derived observations of surfacedisplacements to derive the first probabilistic estimates of fault coupling along the creeping section of theSan Andreas Fault, in between the terminations of the 1857 and 1906 magnitude 7.9 earthquakes. Usinga fully Bayesian approach enables unequaled resolution and allows us to infer a high probability ofsignificant fault locking along the creeping section. The inferred discreet locked asperities are consistentwith evidence for magnitude 6+ earthquakes over the past century in this area and may be associated withthe initiation phase of the 1857 earthquake. As creeping segments may be related to the initiation andtermination of seismic ruptures, such distribution of locked and creeping asperities highlights the centralrole of the creeping section on the occurrence of major earthquakes along the San Andreas Fault

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.

No abstract