USGS Science⌕ Search

SEARCH · USGS Science

Results for “Permafrost and Periglacial Processes”

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,660 records · Page 79Linked to original sources
USGS1998

Ground-water flow in the surficial aquifer system and potential movement of contaminants from selected waste-disposal sites at Cecil Field Naval Air Station, Jacksonville, Florida

As part of the Installation Restoration Program, Cecil Field Naval Air Station, Jacksonville, Florida, is considering remedialaction alternatives to control the possible movement of contaminants from sites that may discharge to the surface. This requires a quantifiable understanding of ground-water flow through the surficial aquifer system and how the system will respond to any future stresses. The geologic units of interest in the study area consist of sediments of Holocene to Miocene age that extend from land surface to the base of the Hawthorn Group. The hydrogeology within the study area was determined from gamma-ray and geologists' logs. Ground-water flow through the surficial aquifer system was simulated with a seven-layer, finite-difference model that extended vertically from the water table to the top of the Upper Floridan aquifer. Results from the calibrated model were based on a long-term recharge rate of 6 inches per year, which fell in the range of 4 to 10 inches per year, estimated using stream hydrograph separation methods. More than 80 percent of ground-water flow circulates within the surficial-sand aquifer, which indicates that most contaminant movement also can be expected to move through the surficial-sand aquifer alone. The surficial-sand aquifer is the uppermost unit of the surficial aquifer system. Particle-tracking results showed that the distances of most flow paths were 1,500 feet or less from a given site to its discharge point. For an assumed effective porosity of 20 percent, typical traveltimes are 40 years or less. At all of the sites investigated, particles released 10 feet below the water table had shorter traveltimes than those released 40 feet below the water table. Traveltimes from contaminated sites to their point of discharge ranged from 2 to 300 years. The contributing areas of the domestic supply wells are not very extensive. The shortest traveltimes for particles to reach the domestic supply wells from their respective contributing areas ranged from 70 to 200 years.

Florida↗
USGS2026

Ice sheet dynamics drive pronounced changes in the subsurface freshwater-saltwater interface

Saltwater is migrating into freshwater aquifers globally with water quality and biogeochemical implications, yet saltwater intrusion in glaciated regions is sparsely investigated. Field observations suggest that groundwater head in glaciated systems is influenced by ice sheet forcings and provides evidence that seawater infiltrated into offshore aquifers during past deglaciation events. To understand links between ice sheet dynamics, groundwater head, and saltwater intrusion, we use numerical models to explore the effects of deglaciation on nearshore head and salinity distributions. We find that ice sheet thinning diminishes groundwater head, and the resulting shift in subsurface pressure gradients drives rapid landward movement of the subsurface freshwater-saltwater interface up to 4.0 km or 1.3 m per m ice sheet loss. Results highlight an overlooked saltwater intrusion mechanism that aligns with field observations and affects glaciated coastlines undergoing ice sheet retreat, underscoring the need to consider this mechanism in studies of contemporary coastal water quality.

Geophysical Research Letters↗
USGS2024

On connecting hydro-social parameters to vegetation greenness differences in an evolving groundwater-dependent ecosystem

Understanding groundwater-dependent ecosystems (i.e., areas with a relatively shallow water table that plays a major role in supporting vegetation health) is key to sustaining water resources in the western United States. Groundwater-dependent ecosystems (GDEs) in Colorado have non-pristine temporal and spatial patterns, compared to agro-ecosystems, which make it difficult to quantify how these ecosystems are impacted by changes in water availability. The goal of this study is to examine how key hydrosocial parameters perturb GDE water use in time and in space. The temporal approach tests for the additive impacts of precipitation, surface water discharge, surface water mass balance as a surrogate for surface–groundwater exchange, and groundwater depth on the monthly Landsat normalized difference vegetation index (NDVI). The spatial approach tests for the additive impacts of river confluences, canal augmentation, development, perennial tributary confluences, and farmland modification on temporally integrated NDVI. Model results show a temporal trend (monthly, 1984–2019) is identifiable along segments of the Arkansas River at resolutions finer than 10 km. The temporal impacts of river discharge correlate with riparian water use sooner in time compared to precipitation, but this result is spatially variable and dependent on the covariates tested. Spatially, areal segments of the Arkansas River that have confluences with perennial streams have increased cumulative vegetation density. Quantifying temporal and spatial dependencies between the sources and effects of GDEs could aid in preventing the loss of a vulnerable ecosystem to increased water demand, changing climate, and evolving irrigation methodologies.

Colorado↗
USGS2026

Groundwater-level elevations in the bedrock aquifers of the Denver Basin aquifer system, Elbert County, Colorado, 2015–23

Water users in Elbert County, Colorado, rely on groundwater from bedrock aquifers in the Denver Basin aquifer system (upper Dawson, lower Dawson, Denver, Arapahoe, and Laramie-Fox Hills aquifers) for approximately half of their water uses. Withdrawals from the bedrocks aquifers have increased to meet the water use needs of expanding regional population growth and development. The U.S. Geological Survey, in cooperation with the Elbert County Board of County Commissioners, began a study in 2015 to monitor groundwater levels within Elbert County. The primary purpose of this report is to present a summary of groundwater levels measured during the study period (2015–23) and present results from statistical analyses of changes in groundwater-level elevations through time. Discrete groundwater levels were measured at 36 wells within Elbert County. Seven of those wells contained equipment to make and record continuous groundwater-level measurements at hourly intervals. All aquifers, except the lower Dawson aquifer, had only declining groundwater-level elevations in discrete measurements for wells with statistically significant trends. Of the eight statistically significant trends in the lower Dawson aquifer, two wells indicated increasing groundwater-level elevation from discrete measurements. The groundwater-level elevation trend medians in the upper Dawson, lower Dawson, Denver, Arapahoe, and Laramie-Fox Hills aquifers were −0.23, −0.66, −0.64, −0.39, and −0.63 feet per year, respectively, for discrete groundwater-level elevation measurements. Trends in continuous groundwater-level elevations were in agreement with statistically significant trends in discrete groundwater-level elevations for all wells. The groundwater-level elevation trend medians in this study, compared to the overall trends in a 2015−2018 study, both indicated declining groundwater-level elevations except in the upper Dawson aquifer, where the trend direction was opposite, a positive trend from 2015 to 2018 and a negative trend (declining groundwater elevations) from 2015 to 2023. The change in trends within the upper Dawson aquifer may be affected by differences in the study period and the trend analysis applied. Trends during the 2015–23 study period were compared to departures from the median 2015 groundwater-level elevation for each site in each aquifer. In general, the departures from the 2015 median supported trends observed at each site and correlated spatially with greater departures near the western border of Elbert County. Additionally, 30-year precipitation data showing wet and dry periods were overlaid with the departure from the 2015 median to assess groundwater-level patterns in wells in the five aquifers. Departures from the 2015 median groundwater-level elevations appeared greatest during the dry period between 2020 and 2023. Potentiometric-surface maps of the upper and lower Dawson aquifers created from static April 2023 groundwater elevations indicated groundwater-flow direction is generally from the south to the north. Results of this study could be used to guide additional groundwater monitoring in Elbert County and could aid in long-term planning of water resources.

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

A synthesis engine for constructing geologic maps of the United States

The geologic history of the United States is cataloged in thousands of geologic maps produced during many decades. However, the disparate nature of these individual maps makes it challenging to assess resources, research geologic histories, or characterize natural hazards holistically across the Nation. The U.S. House of Representatives 2020 appropriations bill for the U.S. Department of the Interior (H.R. 116-100) requires the U.S. Geological Survey to “bring together detailed national and continental-resolution [two-dimensional] and [three-dimensional] information produced throughout the Survey and by [F]ederal and [S]tate partners.” In response to this directive, this report presents a compilation and synthesis of geologic maps across the United States in the form of a relational database. The synthesis database includes thematic maps that synthesize the Nation’s geology, and retains the original input maps as well as linkages to standardized vocabularies to aid the discoverability of geologic information. Specifically, the synthesis database is targeted toward producing four National-resolution maps for the conterminous United States: Quaternary geology, the geology at the Earth’s surface, pre-Quaternary geology, and Precambrian geology. In addition, the synthesis database includes the infrastructure necessary to expand to additional resolutions in the future.

conterminous United States↗
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

Estimating earthquake source depth using teleseismic broadband waveform modeling at the USGS National Earthquake Information Center

The U.S. Geologic Survey National Earthquake Information Center (NEIC) monitors global seismicity, producing a catalog of earthquake source parameters in near-real-time to provide information that can help mitigate the societal impact of earthquakes. The NEIC commonly relies on teleseismic observations to constrain earthquake source parameters (e.g., location, depth, magnitude, and mechanism) due to a lack of local and regional observations. For these ‘teleseismic’ events, depth phase (i.e., pP , sP ) arrival time observations provide the best estimate on source depth. However, depth phases are often difficult to accurately identify and/or pick. Therefore, NEIC relies on waveform modeling, such as those determined from W-phase ( M ww ), body wave ( M wb ), and regional ( M wr ) moment tensor estimations, to provide constraints on source depth. While depth estimates from these approaches are informative, higher frequency observations provide more precise estimates because depth phases are more prominently observed at higher frequencies. Here, we present NEIC’s relatively high-frequency (~0.04 to 1 Hz) teleseismic waveform modeling approach, termed Synthetic Depth Phase Modeling (SynDepth), for determining source depth. SynDepth was developed to provide NEIC with a tool that enables rapid, accurate, and quantifiable estimates of earthquake source depth in cases where locator depths are not reliable. This relatively simple and fast procedure searches over 1 km-incremented source depths and an expanding triangular source-time function to find the best-fitting solution. We compare automatic SynDepth solutions for a dataset of 1,216 earthquakes (M5.5-M7.6) between 2017 and 2021 to NEIC-derived depth estimates from other methods. Our approach provides a robust depth estimate for earthquakes lacking local arrival time data, and it minimizes the need for analyst review of depth-phase picks ( pP , sP ) or using predefined ‘fixed’ depths.

Seismological Research Letters↗
USGS2024

Debris avalanches in the northern California Coast Range triggered by plate boundary earthquakes

Determining the timing and cause for ancient hillslope failures proves difficult in the western United States, yet critical as it ties directly into groundmotion estimates for hazardous events. This knowledge gap is important to confront as hillslope failures are candidates to be triggered by earthquakes along active plate boundaries. We identify two prehistoric, i.e., preinstrumental history, debris avalanches (3 10^6 and 6 10^6 m^3) in the Coast Ranges of northern California. These debris avalanches are well suited for studying the timing (to the exact year) and trigger as the densely forested landscape enables effective use of dendrochronology and high-resolution radiocarbon to pinpoint time of failure, and the legacy of landslide studies in the region provides context for evaluating a seismic trigger as the most likely triggering mechanism. Neither debris avalanche was triggered by the largest precipitation event of the instrumental record in 1964 CE, and the two debris avalanche sites are physiographically suited to accommodate topographic amplification of seismic shaking. Through a suite of geologic and dendrochronologic evidence, we establish the likely time of failure of the two debris avalanches as 1906 CE and 933 CE. The failure of the younger landslide coincides with the 1906 CE San Francisco earthquake on the Pacific/North American plate boundary, while the older landslide dates to 933 CE falling within the broad age range (850-966 CE) of a likely megathrust earthquake on the Cascadia subduction zone as recorded in coastal marshes in Humboldt Bay, California. The precise age on the 933 CE debris avalanche could trim the uncertainty on a Cascadia subduction zone event from 116 years to one year. Utilizing debris avalanche records from sites suitable to record seismic shaking improves understanding of plate-boundary earthquake timing and extent of shaking.

Bulletin of the Seismological Society of America↗
USGS2025

Foundational uncertainties in terminal Ediacaran chronostratigraphy revealed by high-precision zircon U-Pb geochronology of the Nama Group, Namibia

The Nama Group of southern Namibia and northwestern South Africa hosts the best-dated mixed carbonate-siliciclastic foreland basin succession of the terminal Ediacaran [ca. 551 million years (Ma) ago to <538 Ma] and is key for resolving the chronology of early metazoan evolution. Numerous silicified volcanic tuff interbeds are present, but differing interpretations regarding the fidelity of their ages lead to different regional stratigraphic correlations, especially for the Urusis Formation of the Schwarzrand Subgroup. An expanded record of the Urusis Formation is found in the Swartpunt area of southern Namibia, which has yielded an important metazoan biota. But the succession in this area is preserved as a series of thrusts at the leading edge of the Gariep orogenic belt and zircon U-Pb data show systematic age repetition. We use regional stratigraphic and structural mapping, integrated with carbonate carbon isotope (δ 13 C carb ) chemostratigraphy and high-precision radioisotope U-Pb zircon geochronology from outcrop and recently acquired drill core to develop a temporally calibrated basin-wide depositional model. This integrated dataset either reflects complex zircon reworking, inheritance, or potential analytical biases (Scenario 1) or the presence of a Gariep-related cryptic décollement within the Spitskop Member that has resulted in stratigraphic repetition (Scenario 2). We investigate the evidence for and against both scenarios and consider their implications for stratigraphic and δ 13 C carb correlations between the Swartpunt area and coeval exposures along the Orange River border with South Africa. Given that these issues are in an area that hosts numerous silicified ash beds and extensive exposure, an inability to confidently discount either scenario highlights a level of compounding uncertainty in zircon U-Pb geochronology that must be considered when attempting to build global chronostratigraphic frameworks. Scenario 1 implies that some of the weighted mean ages and Bayesian eruption ages from the Swartpunt area may be >1 Myr older than the depositional age of their respective ash beds when assuming existing stratigraphic correlations. If this scenario is preferred, then a cautious approach would be to consider all weighted mean zircon U-Pb ages from ash beds to reflect maximum depositional ages. Both scenarios support deposition of the Huns Member >540 Ma in the Swartpunt area if the oldest weighted mean age reported here represents a near-depositional age, which has significant implications for the temporal calibration of important terminal Ediacaran ichnofossil assemblages and future cyclostratigraphic studies. Stratigraphic correlations common to both scenarios allow us to temporally calibrate a basin evolution model for the Nama Group. Temporal trends in initial hafnium isotope (εHf) compositions of zircon grains from ash beds throughout the succession may support progressive crustal thickening associated with underplating of the Damara orogenic belt along the northern periphery of the Kalahari craton from ca. 547 Ma to ca. 538 Ma. The compilation of new and published zircon U-Pb ages may also imply that the locus of carbonate platform development migrated from north to south (present co-ordinates), tracking the migration of foredeep subsidence.

Earth-Science Reviews↗
USGS2025

Characterizing the scale of regional landslide triggering from storm hydrometeorology

Rainfall strongly affects landslide triggering; however, understanding how storm characteristics relate to the severity of landslides at the regional scale has thus far remained unclear, despite the societal benefits that would result from defining this relationship. As mapped landslide inventories typically cover a small region relative to a storm system, here we develop a dimensionless index for landslide-inducing rainfall, A * , based on extremes of modeled soil water relative to its local climatology. We calibrate A * using four landslide inventories, comprising over 11 000 individual landslides over four unique storm events, and find that a common threshold can be applied to estimate regional shallow-landslide-triggering potential across diverse climatic regimes in California (USA). We then use the spatial distribution of A * , along with topography, to calculate the landslide potential area (LPA) for nine landslide-inducing storm events over the past 20 years, and we test whether atmospheric metrics describing the strength of landfalling storms, such as integrated water vapor transport, correlate with the magnitude of hazardous landslide-inducing rainfall. We find that although the events with the largest LPA do occur during exceptional atmospheric river (AR) storms, the strength of landfalling atmospheric rivers does not scale neatly with landslide potential area, and even exceptionally strong ARs may yield minimal landslide impacts. Other factors, such as antecedent soil moisture driven by storm frequency and mesoscale precipitation features within storms, are instead more likely to dictate the patterns of landslide-generating rainfall throughout the state.

California↗
USGS2025

A soil velocity model for improved ground motion simulations in the U. S. Pacific Northwest

Near-surface seismic velocity structure may significantly impact the intensity, duration, and frequency content of ground shaking during an earthquake. In this study, we compile 649 shear wave velocity (Vs) profiles throughout the U.S. Pacific Northwest and southern British Columbia (PNW) and use these measured profiles to develop a representative soil velocity model for four major Holocene soil provinces: Puget Lowlands, Willamette Valley, fill and alluvium, and `other' soils. The resulting soil velocity model shows good agreement to measured data for a wide range of site conditions, with variability between different geologic domains reflecting fundamental differences in depositional environments. We then show that using this regional soil velocity model in simulations of the 2001 M6.8 Nisqually, Washington earthquake improves the fit to observed high-frequency (≥ 0.5 Hz) ground motions in the Puget Sound region compared to simulations that do not incorporate shallow (≤ 200 m) seismic velocity structure. Overall, this work shows that incorporating localized soil velocity profiles into seismic velocity models is important for accurately estimating high-frequency ground motion and regional seismic hazard in earthquake simulations. Future earthquake simulations and hazard studies in the PNW could incorporate these soil velocity profiles to capture the region's distinct site response characteristics.

Washington↗
USGS2024

Turbidite correlation for paleoseismology

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

Oregon, Washington↗
USGS2024

Parsimonious high-resolution landslide susceptibility modeling at continental scales

Landslide susceptibility maps are fundamental tools for risk reduction, but the coarse resolution of current continental-scale models is insufficient for local application. Complex relations between topographic and environmental attributes characterizing landslide susceptibility at local scales are not transferrable across areas without landslide data. Existing maps with multiple susceptibility classifications under-represent landslide potential in moderate and gently sloping terrain. We leverage an extensive landslide database ( N = 613,724), a high-resolution digital elevation model (10-m), and high-performance computing resources, to develop a new nationwide susceptibility map for the contiguous United States, Hawaii, Alaska, and Puerto Rico. We calculate four alternative linear and nonlinear thresholds of topographic slope and relief using an objective split-sample calibration. We down-sample our results to a 90-m grid to account for uncertainty in the digital elevation model and landslide position, and evaluate these thresholds' ability to differentiate areas of greater susceptibility. The less conservative nonlinear model optimally balances our priorities of capturing observed landslides (99%) while minimizing area covered by susceptible terrain (43%). Independent evaluation with four statewide landslide inventories ( N = 172,367) reinforces our model selection but highlights spatially variable performance. Therefore, we propose a novel approach to susceptibility classification using the concentration of landslide-prone terrain within each down-sampled grid. While landslides are possible within any cells containing susceptible terrain, those with the highest concentration capture the majority of observed landslides. Our new map characterizes landside susceptibility more consistently than prior models; our transparent classification approach also provides flexibility for accommodating different tolerances in risk reduction measures.

AGU Advances↗
USGS2024

Parameter ESTimation with the Gauss–Levenberg–Marquardt algorithm: An intuitive guide

In this paper, we review the derivation of the Gauss–Levenberg–Marquardt (GLM) algorithm and its extension to ensemble parameter estimation. We explore the use of graphical methods to provide insights into how the algorithm works in practice and discuss the implications of both algorithm tuning parameters and objective function construction in performance. Some insights include understanding the control of both parameter trajectory and step size for GLM as a function of tuning parameters. Furthermore, for the iterative Ensemble Smoother (iES), we discuss the importance of noise on observations and show how iES can cope with non-unique outcomes based on objective function construction. These insights are valuable for modelers using PEST, PEST++, or similar parameter estimation tools.

Groundwater↗
USGS2026

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

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

Earthquake Spectra Journal↗
USGS2025

Updating regional‐scale geospatial liquefaction models with locally available geotechnical data

We present a method to update the geospatial liquefaction model used by the U.S. Geological Survey’s near‐real‐time ground failure product with subsurface geotechnical data. The geospatial model estimates liquefaction probability from peak ground velocity (via ShakeMap) and geospatial susceptibility proxies. In many regions, additional information relevant to constraining liquefaction likelihood is also available, including surface geology maps and subsurface geotechnical measurements. There is currently no mechanism to use these data in the ground failure product liquefaction model, even though these data could provide more precise constraints on spatial variations in the lithologic character of the soil (surface geology) and direct measurements of the subsurface mechanical properties that affect liquefaction occurrence and severity (geotechnical measurements). In this study, we develop a method to integrate these data with the geospatial model and assess how these data can improve regional‐scale predictions. We develop a Bayesian updating framework and apply it to the 1989 magnitude 6.9 Loma Prieta, California, earthquake, for which mapped observations are available to evaluate performance. We constrain the Bayesian framework with 373 Northern California cone penetration tests and liquefaction susceptibility classes based on the mapped surface geology. This Bayesian model incorporates geotechnical information into the geospatial model and more accurately predicts liquefaction occurrences than the geospatial model, while sacrificing less accuracy in terms of predicting the absence of liquefaction than the geotechnical model. In future applications, this approach could be adapted to update other geospatial models using locally available subsurface data.

California↗
USGS2024

Agricultural return flow dynamics on a reach of the East River, Colorado, as assessed by mass balance

The U.S. Geological Survey, in cooperation with the Upper Gunnison River Water Conservancy District, studied historical streamflow in a reach of the East River, Colorado, to gain a preliminary understanding of return flow dynamics. Return flow is agricultural irrigation water that is not consumed by evapotranspiration and instead reaches streams by surface and subsurface flow paths. The study reach had a contributing area of 50 square miles and contained 5.23 square miles of pastures irrigated with water diverted from the East River and its tributaries. By comparing upstream inflows to downstream outflows, the net water balance of the study reach from 1994 to 2023 was assessed. Two general hydrologic conditions for the study reach were identified. One hydrologic condition was characterized by a net loss or consumption of water, termed here as general deficit. This general deficit condition extended about 16 years, from 1997 to 2012. During general deficit years, there was usually a notable net loss of streamflow from April through July, and a small net gain, possibly related to return flows, occurred in August about 75 days after the minimums for losses. The second hydrologic condition was characterized by a net gain of water, termed here as general surplus. This second condition extended about 10 years, from 2014 to 2023. During general surplus years, two separate transitions from net loss to net gain commonly occurred during June through August. Losses during general surplus years were smaller than losses during general deficit years, the respective gains were larger, and times between losses and gains were about 18 and 22 days. Differences between the two hydrologic conditions could reflect interactions among irrigation water, available capacity to store additional shallow groundwater, and streamflow. However, deciphering the causes for the shifts between the two general hydrologic conditions was beyond the scope of this report.

Colorado↗