USGS Science⌕ Search

SEARCH · USGS Science

Results for “Philosophical Transactions of the Royal Society B: Biological Sciences”

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,687 records · Page 53Linked to original sources
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 ↗
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↗
USGS2024

Status of water-level altitudes and long-term and short-term water-level changes in the Chicot and Evangeline (undifferentiated) and Jasper aquifers, greater Houston area, Texas, 2023

Since the early 1900s, groundwater withdrawn from the primary aquifers that compose the Gulf Coast aquifer system—the Chicot, Evangeline, and Jasper aquifers—has been an important source of water in the greater Houston area, Texas. This report, prepared by the U.S. Geological Survey in cooperation with the Harris-Galveston Subsidence District, City of Houston, Fort Bend Subsidence District, Lone Star Groundwater Conservation District, and Brazoria County Groundwater Conservation District, is one in an annual series of reports depicting the status of water-level altitudes and water-level changes in these aquifers in the greater Houston area. In this report, the Chicot and Evangeline aquifers are treated as a single aquifer for the purposes of providing annual assessments of regional-scale water-level altitudes and water-level changes over time. In 2023, shaded depictions of water-level altitudes for the Chicot and Evangeline aquifers (undifferentiated) ranged from about 286 feet (ft) below the North American Vertical Datum of 1988 (NAVD 88) to about 169 ft above NAVD 88. The largest decline in water-level altitudes indicated by the 1977–2023 long-term water-level-change map was in south-central Montgomery County southeast of The Woodlands. In comparison, the 1990–2023 long-term water-level-change map depicts the largest declines in water-level altitudes in localized areas at or near certain wells in parts of northwestern Harris County and south-central Montgomery County. The largest rise in water-level altitudes for 1977–2023 is depicted in a relatively large area in southeastern Harris County, whereas the largest rise in water-level altitudes for 1990–2023 is depicted in a relatively large area in central Harris County. The 5-year short-term water-level-change map depicts the largest declines at three wells in northern Fort Bend County, one well in western Harris County, and three wells in south-central Montgomery County and the largest rise at one well in central Harris County. The 1-year short-term water-level-change map depicts the largest declines at one well in northern Fort Bend County and two wells in southwestern Harris County and the largest rises at one well in northern Brazoria County and one well in south-central Montgomery County. In 2023, shaded depictions of water-level altitudes for the Jasper aquifer ranged from about 242 ft below NAVD 88 to about 218 ft above NAVD 88. The 2000–23 long-term water-level-change map depicts water-level declines throughout the study area where water-level-measurement data from the aquifer were collected, with the largest declines in north-central Harris County and south-central Montgomery County south of The Woodlands. The 5-year short-term water-level-change map depicts the largest declines at two wells in central Montgomery County near Conroe and two wells in south-central Montgomery County southeast of The Woodlands and the largest rise at one well in western Montgomery County. The 1-year short-term water-level-change map depicts the largest declines at four wells in south-central Montgomery County southeast of The Woodlands and one well in central Montgomery County near Conroe and the largest rises at two wells in western Montgomery County.

Texas↗
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↗
USGS2026

A temporal look at the influence of topographic amplification on earthquake-triggered landslides in 3D seismic simulations

Earthquakes are a primary trigger for landslides, often leading to catastrophic consequences. While numerous studies have explored the spatial distribution of earthquake-triggered landslides, understanding the interaction between seismic waves and topography remains a critical challenge. Topographic irregularities can cause seismic wave amplification, altering ground shaking, and can trigger landslides that are challenging for predictive models to anticipate. This study investigates the spatial and temporal evolution of topographically amplified landslides, focusing on coseismic landslides triggered by the Mw 7.5 mainshock of the 2018 Papua New Guinea earthquake and post-seismic landslides associated with its four aftershocks, each exceeding Mw 6.0. We employ low-frequency, three-dimensional numerical ground shaking simulations and data-driven multivariate analyses to examine how landslides evolved from the coseismic to post-seismic periods. Our findings reveal a spatial shift in landslide distribution, in which the mainshock triggered slope failures predominantly on steep hillslopes, and the aftershocks triggered landslides on gentler slopes, often near geologic boundaries. We attribute this transition partly to the earthquake legacy effect of the mainshock, where the mainshock caused weakening of these hillslopes, making them more prone to failure when aftershocks occur. Additionally, the concentration of failures along geologic contacts in the post-seismic phase suggests that site amplification, stemming from contrasts in subsurface materials, exerts a key influence on landslide occurrence. Although not explicitly captured in our current numerical simulations, this mechanism warrants further investigation for more accurate hazard modeling.

JGR Solid Earth↗
USGS2025

Groundwater-storage change in the north Phoenix aquifer, Arizona, 2020–23

The city of Phoenix, Arizona, relies primarily on surface water for municipal water supply. The city also maintains wells to withdraw groundwater, particularly in times of drought and reduced surface-water supply, and to recharge groundwater when excess surface water is available. As of 2023, withdrawals from the aquifer in the northeastern part of the city are a small volume of water, less than 3,000 acre-feet in most years. Each year a similar volume of water is recharged through injection wells. The withdrawal wells are permitted to produce in total more than 20,000 acre-feet per year; increased pumping could lead to future groundwater-storage declines. To better understand groundwater-storage change in the north Phoenix aquifer, a repeat microgravity monitoring network was established in 2020. Measurements of changes in Earth’s gravity provide a direct, non-invasive measurement of subsurface mass change. Groundwater-storage changes were small during the 2020–23 study period, consistent with the relatively small volumes of pumping and recharge. Groundwater levels measured in monitoring wells were stable or increased slightly during this period, although the number of monitoring wells within the area of the gravity monitoring network is sparse. In total, about 15,000 acre-feet of water were pumped and 31,000 acre-feet recharged through injection in the north Phoenix aquifer during the 2020–22 period. Within the monitored area, groundwater storage increased by a small amount, about 1,000 acre-feet, in 2020, and decreased by a larger amount, about 6,000 acre-feet, each year in 2021 and 2022. Storage decreased at 89 of 102 stations from 2020 to 2023. Groundwater-storage decreases greater than the volume of net pumping indicate down-gradient subflow from the aquifer is greater than recharge plus incoming subflow, drying of the unsaturated zone resulting from decreased land-surface recharge, or both. At present (2023), the aquifer appears able to store and supply the relatively small amounts of water needed without excessive drawdown or harmful effects, such as land subsidence and surface inundation. If pumping or recharge increases in the future, the established repeat microgravity network is well suited to capture the spatial extent and magnitude of groundwater-storage changes in the aquifer.

Arizona↗