SEARCH · Search USGS

Results for “Computers & Geosciences”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

136 records · Page 4Linked to original sources

USGS2026

Ground-motion characterization for the 2025 U.S. National Seismic Hazard Model for Puerto Rico and the U.S. Virgin Islands

We develop the ground-motion characterization (GMC) for the 2025 U.S. National Seismic Hazard Model for Puerto Rico and the U.S. Virgin Islands (NSHM-PRVI) for earthquakes in active crustal, subduction interface, and subduction intraslab regimes. Using ground-motion models (GMMs) from the Next-Generation Attenuation (NGA)-West2 and NGA-Subduction projects, the GMC is parameterized by scaled-backbone models for median ground motions and by independent logic trees of aleatory variability. We introduce several novel GMC features into the U.S. NSHM: (1) use of regional ground-motion data for modeling median ground motions; (2) development of scaled-backbone models for median ground motions; and (3) development of independent logic trees of aleatory variability from variance components of GMMs and computed from multiple ground-motion datasets, and incorporating regional ground-motion data effects on variability. We compute probabilistic seismic hazard curves and maps to evaluate the contributions from the GMC components and for comparison with the 2003 GMC. Contributions of the GMC to epistemic uncertainty in seismic hazard are evaluated through spatial variations in epistemic uncertainty in hazard maps, comparisons between mean hazard curves and fractiles, as well as investigations of the contributions of logic-tree branches to hazard maps and curves. Comparisons between seismic hazard from the 2025 and 2003 GMCs allow examination of the changes in hazard curves and mapped values with 2% and 10% probabilities of exceedance in 50 years. The 2025 GMC exhibits modest changes in median predictions, relative to the 2003 GMC; overall, values of aleatory variability are higher, except at long periods ( ), resulting in changes to probabilistic ground motions at low probabilities of exceedance (<10% probability of exceedance in 50 years). Changes in hazard at 2% and 10% probabilities of exceedance in 50 years are also relatively modest (within 20%) at most sites, with the impacts from the 2025 GMC exhibiting minor reductions and increases at 0.2 and 1.0 s periods, respectively, relative to the 2003 GMC.

Puerto Rico, U.S. Virgin Islands↗

USGS1986

Water-quality appraisal of NASQAN stations below impoundments, eastern Tennessee

The National Stream Quality Accounting Network (NASQAN) is a network of stations at which systematic and continuing water-quality data are collected. Major objectives of this U.S. Geological Survey program are (1) to depict areal variability of streamflow and water-quality conditions nationwide on a year-by-year basis and (2) to detect long-term changes in streamflow and stream quality. Several NASQAN stations in East Tennessee are downstream from impoundments which have a significant effect on water quality. NASQAN data obtained from the Tennessee River below Watts Bar Dam and the Clinch River below Melton Hill Dam were compared to water-quality data from the basins upstream. The comparison indicates that NASQAN data obtained below impoundment may not be adequate to describe a composite picture of water quality in the accounting unit. Detention time of storage in the impoundments is believed to moderate the range of constituent values observed at the NASQAN stations. Data obtained upstream and downstream from Watts Bar Dam indicate that the water sampled at the NASQAN station comes from stratified layers of the impoundment and is not representative of an integrated sample of water from the impoundment. Values of total recoverable iron suggest that, because of adsorption to sediments in impoundments, some constituents are not accurately described by sampling below impoundments. Relations between water-quality constituents and flow at stations on the Clinch River and Tennessee River are not well defined due to regulation. Direct load computations for many constituents were therefore not possible, which diminished the utility of data from these NASQAN stations to account for quantity versus quality of the water. Load computations were only possible for ionic constituents through use of a continuous specific-conductance record as an intermediary. Compensation for the effects of discharge prior to application of the Seasonal Kendall test for trends could not be done and identification of trends in water-quality constituents caused by some process (source) change was not possible. Some water-quality trends indicated by data from the Clinch and Tennessee Rivers might reflect the decreasing trend in discharge during the 1972-82 water years. Thus the stations below Watts Bar Dam and below Melton Hill Dam do not adequately meet the NASQAN objective to detect and assess long-term changes in stream quality.

Tennessee↗

USGS2026

Estimating aftershock risk for entry into earthquake-damaged buildings

We present a simple method to estimate the risk of experiencing strong shaking from aftershocks during entry into earthquake-damaged buildings. We compute wait times until the probability of strong ground shaking from aftershocks reaches a predefined risk threshold; for example, a 0.4 percent probability of experiencing Modified Mercalli Intensity 7 or greater shaking during the planned building entry. We also develop a relation between aftershock probability and the probability of strong shaking, so that users can reference the U.S. Geological Survey aftershock forecast during an ongoing aftershock sequence to determine if the risk threshold has been met. We apply our method to active continental regions (for example, the Western United States), stable continental regions (for example, the Central and Eastern United States), and subduction zones (for example, Cascadia or Alaska).

Open-File Report↗

USGS2024

StreamStats—A quarter century of delivering web-based geospatial and hydrologic information to the public, and lessons learned

StreamStats is a U.S. Geological Survey (USGS) web application that provides streamflow statistics, such as the 1-percent annual exceedance probability peak flow, the mean flow, and the 7-day, 10-year low flow, to the public through a map-based user interface. These statistics are used in many ways, such as in the design of roads, bridges, and other structures; in delineation of floodplains for land-use zoning and setting of insurance rates; for regulatory purposes, such as the permitting of wastewater discharges; and for hydrologic and climate change studies. StreamStats was first developed for Massachusetts and released in 2001. The application provided users with the ability to obtain streamflow statistics computed from data collected at USGS streamgages and to obtain estimates of streamflow statistics for user-selected ungaged sites. Massachusetts StreamStats used geographic information system software and digital mapping to compute drainage-basin characteristics, which were then used in statistical models to estimate streamflow statistics for the user-selected sites. The statistical models were in the form of equations that were developed through a process known as regression analysis. StreamStats was the first known web application with the ability to do interactive geoprocessing. The utility of Massachusetts StreamStats was instantly apparent, leading the USGS to develop a version of StreamStats that could be implemented nationally. USGS State offices normally were required to develop custom regression equations and prepare local digital mapping data needed for implementing StreamStats for their States. Funding needed to complete this work usually was provided through cooperative agreements between the USGS and State agencies. In 2004, Idaho became the first to be released in the national version of StreamStats. By 2023, 44 States were fully implemented and six were undergoing implementation. StreamStats has undergone many modifications over the years to keep up with changes to the underlying software and to add functionality. Customized functionality and separate linked StreamStats applications were developed for several States. Meeting the high demand for additions and improvements to StreamStats while also adhering to budgetary constraints has, at times, been challenging. The StreamStats development team has identified numerous additional improvements that could be made to provide better performance and more functionality. The lessons learned from the experience of building and operating StreamStats for nearly 25 years could be relevant to others interested in pursuing efforts of a similar scale.

Circular↗

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

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↗

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

USGS1990

Recharge rates and aquifer hydraulic characteristics for selected drainage basins in middle and east Tennessee

Quantitative information concerning aquifer hydrologic and hydraulic characteristics is needed to manage the development of ground-water resources. These characteristics are poorly defined for the bedrock aquifers in Middle and East Tennessee where demand for water is increasing. This report presents estimates of recharge rate, storage coefficient, diffusivity, and transmissivity for representative drainage basins in Middle and East Tennessee, as determined from analyses of stream-aquifer interactions. The drainage basins have been grouped according to the underlying major aquifer, then statistical descriptions applied to each group, in order to define area1 distribution of these characteristics. Aquifer recharge rates are estimated for representative low, average, and high flow years for 63 drainage basins using hydrograph analysis techniques. Net annual recharge during average flow years for all basins ranges from 4.1 to 16.8 in/yr (inches per year), with a mean value of 7.3 in. In general, recharge rates are highest for basins underlain by the Blue Ridge aquifer (mean value11.7 in/yr) and lowest for basins underlain by the Central Basin aquifer (mean value 5.6 in/yr). Mean recharge values for the Cumberland Plateau, Highland Rim, and Valley and Ridge aquifers are 6.5, 7.4, and 6.6 in/yr, respectively. Gravity drainage characterizes ground-water flow in most surficial bedrock aquifer in Tennessee. Accordingly, a gravity yield analysis, which compares concurrent water-level and streamflow hydrographs, was used to estimate aquifer storage coefficient for nine study basins. The basin estimates range from 0.002 to 0.140; however, most estimates are within a narrow range of values, from 0.01 to 0.025. Accordingly, storage coefficient is estimated to be 0.01 for all aquifers in Middle and East Tennessee, with the exception of the aquifer in the inner part of the Central Basin, for which storage coefficient is estimated to be 0.002. Estimates of aquifer hydraulic diffusivity are derived from estimates of the streamflow recession index and drainage density for 75 drainage basins; values range from 3,300 to 130,000 ft^2/d (feet squared per day). Basin-specific and site-specific estimates of transmissivity are computed from estimates of hydraulic diffusivity and specific-capacity test data, respectively. Basin-specific, or areal, estimates of transmissivity range from 22 to 1,300 ft^2/d, with a mean of 240 ft^2/d In general, areal transmissivity is highest for basins underlain by the Cumberland Plateau aquifer (mean value 480 ft^2/d) and lowest for basins underlain by the Central Basin aquifer (mean value 79 ft^2/d). Mean transmissivity values for the Highland Rim, Valley and Ridge, and Blue Ridge aquifer are 320,140, and 120 ft^2/d respectively. Site-specific estimates of transmissivity, computed from specific-capacity data from 118 test wells in Middle and East Tennessee range from 2 to 93,000 ft^2/d with a mean of 2,600 ft^2/d Mean transmissivity values for the Cumberland Plateau, Highland Rim, Central Basin, Valley and Ridge, and Blue Ridge aquifers are 2,800,1,200, 7,800, 390, and 65Oft Id, respectively.

Tennessee↗

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

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↗

USGS2024

Turbidite correlation for paleoseismology

Marine turbidite paleoseismology relies on the assumption of synchronous triggering of turbidity currents by earthquake shaking to infer rupture extent and recurrence. Such inference commonly depends on age dating and correlation of the physical stratigraphy of deposits carried by turbidity currents (i.e., turbidites) across great distances. Along the Cascadia subduction zone, which lies offshore the Pacific Northwest, USA, turbidite facies in core photographs, X-ray computed tomography images, and magnetic susceptibility (MS) data exhibit differences in character over relatively short distances, which implies that not all deposits can be correlated with confidence. Thus, subjective correlation based on expected similarity over great distances and weak age constraints does not independently support paleoseismic models. We present a new method for correlating turbidites along the Cascadia margin that can yield a more objective and repeatable stratigraphic framework to underpin earthquake recurrence. We use dynamic time warping to correlate MS logs and measure correlation coefficients of core pairs to evaluate correlation strength. We then compare these measures to a distribution of correlation coefficients of randomly generated turbidite sequences and find that only a small number of core pairs can be correlated more confidently than randomly stacked turbidites. This methodology promises a more robust correlation strategy for future stratigraphic studies.

Oregon, Washington↗

USGS2026

Opportunities for the U.S. Geological Survey’s National Seismic Hazard Model to improve seismic risk assessment of critical infrastructure.

As fragility and risk modeling techniques and computational capabilities evolve, complemented by moving toward more routine and systematic seismic risk assessment of all buildings and critical infrastructure, the authors pose a few critical questions to investigate how the U.S. Geological Survey (USGS) National Seismic Hazard Models (NSHMs) can be used and enhanced further to serve such issues. In this paper, we use three examples from multiple sectors to (1) identify the role of USGS NSHMs in evaluating seismic risks to critical infrastructure, (2) quantify potential impacts from NSHM enhancements (i.e., [i] hazard curves for the vertical component of ground motion, [ii] stochastic event sets, and [iii] maps of probabilistic ground failure hazards), and (3) clarify the feasibility of relevant NSHM improvements. We illuminate that NSHMs are commonly used in location-specific performance assessments, whereas earthquake effects on critical infrastructure can be widespread across large geospatial regions. Further, we found that without the NSHM extensions considered here, risk can be severely underestimated, e.g., neglecting ground failure hazards can underestimate regional loss by a factor of two or more. Although many challenges remain, we developed example prototypes to clarify the feasibility of the NSHM extensions, which can facilitate improved management of risks to critical infrastructure.

Earthquake Spectra Journal↗