USGS ScienceSearch

SEARCH · USGS Science

Results for “Journal of Water Process Engineering”

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,711 records · Page 23Linked to original sources

Groundwater structures fish growth and production across a riverscape

Landscapes are composed of habitat patches and conditions that vary across space and time. While habitat variability and complexity can support important ecological processes and ecosystem services, the dynamic nature of habitats can also constrain organismal growth and production as optimal conditions are fleeting. In riverine ecosystems, groundwater discharge to streams stabilises water temperature and flow regimes, thus mediating how habitat complexity is expressed. Yet, how stable habitats structure growth and production within the broader landscape matrix is not well understood. In this study, we explored the effects of groundwater on spatiotemporal variation in growth and production for juvenile Yellowstone cutthroat trout ( Oncorhynchus virginalis bouvieri ) across the upper Snake River catchment, Wyoming, USA. We combined machine learning techniques and remotely sensed landscape data to estimate groundwater availability across the river network, which we linked to stream temperature regimes and conspecific density. We then used Bayesian hierarchical models to quantify the effects of temperature, density and groundwater on spatiotemporal variation in fish growth and production in 52 focal reaches. Finally, we predicted body size trajectories and trends in total production continuously over both space and time to understand the effect of groundwater at the riverscape scale. Groundwater discharged to streams where topography changes abruptly in valley-bottom areas underlain by coarse glacial deposits. Groundwater stabilised temperature regimes and was associated with high trout densities. Temperature and density, in turn, interacted to influence growth rates: growth increased strongly with temperature, but this effect was reduced when density was high. Accordingly, variation in groundwater availability among stream reaches diversified growth and production regimes. In reaches with low groundwater availability, growth and production declined over time from summer maxima. In contrast, in reaches with high groundwater availability, temporal trends in growth and production were hump-shaped—peaking in autumn—and mean production was greater. At the riverscape scale, temporal asynchrony in growth rates generated convergent spatial variation in growth capacity, but—when combined with density—led to the formation of distinct hotspots of production. Our results demonstrate how groundwater, an important driver of aquatic ecosystem heterogeneity, structures trout growth and production across space and time. Importantly, rare, but stable habitats may disproportionately affect ecological processes and serve as key sources of population diversity at larger spatial scales.

Wyoming

Forecasting volcanic activity in Germany—A multi-criteria approach

Igneous activity, including shallow intrusions and volcanism, has the potential to disrupt underground critical infrastructure. Notably, future underground infrastructure projects like high-level radioactive waste repositories must be sited in areas of extremely low disruption probability by igneous activity. In Germany, according to the Repository Site Selection Act of 2017 (Standortauswahlgesetz, or StandAG), areas in which Quaternary volcanism is either present or future volcanic activity is expected within the next 1 million years (m.y.) must be excluded from the site selection process. Although the locations of regions with Quaternary volcanism are reasonably well known in Germany, forecasting potential igneous activity at intraplate volcanic fields is challenging, as many processes and their interactions control the spatial distribution of volcanic centers. Here, a semi-quantitative, multi-criteria approach is proposed for a regional evaluation of the relative potential of future igneous activity in Germany. A variety of geoscientific indicators are used, including seismic anomalies in Earth’s mantle, gravity data, tectonic activity, sutures, ground motion, earthquakes, mantle degassing centers, and geochronological data of volcanic rocks. The indicators describe the sequence of processes from potential melt generation in Earth’s mantle, through ascent and accumulation of melt within the lithosphere, to eruption at Earth’s surface. In total, 15 out of 30 proposed geoscientific indicators are selected and quantified using 20 total assigned parameters. Defined threshold values are used to spatially delimit relevant parameter properties to focus on areas with higher potential of future magmatic activity. To consider uncertainties of parameters and their underlying processes, which are usually more spatially extensive below ground, buffer zones are defined in which values of relevance decrease with increasing distance from the initial lateral shape of a parameter. Normalized parameters are combined into an index, whose spatial value distribution is used to differentiate the relative potential of future igneous activity (within the next 1 m.y.). The sensitivity of the results is shown by varying the weighting factors for the relevant parameters in country-wide index maps. Thereby, profiles illustrate the distribution of the resulting index values and respective index fractions of various parameters. Different index maps for the relative potential of future igneous activity are presented and can be used for hazard assessments.

Professional Paper

Evaluation of daily stream temperature predictions (1979-2021) across the contiguous United States using a spatiotemporal aware machine learning algorithm

Stream temperature controls a variety of physical and biological processes that affect ecosystems, human health, and economic activities. We used 42 years (1979–2021) of data to predict daily summary statistics of stream temperature across >50,000 stream reaches in the contiguous United States using a recurrent graph convolution network. We comprehensively documented the performance – both across all reaches and by stream type (e.g., reservoir or groundwater influence) – as a baseline for future improvement. The model showed reach-level RMSE of <2 °C with 90 % prediction intervals that contain 90.7 % of observations. We also assessed how the model captured variability in ecologically relevant metrics (e.g., R 2 for annual 7-day maximum = 0.76; R 2 for days exceeding 25 °C = 0.75). This model does not outperform state-of-the-art machine learning efforts (e.g., RMSE ≤1.5 °C) due to a limited input set but does provide the most spatially complete modeling to date to support water availability assessments.

contiguous United States

Global geologic map of Europa

Discovered by Galileo Galilei more than 400 years ago and imaged in detail by the Voyager 2 Galileo spacecraft, Jupiter’s icy moon Europa has been a source of intrigue. A range of science investigations indicate that it contains the key ingredients for habitability, notably energy, chemistry, and liquid water. Europa’s surface is geologically complex and, based on the dearth of impact craters, interpreted to be as young as ~60 Ma. The array of geologic features that characterize the surface include extensive ridged plains, regions of broad disruption termed chaotic terrain, long, quasi-linear ridges that span thousands of kilometers, and bands as much as 60 kilometers wide and that extend 100s of kilometers. These features, along with other geophysical measurements, indicate the presence of a global briny liquid water ocean beneath the ice shell. It was not until the arrival of the Galileo spacecraft in 1995 that the true nature and level of complexity of the surface was revealed. Although image data returned by Galileo provided insight into the structure of a variety of regions, the entire satellite has yet to be observed at a regional scale (less than 250 meters per pixel) and the detailed geologic nature of much of its surface remains a mystery. Establishing the global context of the distribution and timing of Europan geologic units forms a basis to understand regional and local scale processes, serves as a tool for the planning of future missions, and most of all is essential to gaining insight into the potential habitability of this icy world.

Scientific Investigations Map

Preliminary depth to basement modeling at Salton Sea, California

The San Andreas Fault – Imperial Fault (SAF-IF) transtensional step-over zone along the southern margin of the Salton Sea hosts substantial geothermal production and lithium brine resources. Recent volcanism at the Salton Buttes and active seismicity along the SAFIF fault system highlight active tectonic and magmatic processes that pose natural hazards and may impact energy and mineral production. Characterizing the subsurface architecture and extent of concealed alteration associated with this tectono-magmatic system enhances understanding of these active processes, associated hazards, and resources. We have compiled a gravity database, consisting of new and re-processed existing data, from which we have constructed a new isostatic residual gravity anomaly map of the Salton trough. We have used this new gravity dataset together with a compilation of publicly available borehole data to develop new depth to basement inversion models for the region. These depth to basement models help to constrain basin geometries, inform alteration mapping, and reveal variations in basement rocks. Due to the concealed nature of the complex tectonic framework at the Salton trough, it is necessary to utilize geophysical methods for subsurface characterization. These new depth to basement models are a first step toward constructing 2D and 3D geophysical and geologic models of the Imperial Valley and Salton Sea geothermal area. This analysis complements other geophysical initiatives, including magnetotelluric (MT) modeling (Tokmakoff et al., 2024), magnetic mapping (Glen and Earney, 2023, 2024) and potential field modeling, and seismic studies focused on hazard and resource investigations in the Imperial Valley.

California

A new groundwater energy transport model for the MODFLOW hydrologic simulator

Heat transport in the subsurface is an important aspect of research related to the effects of a warming climate on ecological services (i.e., cold-water refugia); the development of geothermal resources for energy banking schemes (i.e., aquifer thermal energy storage [ATES]); and the effects of temperature on other aspects of groundwater quality, such as nutrient cycling. Historically, simulation of heat transport using the MODFLOW groundwater simulator and related codes was performed by scaling the input parameters of a solute-transport model to emulate heat transport. However, that approach required additional pre- and post-processing of input and output and could not account for the variation in effective thermal storage and transport properties during transient, unsaturated flow, for example. True heat-transport capabilities in the context of MODFLOW were first introduced in a variant called USG-Transport. More recently, a new groundwater energy-transport (GWE) model type has been added to MODFLOW 6, the core version of the MODFLOW hydrologic simulator. GWE supports the simulation of heat transport on structured or unstructured grids as well as within and between features of advanced packages that represent streams, lakes, multi-aquifer wells, and the unsaturated zone. GWE is integrated within MODFLOW 6 and is accessible through the FloPy Python package and the MODFLOW 6 application programming interface (API). An example simulation demonstrates conduction between grid cells through both the water and the solid aquifer material, including thermal bleeding from saturated overburden cells into a groundwater flow field.

Groundwater

Analyzing spatial distributions and alignments of pitted cone features in Utopia Planitia on Mars

Martian geomorphology and surface features provide links to understanding past geologic processes such as fluid movement, local and regional tectonics, and feature formation mechanisms. Pitted cones are common features in the northern plains basins of Mars. They have been proposed to have formed from upwelling volatile-rich fluids, such as magma or water-sediment slurries. In this study, we map the spatial distributions of pitted cone features across the Utopia Planitia (UP) basin. Using the Average Nearest Neighbor technique, we find that pitted cone features appear to be clustered across the basin, occurring in a narrow band around the circumferential rim of UP. Additionally, we find that pitted cone features also appear to be aligned in chains that are sub-parallel with the basin margin. Parallel bands of cones generally follow elevation contours of the UP basin, which suggests elevation or a correlated factor plays a major role in pitted cone and cone chain formation. We propose that pitted cones and cone chains may be related to vertical fractures formed around the UP basin rim from subsidence of infilling basin material.

Icarus

Geophysical identification of potential groundwater – surface water interactions in the Bonita Peak Mining District, San Juan County, Colorado

This report includes the results and interpretation of geophysical investigations to aid in the location, identification, and conceptual site model (CSM) development of the near surface groundwater and surface water hydrology of the Upper Cement Creek and California Gulch of the Animas Rivers within the Bonita Peak Mining Districts (BPMD). Mine influenced water (MIW) from numerous mine adits within the BPMD has altered the geochemistry and ecology of these ground and surface waters. Many studies of varying scales have investigated the geology, hydrogeology, geochemistry, and aqueous chemistry of these adits and impacted areas. Field investigations using non-invasive geophysical methods were conducted in the summer of 2019 to aid in identifying potential groundwater inputs to these streams and further the CSM of the near surface within and along these streams. Two field deployments were completed in August and September of 2019. In August 2019, Phase 1 was completed using multifrequency electromagnetic induction (EMI) and thermal infrared imaging (TIR) to measure the bulk electrical conductivity, magnetic susceptibility, and the temperature response at known and suspected groundwater seeps into the Upper Cement Creek and California Gulch stretches of these alpine streams. Results and interpretation from Phase 1 were used to co-locate stretches of these streams for further detailed fiber optic distributed temperature system (FODTS) monitoring of the temperature dynamics for submeter identification of potential groundwater seeps during Phase 2. Phase 2 included additional EMI measurements, FODTS measurements, and repeated some Phase 1 measurements for improved data quality. All these data were used for the combined qualitative interpretation identifying seep locations. The data from both Phase 1 and Phase 2 are available in the USGS data release at: https://www.sciencebase.gov/catalog/item/5f18a30082cef313ed84879b. These data include detailed near-surface electrical conductivity and magnetic susceptibility mapping, discrete TIR imaging, and several kilometers of FODTS data. These data were georeferenced and used to identify potential groundwater seeps into these streams. The results indicate a complex near surface geophysical and hydrogeologic regime. Seeps are classified and spatially correlated to the EMI, TIR, and geochemical sampling results. This report is a thorough compilation of the geophysical field efforts in the summer of 2019, the data processing, and interpretations of these data. These data alone should not be used for final interpretations of the complex hydrogeology within these alpine systems. Rather, these data should be combined with other known data from various investigations to further refine the CSM and guide decisions for future investigations, watershed protective measures, and other site decisions. Geographical information system (GIS) layers of these results are available to project managers for such integration and combined interpretations.

Colorado

Patterns and drivers of cliff erosion in Big Sur, California, USA using repeat photogrammetry, 2017–2023

Seacliff erosion in steep terrain poses major risks to transportation and critical infrastructure. In Big Sur, California, USA, seacliff erosion threatens the sustainability of the central coast stretch of California State Route 1, a transportation corridor that is critical to the region's economy. Published cliff retreat rates for the region range from 1 to 40 cm yr −1 , highlighting that high-resolution, process-based studies could enhance understanding of the causes of spatial and temporal variability. We quantified cliff erosion and investigated its drivers along ∼13 km of the Big Sur coastline at week–month timescales during the late fall to early spring wet seasons between January 2017 and June 2023 by analyzing 3D point clouds developed from aerial imagery using four-dimensional structure-from-motion (4D SfM) photogrammetry techniques. We calculated cliff face retreat rates of 2.23 ± 3.06 cm yr −1 (mean ±1 σ ), an order of magnitude lower than long-term estimated rates for the region (which included large deep-seated landslides), but in line with short-term rates reported across California. Change detection imagery comparison, cliff profiles through time, and statistical analysis reveal a cyclical cliff evolution process in which erosion by wave action at the cliff base destabilizes the cliff and primes it for subsequent failure during precipitation events. Although more erosion by volume could be attributed to precipitation-induced increases in soil moisture (784 m 3 km −1 yr −1 ) compared with erosion attributed to wave power (282 m 3 km −1 yr −1 ), our observations underscore the coupled nature of these processes in driving cliff evolution, consistent with established theory and observations.

California

Evaluation of water quality in the Langford Valley–Irwin Groundwater Subbasin, Fort Irwin National Training Center, California, 1993–2019

The U.S. Army Fort Irwin National Training Center (NTC), 120 miles northeast of Los Angeles in the Mojave Desert of California, obtains its potable water supply from the Bicycle Valley and Langford Valley groundwater basins; Langford Valley groundwater basin is further subdivided into the Langford Well Lake and Irwin groundwater subbasins. The Irwin groundwater subbasin contains younger, unconsolidated deposits that have a saturated thickness of as much as 200 feet (ft) and a lower aquifer within older unconsolidated deposits as thick as 650 ft. Groundwater recharge under predevelopment conditions (before 1941) occurred primarily from infiltration of intermittent streamflow in small washes that cross the Irwin groundwater subbasin. Since that time, groundwater recharge has increased because of growth of the NTC in recent years and as a result of other processes, including (1) infiltration of treated wastewater into the aquifer through ponds near the NTC wastewater treatment facility (WWTF) and (2) infiltration of imported water and treated wastewater used for landscape irrigation at base housing and athletic fields. Water samples were collected from 17 wells and analyzed for field parameters, chemical constituents, and isotope composition in the Irwin groundwater subbasin between 2014 and 2019. These data were supplemented with water-chemistry data collected during 1993–95 and at other times if available. Between 1993–95 and 2015–19, median dissolved solids and nitrate concentrations in water from wells in the Irwin groundwater subbasin increased from 620 to 1,030 milligrams per liter (mg/L) and from 2.8 to 4.5 mg/L as nitrogen, respectively. After 2014, dissolved solids and nitrate concentrations in water from wells near the NTC WWTF decreased as a result of changes in source water quality attributable to reverse osmosis of treated drinking water delivered within the Irwin groundwater subbasin and to increased levels of treatment at the NTC WWTF. Based on delta oxygen-18 and delta deuterium isotope data, increases in dissolved solids concentrations in water from most wells were consistent with evaporation prior to recharge and mobilization of soluble salts from the unsaturated zone. Arsenic and fluoride concentrations in water from wells decreased between 1993–95 and 2015–19 as the basin filled with treated wastewater, but 2015–19 concentrations generally exceeded the California State Water Resources Control Board maximum contaminant levels of 10 micrograms per liter for arsenic and 2 mg/L for fluoride. Most groundwater in the Irwin groundwater subbasin has unadjusted carbon-14 ages ranging from 18,400 to 12,350 years before present. However, water from well 10E3, located along the wash near the subbasin outflow in the southeastern part of the Irwin groundwater subbasin, contained measurable tritium, which is consistent with infiltration of intermittent streamflow and groundwater recharge from the wash after 1952. Chemical and isotopic data indicate that treated wastewater is present in water from most wells in the upper aquifer that underlies the Irwin groundwater subbasin. Wells were not sampled to adequately determine the extent of treated wastewater and changes in water quality within the lower aquifer that underlies the Irwin groundwater subbasin.

California

Anaerobic biodegradation of perfluorooctane sulfonate (PFOS) and microbial community composition in soil amended with a dechlorinating culture and chlorinated solvents

Perfluorooctane sulfonate (PFOS), one of the most frequently detected per - and polyfluoroalkyl substances (PFAS) occurring in soil, surface water, and groundwater near sites contaminated with aqueous film-forming foam (AFFF), has proven to be recalcitrant to many destructive remedies, including chemical oxidation. We investigated the potential to utilize microbially mediated reduction (bioreduction) to degrade PFOS and other PFAS through addition of a known dehalogenating culture, WBC-2, to soil obtained from an AFFF-contaminated site. A substantial decrease in total mass of PFOS (soil and water) was observed in microcosms amended with WBC-2 and chlorinated volatile organic compound (cVOC) co-contaminants — 46.4 ± 11.0 % removal of PFOS over the 45-day experiment. In contrast, perfluorooctanoate (PFOA) and 6:2 fluorotelomer sulfonate (6:2 FTS) concentrations did not decrease in the same microcosms. The low or non-detectable concentrations of potential metabolites in full PFAS analyses, including after application of the total oxidizable precursor assay, indicated that defluorination occurred to non-fluorinated compounds or ultrashort-chain PFAS. Nevertheless, additional research on the metabolites and degradation pathways is needed. Population abundances of known dehalorespirers did not change with PFOS removal during the experiment, making their association with PFOS removal unclear. An increased abundance of sulfate reducers in the genus Desulfosporosinus (Firmicutes) and Sulfurospirillum (Campilobacterota) was observed with PFOS removal, most likely linked to initiation of biodegradation by desulfonation. These results have important implications for development of in situ bioremediation methods for PFAS and advancing knowledge of natural attenuation processes.

Science of the Total Environment

Cascading land surface hazards as a nexus in the Earth system

Earth’s surface is sculpted by numerous processes that move sediment, ranging from gradual and benign to abrupt and catastrophic. Although infrequent, high-magnitude sediment mobilization events can be hazardous to people and infrastructure, leaving topographic imprints on the landscape and remarkable narratives in the historical record. Hazardous events such as fires, storms, and earthquakes accelerate erosion and sediment transport, increasing landscape sensitivity to subsequent perturbations, thus forming a cascading hazard. Although the redistribution of sediment across Earth’s landscape can result in higher risks to vulnerable populations, cascading processes are commonly unaccounted for in hazard assessments. Cascading hazards can occur almost immediately after triggering events, such as coseismic landslides, or over months, years, or even decades after an initial perturbation, such as debris flows after wildfires or flooding in channels alluviated by volcanic debris. Sediment cascades span Earth’s surface, from mountaintops to river valleys, where erosion, deposition, and aggradation can lead to a myriad of hazardous processes, including decreased river conveyance capacity, which increases the likelihood of downstream flooding. An improved understanding of the magnitude, frequency, and persistence of cascading hazards is critical given the rapid changes in the frequency and severity of storms, fires, sea-level change, and cryospheric melting, as well as the expansion of high-population-density urban footprints in regions susceptible to solid Earth hazards. Understanding the full consequences and underlying physics of Earth’s cascading land surface hazards can help minimize future human and economic losses.

California

An investigation of shallow ground-water quality near East Fork Poplar Creek, Oak Ridge, Tennessee

Alluvial soils and fill materials in and near the floodplain of East Fork Poplar Creek in Oak Ridge, Tennessee, are contaminated with various trace metals (primarily mercury), organic compounds, and radionuclides that were lost to the stream as a result of past operations at the Y-12 Plant, a nuclear-processing facility located within the U.S. Department of Energy’s Oak Ridge Reservation. Observation wells were installed in the shallow (above bedrock) aquifer in and near the floodplain of East Fork Poplar Creek and water-quality samples were collected to determine if contaminants found in the soils and fill are also present in the shallow ground water. The shallow aquifer in the East Fork Poplar Creek flood plain consists primarily of alluvial silt and clay with lesser amounts of sand and gravel. Thickness of the shallow aquifer ranges from essentially zero to as much as 20 feet. A silty-clay glei horizon is present between the base of the alluvium and the top of bedrock at most flood-plain locations and, where present, likely impedes downward ground-water movement. Water in the shallow aquifer near East Fork Poplar Creek occurs under water-table conditions. Recharge to the shallow aquifer is principally from precipitation, and discharge is through springs and seeps to East Fork Poplar Creek and its tributaries. During spring, summer; and fall, evapotranspiration also accounts for the removal of water in storage in the shallow aquifer. Water levels in the shallow aquifer fluctuate seasonally in response to variations in recharge and evapotranspiration. Generally, the depth to water in the observation wells ranged from about 1 to 4 feet below land surface in late winter, and from about 2 to 7 feet below land surface in late fall. During extremely dry periods, the water table recedes below the top of bedrock in some flood-plain areas, possibly causing East Fork Poplar Creek to lose water to the shallow aquifer along some reaches. Contaminants found in water samples collected from several of the observation wells in concentrations (total and (or) total-recoverable) which equaled or exceeded drinking-water standards established by the U.S. Environmental Protection Agency are antimony, chromium, lead, mercury, selenium, total phenols, and strontium-90. Total and dissolved uranium concentrations exceeded 1.0 microgram per liter in samples from nearly 70 percent of the wells in the East Fork Poplar Creek flood plain. Organic compounds that were identified in low concentrations in samples from a few wells in the flood plain are: Arochlor 1260, benzo(a)anthracene, benzo(b)fluoranthene, chrysene, 3,3-dichlorobenzidine, di-n-butylphthalate, N-nitrosodiphenylamine, and pyrene. Water from one well in the East Fork Poplar Creek flood plain at a contaminated fill site contained 37and 8 micrograms per liter of trichloroethene and trans-1,2-dichloroethene, respectively. Comparison of the results of total and (or) total-recoverable trace-metal determinations with those from dissolved determinations demonstrates that elevated concentrations of these substances in water collected from several of the wells in the East Fork Poplar Creek flood plain resulted from sorption of trace metals (and possibly organic compounds and radionuclides) by fine sediment suspended in the samples. The occurrence of contaminated sediment in these samples is suspected to be the result of borehole contamination during well installation.

Tennessee

lasertram: A Python library for time resolved analysis of laser ablation inductively coupled plasma mass spectrometry data

Laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS) data has a wide variety of uses in the geosciences for in-situ chemical analysis of complex natural materials. Improvements to instrument capabilities and operating software have drastically reduced the time required to generate large volumes of data relative to previous methodologies. Raw data from LA-ICP-MS, however, is in counts per unit time (typically counts per second), not elemental concentrations and converting these count ratesto concentrations requires additional processing. For complex materials where the ablated volume may contain a range of material compositions, a moderate amount of user input is also required if appropriate concentrations are to be accurately calculated. In geologic materials such as glasses and minerals that potentially have numerous heterogeneities (e.g., microlites or other inclusions) within them, this is typically determiningwhether the total ablation signal should be filtered to remove these heterogeneities. This necessitates that the LA-ICP-MS data processing pipeline is one that is not automated, but is also designed to enable rapid and efficient processing of large volumes of data. Here we introduce , a Python library for the time resolved analysis of LA-ICP-MS data. We outline its mathematical theory, code structure, and provide an example of how it can be used to provide the time resolved analysis necessitated by LA-ICP-MS data of complex geologic materials. Throughout the pipeline we show how metadata and data are incrementally added to the objects created such that virtually any aspect of an experiment may be interrogated and its quality assessed. We also show, that when combined with other Python libraries for building graphical user interfaces, it can be utilized outside of a pure scripting environment. can be found at https://doi.org/10.5066/P1DZUR3Z

Applied Computing and Geosciences

Teach me how to pycap: A high-capacity well decision support tool using analytical solutions in Python

Regulatory agencies in humid temperate environments rely on timely evaluations of streamflow depletion and drawdown to protect aquatic ecosystems and existing water users. Numerical models offer detailed insights, but their complexity and time demands often preclude their practical use in rapid decision-making. We present pycap-dss, an open-source Python package that implements a suite of analytical solutions for estimating streamflow depletion and drawdown. The tool supports superposition of multiple wells and time-varying pumping, enabling cumulative impact assessments in situations with multiple wells and streams. The software is modular and extensible, allowing users to interchange solutions or add new analytical methods. A YAML-based configuration supports batch processing of multiple wells, and an optional AnalysisProject class facilitates integration with regulatory workflows. Rigorous unit and regression testing ensures computational reliability, and continuous integration supports ongoing development. We demonstrate deterministic examples of drawdown where multiple solutions are readily compared and streamflow depletion with multiple wells in the Central Sands region of Wisconsin. We also show the value of Monte Carlo analyses of streamflow depletion in the same Central Sands example, leveraging computational efficiency to evaluate the uncertainty of individual and cumulative streamflow depletion calculations from over 200 high-capacity wells.

Wisconsin

Neogene hydrothermal Fe- and Mn-oxide mineralization of Paleozoic continental rocks, Amerasia Basin, Arctic Ocean

Rocks dredged from water depths of 1,605, 2,500, 3,300, and 3,400 m in the Arctic Ocean included Paleozoic continental rocks pervasively mineralized during the Neogene by hydrothermal Fe and Mn oxides. Samples were recovered in three dredge hauls from the Chukchi Borderland and one from Mendeleev Ridge north of Alaska and eastern Siberia, respectively. Many of the rocks were so pervasively altered that the protolith could not be identified, while others had volcanic, plutonic, and metamorphic protoliths. The mineralized rocks were cemented and partly to wholly replaced by the hydrothermal oxides. The Amerasia Basin, where the Chukchi Borderland and Mendeleev Ridge occur, supports a series of faults and fractures that serve as major zones of crustal weakness. We propose that the stratabound hydrothermal deposits formed through the flux of hydrothermal fluids along Paleozoic and Mesozoic faults related to block faulting along a rifted margin during minor episodes of Neogene tectonism and were later exposed at the seafloor through slumping or other gravity processes. Tectonically driven hydrothermal circulation most likely facilitated the pervasive mineralization along fault surfaces via frictional heating, hydrofracturing brecciation, and low- to moderate temperature Fe- and Mn-rich hydrothermal fluids, which mineralized the crushed, altered, and brecciated rocks.

Geochemistry, Geophysics, Geosystems

Persistent deep long-period seismicity near the Lassen Volcanic Center

Deep long-period (DLP) earthquakes have been observed at many volcanic settings around the world and linked to the magmatic processes that drive volcanic unrest. At the Lassen Volcanic Center (LVC) of the Cascade arc, limited detection of DLP activity hinders classification of anomalous seismic behavior and its relationship to the LVC magmatic system. This study uses a template matching approach with seismic data from temporary nodal and permanent stations to detect and locate DLP earthquakes near the LVC between 2017 and 2024. Within the DLP catalog of 611 events, a transition occurs from scattered (2017–2020) to oscillatory (2021–2024) occurrence rates. During the oscillatory period, regional earthquakes with large amplitude velocity waveforms observed near the LVC are associated with abrupt changes in DLP occurrence rates. Continued monitoring of DLP activity has the potential to better define the processes that drive volcanic unrest at the LVC in the future.

California

Three-dimensional temperature maps of the Williston Basin, USA: Implications for deep hot sedimentary and enhanced geothermal resources

As part of U.S. Geological Survey's (USGS) efforts to identify and assess geothermal energy resources of the US, a three-dimensional (3D) geologic and thermal model has been constructed for the Williston Basin, USA. The geologic model consists of all sedimentary units above the Proterozoic and Archean crystalline rock (called basement herein), with a total sedimentary thickness of up to 5 km near the basin center. Twenty-nine geologic units were mapped from interpreted formation tops from 16,465 wells. A 3D temperature model was constructed to a depth of 7 km by constructing a 3D heat flow model for the sedimentary units, followed by estimating underlying temperature using a one-dimensional (1D) analytic solution for heat flow within the underlying crystalline basement. Using the sedimentary basin model, heat flow was simulated in 3D and was calibrated using three temperature datasets: 1) 24 high-confidence static temperature logs (equilibrium thermal profiles), 2) more than15,000 drill stem test (DST) measurements from >7,000 wells, and 3) more than 45,000 bottomhole temperature (BHT) measurements from >14,000 wells. The DST and BHT datasets provide broad spatial coverage, but are lower confidence, primarily because measurements were made prior to attaining thermal equilibrium. DST and BHT measurements were binned regionally to develop representative thermal profiles that generally agree with these lower quality data (hereafter called pseudowell temperature profiles). Layer properties (primarily thermal conductivity and compaction curves) were set to best estimate values, then the heat flow model was calibrated to fit pseudowell and static temperature logs primarily by adjusting basal heat flow to approximate the overall temperature profile. Minor adjustments to thermal conductivity allowed adjusting changes in slope at lithologic contacts. Resulting maps include 3D temperature and basal (bottom of sedimentary units) heat flow estimates, which are used as input for the temperature model of the basement. The crystalline basement temperature model uses an analytic 1D solution to the heat flow equation that requires estimates of heat flow and temperature at the upper boundary (i.e., the sediment/basement contact), radiogenic heat production within the crystalline basement, and reference thermal conductivity (i.e., uncorrected for temperature). Two regions of high heat flow are identified: 1) in western North Dakota along the North American Central Plains Conductivity Anomaly and 2) in eastern Montana near the Poplar dome. Within the sedimentary column in the center of the basin of the basin, an area of approximately 100,000 km2 is predicted to have moderate- to high-temperature geothermal resources (>90 °C) under the thickest sequences of sediments. Where thick insulation and high heat flow coincide, electric-grade resources can be less than 4 km deep. Assuming a maximum feasible drilling depth of 7 km, temperatures are predicted to be as high as 175 °C. The geologic model may be used to identify strata at sufficient temperatures that may have natural permeability or that may have conditions that favor development of enhanced/engineered geothermal systems resources.

Montana, North Dakota, South Dakota