Kenneth W. Hudnut

Near-field investigations of the Landers earthquake sequence, April to July 1992.

The Landers earthquake, which had a moment magnitude (Mw) of 7.3, was the largest earthquake to strike the contiguous United States in 40 years. This earthquake resulted from the rupture of five major and many minor right-lateral faults near the southern end of the eastern California shear zone, just north of the San Andreas fault. Its Mw 6.1 preshock and Mw 6.2 aftershock had their own aftershocks and foreshocks. Surficial geological observations are consistent with local and far-field seismologic observations of the earthquake. Large surficial offsets (as great as 6 meters) and a relatively short rupture length (85 kilometers) are consistent with seismological calculations of a high stress drop (200 bars), which is in turn consistent with an apparently long recurrence interval for these faults.

DETECTION OF AQUIFER SYSTEM COMPACTION AND LAND SUBSIDENCE USING INTERFEROMETRIC SYNTHETIC APERTURE RADAR, ANTELOPE VALLEY, MOJAVE DESERT, CALIFORNIA

Interferometric synthetic aperture radar (InSAR) has great potential to detect and quantify land subsidence caused by aquifer system compaction. InSAR maps with high spatial detail and resolution of range displacement (±10 mm in change of land surface elevation) were developed for a groundwater basin (∼103 km2) in Antelope Valley, California, using radar data collected from the ERS-1 satellite. These data allow comprehensive comparison between recent (1993–1995) subsidence patterns and those detected historically (1926–1992) by more traditional methods. The changed subsidence patterns are generally compatible with recent shifts in land and water use. The InSAR-detected patterns are generally consistent with predictions based on a coupled model of groundwater flow and aquifer system compaction. The minor inconsistencies may reflect our imperfect knowledge of the distribution and properties of compressible sediments. When used in conjunction with coincident measurements of groundwater levels and other geologic information, InSAR data may be useful for constraining parameter estimates in simulations of aquifer system compaction.

Tectonic contraction across Los Angeles after removal of groundwater pumping effects.

After the 1987 Whittier Narrows and 1994 Northridge earthquakes revealed that blind thrust faults represent a significant threat to metropolitan Los Angeles, a network of 250 continuously recording global positioning system (GPS) stations was deployed to monitor displacements associated with deep slip on both blind and surface faults. Here we augment this GPS data with interferometric synthetic aperture radar imagery to take into account the deformation associated with groundwater pumping and strike-slip faulting. After removing these non-tectonic signals, we are left with 4.4 mm yr-1 of uniaxial contraction across the Los Angeles basin, oriented N 36° E (perpendicular to the major strike-slip faults in the area). This indicates that the contraction is primarily accommodated on thrust faults rather than on the northeast-trending strike–slip faults. We have found that widespread groundwater and oil pumping obscures and in some cases mimics the tectonic signals expected from the blind thrust faults. In the 40-km-long Santa Ana basin, groundwater withdrawal and re-injection produces 12 mm yr-1 of long-term subsidence, accompanied by an unprecedented seasonal oscillation of 55 mm in the vertical direction and 7 mm horizontally.

Slip pulse and resonance of the Kathmandu basin during the 2015 Gorkha earthquake, Nepal

The bigger they are, the harder they fall The magnitude 7.8 Gorkha earthquake hit Nepal on 25 April 2015. The earthquake killed thousands and caused great damage. Galetzka et al. determined how the fault that caused this earthquake ruptured. The rupture showed a smooth slip pulse 20 km wide that moved eastward along the fault over about 6 s. The nature of the rupture limited damage to regular dwellings but generated shaking that collapsed taller structures. Science, this issue p. 1091 Continuous GPS and InSAR measurements record slip on the fault responsible for the 2015 Mw 7.8 Gorkha earthquake in Nepal. Detailed geodetic imaging of earthquake ruptures enhances our understanding of earthquake physics and associated ground shaking. The 25 April 2015 moment magnitude 7.8 earthquake in Gorkha, Nepal was the first large continental megathrust rupture to have occurred beneath a high-rate (5-hertz) Global Positioning System (GPS) network. We used GPS and interferometric synthetic aperture radar data to model the earthquake rupture as a slip pulse ~20 kilometers in width, ~6 seconds in duration, and with a peak sliding velocity of 1.1 meters per second, which propagated toward the Kathmandu basin at ~3.3 kilometers per second over ~140 kilometers. The smooth slip onset, indicating a large (~5-meter) slip-weakening distance, caused moderate ground shaking at high frequencies (>1 hertz; peak ground acceleration, ~16% of Earth’s gravity) and minimized damage to vernacular dwellings. Whole-basin resonance at a period of 4 to 5 seconds caused the collapse of tall structures, including cultural artifacts.

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.

Uplift and subsidence associated with the great Aceh-Andaman earthquake of 2004

Rupture of the Sunda megathrust on 26 December 2004 produced broad regions of uplift and subsidence. We define the pivot line separating these regions as a first step in defining the lateral extent and the downdip limit of rupture during that great M_w ≈ 9.2 earthquake. In the region of the Andaman and Nicobar islands we rely exclusively on the interpretation of satellite imagery and a tidal model. At the southern limit of the great rupture we rely principally on field measurements of emerged coral microatolls. Uplift extends from the middle of Simeulue Island, Sumatra, at ~2.5°N, to Preparis Island, Myanmar (Burma), at ~14.9°N. Thus the rupture is ~1600 km long. The distance from the pivot line to the trench varies appreciably. The northern and western Andaman Islands rose, whereas the southern and eastern portion of the islands subsided. The Nicobar Islands and the west coast of Aceh province, Sumatra, subsided. Tilt at the southern end of the rupture is steep; the distance from 1.5 m of uplift to the pivot line is just 60 km. Our method of using satellite imagery to recognize changes in elevation relative to sea surface height and of using a tidal model to place quantitative bounds on coseismic uplift or subsidence is a novel approach that can be adapted to other forms of remote sensing and can be applied to other subduction zones in tropical regions.

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.

Southern California Permanent GPS Geodetic Array: Continuous measurements of regional crustal deformation between the 1992 Landers and 1994 Northridge earthquakes

The southern California Permanent GPS Geodetic Array (PGGA) was established in 1990 across the Pacific-North America plate boundary to continuously monitor crustal deformation. We describe the development of the array and the time series of daily positions estimated for its first 10 sites in the 19-month period between the June 28, 1992 (Mw=7.3), Landers and January 17, 1994 (Mw=6.7), Northridge earthquakes. We compare displacement rates at four site locations with those reported by Feigl et al. [1993], which were derived from an independent set of Global Positioning System (GPS) and very long baseline interferometry (VLBI) measurements collected over nearly a decade prior to the Landers earthquake. The velocity differences for three sites 65–100 km from the earthquakes epicenter are of order of 3–5 mm/yr and are systematically coupled with the corresponding directions of coseismic displacement. The fourth site, 300 km from the epicenter, shows no significant velocity difference. These observations suggest large-scale postseismic deformation with a relaxation time of at least 800 days. The statistical significance of our observations is complicated by our incomplete knowledge of the noise properties of the two data sets; two possible noise models fit the PGGA data equally well as described in the companion paper by Zhang et al. [this issue]; the pre-Landers data are too sparse and heterogeneous to derive a reliable noise model. Under a fractal white noise model for the PGGA data we find that the velocity differences for all three sites are statistically different at the 99% significance level. A white noise plus flicker noise model results in significance levels of only 94%, 43%, and 88%. Additional investigations of the pre-Landers data, and analysis of longer spans of PGGA data, could have an important effect on the significance of these results and will be addressed in future work.

Slip history of the 2003 San Simeon earthquake constrained by combining 1--Hz GPS, strong motion, and teleseismic data

[1] The slip history of the 2003 San Simeon earthquake is constrained by combining strong motion and teleseismic data, along with GPS static offsets and 1-Hz GPS observations. Comparisons of a 1-Hz GPS time series and a co-located strong motion data are in very good agreement, demonstrating a new application of GPS. The inversion results for this event indicate that the rupture initiated at a depth of 8.5 km and propagated southeastwards with a speed 3.0 km/sec, with rake vectors forming a fan structure around the hypocenter. We obtained a peak slip of 2.8 m and total seismic moment of 6.2 10 18 Nm. We interpret the slip distribution as indicating that the hanging wall rotates relative to the footwall around the hypocenter, in a sense that appears consistent with the shape of the mapped fault trace. INDEX TERMS: 1242 Geodesy and Gravity: Seismic deformations (7205); 7212 Seismology: Earthquake ground motions and engineering; 8123 Tectonophysics: Dynamics, seismotectonics. Citation: Ji, C., K. M. Larson, Y. Tan, K. W. Hudnut, and K. Choi (2004), Slip history of the 2003 San Simeon earthquake constrained by combining 1-Hz GPS, strong motion, and teleseismic data, Geophys. Res. Lett., 31, L17608, doi:10.1029/2004GL020448.

Geomorphic and geologic controls of geohazards induced by Nepal’s 2015 Gorkha earthquake

Nepals quake-driven landslide hazards Large earthquakes can trigger dangerous landslides across a wide geographic region. The 2015 Mw 7.8 Gorhka earthquake near Kathmandu, Nepal, was no exception. Kargal et al. used remote observations to compile a massive catalog of triggered debris flows. The satellite-based observations came from a rapid response team assisting the disaster relief effort. Schwanghart et al. show that Kathmandu escaped the historically catastrophic landslides associated with earthquakes in 1100, 1255, and 1344 C.E. near Nepals second largest city, Pokhara. These two studies underscore the importance of determining slope stability in mountainous, earthquake-prone regions. Science, this issue p. 10.1126/science.aac8353; see also p. 147 Satellite imaging isolated hazard potential for earthquake-triggered landslides after the 2015 Gorkha earthquake in Nepal. INTRODUCTION On 25 April 2015, the Gorkha earthquake [magnitude (M) 7.8] struck Nepal, followed by five aftershocks of ≥M 6.0 until 10 June 2015. The earthquakes killed ~9000 people and severely damaged a 550 by 200 km region in Nepal and neighboring countries. Some mountain villages were completely destroyed, and the remote locations, blocked roads, and landslide-dammed rivers prevented ground access to many areas. RATIONALE Our “Volunteer Group” of scientists from nine nations, motivated by humanitarian needs, focused on satellite-based systematic mapping and analysis of earthquake-induced geohazards. We provided information to relief and recovery officials as emergency operations were occurring, while supported by one of the largest-ever NASA-led campaigns of responsive satellite data acquisitions over a vast disaster zone. Our analysis of geohazards distribution allowed evaluation of geomorphic, tectonic, and lithologic controls on earthquake-induced landsliding, process mechanisms, and hazard process chains, particularly where they affected local populations. RESULTS We mapped 4312 coseismic and postseismic landslides. Their distribution shows positive associations with slope and shaking intensity. The highest areal densities of landslides are developed on the downdropped northern tectonic block, which is likely explained by momentary reduction of the normal stress along planes of weakness during downward acceleration. The two largest shocks bracket the high-density landslide distribution, the largest magnitudes of the surface displacement field, and highest peak ground accelerations (PGAs). Landslides are heavily concentrated where PGA was >0.6g and slope is >30°. Additional controls on landslide occurrence are indicated by their clustering near earthquake epicenters and within specific lithologic units. The product of PGA and the sine of surface slope (defined as the landslide susceptibility index) is a good indicator of where most landslides occurred. A tail of the statistical distributions of landslides extends to low values of the landslide susceptibility index. Slight earthquake shaking affected vulnerable materials hanging on steep slopes—such as ice, snow, and glacial debris—and moderate to strong shaking affected poorly consolidated sediments deposited in low-sloping river valleys, which were already poised near a failure threshold. In the remote Langtang Valley, some of the most concentrated destruction and losses of life outside the Kathmandu Valley were directly due to earthquake-induced landslides and air blasts. Complex seismic wave interactions and wave focusing may have caused ridgetop shattering and landslides near Langtang but reduced direct shaking damage on valley floors and at glacial lakes. CONCLUSION The Gorkha earthquake took a tremendous, tragic toll on human lives and culture. However, fortunately no damaging earthquake-caused glacier lake outburst floods were observed by our satellite analysis. The total number of landslides was far fewer than those generated by comparable earthquakes elsewhere, probably because of a lack of surface ruptures, the concentration of deformation along the subsurface thrust fault at 10 to 15 km depth, and the regional dominance of competent high-grade metamorphic and intrusive igneous rock types. Landslide distribution and effects of a huge landslide. (A) Landslides (purple dots) are concentrated mostly north of the tectonic hinge-line. Also shown are the epicenters of the main shock and largest aftershock. Displacements are from the JAXA ALOS-2 ScanSAR interferogram (21 Feb and 2 May 2015 acquisitions). (B and C) Before-and-after photographs obtained by D. Breashears in Langtang Valley showing complete destruction of a large part of Langtang village by a huge landslide. The Gorkha earthquake (magnitude 7.8) on 25 April 2015 and later aftershocks struck South Asia, killing ~9000 people and damaging a large region. Supported by a large campaign of responsive satellite data acquisitions over the earthquake disaster zone, our team undertook a satellite image survey of the earthquakes’ induced geohazards in Nepal and China and an assessment of the geomorphic, tectonic, and lithologic controls on quake-induced landslides. Timely analysis and communication aided response and recovery and informed decision-makers. We mapped 4312 coseismic and postseismic landslides. We also surveyed 491 glacier lakes for earthquake damage but found only nine landslide-impacted lakes and no visible satellite evidence of outbursts. Landslide densities correlate with slope, peak ground acceleration, surface downdrop, and specific metamorphic lithologies and large plutonic intrusions.