Cécile Lasserre

High-Resolution Satellite Imagery Mapping of the Surface Rupture and Slip Distribution of the Mw 7.8, 14 November 2001 Kokoxili Earthquake, Kunlun Fault, Northern Tibet, China

The Mw 7.8 Kokoxili earthquake of 14 November 2001, which ruptured a 450-km-long stretch of the sinistral Kunlun strike-slip fault, at the northeastern edge of the Tibet plateau, China, ranks as the largest strike-slip event ever recorded instrumentally in Asia. Newly available high-resolution satellite HRS images (pixel size 1 m) acquired in the months following the earthquake proved a powerful tool to complement field investigations and to produce the most accurate map to date of the coseismic displacements along the central Kusai Hu segment of the rupture. The coseismic rupture geometry south and west of Buka Daban Feng, near the earthquake epicenter, was also investigated in detail. Along the Kusai Hu segment, slip parti- tioning is for the first time observed to occur simultaneously during a single event, with two parallel strands, 2 km apart, localizing almost pure strike-slip and normal faulting. In all, 83 new HRS coseismic offset measurements, some of which calibrated by field work, show large, well-constrained variations (100%) of the slip function over distances of only25 km. Tension cracks opening ahead of the shear dislocation and later offset by the upward propagating strike-slip rupture were observed, dem- onstrating that the rupture front propagated faster at depth than near the surface. The triple junction between the central Kusai Hu segment, the Kunlun Pass fault, where the rupture ended, and the Xidatan-Dongdatan segment, which could be the next segment to fail along the main Kunlun fault, acted as a strong barrier, implying that direct triggering of earthquake rupture on the Xidatan-Dongdatan segment by Kokoxili-type earthquakes may not be the rule.

Postglacial left slip rate and past occurrence of M≥8 earthquakes on the Western Haiyuan Fault, Gansu, China

High-resolution (HR) air photographs and a 1-m horizontal and 2-m vertical resolution digital elevation model derived from them by stereophotogrammetry provide new constraints on the behavior of the western stretch of the active Haiyuan fault, in Gansu province, China. The photographs cover three swaths along the fault, each about 2-km-long and at least 500-m-wide, near the village of Songshan, at 103.5oE. This high-resolution data set is used to map and measure cumulative horizontal offsets of alluvial terraces and risers that range between 115 and la 135 m, and 70 and 90 m, at two sites. Dating these terraces with 14C yields minimum and maximum ages of 8400 and 7600, and 14,200 years B.P., respectively. This leads to a postglacial slip rate of 12±4 mm/yr, with a most likely minimum value of 11.6±1.1 mm/yr. The smallest stream offsets observed on the HR photographs range between 8 and 16 m and are interpreted as coseismic displacements of the last few earthquakes with M ≥ 8 that ruptured the 220-km-long Tianzhu gap of the fault, west of the Yellow River. Earthquakes of that size within this gap, which has been quiescent for at least 800 years, would recur at intervals of 1050±450 years.

Long-term slip rates and characteristic slip: keys to active fault behaviour and earthquake hazard

Over periods of thousands of years, active faults tend to slip at constant rates. Pioneer studies of large Asian faults show that cosmogenic radionuclides ( 10 Be, 26 Al) provide an unparalleled tool to date surface features, whose offsets yield the longest records of recent cumulative movement. The technique is thus uniquely suited to determine long-term (10-100 ka) slip rates. Such rates, combined with coseismic slip-amounts, can give access to recurrence times of earthquakes of similar sizes. Landform dating - morphochronology - is therefore essential to understand fault-behaviour, evaluate seismic hazard, and build physical earthquake models. It is irreplaceable because long-term slip- rates on interacting faults need not coincide with GPS-derived, interseismic rates, and can be difficult to obtain from paleo-seismological trenching.  2001 Academie des sciences / Editions scientifiques et medicales Elsevier SAS

Unsupervised Spatiotemporal Mining of Satellite Image Time Series Using Grouped Frequent Sequential Patterns

An important aspect of satellite image time series is the simultaneous access to spatial and temporal information. Various tools allow end users to interpret these data without having to browse the whole data set. In this paper, we intend to extract, in an unsupervised way, temporal evolutions at the pixel level and select those covering at least a minimum surface and having a high connectivity measure. To manage the huge amount of data and the large number of potential temporal evolutions, a new approach based on data-mining techniques is presented. We have developed a frequent sequential pattern extraction method adapted to that spatiotemporal context. A successful application to crop monitoring involving optical data is described. Another application to crustal deformation monitoring using synthetic aperture radar images gives an indication about the generic nature of the proposed approach.

Millennial recurrence of large earthquakes on the Haiyuan fault near Songshan, Gansu Province, China

The Haiyuan fault is a major active left-lateral fault along the northeast edge of the Tibet-Qinghai Plateau. Studying this fault is important in understanding current deformation of the plateau and the mechanics of continental deformation in general. Previous studies have mostly focused on the slip rate of the fault. Paleo- seismic investigations on the fault are sparse, and have been targeted mostly at the stretch of the fault that ruptured in the 1920 M 8.6 earthquake in Ningxia Province. To investigate the millennial seismic history of the western Haiyuan fault, we opened two trenches in a small pull-apart basin near Songshan, in Gansu Province. The excavation exposes sedimentary layers of alternating colors: dark brown silty to clayey deposit and light yellowish brown layers of coarser-grained sandy deposit. The main fault zone is readily recognizable by the disruption and tilting of the layers. Six paleoseismic events are identified and named SS1 through SS6, from youngest to oldest. Charcoal is abundant, yet generally tiny in the shallowest parts of the trench exposures. Thirteen samples were dated to constrain the ages of paleoseismic events. All six events have occurred during the past 3500-3900 years. The horizontal offsets associated with these events are poorly known. However, events SS3 to SS6 appear to be large ones, judging from comparison of vertical separations and widths of fault zones. The youngest event SS1 instead seems to be a minor one, probably the 1990 Mw 5.8 earthquake. Thus, four large events in 3500-3900 years would imply a re- currence interval of about 1000 years. Three events SS2 to SS4 prior to 1990 occurred sometime during 1440-1640 A.D., shortly after 890-1000 A.D. and 0-410 A.D., re- spectively. We tentatively associate them with the 1514 A.D., 1092 A.D., and 143 or 374 A.D. historical earthquakes. Taking 10 2 m of slip for large events (SS3 and SS4), comparable to the 1920 M 8 Haiyuan earthquake, their occurrence times would be consistent with the long-term 12 4 mm/yr estimate of Lasserre et al. (1999). However, a more realistic evaluation of slip rate and its possible change with time requires a more rigorous determination of coseismic slip amounts of past earth- quakes.

New Radar Interferometric Time Series Analysis Toolbox Released

Interferometric synthetic aperture radar (InSAR) has become an important geodetic tool for measuring deformation of Earth’s surface due to various geophysical phenomena, including slip on earthquake faults, subsurface migration of magma, slow‐moving landslides, movement of shallow crustal fluids (e.g., water and oil), and glacier flow. Airborne and spaceborne synthetic aperture radar (SAR) instruments transmit microwaves toward Earth’s surface and detect the returning reflected waves. The phase of the returned wave depends on the distance between the satellite and the surface, but it is also altered by atmospheric and other effects. InSAR provides measurements of surface deformation by combining amplitude and phase information from two SAR images of the same location taken at different times to create an interferogram. Several existing open‐source analysis tools [Rosen et al., 2004; Rosen et al., 2011; Kampes et al., 2003 ; Sandwell et al., 2011] enable scientists to exploit observations from radar satellites acquired at two different epochs to produce a surface displacement map.

An aseismic slip transient on the North Anatolian Fault

Constellations of Synthetic Aperture Radar (SAR) satellites with short repeat time acquisitions allow exploration of active faults behavior with unprecedented temporal resolution. Along the North Anatolian Fault (NAF) in Turkey, an 80 km long section has been creeping at least since the 1944, M_w 7.3 earthquake near Ismetpasa, with a current Interferometric Synthetic Aperture Radar (InSAR)-derived average creep rate of 8 ± 3 mm/yr (i.e., a third of the NAF long-term slip rate). We use a dense set of SAR images acquired by the COSMO-SkyMed constellation to quantify the spatial distribution and temporal evolution of creep over 1 year. We identify a major burst of aseismic slip spanning 31 days with a maximum slip of 2 cm, between the surface and 4 km depth. This result shows that fault creep along this section of the NAF does not occur at a steady rate as previously thought, highlighting a need to revise our understanding of the underlying fault mechanics.

Slow Slip History for the MEXICO Subduction Zone: 2005 Through 2011

AbstractTo further our understanding of the seismically hazardous Mexico subduction zone, we estimate the first time-dependent slip distributions and Coulomb failure stress changes for the six major slow slip events (SSEs) that occurred below Mexico between late 2005 and mid-2011. Slip dist ributions are the first to be estimated from all continuous GPS data in central and southern Mexico, which better resolves slow slip in space and time than was previously possible in this region. Below Oaxaca, slip during previously un-modeled SSEs in 2008/9 and 2010/11 extended farther to the west than previous SSEs. This constitutes the first evidence that slow slip accounts for deep slip within a previously noted gap between the Oaxaca and Guerrero SSE source regions. The slip that we estimate for the two SSEs that originated below Guerrero between 2005 and 2011 agrees with slip estimated in previous, mostly static-offset SSE modeling studies; however, we show that both SSEs migrated eastward toward the Oaxaca SSE source region. In accord with previous work, we find that slow slip below Guerrero intrudes up-dip into the potentially seismogenic region, presumably accounting for some of the missing slip within the well-described Guerrero seismic gap. In contrast, slow slip below Oaxaca between 2005 and 2011 occurred mostly down-dip from the seismogenic regions defined by the rupture zones of large thrust earthquakes in 1968 and 1978 and released all of the slip deficit that accumulated in the down-dip region during this period.

DEM Corrections Before Unwrapping in a Small Baseline Strategy for InSAR Time Series Analysis

Synthetic aperture radar interferometry (InSAR) is limited by temporal decorrelation and topographic errors, which can result in unwrapping errors in partially incoherent and mountainous areas. In this paper, we present an algorithm to estimate and remove local digital elevation model (DEM) errors from a series of wrapped interferograms. The method is designed to be included in a small baseline subset (SBAS) approach for InSAR time series analysis of ground deformation in natural environment. It is easy to implement and can be applied to all pixels of a radar scene. The algorithm is applied to a series of wrapped interferograms computed from ENVISAT radar images acquired across the Himalayan mountain range. The DEM error correction performance is quantified by the reduction of the local phase dispersion and of the number of residues computed during the unwrapping procedure. It thus improves the automation of the spatial unwrapping step.

Lateral Variations of Interplate Coupling along the Mexican Subduction Interface: Relationships with Long-Term Morphology and Fault Zone Mechanical Properties

Although patterns of interseismic strain accumulation above subduction zones are now routinely characterised using geodetic measurements, their physical origin, persistency through time, and relationships to seismic hazard and long-term deformation are still debated. Here, we use GPS and morphological observations from southern Mexico to explore potential mechanical links between variations in inter-SSE (in between slow slip events) coupling along the Mexico subduction zone and the long-term topography of the coastal regions from Guerrero to Oaxaca. Inter-SSE coupling solutions for two different geometries of the subduction interface are derived from an inversion of continuous GPS time series corrected from slow slip events. They reveal strong along-strike variations in the shallow coupling (i.e. at depths down to 25 km), with high-coupling zones (coupling >0.7) alternating with low-coupling zones (coupling <0.3). Coupling below the continent is typically strong (>0.7) and transitions to uncoupled, steady slip at a relatively uniform