USGS Science⌕ Search

SEARCH · USGS Science

Results for “PLoS ONE”

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,186 records · Page 54Linked to original sources
USGS2025

Machine learning provides reconnaissance-type estimates of carbon dioxide storage resources in oil and gas reservoirs

Oil and gas reservoirs represent suitable containers to sequester carbon dioxide (CO 2 ) in a supercritical state because they are accessible, reservoir properties are known, and they previously contained stored buoyant fluids. However, planners must quantify the relative magnitude of the CO 2 storage resource in these reservoirs to formulate a comprehensive strategy for CO 2 mitigation. Even reconnaissance-type estimates of CO 2 storage resources of known oil and gas reservoirs may require complicated calculations involving 1) estimates of recoverable oil and gas, 2) reservoir properties (depth, temperature, pressure, etc.), and 3) the physical qualities of the retained fluids. We demonstrate the application of machine learning (ML) algorithms to bypass these computations to yield more rapid estimates of CO 2 storage resources in reservoirs capable of hosting CO 2 in a supercritical state. ML algorithms are computationally efficient because they do not impose the strong assumptions on the data-generating process that standard statistical or engineering procedures require. Further, ML algorithms can capture highly complex, particularly nonlinear, relationships among predictor variables. We demonstrate the application of four different ML algorithms using data from onshore and offshore oil and gas reservoirs in Europe, and show they perform well when predictions are compared to engineering estimates. The proposed methods and models provide an effective and novel way to more rapidly and directly determine the subsurface CO 2 storage capacity of oil and gas reservoirs around the world, information that operators, researchers, and policymakers alike require to meet energy transition and decarbonization goals.

Frontiers in Enviornmental Science↗
USGS2026

Reconstructing the Quaternary depositional history using geologic mapping and three-dimensional modeling of the subsurface near Fort Morgan, northeastern Colorado

Centered on Fort Morgan, Colorado, this study is intended to build from previous work by adding a three-dimensional (3D) view of the subsurface to better understand the depositional history of Quaternary deposits. A 1:100,000 scale geologic map was made by combining previous geologic maps, regional soil maps, and recent field investigations. In addition to the geologic mapping, drill hole lithologic data from water wells and oil and gas exploration were compiled and lithologic units simplified to best represent the stratigraphy of the Quaternary deposits. From these subsurface data, a 3D subsurface model was constructed, trimmed at the surface by a digital elevation model, and a bedrock surface foundation gridded from drill hole data was added. The surface of the 3D model was then compared visually to the surficial geologic map. Cross sections were constructed from the 3D model and compared to site-specific drilling that was done as part of this project. Finally, the model was examined in detail to reconstruct the depositional history of the subsurface alluvial and eolian units. Alluvial and fluvial drainage basins exposed in the subsurface have a greater areal extent than the present-day narrow drainages. Older eolian sand in the subsurface tends to be interbedded with loess indicating coeval deposition. Holocene sand, both eroded from bedrock exposed at the surface north of the study area and reworked from the South Platte River, buries most of the interbedded older sand and loess.

Colorado↗
USGS2025

False positives in the identification of dynamic earthquake triggering

Dynamic earthquake triggering is commonly identified through the temporal correlation between increased seismicity rates and global earthquakes that are possible triggering events. However, correlation does not imply causation. False positives may occur when unrelated seismicity rate changes coincidently occur at around the time of candidate triggers. We investigate the expected false positive rate in Southern California with global M ≥ 6 earthquakes as candidate triggers. We compute the false positive rate by applying the statistical tests used by DeSalvio and Fan (2023), https://doi.org/10.1029/2023jb026487 to synthetic earthquake catalogs with no real dynamic triggering. We find a false positive rate of ∼3.5%–8.5% when realistic earthquake clustering is present, consistent with the 95% confidence typically used in seismology. However, when this false positive rate is applied to the tens of thousands of spatial-temporal windows in Southern California tested in DeSalvio and Fan (2023), https://doi.org/10.1029/2023jb026487 , thousands of false positives are expected. The expected false positive occurrence is large enough to explain the observed apparent triggering following 70% of large global earthquakes (DeSalvio & Fan, 2023, https://doi.org/10.1029/2023jb026487 ), without requiring any true dynamic triggering. Aside from the known triggering from the nearby El Mayor-Cucapah, Mexico, earthquake, the spatial and temporal characteristics of the reported triggering are indistinguishable from random false positives. This implies that best practice for dynamic triggering studies that depend on temporal correlation is to estimate the false positive rate and investigate whether the observed apparent triggering is distinguishable from the correlations that may occur by chance.

JGR Solid Earth↗
USGS2025

The Mammoth magnetic anomaly, Pinal County, Arizona

A high-resolution Earth Mapping Resources Initiative airborne geophysical survey was flown in the southwest North American porphyry copper province to improve bedrock geologic maps and to identify areas that have unrecognized critical mineral resource potential. During the review of the aeromagnetic data, a distinctly monopolar-shaped, negative magnetic anomaly was observed at a flight elevation of 200 m above the ground with a maximum amplitude of –9500 nT. We have named this the Mammoth magnetic anomaly (MMA) because it is centered 12 km northeast of the town of Mammoth, Arizona, USA. The total field anomaly (TFA) contour of –500 nT enclosing the MMA defines an elongate shape measuring 2.5 km long by 1 km wide that trends northwest–southeast. Given the striking nature of this negative, monopolar-shaped magnetic anomaly, we conducted a ground campaign in May 2025 to determine its authenticity and potential relationship to critical mineral endowment. The MMA was confirmed on the ground with a TFA approaching –46,000 nT. Total magnetic intensity (TMI) observations routinely fell below the 18,000 nT operating floor of an industry-standard cesium-vapor total field magnetometer, and extremely low TMI measurements were corroborated along coincident traverse lines using two high dynamic range, but lower sensitivity, smartphone vector magnetometers. The lowest TMI values recorded by both smartphone magnetometers were 1000 nT and confirmed with multiple adjacent and crossing lines. Field observations suggest that this magnetic feature is caused by strong remanent magnetization within fine-grained magnetite hosted within locally altered Pinal Schist.

Arizona↗
USGS2025

Sand provenance boundary in the Mu Us Sandy Land of northern China

Desert dunes are often assumed to have uniform mineral compositions due to extensive mixing during lateral transport, which complicates provenance studies. The Mu Us Sandy Land in north-central China, near the East Asian summer monsoon precipitation boundary, experiences a wetter climate than most deserts. Climate wetting as a result of a warming climate, and the ‘Sand Control Project’ implemented by the Chinese government to decrease the lateral movement of sand dunes in this area provide an opportunity to study surface processes of sand production and transport. Previous studies using zircon U Pb geochronology and heavy mineral composition indicate distinct sand sources for the Mu Us Sandy Land: local basement-derived middle Yellow River sediments and recycled dried-up lacustrine sediments for the eastern part, and northeastern Tibetan Plateau-derived upper Yellow River sediments for the western part. However, zircons and heavy minerals only represent trace amounts of the bulk mineralogy within a sand dune, so broader provenance analysis targeting common minerals is essential. We focus on the area near the proposed provenance difference boundary between the western and eastern Mu Us Sandy Land, using comprehensive sampling and provenance techniques to confirm distinct provenance characteristics and delineate the provenance boundary. Our findings reveal that sand from most of the Mu Us Sandy Land originated from the erosion of local basement by the middle Yellow River and recycled local dried-up lacustrine sediments, whereas the southwestern corner and the neighboring western-central Chinese Loess Plateau received sediments from the distal northeastern Tibetan Plateau.

Mu Us Sandy Land↗
USGS2024

Seismic response comparison of a historical masonry church subject to real and simulated ground motions

In recent years, advanced numerical models and high-performance computing have facilitated the utilization of ground motion time series in the assessment of the non-linear dynamic behavior of historic masonry structures. Since recorded accelerograms can be sparse for specific analysis conditions, stochastic ground motion simulations have become a viable alternative to overcome this limitation. This study simulates the recorded acceleration time series of the Central Italy 2016 earthquake event at the closest station to the town of Macerata using a site-based stochastic approach. The simulated motions are seismologically evaluated using a goodness-of-fit method in terms of various intensity measures. The simulated records, in conjunction with real records, are used to study the non-linear dynamic behavior of San Filippo Neri church located in Macerata. The church of San Filippo represents an important example of Baroque religious architecture in central Italy, which was damaged and closed off to the public after the 2016 earthquake events. The construction was investigated with a vast diagnostic campaign which included on-site testing and dynamic identification tests. The collected data is used to calibrate the dynamic response of a three-dimensional finite element model of the church. The model is finally used to compare the non-linear seismic responses under real and simulated ground motions with the site recorded damage. The results of structural responses demonstrate a strong agreement between the real and simulated records, providing evidence to support the validation of the site-based stochastic simulation.

Macerata↗
USGS2025

Beyond the mangroves: A global synthesis of tidal forested wetland types, drivers and future information opportunities

There is increasing awareness of the global diversity of tidal forested wetlands (TFWs) and their significance in the provision of ecosystem services. These ecosystems, including mangrove forests, tidal freshwater forested wetlands, supratidal forests and transitional forests together span multiple climatic zones, geomorphic settings, and inundation and salinity regimes. We utilise case studies across five continents to demonstrate the state of knowledge among TFWs. Intertidal mangroves are the best-defined of the TFWs thanks to decades of research on their geomorphology, hydrology and ecology across their broad distribution. Non-mangrove forest settings, however, demonstrate more diverse hydrological, biochemical and vegetation conditions. In many cases, non-mangrove forests are situated at upper intertidal or supratidal elevations, where surface waters and groundwater are subject to interactions between tides freshwater inputs. Salinity datasets show variations ranging from tidal freshwater forested wetlands and ‘low-salinity mangroves’ to mesohaline or marine salinities, often with high temporal variability. While the floristic composition of non-mangrove forests vary among biogeographic regions, locally dominant TFW species are commonly distributed beyond the tidal niche into non-tidal wetland and upland forests. This presents challenges for traditional remote sensing approaches to ecosystem mapping, which are mostly lacking for non-mangrove forests. Geomorphic approaches and developments in machine learning offer opportunities to address this.

Preprint↗
USGS2026

Small earthquake moment magnitude and implications for frequency–magnitude scaling of injection induced earthquakes of the Raton Basin

Accurate estimation of earthquake source parameters—such as moment magnitudes, corner frequencies, and stress drops—is essential for improving seismic hazard assessments and understanding earthquake physics. In this study, moment magnitudes ( M W ) are calculated for 31,581 earthquakes associated with wastewater injection in the Raton Basin (located along the border between northern New Mexico and southern Colorado) between 2016 and 2024 using radiative transfer theory to fit coda decay envelopes. Our results show that it is feasible to estimate moment magnitudes down to M W ~1 with coda envelopes from a small local monitoring network. Significant differences were found between M W and local magnitudes ( M L ) for small earthquakes ( M < 3.0). A linear relationship was optimized to convert M L to M W : M W = 0.7 M L + 0.96 and M W = 0.73 M L + 0.99 (for the events reported by the U.S. Geological Survey), which can be applied in future studies of Raton Basin seismicity. We find that b -values calculated employing different methods and using M L are approximately 1.0, while those using M W range from 1.2 to 1.4. A larger estimate of the b -value could influence interpretations of the statistical behavior of earthquakes associated with injection and consequently seismic hazard assessments based on a magnitude–frequency distribution. The potential differences between local versus moment magnitude-based earthquake statistics should be considered in other seismically active regions.

Colorado, New Mexico↗
USGS2026

Hydrothermal manganese-oxide mineralization of a carbonate ooze, Samoan hotspot region, South Pacific Ocean

Low-temperature hydrothermal manganese oxides occur throughout the global oceans. However, the hydrothermal replacement of a carbonate ooze by manganese oxides is described here for the first time. The 24 samples dredged from three locations around the Territory of American Samoa in the South Pacific Ocean include six samples with remnant carbonate and volcaniclastic sediments, which we refer to as “low Mn,” and 18 samples in which the mineralization is pervasive and has replaced most original sediment, termed “high Mn.” The 18 high-Mn samples exhibit a mean Mn content of 51 wt.%. Higher Li and Mn contents and lower Fe contents in the high-Mn samples indicate a hydrothermal origin and distinguish these samples from hydrogenetic ferromanganese crusts. Mn-oxide layers are up to 90 mm thick, some with columns to 44 mm long and 10 mm wide, the magnitude of which has not been described previously. Samples are composed of birnessite and 10 Å phyllomanganate minerals. Textures of the thickest Mn layers indicate mineralization below the seabed from ascending fluids during multiple hydrothermal pulses. Mineralization took place by complete to partial replacement and cementation of foraminiferal sediments intermixed with volcaniclastic sediments in varying amounts. Our results highlight the production of carbonate sediment-hosted hydrothermal Mn oxides from multiple hydrothermal sources within the Samoan volcanic chain. The potential extensive distribution on the regional scale of this newly described mineralization process and unique element enrichments raise questions about its broader distribution globally and potential importance to hydrothermal processes, element mass balance, and seabed mineral resources.

Geochemistry, Geophysics, Geosystems↗
USGS2024

Seismic velocity changes from repetitive seismicity at Mauna Loa prior to and during its 2022 eruption

Mauna Loa’s short-lived eruption from late November to early December 2022 marked the culmination of nearly a decade of elevated seismic activity and geodetic inflation. The volcano has been monitored by a network of permanent, short period and broadband seismometers. I used the continuous waveform data from that network starting in 2012 to generate a catalog of seismicity that enhances the US Geological Survey Hawaiian Volcano Observatory’s public seismic catalog with four times the number of earthquakes, which were then grouped by waveform similarity. Analysis of subtle delays in the timing of arrivals of scattered waves between pairs of earthquakes in this catalog yields a history of small changes in the shallow seismic velocity structure of the volcano. Seismic velocities have been shown at other volcanoes to change during unrest and eruption. My results show a decrease in seismic velocity centered on the summit beginning in September 2022, corresponding to the onset of a vigorous precursory swarm of seismic activity and shallow inflation. During the eruption itself, I observe large changes due likely to dike opening along the northeast rift zone and deflation of the summit reservoir. However, seismic velocity changes associated with non-volcanic sources such as ground shaking from large earthquakes and meteorological influences at seasonal and diurnal time scales are also observed, and these dominate the velocity changes prior to the eruption. Proper accounting of these effects will be a requirement for use in real-time monitoring, and this work serves as a starting point in that endeavor for Mauna Loa.

Hawaii↗
USGS2026

Hazard potential of compound flooding from rainfall, storm surge, and groundwater in coastal New York and Connecticut

Compound flood events, the co-occurrence of multiple flood drivers, can result in flood hazard potential exceeding that of any single driver alone. To evaluate compound flooding in a semi-urbanized coastal area, historical records dating back to 1970 are used to study the co-occurrences of high precipitation, storm surge, and shallow groundwater conditions along the coastlines of New York and Connecticut. Joint return periods for coincident precipitation-surge events were computed using statistical dependence models and compared to the assumption of independence as a ratio, referred to here as a return period adjustment. Results indicate distinct seasonality where compound events in the area disproportionately occur in the cold season between October and April. Return period adjustments range from a factor of 1 to almost 9, demonstrating the range in precipitation-storm surge dependence across the study area. Across all 24 station triad locations, groundwater levels were elevated during times of precipitation- surge co-occurrence, reflecting the tendency for coastal storms and shallow groundwater conditions to co-occur seasonally. The result is a pseudo-trivariate compound flood hazard score and corresponding hazard map that integrates dependence between daily precipitation-surge events and overall monthly groundwater levels (as a precondition) into a relative compound hazard score. The location with the highest compound flood hazard score is on the south shore of Long Island, as well as locations across coastal Connecticut where groundwater levels compound the co-occurrence of heavy precipitation and storm surge.

Connecticut, New York↗
USGS2025

Do bighorn sheep use desert tortoise burrow spoil piles as mineral licks in southern California?

In a previous study, bighorn sheep ( Ovis canadensis ) were periodically photographed by trail cameras when they visited desert tortoise ( Gopherus agassizii ) burrows in the San Bernardino Mountains near Palm Springs, California, USA. The authors suggested that bighorn sheep may utilize the excavated spoil material from burrows as a mineral lick. To test that supposition, we collected soil samples from the spoil piles and undisturbed surface soil up to 1 m away at two desert tortoise study sites with bighorn sheep at the aforementioned site and another nearby site. We hypothesized that sodium concentrations would be greater in spoil than in nearby surface soil samples. Common major cations were quantified in soil samples using inductively coupled plasma-mass spectrometry following microwave-assisted digestion to solubilize the fraction of elements that could become biologically or environmentally available. Mean sodium concentrations were greater in burrow spoil samples than in undisturbed soil samples at both sites, likely due to aqueous leaching of sodium from surface to subsurface soil horizons. Although we did not collect evidence of geophagy by bighorn sheep, the behavior is well known. However, trail camera photographs of bighorn sheep with their head next to the entrance of tortoise burrows suggests that they may periodically utilize sodium and other elements brought to the surface by desert tortoises that are otherwise unavailable to bighorn sheep on the surface.

California↗
USGS2024

From exploration to production: Understanding the development dynamics of lithium mining projects

Recently, there has been considerable recent controversy whether current and new lithium mines will be able to supply the rapidly growing needs of the electromobility transition. Mineral exploration projects are typically active for many years, and only some become operational mines. From exploration to production, the projects go through several stages of characterisation and evaluation. At each stage, decisions are made by companies and stakeholders to advance, continue or stop the project. This is a complex process, and even projects with very similar geological and technical characteristics may take very different trajectories, depending on external factors such as global market conditions and local regulatory environments. The present study investigates the dynamics of this process for lithium exploration projects. A global database of 397 lithium projects was compiled, covering their progression through major development stages between 2004 and 2022. Ordinal logistic regression was used for the statistical analysis of this data. Different explanatory variables were tested, including economic, geological, technical, and geographic factors, to identify the best predictors for project progress at each development stage. The results suggest an essential role for lithium carbonate prices, and a variable role for other factors at each stage. Critically, the already elapsed lead time and project economics, which are traditionally considered important for the prediction of the start-up of individual mines, do not appear to be relevant in all cases. The results provide important insights into the dynamics of lithium supply and may eventually allow more realistic forecasts to be made for future lithium market dynamics.

Resources Policy↗
USGS2026

The geologic history of the Chehalis Forearc Basin, Washington State, USA

The Chehalis basin is located between the Cascade arc and the Coast Range in southwest Washington State. It consists of sedimentary and volcanic rocks deposited throughout the Cenozoic and is underlain by the Siletzia terrane, a thick, large igneous province accreted in the Eocene. Here, we constrain evolution of the Chehalis basin depocenter using isochore maps derived from isostatic gravity anomalies, borehole data, and stratigraphy for several time periods: the Eocene (47.6–36.8 Ma), the Oligocene (36.8–20 Ma), and the Neogene (20–0 Ma). Our results suggest that local subsidence is driven by deformation on faults that bound and intersect the basin. We see northward depocenter migration, interpreted as a shift in deformation from the northwest striking Cedar Creek fault to the west striking Doty fault. We interpret these data in terms of the long-term north-south shortening and clockwise rotation of the Cascadia forearc. During the Eocene, the Cedar Creek fault was preferentially aligned with north-south shortening, but became less active when clockwise rotation brought the Doty fault into east-west alignment with the stress field. Sediment accumulation rates decreased from 196 m/Myr in the Paleogene (∼40–20 Ma) to 27 m/Myr in the Neogene (∼20–0 Ma) as Miocene uplift of the Coast Range limited accommodation space and shifted the depositional environment from marine to fluvial. Our results are consistent with the geologic evolution of the Portland and Tualatin basins to the south and reveal a depositional and structural history uniquely shaped by clockwise rotation of the Cascadia forearc.

Washington↗
USGS2026

Egg thiamine concentrations reflect dietary differences between four hatchery strains, including two morphotypes, of Lake Trout in Lake Ontario

Objective Test whether dietary differences among hatchery strains of Lake Trout Salvelinus namaycush were associated with variation in egg thiamine concentrations and assess the potential for using strain-specific traits to mitigate thiamine deficiency complex (TDC) in Lake Ontario restoration efforts. Methods Stomach content analysis and stable isotope mixing models were used to compare diets among four hatchery strains of Lake Trout, representing lean and humper morphotypes. Egg thiamine concentrations were measured to test whether dietary differences correlated with susceptibility to TDC. Results Significant differences in diet and egg thiamine concentrations were observed among Lake Trout strains. A humper-like morphotype strain exhibited elevated thiamine levels, which were associated with reduced consumption of prey known to degrade thiamine levels when consumed, Alewife Alosa pseudoharengus and Rainbow Smelt Osmerus mordax , and increased reliance on Round Goby Neogobius melanostomus or sculpin species (Slimy Sculpin Cottus cognatus and Deepwater Sculpin Myoxocephalus thompsonii ). These findings support the link between the consumption of Alewife and Rainbow Smelt and thiamine deficiency. Conclusions Our findings support a link between strain-specific dietary preferences and susceptibility to TDC in Lake Trout. Preferential stocking of strains less prone to TDC, such as the strain derived from a humper morph from Lake Superior, could enhance restoration success in Lake Ontario. This strategy may be especially effective given the continued dominance of Alewife in the lake's prey base and the current management priorities supporting fisheries for Pacific salmon Oncorhynchus spp.

North American Journal of Fisheries Management↗
USGS2025

A regional model comparison between MODPATH and MT3D of groundwater travel time distributions

Groundwater quality changes in wells and streams lag behind changes to land use due to groundwater travel times. Two contaminant transport methods were compared to assess differences in their simulated travel time distributions (TTDs) to streams and wells in the Wisconsin Central Sands. MODPATH simulates advective groundwater flow with particle tracking, while MT3D simulates age-mass using a finite difference solution without dispersion to allow for direct comparison of the two methods. MODPATH appropriately simulates groundwater TTDs from the water table to surface discharge but is subject to inaccuracies at weak-sink well cells due to the flow-model grid discretization and imprecise location of well discharge within well cells. MT3D better represents weak-sink well cells since it removes mass in proportion to the prescribed pumping rate, although travel time within well cells is neglected. Conversely, MT3D's treatment of surface water boundary cells is not as accurate as MODPATH because mass should be removed from the water table rather than the full cell volume. MT3D simulations of TTDs can also be confounded by the instantaneous vertical distribution of mass introduced throughout recharge cells instead of at the water table, which initiates mass along deeper flow paths. We evaluated 9 MODPATH and 13 MT3D implementations, generating differences in median travel times of up to 18 years. Both methods have strengths and weaknesses, with MT3D better representing weak-sink well cell behavior and MODPATH better representing surficial recharge and discharge. The effect of these characteristics on simulated TTDs, along with ideas for ameliorating method weaknesses, is discussed.

Wisconsin↗
USGS2026

Potential warming-induced changes in stream water sources between proglacial and non-glacial streams, south-central Alaska

Snow and ice-melt are essential for global water resources. Rising air temperatures are causing vegetation encroachment and shrubification of previously snow- and ice-covered landscapes, and precipitation regimes are changing. Collectively, these changes will likely affect the hydrology of mountain environments, although the exact ways in which the hydrology will change are poorly understood. Changing source waters (i.e., the proportion of ice-melt, snowmelt, rain and groundwater) affect the timing and magnitude of streamflow and affect stream temperature, sediment, solute and nutrient fluxes. We investigate spatial and seasonal source water contributions across sub-watersheds of the Nellie Juan River watershed on Alaska's Kenai Peninsula, which differ in glacier cover, elevation and land cover, using a space-for-time approach. To do so, we use a three-endmember mixing model to quantify the contribution of rain, snow/ice-melt and deep groundwater. We show that snow/ice-melt peak later in proglacial streams compared to non-glacial streams. Endmember contributions in non-glacial streams generally vary with elevation and season. We also found that groundwater is a major contributor to streamflow across the study area, contributing between 12% and 47% seasonally, and may become increasingly important as snowpack and ice-melt decrease. This research characterises the shift from a glacierised to a deglaciating landscape through the lens of source water contribution. Broadly, these findings could help improve understanding of how water resources in glacierised watersheds are affected by warming air temperatures.

Alaska↗
USGS2024

Groundwater quality and groundwater levels in Dougherty County, Georgia, April 2020 through January 2023

The Upper Floridan aquifer is the uppermost reliable groundwater source in southwest Georgia. The aquifer lies on top of the Claiborne, Clayton, and Cretaceous aquifers, all of which exhibited water-level declines in the 1960s and 1970s. The U.S. Geological Survey has been working cooperatively with Albany Utilities to monitor groundwater quality and availability in these aquifers since 1977. Flow direction in the Upper Floridan aquifer is to the south and toward the Flint River. During the past 3 years, water levels varied above and below period-of-record median values. Water levels in the Upper Floridan aquifer were primarily above or at median levels during 2020 and 2021 and at or below median levels during 2022. Water levels in the Claiborne aquifer were above median levels, whereas water levels in the Clayton aquifer were at or below median levels, and in the Cretaceous aquifer system were close to median levels. During January 2021, eight wells were sampled for major ions, including nitrate plus nitrite as nitrogen (N). Nitrate plus nitrite as N concentrations ranged from 2.3 to 10.5 milligrams per liter (mg/L). During December 2021, seven wells were sampled for major ions, including nitrate plus nitrite as N. Nitrate plus nitrite as N concentrations ranged from 3.9 to 9.9 mg/L. During November 2022, eight wells were sampled for major ions, including nitrate plus nitrite as N. Nitrate plus nitrite as N concentrations ranged from 3.9 to 10.0 mg/L. Two wells were also sampled for per- and polyfluoroalkyl substances during November 2022.

Georgia↗