Geology topics

Charles Mueller

Publications and source records attributed to Charles Mueller.

At least 19 recordsLinked to original sources

USGS2024

The 2023 Alaska National Seismic Hazard Model

US Geological Survey (USGS) National Seismic Hazard Models (NSHMs) are used extensively for seismic design regulations in the United States and earthquake scenario development, as well as risk assessment and mitigation for both buildings and infrastructure. This 2023 update of the long-term, time-independent Alaska NSHM includes substantial changes to both the earthquake rupture forecast (ERF) and ground motion models (GMMs). The ERF includes numerous additions to the finite-fault model, considers two deformation models, and introduces updated declustering and smoothing algorithms in the gridded background seismicity model. For the Alaska–Aleutian subduction zone, megathrust earthquakes occur on an updated structural and segmentation model, and the moment magnitude (M) 8+ rupture and rate model include a logic tree branch that considers slip rates derived from geodetic models of interface coupling. The megathrust model considers multiple models of down-dip width, and magnitudes are computed using newly developed scaling relations. For subduction intraslab events and subduction interface events with M < 7, the 2023 update uses a smoothed seismicity model with rupture depths derived from Slab2. The 2023 model updates GMMs in all tectonic settings using the recently published Next Generation Attenuation Subduction (NGA-Sub) GMMs for subduction interface and intraslab events, and the NGA-West2 GMMs for active crustal settings. Collectively, additions and updates to the Alaska NSHM result in hazard increases across most of south-central Alaska relative to the previous model, published in 2007. These changes are primarily due to the adoption of updated rate models for the large-magnitude interface events and the NGA-Sub GMMs that have much higher aleatory variability (sigma), consistent with global observations, and that include models of epistemic uncertainty.

Alaska↗

USGS2023

The 2023 US 50-State National Seismic Hazard Model: Overview and implications

The US National Seismic Hazard Model (NSHM) was updated in 2023 for all 50 states using new science on seismicity, fault ruptures, ground motions, and probabilistic techniques to produce a standard of practice for public policy and other engineering applications (defined for return periods greater than ∼475 or less than ∼10,000 years). Changes in 2023 time-independent seismic hazard (both increases and decreases compared to previous NSHMs) are substantial because the new model considers more data and updated earthquake rupture forecasts and ground-motion components. In developing the 2023 model, we tried to apply best available or applicable science based on advice of co-authors, more than 50 reviewers, and hundreds of hazard scientists and end-users, who attended public workshops and provided technical inputs. The hazard assessment incorporates new catalogs, declustering algorithms, gridded seismicity models, magnitude-scaling equations, fault-based structural and deformation models, multi-fault earthquake rupture forecast models, semi-empirical and simulation-based ground-motion models, and site amplification models conditioned on shear-wave velocities of the upper 30 m of soil and deeper sedimentary basin structures. Seismic hazard calculations yield hazard curves at hundreds of thousands of sites, ground-motion maps, uniform-hazard response spectra, and disaggregations developed for pseudo-spectral accelerations at 21 oscillator periods and two peak parameters, Modified Mercalli Intensity, and 8 site classes required by building codes and other public policy applications. Tests show the new model is consistent with past ShakeMap intensity observations. Sensitivity and uncertainty assessments ensure resulting ground motions are compatible with known hazard information and highlight the range and causes of variability in ground motions. We produce several impact products including building seismic design criteria, intensity maps, planning scenarios, and engineering risk assessments showing the potential physical and social impacts. These applications provide a basis for assessing, planning, and mitigating the effects of future earthquakes.

Earthquake Spectra↗

USGS2023

The USGS 2023 Conterminous U.S. time‐independent earthquake rupture forecast

We present the 2023 U.S. Geological Survey time‐independent earthquake rupture forecast for the conterminous United States, which gives authoritative estimates of the magnitude, location, and time‐averaged frequency of potentially damaging earthquakes throughout the region. In addition to updating virtually all model components, a major focus has been to provide a better representation of epistemic uncertainties. For example, we have improved the representation of multifault ruptures, both in terms of allowing more and less fault connectivity than in the previous models, and in sweeping over a broader range of viable models. An unprecedented level of diagnostic information has been provided for assessing the model, and the development was overseen by a 19‐member participatory review panel. Although we believe the new model embodies significant improvements and represents the best available science, we also discuss potential model limitations, including the applicability of logic tree branch weights with respect different types of hazard and risk metrics. Future improvements are also discussed, with deformation model enhancements being particularly worthy of pursuit, as well as better representation of sampling errors in the gridded seismicity components. We also plan to add time‐dependent components, and assess implications with a wider range of hazard and risk metrics.

Bulletin of the Seismological Society of America↗

USGS2023

Induced seismicity and its impact on existing seismic hazard analysis

We develop a scheme for mapping changes in earthquake rates within a region in near-real-time. A specific goal of the work is to track recent changes in the rates of induced earthquakes in the central and eastern United States. We map rates in a time window of interest, map rates in a preceding time window, and then ratio the two maps to show changes. A proof-of-concept map is first prepared, comparing rates during the first six months of 2010 with rates from the preceding five years; this demonstration map shows, among other things, the growth of induced seismicity in Oklahoma during 2010. We then present a series of ten ratio maps: the first six months of 2018 compared with the preceding five years, and so on in six-month increments through the end of 2022. These maps show changes in the rates and locations of induced earthquakes, as well as other seismicity trends. Map regions, time windows, and other model parameters are easily adaptable for other applications.

Technical Letter Report↗

USGS2021

2021 U.S. National Seismic Hazard Model for the State of Hawaii

The 2021 U.S. National Seismic Hazard Model (NSHM) for the State of Hawaii updates the two-decades-old former model by incorporating new data and modeling techniques to improve the underlying ground shaking forecasts of tectonic-fault, tectonic-flexure, volcanic, and caldera collapse earthquakes. Two earthquake ground shaking hazard models (public policy and research) are produced that differ in how they account for declustered catalogs. The earthquake source model is based on (1) declustered earthquake catalogs smoothed with adaptive smoothing methods, (2) earthquake rate forecasts based on three temporally varying 60-year time periods, (3) maximum magnitude models that extend to larger earthquakes than previously considered, (4) a separate Kīlauea-specific seismogenic caldera collapse model which accounts for clustered event behavior observed during the 2018 eruption, and (5) fault ruptures that consider historic seismicity, GPS-based strain rates, and a new Quaternary fault database. Two new Hawaii-specific ground motion models (GMMs) and five additional global models consistent with Hawaii shaking data are used to forecast ground shaking at 23 spectral periods and peak parameters. Site effects are modeled using western U.S. and Hawaii specific empirical equations and provide shaking forecasts for eight site classes. For most sites the new model results in similar spectral accelerations as those in the 2001 NSHM, with a few exceptions caused mostly by GMM changes. Ground motions are highest in the southern portion of the Island of Hawai‘i due to high rates of forecasted earthquakes on décollement faults. Shaking decays to the northwest where lower earthquake rates result from flexure of the tectonic plate. Large epistemic uncertainties in source characterizations and GMMs lead to an overall high uncertainty (more than a factor of 3) in ground shaking at Honolulu and Hilo. The new shaking model indicates significant chances of slight or greater damaging ground motions across most of the island chain.

Hawai'i↗

USGS2021

Changing impacts of Alaska-Aleutian subduction zone tsunamis in California under future sea-level rise

The amplification of coastal hazards such as distant-source tsunamis under future relative sea-level rise (RSLR) is poorly constrained. In southern California, the Alaska-Aleutian subduction zone has been identified as an earthquake source region of particular concern for a worst-case scenario distant-source tsunami. Here, we explore how RSLR over the next century will influence future maximum nearshore tsunami heights (MNTH) at the Ports of Los Angeles and Long Beach. Earthquake and tsunami modeling combined with local probabilistic RSLR projections show the increased potential for more frequent, relatively low magnitude earthquakes to produce distant-source tsunamis that exceed historically observed MNTH. By 2100, under RSLR projections for a high-emissions representative concentration pathway (RCP8.5), the earthquake magnitude required to produce >1 m MNTH falls from ~M w 9.1 (required today) to M w 8.0, a magnitude that is ~6.7 times more frequent along the Alaska-Aleutian subduction zone.

Alaska, California↗

USGS2020

Evaluation of ground‐motion models for U.S. Geological Survey seismic hazard forecasts: Hawaii tectonic earthquakes and volcanic eruptions

The selection and weighting of ground‐motion models (GMMs) introduces a significant source of uncertainty in U.S. Geological Survey (USGS) National Seismic Hazard Modeling Project (NSHMP) forecasts. In this study, we evaluate 18 candidate GMMs using instrumental ground‐motion observations of horizontal peak ground acceleration (PGA) and 5%‐damped pseudospectral acceleration (0.02–10 s) for tectonic earthquakes and volcanic eruptions, to inform logic‐tree weights for the update of the USGS seismic hazard model for Hawaii. GMMs are evaluated using two methods. The first is a total residual visualization approach that compares the probability density function (PDF), mean and standard deviations σ "> σ ⁠ , of the observed and predicted ground motion. The second GMM evaluation method we use is the common total residual probabilistic scoring method (log likelihood [LLH]). The LLH method provides a single score that can be used to weight GMMs in the Hawaii seismic hazard model logic trees. The total residual PDF approach provides additional information by preserving GMM over‐ and underprediction across a broad spectrum of periods that is not available from a single value LLH score. We apply these GMM evaluation methods to two different data sets: (1) a database of instrumental ground motions from historic earthquakes in Hawaii from 1973 to 2007 ( ⁠ M w "> M w 4–7.3) and (2) available ground motions from recent earthquakes ( ⁠ M w "> M w 4–6.9) associated with 2018 Kilauea eruptions. The 2018 Kilauea sequence contains both volcanic eruptions and tectonic earthquakes allowing for statistically significant GMM comparisons of the two event classes. The Kilauea ground observations provide an independent data set allowing us to evaluate the predictive power of GMMs implemented in the new USGS nshmp‐haz software system. We evaluate GMM performance as a function of earthquake depth and we demonstrate that short‐period volcanic eruption ground motions are not well predicted by any candidate GMMs. Nine of the initial 18 candidate GMMs fit the observed ground motions and meet established criteria for inclusion in the update of the Hawaii seismic hazard model. A weighted mean of four top performing GMMs in this study (NGAsubslab, NGAsubinter, ASK14, A10) is 50% lower for PGA than for GMMS used in the previous USGS seismic hazard model for Hawaii.

Hawaii↗

USGS2020

2018 and 2023 U. S. National Seismic Hazard Models

During 2017-2018, the National Seismic Hazard Model (NSHM) was updated by incorporating (1) new median ground motion models, new estimates of their epistemic uncertainties and aleatory variabilities, and new soil amplification factors for the central and eastern U.S., (2) amplification of long-period ground motions in deep sedimentary basins in the Los Angeles, San Francisco, Seattle, and Salt Lake City areas, (3) an updated seismicity catalog, which includes new earthquakes that occurred between 2012 and 2017, and (4) improved computer code and implementation details. Results show significantly increased ground shaking in many (but not all) locations across the central and eastern U.S., and increased ground shaking in the four urban areas (listed above) that overlie deep sedimentary basins in the western U.S. During 2019-2023, the NSHM will consider additional updates for the 2023 NSHM including: 3D simulations in urban areas of Los Angeles and Seattle, additional soil amplification models for other urban areas (e.g., eastern and southern coastal plains of the U.S., Reno, and Las Vegas), non-ergodic ground motion models that reconsider epistemic and aleatory uncertainty, topographic amplification, and new geological and geodetic models for faults. These maps will be considered by the Building Seismic Safety Council, American Society of Civil Engineers, and International Building Code committees for inclusion in upcoming building codes. Due to population growth, more people live and work in areas of moderate or high seismic hazard than ever before, leading to higher risk of undesirable consequences from future ground shaking.

conterminous United States↗

USGS2019

Evaluation of ground‐motion models for U.S. Geological Survey seismic hazard models: 2018 Anchorage, Alaska, Mw 7.1 subduction zone earthquake sequence

Instrumental ground‐motion recordings from the 2018 Anchorage, Alaska ( ⁠ M w "> M w 7.1), earthquake sequence provide an independent data set allowing us to evaluate the predictive power of ground‐motion models (GMMs) for intraslab earthquakes associated with the Alaska subduction zone. In this study, we evaluate 15 candidate GMMs using instrumental ground‐motion observations of peak ground acceleration and 5% damped pseudospectral acceleration (0.02–10 s) to inform logic‐tree weights for the update of the U.S. Geological Survey seismic hazard model for Alaska. GMMs are evaluated using two methods. The first is a total residual visualization approach that compares the probability density function, mean, and standard deviations σ "> σ of the observed and predicted ground motion. The second GMM evaluation method we use is the common total residual probabilistic scoring method (log likelihood [LLH]). The LLH method provides a single score that can be used to weight GMMs in the Alaska seismic hazard model logic trees. To test logic branches in previous seismic hazard models, we evaluate GMM performance as a function of depth and we demonstrate that some GMMs show improved performance for earthquakes with focal depths greater than 50 km. Ten of the initial 15 candidate GMMs fit the observed ground motions and meet established criteria for inclusion in the next update of the Alaska seismic hazard model.

Alaska↗

USGS2019

The 2018 update of the US National Seismic Hazard Model: Overview of model and implications

During 2017–2018, the National Seismic Hazard Model for the conterminous United States was updated as follows: (1) an updated seismicity catalog was incorporated, which includes new earthquakes that occurred from 2013 to 2017; (2) in the central and eastern United States (CEUS), new ground motion models were updated that incorporate updated median estimates, modified assessments of the associated epistemic uncertainties and aleatory variabilities, and new soil amplification factors; (3) in the western United States (WUS), amplified shaking estimates of long-period ground motions at sites overlying deep sedimentary basins in the Los Angeles, San Francisco, Seattle, and Salt Lake City areas were incorporated; and (4) in the conterminous United States, seismic hazard is calculated for 22 periods (from 0.01 to 10 s) and 8 uniform V S30 maps (ranging from 1500 to 150 m/s). We also include a description of updated computer codes and modeling details. Results show increased ground shaking in many (but not all) locations across the CEUS (up to ~30%), as well as near the four urban areas overlying deep sedimentary basins in the WUS (up to ~50%). Due to population growth and these increased hazard estimates, more people live or work in areas of high or moderate seismic hazard than ever before, leading to higher risk of undesirable consequences from forecasted future ground shaking.

Earthquake Spectra↗

USGS2018

Earthquake catalogs for the USGS National Seismic Hazard Maps

We describe a methodology that has been developed at the U.S. Geological Survey for making earthquake catalogs for seismic hazard analysis and review the status of the catalogs for the conterminous United States. A new catalog is assembled from several pre‐existing catalogs. Uniform moment magnitudes and related parameters for estimating unbiased seismicity rates are calculated. Duplicates, explosions, mining‐related earthquakes, and induced earthquakes are flagged, and the catalog is declustered. Distinct catalogs are made for the central and eastern United States and the western United States.

Seismological Research Letters↗

USGS2018

2018 one‐year seismic hazard forecast for the central and eastern United States from induced and natural earthquakes

This article describes the U.S. Geological Survey (USGS) 2018 one‐year probabilistic seismic hazard forecast for the central and eastern United States from induced and natural earthquakes. For consistency, the updated 2018 forecast is developed using the same probabilistic seismicity‐based methodology as applied in the two previous forecasts. Rates of earthquakes across the United States M ≥ 3.0 "> M ≥ 3.0 grew rapidly between 2008 and 2015 but have steadily declined over the past 3 years, especially in areas of Oklahoma and southern Kansas where fluid injection has decreased. The seismicity pattern in 2017 was complex with earthquakes more spatially dispersed than in the previous years. Some areas of west‐central Oklahoma experienced increased activity rates where industrial activity increased. Earthquake rates in Oklahoma (429 earthquakes of M ≥ 3 "> M ≥ 3 and 4 M ≥ 4 "> M ≥ 4 ), Raton basin (Colorado/New Mexico border, six earthquakes M ≥ 3 "> M ≥ 3 ), and the New Madrid seismic zone (11 earthquakes M ≥ 3 "> M ≥ 3 ) continue to be higher than historical levels. Almost all of these earthquakes occurred within the highest hazard regions of the 2017 forecast. Even though rates declined over the past 3 years, the short‐term hazard for damaging ground shaking across much of Oklahoma remains at high levels due to continuing high rates of smaller earthquakes that are still hundreds of times higher than at any time in the state’s history. Fine details and variability between the 2016–2018 forecasts are obscured by significant uncertainties in the input model. These short‐term hazard levels are similar to active regions in California. During 2017, M ≥ 3 "> M ≥ 3 earthquakes also occurred in or near Ohio, West Virginia, Missouri, Kentucky, Tennessee, Arkansas, Illinois, Oklahoma, Kansas, Colorado, New Mexico, Utah, and Wyoming.

Seismological Research Letters↗

USGS2018

Seismic hazard, risk, and design for South America

We calculate seismic hazard, risk, and design criteria across South America using the latest data, models, and methods to support public officials, scientists, and engineers in earthquake risk mitigation efforts. Updated continental scale seismic hazard models are based on a new seismicity catalog, seismicity rate models, evaluation of earthquake sizes, fault geometry and rate parameters, and ground‐motion models. Resulting probabilistic seismic hazard maps show peak ground acceleration, modified Mercalli intensity, and spectral accelerations at 0.2 and 1 s periods for 2%, 10%, and 50% probabilities of exceedance in 50 yrs. Ground shaking soil amplification at each site is calculated by considering uniform soil that is applied in modern building codes or by applying site‐specific factors based on V S 30 "> V S 30 shear‐wave velocities determined through a simple topographic proxy technique. We use these hazard models in conjunction with the Prompt Assessment of Global Earthquakes for Response (PAGER) model to calculate economic and casualty risk. Risk is computed by incorporating the new hazard values amplified by soil, PAGER fragility/vulnerability equations, and LandScan 2012 estimates of population exposure. We also calculate building design values using the guidelines established in the building code provisions. Resulting hazard and associated risk is high along the northern and western coasts of South America, reaching damaging levels of ground shaking in Chile, western Argentina, western Bolivia, Peru, Ecuador, Colombia, Venezuela, and in localized areas distributed across the rest of the continent where historical earthquakes have occurred. Constructing buildings and other structures to account for strong shaking in these regions of high hazard and risk should mitigate losses and reduce casualties from effects of future earthquake strong ground shaking. National models should be developed by scientists and engineers in each country using the best available science.

Bulletin of the Seismological Society of America↗

USGS2018

Preliminary 2018 national seismic hazard model for the conterminous United States

The 2014 U.S. Geological Survey national seismic hazard model for the conterminous U.S. will be updated in 2018 and 2020 to coincide with the Building Seismic Safety Council’s Project 17 timeline for development of new building code design criteria. The two closely timed updates are planned to allow more time for the Provisions Update Committee to analyze the consequences of the hazard model changes in the design criteria. To prepare the 2018 update we held a workshop (March 7-8, 2018) with scientists and engineers to solicit feedback on the model. The 2018 model will be available for public comment during the summer of 2018. The purpose of this paper is to solicit feedback on the modeling choices and results. The 2018 NSHM considers an updated seismicity catalog and certain key changes in the way ground motions are calculated. First, we implemented new Next Generation Attenuation Relationships for the Central and Eastern North America Region and other published models that allow for the calculation of ground motions at additional periods and site classes in the central and eastern U.S. (CEUS). Second, basin depth terms were implemented in the ground motion models in select regions of the western U.S. (WUS) to account for enhanced long-period ground motions at softer soil sites overlying sedimentary basins. Preliminary results indicate higher ground motions for all periods in parts of the CEUS and for long-periods and soft soils in urban areas overlying sedimentary basins in the WUS.

Conference Paper↗

USGS2017

2017 One‐year seismic‐hazard forecast for the central and eastern United States from induced and natural earthquakes

We produce a one‐year 2017 seismic‐hazard forecast for the central and eastern United States from induced and natural earthquakes that updates the 2016 one‐year forecast; this map is intended to provide information to the public and to facilitate the development of induced seismicity forecasting models, methods, and data. The 2017 hazard model applies the same methodology and input logic tree as the 2016 forecast, but with an updated earthquake catalog. We also evaluate the 2016 seismic‐hazard forecast to improve future assessments. The 2016 forecast indicated high seismic hazard (greater than 1% probability of potentially damaging ground shaking in one year) in five focus areas: Oklahoma–Kansas, the Raton basin (Colorado/New Mexico border), north Texas, north Arkansas, and the New Madrid Seismic Zone. During 2016, several damaging induced earthquakes occurred in Oklahoma within the highest hazard region of the 2016 forecast; all of the 21 moment magnitude ( M ) ≥4 and 3 M ≥5 earthquakes occurred within the highest hazard area in the 2016 forecast. Outside the Oklahoma–Kansas focus area, two earthquakes with M ≥4 occurred near Trinidad, Colorado (in the Raton basin focus area), but no earthquakes with M ≥2.7 were observed in the north Texas or north Arkansas focus areas. Several observations of damaging ground‐shaking levels were also recorded in the highest hazard region of Oklahoma. The 2017 forecasted seismic rates are lower in regions of induced activity due to lower rates of earthquakes in 2016 compared with 2015, which may be related to decreased wastewater injection caused by regulatory actions or by a decrease in unconventional oil and gas production. Nevertheless, the 2017 forecasted hazard is still significantly elevated in Oklahoma compared to the hazard calculated from seismicity before 2009.

Seismological Research Letters↗

USGS2016

A rare moderate‐sized (Mw 4.9) earthquake in Kansas: Rupture process of the Milan, Kansas, earthquake of 12 November 2014 and its relationship to fluid injection

The largest recorded earthquake in Kansas occurred northeast of Milan on 12 November 2014 ( M w 4.9) in a region previously devoid of significant seismic activity. Applying multistation processing to data from local stations, we are able to detail the rupture process and rupture geometry of the mainshock, identify the causative fault plane, and delineate the expansion and extent of the subsequent seismic activity. The earthquake followed rapid increases of fluid injection by multiple wastewater injection wells in the vicinity of the fault. The source parameters and behavior of the Milan earthquake and foreshock–aftershock sequence are similar to characteristics of other earthquakes induced by wastewater injection into permeable formations overlying crystalline basement. This earthquake also provides an opportunity to test the empirical relation that uses felt area to estimate moment magnitude for historical earthquakes for Kansas.

Seismological Research Letters↗

USGS2016

Seismic‐hazard forecast for 2016 including induced and natural earthquakes in the central and eastern United States

The U.S. Geological Survey (USGS) has produced a one‐year (2016) probabilistic seismic‐hazard assessment for the central and eastern United States (CEUS) that includes contributions from both induced and natural earthquakes that are constructed with probabilistic methods using alternative data and inputs. This hazard assessment builds on our 2016 final model ( Petersen et al. , 2016 ) by adding sensitivity studies, illustrating hazard in new ways, incorporating new population data, and discussing potential improvements. The model considers short‐term seismic activity rates (primarily 2014–2015) and assumes that the activity rates will remain stationary over short time intervals. The final model considers different ways of categorizing induced and natural earthquakes by incorporating two equally weighted earthquake rate submodels that are composed of alternative earthquake inputs for catalog duration, smoothing parameters, maximum magnitudes, and ground‐motion models. These alternatives represent uncertainties on how we calculate earthquake occurrence and the diversity of opinion within the science community. In this article, we also test sensitivity to the minimum moment magnitude between M 4 and M 4.7 and the choice of applying a declustered catalog with b =1.0 rather than the full catalog with b =1.3. We incorporate two earthquake rate submodels: in the informed submodel we classify earthquakes as induced or natural, and in the adaptive submodel we do not differentiate. The alternative submodel hazard maps both depict high hazard and these are combined in the final model. Results depict several ground‐shaking measures as well as intensity and include maps showing a high‐hazard level (1% probability of exceedance in 1 year or greater). Ground motions reach 0.6 g horizontal peak ground acceleration (PGA) in north‐central Oklahoma and southern Kansas, and about 0.2 g PGA in the Raton basin of Colorado and New Mexico, in central Arkansas, and in north‐central Texas near Dallas–Fort Worth. The chance of having levels of ground motions corresponding to modified Mercalli intensity (MMI) VI or greater earthquake shaking is 2%–12% per year in north‐central Oklahoma and southern Kansas and New Madrid similar to the chance of damage at sites in high‐hazard portions of California caused by natural earthquakes. Hazard is also significant in the Raton basin of Colorado/New Mexico; north‐central Arkansas; Dallas–Fort Worth, Texas; and in a few other areas. Hazard probabilities are much lower (by about half or more) for exceeding MMI VII or VIII. Hazard is 3‐ to 10‐fold higher near some areas of active‐induced earthquakes than in the 2014 USGS National Seismic Hazard Model (NSHM), which did not consider induced earthquakes. This study in conjunction with the LandScan TM Database (2013) indicates that about 8 million people live in areas of active injection wells that have a greater than 1% chance of experiencing damaging ground shaking (MMI≥VI) in 2016. The final model has high uncertainty, and engineers, regulators, and industry should use these assessments cautiously to make informed decisions on mitigating the potential effects of induced and natural earthquakes.

Seismological Research Letters↗

USGS2016

Likelihood testing of seismicity-based rate forecasts of induced earthquakes in Oklahoma and Kansas

Likelihood testing of induced earthquakes in Oklahoma and Kansas has identified the parameters that optimize the forecasting ability of smoothed seismicity models and quantified the recent temporal stability of the spatial seismicity patterns. Use of the most recent 1-year period of earthquake data and use of 10–20-km smoothing distances produced the greatest likelihood. The likelihood that the locations of January–June 2015 earthquakes were consistent with optimized forecasts decayed with increasing elapsed time between the catalogs used for model development and testing. Likelihood tests with two additional sets of earthquakes from 2014 exhibit a strong sensitivity of the rate of decay to the smoothing distance. Marked reductions in likelihood are caused by the nonstationarity of the induced earthquake locations. Our results indicate a multiple-fold benefit from smoothed seismicity models in developing short-term earthquake rate forecasts for induced earthquakes in Oklahoma and Kansas, relative to the use of seismic source zones.

Kansas, Oklahoma↗