USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of the American Water Resources Association”

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,691 records · Page 76Linked to original sources
USGS2025

Arctic fold-and-thrust belts

The modern Arctic has been formed through a series of continent–continent collisions, accretion of terranes and phases of crustal extension. The Neoproterozoic Timanian, Paleozoic Caledonian and Uralian, and late Mesozoic Verkhoyansk–Kolyma, Chukotkan and Brookian orogenies formed several large fold-and-thrust belts (FTBs). The FTBs are exposed across vast areas of continents and continue offshore to form a complex tectonic basement for thick sedimentary basins, playing an important role in the history of accumulation and deformation of younger unmetamorphosed sedimentary successions that are the subject of this volume. Recognition of the importance of FTBs in the Arctic geological history and their role as a controlling factor of development of Arctic sedimentary basins resulted in this chapter, in which we review the current state of knowledge about Arctic FTBs and highlight questions that remain to be addressed. Enclosure D , a map showing boundaries of the FTB and their internal first-order structural fabric, is a part of the overview.

Book chapter↗
USGS2025

Earthquake magnitude and source parameter estimation with a distributed acoustic sensing dataset in the Gorda subduction zone

Distributed acoustic sensing (DAS) systems offer a cost‐effective way to create large‐scale strainmeter arrays for seismological applications using fiber‐optic cables. DAS‐based strain measurements are known to be influenced by various factors, bringing into question their general reliability for accurate earthquake characterization. A 15‐km‐long DAS deployment in northern California was operational within 3 days of the 2022 M w 6.4 Ferndale earthquake and ran continuously throughout the aftershock sequence. We utilize these aftershock data to validate DAS‐based strain measurements in two ways. We first test the accuracy of DAS‐based magnitude estimates from peak dynamic strains by comparing them with magnitude and attenuation scaling relations derived independently from traditional borehole strainmeter (BSM) data. We demonstrate that DAS‐based magnitudes are comparable to BSM‐based magnitudes when corrections for variations in site response along the fiber‐optic cable are properly made. Magnitude errors are spatially correlated, potentially because of factors such as finite‐fault effects (e.g., stress drop) or more complex, unmodeled path attenuation or because of wave propagation effects in heterogeneous media. We then apply more advanced source characterization methodology to the DAS data using a time‐domain empirical Green’s function (EGF) deconvolution approach to measure details of the moment rate history. The EGF approach using DAS data depends on careful treatment of distorting factors such as anthropogenic sources of noise and optical phase wrapping but successfully isolates source spectra for moderate‐magnitude earthquakes: source spectral ratios obtained from DAS data, broadband seismometer data, and BSM data in the same region show consistent results, revealing differences in directivity and spectral shape among earthquakes. Although further research is needed to refine source‐time‐function estimation techniques for DAS data, particularly for larger magnitude events, these case studies demonstrate the clear potential of DAS for earthquake source characterization.

California↗
USGS2025

Quantitative evaluations of earthquake early warning performance using “Did You Feel It?” and post-alert surveys

We examine responses to the U.S. Geological Survey’s “Did You Feel It?” (DYFI) survey and its companion earthquake early warning (EEW) questionnaire to assess the performance of the U.S. ShakeAlert EEW system directly from the alert recipients’ perspectives. ShakeAlert rapidly detects earthquakes and develops alert information, but as official alert delivery partners issue these alerts, it is thus difficult to determine how many people were alerted and when. We investigate DYFI reports for six California earthquakes that had EEW alerts and substantial responses to the DYFI EEW questionnaire. Comparisons of ShakeAlert predictions to reported intensities demonstrate that magnitude estimation accuracy is not necessarily indicative of ground-motion prediction accuracy. Perceived warning time distributions indicate that estimating maximum-expected warning times using the S-wave arrival is a reasonable assumption when discussing public EEW performance. However, we also find many reports of shorter warning times, late alerts, and missed alerts than expected based on ShakeAlert publication times, suggesting alert delivery latencies are substantial and highly variable. The novelty of our analysis is that we demonstrate that the DYFI EEW survey provides useful EEW efficacy information—independent of the specific alerting pathway—that can be used to inform our choices for conveying EEW performance.

California↗
USGS2025

Case study of deep learning image segmentation for the purposes of rapid 2D petrographic analysis in volcanic rocks

Automation using deep learning methods is a useful alternative to manual methods of petrographic segmentation, but often requires user familiarity with coding and/or algorithms. We examine the Dragonfly TM program's deep learning tools for application by users with a variety of skill levels as a method for petrographic image segmentation. An image processing methodology, bimodal image stacking, was created for low-input-data, high-efficacy training of models which can then be applied to varied samples. Using backscatter electron images we show that the resulting model segmentations agree with manual segmentation total and modal crystallinity values within 5%, and calculated plagioclase crystal size distribution (CSD) values within 2σ, despite limitations in discriminating mafic phases. Model creation and training takes <24 hours, 1–3 hours of which are supervised, and the resultant model can then be applied to new uncharacterized samples in <15 minutes per image. This allows for non-experts to create and utilize deep learning models to segment images of variable brightness and texture, at low user-time cost and resulting in size and shape data which are within uncertainty of manual segmentation. While some limitations are noted (for example, sieve-textured phases may need manual correction, and different minerals with similar BSE intensity may not be resolved as separate phases), this methodology can be utilized for general application of models to wide ranges of volcanic crystalline and bubble textures, and to create a library of models for rapid petrological analysis during volcanic eruptions.

Alaska↗
USGS2026

2025 USGS National Seismic Hazard Model for Puerto Rico and the U.S. Virgin Islands: Overview of model and hazard results

The U.S. Geological Survey recently updated the National Seismic Hazard Model (NSHM) for Puerto Rico and the U.S. Virgin Islands (PRVI). The first version of the PRVI NSHM was released in 2003, and therefore this 2025 update includes over 20 years of new geologic, geophysical, and engineering data, methods, and models. Updates follow similar efforts performed in the recent 2023 50-state NSHM. However, this is the first NSHM in which we: (1) apply an inversion methodology to subduction interface fault sources in the earthquake rupture forecast (ERF) model; (2) develop scaled backbone median ground-motion models and independent aleatory variability models that are applied in the ground-motion characterization (GMC) model; and (3) calculate epistemic uncertainty related to alternative scenarios in the ERF and GMC models for all grid points in the study region. Long-term time-independent mean hazard calculations were performed for peak ground acceleration and 5%-damped pseudospectral acceleration at 21 spectral periods from 0.01- to 10.0-s, for eight National Earthquake Hazards Reduction Program site conditions ranging from V S30 = 150 to 1500 m/s, and for 2%, 5%, and 10% in 50-year probabilities of exceedance (return periods of 2475, 975, and 475 years, respectively). Epistemic uncertainty, in the form of selected percentiles, is also provided for a suite of test sites and all grid points in the study region for limited periods, site conditions, and probabilities of exceedance. Selected results, including comparisons with the 2003 PRVI NSHM, are shown and discussed for selected periods, site conditions, and probabilities of exceedance. When comparing the 2025 PRVI NSHM with the 2003 PRVI NSHM, hazard is generally higher at shorter periods and lower at longer periods, as a result of updates in both ERF and GMC models. The 2025 PRVI NSHM is applicable for return periods greater than ∼475 or less than ∼10,000 years.

Puerto Rico, U.S. Virgin Islands↗
USGS2025

Probing lower-crustal fault properties with frequency-dependent tidal tremor triggering

The way seismicity responds to periodic stress perturbations offers crucial insights into the processes that can trigger an earthquake. Laboratory and theoretical analyses have shown that the period of imposed forcing and source properties affect the sensitivity to triggering, but frequency-dependent triggering of tectonic faults is poorly understood. The rate of low-frequency earthquakes (LFEs) near Parkfield, California has been found to be strongly correlated with solid earth tides. Tidal forcing acts over multiple frequencies, and the sensitivity to tidal triggering of LFEs therefore provides a unique opportunity to probe the physics of earthquake triggering and underlying fault properties. Here, we find that the response of LFEs to solid earth tides at diurnal and semi-diurnal frequencies is highly variable but spatially coherent along the San Andreas Fault. Using rate-state friction modeling, we find that the variation of the amplitude of tidal modulation is mainly affected by the spatial variation of the background effective normal stress, whereas the spatially varied tidal modulation at diurnal and semi-diurnal frequencies is mainly affected by the fault frictional property. The spatial variations of the LFE response reveal a heterogeneous lower-crustal geologic structure and complex physical faulting processes below the rupture zone of eventual great earthquakes.

California↗
USGS2024

Field evidence and indicators of rockfall fragmentation and implications for mobility

Rockfall fragmentation can play an important role in hazard studies and the design of protective measures. However, the current lack of modeling tools that incorporate rock fragmentation mechanics is a limitation to enhancing studies and design. This research investigates the fragmentation patterns of rockfalls and analyzes the resulting distribution of fragment sizes within corresponding rockfall deposits. We focus on small rock fragments, which provide insights into the dynamics of the rockfall event and can be used as input for numerical modeling. We analyzed multiple rockfall events from locations worldwide, each exhibiting different degrees of fragmentation. Using image analysis techniques, we mapped all visible blocks, determined their volumes, and measured the distances they travelled from the initial point of impact. A key finding is the identification of three indicators of fragmentation. First, in cases where fragmentation was largely absent, we observed a trend of increasing block size with distance from the impact point or source area, which aligns with previously published findings. However, for energetic rockfall events characterized by intense fragmentation, we observed that small fragments exhibited longer travel distances compared to larger fragments. This distinction allowed us to differentiate blocks primarily resulting from the disaggregation process from those primarily resulting from dynamic fragmentation, with implications for rockfall mobility. Second, although the size distribution of rockfall deposits exhibits a power-law scaling for volumes larger than a minimum size threshold corresponding to a rollover of the distribution, in some case studies a deviation from power-law scaling is observed, indicating a process of larger block comminution due to fragmentation. Third, we found that rockfalls with fragmentation experience reduced mobility, indicated by higher reach angles, and higher lateral dispersion showing a wider distribution of trajectories. We interpret these findings as being directly related to the energy-consuming nature of fragmentation, which prevents farther deposition of fragmented rock blocks.

Albacete province, Lombardy and Aosta Valley, Yose↗
USGS2025

The influence of postfire root strength decay on shallow landslide susceptibility in western Oregon

Severe wildfires can increase the potential for postfire landslides, partly due to the loss of vegetation and root reinforcement. This study investigates the multi-year dynamics of root strength following wildfire in coniferous forests of western Oregon and evaluates implications for postfire slope stability. Root tensile strength was measured through laboratory testing of 969 roots collected from 26 test pits over four postfire timeframes: unburned conditions, 4–12 months postfire (Cedar Creek, 2022), 24–35 months postfire (Holiday Farm, 2020), and 59–67 months postfire (Eagle Creek, 2017). Variables analyzed include root diameter, tensile thread strength, stiffness, progressive tensile strength with displacement, time-dependent strength and ductility, root area ratio, and cohesion. Results indicate a 50% reduction in root strength after wildfire, with the lowest values occurring around 4 years postfire. Projections indicate it could take 10–22 years following fire for root strength to recover to 70% and 90% of the prefire original root strength, respectively. This timeline indicates a multi-year window of vulnerability for shallow landslides that reaches a maximum approximately 4 years following fire. Brittle failure behavior was observed in burned root systems compared to unburned root systems, signaling a more abrupt loss of strength at yield. Simple slope stability analyses show that reduced root strength can lead to instability on more gentle slopes and under less saturation in comparison to unburned conditions. These findings highlight the critical role of root reinforcement in postfire slope stability and the long-term implications of wildfire disturbance on landslide susceptibility.

Oregon↗
USGS2026

Hazard potential of compound flooding from rainfall, storm surge, and groundwater in coastal New York and Connecticut

Compound flood events, the co-occurrence of multiple flood drivers, can result in flood hazard potential exceeding that of any single driver alone. To evaluate compound flooding in a semi-urbanized coastal area, historical records dating back to 1970 are used to study the co-occurrences of high precipitation, storm surge, and shallow groundwater conditions along the coastlines of New York and Connecticut. Joint return periods for coincident precipitation-surge events were computed using statistical dependence models and compared to the assumption of independence as a ratio, referred to here as a return period adjustment. Results indicate distinct seasonality where compound events in the area disproportionately occur in the cold season between October and April. Return period adjustments range from a factor of 1 to almost 9, demonstrating the range in precipitation-storm surge dependence across the study area. Across all 24 station triad locations, groundwater levels were elevated during times of precipitation- surge co-occurrence, reflecting the tendency for coastal storms and shallow groundwater conditions to co-occur seasonally. The result is a pseudo-trivariate compound flood hazard score and corresponding hazard map that integrates dependence between daily precipitation-surge events and overall monthly groundwater levels (as a precondition) into a relative compound hazard score. The location with the highest compound flood hazard score is on the south shore of Long Island, as well as locations across coastal Connecticut where groundwater levels compound the co-occurrence of heavy precipitation and storm surge.

Connecticut, New York↗
USGS2025

Kiloyear cycles of carbonate and Mg-silicate replacement at Von Damm hydrothermal vent field

The Von Damm vent field (VDVF) on the Mid-Cayman Rise in the Caribbean Sea is unique among modern hydrothermal systems in that the chimneys and mounds are almost entirely composed of talc. We analyzed samples collected in 2020 and report that in addition to disordered talc of variable crystallinity, carbonates are a major class of mineral at VDVF. The carbonate minerals include aragonite, calcite, magnesium-rich calcite, and dolomite. Talc and carbonate mineral textures indicate that, rather than replacing volcanic host rock, they precipitate from the mixing of hydrothermal fluids and seawater at the seafloor, occurring in chimneys and surrounding rubble. Alternating precipitation of this mineral assemblage is pervasive, with carbonate minerals typically being succeeded by talc, and with indications that in some cases talc and carbonate minerals replace one another. Stable carbon isotopic data indicate the carbonate minerals originate from the mixing of seawater and hydrothermal fluid, which is supported by U-Th data. Radiocarbon calcite ages and talc 234 U- 230 Th isochron ages indicate mineral ages spanning over thousands to tens of thousands of years. Analyses of these samples illustrate a dynamic system that transitions from carbonate-dominated to Mg-silicate−dominated precipitation over time scales of thousands of years. Our observations raise questions regarding the eventual fate of seafloor precipitates and whether carbonate and silicate minerals in such settings are sequestered and represented in the rock record.

Geology↗
USGS2025

A glimpse into the future of tectonic tremor monitoring

Tectonic tremor is a weak, long-duration seismic signal often observed in subduction zones and on some other plate-bounding faults. Because of tremor's characteristically low amplitude (and low signal-to-noise) and lack of clear phase arrivals, detecting and locating tremor usually requires techniques distinct from those applied to typical earthquakes. Major advances in detection and understanding of tremor have derived in the past from a powerful combination of new data and new analysis techniques. In a recent study, Sagae et al. (2025, https://doi.org/10.1029/2025jb031348 ) exploit that combination again, developing a new machine-learning based workflow and applying it to the S-net cabled seismic network in the Japan trench offshore northern Honshu. Their approach, although complex, succeeds in detecting several times more tremor activity than earlier studies, resulting in new insights and providing a blueprint for similar approaches that could be applied elsewhere. As real-time earthquake monitoring adopts similar tools, it may present an opportunity to bring tremor monitoring into operational workflows. In turn, this could solidify tremor monitoring as a component of future operational earthquake forecasting.

JGR Solid Earth↗
USGS2024

Bathymetric and velocimetric surveys at highway bridges crossing the Missouri River between Kansas City and St. Louis, Missouri, May 19–26, 2021

Bathymetric and velocimetric data were collected by the U.S. Geological Survey, in cooperation with the Missouri Department of Transportation, near nine bridges at eight highway crossings of the Missouri River between Kansas City and St. Louis, Missouri, from May 19 to 26, 2021. A multibeam echosounder mapping system was used to obtain channel-bed elevations for river reaches about 1,640 to 1,840 feet (ft) longitudinally and generally extending laterally across the active channel from bank to bank during low to moderate flood-flow conditions. These surveys provided channel geometry and hydraulic conditions at the time of the surveys and provided characteristics of scour holes that may be useful in developing or verifying predictive guidelines or equations for computing potential scour depth. These data also may be useful to the Missouri Department of Transportation as a low to moderate flood-flow assessment of the bridges for stability and integrity issues with respect to bridge scour during floods. Bathymetric data were collected around every in-channel pier. Scour holes were present at most piers for which bathymetry could be obtained, except those on banks or surrounded by riprap. Occasionally, scour holes were minor and difficult to discern from nearby dunes and ripples. All the bridge sites in this study were previously surveyed and documented in previous studies. Comparisons between bathymetric surfaces from the previous surveys and those of the current (2021) study do not indicate any consistent correlation between channel-bed elevations and streamflow conditions. The average difference between the bathymetric surfaces varied from 1.59 ft higher to 0.95 ft lower in 2021 than 2017, which corresponds to a gain of 100,200 cubic yards and a loss of 55,800 cubic yards, respectively. The average difference between the bathymetric surfaces varied from 2.74 ft higher to 3.05 ft lower in 2021 than 2013, which corresponds to a gain of 111,500 cubic yards and a loss of 169,200 cubic yards, respectively. The average difference between the bathymetric surfaces varied from 4.52 ft higher to 1.38 ft lower in 2021 than 2011, which corresponds to a gain of 221,100 cubic yards and a loss of 90,300 cubic yards, respectively. The most substantial overall net gain was 221,100 cubic yards between 2011 and 2021 at structures L0550 and A4497 at Jefferson City (site 20). The large net gain likely results from a combination of the mitigation of the scour holes near pier 4 of both bridges and the substantially lower flow in 2021 than in 2011. Alternatively, the most substantial overall net loss was 169,200 cubic yards between 2013 and 2021 at structure A6288 at Hermann (site 21), despite comparable streamflows. Pier size, nose shape, and skew to approach flow had a substantial effect on the size of the scour hole observed at a given pier. Larger and deeper scour holes were present at piers with wide or blunt noses caused by exposed footings or caissons. When a pier was skewed to primary approach flow, the scour hole was generally deeper and larger than at a similar pier without skew; furthermore, the shape of the scour hole near skewed piers in this study generally was longer and deeper on the side with impinging flow. At structure A6288 at Hermann (site 21), the scour hole near pier 5 was difficult to discern from nearby dunes and ripples, whereas the upstream edge of the footing was visible at pier 4, which likely contributes to the larger scour hole near that pier; the top of the footing may blunt the horseshoe vortex at pier 5, but the exposed front of the footing may exacerbate the vortex at pier 4.

Missouri↗
USGS2025

Groundwater age estimates for the Mississippi River Valley alluvial aquifer based on tracer data collected during 2018–20

This study characterized groundwater age across the Mississippi River Valley alluvial aquifer (MRVA). Groundwater samples from 69 MRVA wells and 19 wells in Tertiary units of the Mississippi embayment aquifer system (MEAS) were analyzed for sulfur hexafluoride (SF 6 ), tritium ( 3 H), helium (He), and (or) carbon-14 of dissolved inorganic carbon ( 14 C). The age distributions of 89 samples were estimated by fitting lumped parameter models to processed tracer concentrations with the U.S. Geological Survey software TracerLPM. Mean ages of MRVA groundwater samples ranged from 12 to 22,000 years, with a median of 140 years. Mean ages of MEAS groundwater samples ranged from 230 to 52,000 or more years, with a median of 13,500 years. The spatial distribution of MRVA groundwater ages was found to be influenced by depth, inflow of groundwater from deeper units, and soil saturated hydraulic conductivity. In parts of the MRVA, the spatial distribution of MRVA groundwater ages was found to be influenced by annual recharge and (or) annual groundwater pumpage.

Alabama, Arkansas, Illinois, Kentucky, Louisiana, ↗
USGS2024

A global assessment of SAOCOM-1 L-band stripmap data for InSAR characterization of volcanic, tectonic, cryospheric, and anthropogenic deformation

SAOCOM-1 is an L-band (23.5 cm) synthetic aperture radar (SAR) constellation made up of two satellites launched in 2018 and 2020 by Comisión Nacional de Actividades Espaciales (CONAE, Argentina). In this contribution, we present a global summary of interferometric SAR (InSAR) observations of ground deformation with SAOCOM-1 stripmap data for tracking volcanic, tectonic, glacier, and anthropogenic deformation. These examples include: 1) episodes of unrest at volcanoes in the Aleutian Islands, Southern Andes, and Italy, with line-of-sight (LOS) deformation from 4 cm/yr in InSAR time series to ~70 cm in interferograms; 2) dike intrusions in Hawai’i; 3) earthquakes in the Andean fold and thrust belt and the East Anatolian fault; 4) ice flow of the Southern Patagonia icefield; and 5) subsidence due to lithium brine extraction in the Salar de Atacama basin (northern Chile). Comparisons between SAOCOM-1, ALOS-2 SM3, Sentinel-1, and TerraSAR-X/ TanDEM-X/PAZ (TSX/TDX/PAZ) mean velocities from InSAR time series show a 1:1 ± 3% correlation in the LOS velocity, which highlights the high accuracy of SAOCOM-1 data. The minimum deformation that we measured in individual interferograms is 4 ± 0.6 cm. One limitation of SAOCOM-1 is the lack of a global acquisition program, which reduces its global and broader applications. Considering the repeat periods, background observation program, and lack of a controlled orbital tube, the best suited targets for SAOCOM-1 InSAR are two. First, volcanoes that deform with secular rates located in vegetated regions in mid- and high-latitudes, and/or that undergo transient episodes of fast deformation in which C-band coherence is lost quickly. Second, glaciers where coherence can be sustained during the repeat period of eight days.

IEEE Transactions on Geoscience and Remote Sensing↗
USGS2023

Evaluation of models for estimating hydraulic conductivity in glacial aquifers from NMR logging

Nuclear magnetic resonance (NMR) logging is a promising method for estimating hydraulic conductivity ( K ). During the past ∼60 years, NMR logging has been used for petroleum applications, and different models have been developed for deriving estimates of permeability. These models involve calibration parameters whose values were determined through decades of research on sandstones and carbonates. We assessed the use of five models to derive estimates of K in glacial aquifers from NMR logging data acquired in two wells at each of two field sites in central Wisconsin, USA. Measurements of K , obtained with a direct push permeameter (DPP), K DPP , were used to obtain the calibration parameters in the Schlumberger-Doll Research, Seevers, Timur-Coates, Kozeny-Godefroy, and sum-of-echoes (SOE) models so as to predict K from the NMR data; and were also used to assess the ability of the models to predict K DPP . We obtained four well-scale calibration parameter values for each model using the NMR and DPP measurements in each well; and one study-scale parameter value for each model by using all data. The SOE model achieved an agreement with K DPP that matched or exceeded that of the other models. The Timur-Coates estimates of K were found to be substantially different from K DPP . Although the well-scale parameter values for the Schlumberger-Doll, Seevers, and SOE models were found to vary by less than a factor of 2, more research is needed to confirm their general applicability so that site-specific calibration is not required to obtain accurate estimates of K from NMR logging data.

Wisconsin↗
USGS2024

Projections of multiple climate-related coastal hazards for the US Southeast Atlantic

Faced with accelerating sea level rise and changing ocean storm conditions, coastal communities require comprehensive assessments of climate-driven hazard impacts to inform adaptation measures. Previous studies have focused on flooding but rarely on other climate-related coastal hazards, such as subsidence, beach erosion and groundwater. Here, we project societal exposure to multiple hazards along the Southeast Atlantic coast of the United States. Assuming 1 m of sea level rise, more than 70% of the coastal residents and US \ $ 1 trillion in property are in areas projected to experience shallow and emerging groundwater, 15 times higher than daily flooding. Storms increase flooding exposure by an order of magnitude over daily flooding, which could impact up to ~50% of all coastal residents and US \ $ 770 billion in property value. The loss of up to ~80% of present-day beaches and high subsidence rates that currently affect over 1 million residents will exacerbate flooding and groundwater hazard risks.

Florida, Georgia, North Carolina, South Carolina, ↗
USGS2025

Site-specific, extended ShakeMaps for earthquake engineering applications

The U.S. Geological Survey (USGS) routinely produces ShakeMaps of shaking intensity across the globe. Due to practical constraints, the number of response spectral periods was limited to three standard periods (0.3, 1.0, and 3.0 sec). We have recently developed the tools that are necessary to expand this functionality to include 22 periods (matching the current U.S. National Seismic Hazard Model periods) as well as the orientation-independent components (e.g., “RotD50”). We refer to ShakeMap products that include these extensions as “extended ShakeMaps.” The added level of complexity motivated us to also develop a user-friendly tool called the “ShakeMap Sampling Tool” (SST) that gives all the estimated shaking metrics for a specific location (or list of locations). Additionally, we develop a web application where users can input locations of interest and view/download the SST results. We further familiarize users with the concept of “Composite ShakeMaps.” For earthquakes sequences such as a mainshock and larger foreshocks and aftershocks, this provides a map of the maximum value of each shaking metric, which is useful for overall loss estimates, the full extent of ground failure triggering potential, and a better portrayal of the repeated shaking levels at a given point for a series of earthquakes. Such a site-specific shaking history facilitates earthquake forensics at building or infrastructure sites for which damage may be of concern, as described in the Disproportionate Damage Earthquake trigger specified in the IEBC (2018, Section 405.2.2) and in developing ATC-145 guidelines (Guidelines for Post-Earthquake Assessment, Repair, and Retrofit of Buildings). The composite ShakeMap can be combined with the SST for a variety of earthquake-hazard applications, such as systematically inferring triggering shaking estimates at specific sites of geotechnical interest for landsliding, liquefaction, and lateral-spreading hazards.

Kahramanmaraş↗