Volcano hazards; lessons learned in the eastern Caribbean
Explore the source record for details and available documents.
SEARCH · USGS Science
Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Explore the source record for details and available documents.
The 2018 US Geological Survey National Seismic Hazard Model (NSHM) incorporates new data and updated science to improve the underlying earthquake and ground motion forecasts for the conterminous United States. The NSHM considers many new data and component input models: (1) new earthquakes between 2013 and 2017 and updated earthquake magnitudes for some earlier earthquakes; (2) two updated smoothed seismicity models to forecast earthquake rates; (3) two suites of new central and eastern US (CEUS) ground motion models (GMMs) to translate ground shaking for various earthquake sizes and source-to-site distances considered in the model; (4) two CEUS GMMs for aleatory variability; (5) two CEUS site-effect models that modify ground shaking based on alternative shallow site conditions; (6) more advanced western US (WUS) lithologic and structural information to assess basin site effects for selected urban regions; and (7) a more comprehensive range of outputs (22 periods and 8 site classes) than in previous versions of the NSHMs. Each of these new datasets and models produces changes in the probabilistic ground shaking levels that are spatially and statistically analyzed. Recent earthquakes or changes to some older earthquake magnitudes and locations mostly result in probabilistic ground shaking levels that are similar to previous models, but local changes can reach up to +80% and −60% compared to the 2014 model. Newly developed CEUS models for GMMs, aleatory variability, and site effects cause overall changes up to ±64%. The addition of the WUS basin amplifications causes changes of up to +60% at longer periods for sites overlying deep soft soils. Across the conterminous United States, the hazard changes in the model are mainly caused by new GMMs in the CEUS, by sedimentary basin effects for long periods (≥1 s) in the WUS, and by seismicity changes for short (0.2 s) and long (1 s) periods for both areas.
Earthquakes can impart varying degrees of damage and permanent, inelastic strain on materials, potentially resulting in ruptures that may promote hazards such as landslides and other collapse events. However, the accumulation of damage in rocks under the frequency and amplitude of shaking experienced during earthquake events is rarely systematically measured due to technical limitations. Here, we characterize damage evolution during laboratory experiments on a suite of dacitic rocks from Unzen volcano, Japan, to help resolve accumulated damage and landslide susceptibility of lava domes during regional earthquake events. Damage was imparted during slow (time-dependent creep) and fast (stress-oscillation earthquake simulations) uniaxial loading in compression and tension. Damage evolution is approximated from strain during experiments; all samples accumulate strain during earthquake events, but microfracture-dominated samples tend to be more susceptible to damage than vesicle-dominated samples. The orientation of existing fabrics with respect to loading direction dictates the magnitude of strain accumulation under load oscillations. During each “earthquake” experiment of multiple dynamic stress-oscillations, samples accumulate inelastic strain. The strain imparted during each successive event is initially high and then reduces after 5-7 events, except when stressing results in failure. The strain rate during phases of intermittent stressing tends to be higher than prior to them. Understanding the accumulation of damage and the potential for brittle failure of rocks subjected to earthquakes can help define the origin and timing of certain landslides, rockfalls, lava dome collapses, and other failure events.
The Seattle, Washington, area is known for its livability and its magnificent natural setting. The city and nearby communities are surrounded by an abundance of rivers and lakes and by the bays of Puget Sound. Two majestic mountain ranges, the Olympics and the Cascades, rim the region. These splendid natural features are products of dynamic forces -- landslides, earthquakes, tsunamis, glaciers, volcanoes, and floods. The same processes that formed this beautiful landscape pose hazards to the ever-growing population of the region. To maintain the Seattle area's livability, public and private policymakers must learn to manage the area's vulnerability to natural hazards to protect its three million residents from loss and damage from future disasters. The U.S. Geological Survey (USGS) is working with other Federal and State agencies, the city of Seattle, and other local governments to provide necessary scientific information that will help communities manage the natural hazards. This information will be useful in planning future development, siting public facilities and businesses, and developing effective emergency plans. -- Gori, et.al., 1999
Using the high-performance computing resources of the Southern California Earthquake Center, we simulate broadband (0-10 Hz) ground motions for three Mw 7.8 rupture scenarios of the southern San Andreas fault. The scenarios incorporate a kinematic rupture description with the average rupture speed along the large slip portions of the fault set at 0.96, 0.89, and 0.84 times the local shear wave velocity. Consistent with previous simulations, a southern hypocenter efficiently channels energy into the Los Angeles region along the string of basins south of the San Gabriel Mountains. However, we find the basin ground motion levels are quite sensitive to the prescribed rupture speed, with peak ground velocities at some sites varying by over a factor of two for variations in average rupture speed of about 15%. These results have important implications for estimating seismic hazards in Southern California and emphasize the need for improved understanding of earthquake rupture processes. Copyright 2008 by the American Geophysical Union.
The devastating consequences of the 2004 Indian Ocean and 2011 Tohoku-oki tsunamis have led to increased research into many different aspects of the tsunami phenomenon. In this paper, we review research related to the observed complexity and uncertainty associated with tsunami generation, propagation, and occurrence described and analyzed using a variety of stochastic models. In each case, tsunamis generated by earthquakes are primarily considered. Stochastic models are developed from the physical theories that govern tsunami evolution combined with empirical models fitted to seismic and tsunami observations, as well as tsunami catalogs. These stochastic models are key to providing probabilistic forecasts and hazard assessments for tsunamis. The stochastic methods described here are similar to those described for earthquakes (Vere-Jones, 2013) and volcanoes (Bebbington, 2013) in this Encyclopedia.
Strong aftershocks following major earthquakes present significant challenges for infrastructure recovery as well as for emergency rescue efforts. A tragic instance of this is the 22 February 2011 M w 6.3 Christchurch aftershock in New Zealand, which caused more than 100 deaths while the 2010 M w 7.1 Canterbury mainshock did not cause a single fatality (Figure 1). Therefore, substantial efforts have been directed toward understanding the generation mechanisms of aftershocks as well as mitigating hazards due to aftershocks. Among these efforts are the prediction of strong aftershocks, earthquake early warning, and aftershock probability assessment. Zhang et al. (1999) reported a successful case of strong aftershock prediction with precursory data such as changes in seismicity pattern, variation of b-value, and geomagnetic anomalies. However, official reports of such successful predictions in geophysical journals are extremely rare, implying that deterministic prediction of potentially damaging aftershocks is not necessarily more scientifically feasible than prediction of mainshocks.
In 1912 the Hawaiian Volcano Observatory (HVO) was established by Massachusetts Institute of Technology professor Thomas A. Jaggar Jr. on the island of Hawaii. Driven by the devastation he observed while investigating the volcanic disasters of 1902 at Montagne Pelée in the Caribbean, Jaggar conducted a worldwide search and decided that Hawai‘i provided an excellent natural laboratory for systematic study of earthquake and volcano processes toward better understanding of seismic and volcanic hazards. In the 100 years since HVO’s founding, surveillance and investigation of Hawaiian volcanoes have spurred advances in volcano and seismic monitoring techniques, extended scientists’ understanding of eruptive activity and processes, and contributed to development of global theories about hot spots and mantle plumes.
A new history of great earthquakes (and their tsunamis) for the central and southern Cascadia subduction zone shows more frequent (17 in the past 6700 yr) megathrust ruptures than previous coastal chronologies. The history is based on along-strike correlations of Bayesian age models derived from evaluation of 554 radiocarbon ages that date earthquake evidence at 14 coastal sites. We reconstruct a history that accounts for all dated stratigraphic evidence with the fewest possible ruptures by evaluating the sequence of age models for earthquake or tsunami contacts at each site, comparing the degree of temporal overlap of correlated site age models, considering evidence for closely spaced earthquakes at four sites, and hypothesizing only maximum-length megathrust ruptures. For the past 6700 yr, recurrence for all earthquakes is 370–420 yr. But correlations suggest that ruptures at ∼1.5 ka and ∼1.1 ka were of limited extent (<400 km). If so, post-3-ka recurrence for ruptures extending throughout central and southern Cascadia is 510–540 yr. But the range in the times between earthquakes is large: two instances may be ∼50 yr, whereas the longest are ∼550 and ∼850 yr. The closely spaced ruptures about 1.6 ka may illustrate a pattern common at subduction zones of a long gap ending with a great earthquake rupturing much of the subduction zone, shortly followed by a rupture of more limited extent. The ruptures of limited extent support the continued inclusion of magnitude-8 earthquakes, with longer ruptures near magnitude 9, in assessments of seismic hazard in the region.
Two recent earthquakes left their mark on Santiago de Chile and Tokyo, well beyond the rupture zones, raising questions about the future vulnerability of these and other cities that lie in seismically active regions. Though spared strong shaking, the megacities nevertheless lit up in small quakes, perhaps signaling an abrupt change in the condition for failure on the faults beneath the cities. To detect such changes in earthquake rate requires good seismic monitoring networks; to respond to such hazard increases with civic preparations requires good government.
ASCE 7-10 introduced new seismic design maps that define risk-targeted ground motions such that buildings designed according to these maps will have 1% chance of collapse in 50 years. These maps were developed by iterative risk calculation, wherein a generic building collapse fragility curve is convolved with the U.S. Geological Survey hazard curve until target risk criteria are met. Recent research shows that this current approach may be unconservative at locations where the tectonic environment is much different than that used to develop the generic fragility curve. This study illustrates how risk-targeted ground motions at selected sites would change if generic building fragility curve and hazard assessment were modified to account for seismic risk from subduction earthquakes and near-fault pulses. The paper also explores the difficulties in implementing these changes.
Building a new model, especially one used for policy purposes, takes considerable time, effort, and resources. In justifying such expenditures, one inevitably spends a lot of time denigrating previous models. For example, in pitching the third Uniform California Earthquake Rupture Forecast (UCERF3) ( http://www.WGCEP.org/UCERF3 ), criticisms of the previous model included fault‐segmentation assumptions and the lack of multifault ruptures. In the context of including spatiotemporal clustering for operational earthquake forecasting (e.g., Jordan et al. , 2011 ), another criticism has been that previous candidate models not only ignore elastic rebound but also produce results that are antithetical to that theory. For instance, the short‐term earthquake probabilities model ( Gerstenberger et al. , 2005 ), which provided California aftershock hazard maps at the U.S. Geological Survey web site between 2005 and 2010, implies that the time of highest likelihood for any rupture will be the moment after it occurs, even for a big one on the San Andreas fault. Furthermore, Monte Carlo simulations imply that excluding elastic rebound in such models also produces unrealistic triggering statistics ( Field, 2012 ).
The U.S. Geological Survey recently updated the National Seismic Hazard Maps by incorporating new seismic, geologic, and geodetic information on earthquake rates and associated ground shaking. The 2008 versions supersede those released in 1996 and 2002. These maps are the basis for seismic design provisions of building codes, insurance rate structures, earthquake loss studies, retrofit priorities, and land-use planning. Their use in design of buildings, bridges, highways, and critical infrastructure allows structures to better withstand earthquake shaking, saving lives and reducing disruption to critical activities following a damaging event. The maps also help engineers avoid costs from over-design for unlikely levels of ground motion.
The San Francisco Public Utilities Commission is seismically retrofitting the water delivery system at San Andreas Lake, San Mateo County, California, where the reservoir intake system crosses the San Andreas Fault (SAF). The near-surface fault location and geometry are important considerations in the retrofit effort. Because the SAF trends through highly distorted Franciscan mélange and beneath much of the reservoir, the exact trace of the 1906 surface rupture is difficult to determine from surface mapping at San Andreas Lake. Based on surface mapping, it also is unclear if there are additional fault splays that extend northeast or southwest of the main surface rupture. To better understand the fault structure at San Andreas Lake, the U.S. Geological Survey acquired a series of seismic imaging profiles across the SAF at San Andreas Lake in 2008, 2009, and 2011, when the lake level was near historical lows and the surface traces of the SAF were exposed for the first time in decades. We used multiple seismic methods to locate the main 1906 rupture zone and fault splays within about 100 meters northeast of the main rupture zone. Our seismic observations are internally consistent, and our seismic indicators of faulting generally correlate with fault locations inferred from surface mapping. We also tested the accuracy of our seismic methods by comparing our seismically located faults with surface ruptures mapped by Schussler (1906) immediately after the April 18, 1906 San Francisco earthquake of approximate magnitude 7.9; our seismically determined fault locations were highly accurate. Near the reservoir intake facility at San Andreas Lake, our seismic data indicate the main 1906 surface rupture zone consists of at least three near-surface fault traces. Movement on multiple fault traces can have appreciable engineering significance because, unlike movement on a single strike-slip fault trace, differential movement on multiple fault traces may exert compressive and extensional stresses on built structures within the fault zone. Such differential movement and resulting distortion of built structures appear to have occurred between fault traces at the gatewell near the southern end of San Andreas Lake during the 1906 San Francisco earthquake (Schussler, 1906). In addition to the three fault traces within the main 1906 surface rupture zone, our data indicate at least one additional fault trace (or zone) about 80 meters northeast of the main 1906 surface rupture zone. Because ground shaking also can damage structures, we used fault-zone guided waves to investigate ground shaking within the fault zones relative to ground shaking outside the fault zones. Peak ground velocity (PGV) measurements from our guided-wave study indicate that ground shaking is greater at each of the surface fault traces, varying with the frequency of the seismic data and the wave type (P versus S). S-wave PGV increases by as much as 5–6 times at the fault traces relative to areas outside the fault zone, and P-wave PGV increases by as much as 3–10 times. Assuming shaking increases linearly with increasing earthquake magnitude, these data suggest strong shaking may pose a significant hazard to built structures that extend across the fault traces. Similarly complex fault structures likely underlie other strike-slip faults (such as the Hayward, Calaveras, and Silver Creek Faults) that intersect structures of the water delivery system, and these fault structures similarly should be investigated.
Yucca Mountain is the designated site of the underground repository for the United States' high-level radioactive waste (HLW), consisting of commercial and military spent nuclear fuel, HLW derived from reprocessing of uranium and plutonium, surplus plutonium, and other nuclear-weapons materials. Yucca Mountain straddles the western boundary of the Nevada Test Site, where the United States has tested nuclear devices since the 1950s, and is situated in an arid, remote, and thinly populated region of Nevada, ~100 miles northwest of Las Vegas. Yucca Mountain was originally considered as a potential underground repository of HLW because of its thick units of unsaturated rocks, with the repository horizon being not only ~300 m above the water table but also ~300 m below the Yucca Mountain crest. The fundamental rationale for a geologic (underground) repository for HLW is to securely isolate these materials from the environment and its inhabitants to the greatest extent possible and for very long periods of time. Given the present climate conditions and what is known about the current hydrologic system and conditions around and in the mountain itself, one would anticipate that the rates of infiltration, corrosion, and transport would be very low—except for the possibility that repository integrity might be compromised by low-probability disruptive events, which include earthquakes, strong ground motion, and (or) a repository-piercing volcanic intrusion/eruption. Extreme ground motions (ExGM), as we use the phrase in this report, refer to the extremely large amplitudes of earthquake ground motion that arise at extremely low probabilities of exceedance (hazard). They first came to our attention when the 1998 probabilistic seismic hazard analysis for Yucca Mountain was extended to a hazard level of 10 -8 /yr (a 10-4/yr probability for a 10 4 -year repository “lifetime”). The primary purpose of this report is to summarize the principal results of the ExGM research program as they have developed over the past 5 years; what follows will be focused on Yucca Mountain, but not restricted to it.
Stress levels on and adjacent to megathrust faults at seismogenic depths remain a key but difficult to constrain parameter for assessing seismic hazard in subduction zones. Although strong ground motions have been observed to be generated from distinct, high-stress regions on the downdip end of the megathrust rupture areas in many great earthquakes, we lack direct constraints on the stress level in the lower seismogenic portion of the Cascadia megathrust. On 2022 December 20, a M w 6.4 strike-slip earthquake occurred near Ferndale, California in southern Cascadia and likely ruptured the Gorda slab crust in the lower seismogenic portion, providing an opportunity to assess the stress level in this region. Here, we relocate the Ferndale mainshock and the first two weeks of aftershocks using a high-resolution 3-D velocity model and estimate rupture dimensions, directivity, and stress drop for several M w 4-5 aftershocks and recent earthquakes. The aftershocks define a strike-slip fault in the slab crust striking ENE, consistent with the mainshock focal mechanism. The orientation of this fault is about 45 ° off the ideally oriented fault plane given the stress state in the slab. The aftershock zone is extensive and broad in the forward direction of the mainshock rupture but still constrained within the volume of high Vp/Vs within the slab crust. Our stress drop estimates are generally lower for M w 4-5 earthquakes located within the slab crust compared to those a few km deeper in the slab mantle. Combined, our results support a relatively low effective stress level in the vicinity of the megathrust in the lower portion of the seismogenic zone in southern Cascadia, likely due to elevated fluid pressures. Consequently, the ground motion in the onshore region above this low-stress seismogenic portion in southern Cascadia may not be as intense as that observed during great earthquakes in other subduction zones.
The size of the logic tree within the Uniform California Earthquake Rupture Forecast Version 3, Time-Dependent (UCERF3-TD) model can challenge risk analyses of large portfolios. An insurer or catastrophe risk modeler concerned with losses to a California portfolio might have to evaluate a portfolio 57,600 times to estimate risk in light of the hazard possibility space. Which branches of the logic tree matter most, and which can one ignore? We employed two model-order-reduction techniques to simplify the model. We sought a subset of parameters that must vary, and the specific fixed values for the remaining parameters, to produce approximately the same loss distribution as the original model. The techniques are (1) a tornado-diagram approach we employed previously for UCERF2, and (2) an apparently novel probabilistic sensitivity approach that seems better suited to functions of nominal random variables. The new approach produces a reduced-order model with only 60 of the original 57,600 leaves. One can use the results to reduce computational effort in loss analyses by orders of magnitude.
The seismic potential of the Lesser Antilles megathrust remains poorly known, despite the potential hazard it poses to numerous island populations and its proximity to the Americas. As it has not produced any large earthquakes in the instrumental era, the megathrust is often assumed to be aseismic. However, historical records of great earthquakes in the 19 th century and earlier, which were most likely megathrust ruptures, demonstrate that the subduction is not entirely aseismic. Recent occurrences of giant earthquakes in areas where such events were previously thought to be improbable have illustrated the importance of critically evaluating the seismic potential of other “low-hazard” subduction zones, such as the Lesser Antilles. Using the method of coral microatoll paleogeodesy developed in Sumatra, we examine 20 th -century vertical deformation on the forearc islands of the Lesser Antilles and model the underlying strain accumulation on the megathrust. Our data indicate that the eastern coasts of the forearc islands have been subsiding by up to ∼8 mm/yr relative to sites closer to the arc, suggesting that on the time scale of the 20 th century, a portion of the megathrust just east of the forearc islands has been locked. Our findings are in contrast to recent models based on satellite geodesy that suggest little or no strain accumulation anywhere along the Lesser Antilles megathrust. This discrepancy is potentially explained by the different time scales of measurement, as recent studies elsewhere have indicated that interseismic coupling patterns may vary on decadal time scales and that century-scale or longer records are required to fully assess seismic potential. The accumulated strain we have detected will likely be released in future megathrust earthquakes, uplifting previously subsiding areas and potentially causing widespread damage from strong ground motion and tsunami waves.