Gilles Peltzer

Updated repeat orbit interferometry package released

RO1_PAC V2.3, a Repeat Orbit Interferometry package that allows topographic and surface change researchers to apply Interferometric Synthetic Aperture Radar (InSAR) methods, is now freely available to the community InSAR is the synthesis of conventional SAR and interferometry techniques that have been developed over several decades in radio astronomy and radar remote sensing. In recent years, it has opened entirely new application areas for radar in the Earth system sciences, including topographic mapping and geodesy. RO1_PAC, developed primarily to work with European Remote Sensing (ERS) satellite radar data, currently supports ERS-1, ERS-2, and Japanese Earth Resources Satellite (JERS) radar data, and is configurable to work with “strip-mode” data from all existing satellite radar instruments. The first release of RO1_ PAC (V1.0) was made quietly in 2000, and roughly 30 groups in the academic and research community currently use it.

Active thrusting and folding in the Qilian Shan, and decoupling between upper crust and mantle in northeastern Tibet

Fieldwork south of the city of Gaotai (Gansu province, China) shows that active shortening of surface sediments in the foothills of the Yumu Shan, a large fore-mountain of the Qilian Shan, at the northeastern edge of Tibet, involves both overthrusting and flexural-slip folding. North of this mountain, we found and mapped a prominent north-facing thrust scarp that offsets a Holocene fan sloping gently (3.4°) to the north. Part of this scarp appears to be related to the M ≈ 7.5, 180 A.D. earthquake that may have led to the demise of the Han Dynasty city of Luo Tuo Chen, in the Hexi corridor. A set of 10, 100–150 m long profiles measured across this scarp, 3.2 m high on the average, can be made to fit the diffusion-degraded morphology of a surface break related to the 180 A.D. event using a value of about 3.3 m^2/10^3 yr for the mass diffusivity ϰ of fanglomerates in this part of Gansu province. Smaller mountain-facing scarps on a terrace-capped foothill result from bedding slip concurrent with active folding of underlying, steeply northdipping, Plioquaternary sandstone and conglomerate beds. Holocene uplift rates along the Yumu Shan, which is only one of the Qilian Shan ranges, are estimated to be between 0.4 and 1.9 mm/yr, which implies that much of the mountain formed in the Quaternary. The periclinal structure of the Plioquaternary envelope under which the Paleozoic core of the Yumu Shan plunges towards the west suggests that the whole 3200 m high mountain is a basement ramp anticline. Mountains striking parallel to the Yumu Shan, with similar structure and comparable or greater sizes north and south of the Hexi corridor probably also correspond to recent, crustal ramp anticlines. This implies that the wide, mountainous upper crustal wedge making the northeastern edge of the Tibet-Qinghai plateau is detached from the underlying lower crust and upper mantle.

Magnitude of Late Quaternary Left-Lateral Displacements Along the North Edge of Tibet

Images taken by the earth observation satellite SPOT of the Quaternary morphology at 18 sites on the 2000-kilometer-long Altyn Tagh fault at the north edge of Tibet demonstrate that it is outstandingly active. Long-term, left-lateral strike-slip offsets of stream channels, alluvial terrace edges, and glacial moraines along the fault cluster between 100 and 400 meters. The high elevation of the sites, mostly above 4000 meters in the periglacial zone, suggests that most offsets resulted from slip on the fault since the beginning of the Holocene. These data imply that slip rates are 2 to 3 centimeters per year along much of the fault length and support the hypothesis that the continuing penetration of India into Asia forces Tibet rapidly toward the east.

Mobility of continental mantle: Evidence from postseismic geodetic observations following the 1992 Landers earthquake

The crust around the rupture zone of the 1992 Landers earthquake has continued to deform in the years following the earthquake at rates ∼3 times greater than pre-earthquake rates. We use a combination of Global Positioning System (GPS) and synthetic aperture radar (InSAR) data collected during a ∼3-year epoch following the earthquake in order to investigate postseismic mechanisms responsible for the high transient velocities. In order to maximize the potential signal from viscoelastic relaxation we evaluate and model postseismic relaxation following the first few months of documented accelerated deformation. The combination of GPS and InSAR data allows us to establish viscoelastic relaxation of the lower crust and upper mantle as the dominant postseismic process and to discriminate among possible viscoelastic models. The data particularly require the presence of a highly ductile uppermost mantle beneath the central Mojave Domain, with temperature between the wet and dry basalt solidus. This is consistent with independent seismic and geochemical inferences of a regionally warm uppermost mantle. Further consideration of seismic velocity variations in conjunction with faulting patterns within the Mojave Desert suggests that the primary faulting characteristics of the Mojave Desert, namely, the pervasive late Cenozoic deformation within the Eastern California Shear Zone versus the near absence of faults in the Western Mojave Domain, are controlled by the rheology of the uppermost mantle.

Postseismic Rebound in Fault Step-Overs Caused by Pore Fluid Flow

Near-field strain induced by large crustal earthquakes results in changes in pore fluid pressure that dissipate with time and produce surface deformation. Synthetic aperture radar (SAR) interferometry revealed several centimeters of postseismic uplift in pull-apart structures and subsidence in a compressive jog along the Landers, California, 1992 earthquake surface rupture, with a relaxation time of 270 ± 45 days. Such a postseismic rebound may be explained by the transition of the Poissons ratio of the deformed volumes of rock from undrained to drained conditions as pore fluid flow allows pore pressure to return to hydrostatic equilibrium.

Active tectonics in southern Xinjiang, China: Analysis of terrace riser and normal fault scarp degradation along the Hotan-Qira fault system

The northern piedmont of the western Kunlun mountains (Xinjiang, China) is marked at its easternmost extremity, south of the Hotan-Qira oases, by a set of normal faults trending N50E for nearly 70 km. Conspicuous on Landsat and SPOT images, these faults follow the southeastern border of a deep flexural basin and may be related to the subsidence of the Tarim platform loaded by the western Kunlun northward overthrust. The Hotan-Qira normal fault system vertically offsets the piedmont slope by 70 m. Highest fault scarps reach 20 m and often display evidence for recent reactivations about 2 m high. Successive stream entrenchments in uplifted footwalls have formed inset terraces. We have leveled topographic profiles across fault scarps and transverse abandoned terrace risers. The state of degradation of each terrace edge has been characterized by a degradation coefficient τ, derived by comparison with analytical erosion models. Edges of highest abandoned terraces yield a degradation coefficient of 33 ± 4 m^2. Profiles of cumulative fault scarps have been analyzed in a similar way using synthetic profiles generated with a simple incremental fault scarp model. The analysis shows that (1) rate of fault slip remained essentially constant since the aggradation of the piedmont surface and (2) the occurrence of inset terraces was synchronous at all studied sites, suggesting a climate-driven terrace formation. Observation of glacial and periglacial geomorphic features along the northern front of the western Kunlun range indicates that the Qira glaciofluvial fan emplaced after the last glacial maximum, during the retreat of the Kunlun glaciers (12–22 ka). The age of the most developed inset terrace in uplifted valleys is inferred to be 10 ± 3 ka, coeval with humid climate pulses of the last deglaciation. The mass diffusivity constant (k=τ/T, being time B.P.) in the Hotan region is determined to be 3.3 ± 1.4 m^2/10^3 years, consistent with other estimates in similar climatic and geologic environments of western China. These results imply a minimum rate for the Tarim subsidence of 3.5 ± 2 mm/yr. If Western Kunlun overthrusts the Tarim platform on a crustal ramp dipping 40°–45° to the south, it would absorb at least 4.5 ± 3 mm/yr of convergence between western Tibet and Tarim.

Shallow creep on the Haiyuan Fault (Gansu, China) revealed by SAR Interferometry

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

Improving InSAR geodesy using Global Atmospheric Models

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.

Current slip rates on conjugate strike-slip faults in central Tibet using synthetic aperture radar interferometry

This is the published version. Copyright 2006 American Geophysical Union. All Rights Reserved.

Analysis of coseismic surface displacement gradients using radar interferometryc New insights into the Landers earthquake

The map of the coseismic displacement field generated by interferometric processing of synthetic aperture radar lSARr images taken before and after the June 28, 1992, Landers earthquake sequence brings new insights into the nature of deformation caused by these earthquakes. We use the interferometric map generated by Massonnet et al. l1993r to analyze the surface displacement field in the vicinity of the fault trace. Complexities in the fringe pattern near the fault reflect shorthwavelength variations of the surface rupture and slip distribution, and attest to large displacement gradients. Along two sections of the fault, characteristic fringe patterns can be recognized, contrasting in density and direction with patterns observed away from the rupture. In order to understand the observed fringe patterns, we compute synthetic interferograms in three simple cases: l1r rigidhbody rotations about a vertical axis, l2r about a horizontal axis ltiltr, and l3r distributed, simple shear. The orientation and spatial separation of interferometric fringes predicted by these models help constrain nearhfield deformation and rupture parameters. Where the Kickapoo fault connects with the Homestead Valley fault, the interferogram shows a clear pattern of parallel N20°W fringes separated by about 160 m. This pattern and vertical offsets measured along the Kickapoo fault suggest that the block between this fault and the Johnson Valley fault may have been tilted, down to the west. A 5hkm block lifted by l m on one side would be tilted by an angle of 0.01° l190 μradr, producing fringes separated by about 160 m, parallel to the tilt axis. Such a tilt, parallel to a N20°W direction, would account for the gradual, northward increase of the vertical slip component observed along the Kickapoo fault. This tilt may also explain the 1 m of reverse slip observed along the “slip gap” section of the Homestead Valley break. Between the southern end of the Johnson Valley fault and the Eureka Peak fault, where no surface rupture has been mapped, the dense pattern of fringes implies distributed shear, probably resulting from fault slip at depth. The density and direction of the fringes in the gap are consistent with a righthlateral slip of 1.2–3.8 m on a blind fault locked above the depth of 1.5–2 km. Such observations of small wavelength features in the SAR interferogram bring new insights into the nearhfield displacement gradient and thus on response of the uppermost crust to seismic rupture.