Geology topics

Research about United States

Source-linked reports with geographic coverage including United States.

1,087 records · Page 41Linked to original sources

USGS2024

Imaging of seismic discontinuities using an adjoint method

For imaging of seismic discontinuities at depth, reverse time migration (RTM) is a powerful method to apply to recordings of seismic events. It is especially powerful when an extensive receiver array, numerous seismic sources, or both, permit adequate reconstruction of incident and scattered wavefields at depth. Reconstructing either the incident or scattered wavefield at depth becomes less accurate when relatively few recordings of seismic events are available. Here we explore an inverse scattering approach to imaging discontinuities based on an adjoint method, employing sensitivity kernels (Frechet derivatives) that represent jumps in material properties across seismic-discontinuity surfaces. When combined with ray-based requirements on scattering geometry, it constitutes a powerful approach to determining the locations and amplitudes of the discontinuities, recovering only those properties that can be resolved by a spatially limited source and/or receiver distribution. This is illustrated by synthetic examples with local sources followed by a field example in a subduction zone setting

Washington↗

USGS2024

The geometry of fault reactivation and uplift along the central part of the Maacama fault zone, northern California Coast Ranges (USA)

Fault reactivation of bedrock structures in active fault zones influences stress state and earthquake rupture phenomena through the introduction of weak slip surfaces that impact fault zone geometry and width. Yet, geometric relationships between modern faults and older reactivated faults are difficult to quantify in rocks that have experienced multiple deformation episodes. We used new geologic mapping, geomorphic tools, and structural modeling to quantify rock uplift and subsurface fault geometry of the central part of the Maacama Fault Zone near Ukiah, California, USA, and the surrounding area. Results suggest that the northern Mayacamas Mountains are in a tectonically driven disequilibrium, with differential rock uplift focused on the western side of the range. Steeply east-dipping fault surfaces and splays characterize the geometry of the Maacama Fault Zone. We mapped two newly identified faults to the east of the main Maacama Fault, the Cow Mountain–Mill Creek Fault, and Willow Creek Fault, which align with a moderately east-dipping cluster of microseismicity between 4–10 km depth beneath the Mayacamas Mountains. Static stress modeling on the Maacama Fault Zone and newly identified faults to the east quantify slip tendency values of 0.5–0.4, which suggests that the faults are moderately to poorly suited for slip in the modern stress field and may be weak. We infer that modern uplift is driven by oblique reverse, up-to-the-east, dip-slip motion on the reactivated Cenozoic Cow Mountain–Mill Creek and Willow Creek Faults as material is advected through a restraining bend on the Maacama Fault. This study shows that reactivated bedrock faults increase the fault zone width and introduce fault surfaces that contribute a component of vertical deformation and uplift in major strike-slip fault zones. Deformation is accommodated on an interconnected network of new and reactivated faults that delineate a complex seismic hazard.

California↗

USGS2024

Use of a numerical groundwater-flow model and projected climate scenarios to simulate the effects of future climate conditions on base flow for reach 1 of the Washita River alluvial aquifer and Foss Reservoir storage, western Oklahoma

To better understand the relation between climate variability and future groundwater resources in reach 1 of the Washita River alluvial aquifer and Foss Reservoir in western Oklahoma, the U.S. Geological Survey, in cooperation with the Bureau of Reclamation, used a previously published numerical groundwater-flow model and climate-model data to investigate changes in base flow and reservoir storage by evaluating three scenarios. The three projected climate scenarios were (1) a central-tendency scenario, (2) a warmer/drier scenario, and (3) a less-warm/wetter scenario. To estimate future base flow and groundwater availability in western Oklahoma, specifically in reach 1 of the Washita River alluvial aquifer, downscaled climate-model data from 231 Coupled Model Intercomparison Project phase 5 (CMIP5) projections coupled with a previously published numerical groundwater-flow model were used to compare the effects of different climate scenarios on the aquifer. Changes in base flow and groundwater-level elevations during a 30-year baseline scenario (1985–2014) and the three 30-year projected climate scenarios (2050–79) under central-tendency, warmer/drier, and less-warm/wetter climatic conditions were assessed by using the calibrated model. In the simulations, the amount of base flow and reservoir storage declined in the central-tendency and warmer/drier scenarios compared to the amount of base flow and reservoir storage under historical climatic conditions (baseline scenario). Mean annual change in reservoir storage decreased from the baseline scenario the most in the warmer/drier scenario, followed by the central-tendency scenario, but increased in the less-warm/wetter scenario compared to the baseline scenario. At the end of the simulation period (2079), the largest magnitude differences in groundwater-level elevations in all three projected climate scenarios relative to the baseline scenario occurred upstream from Foss Reservoir. Results from incorporating downscaled climate projections into localized numerical groundwater-flow models can highlight potential future changes in and implications for groundwater resources and availability.

Oklahoma↗

USGS2024

Analysis of factors affecting plume remediation in a sole-source aquifer system, southeastern Nassau County, New York

Several plumes of dissolved, chlorinated solvents, including trichloroethylene, have been identified in a sole-source aquifer near the former Northrop Grumman Bethpage Facility and Naval Weapons Industrial Reserve Plant sites in southeastern Nassau County, New York. Past investigations have documented that the groundwater contamination originated from this industrial area and now extends to the south, in the direction of groundwater flow. The intermixed plumes are commonly referred to as the “Navy Grumman groundwater plume.” Detailed groundwater-flow modeling was needed for the New York State Department of Environmental Conservation (NYSDEC) to evaluate design options necessary for the construction, operation, optimization, maintenance, and monitoring of a groundwater extraction and treatment cleanup plan selected in a December 2019 Amended Record of Decision by the NYSDEC to comprehensively address these plumes. Consequently, the NYSDEC began a cooperative study with the U.S. Geological Survey in 2020 to better understand the local hydrogeologic framework using two independent approaches to characterize aquifer heterogeneity and update an existing regional groundwater-flow model to provide transient boundary conditions for new inset groundwater-flow models of the plume area. We developed these detailed inset models for the two independent aquifer characterizations using history-matching techniques coupled with a novel approach to risk-based management optimization of the remedial design. We also used the updated regional model to assess this optimized groundwater extraction and treatment design for potential saltwater intrusion. The ensembles of parameters resulting from history matching provided a platform with which to evaluate capture by water-supply and remedial wells using particle-tracking techniques. Using the ensemble to select a risk stance, we performed multiobjective optimization to identify various configurations of remedial pumping that are consistent with external constraints and that favor potentially competing objectives. Multiple solutions provide tradeoffs that NYSDEC can consider. In general, pumping redistribution may help to prevent further contamination migration downgradient. These and other study results are intended to support decisions for the remedial design focused on the local area encompassing the full extent of the Navy Grumman groundwater plume.

New York↗

USGS2024

New occurrences of the rare, REE minerals daqingshanite, törnebohmite, biraite, sahamalite, and ferriperbøeite from the Sheep Creek area, Montana, USA

Over 30 small, discontinuous, tabular carbonatite bodies are located in the Sheep Creek area, Ravalli County, southwest Montana. The age and origin of these REE-Nb-rich deposits are currently being investigated. The purpose of this paper is to document the occurrence of several rare minerals, including daqingshanite, törnebohmite, biraite, sahamalite, and ferriperbøeite, in two of the carbonatite bodies. These minerals are found in association with monazite, hydroxylbastnäsite, ferriallanite, calcite, dolomite, baryte, quartz, actinolite, apatite, celsian, and Sr-rich aragonite. Automated SEM-EDS was used to target the areas of interest in polished specimens for more detailed spot SEM-EDS and electron probe microanalysis. Raman spectra were also acquired for each of the rare minerals. The complex mineralogy of the Sheep Creek carbonatites is most likely due to several overlapping thermal events, including primary magmatic, overprinting hydrothermal, and supergene weathering stages. The rare minerals described in this study are believed to be hydrothermal and/or carbothermal in origin, although no estimates of temperature are available at this time.

Montana↗

USGS2024

Temporal, spatial, and chemical evolution of Quaternary high-silica rhyolites in the Mineral Mountains, Utah

The Mineral Mountains in southwestern Utah are a structurally controlled core complex at the confluence of the Colorado Plateau and the Basin and Range physiographic provinces. Aside from hosting Utah’s largest batholith, the Mineral Mountains host some of the State’s youngest high-silica rhyolites, which have been linked to a magma source that is presently being utilized as an enhanced geothermal system. The high-silica rhyolites take the form of effusive lavas and domes, and explosive products are rare. Previous K-Ar dating of these Pleistocene rhyolites placed eruptions between about 790 and 500 kilo-annum (ka) with contemporaneous basalts erupting in the valley to the east of the Mineral Mountains. Large uncertainties on these ages obscured the tempo of eruptions and thus hindered attempts to constrain the timescales of the petrogenetic processes that produced the rhyolites. In this study, we build on previous studies conducted in the 1970s and 1980s by using new geochronologic and geochemical data to investigate the temporal and spatial evolution of the youngest phase of volcanism in the Mineral Mountains. We identify two major eruptive periods, from approximately 850 to 750 ka and from approximately 590 to 480 ka. The older phase is characterized by the eruption of several basaltic lavas, two obsidian flows, and a series of coalescing porphyritic rhyolite domes. The younger phase included the eruption of six evolved high-silica rhyolite domes and one pyroclastic deposit, followed by the eruption of trachyandesite in the adjacent valley to the east. Whole-rock geochemical data indicate that the rhyolites can be divided into three chemical groups, with more evolved compositions erupting through time. The youngest rhyolites along the range crest have the lowest total iron and TiO 2 concentrations and the highest incompatible element concentrations, indicative of increasing differentiation with time and elevation. Improved precision on the eruption ages indicates a recurrence interval of approximately 20 thousand years. The eruptive flux for both periods of rhyolitic volcanism is about 0.01 cubic kilometers per thousand years, which is less than the magma resurgence flux rates for syn-caldera and post-caldera eruptions of the Valles Caldera and Yellowstone Caldera volcanic systems. Collectively, these geochemical, geochronological, and volumetric data may facilitate a better understanding of heat flux and the longevity of magmatic sources related to geothermal resources in similar small-volume, silicic systems.

Utah↗

USGS2024

Conceptualization and simulation of groundwater flow and groundwater availability in the Boone and Roubidoux aquifers in northeastern Oklahoma, 1980–2017

Oklahoma Groundwater Law (Oklahoma Statute § 82-1020.5) requires that the Oklahoma Water Resources Board conduct hydrologic investigations to determine the maximum annual yield for the State’s groundwater basins. The Boone and Roubidoux aquifers (also known as the Springfield Plateau aquifer and Ozark aquifer, respectively) are bedrock aquifers that extend from northeastern Oklahoma into Kansas, Arkansas, and Missouri. At present (2024), the Oklahoma Water Resources Board has yet to legally issue orders for the final determination of maximum annual yields for the Boone and Roubidoux aquifers. To support determination of a maximum annual yield, the U.S. Geological Survey, in cooperation with the Oklahoma Water Resources Board, developed a hydrogeologic framework, a conceptual groundwater-flow model, and a calibrated numerical groundwater-flow model for the Boone and Roubidoux aquifers. Three types of groundwater-availability scenarios were simulated by using the calibrated numerical model. These scenarios were used to (1) estimate equal-proportionate-share groundwater withdrawal rates (groundwater withdrawal applied equally over the aquifer), (2) quantify the potential effects of projected groundwater withdrawals on groundwater storage over a 50-year period, and (3) simulate the potential effects of a hypothetical 10-year drought. For the Boone aquifer, equal-proportionate-share groundwater withdrawal rates were 1.10, 0.98, and 0.96 acre-feet per acre per year for the 20-, 40-, and 50-year scenarios, respectively. For the Roubidoux aquifer, equal-proportionate-share groundwater withdrawal rates were 1.76, 1.34, and 1.25 acre-feet per acre per year for the 20-, 40-, and 50-year simulations, respectively. For the 50-year scenarios, stream seepage was minimally affected. Over the 10-year drought scenario, groundwater storage in the Boone and Roubidoux aquifers decreased by 660,451 acre-feet (6.7 percent) and 508,472 acre-feet (1.0 percent), respectively.

Oklahoma↗

USGS2024

State of science, gap analysis, and prioritization for southeastern United States water-quality impacts from coastal storms—Fiscal year 2023 program report to the Water Resources Mission Area from the Water Availability Impacts of Extreme Events Program—Hurricanes

Tropical cyclones (coastal storm events that include tropical depressions, tropical storms, and hurricanes) cause landscape-scale disturbances that can lead to impaired water quality and thus reduce water availability for use. Stakeholders and scientists at local and national scales have illustrated a need for understanding these risks to water quality. A regional and comprehensive understanding of the impacts of tropical storms and hurricanes on surface-water and groundwater quality—and thus water availability—is lacking for potentially impacted coastal and inland areas. As the U.S. Geological Survey considers development of tools to predict the extent to which water-quality impacts of hurricanes affect water availability, an assessment of the state of the science of hurricane impacts is needed, including a gap analysis and prioritization of data and science needs. This assessment focuses on the southeastern coastal States.

southeastern United States↗

USGS2024

Geologic input databases for the 2025 Puerto Rico – U.S. Virgin Islands National Seismic Hazard Model update: Crustal faults component

The last National Seismic Hazard Model (NSHM) for Puerto Rico and the U.S. Virgin Islands (PRVI) was published in 2003. In advance of the 2025 PRVI NSHM update, we created three geologic input databases to summarize new onshore and offshore fault source information in the northern Caribbean region between 62°–70° W and 16°–21° N. These databases, of fault sections, fault‐zone polygons, and geologic estimates of fault activity (fault‐slip rate and earthquake recurrence intervals) at specific sites, document updates to fault parameters used in prior seismic hazard models in PRVI. Fault sources were reviewed from published studies since 2003, which document substantial changes to the understanding of fault location, geometry, or activity. New fault section sources were added for features that meet the criteria of (1) length ≥7 km, (2) unequivocal evidence of recurrent tectonic Quaternary activity, and (3) documentation that is publicly available in a peer‐reviewed source. In addition, we revised several broad areal sources, such as the Mona and Anegada extensional zones. The 2003 model included three fault sections and two fault‐zone polygons (areal sources). These databases include 35 fault sections, 6 fault‐zone polygons, and 51 earthquake geology sites. To characterize fault activity rates, slip‐rate bins were assigned based on landscape expression and paleoseismic trench observations for faults without published slip‐rate sites. Additional fault sources were evaluated but not included in these databases due to a lack of published information about fault location, geometry, or recurrent Quaternary activity. The PRVI NSHM 2025 geologic input databases describe crustal faulting; the geometries and coupling of Puerto Rico subduction zone and Muertos Trough models are considered in a separate database. Updates to the fault sections, fault‐zone polygons, and earthquake geology databases can help inform the location and recurrence rate of damaging earthquakes in the PRVI NSHM implementation.

Puerto Rico, U.S. Virgin Islands↗

USGS2024

Correlation analysis of groundwater and hydrologic data, Kaloko-Honokōhau National Historical Park, Hawai‘i

Designated in 1978, Kaloko-Honokōhau National Historical Park is located on the west coast of the Island of Hawaiʻi. The Kaloko-Honokōhau National Historical Park encompasses about 1,200 acres of coastal land and nearshore ecosystems, which include wetlands, anchialine pools (landlocked bodies of brackish water with hydrologic connections to the ocean), fishponds, a fishtrap, and coral reefs. These nearshore ecosystems are dependent on groundwater discharge with a freshwater component and provide habitat for threatened and endangered, endemic species, such as the orangeblack Hawaiian damselfly ( Megalagrion xanthomelas ) and the Hawaiian coot (ʻAlae keʻokeʻo, Fulica alai ). The populations of these native species, however, are threatened because of habitat loss related to urban development and environmental changes. Kaloko-Honokōhau National Historical Park is within the Keauhou aquifer system and the North Kona District, which experienced a 52 percent resident-population increase between 2000 and 2020 and a 41 percent visitor increase between 2008 and 2019. To support the current water demand associated with this growing population, groundwater is the primary source of freshwater used in the North Kona District, with about 15 million gallons of groundwater withdrawn from the Keauhou aquifer system per day since 2009. With anticipated development, future (2015–35) groundwater withdrawal from the Keauhou aquifer system is projected to be about 55 percent greater than recent (2012–14) withdrawal. Because Kaloko-Honokōhau National Historical Park is located within a coastal aquifer, natural and human-induced changes can affect the quality and quantity of groundwater, which can threaten groundwater-dependent ecosystems. To improve understanding of recent groundwater conditions, the U.S. Geological Survey, in cooperation with the National Park Service, undertook this study to document correlations between hydrologic time-series datasets from sites in and near Kaloko-Honokōhau National Historical Park using the nonparametric (distribution-free) Kendall’s tau statistical test. For the statistical analyses, dependent variables representing the groundwater system include groundwater level, the groundwater-level difference between pairs of sites, and specific conductance, and independent variables include datasets of sea level, rainfall, and groundwater withdrawal. About 34 percent of the 140 non-time-lagged Kendall’s tau statistical tests evaluated in this report are statistically significant ( p -value ≤ 0.050) with generally weak (0.1 ≤ tau ≤ 0.2) to moderate (0.2 ≤ tau ≤ 0.3) correlations. Groundwater levels measured at monitoring sites have the strongest correlation with the multivariate El Niño–Southern Oscillation index and withdrawal from production wells at the nearby Kohanaiki Private Club Community. Specific conductance is not consistently and significantly correlated with the independent hydrologic variables investigated in this report. Because the relations between hydrologic variables are commonly not instantaneous, a second set of correlations was evaluated after applying a range of time lags to the independent variable datasets. Relative to the non-time-lagged case (the set of correlations that did not use time-lagged independent variables), some of the time-lagged independent variables improved correlations with some of the dependent variables. For a particular independent variable, similar time lags were expected between the independent variable and dependent variable at all four monitoring sites. However, different time lags among the four sites sometimes produced the strongest correlations. This study identified several correlations that are statistically significant and hydrologically plausible, but the correlations could indicate that multiple concurrent factors are controlling the observed groundwater-system response, which might be better addressed using multivariate analyses. This study only investigates bivariate correlations, which may not explain all the variance in the data. The correlations analyzed in this report are limited by the quantity of available hydrologic data in the area near Kaloko-Honokōhau National Historical Park and are based on 14 years of time-series data, which were aggregated to a relatively coarse monthly temporal resolution that represents the minimum resolution common to all datasets.

Hawaii↗

USGS2024

Pesticide occurrence in shallow groundwater in three regions of agricultural land use: Baldwin County, the Wiregrass region, and the Tennessee River valley region of Alabama, 2009–20

As part of a cooperative investigation between the U.S. Geological Survey and the Alabama Department of Agriculture and Industries, a network of 22 groundwater wells were sampled from 2014 through 2020 for about 230 pesticide and pesticide degradate compounds. Wells were located in three regions of intensive agricultural land use in Alabama: Baldwin County, the Wiregrass region, and the Tennessee River valley region. Metolachlor sulfonic acid, a degradate of the herbicide metolachlor, was the most frequently detected compound, occurring in about 70 percent of the samples. Three other compounds, metolachlor, atrazine, and 2-chloro-4-isopropylamino-6-amino-s-triazine, were also detected in more than half of the samples. Metolachlor and its degradates accounted for 33 of the 50 greatest compound concentrations study-wide, including the maximum pesticide concentration across all compounds (62,500 nanograms per liter). The frequency and magnitude of detections of many specific pesticide compounds varied among the three regions, but all detected pesticide concentrations were well below the U.S. Environmental Protection Agency maximum contaminant levels and applicable human health benchmarks. Sample results were combined with results of previous (2009–13) sampling to provide a continuous time-series of data for 2009–20. More than half of the 289 pesticide compounds analyzed during 2009–20 were not detected in any samples. Only four compounds were detected at great enough frequency throughout the 10 sampling years to evaluate patterns of change through time. Metolachlor and its degradate, metolachlor sulfonic acid, were frequently detected in all regions. Atrazine and its degradate, 2-chloro-4-isopropylamino-6-amino-s-triazine, were also detected in wells from all regions, but the variability and magnitude of concentrations were greatest in the Tennessee River valley region. No apparent temporal pattern in concentrations was found.

Alabama↗

USGS2024

Preliminary observations of the April 5th, 2024, Mw4.8 New Jersey earthquake

On 5 April 2024, 10:23 a.m. local time, a moment magnitude 4.8 earthquake struck Tewksbury Township, New Jersey, about 65 km west of New York City. Millions of people from Virginia to Maine and beyond felt the ground shaking, resulting in the largest number (>180,000) of U.S. Geological Survey (USGS) “Did You Feel It?” reports of any earthquake. A team deployed by the Geotechnical Extreme Events Reconnaissance Association and the National Institute of Standards and Technology documented structural and nonstructural damage, including substantial damage to a historic masonry building in Lebanon, New Jersey. The USGS National Earthquake Information Center reported a focal depth of about 5 km, consistent with a lack of signal in Interferometric Synthetic Aperture Radar data. The focal mechanism solution is strike slip with a substantial thrust component. Neither mechanism’s nodal plane is parallel to the primary northeast trend of geologic discontinuities and mapped faults in the region, including the Ramapo fault. However, many of the relocated aftershocks, for which locations were augmented by temporary seismic deployments, form a cluster that parallels the general northeast trend of the faults. The aftershocks lie near the Tewksbury fault, north of the Ramapo fault.

New Jersey↗

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

Beyond the wedge: Impact of tidal streams on salinization of groundwater in a coastal aquifer stressed by pumping and sea-level rise

Saltwater intrusion (SWI) is a well-studied phenomenon that threatens the freshwater supplies of coastal communities around the world. The development and advancement of numerical models has led to improved assessment of the risk of salinization. However, these studies often fail to include the impact of surface waters as potential sources of aquifer salinity and how they may impact SWI. Based on field-collected data, we developed a regional, variable-density groundwater model using SEAWAT for east Dover, Delaware. In this location, major users of groundwater from the surficial aquifer are the City of Dover and irrigation for agriculture. Our model includes salinized marshland and tidal streams, along with irrigation and municipal pumping wells. Model scenarios were run for 100 years and included changes in pumping rates and sea-level rise (SLR). We examined how these drivers of SWI affect the extent and location of salinization in the surficial aquifer by evaluating differences in chloride concentration near surface waters and the subsurface freshwater-saltwater interface. We found the presence of the marsh inverts the typical freshwater-saltwater wedge interface and that the edge of the interface did not migrate farther inland. Additionally, we found that tidal streams are the dominant pathways of SWI at our site with salinization from streams being exacerbated by SLR. Our results also show that spatial distribution of pumping affects both the magnitude and extent of salinization, with an increase in concentrated pumping leading to more intensive salinization than a more widely distributed increase of the same total pumping volume.

Delaware↗

USGS2024

Estimating groundwater level records using MOVE.1 and computing monthly percentiles from estimated groundwater records in Massachusetts

The U.S. Geological Survey, in cooperation with the Massachusetts Department of Environmental Protection, performed record extensions on groundwater levels at select wells using the Maintenance of Variance Extension type 1 (MOVE.1) method. The groundwater levels estimated from these record extensions were used to compute monthly percentiles to improve future determinations of a groundwater index. In Massachusetts, 27 of 29 short-record study wells with continuous groundwater levels between 0.8 and 8.1 years were suitable for record extensions; 37 long-record index wells were used to extend the groundwater level records at the study wells. The index well selected to pair with a study well was chosen based on Pearson correlation coefficient values; cross-correlation between the two wells; geologic and topographic similarity; and smallest distance spanning the wells. Each study well and its corresponding index well have 1 or more years of concurrent, overlapping data; a Pearson correlation coefficient that exceeded a threshold value of 0.8; and a similar aquifer type and hydrologic characteristics. Of the 29 study wells, 2 showed poor correlations with all index wells and were not considered for record extensions. Performance metrics used to assess the accuracy of the MOVE.1 models indicated that most models provided reasonable estimates of groundwater levels. Root mean square error values ranged from 0.097 to 2.292 feet, with a median of 0.536 foot. Nash-Sutcliffe efficiency coefficient values ranged from 0.623 to 0.996, with a median value of 0.759. Generally, study wells in close geographical proximity to their index well resulted in stronger model performance. The average length of groundwater level records was extended by 14.1 years to a new average of 18.1 years. The estimated groundwater level records from the MOVE.1 models resulted in an increase in the range of highest and lowest groundwater levels at 23 of 27 wells. The increase in range of groundwater levels was between 0.08 to 7.95 feet. Monthly percentiles for State drought indices were computed from the estimated MOVE.1 records and observed records through December 31, 2021. Percentiles computed from estimated records show an average groundwater level about 1.0 foot lower than observed data at the 2d percentile and 0.1 foot lower at the 30th percentile.

Massachusetts↗

USGS2024

Stream nitrate dynamics driven primarily by discharge and watershed physical and soil characteristics at intensively monitored sites: Insights from deep learning

We developed a suite of models using deep learning to make hindcast predictions of the 7‐day average backward‐looking nitrate concentration at 46 predominantly agricultural sites across the midwestern and eastern United States. The models used daily observations of discharge and meteorological variables and watershed attributes describing anthropogenic modification to hydrology, nitrogen application, climate, groundwater, land use, watershed physiographic attributes, and soils. Across all sites, discharge and watershed soil and physiographic attributes showed a strong influence on model performance. Analysis of drivers across sites revealed considerable regional differences related to controlling processes such as groundwater contributions. We tested several ways to pool data across sites to develop accurate models and make the most effective use of available data. Single‐site models, in which models are trained and tested at a single location, showed generally strong predictive performance (median Kling‐Gupta Efficiency = 0.66), and accuracy at poorly performing sites could be improved by grouping sites with similar characteristics. Developing a single model for all sites reduced performance at several locations with distinct characteristics, suggesting that there is a threshold of dissimilarity beyond which more data does not improve the model. While many deep learning studies have shown that national or even global models can outperform local models, it is not clear that this is true for water quality constituents. This study demonstrates how data can be combined effectively, using deep learning to develop accurate and interpretable models of instream nitrate at sites where varying processes are responsible for changes in nitrate concentration.

Water Resources Research↗

USGS2024

Groundwater quality near the Placerita Oil Field, California, 2018

Groundwater-quality data and potential fluid-migration pathways near the Placerita Oil Field in Los Angeles County, California, were examined by the U.S. Geological Survey to determine if oil-field fluids (water and gas from oil-producing and non-producing zones) have mixed with groundwater resources. Six of the 13 new groundwater samples collected for this study contained petroleum hydrocarbons, thermogenic gas, inorganic chemical signatures, and (or) isotopic values consistent with potential mixing with fluids from hydrocarbon-bearing formations. For historical groundwater samples, benzene was the most detected petroleum hydrocarbon. The historical groundwater samples with a benzene concentration greater than 0.5 micrograms per liter were from environmental monitoring wells at industrial or commercial facilities unrelated to oil and gas development that, in many cases, have identified soil or groundwater contamination and were not typically analyzed for other constituents that could provide additional lines of evidence for potential mixing with oil-field fluids. Methane was not detected in any of the 12 historical samples with a reported measurement. Reviewing historical data revealed factors that could potentially adversely affect groundwater quality in the study area. These factors include modified hydraulic gradients caused by large volumes of water extracted from the main production area and reinjected downgradient into nonproducing zones, well-barrier failures in wells constructed in the northern part of the oil field before the 1970s, well-barrier failures in produced-water disposal wells downgradient from the main production area, and naturally occurring hydrocarbons at shallow intervals. The groundwater samples most geochemically similar to samples from hydrocarbon-bearing formations were in areas where hydrocarbons are naturally occurring at shallow intervals and where oil development is at shallow depths. Additional data for hydraulic heads, water quality, and formation temperatures at multiple depths in areas with large injection volumes and well-integrity issues are needed to evaluate whether those factors have contributed to mixing between fluids from oil-producing or injection formations and groundwater resources.

California↗

USGS2024

Simulated mean monthly groundwater-transported nitrogen loads in watersheds on the north shore of Long Island Sound, 1993–2022

Elevated nitrogen loads are pervasive in the Long Island Sound, an estuary that receives freshwater and nutrients from both surface-water and groundwater discharge. Surface-water nitrogen loads to the Long Island Sound are relatively well characterized, but less is known about groundwater-transported nitrogen loads. Prior work on the northern shore of Long Island Sound (Connecticut and areas of New York and Rhode Island) suggested that groundwater travel times are relatively short (median less than 2 years) and that decade-long nutrient legacies are not widespread. Because the travel times are short, groundwater flow and nutrient loads likely vary substantially between months. In the current study, the U.S. Geological Survey, in cooperation with the U.S. Environmental Protection Agency’s Long Island Sound Study and the Connecticut Department of Energy and Environmental Protection, developed a set of models to better characterize spatial and temporal patterns of groundwater-transported nitrogen loading from atmospheric deposition, septic systems, and fertilizers within the study area. The models provide an estimate, with uncertainty, of groundwater-transported nitrogen loads in the study area, filling a key gap in the nitrogen budget for Long Island Sound. The models also highlight the spatial and temporal variation in nitrogen loading throughout the study area. The modeling workflow involved four models. (1) A soil-water-balance model was developed by using the Soil-Water-Balance software to simulate groundwater recharge across the study area for water years 2005 through 2022. The simulated mean monthly recharge from the soil-water-balance model was used as input into a groundwater-flow model. (2) The groundwater-flow model was developed by using the MODFLOW 6 software and data for water years 1993 through 2022 and simulates average monthly hydrologic conditions. The groundwater-flow model was calibrated by using the Iterative Ensemble Smoother method within the PEST++ software. The Iterative Ensemble Smoother method generates an ensemble of sets of parameter values, with each set producing reasonable simulated hydrologic parameter values. (3) An ensemble of MODPATH particle-tracking simulations were run to generate particle flow paths and travel times, with each simulation using a different set of the flow model parameters. (4) A nitrogen load model uses the MODPATH simulation outputs to track nitrogen from the land surface through multiple attenuation zones until it discharges into fresh or saline surface water. As with the groundwater-flow model, the nitrogen model simulated average monthly groundwater-transported nitrogen loads for water years 1993 through 2022. One novel aspect of the nitrogen load model is that the nitrogen attenuation parameters were calibrated to observed nitrogen loads. Across the ensemble of simulated nitrogen loads, the median study-area-wide monthly simulated nitrogen loads from the aquifer to Long Island Sound throughout the year ranged from 900 to 18,600 kilograms of nitrogen per day, with a median load of 5,100 kilograms of nitrogen per day. The simulated loads were based on average monthly conditions for water years 1993 through 2022. Loads were highest during the winter and early spring and lowest during the late summer. However, simulated travel times for groundwater and nitrogen loads discharged to Long Island Sound during summer were longer than travel times for groundwater and loads discharged during the winter, indicating that, on average, groundwater discharged during summer traveled along different, and longer, flow paths, than groundwater discharged during winter. This indicates that summer loads would respond more slowly to changes in nitrogen inputs at the water table than winter loads. Over the entire study area, approximately 15 percent of the simulated load is from atmospheric deposition sources, 30 to 40 percent is from fertilizer, and 50 to 60 percent is from septic systems. The final analysis of the study involved simulating the change in groundwater-transported nitrogen load in response to upgrading septic systems or reducing fertilizing inputs to areas of turf grass. Both management interventions reduced the groundwater-transported nitrogen load, and reductions were greater in areas with greater loads from septic systems or turf-grass fertilizers. The delay between management actions and substantial reductions in groundwater-transported nitrogen loads varied seasonally; loads during the late summer months remained elevated longer than the winter loads.

Connecticut, Rhode Island↗