Sébastien Leprince

Earth-like sand fluxes on Mars

Strong and sustained winds on Mars have been considered rare, on the basis of surface meteorology measurements and global circulation models, raising the question of whether the abundant dunes and evidence for wind erosion seen on the planet are a current process. Recent studies showed sand activity, but could not determine whether entire dunes were moving—implying large sand fluxes—or whether more localized and surficial changes had occurred. Here we present measurements of the migration rate of sand ripples and dune lee fronts at the Nili Patera dune field. We show that the dunes are near steady state, with their entire volumes composed of mobile sand. The dunes have unexpectedly high sand fluxes, similar, for example, to those in Victoria Valley, Antarctica, implying that rates of landscape modification on Mars and Earth are similar.

Monitoring Earth Surface Dynamics With Optical Imagery

Despite the increasing availability of high-quality optical satellite images, continuous monitoring of Earths surface changes is still of limited use due to technical limitations. To overcome these limitations, this thesis presents a processing chain to accurately orthorectify and co-register sets of satellite and aerial images, which, associated with a precise correlation technique, allow for the measurement of horizontal ground deformations with accuracy better than 1/10 of the pixel size. The irregular resampling problem is addressed to avoid introducing aliasing in the orthorectified images. Image registration and correlation is achieved with an iterative, unbiased processor that estimates the phase plane in the Fourier domain for sub-pixel shift detection. Errors due to the imaging system are calibrated and modeled, topography artifacts are characterized and solutions are proposed to compensate or to filter them. A software package implementing these procedures, Co-registration of Optically Sensed Images and Correlation (COSI-Corr), is available from the Caltech Tectonics Observatory website. The procedure is validated in several different contexts, and applied to seismo-tectonics and glaciology studies. Accurate measurements of horizontal co-seismic displacements in the near fault zone allow unambiguous imaging of surface ruptures. It is shown that measurements of surface ruptures from optical aerial and satellite images compare well with field measurements, and that in addition they have the potential of densely measuring the fault perpendicular component, and the off-fault distributed slip. When combined with seismic waveform modeling, fault geometry and surface offsets add crucial constraints to describe in details the seismic faulting process. Dense maps of glacier velocity are reported for several glaciers in Europe and in the Himalayas. Optical image correlation proves robust even in challenging mountainous areas, allowing accurate measurements of glacier flow velocity. Seasonal variations of glacier flow velocity are well identified, suggesting that such measurements can be used to better study the effects of climate change, and to refine the tuning of numerical glacier models.

Assembly of a large earthquake from a complex fault system: Surface rupture kinematics of the 4 April 2010 El Mayor–Cucapah (Mexico) Mw 7.2 earthquake

The 4 April 2010 moment magnitude (M_w) 7.2 El Mayor–Cucapah earthquake revealed the existence of a previously unidentified fault system in Mexico that extends ∼120 km from the northern tip of the Gulf of California to the U.S.–Mexico border. The system strikes northwest and is composed of at least seven major faults linked by numerous smaller faults, making this one of the most complex surface ruptures ever documented along the Pacific–North America plate boundary. Rupture propagated bilaterally through three distinct kinematic and geomorphic domains. Southeast of the epicenter, a broad region of distributed fracturing, liquefaction, and discontinuous fault rupture was controlled by a buried, southwest-dipping, dextral-normal fault system that extends ∼53 km across the southern Colorado River delta. Northwest of the epicenter, the sense of vertical slip reverses as rupture propagated through multiple strands of an imbricate stack of east-dipping dextral-normal faults that extend ∼55 km through the Sierra Cucapah. However, some coseismic slip (10–30 cm) was partitioned onto the west-dipping Laguna Salada fault, which extends parallel to the main rupture and defines the western margin of the Sierra Cucapah. In the northernmost domain, rupture terminates on a series of several north-northeast–striking cross-faults with minor offset (<8 cm) that cut uplifted and folded sediments of the northern Colorado River delta in the Yuha Desert. In the Sierra Cucapah, primary rupture occurred on four major faults separated by one fault branch and two accommodation zones. The accommodation zones are distributed in a left-stepping en echelon geometry, such that rupture passed systematically to structurally lower faults. The structurally lowest fault that ruptured in this event is inclined as shallowly as ∼20°. Net surface offsets in the Sierra Cucapah average ∼200 cm, with some reaching 300–400 cm, and rupture kinematics vary greatly along strike. Nonetheless, instantaneous extension directions are consistently oriented ∼085° and the dominant slip direction is ∼310°, which is slightly (∼10°) more westerly than the expected azimuth of relative plate motion, but considerably more oblique to other nearby historical ruptures such as the 1992 Landers earthquake. Complex multifault ruptures are common in the central portion of the Pacific North American plate margin, which is affected by restraining bend tectonics, gravitational potential energy gradients, and the inherently three-dimensional strain of the transtensional and transpressional shear regimes that operate in this region.

The 2001 Mw 7.6 Bhuj earthquake, low fault friction, and the crustal support of plate driving forces in India

We present a source model for the 2001 M_w 7.6 Bhuj earthquake of northwest India. The slip distribution suggests a high stress drop (~35 MPa) and, together with the depth distribution of aftershocks, that the entire crust is seismogenic. We suggest that the active faults have an effective coefficient of friction of ~0.08, which is sufficient for the seismogenic crust to support the majority of the compressive force transmitted through the Indian lithosphere. This model is consistent with the midcrustal depth of the transition from extension to compression beneath the Ganges foreland basin where India underthrusts southern Tibet. If the coefficient of friction were the more traditional value of 0.6, the lithosphere would be required to support a net force roughly an order of magnitude higher than current estimates in order to match the observed depth of the neutral fiber.

Co-Registration of Optically Sensed Images and Correlation (COSI-Corr): an operational methodology for ground deformation measurements

Recent methodological progress, co-registration of optically sensed images and correlation, outlined here, makes it possible to measure horizontal ground deformation from optical images on an operational basis, using the COSI-Corr software package. In particular, its sub-pixel capabilities allow for accurate mapping of surface ruptures and measurement of co-seismic offsets. We retrieved the fault rupture of the 2005 Mw 7.6 Kashmir earthquake from ASTER images, and we also present a dense mapping of the 1992 Mw 7.3 Landers earthquake of California, from the mosaicking of 30 pairs of aerial images.

Erosion by an Alpine glacier

The glacial way of wearing away The rate at which glaciers erode landscapes is an important but poorly constrained relationship. Herman et al. tackle this issue by considering the Franz Josef alpine glacier in New Zealand. The amount of sediment piling up at the edge of the glacier provided erosion rates, whereas remote sensing allowed for simultaneous tracking of glacial motion. The result was a nonlinear relationship, suggesting that fast glaciers are much more effective at gouging landscapes. This could explain the paradox of why long-term erosion rates are so much lower in polar regions with more permanent glaciers. Science, this issue p. 193 Erosion rate and sliding velocity have a nonlinear relationship in the Franz Josef alpine glacier. Assessing the impact of glaciation on Earth’s surface requires understanding glacial erosion processes. Developing erosion theories is challenging because of the complex nature of the erosion processes and the difficulty of examining the ice/bedrock interface of contemporary glaciers. We demonstrate that the glacial erosion rate is proportional to the ice-sliding velocity squared, by quantifying spatial variations in ice-sliding velocity and the erosion rate of a fast-flowing Alpine glacier. The nonlinear behavior implies a high erosion sensitivity to small variations in topographic slope and precipitation. A nonlinear rate law suggests that abrasion may dominate over other erosion processes in fast-flowing glaciers. It may also explain the wide range of observed glacial erosion rates and, in part, the impact of glaciation on mountainous landscapes during the past few million years.

The 2013 Mw 7.7 Balochistan Earthquake: Seismic Potential of an Accretionary Wedge

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 decollement. 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.

Deformation during the 1975–1984 Krafla rifting crisis, NE Iceland, measured from historical optical imagery

We measure the displacement field resulting from the 1975–1984 Krafla rifting crisis, NE Iceland, using optical image correlation. Images are processed using the COSI-Corr software package. Surface extension is accommodated on normal faults and fissures which bound the rift zone, in response to dike injection at depth. Correlation of declassified KH-9 spy and SPOT5 satellite images reveals extension between 1977–2002 (2.5 m average opening over 80 km), while correlation of aerial photos between 1957–1990 provide measurements of the total extension (average 4.3 m opening over 80 km). Our results show ∼8 m of opening immediately north of Krafla caldera, decreasing to 3–4 m at the northern end of the rift. Correlation of aerial photos from 1957–1976 reveal a bi-modal pattern of opening along the rift during the early crisis, which may indicate either two different magma sources located at either end of the rift zone (a similar pattern of opening was observed in the 2005 Afar rift crisis in East Africa), or variations in rock strength along the rift. Our results provide new information on how past dike injection events accommodate long-term plate spreading, as well as providing more details on the Krafla rift crisis. This study also highlights the potential of optical image correlation using inexpensive declassified spy satellite and aerial photos to measure deformation of the Earths surface going back many decades, thus providing a new tool for measuring Earth surface dynamics, e.g. glaciers, landsliding, coastal erosion, volcano monitoring and earthquake studies, when InSAR and GPS data are not available.

Quantifying near‐field and off‐fault deformation patterns of the 1992 Mw 7.3 Landers earthquake

Coseismic surface deformation in large earthquakes is typically measured using field mapping and with a range of geodetic methods (e.g., InSAR, lidar differencing, and GPS). Current methods, however, either fail to capture patterns of near-field coseismic surface deformation or lack preevent data. Consequently, the characteristics of off-fault deformation and the parameters that control it remain poorly understood. We develop a standardized method to fully measure the surface, near-field, coseismic deformation patterns at high resolution using the COSI-Corr program by correlating pairs of aerial photographs taken before and after the 1992 M_w 7.3 Landers earthquake. COSI-Corr offers the advantage of measuring displacement across the entire zone of surface deformation and over a wider aperture than that available to field geologists. For the Landers earthquake, our measured displacements are systematically larger than the field measurements, indicating the presence of off-fault deformation. We show that 46% of the total surface displacement occurred as off-fault deformation, over a mean deformation width of 154 m. The magnitude and width of off-fault deformation along the rupture is primarily controlled by the macroscopic structural complexity of the fault system, with a weak correlation with the type of near-surface materials through which the rupture propagated. Both the magnitude and width of distributed deformation are largest in stepovers, bends, and at the southern termination of the surface rupture. We find that slip along the surface rupture exhibits a consistent degree of variability at all observable length scales and that the slip distribution is self-affine fractal with dimension of 1.56.

Fault kinematics and surface deformation across a releasing bend during the 2010 MW 7.1 Darfield, New Zealand, earthquake revealed by differential LiDAR and cadastral surveying

Dextral slip at the western end of the east-west–striking Greendale fault during the 2010 M_W 7.1 Darfield earthquake transferred onto a northwest-trending segment, across an apparent transtensional zone, here named the Waterford releasing bend. We used detailed surface mapping, differential analysis of pre- and postearthquake light detection and ranging (LiDAR), and property boundary (cadastral) resurveying to produce high-resolution (centimeter-scale) estimates of coseismic ground-surface displacements across the Waterford releasing bend. Our results indicate that the change in orientation on the Greendale fault incorporates elements of a large-scale releasing bend (from the viewpoint of westward motion on the south side of the fault) as well as a smaller-scale restraining stepover (from the viewpoint of southeastward motion on the north side of the fault). These factors result in the Waterford releasing bend exhibiting a decrease in displacement to near zero at the change in strike, and the presence within the overall releasing bend of a nested, localized restraining stepover with contractional bulging. The exceptional detail of surface deformation and kinematics obtained from this contemporary surface-rupture event illustrates the value of multimethod investigations. Our data provide insights into strike-slip fault bend kinematics, and into the potentially subtle but important structures that may be present at bends on historic and prehistoric rupture traces.