Search USGS⌕ Search

SEARCH · Search USGS

Results for “Applied Computing and Geosciences.”

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

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

1,656 records · Page 3Linked to original sources
USGS2026

SlideDetect: Spatio-temporal landslide detection using a three-dimensional convolutional neural network

Landslides pose a serious and ongoing threat to both human lives and infrastructure worldwide; therefore, it is of interest to predict where and when landslides are likely to occur. Advances in machine learning techniques have spurred numerous studies aimed at estimating relative landslide propensity, but are limited to spatial (as opposed to temporal) prediction due to the sparsity of landslide timing data. We address this data gap by training SlideDetect, a 3-dimensional convolutional neural network (3D CNN), to identify landslides based on their spatial and temporal occurrence within multitemporal image stacks. We use an inventory of landsides triggered by the 2018 Hokkaido earthquake and two years of monthly composite optical imagery spanning this event. The model can identify not only landslide location but also landslide date with an area under the precision-recall curve (PR-AUC) of 0.84. We further present a new standard for presenting PR curve results that explicitly compares model performance at different confidence thresholds, allowing for clearer model evaluation and comparison. Our new approach to constraining landslide timing paired with this more consistent and objective method for evaluating model performance shows considerable promise, and with further application and testing, SlideDetect could enhance the data availability and tools needed to advance landslide hazard and risk assessments.

JGR Machine Learning and Computation↗
USGS2025

Identification of representative earthquakes for probabilistic tsunami hazard analysis (PTHA) using earthquake rupture forecasts and machine learning

As probabilistic tsunami hazard analysis (PTHA) focuses more on assessments for localized, populous regions, techniques are needed to identify a subsample of representative earthquake ruptures to make the computational requirements for producing high-resolution hazard maps tractable. Moreover, the greatest epistemic uncertainty in seismic PTHA is related to source characterization, which is often poorly defined and subjective. We address these two salient issues by applying streamlined earthquake rupture forecasts (ERFs), based on combinatorial optimization methods, to an unsupervised machine learning workflow for identifying representative ruptures. ERFs determine the optimal distribution of a millennia-scale sample of earthquakes by inverting the observed slip rate on major faults. We use two previously developed combinatorial optimization ERFs, integer programming and greedy sequential, to produce the optimal location of ruptures with seismic moments sampled from a regional Gutenberg–Richter magnitude–frequency distribution. These ruptures in turn are used to calculate peak nearshore tsunami amplitude, using computationally efficient tsunami Green's functions. An unsupervised machine learning workflow is then used to identify a small subsample of the earthquakes input to ERFs for onshore PTHA analysis. We eliminate epistemic uncertainty related to source distribution under traditional PTHA analysis; in its place, a quantifiable, less subjective and generally smaller uncertainty related to the input to ERFs is included. The Nankai subduction zone is used as a test case, where previous ERFs have been conducted. Results indicate that the locations of representative earthquakes are sensitive to choice of magnitude–area relation and to whether a minimum cumulative stress objective is imposed on the fault. In general, incorporating ERFs into PTHA provide a physically self-consistent method to incorporate fault slip information in determining representative earthquakes for onshore PTHA, eliminating a major source of epistemic uncertainty.

Nankai subduction zone↗
USGS2026

Computation of regional groundwater budgets for the Virginia Coastal Plain aquifer system

Computation of detailed groundwater flow budgets for subdivisions of Virginia’s Coastal Plain aquifer system has enabled quantification and more thorough understanding of groundwater flow within this important water resource. A zone budget analysis conducted on previously published groundwater models of the Virginia Coastal Plain and Virginia Eastern Shore shows that groundwater conditions vary substantially throughout the Coastal Plain aquifer system due to local variations in hydrogeology and historical and ongoing variations in groundwater use and management. Decades of substantial groundwater withdrawal from the Coastal Plain aquifer system have fundamentally altered groundwater flow from pre-development conditions. Rates of sustainable withdrawal are limited because the downward groundwater flow rate into confined aquifers supplying groundwater is a relatively small portion of the total groundwater water budget for the aquifer system. Analyses of groundwater budgets from the Virginia Coastal Plain model show that groundwater flow is generally outward from the surficial aquifer to rivers and coastal water bodies and downward through a series of underlying aquifers and confining units to the Potomac aquifer, which is the deepest aquifer and the source of most groundwater withdrawals. Downward flow into the Potomac aquifer currently is estimated to be only 7 percent of total net precipitation-derived net recharge at the land surface but makes up about 66 percent of inflow to the aquifer in Virginia, with much of the remaining inflow occurring laterally from areas outside of defined groundwater budget regions in Virginia. For several decades prior to 2010, high rates of withdrawal from the Potomac aquifer resulted in substantial decline in groundwater storage in the aquifer and in most overlying aquifers and confining units. From 2010 to 2025, rates of withdrawal substantially lower than the historical maximum have resulted in small net increases in groundwater storage in the confined aquifer system for most regions of the Virginia Coastal Plain. Nevertheless, for the same period, groundwater storage for the entire model domain continues to incrementally decline, indicating that storage recovery in Virginia is offset by a continued decrease in storage in areas beneath the Chesapeake Bay or in adjacent areas of Maryland and North Carolina. Withdrawals from the Potomac aquifer have induced substantial downward flow which is a large part of groundwater budgets for confined aquifers such as the Potomac. Downward groundwater flow continues under current conditions, but because vertical flow rates are a function of the difference between water pressure in the upper surficial systems and lower confined units, those rates are lower than those in earlier decades as the confined water levels partially recover from larger groundwater withdrawals in the past. Geographically, groundwater flow is generally inward from perimeter regions of the Virginia Coastal Plain toward central regions with the largest withdrawal rates. Estimated groundwater inflow from coastal regions could be contributing to saltwater intrusion, though that was not measured directly in this study. Analyses of groundwater budgets from the Virginia Eastern Shore peninsula, a geographic region of the Virginia Coastal Plain, show that groundwater flow for that isolated aquifer system is generally outward from the surficial aquifer to coastal water bodies and downward into the confined Yorktown-Eastover aquifer system, which is the source of most withdrawals. Downward groundwater flow into the confined Yorktown-Eastover aquifer system is estimated to be less than 2 percent of total recharge and less than 9 percent of net recharge at the water table but makes up over 93 percent of all inflow to the confined aquifer system. Decades of substantial but relatively consistent groundwater withdrawals have induced greater downward flow rates into the confined aquifer system but also have resulted in loss of groundwater from storage. Currently, estimated storage loss accounts for slightly under 7 percent of withdrawals from the confined aquifer system. The current withdrawal rate from the confined Yorktown-Eastover system is near the highest reported rate for the Eastern Shore, which means that the storage depletion is expected to continue, even though groundwater levels appear to be relatively stable. Estimated groundwater flow rates upward from the confining unit underlying the Yorktown-Eastover system and small rates of inflow from coastal water bodies underscore ongoing concerns about up-coning and lateral intrusion of salty groundwater.

EarthArXiv↗
USGS2024

Representation of surface-water flows using Gradient-Related Discharge in an Everglades Network

The Everglades Depth Estimation Network interpolates water-level gage data to produce daily water-level elevations for the Everglades in south Florida. These elevations were used to estimate flow vectors (gradients and directions) and volumetric flow rates using the Gradient-Related Discharge in an Everglades Network (GARDEN) application developed by the U.S. Geological Survey in cooperation with the U.S. Army Corps of Engineers. Flow rates in both the east-west and north-south directions were computed on a 400-meter square grid using modified parameters in the Manning’s equation. The frictional resistance parameter in the Manning’s equation was calibrated to measured flow rates at coastal creeks fed by Everglades Depth Estimation Network boundary flows. Levees and other features that act as barriers to flow were defined as “no-flow” grid cells where vectors were set to zero. The flow volume magnitudes were calibrated with 2020 daily values of coastal river flows, and verification was performed using 2021 data. Within a given day, the measured coastal river flows fluctuate more than the GARDEN boundary flows because of tidal and wind forcings. Because the GARDEN boundary flows were the upstream water source for the coastal rivers, calibration focused on matching average daily flow volumes rather than daily fluctuations. The Pearson’s correlation coefficient is 0.766 for the 2020 calibration period and 0.566 for the 2021 verification period. Applying GARDEN to periods with hydraulic-control-structure releases allows the propagation of structure flows to be seen in the daily flow-vector maps along with the multiday response of flows farther downgradient. Flow vectors may be overestimated near control structures because of difficulties in resolving the water gradient downstream from the structure. Flow vectors farther from the structure are more accurate than those near the structure.

Florida↗
USGS2026

Computation of regional groundwater budgets for the Virginia Coastal Plain aquifer system

Computation of detailed groundwater flow budgets for subdivisions of the Virginia Coastal Plain aquifer system has enabled quantification and more thorough understanding of groundwater flow within this important water resource. A zone budget analysis based on previously published groundwater models of the Virginia Coastal Plain and Virginia Eastern Shore indicates that groundwater conditions vary substantially throughout the Coastal Plain aquifer system because of local variations in hydrogeology and historical and ongoing variations in groundwater use and management. Decades of substantial groundwater withdrawal from the Coastal Plain aquifer system have altered groundwater flow from predevelopment conditions. Rates of sustainable withdrawal are limited because the downward groundwater flow rate into confined aquifers is a relatively small part of the total groundwater budget for the aquifer system compared to the rate of recharge at the land surface. Analyses of groundwater budgets from the Virginia Coastal Plain model indicate that groundwater flow is generally outward from the surficial aquifer to rivers and coastal waterbodies and downward through a series of underlying aquifers and confining units to the Potomac aquifer, which is the deepest aquifer and the source of most groundwater withdrawals. Downward flow into the Potomac aquifer is estimated to be only 7 percent of total net precipitation-derived net recharge at the land surface but makes up about 66 percent of inflow to the aquifer in Virginia, with much of the remaining inflow occurring laterally from outside of defined groundwater budget regions in Virginia. For several decades prior to 2010, high rates of withdrawal from the Potomac aquifer resulted in substantial decline in groundwater storage in the aquifer and in most overlying aquifers and confining units. From 2010 to 2023, rates of withdrawal substantially lower than the historical maximum resulted in small net increases in groundwater storage in the confined aquifer system for most regions of the Virginia Coastal Plain. Nevertheless, for the same period, groundwater storage for the entire model domain continues to incrementally decline, indicating that storage recovery in Virginia is offset by a continued decrease in storage in areas beneath the Chesapeake Bay or adjacent areas of Maryland and North Carolina. Withdrawals from the Potomac aquifer have induced substantial downward flow which is a large part of groundwater budgets for confined aquifers such as the Potomac. For the most recent simulated conditions (2023) downward groundwater flow continues, but because vertical flow rates are a function of the difference between water pressure in the upper surficial systems and lower confined units, rates of downward flow are lower than those in earlier decades as the confined water levels partially recover from larger groundwater withdrawals in the past. Geographically, groundwater flow is generally inward from perimeter regions of the Virginia Coastal Plain toward central regions with the largest withdrawal rates. Groundwater inflow from coastal regions could be contributing to saltwater intrusion, even though that was not measured in this study. Analyses of groundwater budgets from the Virginia Eastern Shore peninsula, a geographic region of the Virginia Coastal Plain, indicate that groundwater flow for that isolated aquifer system is generally outward from the surficial aquifer to coastal water bodies and downward into the confined Yorktown-Eastover aquifer system, which is the source of most withdrawals. Downward groundwater flow into the confined Yorktown-Eastover aquifer system is estimated to be less than 2 percent of total recharge and less than 9 percent of net recharge at the water table but makes up more than 93 percent of all inflow to the confined aquifer system. Decades of substantial but relatively consistent groundwater withdrawals have induced greater downward flow rates into the confined aquifer system but also have resulted in loss of groundwater from storage. For the most recent simulated period (2023), estimated storage loss accounts for slightly under 7 percent of withdrawals from the confined aquifer system. The reported withdrawal rate for this period from the confined Yorktown-Eastover system is near the highest reported rate for the Virginia Eastern Shore, which means that the storage depletion is expected to continue, even though groundwater levels appear to be relatively stable. Estimated groundwater flow rates upward from the confining unit underlying the Yorktown-Eastover system and low rates of inflow from coastal water bodies underscore ongoing concerns about up-coning and lateral intrusion of salty groundwater.

Virginia↗
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↗
USGS2025

SURF: An automated method for building nonplanar 3D fault models from earthquake hypocenters

Accurately characterizing 3D fault geometry is vital for improving our understanding of earthquake behavior and informing the development of seismic hazard models. Despite their importance, subsurface fault structures tend to be poorly constrained because of limitations in observational data. Improvements to the seismic networks and earthquake detection algorithms have increased the precision and volume of earthquake catalogs, which help illuminate detailed subsurface fault structure and provide the most direct information available about fault geometries at depth. We present a Python package to automate generating 3D fault geometries directly from hypocentral seismicity patterns. This method begins with clustering events based on their spatial density, identifying coherent patterns. Nearby clusters are then merged based on the similarity of their orientations. We fit nonplanar surfaces using support vector regression to balance surface accuracy with minimal deviations from planarity. The fault models are output as quadrilateral meshes at user‐defined resolution. In the process of generating the 3D fault surfaces, we compute the spatial density of seismicity around the surface and the planarity as quantitative metrics of the model outputs. As a proof of concept, we apply this approach to the San Andreas–Calaveras fault junction region and the 2019 Ridgecrest earthquake sequence, both in California, which contain complex subparallel faults well defined at the Earth’s surface and abundant microseismicity. These case studies demonstrate the method’s ability to model complex fault structures, including long continuous fault surfaces, crossing faults, variably dipping segments, and subparallel faults. We test the method on both standard network catalogs and double‐difference relocated catalogs. We find that our seismicity‐based fault model results align with published 3D models that incorporate additional constraints and interpretations ( Plesch et al. , 2020 ; Aagaard and Hirakawa, 2021 ). This workflow provides a low‐user‐input solution for estimating fault geometries at depth from earthquake catalogs.

California↗
USGS1997

Application of acoustical methods for estimating water flow and constituent loads in Perdido Bay, Florida

Water flow and quality data were collected from December 1994 to September 1995 to evaluate variations in discharge, water quality, and chemical fluxes (loads) through Perdido Bay, Florida. Data were collected at a cross section parallel to the U.S. Highway 98 bridge. Discharges measured with an acoustic Doppler current profiler (ADCP) and computed from stage-area and velocity ratings varied roughly between + or - 10,000 cubic feet per second during a typical tidal cycle. Large reversals in flow direction occurred rapidly (less than 1 hour), and complete reversals (resulting in near peak net-upstream or downstream discharges) occurred within a few hours of slack water. Observations of simultaneous upstream and downstream flow (bidirectional flow) were quite common in the ADCP measurements, with opposing directions of flow occurring predominantly in vertical layers. Continuous (every 15 minutes) discharge data were computed for the period from August 18, 1995, to September 28, 1995, and filtered daily mean discharge values were computed for the period from August 19 to September 26, 1995. Data were not computed prior to August 18, 1995, either because of missing data or because the velocity rating was poorly defined (because of insufficient data) for the period prior to landfall of hurricane Erin (August 3, 1995). The results of the study indicate that acoustical techniques can yield useful estimates of continuous (instantaneous) discharge in Perdido Bay. Useful estimates of average daily net flow rates can also be obtained, but the accuracy of these estimates will be limited by small rating shifts that introduce bias into the instantaneous values that are used to compute the net flows. Instantaneous loads of total nitrogen ranged from -180 to 220 grams per second for the samples collected during the study, and instantaneous loads of total phosphorous ranged from -10 to 11 grams per second (negative loads indicate net upstream transport). The chloride concentrations from the water samples collected from Perdido Bay indicated a significant amount of mixing of saltwater and freshwater. Mixing effects could greatly reduce the accuracy of estimates of net loads of nutrients or other substances. The study results indicate that acoustical techniques can yield acceptable estimates of instantaneous loads in Perdido Bay. However, estimates of net loads should be interpreted with great caution and may have unacceptably large errors, especially when saltwater and freshwater concentrations differ greatly.

Florida↗
USGS2026

The 2025 Puerto Rico and Virgin Islands U.S. National Seismic Hazard Model Update: Ground motion model selection and comparison

We evaluate, select, and describe the ground-motion models (GMMs) used in the 2025 update of the U.S. National Seismic Hazard Model (NSHM) for Puerto Rico and the U.S. Virgin Islands (PRVI). We identify the most appropriate models that align with GMM selection criteria for use in the PRVI region to improve the accuracy of seismic hazard assessments. The update incorporates globally applicable GMMs suited for the active crustal and subduction earthquakes in the region. We include region-specific adjustments to these GMMs derived from local site response analyses derived from ground motion records. The unadjusted and regionally-corrected GMMs are combined to create a robust model for predicting median ground motion. The model integrates epistemic uncertainty through a median ground motion logic tree that accounts for variations in magnitude and distance. This study compares the GMMs selected for the 2025 PRVI NSHM, including both as-provided and regionally adjusted NGA-West2 and NGA-Subduction models, with those used in the 2003 PRVI NSHM. We evaluate how changes in model selection, weighting, aleatory variability, and epistemic uncertainty influence seismic hazard estimates. Trends with distance, magnitude, and spectral period are analyzed to evaluate how the scaling behavior of the newer GMMs differs from that of earlier models. Relative to the GMMs used in the 2003 NSHM for this region, the 2025 models generally predict lower ground motions. Comparisons with additional GMMs indicate that the adjustments applied for PRVI are consistent with regional-specific modifications developed elsewhere globally. The increase in aleatory variability and epistemic uncertainty in the 2025 update results in a notable increase in hazard levels from these wider uncertainty bounds. These changes can result in as much as a 10%–20% variation in probabilistic ground motion at the 2% in 50 years exceedance level for hazard maps computed across the region for representative site classes and periods.

Puerto Rico, Virgin Islands↗
USGS2025

Capturing the uncertainty of seismicity observations in earthquake rate estimates: Implications for probabilistic seismic hazard analysis and the USGS National Seismic Hazard Model

The rate of earthquakes in a region is a fundamental input to Probabilistic Seismic Hazard Analysis. We present a Monte Carlo method for computing that rate from seismicity catalogs while including a range of data and analysis uncertainties. This method is applied to regions for which the b value is assumed to be spatially invariant. Each region is broken down into epochs for which each epoch is estimated to have a uniform magnitude of completeness ( ⁠⁠ M c ). The distribution of earthquake rates for M ≥ M c is determined for each epoch by considering the Poisson likelihood of rates given the number of observed earthquakes with M ≥ M c ⁠ . We use a Monte Carlo process to include the uncertainty in b , ⁠ M c , and individual event magnitudes. The result for each epoch is the joint distribution of the Poisson rate of earthquakes with magnitudes larger than the minimum value used to calculate hazard ( ⁠⁠ M 1 ) and the Gutenberg–Richter b values, which control the extrapolation to other magnitudes. The rate for each region is either the duration‐weighted average over the epochs or, to better capture temporal variations, we also consider mixture models. The mixture models also provide an avenue to allow temporal variations in b values. To implement this joint distribution in a logic tree, we use the mean and 95% confidence branches, each of which is parameterized with an M ≥ M 1 rate and b value. We explore different ways of defining those branches, as well as non‐Gutenberg–Richter branches, and their impact on hazard estimates. The mean hazard, but not the fractiles, is robust with respect to these choices. To illustrate these new methods, we use synthetic data and catalogs from recent U.S. Geological Survey National Seismic Hazard Models for the Central and Eastern United States and for Puerto Rico and the U.S. Virgin Islands.

central and eastern United States, Puerto Rico, U.↗
USGS2026

Bathymetric maps, surface areas, and storage capacities of Council Grove Lake and Marion Reservoir, Kansas, and Pine Creek Lake, Oklahoma, 2024

The U.S. Geological Survey, in cooperation with the U.S. Army Corps of Engineers, completed high-resolution multibeam bathymetric surveys to compute new elevation-area-capacity tables for Council Grove Lake and Marion Reservoir, Kansas, and Pine Creek Lake, Oklahoma. Elevation-area-capacity tables identify the relation between the water-surface elevation, surface area, and storage capacity of the lake. The surface areas and storage capacities of each lake were computed from bathymetric surfaces combining multibeam echo sounder data collected in 2024 and light detection and ranging point-cloud data collected in 2016 and 2018.

Kansas, Oklahoma↗
USGS2026

Perspectives and insights gained from applying landscape ecological approaches to the Upper Mississippi riverscape

Advances in remote sensing, spatial data collection, hydraulic and geospatial models, and the maturation of landscape ecology during the early 2000s spurred the concept of rivers as landscapes or ‘riverscapes’. Since that time, the Upper Mississippi River Restoration (UMRR) Program - one of the world’s longest running large-river ecosystem restoration and monitoring programs - has applied landscape ecological principles to the Upper Mississippi River System (UMRS). I examine spatial patterns of aquatic nutrient concentrations, submersed aquatic vegetation, and fish and mussel communities, finding that within-river patterns are often patchy in nature, likely resulting from how nutrient processes and community traits interact with spatial variability in hydraulic exchange patterns. Local-scale patchy spatial patterns in the UMRS are embedded in system-scale longitudinal gradients and a hierarchical view of the river has provided a foundation for strategic river restoration planning. Key Words: Aquatic Habitat, hydraulic connectivity, patch, restoration, riverscape

BioScience↗
USGS2024

Probabilistic assessment of postfire debris-flow inundation in response to forecast rainfall

Communities downstream of burned steep lands face increases in debris-flow hazards due to fire effects on soil and vegetation. Rapid postfire hazard assessments have traditionally focused on quantifying spatial variations in debris-flow likelihood and volume in response to design rainstorms. However, a methodology that provides estimates of debris-flow inundation downstream of burned areas based on forecast rainfall would provide decision-makers with information that directly addresses the potential for downstream impacts. We introduce a framework that integrates a 24 h lead-time ensemble precipitation forecast with debris-flow likelihood, volume, and runout models to produce probabilistic maps of debris-flow inundation. We applied this framework to simulate debris-flow inundation associated with the 9 January 2018 debris-flow event in Montecito, California, USA. When the observed debris-flow volumes were used to drive the probabilistic forecast model, analysis of the simulated inundation probabilities demonstrates that the model is both reliable and sharp. In the fully predictive model, however, in which debris-flow likelihood and volume were computed from the atmospheric model ensemble's predictions of peak 15 min rainfall intensity, I 15 , the model generally under-forecasted the inundation area. The observed peak I 15 lies in the upper tail of the atmospheric model ensemble spread; thus a large fraction of ensemble members forecast lower I 15 than observed. Using these I 15 values as input to the inundation model resulted in lower-than-observed flow volumes which translated into under-forecasting of the inundation area. Even so, approximately 94 % of the observed inundated area was forecast to have an inundation probability greater than 1 %, demonstrating that the observed extent of inundation was generally captured within the range of outcomes predicted by the model. Sensitivity analyses indicate that debris-flow volume and two parameters associated with debris-flow mobility exert significant influence on inundation predictions, but reducing uncertainty in postfire debris-flow volume predictions will have the largest impact on reducing inundation outcome uncertainty. This study represents a first step toward a near-real-time hazard assessment product that includes probabilistic estimates of debris-flow inundation and provides guidance for future improvements to this and similar model frameworks by identifying key sources of uncertainty.

California↗
USGS2026

The damability function: A probabilistic approach to regional landslide dam susceptibility analysis applied to the Oregon Coast Range, USA

Landslides can dam rivers and require rapid response to mitigate catastrophic outburst floods. Here, we present a workflow to map landslide dam formation susceptibility at a regional scale. We define a probabilistic function that combines river valley width and landslide volume to efficiently determine the likelihood of a landslide dam or “damability”. We combine damability values with landslide susceptibility to estimate landslide dam susceptibility. The valley width measurements are automated using a new elevation threshold-based algorithm. Landslide volume is represented as a statistical distribution from mapped landslides. We validate and apply our approach to the Oregon Coast Range, USA and find that 36 % of river stretches exceed a dam potential threshold; these are in river headwaters and steeper terrain, which in this case correlate with more resistant lithologies. We also estimate volumes of the potential dammed lakes and find that most rivers with high dam susceptibility are less likely to impound large lakes because they have low drainage areas. However, widespread susceptibility, and the potential impacts from exceptionally large landslides, suggest that this hazard should be considered in the Pacific Northwest. The damability function workflow can ingest new data and be applied more broadly to assess future landslide dam hazards.

Oregon↗
USGS2026

Top-down targeted network analysis of critical mineral commodities applied to international geochemistry database

The global demand for critical mineral commodities is rapidly increasing, making domestic production an important factor in supporting the economy and national security. Large scale, publicly available geochemical databases allow for the application of data informatics methods to interrogate critical mineral commodities data for correlations in deposit formation and distribution, particularly for identifying enrichment of multiple critical mineral commodities at the same deposit. In this study, we applied network analysis to the Critical Minerals Mapping Initiative (CMMI) ore geochemistry (Critical Minerals in Ores, CMiO) database to identify the high concentration (defined as 10× bulk crustal abundance) co-occurrence of different critical mineral commodities across a mineral system hierarchy from deposit environments to individual deposits. Identifying patterns or unique outliers in enrichment in network communities will allow for the location of secondary critical mineral commodity resources from under-utilized deposits. We find trends in the enrichment of critical mineral commodities in network-communities between the elements praseodymium (Pr), neodymium (Nd), terbium (Tb), and dysprosium (Dy) across multiple CMiO database deposit environments and groups down to specific deposit types and sites. A separate trend in network community deposition is observed as well between iridium (Ir) and platinum (Pt) in deposit environments, groups, types, and sites. Network analysis focused on critical minerals in magmatic-hydrothermal deposits identified multiple deposit sites from different deposit types within the CMiO database with concentrations of Dy, Nd, Tb, Pr, Ir, and Pt that are at least ten times greater than the crustal average. This approach can be applied to any target element(s) or deposit(s) of interest, allowing broad investigation of co-enriched critical mineral commodities.

Journal of Geochemical Exploration↗
USGS1997

An optimized network for phosphorus load monitoring for Lake Okeechobee, Florida

Phosphorus load data were evaluated for Lake Okeechobee, Florida, for water years 1982 through 1991. Standard errors for load estimates were computed from available phosphorus concentration and daily discharge data. Components of error were associated with uncertainty in concentration and discharge data and were calculated for existing conditions and for 6 alternative load-monitoring scenarios for each of 48 distinct inflows. Benefit-cost ratios were computed for each alternative monitoring scenario at each site by dividing estimated reductions in load uncertainty by the 5-year average costs of each scenario in 1992 dollars. Absolute and marginal benefit-cost ratios were compared in an iterative optimization scheme to determine the most cost-effective combination of discharge and concentration monitoring scenarios for the lake. If the current (1992) discharge-monitoring network around the lake is maintained, the water-quality sampling at each inflow site twice each year is continued, and the nature of loading remains the same, the standard error of computed mean-annual load is estimated at about 98 metric tons per year compared to an absolute loading rate (inflows and outflows) of 530 metric tons per year. This produces a relative uncertainty of nearly 20 percent. The standard error in load can be reduced to about 20 metric tons per year (4 percent) by adopting an optimized set of monitoring alternatives at a cost of an additional $200,000 per year. The final optimized network prescribes changes to improve both concentration and discharge monitoring. These changes include the addition of intensive sampling with automatic samplers at 11 sites, the initiation of event-based sampling by observers at another 5 sites, the continuation of periodic sampling 12 times per year at 1 site, the installation of acoustic velocity meters to improve discharge gaging at 9 sites, and the improvement of a discharge rating at 1 site.

Florida↗
USGS2024

Methods to evaluate and improve the modeling of rupture directivity in assessment of seismic hazard

In recent years, there have been several advancements related to the modelling of near-source effects of earthquake rupture on strong ground shaking, leading to an improved characterization of ground motions and resulting seismic hazard. Some of these modifications have stemmed from physics-based numerical modelling of the earthquake rupture process, using physics-based dynamic rupture simulations. These contributions have led to a better understanding of how fault rupture characteristics, geometry, and the style of faulting can interact with the hypocenter-dependence on the path from source to site that may ultimately guide the development of seismic directivity models. Moving forward, the application of modern techniques can be used to incorporate these source characteristics and near-fault ground motion behavior that contribute to the azimuthally varying effects that result in rupture directivity. One example is the application of machine learning methods to support more automated integration of new predictor variables in model development and open more evaluation opportunities to access residuals. Here, we utilize several techniques to take advantage of the plethora of synthetic data and its ability to supplement preexisting trends observed in data. We showcase two examples of how models can be either developed, expanded upon, or constrained using artificial neural network model (ANNs). We evaluate the performance of the ANN with existing methods, comparing misfit, potential limitations, and ability to continue to improve upon these methods in the future. One approach uses a set of simulations with corresponding synthetic ground motions from the Southern California Earthquake Center (SCEC) CyberShake study to develop a ground motion model adapted to incorporate seismic directivity information using an ANN. This large database (TBs) enables us to train the model to capture magnitude, period, and distance variations and how these parameters relate to amplification from hypocenters located along finite-faults. In some cases, there is reduced misfit from better representing source features that aren’t included in base ground motion models that neglect hypocenter location (e.g. azimuthal variation, source-to-site terms). Another ANN method uses a shallow-layered neural network model to better fit a hypocenter-independent model. This method adjusts the median and aleatory variability to account for the averaged impact of various hypocenter distributions to fit the underlying directivity adjustment model. This method serves as a template to apply to other directivity models, improving computational efficiency and more readily enabling integration in hazard codes.

California↗
USGS2025

Groundwater salinity: Applying the specific conductance and water type proxy

Groundwater is increasingly needed for water supplies but may have limited utility in some locations because of its salinity. Salinity, often expressed as total dissolved solid (TDS), is frequently estimated using specific conductance (SC) measurements. However, the commonly used proxy (0.65 multiplied by SC to indicate TDS, common in many handheld meters) can result in inaccurate TDS estimates. First, the TDS–SC relationship is not linear over the entire concentration range of groundwater. Furthermore, the TDS (and salinity)–SC relationships vary substantially depending on the major-ion composition. Here we develop a proxy method utilizing SC and major-ion water type to estimate TDS and salinity specifically for groundwaters. Compared to most surface waters, groundwater tends to have a wider range of salinity (fresh to highly saline) and higher concentrations of bedrock-derived solutes such as carbonate ions, silica, and many other ions. The dataset used to develop the proxies includes water chemistry data from 149,059 discrete groundwater samples. The groundwater proxies, which employ nonlinear log–log relations, utilize five water types (HCO 3 , Cl, Ca-Mg-SO 4 , Na-K-SO 4 , and mixed waters), are accurate (median percent difference between TDS and salinity determined using the proxy compared to discrete measurements was <±0.8%) over a wide range of SC (up to 200 mS/cm), rapid, cost-effective, and can be measured on-site.

Groundwater↗