USGS Science⌕ Search

FIND YOUR NEXT DISCOVERY

Results for “Applied Computing and Geosciences”

Original records, connected by a shared subject.

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

1,223 records · Page 2Linked to original sources
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↗
USGS2025

False positives in the identification of dynamic earthquake triggering

Dynamic earthquake triggering is commonly identified through the temporal correlation between increased seismicity rates and global earthquakes that are possible triggering events. However, correlation does not imply causation. False positives may occur when unrelated seismicity rate changes coincidently occur at around the time of candidate triggers. We investigate the expected false positive rate in Southern California with global M ≥ 6 earthquakes as candidate triggers. We compute the false positive rate by applying the statistical tests used by DeSalvio and Fan (2023), https://doi.org/10.1029/2023jb026487 to synthetic earthquake catalogs with no real dynamic triggering. We find a false positive rate of ∼3.5%–8.5% when realistic earthquake clustering is present, consistent with the 95% confidence typically used in seismology. However, when this false positive rate is applied to the tens of thousands of spatial-temporal windows in Southern California tested in DeSalvio and Fan (2023), https://doi.org/10.1029/2023jb026487 , thousands of false positives are expected. The expected false positive occurrence is large enough to explain the observed apparent triggering following 70% of large global earthquakes (DeSalvio & Fan, 2023, https://doi.org/10.1029/2023jb026487 ), without requiring any true dynamic triggering. Aside from the known triggering from the nearby El Mayor-Cucapah, Mexico, earthquake, the spatial and temporal characteristics of the reported triggering are indistinguishable from random false positives. This implies that best practice for dynamic triggering studies that depend on temporal correlation is to estimate the false positive rate and investigate whether the observed apparent triggering is distinguishable from the correlations that may occur by chance.

JGR Solid Earth↗
USGS1998

Methods and guidelines for effective model calibration; with application to UCODE, a computer code for universal inverse modeling, and MODFLOWP, a computer code for inverse modeling with MODFLOW

This report documents methods and guidelines for model calibration using inverse modeling. The inverse modeling and statistical methods discussed are broadly applicable, but are presented as implemented in the computer programs UCODE, a universal inverse code that can be used with any application model, and MODFLOWP, an inverse code limited to one application model. UCODE and MODFLOWP perform inverse modeling, posed as a parameter-estimation problem, by calculating parameter values that minimize a weighted least-squares objective function using nonlinear regression. Minimization is accomplished using a modified Gauss-Newton method, and prior, or direct, information on estimated parameters can be included in the regression. Inverse modeling in many fields is plagued by problems of instability and nonuniqueness, and obtaining useful results depends on (1) defining a tractable inverse problem using simplifications appropriate to the system under investigation and (2) wise use of statistics generated using calculated sensitivities and the match between observed and simulated values, and associated graphical analyses. Fourteen guidelines presented in this work suggest ways of constructing and calibrating models of complex systems such that the resulting model is as accurate and useful as possible.

Water-Resources Investigations Report↗
USGS2026

pySATSI: A Python package for computing focal mechanism stress inversions

We introduce pySATSI, a Python package for computing earthquake focal mechanism stress inversions. This algorithm can handle a wide variety of types of stress inversion problems with a single script and can duplicate many capabilities of preceding methodologies. We also add new capabilities that include spatiotemporally variable inversion grids, damped stress estimates for clusters with few or no focal mechanisms, and variable fault‐plane ambiguities that the user can assign to individual events. In addition, we added the ability to use damped stress inversions with fault‐plane ambiguity probabilities that are weighted by fault instabilities. Our algorithm is computationally efficient with faster runtimes than previous algorithms, scales well for large datasets, and can be easily parallelized.

Seismological Research Letters↗
USGS2010

Divisions of geologic time—Major chronostratigraphic and geochronologic units

Effective communication in the geosciences requires consistent uses of stratigraphic nomenclature, especially divisions of geologic time. A geologic time scale is composed of standard stratigraphic divisions based on rock sequences and is calibrated in years. Over the years, the development of new dating methods and the refinement of previous methods have stimulated revisions to geologic time scales. Advances in stratigraphy and geochronology require that any time scale be periodically updated. Therefore, Divisions of Geologic Time, which shows the major chronostratigraphic (position) and geochronologic (time) units, is intended to be a dynamic resource that will be modified to include accepted changes of unit names and boundary age estimates. This fact sheet is a modification of USGS Fact Sheet 2007-3015 by the U.S. Geological Survey Geologic Names Committee.

Fact Sheet↗
USGS2026

Bascom Laser Diffraction Sedimentology Laboratory, Reston, Virginia

Introduction At the Bascom Laser Diffraction Sedimentology Laboratory, which is located in the Florence Bascom Geoscience Center at U.S. Geological Survey (USGS) headquarters in Reston, Virginia, scientists use physical sedimentology and particle characterization techniques to conduct detailed sediment characterization. Scientists address research problems in collaboration with other USGS science centers, State geological surveys, commercial industry, universities, and other partners. Laboratory capabilities include laser diffractometry for quantitative particle-size analysis, portable x-ray fluorescence (XRF) analysis for determining elemental abundances in rock or sediment samples, petrographic analysis of geologic media, and mechanical sieve analysis. These methods are used to analyze soil and sediment core material from terrestrial, marine, and lacustrine environments, surface sediments from coastal regions, and calcareous materials. Work done by the laboratory supports geologic mapping, resource assessments, land change studies, and geohazard analyses.

Fact Sheet↗
USGS2026

Orientation dependence of probabilistic seismic hazard estimates from CyberShake physics-based simulations

Earthquake ground‐motion intensities, such as pseudospectral accelerations (SAs), can vary significantly with horizontal orientation. However, conducting probabilistic seismic hazard analysis (PSHA) for each horizontal orientation is challenging because current ground‐motion models used in PSHA consider only a single horizontal intensity value, usually the median across all orientations, known as RotD50. To address this limitation, we employ physics‐based simulations for PSHA, which contain full waveforms from which ground‐motion intensities can be computed for all horizontal orientations to study directional seismic hazard. We apply our approach to the latest CyberShake study of the Greater Los Angeles metropolitan area, developed by the Statewide California Earthquake Center, finding that seismic hazard at a 2475‐yr return period, a common value used for earthquake‐resistant design, varies significantly with horizontal orientation. For instance, for SAs at 3 s, the maximum seismic hazard across all horizontal orientations is, on average, 15% higher than the median RotD50 hazard, with these differences becoming more pronounced at longer periods. These observed variations can generally be attributed to physical mechanisms that polarize seismic waves, such as the radiation pattern of the earthquake source and the influence of the subsurface structure. These results may have important implications for earthquake engineering applications, particularly for long‐period structures in areas with substantial horizontal variations in seismic hazard.

California↗
USGS2025

Critical Minerals in Ores (CMiO) database

Critical minerals are commodities essential to modern industrial and strategic technologies and are highly vulnerable to supply chain disruption. The Critical Minerals Mapping Initiative (CMMI) is a collaboration among the U.S. Geological Survey (USGS), the Geological Survey of Canada, and Geoscience Australia that aims to deepen global understanding of where critical minerals are located. A key output of this initiative is the Critical Minerals in Ores (CMiO) database that is advancing our collective understanding of critical minerals distributions. For instance, publicly available data on the concentrations of many critical minerals are sparse because these commodities can only be produced in small, yet essential, quantities compared to the primary commodities like copper and zinc. The CMiO database helps bridge this gap by offering high-quality, multielement geochemical data from a wide variety of critical mineral-bearing deposits around the world. Importantly, it uses a novel consensus deposit environment, group, and type classification scheme developed by the agencies that allows comparisons among ore deposits from different regions. The CMiO database contains geochemical data for more than 20,000 samples from more than 100 deposit types comprising 10 deposit environments.

Fact Sheet↗
USGS2025

Total uncertainty quantification in inverse solutions with deep learning surrogate models

We propose an approximate Bayesian method for quantifying the total uncertainty in inverse partial differential equation (PDE) solutions obtained with machine learning surrogate models, including operator learning models. The proposed method accounts for uncertainty in the observations, PDE, and surrogate models. First, we use the surrogate model to formulate a minimization problem in the reduced space for the maximum a posteriori (MAP) inverse solution. Then, we randomize the MAP objective function and obtain samples of the posterior distribution by minimizing different realizations of the objective function. We test the proposed framework by comparing it with the iterative ensemble smoother and deep ensembling methods for a nonlinear diffusion equation with an unknown space-dependent diffusion coefficient. Among other applications, this equation describes the flow of groundwater in an unconfined aquifer. Depending on the training dataset and ensemble sizes, the proposed method provides similar or more descriptive posteriors of the parameters and states than the iterative ensemble smoother method. Deep ensembling underestimates uncertainty and provides less-informative posteriors than the other two methods. Our results show that, despite inherent uncertainty, surrogate models can be used for parameter and state estimation as an alternative to the inverse methods relying on (more accurate) numerical PDE solvers.

Journal of Computational Physics↗
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↗
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↗
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 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↗
USGS2026

USGS Geochron—A database of geochronological and thermochronological dates and data—Technical documentation

Geochronological and thermochronological data are essential for constraining the timing and rates of geological processes, supporting geologic mapping, natural hazard assessment, and resource exploration. The U.S. Geological Survey (USGS) Geochron Database is a centralized, relational database that integrates USGS and State geological survey data in accordance with the National Geologic Mapping Act of 1992 and its reauthorizations. This report documents the structure and development of the database, including its standardized schema, controlled vocabularies, and protocols for compiling legacy and newly published data. The database is designed to adhere to FAIR (Findable, Accessible, Interoperable, and Reusable) data principles and currently (2026) includes data from a wide range of geochronological and thermochronological methods. Public access is provided through the USGS Geochron Database Explorer, a browser-based geographic information system (GIS) interface, and through versioned data releases available on ScienceBase, which provide sample-level summary data and detailed analytical information. The USGS Geochron Database is structured to be extensible, supporting the inclusion of additional methods and data types as they are compiled. This format ensures long-term utility and aligns with the Mapping Act’s directive to create a national archive of geochronological information that adds interpretive value to geologic map data. Its design reflects a commitment to data transparency, scientific reproducibility, and national geoscience priorities.

Data Report↗
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.↗