Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Earth's Future”

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,101 records · Page 10Linked to original sources
USGS2026

Regional models for postfire debris-flow likelihood and rainfall thresholds across the western United States

The U.S. Geological Survey (USGS) uses an empirical model developed with logistic regression (the ‘M1’ model) to rapidly assess debris-flow likelihood and to identify quantitative rainfall thresholds for debris flows after wildfire in the western United States. The M1 model was calibrated to a debris-flow inventory from southern California (United States) and has been applied throughout the western United States. Limited spatial coverage in the calibration dataset has motivated evaluation of M1 model accuracy outside the calibration region (e.g., the Sierra Nevada or the eastern Cascade Range, United States). Previous test cases showed that M1 overpredicts debris-flow likelihood and underpredicts rainfall thresholds for some locations (e.g., Arizona, northern California, Colorado, New Mexico, United States). We sought to improve the regional applicability of a debris-flow likelihood model by expanding the debris-flow inventory used for calibration, testing multiple potential models and generating an updated model framework. The updated inventory includes 3788 observations from 67 burned areas paired with short duration rainfall ratios. The updated model framework consists of a modified model structure and sets of coefficients calibrated separately to the entire updated inventory and to subsets of the inventory that intersect three Environmental Protection Agency (EPA) Level 2 ecoregions (Mediterranean California, Upper Gila Mountains and Western Cordillera). Comparisons of predictions from the updated models with observed rainfall and debris-flow activity show that the updated models outperform the M1 model by ~15%–60% and improve the uniformity of predictive performance across the western United States. The updated models also reduce false positive rates relative to M1 and generate rainfall thresholds that are better aligned with relative differences in regional climatology and debris-flow activity.

Arizona, California, Colorado, Idaho, Montana, Nev↗
USGS2021

Remote sensing-based actual evapotranspiration assessment in a data-scarce area of Brazil: A case study of the Urucuia Aquifer System

The large groundwater reserves of the Urucuia Aquifer System (UAS) enabled agricultural development and economic growth in the western Bahia State, in northeastern Brazil. Over the last several years, concern has grown around the aquifer’s diminishing water levels, and water balance (WB) studies are in demand. Considering the lack of measured actual evapotranspiration (ET a ), a major component of the water cycle, this work uses the Operational Simplified Surface Energy Balance (SSEBop) model to estimate ET a , and compares it to basin-scale estimates from the Soil Moisture Accounting Procedure (SMAP) monthly model and from an annual WB closure method, based on gridded meteorological data and the Gravity Recovery and Climate Experiment (GRACE) product. Additionally, a comparative assessment of different versions of the SSEBop parameterization was performed. Moderate Resolution Imaging Spectroradiometer (MODIS) imagery was used to implement eight different versions of the SSEBop algorithm over the UAS between 2000 and 2013. SSEBop and SMAP ET a yielded similar seasonal patterns, with correlation coefficient (r) up to 0.65, mean difference (MD) of 0.8 mm/month and mean absolute difference (MAD) of 18.5 mm/month. Comparison of SSEBop annual ET a estimates to annual SMAP and WB closure estimates yielded low MD (12.1 and −7.3 mm/year, respectively) and MAD (82.5 and 82.8 mm/year, respectively), but also low r values (0.00 and 0.37, respectively). The comparison of the different SSEBop versions indicated the need to incorporate a calibration step of the aerodynamic heat resistance (r ah ) parameter. SSEBop results were also used for land cover and drought monitoring. Analysis indicates that agriculture, associated with an increasing trend of atmospheric evaporative demand, is responsible for the decrease in groundwater levels and streamflow in the studied time period.

Urucuia Aquifer System↗
USGS2025

Simulation of groundwater flow in Wake County, North Carolina, 2000 through 2070

In 2019, the U.S. Geological Survey and Wake County Environmental Services began a collaborative study to evaluate groundwater resources and long-term groundwater availability in the county’s fractured-rock groundwater system. Wake County, in central North Carolina, is experiencing rapid population growth, associated land development, and changing water use. Hydrogeologic data including groundwater levels, aquifer testing, borehole fracture flow measurements, water-quality samples, and groundwater age-dating tracers were collected, along with findings from previous investigations, to help inform a conceptual model of the flow system used to develop a modular three-dimensional finite-difference groundwater-flow model (MODFLOW) for simulating historical and future groundwater conditions from 2000 to 2070. Hydraulic conductivity and transmissivity ranges were estimated from 17 slug tests and 21 borehole-flow measurements. Groundwater-quality analytical results from 19 sampling sites indicate that oxidation-reduction (redox) conditions varied within the regolith and bedrock and that minimal evaporation occurred before recharge entered the groundwater system. Age dating revealed mixtures of older and younger water, ranging from the 1940s to the 1990s—indicating variable flow pathways of recharge within permeable bedrock fracture zones. To simplify the complex fractured-rock groundwater system, two layers representing the regolith and the fractured bedrock were used in the MODFLOW model. Model calibration included parameter estimation and provided a reasonable fit to observed groundwater levels and estimated stream base flows. The model forecast scenarios incorporated future climate-model data for two emissions scenarios with land cover change projections to simulate potential impacts to future groundwater levels, recharge, and base flows. Recharge and base flow projections were largely within historical ranges, with no apparent long-term trends, but did indicate a slight downward shift in median values—likely, in part, because of differences in spatial resolution of input climate datasets. Seasonal patterns were consistent with historical data, with projections of possible increases in future winter recharge. Model limitations are discussed, and additional monitoring and model refinement needs are highlighted to support decision making for local groundwater management.

North Carolina↗
USGS2025

Assessment of western Oregon debris-flow hazards in burned and unburned environments

In the steep and mountainous environment of western Oregon, debris flows pose a considerable threat to property, infrastructure and life. Wildfire is commonly known to increase the susceptibility of steep slopes to debris flows, but the extent of this process in the western Cascades is not well understood. The US Geological Survey (USGS) currently estimates postfire debris-flow likelihood and triggering rainfall thresholds using a model calibrated to a southern California inventory of debris flows generated by excess runoff within the first year after fire. Because of a lack of available data, this model has not been tested in western Oregon, or in locations where postfire debris flows initiate via other mechanisms (e.g., shallow landslides or in-channel failures). Using repeat field observations and aerial imagery, we developed two new debris-flow inventories within and adjacent to the perimeters of five 2020 wildfires in western Oregon: Archie Creek, Holiday Farm, Beachie Creek, Lionshead and Riverside. The first inventory focuses on postfire debris flows (2020–2022); the second focuses on debris flows prior to fires (1995–2020). Our inventories of prefire and postfire debris flows were used to document initiation mechanisms in Oregon's western Cascades and to evaluate the effects of wildfire. We found that wildfire changed the distribution of debris-flow initiation mechanisms in the western Cascades. After the wildfires, annual rates of runoff-generated debris flows increased by 22% and the number of shallow landslide-initiated debris flows decreased by 17% relative to before the wildfires. Despite this shift, shallow landsliding was the dominant debris-flow initiation mechanism in both unburned and burned environments. We found the performance of the current USGS debris-flow likelihood model was degraded relative to other previously tested locations across the intermountain western United States. Our results highlight the need for improved postfire hazard assessment in western Oregon based on regional model calibration that is tuned to the dominant debris-flow initiation mechanisms.

Oregon↗
USGS2025

Comparison of creek and bay influences on salt marsh sediment budget and deposition patterns

The resilience of salt marshes with low organic production depends on their effective capture and retention of mineral sediment from adjacent waters. Little prior work has directly compared mechanisms of sediment import from wave-influenced marsh boundaries against those of tidal creeks. We used simultaneous deployment of net-deposition tiles and oceanographic sensors to identify the timing and magnitude of sediment import/export to, and redistribution within, a marsh in south San Francisco Bay. As the marsh has both an eroding bay-exposed scarp and a prominent tidal creek, we investigated the mechanisms and magnitudes of sediment import from the marsh-bay versus the marsh-creek interface. The strong daily sea breezes of the summer season produced most of the wave-driven erosion of the marsh scarp and controlled suspended sediment concentrations; the winter season had weaker winds punctuated by a few storms. A large seasonal difference in suspended sediment concentrations influenced both flood and ebb sediment fluxes to the marsh and led to much higher rates of import in the summer. Both bay-side and creek-side processes were important to total marsh sediment budget. Bay-side sediment contributions were more variable in time due to the bay-influenced environment, and creek-side contributions were overall larger, reflecting the large proportion of the marsh fed by creek water. Sediment was redistributed throughout the system, with erosion near the bay-edge, accretion near the creek-edge and slow import to the marsh interior. The marsh was net importing sediment in the summer and exporting in the winter from different rates of these processes; on an annual scale, the marsh was net importing despite rapid lateral marsh loss. These findings emphasize that a positive sediment budget does not imply a stable marsh and that both creek- and edge-side dynamics are important for marsh sedimentation and geomorphic trajectories. Further, we expand understandings of non-storm and seasonal controls on marsh sedimentation.

California↗
USGS2025

Wetter winters, drier summers: Quantifying the change in hydrological response around the Puget Sound area using the wflow_sbm hydrological model and CMIP6 projections

Climate change is expected to impact hydrological regimes worldwide, including the Pacific Northwest of the United States. This study investigates how climate change will affect river discharge in the Puget Sound region of the State of Washington, with a focus on King and Pierce Counties. We simulated river discharge under historical and future conditions using the physically based, spatially distributed wflow_sbm hydrological model, which was calibrated and validated against U.S. Geological Survey discharge records. Future forcing was based on an ensemble of six high-resolution CMIP6 climate models, which were bias corrected using empirical quantile mapping. The results indicate a decrease in summer discharges (5–10%) and an increase in winter discharges (5–10%) across the study region. The high discharges (90th percentile) are projected to increase in winter, and the low discharges are projected to decrease in summer, due to shifts in precipitation regimes, snowpack hydrology, and evapotranspiration. However, variability between individual CMIP6 models often exceeds the magnitude of ensemble mean changes, underscoring substantial uncertainty in climate projections and the importance of including multiple climate models in climate change analysis. Furthermore, model consensus increased with elevation, which could be the result of the higher elevation areas being driven by less diverse hydrological processes. These findings highlight potential challenges for regional water management, ecosystem health, and flood risk mitigation in the Puget Sound region under future climate conditions.

EarthArXiv↗
USGS2024

Influence of inherited structure on flexural extension in foreland basin systems: Evidence from the northern Arkoma basin and southern Ozark dome, USA

Extensional faults are key components of foreland basin systems. They form within the upper crust in response to flexure of the lithosphere and accommodate subsidence within the foredeep and forebulge depozones. Such faults are excellent proxies for orogenic system evolution and control the distribution of natural resources and hazards. However, the spatiotemporal evolution of flexural extension has not been documented previously at a regional scale, thereby limiting our understanding of underlying geodynamic controls. Here, we resolve late Paleozoic flexural extension in the northern Arkoma basin and southern Ozark dome, USA. We synthesize a large database of previous mapping, existing research, subsurface data, and geophysical data into 3D geologic and 2D kinematic models. Mesh surfaces representing several key horizons from the Carboniferous Period (ca. 335-306 Ma) were constructed. These surfaces were built from oil and gas well tops (n = ∼10,000) and surface geologic map contacts using an advanced kriging method. The mesh surfaces are offset by a complex 3D fault network, allowing detailed analysis of along-strike and down-dip variations in fault displacement. Analysis of the 3D model reveals a regular and repeated fault segmentation pattern wherein E -W striking, left- and foreland-stepping en échelon normal faults are segmented by inherited NE striking basement faults. Maximum vertical separation along the E -W normal faults is generally focused between the inherited NE-trending faults. This suggests that the inherited basement faults delocalized extensional strain during late Paleozoic normal faulting. Maximum vertical separation and fault localization may correlate to areas with high-amplitude positive magnetic anomalies interpreted as Mesoproterozoic granitic rocks. Speculative covariance of magnetic anomalies and fault displacements implies that the relatively strong basement granite concentrated stress, leading to localized faulting within the relatively thin sedimentary cover. Lastly, we show that flexural extension migrated southeast to northwest from the Chesterian-Morrowan (ca. 335-319 Ma) to the Desmoinesian (ca. 306 Ma). The migratory flexural extension may be explained by diachronous loading during Pangean assembly, or by synchronous loading but variable load compensation due to inherent factors.

Arkansas, Oklahoma↗
USGS2024

Isotopic evaluation of the National Water Model reveals missing agricultural irrigation contributions to streamflow across the western United States

The National Water Model (NWM) provides critical analyses and projections of streamflow that support water management decisions. However, the NWM performs poorly in lower-elevation rivers of the western United States (US). The accuracy of the NWM depends on the fidelity of the model inputs and the representation and calibration of model processes and water sources. To evaluate the NWM performance in the western US, we compared observations of river water isotope ratios ( 18 O / "> 16 O and 2 H / "> 1 H expressed in δ notation) to NWM-flux-estimated (model) river reach isotope ratios. The modeled estimates were calculated from long-term (2000–2019) mean summer (June, July, and August) NWM hydrologic fluxes and gridded isotope ratios using a mass balance approach. The observational dataset comprised 4503 in-stream water isotope observations in 877 reaches across 5 basins. A simple regression between observed and modeled isotope ratios explained 57.9 % ( δ 18 O) and 67.1 % ( δ 2 H) of variance, although observations were 0.5 ‰ ( δ 18 O) and 4.8 ‰ ( δ 2 H) higher, on average, than mass balance estimates. The unexplained variance suggest that the NWM does not include all relevant water fluxes to rivers. To infer possible missing water fluxes, we evaluated patterns in observation–model differences using δ 18 O diff ( δ 18 O obs − δ 18 O mod ) and d diff ( δ 2 H diff - 8 ⋅ δ 18 O diff "> ). We detected evidence of evaporation in observations but not model estimates (negative d diff and positive δ 18 O diff ) at lower-elevation, higher-stream-order, arid sites. The catchment actual-evaporation-to-precipitation ratio, the fraction of streamflow estimated to be derived from agricultural irrigation, and whether a site was reservoir-affected were all significant predictors of d diff in a linear mixed-effects model, with up to 15.2 % of variance explained by fixed effects. This finding is supported by seasonal patterns, groundwater levels, and isotope ratios, and it suggests the importance of including irrigation return flows to rivers, especially in lower-elevation, higher-stream-order, arid rivers of the western US.

Hydrology and Earth Systems Science↗
USGS2025

Estimation of dynamic geologic CO2 storage resources in the Illinois Basin, including effects of brine extraction, anisotropy, and hydrogeologic heterogeneity

Since the vast majority of carbon dioxide (CO 2 ) storage resources in the United States are in deep saline aquifers, optimizing the use of these saline storage resources could be crucial for efficient development of geologic CO 2 storage (GCS) resources and basin- or larger-scale deployment of GCS in the country. Maximum CO 2 injection rates can be enhanced by extracting brine from the CO 2 storage unit. However, disposal of the extracted brine is both a technological and economic challenge. The lowest-cost option would likely be reinjection of the extracted brine into another formation above or below the CO 2 storage unit. Therefore, it is important to estimate brine injectivity as it will constrain the potential to increase CO 2 injectivity at an injection site that has access to multiple geologic storage units where either CO 2 or brine can be injected. Using a simulation-optimization framework, coupled with a non-isothermal, multiphase CO 2 -water-salt equation-of-state module, we developed a computationally efficient method for evaluating optimization of simultaneous CO 2 injection, brine extraction, and brine (re)injection at hypothetical injection sites deployed across a geologic basin. The Illinois basin is ideal for testing our methodology because it contains multiple geologic storage units with seals in between them to isolate injection of CO 2 in one unit from interfering with the injection of either brine or CO 2 in another unit above or below it. In addition, we investigated the relative effects of variation in key geologic parameters as well as two reservoir structures (hydrogeologic heterogeneity/anisotropy and homogeneity/isotropy) on CO 2 injectivities and enhancement of CO 2 injectivity through extracting brine. Results suggest that permeability, depth, and especially thickness of the storage unit could be the most influential parameters determining CO 2 injectivity. They also suggest that only injecting CO 2 into the storage unit with the greatest injectivity, enhancing that unit’s injectivity by extracting brine, and disposing of the produced brine in other suitable units could maximize total CO 2 injectivity in limited regions of the basin. At the majority of simulated injection sites, however, we found that injecting CO 2 into all of the accessible and suitable storage units was more likely to maximize the CO 2 storage resource.

Illinois, Kentucky, Indiana↗
USGS2025

State shifts in the deep Critical Zone drive landscape evolution in volcanic terrains

Understanding the near-surface environment where atmospheric and solid earth processes interact, often termed the “Critical Zone,” is important for assessing resources and building resilient societies. Here, we examine a volcanic landscape in the Oregon Cascade Range, an understudied Critical Zone setting that is host to major regional water resources, pervasive silicate weathering, and significant geohazards. We leverage a bedrock age chronosequence to show that the volcanic Critical Zone undergoes a structural shift, from depth extents of >1 km to meters, over timescales of ~1 My. We map an active groundwater volume comparable to major continental lakes, stored at the Cascade Range crest. This state shift makes volcanic landscape evolution a unique probe of deep coupling between Earth systems.

PNAS↗
USGS2025

Projecting stream water quality using Weighted Regression on Time, Discharge, and Season (WRTDS): An example with drought conditions in the Delaware River Basin

Future water availability depends on understanding the responses of constituent concentrations to hydrologic change. Projecting future water quality remains a methodological challenge, particularly when using discrete observations with limited temporal resolution. This study introduces Weighted Regression on Time, Discharge, and Season for Projection (WRTDS-P), a novel, computationally efficient method that enables the projection of daily stream water quality under varying hydrologic conditions using commonly available discrete monitoring data. WRTDS-P model performance was validated using 39 sites in the Delaware River Basin (DRB) and four key constituents: specific conductance (SC), nitrate (NO 3 − ), magnesium (Mg 2+ ) and calcium (Ca 2+ ). Projections were tested against holdout data from the final 1 to 5 years of each time series, demonstrating robust predictive capability, with median Nash-Sutcliffe efficiencies of 0.67 for SC, 0.56 for NO 3 − , 0.65 for Ca 2+ , and 0.79 for Mg 2+ . Model uncertainty was correlated with indicators of hydrologic or geochemical mass-sinks, such as groundwater storage and adsorption in wetland soils. Drought scenario analyses for SC used ranges of reduced discharge including flows from the 1965 drought of record. Scenarios predicted widespread increases of SC, especially in southern DRB streams where baseline SC levels are already elevated. Fractional increases of SC were more uniformly distributed, indicating potential risk to sensitive ecosystems. Notably, drought-induced SC increases were positively correlated with interannual SC trends, indicating that hydrologic extremes could exacerbate ongoing salinization. This work provides a transferable and interpretable framework for projecting future water quality and assessing hydrologic risk to water resources and aquatic ecosystems.

Delaware, New Jersey, New York, Pennsylvania↗
USGS2024

Editorial: From cold seeps to hydrothermal vents: Geology, chemistry, microbiology, and ecology in marine and coastal environments

This Research Topic compiles contemporary studies on cold seeps, hydrothermal vents, mud volcanoes, and related seafloor features that are associated with focused fluid emissions and the transfer of carbon, other chemical species, and sometimes heat from the geosphere to the ocean. Because these features sometimes tap fluids and gas originating kilometers below the seafloor, they provide an important window into deep processes that are otherwise inaccessible to scientists. At the shallow portion of their journey, migrating fluids nearing the seafloor contribute to a range of unique biological, physical, and chemical processes within the sediments themselves and at the sediment-water interface. Seafloor fluid emissions play a critical role in global biogeochemical cycles, ocean chemistry, and possibly even climate change. Seafloor leakage points often emit hydrocarbon gases (especially methane and CO 2 ) and are sometimes the loci for deposition of seafloor minerals that have economic value. A burgeoning area of research focuses on natural products generated at these features, seeking compounds with potential pharmaceutical or other applications. Multidisciplinary studies have become routine for characterization of seafloor fluid emission sites, attesting to the inseparability of geologic, physical, chemical, and biological processes in these settings. It is increasingly common for researchers to combine in a single research cruise: subbottom imaging and seafloor mapping; porewater and water column geochemistry and gas sampling; sediment retrieval for lithologic, biostratigraphic, and solid phase analyses; and studies of benthic and subseafloor communities at the microbial to macrofaunal scales. This multidisciplinary approach has the advantage of ensuring the spatial and temporal coincidence of surveys and samples, an important factor at highly dynamic seafloor fluid emission sites. In addition, researchers often use remotely operated vehicles (ROVs), autonomous underwater vehicles (AUVs), or human-occupied vehicles (HOVs) to record video of the seafloor, compile photomosaics, collect targeted samples, and survey with high-resolution geophysical near-seafloor systems, providing a degree of detail about seafloor fluid emission sites that is unprecedented compared to most areas of the deep ocean. While rarer, long-term cabled observatories or shorter-term deployments of portable observatories are also used at some loci for seafloor fluid flux and are particularly helpful for capturing temporal variations at these dynamic features. Here we summarize the Research Topic’s contribution to multidisciplinary seafloor emission studies in the categories of cold seeps, mud volcanoes, and hydrothermal vents. Figure 1 shows the geographic distribution of the studies in this Research Topic and key features referred to in this Introduction.

Frontiers in Earth Science↗
USGS2024

Estimated groundwater recharge for mid-century and end-of-century climate projections, Kaua‘i, O‘ahu, Moloka‘i, Lāna‘i, Maui, and the Island of Hawai‘i

Demand for freshwater in the State of Hawaiʻi is expected to increase by roughly 13 percent from 2020 to 2035. Groundwater availability in Hawaiʻi is affected by a number of factors, including land cover, rainfall, runoff, evapotranspiration, and climate change. To evaluate the availability of fresh groundwater under projected future-climate conditions, estimates of groundwater recharge are needed. A water-budget model with a daily computation interval was used to estimate the spatial distribution of groundwater recharge for Kauaʻi, Oʻahu, Molokaʻi, Lānaʻi, Maui, and the Island of Hawaiʻi for recent climate conditions and three future-climate scenarios. Climate conditions from 1978 to 2007 were used as the reference period for recent climate conditions on each island. The three future-climate scenarios were developed using available high-resolution downscaled climate projections that include (1) a mid-century scenario using projected rainfall conditions for the Representative Concentration Pathway (RCP) scenario during 2041–71 with a total radiative forcing of 8.5 watts per square meter by the year 2100 (RCP8.5 2041–71 scenario), (2) a dry-climate scenario using projected rainfall conditions for the RCP8.5 scenario during 2071–99, and (3) a wet-climate scenario using projected rainfall conditions for the “Special Report on Emissions Scenarios” A1B scenario during 2080–99 for Maui, the RCP4.5 scenario during 2080–99 for Kauaʻi, Lānaʻi, and the Island of Hawaiʻi, and the RCP8.5 scenario during 2080–99 for Oʻahu and Molokaʻi. An additional drought scenario was added for Lānaʻi to assess the effect of extreme drought conditions during 2008–12 on groundwater recharge. All scenarios used 2020 land cover. Mean annual groundwater recharge is estimated to decrease between 5 and 55 percent on all six islands in this study for the mid-century and dry-climate scenarios relative to the reference-period recharge. Recharge is estimated to increase for Kauaʻi, Oʻahu, Molokaʻi, Lānaʻi, and Maui between 2 and 43 percent and decrease for the Island of Hawaiʻi by about 4 percent for the wet-climate scenario. Comparing the mid-century and dry-climate scenarios, all 110 aquifer systems (management areas defined by the State of Hawaiʻi Commission on Water Resource Management) from all six islands show similar direction in drying (104 aquifer systems) or wetting (6 aquifer systems) changes for recharge. However, among the three future scenarios, only 35 of 110 aquifer systems show similar direction in drying (30 aquifer systems) or wetting (5 aquifer systems) changes for recharge.

Hawaii↗
USGS2025

Constraining landslide frequency across the United States to inform county-level risk reduction

Informative landslide hazard estimates are needed to support landslide mitigation strategies to reduce landslide risk across the United States. Whereas existing national-scale landslide susceptibility products assess where landslides are likely to occur, they do not address how often , which is a critical element of landslide hazard and risk assessments. In particular, the U.S. Federal Emergency Management Agency's National Risk Index (NRI) requires landslide frequency estimates to inform expected annual loss estimates. We present county-level landslide frequency (landslides per area per year) estimates for the 50 US states. We applied Bayesian negative binomial regression to estimate both the expected (average) reported landslide frequency and full distribution of annual landslide counts for each county. We compared a suite of models that used combinations of landslide-susceptible area, probability of potentially triggering earthquakes, frequency of potentially triggering precipitation, and ecological region as predictors. We trained our models with landslide inventory data from counties with the most comprehensive records available nationwide and used zero-inflated negative binomial distributions as an incompleteness model to correct for temporal reporting gaps. We selected a preferred frequency model to inform the NRI based on information criteria and physically plausible parameter estimates. The model showed that average annual reported landslide frequencies vary by 5 orders of magnitude across US counties, ranging from 0.002 (0.00015–0.05) landslides 1000 km −2 yr −1 in Kusilvak Census Area, Alaska, to 29 (19–46) landslides 1000 km −2 yr −1 in Lake County, California, reflecting the country's strong variations in landslide susceptibility, earthquake probability, and other factors for which ecological region serves as a proxy. Counties with estimated frequencies in the top 20 % of all counties are predominately along the West Coast of the continental United States, in mountainous regions of the Pacific Northwest and Intermountain West, in locally steep or earthquake-prone regions of the Midwest and Southeast, along the Appalachians, in southern and southeastern Alaska, and on some Hawaiian islands. By examining the number of landslides predicted in 99th percentile years for each county, we identified that 26 % of US counties likely have potential for widespread landsliding with more than 10 landslides 1000 km −2 yr −1 , even when such large events have not been reported in the training data for that county. Overall, our results better represent the range of possible landslide frequencies and spatial variations than previous national-scale estimates reported in the NRI, and our approach can inform other risk-reduction and loss-mitigation efforts across the United States and globally.

Natural Hazards and Earth System Sciences↗
USGS2025

Strong shaking from past Cascadia Subduction Zone earthquakes encoded in coastal landforms

Strong earthquakes along subduction zones are often devastating events, but sparse records along some tectonic margins limit our understanding of seismic hazards. Constraining shaking intensities is critical, especially in subduction zones with infrequent but large-magnitude earthquakes like the Cascadia Subduction Zone (CSZ), where the lack of recorded ground motions has led to uncertainty in the severity and potential impacts of future earthquakes. Here we fill this observational gap with a novel inventory of quantitative estimates of past shaking intensities from geotechnical modeling of coastal landforms. One hundred fifty-four deep-seated landslides and 65 fragile geologic features constrain minimum and maximum peak ground accelerations, respectively. These estimates are broadly consistent with model predictions of M9 ruptures, suggesting strong shaking of 0.4–0.8 g during past CSZ earthquakes. Local discrepancies between our geologic shaking constraints and earthquake simulations may inform past rupture behavior, leading to better predictions of shaking intensity for future earthquakes.

California, Oregon, Washington↗
USGS2024

A benchmark dataset and workflow for landslide susceptibility zonation

Landslide susceptibility shows the spatial likelihood of landslide occurrence in a specific geographical area and is a relevant tool for mitigating the impact of landslides worldwide. As such, it is the subject of countless scientific studies. Many methods exist for generating a susceptibility map, mostly falling under the definition of statistical or machine learning. These models try to solve a classification problem: given a collection of spatial variables, and their combination associated with landslide presence or absence, a model should be trained, tested to reproduce the target outcome, and eventually applied to unseen data. Contrary to many fields of science that use machine learning for specific tasks, no reference data exist to assess the performance of a given method for landslide susceptibility. Here, we propose a benchmark dataset consisting of 7360 slope units encompassing an area of about 4,100 km 2 "> 4,100 km 2 in Central Italy. Using the dataset, we tried to answer two open questions in landslide research: (1) what effect does the human variability have in creating susceptibility models; (2) how can we develop a reproducible workflow for allowing meaningful model comparisons within the landslide susceptibility research community. With these questions in mind, we released a preliminary version of the dataset, along with a “call for collaboration,” aimed at collecting different calculations using the proposed data, and leaving the freedom of implementation to the respondents. Contributions were different in many respects, including classification methods, use of predictors, implementation of training/validation, and performance assessment. That feedback suggested refining the initial dataset, and constraining the implementation workflow. This resulted in a final benchmark dataset and landslide susceptibility maps obtained with many classification methods. Values of area under the receiver operating characteristic curve obtained with the final benchmark dataset were rather similar, as an effect of constraints on training, cross–validation, and use of data. Brier score results show larger variability, instead, ascribed to different model predictive abilities. Correlation plots show similarities between results of different methods applied by the same group, ascribed to a residual implementation dependence. We stress that the experiment did not intend to select the “best” method but only to establish a first benchmark dataset and workflow, that may be useful as a standard reference for calculations by other scholars. The experiment, to our knowledge, is the first of its kind for landslide susceptibility modeling. The data and workflow presented here comparatively assess the performance of independent methods for landslide susceptibility and we suggest the benchmark approach as a best practice for quantitative research in geosciences.

Earth-Science Reviews↗
USGS2026

Guidelines for producing integrated 210Pb and 14C age-models

Accurate reconstructions of past environmental changes are crucial in paleoecological research and require reliable chronologies of sedimentary archives. Establishing robust age-models and obtaining the most appropriate proxies for analysis is a complex scientific endeavor, requiring extensive resources and collaboration among specialists, including radiochronologists. Radiometric dating methods, such as 210 Pb and radiocarbon ( 14 C), are frequently employed to establish chronologies in aquatic sedimentary deposits and peat bogs. In this study, we review key aspects of sampling, analysis, and the principles underlying 210 Pb and 14 C age-models, focusing on methods for developing robust joint chronologies for paleoenvironmental research. Drawing largely from the authors' experiences and group discussions during and after a scientific workshop in 2022, we discuss important considerations for site selection, sampling strategies, and radiometric dating to construct integrated 210 Pb and 14 C age-models. Using expert consensus, this group – called Paleostats – aims to provide a set of best practices for other geochronologists with this methods paper. Among our conclusions, we emphasize the importance of accounting for site-specific factors such as prior information on sedimentation rates to establish appropriate sampling and analytical strategies. The use of appropriate coring devices can minimize disturbance to sediments and ensure the core surface remains intact and preserved until sectioning. Where excess 210 Pb is expected, sectioning at intervals of ≤1 cm provides an adequate sampling resolution for 210 Pb dating. Exceptions are possible, allowing for ~2–3 cm sections in areas with confirmed high sedimentation rates (e.g., > 1 cm yr −1 ). Recovering deeper core sections for 14 C dating with sufficient overlap allows for accounting errors in depth estimates made in the field. Special attention is advised during time intervals where validation proxies, such as the human-made radionuclides 137 Cs or post-bomb 14 C, are expected, and to determine the depth of secular equilibrium between 210 Pb and 226 Ra. Radiocarbon analyses are commonly performed by accelerator mass spectrometry, and age models are constructed mainly using Bayesian statistics with Markov Chain Monte Carlo techniques (e.g., Bacon ). A Bayesian approach ( Plum ) is now available for producing 210 Pb age-models, which infers the 210 Pb ex flux, eliminates the need for selecting an equilibrium depth, and allows dating cores with incomplete 210 Pb ex inventory. Plum offers improved chronologies by integrating raw 210 Pb and 14 C data, and these age-models can be enriched with other dating methodologies, such as identifying tephras and other well-recorded historical events. Harmonized reporting would contribute to making radiometric age-models reproducible, which would benefit from an international effort. Using 210 Pb and 14 C to produce integrated age-models may yield better insights into the interplay between natural and recent anthropogenic forcings on ecosystems. This can enhance our understanding of environmental processes and their impacts on climate change, ultimately supporting science-based assessments and decisions.

Earth-Science Reviews↗
USGS2025

Estimated average annualized tsunami losses for the United States

Tsunami hazards are substantial threats to coastal communities across the United States (U.S.) and its territories. U.S. states and territories collaborate through the National Tsunami Hazard Mitigation Program (NTHMP) to develop their own tsunami-hazard information for outreach and evacuation planning. An effort to curate this tsunami-hazard information to support comprehensive risk analysis at the national level has not yet been completed. In support of this effort, the Federal Emergency Management Agency (FEMA) collaborated with the NTHMP, the National Oceanic and Atmospheric Administration (NOAA) and the U.S. Geological Survey (USGS) starting in 2023. This collaboration included the collection and analysis of existing tsunami hazard data and methods in the U.S. Tsunami subject matter experts identified and selected scientifically defensible methods for estimating the risks to buildings and populations in coastal communities. These efforts may support decision making regarding resilience policies, priorities, strategies and funding levels. Tsunamis can be triggered by earthquakes, subaerial or submarine landslides, volcanic eruptions, glacial calving, near-earth objects, weather or other events. These events can cause severe destruction, injuries, and loss of life due to powerful currents and flooding. Tsunamis pose a substantial threat to the western United States and all U.S. territories, as described below. ■ Hawaii is threatened by distant tsunamis due to its central location in the Pacific Ocean basin and has a history of local events. ■ Alaska, particularly the Aleutian Islands, faces local tsunami threats due to proximity to the Alaska-Aleutian Subduction Zone, as well as distant tsunamis from around the Pacific Ocean basin. ■ The western coast of the U.S. is threatened by distant tsunamis from around the Pacific Ocean basin and local source tsunamis from earthquakes generated within the Cascadia Subduction Zone in the Pacific Northwest. ■ American Samoa faces local tsunami threats from earthquakes generated in the nearby Tonga Trench, as well as distant tsunami threats. ■ Guam and the Commonwealth of the Northern Mariana Islands are threatened by local tsunamis from the nearby Mariana Subduction Zone, as well as distant sources from around the Pacific Ocean Basin. ■ Puerto Rico and the United States Virgin Islands are threatened by multiple local and distant tsunami sources, such as the Puerto Rico Trench (PRT), given their location in the complex seismic region of the Caribbean Sea. Several historical events stand out because of their catastrophic impacts. ■ In the Pacific Northwest, the 1700 Cascadia earthquake caused a tsunami that affected coastal Native American communities, though the extent of the damage is not fully documented (Ludwin, et al., 2005). ■ In Puerto Rico, the 1918 earthquake triggered a tsunami that caused $77 million in damage in 2022 dollars and 116 fatalities, primarily along the western coast (Coffman et al., 1982). ■ The 1946 Aleutian Islands earthquake triggered a massive tsunami that devastated Hilo, Hawaii, killing 158 people and resulting in approximately $375 million in damage (adjusted to 2022 dollars) (Fisher et al., 2023). ■ The 1964 Alaska earthquake (M 9.2) generated tsunamis that caused severe destruction in some communities across Alaska, Oregon, and California. This disaster led to a total of 124 fatalities and approximately $2.9 billion in property damage (adjusted to 2022 dollars) (Brocher et al., 2014) (Alaska Science Center, 2024). ■ In American Samoa, a tsunami generated by the 2009 Samoa earthquake (Mw 8.1) caused widespread devastation, resulting in 34 confirmed fatalities (Apatu et al., 2013) and economic losses exceeding $160 million (adjusted to 2022 dollars) (DHS, 2011). More recent events, including the 2010 Chile earthquake, the 2011 Japan earthquake, and the 2022 Tonga volcanic eruption, resulted in millions of dollars in damage to numerous ports and harbors in the U.S. South Pacific territories, Hawaii, and along the west coast of the U.S. (Lynett, et al., 2022) (Wilson, et al., 2013). Since these events, the expansion of the built environment in lowlying areas along the coast has increased the exposure of buildings and people, thereby further escalating community risk from tsunamis. This report provides a comprehensive national assessment of earthquake-generated tsunami risk. It does not include impacts from tsunamis generated by landslides, volcanic eruptions, glacial calving, near-earth objects, weather, or other events. This study is based on the best available hazard data from the U.S. Pacific Coast (California, Oregon and Washington), Alaska, Hawaii, U.S. Pacific Territories (American Samoa, Guam and Commonwealth of the Northern Mariana Islands) and Caribbean Territories (Puerto Rico and United States Virgin Islands). Tsunami risks associated with states along the East Coast, Gulf Coast, and Great Lakes are not included in this study because Hazus 6.1 software (FEMA 2024a) does not currently include the ability to analyze tsunami risk in those states. Once modeling capabilities and tsunami hazard data become available for additional states, FEMA may incorporate these data into future editions of this study.

Alaska, California, Hawaii Oregon, Washington↗