Search

The 7 January 2020 M 6.4 Puerto Rico earthquake, the mainshock of an extended earthquake sequence, triggered significant ground failure. In this study, we detail the ground failure that occurred based largely on a postearthquake field reconnaissance campaign that we conducted. We documented more than 300 landslides, mainly rock falls that were concentrated in areas where peak ground acceleration (PGA) exceeded 30%g; sparse smaller landslides occurred in highly susceptible areas more than 50 km from the epicenter and at PGA values

We present an open-source Python package, lsforce, for performing single-force source inversions of long-period (tens to hundreds of seconds) seismic signals. Although the software is designed primarily for landslides, it can be used for any single-force seismic source. The package allows users to produce estimates of the three-component time series of forces exerted on the Earth by a landslide with postprocessing options to estimate the trajectory of its center of mass. Green’s functions for a user-selected 1D Earth model are obtained automatically from the Incorporated Research Institutions for Seismology Synthetics Engine webservice or can be computed for custom 1D Earth models using Computer Programs in Seismology. lsforce implements the two most commonly used source parameterizations: a fully flexible, high-resolution approach and a more stable but lower-resolution method of overlapping triangle sources. Regularization options include a blended zeroth-, first-, and second-order semiautomated Tikhonov regularization scheme, as well as additional optional constraints on start times, end times, and on the sum of forces. Uncertainty due to data selection can be assessed using either a leave-one-out approach or a modified jackknife technique that randomly excludes subsets of the data for multiple re-inversions. Numerous built-in plotting methods allow for easy quality control and assessment of results. In this article, we briefly outline the theory and methodology, describe our implementation, and demonstrate the usage of lsforce using the well-studied 28 June 2016 Lamplugh rock avalanche in Alaska. Despite the rapidly increasing prevalence of landslide single-force inversions in the landslide and seismology literature over the past decade, to our knowledge this is the first open-source code for performing such inversions.

In the minutes to hours after a major earthquake, such as the recent 2018 Mw 7.1 Anchorage event, the U.S. Geological Survey (USGS) produces a suite of interconnected earthquake products that provides diverse information ranging from basic earthquake source parameters to loss estimates. The 2018 Anchorage earthquake is the first major domestic earthquake to occur since several new USGS products have been developed, thus providing an opportunity to discuss the newly expanded USGS earthquake product suite, its timeliness, performance, and reception. Overall, the products were relatively timely, accurate, well received, and widely used, including by the media, who used information and visualizations from many products to frame their early reporting. One downside of the codependence of multiple products is that reasonable updates to upstream products (e.g., magnitude and source characterization) can result in significant changes to downstream products; this was the case for the Anchorage earthquake. However, the coverage of strong‐motion stations and felt reports was so dense that the ShakeMap and downstream products were relatively insensitive to changes in magnitude or fault‐plane orientation once the ground‐motion data were available. Shaking and loss indicators initially fluctuated in the first hour or two after the earthquake, but they stabilized quickly. To understand how the products are being used and how effectively they are being communicated, we analyze the media coverage of USGS earthquake products. Most references to USGS products occurred within the first 48 hr after the event. The lack of coverage after 48 hr could indicate that longer‐term products addressing what actions the USGS is taking or what early reconnaissance has revealed might be useful for those people wanting additional information about the earthquake.

Investigation of ground failure triggered by the 2018 Mw 7.1 Anchorage earthquake showed that landslides, liquefaction, and ground cracking all occurred and caused significant damage. Shallow rock falls and rock slides were the most abundant types of landslides, but they occurred in smaller numbers than global models that are based on earthquake magnitude predict; this might result from the 2018 earthquake being an intraslab event. Liquefaction was common in alluvial and intertidal areas; ground deformation probably related to liquefaction damaged numerous houses and port facilities in Anchorage. Ground cracking was pervasive near the edges of slopes in hilly areas and caused perhaps the most significant property damage of all types of ground failure. A complex of slump–earth flows was triggered along coastal bluffs in southern Anchorage where slides also occurred in 1964; the 2018 slides involved both mobilization of new landside material and reactivation of parts of the 1964 landslide deposits. Large translational slides that formed during the 1964 Alaska earthquake showed evidence of deformation along pre‐existing failure surfaces but did not reactivate with new net downslope displacement. Modeling suggests that ground motion in 2018 was of insufficient duration and too high frequency to trigger reactivation of the deep landslides.

We demonstrate the value of utilizing broadband synthetic seismograms to assess regional seismically induced landslide hazard. Focusing on a case study of an M w 7.0 Seattle fault earthquake in Seattle, Washington, we computed broadband synthetic seismograms that account for rupture directivity and 3D basin amplification. We then adjusted the computed motions on a fine grid for 1D amplifications based on the site response of typical geologic profiles in Seattle and used these time‐series ground motions to trigger shallow landsliding using the Newmark method. The inclusion of these effects was critical in determining the extent of landsliding triggered. We found that for inertially triggered slope failures modeled by the Newmark method, the ground motions used to simulate landsliding must have broadband frequency content in order to capture the full slope displacement. We applied commonly used simpler methods based on ground‐motion prediction equations for the same scenario and found that they predicted far fewer landslides if only the mean values were used, but far more at the maximum range of the uncertainties, highlighting the danger of using just the mean values for such methods. Our results indicate that landsliding triggered by a large Seattle fault earthquake will be extensive and potentially devastating, causing direct losses and impeding recovery. The high impact of landsliding predicted by this simulation shows that this secondary effect of earthquakes should be studied with as much vigor as other earthquake effects. Online Material: High‐resolution maps of relative seismically induced landslide hazard for an M w 7.0 Seattle fault earthquake for dry and saturated soil conditions.