SEARCH · USGS Science

Results for “Transaction of the American Fisheries Society”

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

SourceResearch connectionFrom yearThrough year

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

1,688 records · Page 40Linked to original sources

USGS2025

Ultralong, supershear rupture of the 2025 Mw 7.7 Mandalay earthquake reveals unaccounted risk

The 28 March 2025 moment magnitude ( M w ) 7.7 earthquake in Mandalay, Burma (Myanmar), ruptured 475 kilometers of the Sagaing Fault, which was more than twice the length predicted by magnitude scaling relationships. Kinematic slip models and observation of a Rayleigh Mach wave that passed through parts of Thailand confirmed that rupture occurred at supershear velocities of greater than 5 kilometers per second. The anomalous length exposed a vast population to violent near-fault shaking. The Mandalay earthquake is a modern analog for the M w 7.9 1906 San Francisco earthquake, another atypically long and fast rupture. Probabilistic seismic hazard analyses use scaling relations that do not account for such long ruptures at moderate magnitudes. This limitation, in conjunction with a likely increased population and infrastructure exposure for atypically long ruptures, contributes to a potential mischaracterization of seismic risk.

Mandalay↗

USGS2026

Generating geochemical and mineralogy distributions of soil in the conterminous United States using Bayesian hierarchical spatial models

Characterizing geochemical and mineralogical soil distributions across large spatial extents is essential for understanding mineral resources, ecosystem processes, and environmental risks. Rasters of soil geochemical distributions for the conterminous United States, however, are limited. We present a Bayesian modeling workflow and tool for generating predictive geochemical and mineralogy distribution maps for the conterminous United States using integrated nested Laplace approximation (INLA) with the stochastic partial differential equation approach. By modeling soil geostatistical data with environmental covariates (soil properties, topography, climate, and land cover), we generate predictive distributions of soil geochemistry that can be mapped or extracted for further analyses. As an example, we model the spatial distribution of trace elements in soil relevant to vertebrate health (cobalt, copper, iron, manganese, selenium, and zinc) and provide a workflow that can be used to generate and visualize predictive distributions of 39 other major and trace elements and 21 minerals of the soil survey, supporting a variety of ecological, environmental, and agricultural applications.

MethodsX↗

USGS2026

Quantitative mineral resource assessment of lithium pegmatite deposits in the southern Appalachian orogen

The first quantitative mineral resource assessment for undiscovered lithium pegmatite deposits in the southern Appalachian region of the United States was conducted. Permissive tracts for lithium pegmatite deposits were delineated by integrating lithological, tectonic, geochemical, geophysical and mineral occurrence data. Lithium pegmatite prospectivity of the tracts was ranked with simplified mappable criteria, including proximity to Paleozoic felsic intrusions and major lithotectonic structures, stream sediment geochemical anomalies, and pegmatite occurrence data. The geospatial data and permissive tracts were used to estimate the number of undiscovered lithium pegmatite deposits. These estimates were integrated into probabilistic simulations along with a new global lithium pegmatite grade and tonnage dataset to quantify potential contained undiscovered lithium resources. An economic filter was applied to convert the probabilistic estimates of contained lithium into recoverable material. The identified lithium pegmatite resources for the Carolina Lithium and Kings Mountain deposits, North Carolina, contain 1589 thousand tons (kt) of Li 2 O. The median contained undiscovered resource for the southern Appalachian orogen was estimated to be 2240 kt Li 2 O. At 90% confidence, the region contains at least 130 kt Li 2 O, and 10,700 kt at 10% confidence. After applying economic filters, the median recoverable contained resource was 1430 kt Li 2 O, corresponding to approximately 201 years of current lithium imports for consumption in the United States. North and South Carolina are likely to contain most of these resources. Coarse data resolution and intra-state variations in the geological data contribute to uncertainty of undiscovered lithium pegmatite resources. Continued efforts to harmonize disparate geospatial datasets with updated or new information can improve the accuracy and precision of estimated undiscovered lithium pegmatite resources in the study area and at broader scales.

Alabama, Georgia, Maryland, North Carolina, South ↗

USGS2025

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↗

USGS2024

Methodology for inclusion of produced and stored carbon dioxide in the U.S. Geological Survey Federal lands greenhouse gas inventory

The U.S. Geological Survey (USGS) has developed two new carbon dioxide (CO2) emissions and sequestration accounting methods for use in future reports. The first method is a Federal lease-produced CO2 emissions calculation for an update of the report, “Federal Lands Greenhouse Gas Emissions and Sequestration in the United States.” The methodology to incorporate Federal lease CO2 production emissions into the updated report relies on CO2 sales royalty data from the Office of Natural Resources Revenue (ONRR). The end usage points for the gas include enhanced oil recovery with CO2 (CO2-EOR), food and beverage, and chemical production. CO2-EOR is the main end point for natural CO2 production in the United States; it accounted for 94% of usage in 2022 [1]. Federal lands emissions from this sector are estimated at 460 metric tons of CO2 in 2022, a very small amount relative to most other Federal lands emissions sector estimates. The second new method, planned for a separate report, is a calculation of the geologic storage of CO2 on Federal lands. The second method estimates the CO2 stored under Federal surface lands and documents Federal climate change mitigation efforts. Currently, there is no storage of CO2 at an industrial level on Federal lands, however multiple proposals and projects are planned. This method was developed on non-Federal lands datasets in an effort to prepare for when these activities on Federal lands will require accounting. National estimates for CO2 geologic storage using this method, but without a Federal lands filtering step, totaled 8.0 million metric tons (Mt) in 2022. The two methods described here are new benchmark methods in a collection of accounting procedures to document the current state of greenhouse gas emissions and their storage on Federal lands. These benchmarks can then be used to measure any subsequent changes in emissions from or carbon storage beneath Federal lands. While the magnitude of the values is currently non-existent to small, emissions mitigation goals established by decision makers indicate that these values will grow, and their documentation will take on greater value and use.

continental United States↗

USGS2025

Implications of physics-based M9 ground motions on liquefaction-induced damage in the Cascadia Subduction Zone: Looking forward and backward

Given the likelihood of future M9 Cascadia Subduction Zone (CSZ) earthquakes, various estimates of the resulting, regional ground motions have been made, including a suite of 30 physics-based simulations that reflect key modeling uncertainties. However, because the last CSZ interface rupture occurred in 1700 CE, the shaking expected in such an event is especially uncertain, as are the impacts to the built and living environments. Like other coseismic impacts, soil liquefaction poses a significant threat and must be considered by any scenario study used to inform planning and response, or to focus mitigation resources. Liquefaction is also notable for its potential to “ground truth” ground-motion estimates, given that its presence or absence in the geologic record can provide constraint on the intensities of shaking in past events. It is thus an important phenomenon looking both forward and backward. Accordingly, using recent physics-based simulations, this study (1) predicts liquefaction in M9 CSZ ruptures at 400 locations in Oregon, Washington, and British Columbia (BC) using an array of cone-penetration-test based models and (2) uses paleoliquefaction evidence at ten sites spanning from Southern Oregon to Vancouver, BC to constrain possible ground-motion intensities experienced in the 1700 CE earthquake. The forward predictions indicate that liquefaction in M9 events could be pervasive in the region and affect numerous population hubs, with the potential for damage across hundreds of square kilometers. The backward analyses suggest that 1700 CE ground-motion intensities may have been less than expected from M9 simulations in some northern portions of the CSZ (e.g. Seattle), given the paucity of 1700 CE liquefaction evidence in these areas. Ultimately, further discovery and analysis of CSZ paleoliquefaction, or lack thereof, will confirm or modify this possibility and the conclusions drawn herein.

Oregon, Washington↗

USGS2025

Comparing magmatism and hydrothermal alteration using magnetic modelling and stable isotopes at the Clementine porphyry copper prospect, Montana, USA

Recent mapping of the underexplored Clementine prospect in southwestern Montana has revealed evidence of hydrothermal alteration and mineralized breccia vein gossan interpreted to represent the upper expression of a Cretaceous, sediment-hosted copper porphyry system. The prospect is at the nexus of several Cretaceous granites, including the pre-mineralization Butte Granite and Big Hole Canyon plutons and local pre- to post-mineralization granite intrusions. Here, 3D magnetic inversions and stable isotope data are used to evaluate the spatial and genetic relationship of mineralization and Cretaceous magmatism. Magnetic inversions reveal a zone of high magnetic susceptibility beneath the prospect that is related to an exposure of unaltered, likely post-mineralization granite. The granite appears to connect to the Big Hole Canyon pluton at depth, suggesting it is a late-stage differentiate of the cooling, pre-mineralization pluton that may have been one of several potential heat sources for the mineral system at Clementine. This is supported by carbon isotope values of graphite that indicate organic material in shale of the Cretaceous Kootenai Formation was subject to metamorphism. Preliminary light hydrogen isotope values of actinolite hydrothermal orbs suggest formation following magmatic outgassing during thermal advection from an underlying intrusion.

Montana↗

USGS2024

Geophysical characterization of an alkaline‑carbonatite complex using gravity and magnetic methods at Magnet Cove, Arkansas, USA

The Magnet Cove alkaline‑carbonatite complex (MCC), located in the Ouachita Mountains of south-central Arkansas in the United States, hosts an extensive variety of rare rock types and critical mineral resources with physical properties (density and magnetization) that contrast significantly with the sedimentary rocks into which they have intruded. Newly acquired ground-based gravity and magnetic data were used to develop two-dimensional and three-dimensional geophysical models of the Cretaceous-aged Magnet Cove intrusive complex. The models reveal that the MCC: (1) widens out at middle crustal depths to as much 22 km across, and may reach a depth of 20 km; (2) has a total volume (exposed and subsurface) that may be over 800 km 3 ; (3) is likely connected at depth to other intrusions in the Arkansas alkaline province; and (4) has a geometry that is aligned with pre-existing structures such as the Reelfoot rift and the Ouachita orogenic belt, some of which were likely structurally controlled by the Precambrian crystalline basement and the continent-ocean transition zone buried beneath the Ouachita orogen. For the first time, the magnetic models of the MCC account for the presence of strong remanent magnetization. This results in a geophysical workflow necessary to accurately interpret magnetic anomalies over the much larger Arkansas alkaline province, its geologic and structural framework, and critical mineral potential.

Arkansas↗

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↗

USGS1992

Public water-supply systems and water use in Tennessee, 1988

This report summarizes the results of a study conducted by the U.S. Geological Survey, in cooperation with the Tennessee Department of Environment and Conservation (TDEC), Division of Water Supply in 1988. Data gathered during an inventory by the TDEC were collated to determine water use, supply sources, population served, and design and storage capacities of the systems. The inventory was limited to systems that were active on June 30, 1988. Results of a survey of the systems conducted by the Tennessee Department of Health and Environment during 1988 were a primary source of data for this report. Data from computer and manual files maintained by the Tennessee Department of Health and Environment and the U.S. Geological Survey also were used. The Division of Water Supply, TDEC, surveyed 541 public water-supply systems. These systems served 81 percent of the population of the State, or 3.95 million people. The gross per capita use statewide for public-supplied water was 179 gallons per day. Total water withdrawals for public supply increased about 39 percent from 510 million gallons per day (Mgal/d) in 1980, to 708 Mgalld in 1988. During the same period, the population increased about 7 percent. Surface-water withdrawals accounted for 63 percent (446 Mgal/d) of the total water withdrawn in the State. All of these withdrawals occurred in the Tennessee (56 percent or 249 Mgal/d) and the Ohio (44 percent or 197 Mgalld) hydrologic regions. Ground water supplied 262 Mgal/d or 37 percent of the total water withdrawn by public-supply systems statewide. Of that amount, 79 percent, or 208 Mgalld, was used in western Tennessee.

Tennessee↗

USGS2025

Accurate simulation of flow through dipping aquifers with MODFLOW 6 using enhanced cell connectivity

In simulations of groundwater flow through dipping aquifers, layers of model cells are often “deformed” to follow the top and bottom elevations of the aquifers. When this approach is used in MODFLOW, adjacent cells within the same model layer are vertically offset from one another, and the standard conductance-based (two-point) formulation for flow between cells does not rigorously account for these offsets. The XT3D multi-point flow formulation in MODFLOW 6 is designed to account for geometric irregularities in the grid, including vertical offsets, and to provide accurate results for both isotropic and anisotropic groundwater flow. A recent study evaluated the performance of the standard formulation and XT3D using a simple, synthetic benchmark model of a steeply dipping aquifer. Although XT3D generally improved the accuracy of flow simulations relative to the standard formulation as expected, neither formulation produced accurate flows in cases that involved large vertical offsets. In this paper, we explain that the inability of XT3D to produce accurate flows in the steeply dipping aquifer benchmark was not due to an inherent limitation of the flow formulation, but rather to the limited cell connectivity inherent in the most commonly used discretization packages in MODFLOW 6. Furthermore, we demonstrate that XT3D is able to produce the expected accuracy when adequate cell connectivity is introduced using MODFLOW's unstructured grid type and the aquifer is discretized vertically using at least two model layers.

Groundwater↗

USGS2025

Assessment of channel morphology, hydraulics, and bedload transport along the Siletz River, western Oregon

Significant Findings Chinook salmon ( Oncorhynchus tshawytscha ) and Pacific lamprey ( Entosphenus tridentatus ) are native, anadromous fish species in the Siletz River Basin, western Oregon, that face many threats to their survival in freshwater and the ocean. The Confederated Tribes of Siletz Indians of Oregon seek to mitigate freshwater threats to Chinook salmon and Pacific lamprey, where possible, with habitat conservation and restoration efforts. This study was conducted to assist the Confederated Tribes of Siletz Indians of Oregon in documenting and understanding the hydrogeomorphic processes shaping present-day habitat conditions and assessing future habitat implications for Chinook salmon and Pacific lamprey along the main-stem Siletz River. As such, this study focused on understanding geomorphic processes and patterns of channel change, including lateral and vertical adjustments in channel position and changes in bed-material sediment (sands, gravels, and cobbles that mantle the channel bed), which collectively determine overall patterns of channel morphology and fluvial habitats. Objective One was to evaluate lateral changes in channel position, vertical changes in bed elevation, and longitudinal patterns in bed-material particle size along the Siletz River using detailed channel maps developed from aerial photographs collected from 1939 to 2016, long-term records of stage and discharge collected by the U.S. Geological Survey (USGS) near the City of Siletz, and sediment particle size data. Objective Two was to assess hydraulic conditions using one- and two-dimensional hydraulic models and transport capacity of bed-material sediment using bedload transport models and sediment particle size data for a range of discharge conditions. Objective Three was to identify potential burrowing habitat for lamprey larvae (PBH) along the Siletz River network and provide insights in local factors influencing PBH along the main-stem Siletz River. The overall findings are synthesized to describe habitat implications for Chinook salmon and Pacific lamprey under present-day and future conditions. Results of Objective One, an evaluation of changes in channel position and bed elevations and longitudinal patterns in bed-material particle size along the Siletz River, include the following From 1939 to 2016, channel planform and the mapped area of gravel bars did not change considerably along the 97.2-kilometers (km) of Siletz River between Elk Creek and Millport Slough, except for in short sections generally less than 1 km long. This inherent lateral channel stability results from the resistant bedrock and terraces that bound most of the channel and limit lateral changes in channel position. Intermittent sections along the study reach where the active channel widened at channel bends displayed noticeable planform changes and increases in mapped bar area. From water year (WY; a 12-month period from October 1 through September 30 and named for the year in which it ends) 1906 to 2021, changes in the stage-discharge relation interpreted as rising and lowering channel bed levels were observed at the USGS streamgage 14305500 (Siletz River at Siletz, OR) in response to floods (such as high flows in February 1996 and January 2002 that exceeded 0.1 and 0.667 annual exceedance probabilities [AEP] events, respectively). However, the rating curve representing the stage-discharge relation did not change in response to high-magnitude floods between 2007 and 2021. Along the approximately 54-km of the Siletz River between Moonshine Park and the Bulls Bag area, surficial particle distributions varied considerably between sampling sites in response to changes in channel width and gradient, sediment inputs from tributaries, and basin geology. Despite this variability, median particle sizes tended to decrease in size in the downstream direction over the 54 km. Results of Objective Two, an evaluation of hydraulic and bedload transport conditions along the Siletz River, include the following The most substantial increases in maximum and mean water velocity and bedload sediment transport capacity occurred at events between the 0.995 and 0.50 AEPs. Events of these magnitudes occur approximately every 1–2 years. Smaller events (0.995 AEP) are generally contained by the banks of the main channel, whereas larger events (0.50 AEP) generally spill over the banks and inundate high-elevation bar and low-elevation floodplain surfaces. Multiple smaller floods within a WY that exceed the 0.995 and 0.50 AEPs (such as in WY 1996) can transport as much or more bed-material sediment than a single, higher magnitude event (such as the maximum event recorded in WY 2000 with an estimated AEP of around 0.002). Bedload transport capacity generally exceeds sediment supply (greater than 2 millimeters [mm]) for most of the study area from Wildcat Creek to the City of Siletz as evidenced by substantial in-channel bedrock. Despite overall conditions of limited sediment supply, transport capacity still varies considerably within and between years with discharge magnitude and spatially in relation to local hydraulics imposed by bedrock, channel morphology, and human infrastructure (such as bridges). Results of Objective Three, an analysis of PBH for lamprey larvae, include the following About 28 percent of the Siletz River network meets the mean annual suspended sediment loads and channel slope criteria for PBH. Along the main-stem Siletz River, in-channel bedrock outcrops and high transport capacity are expected to further constrain PBH. Together, these results suggest that most of the Siletz River between Wildcat Creek and the City of Siletz has had only modest vertical and lateral change between the 1930s and 2010s because of the bedrock in and along the main channel and the river’s relatively high transport capacity relative to bed-material sediment supply. However, localized sections of the Siletz River where the active channel widens, particularly at channel bends, exhibited some change in channel planform and the locations and area of gravel bars. In the future, moderate increases in autumn-winter discharge may not result in substantial changes in coarse gravel bars along the Siletz River but may result in selective transport of finer bed-material sediment (gravel, sands, and silts) that provide spawning habitats for Chinook salmon and Pacific lamprey and burrowing habitats for lamprey larvae. Assuming no substantial changes in bed-material sediment supply, increased bedload transport capacity may cause frequent entrainment of lamprey larvae that are burrowed in coarse sand deposits, suspension and downstream transport of salmon eggs incubating in gravels, and reductions in the areas of spawning gravels for Chinook salmon and Pacific lamprey. Exact implications of current and future discharge conditions for these species along the Siletz River depends on many factors, including sediment supply, local hydraulics, and the timing of flood events relative to fish life stages.

Oregon↗

USGS2026

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↗

USGS2025

Spatio-temporal evolution of distributed volcanic fields, case studies—Sierra Chichinautzin and Michoacán-Guanajuato, México

An analysis of 1,375 volcanoes in the Michoacán-Guanajuato (1,148 volcanoes in a 26,200 square-kilometer area) and Sierra Chichinautzin (227 volcanoes in a 3,500 square-kilometer area) volcanic fields in central Mexico identified patterns in the spatial and temporal distribution of past eruptions. A cluster agglomerative hierarchical method and kernel analysis confirmed that the Michoacán-Guanajuato volcanic field comprises four volcanic fields (Valle de Santiago, Uruapan, Apatzingán, and Pátzcuaro volcanic fields) controlled by different fault systems, indicating that it is not a single volcanic field but rather a group of volcanic fields (a “superfield”), each of which has distinct characteristics. In the Sierra Chichinautzin volcanic field, well-constrained isotopic ages were used to build a model of how the spatial distribution of the eruptions has changed over time. Two new 40Ar/39Ar ages from a locally recognized volcanic feature near the town of El Cantil, herein called El Cantil volcano (1,537±17 kilo-annum [ka]) and the volcanic feature at Cerro el Elefante (herein called El Elefante dome) (1,485±92 ka) belong to the oldest volcanic group identified in the Sierra Chichinautzin volcanic field, confirming the timing of the beginning of monogenetic volcanism in the region. Based on the volcanic groups identified in the Sierra Chichinautzin volcanic field, the youngest volcanism (less than 35 ka) is found only in the central-western sector of the field. Principal component analysis determined the directional trends of feeder dikes only for vents <10 ka in the Sierra Chichinautzin volcanic field. Possible magma migration paths through the crust were identified using seismic data from both volcanic fields using an earthquake catalog from 1973 to 2023, which includes 9,016 earthquakes in the Michoacán-Guanajuato volcanic field and 841 in the Sierra Chichinautzin volcanic field. The spatial distribution of the hypocenters does not highlight any trend that could be associated with superficial movement of magma in the Sierra Chichinautzin volcanic field. In the Michoacán-Guanajuato volcanic field, however, eight seismic swarms since 1997 have been detected. These swarms are interpreted to result from ascending magma. Strengthening monitoring systems and reinforcing mitigation measures to address volcanic hazards and risk are important means of preparing for future eruptions in both regions. Analysis such as those herein provide insights as to where an eruption might occur and may help mitigate volcanic hazards.

Michoacán-Guanajuato, Sierra Chichinautzin↗

USGS2025

Variable partitioning of lithium in rhyolitic melt during decompression and ascent

The partitioning behavior of Li in magmatic systems is increasingly being investigated due to the economic importance of Li in the transition to sustainable energy resources (e.g., batteries). However, at upper crustal pressures, it remains uncertain whether Li preferentially partitions into the vapor or liquid (brine) phase or remains in the silicate melt. This complicates our ability to determine where Li resides—silicate melt, minerals, or fluid phase—upon eruption, a crucial factor for understanding its postdepositional movement and concentration into a brine or volcano-sedimentary deposit. Here, we present a novel investigation into the behavior of Li within natural evolved melts during continuous magma decompression and ascent using melt embayments (open melt inclusions). Mineral-hosted melt embayments preserve records of the evolving composition of the exterior melt, including degassing pathways and ascent timescales, when paired with appropriate diffusion coefficients. Lithium concentration profiles were measured in quartz-hosted melt embayments from the rapidly quenched eruptive phases of five rhyolitic, caldera-forming eruptions to investigate the behavior of Li during magma decompression and ascent, where vapor partitioning and ascent dynamics were previously established by investigating H 2 O and CO 2 profiles. We find that in four systems, embayments contain lower interior Li concentrations than the coerupted melt inclusions; the fifth system contains the same Li concentrations in embayments and melt inclusions. However, many of these embayments contain gradients, with 84 % preserving Li enrichment near the melt-bubble interface, as compared to their interior concentration. We interpret these characteristics to represent two distinct stages of Li partitioning during magma decompression and ascent, in contrast to existing literature that proposes only one type of partitioning behavior. The first stage is interpreted as melt depletion of Li, likely driven by partitioning into an exsolved supercritical fluid phase, supported by the strong correlation between the extent of Li depletion and Cl concentration in the melt, as well as the decompression rate. This behavior then fundamentally shifts, where Li reenriches in the melt, postulated to be driven by the unmixing of the supercritical fluid phase at shallow pressures. For the one system that did not develop Li gradients through decompression, we attribute this to the lower values of Na and Cl in the melt, potentially inhibiting the partitioning of Li into a fluid phase. Importantly, the behavior of Li during decompression is not consistent within or between volcanic centers, highlighting the need for systematic experimental investigation in variable composition melts at pressures relevant to conduit dynamics. This knowledge would improve our ability to model Li profiles to understand magma decompression, and predict where Li resides (e.g., stored in volcanic glass, gas, or crystals) upon eruption prior to any later extraction.

Economic Geology↗

USGS2024

Comparison of water quality in shallow groundwater near agricultural areas in the Delaware Coastal Plain, 2014 and 2019

The State of Delaware has encouraged agricultural conservation practices to improve nutrient uptake by crops and mitigate nutrient transport to groundwater in the surficial aquifer. To study recent changes in groundwater quality, the U.S. Geological Survey and the Delaware Department of Agriculture (DDA) developed a network of shallow wells near agricultural areas throughout the Delaware Coastal Plain. This network was designed to characterize water quality related to agricultural practices and to detect any recent changes in shallow groundwater quality, in particular groundwater nitrate concentrations. The shallow well network was first sampled in 2014 and resampled in 2019. In 2019, field parameters (including dissolved oxygen, pH, specific conductance, and temperature), major ions, nutrients, stable isotopes of water, and isotopes of nitrate were measured in groundwater samples collected between October and December. Wells were organized into three groups based on their geochemical characteristics measured in 2014: the Agricultural, Urban, and Mixed Groups. Results from the 2019 sampling showed little change in water quality from the 2014 sampling. Land-use factors continued to be the driving influence between groups. Groundwater moves slowly and changes in groundwater quality are likely to respond slowly to changes in conservation practices. Continued sampling of both groundwater quality in this network and monitoring land management practices can help detect groundwater quality trends in the future.

Delaware↗

USGS2024

Hydroacoustic observations reveal drivers of mixing and salinization of a karst subterranean estuary during intense precipitation

Karst subterranean estuaries within globally ubiquitous carbonate aquifers are coastal groundwater ecosystems that provide an essential water resource for human populations. To understand the drivers of salinization within a coastal aquifer in the Yucatan Peninsula (Mexico), we employed hydroacoustics in flooded caves to observe how oceanic and atmospheric events facilitate mixing between the meteoric lens (fresh-brackish groundwater) and the saline groundwater on tidal and episodic timescales. Precipitation during Tropical Storm Carlotta increased the flow and salinity of the meteoric lens without evidence for vertical mixing across the halocline. We postulate that vertical migration of haloclines in the conduit relative to those within the rock matrix during precipitation creates lateral density gradients that drive mixing, and ultimately creates a brackish layer within the meteoric lens. These results provide a mechanistic explanation for vertical and lateral exchange in a coastal carbonate aquifer, which has implications for groundwater response to future climatic change.

Ox Bel Ha cave system, Yucatan peninsula↗

USGS2026

Bathymetric survey and storage capacity of Upper Lake Mary near Flagstaff, Arizona in 2024

The U.S. Geological Survey (USGS), in cooperation with the city of Flagstaff, collected bathymetric, light detection and ranging (lidar), and land-survey data of Upper Lake Mary in Arizona during the months of April and October 2024. The city of Flagstaff uses a combination of groundwater from well fields throughout the Flagstaff area and surface water, mainly from Upper Lake Mary, for its potable water supply. The purpose of the survey is to update previous surveys using new technology and compare the results to previous surveys to determine if there was a decrease in storage capacity that could affect the city’s water supply. The lakebed was mapped in April 2024 using a vessel equipped with a multibeam echosounder (MBES) and mobile lidar scanner with positioning captured using a real-time kinematic global navigation satellite system (RTK GNSS) base and receivers. In October 2024, areas of the reservoir that were too shallow for the boat and shoreline that were not captured by the vessel-based lidar were surveyed on foot using hand-held RTK GNSS receivers. At full pool (spillway elevation of 6,831.82 feet above NAVD 88 [2,082.34 meters (m)], Upper Lake Mary has a storage capacity of 16,449.80 acre-feet (20,290,611.73 cubic meters) and a surface area of 953.57 acres (3,860,926.075 square meters). The reservoir is 5.7 miles (9.7 kilometers) long and varies in width from 326 feet (99.36 m) near the central, narrow portion of the reservoir to 2,613 feet (796.44 m) in the upper portion. Comparisons between this survey and the previous two surveys from the 1950s and 2006 indicate no apparent decrease in reservoir area or storage capacity. Results of the 2024 survey indicate that Upper Lake Mary’s storage capacity increased by 0.9 percent from the 2006 survey and a 1.6 percent increase in surface area from the 2006 survey.

Arizona↗