USGS Science⌕ Search

SEARCH · USGS Science

Results for “Computers & 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.

97 records · Page 5Linked to original sources
USGS2024

A partially nonergodic ground-motion model for Fourier amplitude spectra for the San Francisco Bay area, California, USA

We develop a partially nonergodic ground-motion model (GMM) for Fourier amplitude spectra for the San Francisco Bay Area, California, USA, using the Bayless and Abrahamson (2019) GMM as a reference ergodic GMM and developing location-dependent adjustments to the predicted median and variance. We compile regional ground-motion data from moment magnitude (𝑀 w ) >3 earthquakes occurring during 2000–2022 for which magnitude information is available in the U.S. Geological Survey Comprehensive Catalog (Guy et al., 2015). The data set predominantly consists of records from 𝑀 w 3.5–4.5 earthquakes but includes three well-recorded 𝑀 w > 5 events. Ground-motion residuals are evaluated using the time-averaged shear-wave velocity in the top 30 m (𝑉 S30 ) from the California-specific map of Thompson et al. (2018) and basin-depth site parameters from the seismic velocity model of Aagaard and Hirakawa (2021). The 𝑉 S30 dependence and basin-depth scaling of the reference ergodic GMM of Bayless and Abrahamson (2019) are evaluated and modified with the updated data set. We compute maps of site adjustments using a varying-coefficient model that considers the spatial correlation structure and uncertainties at each observation location. The spatial covariance model is developed using ground-motion residuals that are standardized by the uncertainty model, which allows for consideration of the aleatory variability in developing the site adjustments. The covariance model is fit considering the means and standard deviations of the site terms at all locations. The use of partially nonergodic median adjustments results in modified variance components of the within-event variability. Due to the low number of large-magnitude earthquakes that control seismic hazard in the data set, we do not modify between-event variance; however, we present adjustments to site-to-site variability for use in partially nonergodic hazard assessments.

California↗
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↗
USGS2024

Estimation of reservoir storage capacity and geomorphic change detection analysis from a multibeam bathymetric survey of Randy Poynter Lake, Rockdale County, Georgia

Rockdale County Department of Water Resources has a directive to update estimates of the reservoir storage capacity of Randy Poynter Lake, located in northern Georgia, and to assess recent sedimentation and associated storage capacity loss. In 2022, the U.S. Geological Survey completed a multibeam bathymetric survey of Randy Poynter Lake to update storage capacity estimates and to quantify storage capacity change since the first multibeam bathymetric survey in 2012 in consideration of estimated errors inherent to bathymetric surveys. Data from the 2022 survey were used to generate contours of the reservoir as well as compute storage capacity at regular increments of water-surface elevation. Storage capacity comparisons between 2012 and 2022 at Randy Poynter Lake show minimal changes that are within the estimated uncertainties, with consistent or slightly increased storage capacities observed at most water-surface elevations and reductions observed at the remaining few elevations. Comparison of the multibeam bathymetric data collected in 2012 with data collected in 2022 further allowed for a formal geomorphic change detection analysis to map, quantify, and infer causation of morphological change over time with respect to a level of detectable change. The volume change in Randy Poynter Lake for the decade between 2012 and 2022 was slightly net-depositional and within the estimated uncertainty. The spatial distribution of sediment deposition was primarily concentrated in the northern portion of the lake, where the principal tributary flows into Randy Poynter Lake. The results of the geomorphic change analysis were used to further understand the future implications to storage capacity change. Despite the challenges of confirming systematic biases because of uncertainties exceeding the observed changes, insights from the study help predict long-term reservoir sediment accumulation, indicating a reservoir half-life extending about 650 years from 2022 on the basis of the current sediment yield estimates.

Georgia↗
USGS2024

Framework for mapping liquefaction hazard–Targeted design ground motions

Liquefaction-induced ground failure poses substantial challenges to geotechnical earthquake engineering design. Current approaches for designing against liquefaction hazards, as specified in most seismic provisions, focus on estimating a liquefaction factor of safety ( 𝐹⁢𝑆𝐿 ) and typically characterize earthquake loading using design parameters based on probabilistic or deterministic ground motion levels. Because 𝐹⁢𝑆𝐿 is estimated deterministically, this basis of design neglects considerable uncertainties for estimating liquefaction triggering and its consequences and results in a lack of liquefaction-specific design criteria, particularly as structural design has advanced toward risk-targeted performance objectives. This study presents a framework for developing liquefaction-targeted design criteria based on a minimum acceptable return period of liquefaction, informed by probabilistic liquefaction hazard analysis (PLHA). PLHA quantifies annualized rates of liquefaction by considering contributions from (1) the full ground-motion probability space, and (2) uncertainties in liquefaction triggering using probabilistic models. PLHA is used in this study to characterize the current, effective return periods of 𝐹⁢𝑆𝐿 ( 𝑇𝑅,𝐹⁢𝑆 ) obtained from conventional liquefaction hazard analysis (CLHA) using uniform-hazard ground motions. 𝑇𝑅,𝐹⁢𝑆 is evaluated in a parametric study of nearly 100 sites throughout the conterminous United States. The results indicate large geographic variations in acceptable liquefaction hazard levels, with implied 𝑇𝑅,𝐹⁢𝑆 ranging between approximately 1,000 to 3,000 years. To address these inconsistencies without the computational demands of full PLHA, a framework is proposed for developing a liquefaction-targeted design peak ground acceleration, 𝑃⁢𝐺⁢𝐴𝐿 , for use in liquefaction models that result in consistent liquefaction design levels across all geographic locations. The mapped 𝑃⁢𝐺⁢𝐴𝐿 is shown to be somewhat sensitive to site-specific properties, and adjustment factors are developed and presented. The proposed 𝑃⁢𝐺⁢𝐴𝐿 mapping procedure produces 𝐹⁢𝑆𝐿 estimates that are consistent with those obtained from full PLHA at a target 𝑇𝑅,𝐹⁢𝑆 , providing a promising roadmap to incorporating PLHA concepts into current liquefaction design methods.

Journal of Geotechnical and Geoenvironmental Engin↗
USGS2025

The GorDAS Distributed Acoustic Sensing experiment above the Cascadia locked zone and subducted Gorda Slab

The southernmost portion of the Cascadia Subduction zone in Northern California produces high rates of moderate and large earthquakes owing to subduction of the Gorda slab and deformation associated with the Mendocino Triple Junction. Distributed Acoustic Sensing (DAS) is rapidly advancing as a method for detecting earthquakes and imaging crustal structure. We have begun a long-term DAS monitoring experiment on buried telecom fiber in Arcata, California, with the goal of increasing the available recordings of moderate to large earthquakes as well as imaging seismogenic structures. We have recorded over a year's worth of data, including most aftershocks of the 2022 M w 6.4 Ferndale earthquake, though not the mainshock itself. The dataset includes numerous magnitude 3.5 and larger earthquakes including the 2023/01/01 M w 5.4 Rio Dell earthquake. Here we present initial results comparing an earthquake detection algorithm, run in real-time on the processing unit of the interrogator system, with both the ShakeAlert earthquake early warning system as well as a post-processed earthquake catalog developed with deep-learning phase-picker algorithms. The rapid onboard processing of the detector demonstrates the potential utility of DAS-based edge computing for earthquake early warning. We also verify the quality of the strain waveforms both in terms of peak amplitudes and waveform similarity using about five months of nodal seismometer data. These instruments were deployed roughly every 300 m along the ~15km long cable and validate large variations in peak strain over short distances that are seen in the DAS data. All data from time windows surrounding both the local and teleseismic earthquakes are publicly available, which will improve our understanding of both the performance of DAS systems in moderate earthquakes and earthquake hazards associated with the Gorda subduction zone.

California↗
USGS2026

Suitability of ground-motion models for seismic hazard assessment in Puerto Rico and the U.S. Virgin Islands

We perform linear mixed-effects analyses with a ground-motion dataset to evaluate how well ground-motion models (GMMs) fit active crustal, subduction interface, and subduction intraslab earthquakes in Puerto Rico and the US Virgin Islands. Most of the GMMs reproduce the spatial variation in peak ground motions with earthquake magnitude and rupture distance but predict ground motions 0.3–1.0 natural log units (35%–270%) greater than observations. Two GMMs developed for Puerto Rico that are based on ground-motion records from mostly magnitude 4–5 earthquakes do not perform as well. We attribute the period-dependent overprediction in observed ground motions to differences between the observed site response and the linear site response in the GMMs. Consequently, we developed region-specific GMMs by adjusting the period-dependent linear site response coefficients and period-dependent constant coefficients to remove most of the bias between predicted ground motions and observations. For the analysis, we compile ground-motion records and process waveforms to build a dataset with 10,127 records at 72 stations from 849 magnitude 4.0–6.4 earthquakes between 1 January 2006 and 31 March 2024. The earthquakes include active crustal, subduction interface, and subduction intraslab events. We evaluate the GMMs using the time-averaged shear wave speed in the top 30 m ( ), which we compile from site surveys and proxy values computed from horizontal to vertical spectral ratios. Site terms exhibit strong consistency across GMMs and crustal and subduction earthquakes, indicating that the linear mixed-effects analysis successfully isolates the effects of local site response. The event terms show little spatial correlation and more substantial variability than in other regions, which we attribute to uncertainties in the earthquake magnitudes. This analysis guides the selection of the GMMs for the 2025 update of the National Seismic Hazard Model for Puerto Rico and the US Virgin Islands.

Puerto Rico, U.S. Virgin Islands↗
USGS2026

Groundwater drought in the United States: Spatial and temporal variability

Many communities and ecosystems in the United States that are dependent on groundwater are potentially adversely affected by groundwater drought. We computed yearly groundwater-drought metrics and mean groundwater levels at well locations across the conterminous United States (CONUS), using data from wells and remotely sensed and modeled Gravity Recovery and Climate Experiment Drought Monitor Data Assimilation (GRACE-DADM). We also modeled the probability of low or high human impact at each well location. The spatial distribution of groundwater-drought duration and severity from 2001 to 2020 for 1,510 wells shows longer maximum duration and higher maximum severity events in drier regions like the Southwest than in wetter regions like the Northeast. Based on 613 wells in CONUS from 1981 to 2020, there are many significant decreases in drought duration and severity in the Northeast and many significant increases in annual-mean groundwater levels. In contrast, there are many significant increases in drought metrics and decreases in mean water levels in parts of the Southeast. There are major differences in trends from 2001 to 2020 between well-based and GRACE-DADM-based groundwater metrics in some CONUS regions and a very low correlation between trends at individual locations across CONUS. A potential reason for this disparity is the low GRACE-DADM resolution (∼12 km) and the potential for a large amount of groundwater variation at the local scale. Also, GRACE-DADM represents shallow, unconfined aquifers which may not match the screened interval of the monitoring wells we evaluated. Large spatial gaps in long-term, high frequency, and quality-assured groundwater-well monitoring data present a challenge for understanding groundwater-drought variability across CONUS. Remote sensing tools such as GRACE can help but cannot fully replace well monitoring, as highlighted by our study results. Substantially more long-term monitoring wells would more accurately represent groundwater-drought trends and spatial variability across CONUS, particularly in western regions.

conterminous United States↗
USGS2026

As above, so below? A framework for integrating long-term water quantity trends reveals divergent patterns in groundwater and low streamflow across the United States

Climate, land-use, and disturbance drive long-term global trends in groundwater levels and streamflow. At large scales, these trends are typically considered separately, despite the well-established concept that groundwater and surface water comprise a single resource. Joint trend assessment at national scales is challenging because it requires pairing and aggregating data from spatially disparate streamflow and groundwater monitoring sites for which no established framework exists. Here, we evaluate alternative approaches for integrating groundwater and streamflow data to enable joint trend analysis—a critical step toward understanding how water-budget components respond concurrently and interactively to environmental drivers. Mann–Kendall trends were computed for individual groundwater (annual mean depth) and streamflow (annual low of 7 d averages) sites across the U.S over 21- (2000–2020), 31- (1990–2020), and 41-year (1980–2020) periods. Regional Kendall trends were calculated using five spatially contiguous and noncontiguous regional classifications for aggregation based on subsurface (e.g. aquifer, geology) and surface (e.g. watershed, landscape) characteristics. Site-level results revealed contrasting trends, with tendencies toward increasing low flows (wetting) and increasing groundwater depths (drying). Agreement between streamflow and groundwater trends increased with regional aggregation and longer timeframes, though persistent skew toward streamflow wetting and groundwater drying remained. Results varied by region and trend period, with notable consistencies: unified drying in the West/Southwest and wetting in the Upper Midwest. Directional mismatches in long-term trends were prominent in the High Plains and Mississippi Alluvial Plain, whereas near-term mismatches were most evident in the Northwest. Aggregation by hydrologic landscape regions (HLR) yielded the greatest agreement between groundwater and streamflow trends. These findings indicate that coupled responses may represent combined influences of climate, relief, and geology, as captured by HLR, more strongly than geography or geology alone. Integrated water availability assessments may benefit from a multi-characteristic classification framework to treat groundwater and surface water as a unified resource.

Environmental Research: Water↗
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

End-user needs for remote sensing wetlands of the Prairie Pothole Region of North America

The Prairie Pothole Region (PPR) of North America comprises globally important grassland and wetland ecosystems critical for numerous populations of migratory birds. Due to the importance of this region for migratory birds, and particularly waterfowl, and the threats of habitat loss due to intensifying agriculture, there is a mature and diverse system of conservation organizations, agencies, and partnerships that spends hundreds of millions of dollars annually on habitat conservation to support migratory bird populations. Remote sensing can be a powerful tool for observing and evaluating global change at large scales as well as expanding inferences from field studies to the broader landscape with statistical models. However, development and utilization of these tools has lagged behind their demand for several reasons, including concerns over spatial and temporal resolution and accuracy of products; perception of a misalignment with decision-maker needs; technological barriers such as skill sets of conservation professionals, computing resources, data access, and usability. In this report, we summarize the needs of conservation professionals and scientists who use or want to use remote sensing data products to inform science about wetland change and conservation of wetlands in the PPR. We assembled this information through several methods leading up to, during, and following a January 2026 PPR Wetland Remote Sensing Workshop. The workshop included United States and Canadian scientists, conservation professionals, and policy experts. Our goal was to bring together end-users and remote sensing product developers jointly to explore reducing the lag between product development and utilization of products to inform science and conservation. Specifically, we aimed to identify gaps in wetland remote sensing that limit effective monitoring, management, and conservation in the PPR, and to develop a framework that outlines pathways to address these gaps by fostering collaboration, improving communication networks, encouraging discussion, and building on existing and ongoing efforts. This report summarizes our participants’ descriptions of end-user needs and the outcomes of the workshop.

Prairie Pothole region↗
USGS2025

Near-surface material and topography generate anomalous high-frequency ground motion amplification in Chugiak, Alaska

An ∼3 km long nodal array oriented approximately east–west was deployed in Chugiak, Alaska, by the U.S. Geological Survey during 2021. The array intersects with the permanent NetQuakes station NP.ARTY, where peak ground acceleration (PGA) value of 1.98 g was recorded during the 2018 M w 7.1 Anchorage, Alaska, earthquake, in sharp contrast to the PGA of ∼0.3 g at a site just 4 km to the west. Seismic data for M w 1.8–4.3 aftershocks from the M w 7.1 event recorded by the nodal array confirm the anomalously large ground motions obtained at NP.ARTY as well as similar amplifications at nodes within ∼1 km to the east. Here, we performed 0–10 Hz 3D finite‐difference simulations, including high‐resolution surface topography, to explore the cause of the unexpectedly large amplification. As expected, the simulations computed with a regional 3D tomography velocity model severely underpredict the 0–10 Hz acceleration records at almost all sites. Adding a near‐surface low‐velocity taper to 300 m depth amplifies the accelerations by up to a factor of 5 and enables a reasonable match between the nodal data and simulations at sites to the west of NP.ARTY. However, this model still underpredicts the spectral energy in the area covered by glacial sediments by up to an order of magnitude. The addition of a till layer using a depth‐dependent shear‐wave velocity ( ⁠⁠ V s ) profile along with a homogeneous, 8 m thick low‐velocity layer with V s = 250 m/s representing the kame terraces improves the fit to data to within a factor of 2 at nodes located on top of the glacial sediments. Our study shows that the anomalously large high‐frequency amplification recorded at and near NP.ARTY can be explained by a combination of topographic effects and near‐surface low‐velocity material with amplification effects on the high‐frequency ground motion by up to about 40% and an order of magnitude, respectively.

Alaska↗
USGS2025

Methods for estimating selected low-flow statistics at gaged and ungaged stream sites in Massachusetts

The U.S. Geological Survey, in cooperation with the Massachusetts Department of Conservation and Recreation, Office of Water Resources, computed selected at-site streamflow statistics at U.S. Geological Survey streamgages in and near Massachusetts and developed regional regression equations for estimating selected streamflows at ungaged stream sites in Massachusetts. Two sets of regional regression equations were developed: (1) the “mainland” equations, for mainland Massachusetts excluding the area covered by the second set, and (2) the “southeastern” equations, for the Plymouth-Carver-Kingston-Duxbury aquifer area in southeastern Massachusetts and for Cape Cod. The regression equations and at-site statistics may be used by Federal, State, and local water managers in addressing water-resources issues relevant in Massachusetts. Regional regression analyses for the mainland equations were developed to estimate the following 27 streamflow statistics: 99-, 98-, 95-, 90-, 85-, 80-, 75-, 70-, 60-, and 50-percent flow durations; monthly June, July, August, and September 90- and 50-percent flow durations; February, June, and August median of the monthly means; harmonic mean; and medians of the following annual low-flow frequency statistics: 7-day; 7-day, 2-year; 7-day, 10-year; 30-day, 2-year; and 30-day, 10-year. The analyses used 81 streamgages with minimal to no regulations in and near Massachusetts. The regression analyses determined that four basin characteristics—drainage area, combined hydrologic soils A and B, streamflow variability index, and annual mean temperature—were the only significant explanatory variables for the different mainland equations. Regional regression equations were developed for the Plymouth-Carver-Kingston-Duxbury aquifer area in southeastern Massachusetts and Cape Cod, because surface-water drainage areas and groundwater contributing areas do not always coincide in this area of the State. The regression analyses to estimate 10 flow durations from the 99th to 50th percentiles used 18 streamflow sites with some occasional minor regulations—because there are few unregulated streams in southeastern Massachusetts. The analyses determined that groundwater contributing area and storage (combined water bodies and wetlands) were the only significant explanatory variables in the southeastern equations.

Massachusetts↗
USGS2000

Quality assurance and analysis of water levels in wells on Pahute Mesa and vicinity, Nevada Test Site, Nye County, Nevada

Periodic and continual water-level data from 1963 to 1998 were compiled and quality assured for 65 observation wells on Pahute Mesa and vicinity, Nye County, Nevada. As part of the quality assurance of all water levels, ancillary data pertinent to computing hydraulic heads in wells were compiled and analyzed. Quality-assured water levels that were not necessarily in error but which did not represent static heads in the regional aquifer system, or required some other qualification, were flagged. Water levels flagged include those recovering from recent pumping or well construction, water levels affected by nuclear tests, and measurements affected by borehole deviations. A cursory examination of about 30 wells with available water-level and down-hole temperature data indicate that water levels in most wells on Pahute Mesa would not be significantly affected by temperature if corrected to 95 degrees Fahrenheit. Wells with large corrections (greater than 10 feet) are those with long water columns (greater than 1,500 feet of water above the assumed point of inflow) in combination with mean water-column temperatures exceeding 105 degrees Fahrenheit. Water-level fluctuations in wells on Pahute Mesa are caused by several factors including infiltration of precipitation, barometric pressure, Earth tides, ground-water pumpage, and seismic events caused by tectonic activity and underground nuclear testing. No observed water-level fluctuations were attributed to a naturally occurring earthquake. The magnitude and duration of changes in water levels caused by nuclear tests are affected by the test size and the distance from a well to the test. Identifying water levels that might be affected by past nuclear tests is difficult because pre-testing water-level data are sparse. Hydrologically significant trends were found in 13 of 25 wells with multiple years of water-level record. The largest change in water levels (1,029 feet in 25 years) occurred in well U-19v PS 1D as a result of the Almendro nuclear test. Likely explanations for trends in most of the wells are either changes in precipitation patterns that affect recharge rates to the ground-water system, pumping effects from water-supply well U-20 WW, or a combination of these two factors.

Nevada↗
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

Bathymetric and velocimetric surveys at highway bridges crossing the Missouri River near Kansas City, Missouri, August 8–9, 2023

Bathymetric and velocimetric data were collected by the U.S. Geological Survey, in cooperation with the Missouri Department of Transportation, near 8 bridge crossings of the Missouri River near Kansas City, Missouri, on August 8–9, 2023. A multibeam echosounder mapping system was used to obtain channel-bed elevations for river reaches that extended about 1,550 to 1,640 feet longitudinally and generally extended laterally across the active channel from bank to bank during low floodflow to nonflood conditions. These surveys provided the channel geometry and hydraulic conditions of the river at the time of the surveys and provided characteristics of scour holes, which may be useful in developing or verifying predictive guidelines or equations for computing potential scour depth. The data collected from the surveys may also be useful to the Missouri Department of Transportation as a record of low floodflow conditions in regards to the stability and integrity of the bridges with respect to bridge scour. Bathymetric data were collected around every in-channel pier. Scour holes were at most piers where bathymetry could be obtained, except for those piers on banks or surrounded by riprap. All the bridge sites in this study were surveyed and documented in previous studies. The average difference between the bathymetric surfaces ranged from 0.07 to 4.16 feet higher in 2023 than 2019, which indicates overall deposition between the survey dates, as might be expected based purely on streamflow at the time of the survey. However, the average difference between the bathymetric surfaces ranged from 1.44 feet higher to 1.88 feet lower in 2023 than 2015, which indicates a dynamic equilibrium of scour and deposition overall between those surveys, despite the lower flow conditions in 2023. Similarly, the average difference between the bathymetric surfaces ranged from 3.18 feet higher to 5.19 feet lower in 2023 than 2011, which indicates a relative equilibrium between scour and deposition overall, albeit the trend was toward scour as might be expected because of the substantial flood event in 2011. Riprap blankets and alignment to flow had a substantial effect on the size of the scour hole for a given pier. Piers that were partially or fully surrounded by riprap blankets had scour holes that were substantially smaller (to nonexistent) compared to piers with no rock or riprap and effectively mitigated the scour holes historically observed at these piers. Several of the structures had piers that were skewed to primary approach flow. At most of the structures, the scour hole was deeper and longer on the side of the pier with impinging flow than the leeward side, with some amount of deposition on the leeward side, as typically observed at piers skewed to approach flow.

Kansas, Missouri↗
USGS2026

Conventional hydrothermal power-producing systems of the Great Basin, USA

As part of the update to the electric-grade conventional hydrothermal assessment of the Great Basin, USA, Monte Carlo analyses of identified resources within explored regions will be performed to make estimates of discovered resources and associated uncertainty. Analyses use conditional statistics where estimates are conditioned upon a hydrothermal favorability map, allowing for the likelihood that more resources exist in regions of higher hydrothermal favorability. For these analyses, a dataset of identified hydrothermal systems is compiled, and the new compilation is described herein. Recognizing that a single hydrothermal system may be developed with multiple power plants, and that the hydrothermal upflow zone may be several kilometers across with many measurements characterizing a single hydrothermal system, a procedure was developed and employed to create clusters of points (power plants, measurements, etc.) that are associated with a single system, and a new central point was defined as the best estimator of the center of the hydrothermal system. Hydrothermal systems were uniquely identified by grouping electric-grade hydrothermal measurements and operating power plants within a distance of 10 km. Groups that are >10 km apart are assumed to be different electric-grade hydrothermal systems. While 10 km was used as the threshold, most systems were significantly further apart, and most points within groups were typically within 5 km of each other. A well measurement was considered an electric-grade measurement of a hydrothermal system if it had two properties: a measured temperature of >85 °C and evidence of hydrothermal convection. Other points that were added to the dataset are locations of operating powerplants or locations that have been classified as an electric-grade hydrothermal resource by either the U.S Geological Survey (USGS) or the Great Basin Center for Geothermal Energy. After all points are assigned to systems, new points were computed with the goal of identifying the center of the throat of the hydrothermal upflow zone. If operating powerplants exist for a system, then the arithmetic average of all power plant locations is used. Otherwise, if USGS made an estimate, that location is used. In the absence of both powerplants or USGS estimates, the arithmetic average of all electric-grade measurement locations is used. An example is shown of how these newly compiled locations might be ranked for uncertainty analyses, where higher confidence is assumed if measured temperature is higher and there are many supporting measurements indicating an electric-grade resource. In summary, 28 systems have operating power plants, an additional 78 systems are known identified electric-grade hydrothermal resources, and 99 new systems were identified as probable electric-grade systems with varying levels of confidence. These 205 locations are shown as a function of a recent hydrothermal favorability map, conceptually illustrating the conditional statistics that can be used to make estimates of the undiscovered resources of the Great Basin. An accompanying data release provides summaries of developed capacity by system and USGS estimates of likely total capacity and associated uncertainty.

Arizona, California, Idaho, Nevada, Oregon, Utah↗
USGS2026

Multiple-well monitoring site adjacent to the Midway- Sunset and Buena Vista Oil Fields, Kern County, California

Groundwater quality in and around oil fields in the Southern San Joaquin Valley is of interest to many California residents that rely heavily on groundwater for domestic, commercial, and agricultural use. To help assess the effects of historical oil-field activities and natural geologic sources on groundwater near the southwest margins of the Kern County Groundwater Subbasin, a multiple-well monitoring site was installed near the administrative boundary between the Midway-Sunset and Buena Vista Oil Fields in Kern County, California. The installation of the Midway-Sunset Buena Vista multiple-well monitoring site (MSBV) supports regional analysis of the relations of oil and gas sources to groundwater quality by providing information about the geology, hydrology, geophysical properties, and water quality of the alluvial and upper Tulare aquifers in areas where groundwater data were limited. Data collected from the site included drill cuttings, whole core samples, sidewall core samples, mud-gas analysis, borehole geophysical logs, depth to water measurements, and water quality samples. Whole cores were scanned using dual energy computed tomography. Subsamples of selected cores were analyzed for density, porosity, specific retention, and bulk minerology. Thin sections of the subsamples were prepared, photographed, and examined. Two samples were analyzed using scanning electron microscope technology to examine the microporosity of diatomite laden sediment. Instrumentation installed in the wells collect hourly depth to water measurements. Analysis of the data show there is 355 feet of alluvium overlying the Tulare Formation at the well site. The contact between the two formations is an aquitard resulting in a perched aquifer in the alluvium and unconfined aquifer in the Tulare Formation. The alluvium is more heterogenous and finer grained than the Tulare Formation resulting in markedly higher porosity in the alluvium compared to the Tulare Formation. Higher specific retention observed in the alluvium is attributed to the finer grained sediment and greater abundance of reworked diatomite (as represented by opal-CT [cristobalite-tridymite]) compared to the Tulare Formation. Total dissolved solids (TDS) approached or exceeded 10,000 milligrams per liter (mg/L) in the alluvium from approximately 176 to 242 feet below land surface and at the top of the Amnicola clay at approximately 670 feet below land surface within the Tulare Formation. Elevated TDS, chloride, and boron concentrations in the alluvium and on top of the Amnicola clay likely reflect groundwater that is mixed with oil-field water. Water chemistry and modern-aged groundwater in the alluvial monitoring well (MSBV #3) are consistent with the oil-field water in the alluvium being derived from documented historical surface disposal of oil-field water upslope (northwest) of the site. Water chemistry and pre-modern groundwater age in the deeper Tulare monitoring well (MSBV #1) on top of the Amnicola clay are consistent with oil-field fluids derived from upslope natural geologic sources or old oil wells that leak in the subsurface. Shallow groundwater in the Tulare (MSBV #2) is not affected by mixing with oil-field sources.

EarthArXiv↗