Eric O. Lindsey

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.

Line of Sight Displacement from ALOS-2 Interferometry: Mw 7.8 Gorkha Earthquake and Mw 7.3 Aftershock

Interferometric synthetic aperture radar (InSAR) is a key tool for the analysis of displacement and stress changes caused by large crustal earthquakes, particularly in remote areas. A challenge for traditional InSAR has been its limited spatial and temporal coverage especially for very large events, whose dimensions exceed the typical swath width of 70–100 km. This problem is addressed by the ALOS-2 satellite, whose PALSAR-2 instrument operates in ScanSAR mode, enabling a repeat time of 2 weeks and a swath width of 350 km. Here we present InSAR line-of-sight displacement data from ALOS-2/PALSAR-2 observations covering the Mw 7.8 Gorkha, Nepal earthquake and its Mw 7.3 aftershock that were acquired within 1 week of each event. The data are made freely available and we encourage their use in models of the fault slip and associated stress changes. The Mw 7.3 aftershock not only extended the rupture area of the main shock toward the east but also left a 20 km gap where the fault has little or no coseismic slip. We estimate this unslipped fault patch has the potential to generate a Mw 6.9 event.

Geodetic slip rates in the southern San Andreas Fault system: Effects of elastic heterogeneity and fault geometry

[1] We use high resolution interferometric synthetic aperture radar and GPS measurements of crustal motion across the southern San Andreas Fault system to investigate the effects of elastic heterogeneity and fault geometry on inferred slip rates and locking depths. Geodetically measured strain rates are asymmetric with respect to the mapped traces of both the southern San Andreas and San Jacinto faults. Two possibilities have been proposed to explain this observation: large contrasts in crustal rigidity across the faults, or an alternate fault geometry such as a dipping San Andreas fault or a blind segment of the San Jacinto Fault. We evaluate these possibilities using a two-dimensional elastic model accounting for heterogeneous structure computed from the Southern California Earthquake Center crustal velocity model CVM-H 6.3. The results demonstrate that moderate variations in elastic properties of the crust do not produce a significant strain rate asymmetry and have only a minor effect on the inferred slip rates. However, we find that small changes in the location of faults at depth can strongly impact the results. Our preferred model includes a San Andreas Fault dipping northeast at 60 , and two active branches of the San Jacinto fault zone. In this case, we infer nearly equal slip rates of 18 1 and 19 2mm/yr for the San Andreas and San Jacinto fault zones, respectively. These values are in good agreement with geologic measurements representing average slip rates over the last 10–10 years, implying steady long-term motion on these faults.

Localized and distributed creep along the southern San Andreas Fault

We investigate the spatial pattern of surface creep and off-fault deformation along the southern segment of the San Andreas Fault using a combination of multiple interferometric synthetic aperture radar viewing geometries and survey-mode GPS occupations of a dense array crossing the fault. Radar observations from Envisat during the period 2003-2010 were used to separate the pattern of horizontal and vertical motion, providing a high-resolution image of uplift and shallow creep along the fault trace. The data reveal pervasive shallow creep along the southernmost 50 km of the fault. Creep is localized on a well-defined fault trace only in the Mecca Hills and Durmid Hill areas, while elsewhere creep appears to be distributed over a 1-2 km wide zone surrounding the fault. The degree of strain localization is correlated with variations in the local fault strike. Using a two-dimensional boundary element model, we show that stresses resulting from slip on a curved fault can promote or inhibit inelastic failure within the fault zone in a pattern matching the observations. The occurrence of shallow, localized interseismic fault creep within mature fault zones may thus be partly controlled by the local fault geometry and normal stress, with implications for models of fault zone evolution, shallow coseismic slip deficit, and geologic estimates of long-term slip rates.

The mechanism of partial rupture of a locked megathrust: The role of fault morphology

Assessment of seismic hazard relies on estimates of how large an area of a tectonic fault could potentially rupture in a single earthquake. Vital information for these forecasts includes which areas of a fault are locked and how the fault is segmented. Much research has focused on exploring downdip limits to fault rupture from chemical and thermal boundaries, and along-strike barriers from subducted structural features, yet we regularly see only partial rupture of fully locked fault patches that could have ruptured as a whole in a larger earthquake. Here we draw insight into this conundrum from the 25 April 2015 M w 7.8 Gorkha (Nepal) earthquake. We invert geodetic data with a structural model of the Main Himalayan thrust in the region of Kathmandu, Nepal, showing that this event was generated by rupture of a decollement bounded on all sides by more steeply dipping ramps. The morphological bounds explain why the event ruptured only a small piece of a large fully locked seismic gap. We then use dynamic earthquake cycle modeling on the same fault geometry to reveal that such events are predicted by the physics. Depending on the earthquake history and the details of rupture dynamics, however, great earthquakes that rupture the entire seismogenic zone are also possible. These insights from Nepal should be applicable to understanding bounds on earthquake size on megathrusts worldwide.

Geodetic constraints on frictional properties and earthquake hazard in the Imperial Valley, Southern California

We analyze a suite of geodetic observations across the Imperial Fault in southern California that span all parts of the earthquake cycle. Coseismic and postseismic surface slips due to the 1979 M 6.6 Imperial Valley earthquake were recorded with trilateration and alignment surveys by Harsh (1982) and Crook et al. (1982), and interseismic deformation is measured using a combination of multiple interferometric synthetic aperture radar (InSAR)-viewing geometries and continuous and survey-mode GPS. In particular, we combine more than 100 survey-mode GPS velocities with InSAR data from Envisat descending tracks 84 and 356 and ascending tracks 77 and 306 (149 total acquisitions), processed using a persistent scatterers method. The result is a dense map of interseismic velocities across the Imperial Fault and surrounding areas that allows us to evaluate the rate of interseismic loading and along-strike variations in surface creep. We compare available geodetic data to models of the earthquake cycle with rateand state-dependent friction and find that a complete record of the earthquake cycle is required to constrain key fault properties including the rate-dependence parameter (a − b) as a function of depth, the extent of shallow creep, and the recurrence interval of large events. We find that the data are inconsistent with a high (>30 mm/yr) slip rate on the Imperial Fault and investigate the possibility that an extension of the San Jacinto-Superstition Hills Fault system through the town of El Centro may accommodate a significant portion of the slip previously attributed to the Imperial Fault. Models including this additional fault are in better agreement with the available observations, suggesting that the long-term slip rate of the Imperial Fault is lower than previously suggested and that there may be a significant unmapped hazard in the western Imperial Valley.