USGS Science⌕ Search

SEARCH · USGS Science

Results for “earthquake hazards”

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.

At least 1,117 records · Page 62Linked to original sources

Findings from a decade of ground motion simulation validation research and a path forward

Simulated ground motions have the potential to advance seismic hazard assessments and structural response analyses, particularly for conditions with limited recorded ground motions such as large magnitude earthquakes at short source-to-site distances. However, rigorous validation of simulated ground motions is needed for hazard analysts, practicing engineers, or regulatory bodies to be confident in their use. A decade ago, validation exercises were mainly limited to comparisons of simulated-to-observed waveforms and median values of spectral accelerations for selected earthquakes. The Southern California Earthquake Center (SCEC) Ground Motion Simulation Validation (GMSV) group was formed to increase coordination between simulation modelers and research engineers with the aim of devising and applying more effective methods for simulation validation. Here, we summarize what has been learned in over a decade of GMSV activities, principally reflecting the views of the SCEC research community but also extending our findings and suggestions for a path forward to broader United States and worldwide simulation validation efforts. We categorize different validation methods according to their approach and the metrics considered. Two general approaches are to compare validation metrics from simulations to those from historical records or to those from semi-empirical models. Validation metrics are categorized into ground motion characteristics and structural responses. We discuss example validation studies that have been impactful in the past decade and suggest future research directions. Key lessons learned are that validation is application-specific, our outreach and dissemination need improvement, and much validation-related research remains unexplored.

Earthquake Spectra↗

Landslide inventories: The essential part of seismic landslide hazard analyses

A detailed and accurate landslide inventory is an essential part of seismic landslide hazard analysis. An ideal inventory would cover the entire area affected by an earthquake and include all of the landslides that are possible to detect down to sizes of 1–5 m in length. The landslides must also be located accurately and mapped as polygons depicting their true shapes. Such mapped landslide distributions can then be used to perform seismic landslide hazard analysis and other quantitative analyses. Detailed inventory maps of landslide triggered by earthquakes began in the early 1960s with the use of aerial photography. In recent years, advances in technology have resulted in the accessibility of satellite imagery with sufficiently high resolution to identify and map all but the smallest of landslides triggered by a seismic event. With this ability to view any area of the globe, we can acquire imagery for any earthquake that triggers significant numbers of landslides. However, a common problem of incomplete coverage of the full distributions of landslides has emerged along with the advent of high resolution satellite imagery.

Engineering Geology↗

Timescales of cumulate mobilization and mixing for the 1868 A.D. eruption of Mauna Loa, Island of Hawai‘i

The deadly 1868 A.D. eruption of Mauna Loa’s lower Southwest Rift Zone (Island of Hawai‘i) included a M7.9 earthquake and associated tsunami and landslides, demonstrating the severe hazards posed by Earth’s largest active subaerial volcano. To better understand the relationship between intense seismic activity, dike emplacement, magma storage, transport histories, and mobilization of olivine cumulates at Mauna Loa, we examine compositional zoning of olivine in the 1868 lava flows. Samples range from basalt (< 10% olivine) to picrite (30–40% olivine). The olivine cargo is heterogeneous (Fo 78.2–89.2 ; forsterite = [Mg/(Mg + Fe) × 100]) but dominated by ~ Fo 89 cores that lie above the Fe-Mg equilibrium field of host glasses. Crystal rims < Fo 80 are due to post-eruptive modification in slow cooling lava flows. Minor element compositions fall within the range of other Mauna Loa olivine erupted in the past 200 years. Olivine crystals exhibit both normal and complex Fo zoning patterns that yield timescales of diffusive re-equilibration that range from 3 to 258 days, with 72% of crystals recording 71 days or less. These timescales correspond to magmatic priming of the summit reservoir system ~ 2 months prior to the eruption and the M7.9 earthquake likely facilitated the transport of the crystal-rich summit-derived magmas downrift shortly prior to eruption. If the recently proposed faster Fe-Mg diffusion coefficient is used, timescales instead range from < 1 day to 25 days, with most recording 1 week or less. In this scenario, most of the olivine zoning would have to have been generated after the M7.9 earthquake perturbed the system.

Hawaii↗

Economic resilience lessons from the ShakeOut earthquake scenario

Following a damaging earthquake, “business interruption” (BI)—reduced production of goods and services—begins and continues long after the ground shaking stops. Economic resilience reduces BI losses by making the best use of the resources available at a given point in time (static resilience) or by speeding recovery through repair and reconstruction (dynamic resilience), in contrast to mitigation that prevents damage in the first place. Economic resilience is an important concept to incorporate into economic loss modeling and in recovery and contingency planning. Economic resilience framework includes the applicability of resilience strategies to production inputs and output, demand- and supply-side effects, inherent and adaptive abilities, and levels of the economy. We use our resilience framework to organize and share strategies that enhance economic resilience, identify overlooked resilience strategies, and present evidence and structure of resilience strategies for economic loss modelers. Numerous resilience strategies are compiled from stakeholder discussions about the ShakeOut Scenario (Jones et. al. 2008). Modeled results of ShakeOut BI sector losses reveal variable effectiveness of resilience strategies for lengthy disruptions caused by fire-damaged buildings and water service outages. Resilience is a complement to mitigation and may, in fact, have cost and all-hazards advantages.

Earthquake Spectra↗

Development of a global seismic risk model

Since 2015 the Global Earthquake Model (GEM) Foundation and its partners have been supporting regional programmes and bilateral collaborations to develop an open global earthquake risk model. These efforts led to the development of a repository of probabilistic seismic hazard models, a global exposure dataset comprising structural and occupancy information regarding the residential, commercial and industrial buildings, and a comprehensive set of fragility and vulnerability functions for the most common building classes. These components were used to estimate probabilistic earthquake risk globally using the OpenQuake-engine, an open-source software for seismic hazard and risk analysis. This model allows estimating a number of risk metrics such as annualized average losses or aggregated losses for particular return periods, which are fundamental to the development and implementation of earthquake risk mitigation measures.

Earthquake Spectra↗

A new amphibious earthquake catalog for the Alaska Peninsula: Training, testing, and applying machine learning earthquake picking in offshore and amphibious environments

The dual seismological developments of large-scale offshore observational networks and high-quality machine-learning earthquake detection and picking allow unprecedented seismicity detection in high hazard subduction zones. We develop a high resolution machine learning enhanced seismicity catalog for the Alaska Peninsula stretch of the Alaska-Aleutian Subduction Zone, which was instrumented with both land and ocean bottom seismometers during the Alaska Amphibious Community Seismic Experiment (AACSE) in 2018–2019. We first use an analyst-checked AACSE catalog to develop new local land and ocean-bottom earthquake waveform training data sets, which we then use to fine-tune the EQTransformer machine learning picking model. In the process, we test multiple training strategies for best picking performance on ocean-bottom waveform data. We evaluate the effect of choice of pre-trained starting model, addition of waveform bandpass filtering, and choice of random metadata split. We find no training strategy to be substantially better than the others. Instead, simply having a local waveform data set to fine-tune the picking model substantially improves performance on both the ocean-bottom and land seismic data from AACSE. We finally apply the separately trained land and ocean-bottom models to continuous seismic data for AACSE and use an iterative catalog refinement method to detect new small-magnitude seismicity. The machine learning enhanced catalog has roughly 3x more earthquakes than the analyst-checked catalog, and the magnitude of completeness is regionally variable but decreases 0.1–0.4 magnitude units. With more small seismicity detected, seismicity patterns mirror prior observations in the region. Detection is especially improved far offshore and in the outer rise.

Alaska↗

Liquefaction timing and post-triggering seismic energy: A comparison of crustal and subduction zone earthquakes

The objective of the study is to assess when liquefaction is triggered in a suite of ground motions following simplified approaches and measure the remaining post-triggering energy content of those ground motions. For liquefaction-induced deformations, current simplified analysis procedures do not directly incorporate temporal effects and rely on peak transient intensity measurements. Liquefaction hazard from short-duration, small to moderate-magnitude (M) earthquakes (M4.5–7.5) is adequately expressed using transient intensity measurements. However, subduction-zone interface earthquakes can have magnitudes greater than 9.0, with ground-motion durations exceeding 300 s. Using 525 ground motions from the NGA-Subduction (NGA-Sub) database for subduction-zone earthquakes with M8.25–9.25, the timing of liquefaction was calculated using cyclic counting procedures by assuming a reference stress condition and incorporating cyclic strengths from laboratory element testing. A complementary analysis was completed using 514 crustal ground motion records from the NGA-West2 database for M6.75–7.75. Several trends were identified during this study. First, liquefaction will likely trigger during the first half of the ground motion duration, independent of the earthquake source type. However, subduction-zone motions have larger post-triggering energy content compared to crustal earthquakes. The findings from this work indicate that accurately predicting liquefaction-induced deformations from subduction-zone earthquakes may be substantially improved by using robust time-based liquefaction analysis procedures.

Conference Paper↗

CRESCENT earthquake dynamic rupture, earthquake cycle, and tsunami code verification platform

Physics-based simulations are critical for understanding natural hazards. The increasing complexity of numerical codes requires benchmark exercises to verify that different computational methods yield consistent results when solving the same governing equations. Here, we present an open-access web platform designed for the verification of earthquake dynamic rupture, seismic cycle, and tsunami simulations. The platform architecture utilizes a modular, serverless backend on Amazon Web Services (AWS) to provide scalable file processing and visualization. A lightweight static web application provides a secure interface for uploading and managing results, while the browser-based data visualization enables interactive analysis of time series and surface grid data. By using structured JavaScript Object Notation (JSON) text files to define benchmark structures, the system remains fully extensible, allowing the addition of new scenarios without modifying the underlying software logic. The platform hosts the "The Tsunami Problem Versions" (TTPV) 1 & 2, two benchmarks for 3D fully coupled earthquake dynamic rupture and tsunami generation, and provides a framework for earthquake cycle models. This community resource aims to build trust in numerical simulations and facilitate long-term collaborative code verification as modeling software continues to evolve.

Seismica↗

Surface fault displacement models for strike-slip faults

Fault displacement models (FDMs) are an essential component of the probabilistic fault displacement hazard analyses (PFDHA), much like ground motion models in the probabilistic seismic hazard analyses for ground motion hazards. In this study, we develop several principal surface FDMs for strike-slip earthquakes. The model development is based on analyses of the new and comprehensive fault displacement database developed as part of the Fault Displacement Hazard Initiative project led by the University of California, Los Angeles. The main objective of our study is to update the FDMs that were developed over a decade ago by the U.S. Geological Survey and California Geological Survey, in which a reference trace was drawn manually, FDMs are fixed-effect models for lateral displacement, displacement on multiple subparallel ruptures is not aggregated, magnitude (M) scaling is linear, and natural logarithm of displacement is assumed to be normally distributed.

Report GIRS↗

Efforts to monitor and characterize the recent increasing seismicity in central Oklahoma

The sharp increase in seismicity over a broad region of central Oklahoma has raised concerns regarding the source of the activity and its potential hazard to local communities and energy-industry infrastructure. Efforts to monitor and characterize the earthquake sequences in central Oklahoma are reviewed. Since early 2010, numerous organizations have deployed temporary portable seismic stations in central Oklahoma to record the evolving seismicity. A multiple-event relocation method is applied to produce a catalog of central Oklahoma earthquakes from late 2009 into early 2015. Regional moment tensor (RMT) source parameters were determined for the largest and best-recorded earthquakes. Combining RMT results with relocated seismicity enabled determination of the length, depth, and style of faulting occurring on reactivated subsurface fault systems. It was found that the majority of earthquakes occur on near-vertical, optimally oriented (northeast-southwest and northwest-southeast) strike-slip faults in the shallow crystalline basement. In 2014, 17 earthquakes occurred with magnitudes of 4 or larger. It is suggested that these recently reactivated fault systems pose the greatest potential hazard to the region.

Oklahoma↗

A big problem for small earthquakes: Benchmarking routine magnitudes and conversion relationships with coda-envelope-derived Mw in southern Kansas and northern Oklahoma

Earthquake magnitudes are widely relied upon measures of earthquake size. Although moment magnitude ( ⁠ M w "> M w Mw ⁠ ) has become the established standard for moderate and large earthquakes, difficulty in reliably measuring seismic moments for small (generally M w &lt; 4 "> M w < 4 Mw<4 ⁠ ) earthquakes has meant that magnitudes for these events remain plagued by a patchwork of inconsistent measurement scales. Because of this, magnitudes of small earthquakes and statistics derived from them can be biased. Furthermore, because small earthquakes are much more numerous than large ones, many applications, such as seismic hazard modeling, depend critically on analysis of events characterized by magnitudes other than M w "> M w Mw ⁠ . To assess this problem, we apply coda envelope analysis to reliably determine moment magnitudes for a case study of small earthquakes from northern Oklahoma and southern Kansas. Not surprisingly, we find significant differences among M L "> M L ML ⁠ , m b L g "> m b L g mbLg ⁠ , and M w "> M w Mw for M ∼2–4 earthquakes examined here. More troublingly, we find that relations designed to convert other magnitudes to M w "> M w Mw ⁠ , which are relied upon for important applications such as seismic hazard analysis, often increase rather than decrease this bias for our dataset. In our case study, we find that converted magnitudes can result in a systematic bias sometimes exceeding 0.5 magnitude units, a difference that typically corresponds to a factor of ∼3 in seismicity rate. Moreover, we find a correspondingly large bias in Gutenberg–Richter b ‐values, controlled primarily by inaccurate magnitude scaling in the conversion relationships. Although this study focuses on a relatively small geographic area, we can expect that similar issues exist with varying severity in other regions. Therefore, magnitudes of small earthquakes and their associated statistics, including seismicity rates and b ‐values, should be treated with caution.

Kansas, Oklahoma↗

A time-dependent probabilistic seismic-hazard model for California

For the purpose of sensitivity testing and illuminating nonconsensus components of time-dependent models, the California Department of Conservation, Division of Mines and Geology (CDMG) has assembled a time-dependent version of its statewide probabilistic seismic hazard (PSH) model for California. The model incorporates available consensus information from within the earth-science community, except for a few faults or fault segments where consensus information is not available. For these latter faults, published information has been incorporated into the model. As in the 1996 CDMG/U.S. Geological Survey (USGS) model, the time-dependent models incorporate three multisegment ruptures: a 1906, an 1857, and a southern San Andreas earthquake. Sensitivity tests are presented to show the effect on hazard and expected damage estimates of (1) intrinsic (aleatory) sigma, (2) multisegment (cascade) vs. independent segment (no cascade) ruptures, and (3) time-dependence vs. time-independence. Results indicate that (1) differences in hazard and expected damage estimates between time-dependent and independent models increase with decreasing intrinsic sigma, (2) differences in hazard and expected damage estimates between full cascading and not cascading are insensitive to intrinsic sigma, (3) differences in hazard increase with increasing return period (decreasing probability of occurrence), and (4) differences in moment-rate budgets increase with decreasing intrinsic sigma and with the degree of cascading, but are within the expected uncertainty in PSH time-dependent modeling and do not always significantly affect hazard and expected damage estimates.

California↗

The current unlikely earthquake hiatus at California’s transform boundary paleoseismic sites

Paleoseismic and historical earthquake records used to quantify earthquake recurrence rates can also be used to test the likelihood of seismically quiescent periods. At principal paleoseismic sites in California on the San Andreas, San Jacinto, Elsinore, and Hayward faults, no ground‐rupturing earthquake has occurred in the last 100 yr, yet this interval is about three times the average interearthquake period for the ensemble of sites. We examine long paleoseismic records from these faults, as they carry most of the transform fault slip on the plate boundary, to see if the current hiatus has any precedent in the last 1000 yr. The selection of sites is designed to sample fault sections unlikely to have ruptured together, so their conditional probabilities of a hiatus can be combined as independent events. We find a 100‐yr hiatus is not predicted by common time‐dependent or time‐independent recurrence models. Paleoearthquake dating uncertainties can allow long open intervals at individual sites or subsets of sites, but do not explain the observed gap in the ensemble. After approximately removing redundancies in the full paleoearthquake record, the time‐independent probability of the current 100‐yr gap is of order 0.3%. This raises several questions. Do we live in a statistically exceptional time? Or does some wide‐scale effect modulate earthquake occurrence among sites over longer timescales? Finally, how should we understand seismic hazard estimates in California if the recurrence models on which they rely seem, at minimum, incomplete? Whether due to a statistical anomaly, some longer‐term modulation of earthquake occurrence, or another cause, our results emphasize that the hiatus of the last century has been exceptional.

California↗

Developing seismogenic source models based on geologic fault data

Calculating seismic hazard usually requires input that includes seismicity associated with known faults, historical earthquake catalogs, geodesy, and models of ground shaking. This paper will address the input generally derived from geologic studies that augment the short historical catalog to predict ground shaking at time scales of tens, hundreds, or thousands of years (e.g., SSHAC 1997). A seismogenic source model, terminology we adopt here for a fault source model, includes explicit three-dimensional faults deemed capable of generating ground motions of engineering significance within a specified time frame of interest. In tectonically active regions of the world, such as near plate boundaries, multiple seismic cycles span a few hundred to a few thousand years. In contrast, in less active regions hundreds of kilometers from the nearest plate boundary, seismic cycles generally are thousands to tens of thousands of years long. Therefore, one should include sources having both longer recurrence intervals and possibly older times of most recent rupture in less active regions of the world rather than restricting the model to include only Holocene faults (i.e., those with evidence of large-magnitude earthquakes in the past 11,500 years) as is the practice in tectonically active regions with high deformation rates. During the past 15 years, our institutions independently developed databases to characterize seismogenic sources based on geologic data at a national scale. Our goal here is to compare the content of these two publicly available seismogenic source models compiled for the primary purpose of supporting seismic hazard calculations by the Istituto Nazionale di Geofisica e Vulcanologia (INGV) and the U.S. Geological Survey (USGS); hereinafter we refer to the two seismogenic source models as INGV and USGS, respectively. This comparison is timely because new initiatives are emerging to characterize seismogenic sources at the continental scale (e.g., SHARE in the Euro-Mediterranean, http://www.share-eu.org/; EMME in the Middle East, http://www.emme-gem.org/) and global scale (e.g., GEM, http://www.globalquakemodel.org/; Anonymous 2008). To some extent, each of these efforts is still trying to resolve the level of optimal detail required for this type of compilation. The comparison we provide defines a common standard for consideration by the international community for future regional and global seismogenic source models by identifying the necessary parameters that capture the essence of geological fault data in order to characterize seismogenic sources. In addition, we inform potential users of differences in our usage of common geological/seismological terms to avoid inappropriate use of the data in our models and provide guidance to convert the data from one model to the other (for detailed instructions, see the electronic supplement to this article). Applying our recommendations will permit probabilistic seismic hazard assessment codes to run seamlessly using either seismogenic source input. The USGS and INGV database schema compare well at a first-level inspection. Both databases contain a set of fields representing generalized fault three-dimensional geometry and additional fields that capture the essence of past earthquake occurrences. Nevertheless, there are important differences. When we further analyze supposedly comparable fields, many are defined differently. These differences would cause anomalous results in hazard prediction if one assumes the values are similarly defined. The data, however, can be made fully compatible using simple transformations.

Seismological Research Letters↗

An Atlas of ShakeMaps and population exposure catalog for earthquake loss modeling

We present an Atlas of ShakeMaps and a catalog of human population exposures to moderate-to-strong ground shaking (EXPO-CAT) for recent historical earthquakes (1973-2007). The common purpose of the Atlas and exposure catalog is to calibrate earthquake loss models to be used in the US Geological Survey's Prompt Assessment of Global Earthquakes for Response (PAGER). The full ShakeMap Atlas currently comprises over 5,600 earthquakes from January 1973 through December 2007, with almost 500 of these maps constrained-to varying degrees-by instrumental ground motions, macroseismic intensity data, community internet intensity observations, and published earthquake rupture models. The catalog of human exposures is derived using current PAGER methodologies. Exposure to discrete levels of shaking intensity is obtained by correlating Atlas ShakeMaps with a global population database. Combining this population exposure dataset with historical earthquake loss data, such as PAGER-CAT, provides a useful resource for calibrating loss methodologies against a systematically-derived set of ShakeMap hazard outputs. We illustrate two example uses for EXPO-CAT; (1) simple objective ranking of country vulnerability to earthquakes, and; (2) the influence of time-of-day on earthquake mortality. In general, we observe that countries in similar geographic regions with similar construction practices tend to cluster spatially in terms of relative vulnerability. We also find little quantitative evidence to suggest that time-of-day is a significant factor in earthquake mortality. Moreover, earthquake mortality appears to be more systematically linked to the population exposed to severe ground shaking (Modified Mercalli Intensity VIII+). Finally, equipped with the full Atlas of ShakeMaps, we merge each of these maps and find the maximum estimated peak ground acceleration at any grid point in the world for the past 35 years. We subsequently compare this "composite ShakeMap" with existing global hazard models, calculating the spatial area of the existing hazard maps exceeded by the combined ShakeMap ground motions. In general, these analyses suggest that existing global, and regional, hazard maps tend to overestimate hazard. Both the Atlas of ShakeMaps and EXPO-CAT have many potential uses for examining earthquake risk and epidemiology. All of the datasets discussed herein are available for download on the PAGER Web page ( http://earthquake.usgs.gov/ eqcenter/pager/prodandref/ ). ?? 2009 Springer Science+Business Media B.V.

Bulletin of Earthquake Engineering↗

SCIGN; new Southern California GPS network advances the study of earthquakes

Southern California is a giant jigsaw puzzle, and scientists are now using GPS satellites to track the pieces. These puzzle pieces are continuously moving, slowly straining the faults in between. That strain is then eventually released in earthquakes. The innovative Southern California Integrated GPS Network (SCIGN) tracks the motions of these pieces over most of southern California with unprecedented precision. This new network greatly improves the ability to assess seismic hazards and quickly measure the larger displacements that occur during and immediatelyafter earthquakes.

Fact Sheet↗

Multi-segment earthquakes and tsunami potential of the Aleutian megathrust

Large to great earthquakes and related tsunamis generated on the Aleutian megathrust produce major hazards for both the area of rupture and heavily populated coastlines around much of the Pacific Ocean. Here we use paleoseismic records preserved in coastal sediments to investigate whether segment boundaries control the largest ruptures or whether in some seismic cycles segments combine to produce earthquakes greater than any observed since instrumented records began. Virtually the entire megathrust has ruptured since AD1900, with four different segments generating earthquakes >M8.0. The largest was the M9.2 great Alaska earthquake of March 1964 that ruptured ???800 km of the eastern segment of the megathrust. The tsunami generated caused fatalities in Alaska and along the coast as far south as California. East of the 1964 zone of deformation, the Yakutat microplate experienced two >M8.0 earthquakes, separated by a week, in September 1899. For the first time, we present evidence that earthquakes ???900 and ???1500 years ago simultaneously ruptured adjacent segments of the Aleutian megathrust and the Yakutat microplate, with a combined area ???15% greater than 1964, giving an earthquake of greater magnitude and increased tsunamigenic potential. ?? 2008 Elsevier Ltd. All rights reserved.

Quaternary Science Reviews↗

Seismic hazard in the Istanbul metropolitan area: A preliminary re-evaluation

In 1999, two destructive earthquakes (M7.4 Kocaeli and M7.2 Duzce) occurred in the north west of Turkey and resulted in major stress-drops on the western segment of the North Anatolian Fault system where it continues under the Marmara Sea. These undersea fault segments were recently explored using bathymetric and reflection surveys. These recent findings helped to reshape the seismotectonic environment of the Marmara basin, which is a perplexing tectonic domain. Based on collected new information, seismic hazard of the Marmara region, particularly Istanbul Metropolitan Area and its vicinity, were re-examined using a probabilistic approach. Two seismic source and alternate recurrence models combined with various indigenous and foreign attenuation relationships were adapted within a logic tree formulation to quantify and project the regional exposure on a set of hazard maps. The hazard maps show the peak horizontal ground acceleration and spectral acceleration at 1.0 s. These acceleration levels were computed for 2 and 10 % probabilities of transcendence in 50 years.

Journal of Earthquake Engineering↗