Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Applied Computing and Geosciences.”

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.

1,657 records · Page 6Linked to original sources
USGS2026

Remote compositional analyses of space-weathered lunar maria

Visible-to-shortwave infrared (VSWIR) reflectance spectroscopy has revolutionized our understanding of planetary surface compositions. However, space-weathering processes on airless bodies complicate quantitative compositional analyses. Here, we present a framework to isolate the signatures of space weathering in VSWIR spectra of lunar maria by leveraging radiative transfer modeling under the assumptions that (i) a space-weathered target can be expressed as a mixture of fresh and fully space-weathered components and (ii) remaining signatures can be modeled by including agglutinates as an end-member component. We first validate this approach against laboratory spectra of space-weathered Apollo mare soils of known mineral compositions using a probabilistic Markov Chain Monte Carlo implementation of the Hapke radiative transfer model. Second, we illustrate how this approach can be applied to orbital Moon Mineralogy Mapper data. The proposed space-weathering correction workflow for lunar maria could be expanded to other lunar lithologies and applied to existing and future data sets.

Planetary Science Journal↗
USGS2026

Open-source gravity reduction workflows for geothermal resource assessment

Potential-field geophysical data such as gravity can enhance understanding of geothermal resources at all stages of the resource life cycle, including assessment, exploration, development, and monitoring, and at multiple scales, from the reservoir scale to regional scale. However, to make gravity data useful for geothermal resource characterization, several processing steps are required to isolate the effects of density variations in the Earth’s crust to enable the identification of structural features associated with geothermal resources. Although this process is well-established, standard computational implementations for processing gravity data that are FAIR (Findable, Accessible, Interoperable, and Reproduceable) are still lacking. This paper details ongoing efforts at the U.S. Geological Survey (USGS) to develop a standard set of open-source Python tools for gravity data reduction that align with the FAIR principles. This workflow makes use of existing open-source tools for geophysical data processing with the goal of maximizing opportunities for rapid improvements, interoperability, and adaptability to other types of geophysical data.

Conference Paper↗
USGS2025

Chemostratigraphy of the Cretaceous Hue Shale in Arctic Alaska: Exploring paleoceanographic controls on trace element enrichment, organic matter accumulation, and source-rock evolution

We document chemostratigraphy in an outcrop of late Albian to early Campanian (∼103–82 Ma) marine source rocks to better understand paleoenvironmental controls on trace element (TE) enrichment and organic matter accumulation in the distal Colville foreland basin of Arctic Alaska and how those drivers are linked to arc volcanism and successions of Cretaceous oceanographic and climatic biogeochemical events. This unique, 113-m-thick section of Cretaceous Hue Shale deposited during a series of previously undocumented Arctic Cretaceous oceanic anoxic events (Lease et al., 2024) is the only known exposure of thermally immature (0.48–0.52% R o , random vitrinite reflectance) Hue Shale in Arctic Alaska. Strata comprise mainly clay-rich mudstone with elevated total organic carbon (TOC) and hydrogen index values reaching 26.3 wt% (mean = 7.5 wt%) and 689 mg hydrocarbon (HC)/g TOC (mean = 385 mg HC/g TOC), respectively. Maceral composition consists predominantly of fluorescent amorphous organic matter, with abundant brightly fluorescent alginite, including Tasmanites , acritarchs, and Leiosphaeridia . Discrete layers of volcanic ash (preserved as bentonite) are present throughout the section and provide quantitative age control based on U–Pb dates. Chemostratigraphic trends are investigated to advance our understanding of local oceanographic conditions and controls on stratigraphic and temporal heterogeneity of Brookian source rocks. Concurrent sedimentary enrichment in Mo, U, V, Pb, and Cu across the Albian–Cenomanian boundary of the exposed basal gamma-ray zone, may reflect anoxic to euxinic benthic redox conditions favoring organic matter accumulation and preservation. Fluctuating degrees of anoxia-euxinia are inferred throughout the overlying Hue Shale succession, reflected by varying patterns of TE enrichment and TE–TOC covariation. Elevated C org /P molar ratios (>250) across most of the section, with several values exceeding 690, signify that enhanced biological productivity is sustained throughout deposition. Enhanced productivity, recorded by both C org /P and excess Ba, also parallels increases in source rock richness (elevated TOC and S 2 values) during the late Albian–early Cenomanian and late Cenomanian – Turonian. Enhanced productivity and variations in oceanic circulation/stratification likely both drove changes in benthic redox conditions that favored organic carbon accumulation and preservation. Increased continental arc volcanism (e.g., Okhotsk–Chukotka volcanic belt) and High Arctic Large Igneous Province magmatic eruptions throughout the Cretaceous, inferred to have influenced nutrient cycling and local aqueous nutrient availability, also have been invoked as potential drivers of organic carbon burial and source-rock development across the sedimentary sequence. Results presented here document the organic-rich and oil-prone source-rock quality of the Hue Shale in the distal part of the Colville foreland basin and bolster the potential for a Cretaceous petroleum system beneath the eastern North Slope.

Alaska↗
USGS2025

Flood-Inundation Maps of the Current and Jacks Fork Rivers including the Ozark National Scenic Riverways, Southeast Missouri, 2023

Digital flood-inundation maps for a 131.8-mile reach of the Current River and a 44.6-mile reach of the Jacks Fork River, in southeast Missouri, were created by the U.S. Geological Survey (USGS) in cooperation with the Ozark Foothills Regional Planning Commission and the South Central Ozark Council of Governments. The maps also encompass the 134 miles of the Current and Jacks Fork Rivers within the Ozark National Scenic Riverways, which is the first national park area to protect a river system. The flood-inundation maps, which can be accessed through the USGS Flood Inundation Mapping Program website at https://water.usgs.gov/osw/flood_inundation/ , depict estimates of the areal extent and depth of flooding corresponding to selected water levels (hereafter referred to as “stages”) at eight reference USGS streamgages—five on the Current River (USGS station numbers 07064440, 07064533, 07066510, 07067000, and 07068000) and three on the Jacks Fork River (USGS station numbers 07065200, 07065495, and 07066000). Near-real-time stages at these streamgages may be obtained from the USGS National Water Information System at https://doi.org/10.5066/F7P55KJN or the National Weather Service National Water Prediction Service at http://water.noaa.gov/ , which also forecasts flood hydrographs at four of these sites (USGS station numbers 07067000, 07068000, 07065495, and 07066000). Flood profiles were computed for seven of the eight map reaches by means of two-dimensional hydraulic models and the remaining reach by a one-dimensional hydraulic model. The models were calibrated by using stage-streamflow relations or streamflow measurements at the USGS streamgages and from high-flow stage measurements from water-level loggers distributed throughout the reaches. The hydraulic models were used to compute water-surface profiles for flood stages at 1-foot intervals referenced to the streamgage datums. The profile stages ranged from the National Weather Service “action stage” or near bankfull, to a stage exceeding the highest recorded water level at each streamgage. The simulated water-surface profiles were then combined with a digital elevation model (derived from light detection and ranging data having a nonvegetated vertical accuracy of a maximum 10-centimeter root mean square error) to delineate the area flooded at each water level and the associated water depths. The availability of these maps, along with information regarding current stage from the USGS streamgage and forecasted high-flow stages from the National Weather Service, will provide emergency management personnel, resource managers, and residents with information that is critical for flood-response activities such as evacuations and road closures, as well as for postflood recovery efforts.

Missouri↗
USGS2024

Wave propagation in layered soil deposits

Recent advances in the general theory of viscoelastic waves and rays in layered media provide a rigorous mathematical framework for site-specific, soil-response models used for earthquake resistant design. The advances provide general closed-form anelastic solutions for the classic problems of the response of a stack of soil layers to S and P waves, ray theory for reflected and refracted waves, Rayleigh- and Love-Type surface waves, and head waves. These general solutions valid for anelastic media regardless of the amount of material damping yield new insights regarding the characteristics of seismic waves and their ray paths that are not provided by conventional models. They provide corresponding numerical ground-response models and ray-tracing computation algorithms that account for changes in velocity and attenuation of anelastic waves associated with changes in inhomogeneity of the waves induced by anelastic soil and soil-rock boundaries. These anelastic effects manifest themselves as variations in amplitude response, amplitude attenuation, raypath location, and travel time as observed at the Earth’s surface. Implications of these anelastic effects for soil-response models used for earthquake resistant design are provided herein.

Conference Paper↗
USGS2026

Flood-inundation maps of the Blue River and major tributaries in and near Kansas City, Missouri, 2023–25

Digital flood-inundation maps for 35.5 miles of the Blue River, in and near Kansas City, Missouri, were created by the U.S. Geological Survey (USGS) in cooperation with the City of Kansas City. Maps were also developed for a combined 7.3 miles of Brush Creek (4.4 miles) and Indian Creek (2.9 miles), two primary tributaries to the Blue River in the study reach. The flood-inundation maps, available through the USGS Flood Inundation Mapping Program website at https://water.usgs.gov/osw/flood_inundation/ , depict estimates of the areal extent and depth of flooding, depth of water above first-floor elevations in selected buildings, water velocity, water-surface elevations, and time of travel corresponding to selected water levels or gage heights (hereafter referred to as “stages”) at nine reference USGS streamgages—five on the Blue River (USGS station numbers 06893590, 06893578, 06893500, 06893150, and 06893100), two on Brush Creek (USGS station numbers 06893562 and 06893557), one on Indian Creek (USGS station number 06893390), and one on the Missouri River (USGS station number 06893000). Near-real-time stages at these streamgages may be obtained from USGS Water Data for the Nation at https://doi.org/10.5066/F7P55KJN or the National Weather Service National Water Prediction Service at http://water.noaa.gov/ , which also forecasts flood hydrographs at seven of these sites (USGS station numbers 06893000, 06893590, 06893578, 06893500, 06893150, 06893557, and 06893390). Flood profiles were computed for eight map reaches by means of two-dimensional hydraulic models. The models were calibrated using the USGS-developed stage–streamflow relations at each reference streamgage. Two map reaches on the lower Blue River also include the effects of backwater from the Missouri River on flood inundation extent and hydraulic characteristics. The hydraulic models were used to compute water-surface profiles for flood stages at 1-foot intervals referenced to the streamgage datums. The profile stages ranged from the National Weather Service “Action Stage” or near bankfull to a stage exceeding the highest recorded water level at each streamgage. The simulated water-surface profiles were then combined with a digital elevation model (derived from light detection and ranging data having a nonvegetated vertical accuracy of a maximum 10-centimeter root-mean-square error) to delineate the area flooded at each water level and the associated hydraulic characteristics. The availability of these maps, along with information regarding current stage from the USGS streamgage and forecasted high-flow stages from the National Weather Service, will provide emergency management personnel, resource managers, and residents with information that could be critical for flood-response activities such as evacuations and road closures, as well as for postflood recovery efforts.

Kansas, Missouri↗
USGS2026

A diatom-based quantitative sea-ice proxy for the Bering and Chukchi seas

Sea ice affects Earth's climate system on both regional and global scales. Its incorporation into climate can be used to achieve more accurate predictions of future climate. However, instrumental records of sea-ice concentration do not extend earlier than 1978. In an effort to extend this record, we constructed a proxy using the generalized additive model based on relative abundances of five easy-to-identify diatom species found in sediment samples across the Bering and Chukchi seas. Here we present the first quantitative diatom-based sea-ice proxy developed for Beringia. The developed proxy has been applied to two sediment cores in the Bering Sea ranging from 0 to 25.7 ka (HLY0204 51JPC) and 369 to 430 ka (IODP Exp 323 Site U1345) and one in the Chukchi Sea ranging from 2.7 to 10 ka (HLY0204 24JPC). The obtained reconstructions of sea-ice concentrations are similar, but not identical to previously published qualitative and nearby records based on other proxies. Because our results are quantitative, they can be incorporated into regional climate models. The proxy is publicly available as an R Shiny application (app) and can be applied to any diatom count from marine sediments in the region.

Bering Sea, Chukchi Sea↗
USGS2024

Framework for implementing damping scaling factors in U.S. Geological Survey National Seismic Hazard Models

Traditionally, probabilistic seismic hazard analysis (PSHA) has focused on calculating ground motion hazard curves for elastic, 5%-damped pseudo spectral accelerations, Sa(T,5%), which are used as the basis for engineering design parameters and targets for ground motion selection and modification. However, structures and geotechnical systems can exhibit a wide range of damping ratios both above and below the 5% level, depending on the construction material, structural system, nonstructural elements, or subsurface soil properties. When spectral parameters at such damping levels are required for certain applications, 5%-damped accelerations have traditionally been extracted from PSHA-based hazard curves and adjusted outside of the hazard integral using damping scaling factors (DSF) such as those from Newmark & Hall (1982). Recent advances in the development of more rigorous and comprehensive damping scaling models (e.g., Rezaeian et al., 2014; Rezaeian et al., 2021) have allowed for the modeling of means and standard deviations of DSFs as functions of earthquake source and path properties for crustal, intraslab, and subduction interface tectonic environments. These DSF models can be applied to ground motion model (GMM) estimates of Sa(T,5%) for a given earthquake rupture scenario to produce a corresponding mean and standard deviation Sa at a specified damping ratio β, Sa(T,β). In this study, the DSF models of Rezaeian et al. (2014) and Rezaeian et al. (2021) are implemented within the U.S. Geological Survey National Seismic Hazard Model (NSHM) PSHA framework to calculate probabilistic hazard curves for spectral accelerations at damping ratios from 0.5% to 30%. The DSF models are applied directly to the mean and standard deviation of Sa(T,5%) predictions from each GMM in the NSHM logic tree. Resulting hazard curves and uniform hazard and risk spectra for Sa(T,β) are presented for several geographic locations and compared with corresponding spectra estimated using current design practices by applying the same DSFs outside of the PSHA calculation. Key differences between the two methods for estimating Sa(T,β) are discussed, and potential strategies are presented for the implementation and usage of the hazard-consistent Sa(T,β) in building codes. Comparing the results to those from DSFs used in current design practices that are mainly based on Newmark & Hall (1982) is not explored in this study.

Conference Paper↗
USGS2025

Subduction zone earthquake catalog separation tool: Implementation in the USGS 2025 Puerto Rico and U.S. Virgin Islands National Seismic Hazard Model

The U.S. Geological Survey (USGS) periodically releases updates to National Seismic Hazard Model (NSHM) for the United States and its territories leveraging current scientific knowledge and methodologies to guide public policy, building codes, and risk assessments regarding potential ground shaking due to earthquakes that may result in infrastructure damage. In subduction zones, there is a need to separate the earthquake catalog into tectonic regimes to create specific seismicity models for which the most appropriate ground‐motion models are then applied. Here, we describe newly developed methods and software, called CatSep, that classifies subduction zone events into three primary tectonic regimes: crustal, interface, and intraslab. This method incorporates information about the location of the earthquake relative to the subducting slab, the depth of the Mohorovičić discontinuity, and the earthquake’s moment tensor. Applying this method is a first step in the NSHM workflow for regions covering U.S. subduction zones. Results using this subduction zone earthquake catalog separation tool for the 2025 Puerto Rico and U.S. Virgin Islands NSHM earthquake catalog are presented and analyzed.

Puerto Rico, U.S. Virgin Islands↗
USGS2026

Remote sensing enables basin-scale inventories of coal mine methane

Underground coal mines are important global sources of methane, but emission estimates are uncertain. We show that emission estimates for individual mines from aircraft remote-sensing surveys in the United States agree within 40% with direct measurements used for national emission reporting (IPCC Tier 3 estimate). Such direct measurements are unavailable in most countries, which rely on estimated emission factors (EFs) applied to coal-production rates. We find that EFs from IPCC Tier 1 and the Model for Calculating Coal Mine Methane (MC2M) methods overestimate U.S. emissions 3-fold due to incorrect dependence on mine depth. An IPCC Tier 2 method using measured basin-specific mine gas content agrees with direct emission measurements but does not account for gob well emissions and requires gas content data that are generally unavailable. We show that aircraft remote sensing for a small sample of mines can successfully estimate basin-specific EFs for ventilation shafts and gob wells, enabling estimates of basin- and national-scale emissions. We discuss how the method can be applied with satellite remote sensing to quantify coal emissions worldwide.

Alabama, Colorado, Kentucky, New Mexico, Ohio, Pen↗
USGS2023

WellSTIC: A cost-effective sensor for performing point dilution tests to measure groundwater velocity in shallow aquifers

Many individual measurement points are required to characterize groundwater velocity within an aquifer. Groundwater velocity is most commonly measured using a network of >5 cm diameter monitoring wells, which, if not already present at a site, are expensive and labor-intensive to install. Drive-point piezometers—simple, cost-effective wells that can be installed by hand—are a common tool for sampling groundwater in shallow, alluvial aquifers, but most groundwater velocity measurement techniques require equipment that is too large for these narrow (usually <2 cm inside diameter) piezometers. In this technical note, we introduce a low-cost sensor and well packer system (<$90 USD) for performing point dilution tests in narrow piezometers. Field data show that the magnitude of groundwater velocity measured with this technique agrees with velocities computed from natural gradient tracer tests. Additionally, with proper calibration, these sensors can be used to continuously monitor in-well specific conductance, either during inter-well tracer tests with saline tracers or for water quality monitoring. This system is a viable tool for rapid assessment of the magnitude of groundwater velocity in shallow aquifers.

Water Resources Research↗
USGS2023

A comprehensive fault system inversion approach: Methods and application to NSHM23

We present updated inversion‐based fault‐system solutions for the 2023 update to the National Seismic Hazard Model (NSHM23), standardizing earthquake rate model calculations on crustal faults across the western United States. We build upon the inversion methodology used in the Third Uniform California Earthquake Rupture Forecast (UCERF3) to solve for time‐independent rates of earthquakes in an interconnected fault system. The updated model explicitly maps out a wide range of fault recurrence and segmentation behavior (epistemic uncertainty), more completely exploring the solution space of viable models beyond those of UCERF3. We also improve the simulated annealing implementation, greatly increasing computational efficiency (and thus inversion convergence), and introduce an adaptive constraint weight calculation algorithm that helps to mediate between competing constraints. Hazard calculations show that ingredient changes (especially fault and deformation models) are the primary driver of hazard changes between NSHM23 and UCERF3. Updates to the inversion methodology are also consequential near faults in which the slip rate in UCERF3 was poorly fit or was satisfied primarily using large multifault ruptures that are now restricted by explicit b ‐value and segmentation constraints.

Bulletin of the Seismological Society of America↗
USGS2024

Slope Unit Maker (SUMak): An efficient and parameter-free algorithm for delineating slope units to improve landslide modeling

Slope units are terrain partitions bounded by drainage and divide lines. In landslide modeling, including susceptibility modeling and event-specific modeling of landslide occurrence, slope units provide several advantages over gridded units, such as better capturing terrain geometry, improved incorporation of geospatial landslide-occurrence data in different formats (e.g., point and polygon), and better accommodating the varying data accuracy and precision in landslide inventories. However, the use of slope units in regional ( > 100 km 2 ) landslide studies remains limited due, in part, to the large computational costs and/or poor reproducibility with current delineation methods. We introduce a computationally efficient algorithm for the parameter-free delineation of slope units that leverages tools from within TauDEM and GRASS, using an R interface. The algorithm uses geomorphic laws to define the appropriate scaling of the slope units representative of hillslope processes, avoiding the often ambiguous determination of slope unit size. We then demonstrate how slope units enable more robust regional-scale landslide susceptibility and event-specific landslide occurrence maps.

Natural Hazards and Earth Systems Sciences (NHESS)↗
USGS1976

Digital telemetry system network installation for the National Earthquake Information Service

A Digital Telemetry System Network, developed by the Albuquerque Seismological Laboratory, was installed for the USGS National Earthquake Information Service (NEIS) during January 1975. NEIS is located in Golden, Colorado. The initial Digital Telemetry System Network consists of four independent digital telemetry systems that provide wide area coverage over most of the United States. The four systems were installed between: Albuquerque, New Mexico to NEIS; St. Louis, Missouri to NEIS; Salt Lake City, Utah to NEIS; and Menlo Park, California to NEIS. Each digital telemetry system will transmit fourteen channels of high accuracy seismic data in a digital format over a standard C-1 grade telephone line. Digitized data words will have a resolution of eleven digital bits plus one sign bit. Each digital telemetry receiver, located at NEIS, will have the capacity to provide fourteen channels of seismic data in a digital format suitable for either on-line digital computer analysis or magnetic tape recording. In addition the data can be presented in an analog format suitable for recording on helicorders or develocorders. This network is the first large scale digital telemetry system used by the U.S. Geological Survey for the transmission of real time seismic data.

Open-File Report↗
USGS2025

Methodology for compilation of previously published contour data showing the altitude of the base of Dakota Sandstone on the Colorado Plateau, Arizona, Colorado, New Mexico, and Utah

Structure contours and other geologic information from numerous published geologic maps were digitized and compiled into a digital dataset showing the configuration of a single stratigraphic datum, the base of the Dakota Sandstone and its equivalents across the Colorado Plateau. The principal maps compiled in digital form are a series of 1:250,000-scale 1 degree (°) × 2° quadrangle maps published by the U.S. Geological Survey, augmented by other geologic maps published at various map scales. The compiled digital dataset contains geologic map polygons of the Dakota Sandstone and regional stratigraphic equivalents, the location of faults and fold axes, structure contour lines that define the altitude of the base of the unit and bedding orientation data computed from the structure contour lines. This report provides the scientific rationale for compilation of these data and describes the compilation methodology for each of the data elements. This report provides an extended description of the data compilation in a companion U.S. Geological Survey digital data release of spatial data and attributes associated with the contoured surface and associated geologic data layers.

Arizona, Colorado, New Mexico, Utah↗
USGS2026

Transit time modeling framework for predicting freshwater salinization in urban catchments

The salinity of inland freshwaters is rising globally, particularly in urban watersheds where winter road deicers are widely applied. Attributing stream salinity dynamics to specific sources and transport pathways remains challenging due to episodic salt inputs, engineered drainage, and strong coupling between hydrology and subsurface storage. We present a modeling framework that couples climate-driven deicer build-up and wash-off with transient transit time distribution theory to simulate salt transport through drainage, interflow, and groundwater pathways. Applied to an urban watershed in Northern Virginia (USA), the model reproduces ten years of high-frequency stream salinity measurements across daily-to-decadal timescales. The calibrated model implies an average deicer application of 206 tonnes Cl −1 yr −1 , or roughly one 20 kg bag of rock salt person −1 yr −1 when normalized by the 20,000 people living in the watershed. In winter months, higher infiltration routes a large fraction of snowmelt and deicers into shallow subsurface pathways, enhancing vadose-zone and interflow contributions to stream salinity. Limited subsurface storage capacity and seasonal hydrologic turnover flush excess chloride from the vadose zone and groundwater during subsequent summer storms. By linking climate-driven deicer inputs, hydrologic connectivity, and stream water age, the framework provides a transferable basis for diagnosing and managing freshwater salinization in urban watersheds.

Virginia↗
USGS2026

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↗
USGS2024

Uncertainty in ground-motion-to-intensity conversions significantly affects earthquake early warning alert regions

We examine how the choice of ground‐motion‐to‐intensity conversion equations (GMICEs) in earthquake early warning (EEW) systems affects resulting alert regions. We find that existing GMICEs can underestimate observed shaking at short rupture distances or overestimate the extent of low‐intensity shaking. Updated GMICEs that remove these biases would improve the accuracy of alert regions for the ShakeAlert EEW system for the West Coast of the United States. ShakeAlert uses ground‐motion prediction equations (GMPEs), which calculate spatial distributions of peak ground acceleration (PGA) and peak ground velocity (PGV) from earthquake source estimates, combined with GMICEs to translate GMPE output into modified Mercalli intensity (MMI). We find significant epistemic uncertainty in alert distances; near‐source MMI estimates from different GMICEs can differ by over 1 MMI unit, and MMI extents used for public EEW alerts can differ by hundreds of kilometers for larger magnitude earthquakes ( M ∼6.5+). We use a catalog of “Did You Feel It?” shaking reports to evaluate how well GMICEs predict observed shaking. Our preferred GMICE is the one that computes MMI using PGV for high intensities and transitions to using PGA for nondamaging intensities. These results motivate updating GMICE relationships more generally, including in ShakeMap applications.

The Seismic Record↗