Jing Liu-Zeng

Convergence rate across the Nepal Himalaya and interseismic coupling on the Main Himalayan Thrust: Implications for seismic hazard

We document geodetic strain across the Nepal Himalaya using GPS times series from 30 stations in Nepal and southern Tibet, in addition to previously published campaign GPS points and leveling data and determine the pattern of interseismic coupling on the Main Himalayan Thrust fault (MHT). The noise on the daily GPS positions is modeled as a combination of white and colored noise, in order to infer secular velocities at the stations with consistent uncertainties. We then locate the pole of rotation of the Indian plate in the ITRF 2005 reference frame at longitude = − 1.34° ± 3.31°, latitude = 51.4° ± 0.3° with an angular velocity of Ω = 0.5029 ± 0.0072°/Myr. The pattern of coupling on the MHT is computed on a fault dipping 10° to the north and whose strike roughly follows the arcuate shape of the Himalaya. The model indicates that the MHT is locked from the surface to a distance of approximately 100 km down dip, corresponding to a depth of 15 to 20 km. In map view, the transition zone between the locked portion of the MHT and the portion which is creeping at the long term slip rate seems to be at the most a few tens of kilometers wide and coincides with the belt of midcrustal microseismicity underneath the Himalaya. According to a previous study based on thermokinematic modeling of thermochronological and thermobarometric data, this transition seems to happen in a zone where the temperature reaches 350°C. The convergence between India and South Tibet proceeds at a rate of 17.8 ± 0.5 mm/yr in central and eastern Nepal and 20.5 ± 1 mm/yr in western Nepal. The moment deficit due to locking of the MHT in the interseismic period accrues at a rate of 6.6 ± 0.4 × 10^(19) Nm/yr on the MHT underneath Nepal. For comparison, the moment released by the seismicity over the past 500 years, including 14 M_W ≥ 7 earthquakes with moment magnitudes up to 8.5, amounts to only 0.9 × 10^(19) Nm/yr, indicating a large deficit of seismic slip over that period or very infrequent large slow slip events. No large slow slip event has been observed however over the 20 years covered by geodetic measurements in the Nepal Himalaya. We discuss the magnitude and return period of M > 8 earthquakes required to balance the long term slip budget on the MHT.

Near-Field Deformation from the El Mayor–Cucapah Earthquake Revealed by Differential LIDAR

Earthquakes from Above Preparing for risks and hazards associated with large earthquakes requires detailed understanding of their mechanical properties. In addition to pinpointing the location and magnitude of earthquakes, postmortem analyses of the extent of rupture and amount of deformation are key quantities, but are not simply available from seismological data alone. Using a type of optical remote sensing, Light Detection and Ranging (LiDAR), Oskin et al. (p. 702) surveyed the surrounding area that ruptured during the 2010 Mw 7.2 El Mayor–Cucapah earthquake in Northern Mexico. Because this area had also been analyzed in 2006, a comparative analysis revealed slip rate and strain release on the shallow fault zone and a number of previously unknown faults. As remote imaging becomes cheaper and more common, differential analyses will continue to provide fault-related deformation data that complements modern seismological networks. Optical remote sensing before and after a large earthquake reveals its rupture dynamics. Large [moment magnitude (Mw) ≥ 7] continental earthquakes often generate complex, multifault ruptures linked by enigmatic zones of distributed deformation. Here, we report the collection and results of a high-resolution (≥nine returns per square meter) airborne light detection and ranging (LIDAR) topographic survey of the 2010 Mw 7.2 El Mayor–Cucapah earthquake that produced a 120-kilometer-long multifault rupture through northernmost Baja California, Mexico. This differential LIDAR survey completely captures an earthquake surface rupture in a sparsely vegetated region with pre-earthquake lower-resolution (5-meter–pixel) LIDAR data. The postevent survey reveals numerous surface ruptures, including previously undocumented blind faults within thick sediments of the Colorado River delta. Differential elevation changes show distributed, kilometer-scale bending strains as large as ~103 microstrains in response to slip along discontinuous faults cutting crystalline bedrock of the Sierra Cucapah.

Tectonic control of Yarlung Tsangpo Gorge revealed by a buried canyon in Southern Tibet

The Himalayan mountains are dissected by some of the deepest and most impressive gorges on Earth. Constraining the interplay between river incision and rock uplift is important for understanding tectonic deformation in this region. We report here the discovery of a deeply incised canyon of the Yarlung Tsangpo River, at the eastern end of the Himalaya, which is now buried under more than 500 meters of sediments. By reconstructing the former valley bottom and dating sediments at the base of the valley fill, we show that steepening of the Tsangpo Gorge started at about 2 million to 2.5 million years ago as a consequence of an increase in rock uplift rates. The high erosion rates within the gorge are therefore a direct consequence of rapid rock uplift. Sediment cores from a buried canyon upstream of the Tsangpo Gorge support a rapid uplift event to explain gorge formation. [Also see Perspective by Whipple] Tibetan gorge avoids a tectonic aneurysm Rapid tectonic uplift was responsible for the immense Tsangpo Gorge on the eastern edge of the Tibetan plateau 2.5 million years ago. Wang et al. found a buried canyon upstream from the gorge along the Yarlung Tsangpo River that began filling with sediments after sudden uplift. Drill cores of the buried canyon sediments show the same river gradient as found downstream of the gorge. The constant river gradient strongly suggests a rapid uplift event created the gorge, rather than river incision as previously believed. Science, this issue p. 978

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.

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.

Focused modern denudation of the Longmen Shan margin, eastern Tibetan Plateau

We use river sediment load data to map the pattern of modern denudation across the Longmen Shan margin of the Tibetan Plateau. Suspended sediment load, with corrections of bed load and solute load contributions, is used to calculate watershed-averaged denudation rates. Decadal erosion is spatially heterogeneous, and seasonally modulated by monsoon flows, which account for 80–90% of the sediment load. Enhanced denudation occurs in a ∼50 km wide band on the hanging wall of the Longmen Shan and Huya fault zones, reaching 0.5–0.8 mm/yr. These rates are similar to kyr-scale rates deduced from cosmogenic 10Be and to Myr-scale rates from low-temperature thermochronology. The sediment flux-derived erosion rates decrease with increasing distance plateauward, to less than 0.05 mm/yr at a distance ∼200 km northwest of the foot of the Longmen Shan. The gradient in precipitation across this margin alone cannot explain the one order of magnitude spatial difference in erosion. Rather, the river sediment load data delineates a zone of relatively rapid denudation around active faults that carry the Longmen Shan in their hanging wall. From the similarity of denudation rates measured over Myr, kyr, and decadal time scales, we propose that erosion of the Longmen Shan margin has approached a flux steady state. The erosional efflux is balanced by advection of rock toward the Longmen margin above the ∼20°NW dipping ramp of the margin-bounding fault. Our results suggest that high amounts of landslide material mobilized by earthquakes such as the Mw 7.9 2008 Wenchuan event are gradually removed by rivers, smoothing sediment flux over time. Our results also suggest that caution should be exercised when interpreting young cooling ages as evidence of the initiation of plateau uplift. Advection of an already high plateau into the belt of higher erosion rate at the Longmen Shan could also give rise to an abrupt cooling history.

Incision rate of the Yellow River in Northeastern Tibet constrained by (10)Be and (26)Al cosmogenic isotope dating of fluvial terraces: implications for catchment evolution and plateau building

Abstract Unlike other large rivers flowing out of Tibet, the Yellow River escapes from the plateau towards the NE crossing no less than five NW–SE striking, actively growing ranges and intervening basins. Thick Plio-Quaternary deposits and fluvial terraces testify to a phase of aggradation and sediment infill up to the average surface elevation (3200–3250 m a.s.l.) of the Gonghe, Guide and Qinghai Lake basins. A set of seven main terraces across the Gonghe Basin suggests progressive down-cutting of the Yellow River carving the 500 m deep Longyang gorge at the basin exit. 10Be and 26Al concentrations in quartz of surface and sub-surface samples of four terraces constrain the timing of incision by determining the burial age of the deposit and the exposure age of its surface. Modelling the depth dependence of the 10Be concentration and the 26Al/10Be ratio allows us to constrain the onset of the ongoing phase of incision to 120–250 ka. These ages suggest long-term incision rates between 2–6 mm a−1. Together with the present morphology of the Yellow River terraces across the Gonghe basin and the Longyang gorge, our results imply rapid river catchment evolution and interaction between river dynamics, tectonic and climate in northeastern Tibet.

Erratum to Surficial Slip and Rupture Geometry on the Beichuan Fault near Hongkou during the Mw 7.9 Wenchuan Earthquake, China

Worldwide occurrence and documentation of reverse-type ruptures are sparse. Near Hongkou, the Wenchuan rupture passes through the broad Baisha River valley and provides excellent opportunities to trace the surface faulting in fine details for 13 km distance, one of the longest continuous sections along the entire rupture. In this paper, we present the results of our mapping of the surface rupture in this reach. Based on the discontinuities in slip and geometry, the rupture was divided into four sections for convenience, from west to east: the Shenxi Gou, the Miaoba, the Gaoyuan, and the Bajiao Miao sections, respectively. The vertical offset is large in the Shenxi Gou and the Bajiao Miao sections, locally reaching 5-6 m in maxima, and generally low in the Miaoba section (1-2 m or less in most places). The slip gradient for vertical offset is generally 10 3 , locally up to 10 1 , similar to that in well-documented strike-slip ruptures. Near Gaoyuan village, the surface rupture con- sists of two subparallel branches, with the northern one exhibiting right-lateral slip with minor southeast-side-up thrusting, while the southern one is almost pure south- east-side-up thrusting. This pattern mimics the incomplete slip-partitioning of oblique thrusting on parallel faults but at a local scale. In addition, the sense of vertical throw on these two strands is opposite to the general northwest-side-up thrusting of the Wenchuan rupture. We propose that it is likely due to the inheritance at shallow depth from the southeast-dipping geological faults, and that old fault zone structures can have a strong effect on the dynamic rupture by guiding the rupture propagation onto paths of preexisting, though locally unfavorable, dipping fault planes. We also discuss the cross-cutting slickenside striations observed near Bajiao Miao, which indicated temporal rake rotation during dynamic rupture, and their geological and seismological implications.

Rupture termination at restraining bends: The last great earthquake on the Altyn Tagh Fault

Strike-slip rupture propagation falters where changes in fault strike increase Coulomb failure stress. Numerical models of this phenomenon offer predictions of rupture extent based on bend geometry, but have not been verified with field data. To test model predictions of rupture barriers, we examine rupture extent along a section of the sinistral Altyn Tagh Fault punctuated by three major double bends. We measure 3–8 m offsets and map >95 km of continuous scarps that define the most recent surface rupture. We document the eastern terminus of this rupture within the Aksay bend, where an undeformed Pleistocene alluvial fan we mapped and dated overlaps the fault. We conclude, based on this geomorphologic evidence, that multiple Holocene ruptures have stopped in the Aksay bend. Our field data validate model predictions of rupture termination at a >18° restraining bend and support use of geometric parameters to define expected earthquake sizes in seismic hazard models.

Dating the incision of the Yangtze River gorge at the First Bend using three‐nuclide burial ages

Incision of the Yangtze River gorge is widely interpreted as evidence for lower crustal flow beneath the southeast margin of the Tibetan Plateau. Previous work focused on the onset of incision, but the duration of incision remains unknown. Here we present cosmogenic nuclide burial ages of sediments collected from caves on the walls of the gorge that show the gorge was incised ~1 km sometime between 18 and 9 Ma. Thereafter, incision slowed substantially. We resolve middle Miocene burial ages by using three nuclides and accounting for in situ muogenic production. This approach explains the absolute concentrations of 10Be, 26Al, and 21Ne, as well as 26Al/10Be and 21Ne/10Be ratios. A declining incision rate challenges existing geodynamic interpretations by suggesting that either (1) surface uplift has ceased immediately south of the plateau margin or (2) gorge incision is not a useful proxy for the timing of surface uplift.