USGS Science⌕ Search

SEARCH · USGS Science

Results for “Frontiers in Enviornmental Science”

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,605 records · Page 28Linked to original sources
USGS2026

Deep groundwater total dissolved solids mapping in the Dakota Group, Williston Basin, USA

Growing concern about the quantity of available freshwater around the world has led to interest in surveying groundwater total dissolved solids (TDS) below water well depths. Deep TDS has not been systematically mapped, and there is much to learn about the distribution and controls on deeper groundwater. In sedimentary basins across the United States, groundwater resources often overlie hydrocarbon resources, providing an opportunity to use borehole geophysical data collected for hydrocarbons to characterize groundwater and pore space resources. This study adapts a recently developed subsurface geostatistical and geophysical modeling approach to continuously map groundwater TDS, porosity, and temperature in the Dakota Group of the Williston Basin—an undercharacterized regional aquifer system overlying deeper hydrocarbon reservoirs. Groundwater TDS in the Dakota Group ranges from approximately 4800 to 26,900 mg/L. TDS patterns are stratified with higher TDS in the lower and upper Dakota Group, and relatively lower TDS in the middle Dakota Group. The lower TDS in the middle zone may represent a preferential regional flow path for lower-TDS meteoric recharge from the west. The alternating pattern of TDS may also be evidence of higher-TDS inflows into the Dakota Group from underlying and potentially from overlying aquifers. Porosity is lower near the center of the Williston Basin and tends to be higher to the east, which may be related to grain size distributions. The new regional TDS and porosity modeling serves as a quantitative reference for water users and provides supporting evidence for hypotheses on Dakota Group recharge.

Montana, North Dakota, South Dakota↗
USGS2025

Climate change and future water availability in the United States

The steady rise in global temperature as a result of human activity is causing changes in Earth’s water cycle. The balance of water stored within and moving between vapor, liquid, and frozen states in the water cycle is shifting, with consequences for water availability that include increases in drought, fire weather, flooding, and heavy precipitation, as well as cryosphere decline and sea-level rise. In this chapter of the U.S. Geological Survey Integrated Water Availability Assessment—2010–20, we provide an overview of climate-change observations and projections from Earth-system model simulations that relate to future water availability, from global and national climate assessments and from the published literature. Effects of climate change on primary water-cycle components are discussed in context of how global-scale hydroclimate drivers influence regional processes within the United States. Understanding the major climate drivers impacting the water cycle is crucial to predicting future changes in water availability and developing adaptation strategies to ensure human and ecosystem water supplies. First, we provide background information on the water cycle, the climate-model ensemble simulations developed to produce projections based on warming scenarios, and attribution and certainty levels. Tipping points, self-reinforcing feedbacks, cascading effects, and compound extremes are introduced. The framework of climatic impact drivers (CIDs) outlined in the Intergovernmental Panel on Climate Change Sixth Assessment Report (IPCC AR6) is used to show primary drivers of physical change to the water cycle and to understand and predict changes in future water availability. Specific climate-change related observations and projections are discussed for water cycle components of precipitation, evapotranspiration, soil moisture, streamflow, lakes and wetlands, ice and snow, and groundwater, as well as their implications for future water availability for humans and ecosystems. The chapter concludes with a synthesis discussion of three examples of complex regional-scale hydroclimate processes that influence water availability for populations in the United States, including (1) mountain and coastal precipitation, (2) aridification and drought, and (3) the influence of forest-cover change on terrestrial water-vapor recycling.

Professional Paper↗
USGS2025

Data gap analysis for estimation of agricultural return flows in the Upper Gunnison River Basin, Colorado

The Gunnison River and many tributaries in the Upper Gunnison River Basin provide water to irrigate agricultural crops. The application of irrigation water can recharge some aquifers locally by water percolating below the root zone and eventually flowing back to the stream or river through the subsurface. Diverting surface water for irrigation reduces streamflow during the irrigation season but can provide temporary storage of water and supplement streamflow after the snowmelt runoff season. Understanding the timing and quantity of agricultural return flows could help resource managers make informed decisions and adapt to potential changes in water management and availability that could affect irrigation practices. In 2024, the U.S. Geological Survey, in cooperation with the Upper Gunnison River Water Conservancy District, began a study to characterize agricultural return flows in the Upper Gunnison River Basin by using endmember mixing analysis and developing a groundwater model. Both approaches require data from multiple sources, but data gaps exist in the East River study reach and other reaches of interest (Ohio Creek, Tomichi Creek, and Cochetopa Creek). The East River Basin, which is the initial focus of the study, has fewer data gaps than the other basins. Data gaps could be addressed by installing additional surface water and groundwater monitoring sites, making regular streamflow measurements on tributaries, and completing tests to characterize local aquifer properties.

Colorado↗
USGS2025

2024 Surprise Inlet landslides: Insights from a prototype landslide‐triggered tsunami monitoring system in Prince William Sound, Alaska

Alaska's coastal communities face growing landslide hazards owing to glacier retreat and extreme weather intensified by the warming climate, yet hazard monitoring remains challenging. As part of ongoing experimental monitoring in Prince William Sound, we detected three large landslides (0.5–2.3 M m 3 ) at Surprise Inlet on 20 September 2024, within the span of an hour. These events were identified in near real-time through seismic data and later confirmed using satellite imagery, tidal records, and infrasound. The landslides generated a modest tsunami, and a 4 cm wave was recorded by a tide gauge 18 km away, marking the first recorded landslide to reach water since monitoring began in this region in 2021. Here, we examine the detection and interpretation of these landslides using multiple data sources and modeling. We demonstrate the effectiveness of this regional seismic monitoring system and show how complementary instrumentation, where available, can enhance detection capabilities.

Alaska↗
USGS2025

Timescales of surface faulting preservation in low-strain intraplate regions from landscape evolution modeling and the geomorphic and historical record

Large surface-rupturing intraplate earthquakes in stable continental regions (SCRs) are uncommon globally and have recurrence intervals of thousands to hundreds of thousands of years based on the paleoseismic and geomorphic record, challenging accurate active fault identification in these regions. To constrain the timescales of preservation for scarps created by surface ruptures from dip-slip earthquakes, we use a two-dimensional scarp diffusion model for typical intraplate settings and explore which parameters influence fault scarp preservation. These parameters include the coseismic vertical surface offset, the recurrence interval of similar magnitude earthquakes, diffusivity (as a proxy for mean annual precipitation rate), and the erodibility of the surficial material. We constrain parameter ranges from a compilation of historical surface ruptures in intraplate settings in a variety of climates, including the Central and Eastern United States, Australia, Europe, Central Asia (Mongolia, China), India, and West Africa. The timescales of scarp preservation from landscape evolution modeling agree well with observations of scarp preservation in low-strain SCR and intraplate tectonic settings, with some notable exceptions for Australian scarps. We find that the erodibility of the surficial material and earthquake recurrence interval have a stronger effect on the timescales of scarp preservation than diffusivity or coseismic vertical surface offset. Our model results may aid in identifying and characterizing subtle, slow-moving active faults in low-strain SCR and intraplate tectonic settings for different tectonic, geomorphic, and climatic characteristics. Accurate fault locations and characterization from the landscape record has implications for both probabilistic seismic and fault displacement hazard analyses.

Journal of Geophysical Research Solid Earth↗
USGS2024

A methodology to estimate CO2 and energy gas storage resources in depleted conventional gas reservoirs

Depleted hydrocarbon reservoirs are subsurface geological structures capable of sequestering vast quantities of carbon dioxide (CO 2 ) as well as storing other energy gases for later usage, such as natural gas, and potentially hydrogen (H 2 ). Here we outline a methodology to quantify multi-gas storage resources in depleted conventional gas reservoirs for usage in assessments by the United States Geological Survey (USGS) at the scale of sedimentary basins. The methodology consists first of quantifying accessible pore volume in a depleted reservoir for natural gas storage using up to three equations. Input data are derived from commonly reported or estimated reservoir parameters and natural gas production volumes, and equations may be combined in linear models to improve pore volume estimates. Storage estimates from these equations are tested and validated for 31 reservoirs in the Michigan Basin Province, USA that were previously converted to underground gas storage facilities and have known (federally reported) natural gas storage capacities. Secondly, natural gas storage capacities can be transformed via fluid substitution calculations to estimate the storage resources for non-native fluids, applied here for, CO 2 , H 2 , and methane-H 2 blends, accounting for molecule-specific deviations from ideal gas behavior at reservoir pressures and temperatures as well as differing storage efficiencies. Importantly, the storage of non-native fluids may not be appropriate in all depleted gas reservoir settings due to potential risks like leakage, in particular in the case of H 2 storage, requiring additional knowledge of caprock sealing capacity. Given this caveat, we demonstrate the fluid substitution method for natural gas reservoirs of the Northern Niagaran Reef and Southern Niagaran Reef USGS plays in the Michigan Basin Province, as these trends of Silurian pinnacle reefs are capped with tight-sealing evaporite facies. The deterministic equations outlined from this methodology can be incorporated into future probabilistic USGS gas storage assessments for CO 2 , H 2 , and natural gas resources in the United States.

Michigan↗
USGS2024

Models no not provide proof: An example of model ambiguity and application of isotopic data in a mine pit lake

Geochemical and hydrologic models of pit lakes are commonly used in environmental regulatory decisions to predict future water quality and hydrologic conditions and to understand existing pit lakes. Models may be used to quantify sulfide oxidation, predict thermal/chemical stratification and mixing, and better understand connections between pit lakes and aquifers. One concern related to the hydrologic character of pit lakes is if they are terminal (a groundwater sink with no outflow) or flowthrough (both receiving groundwater inflow and discharging to groundwater). This question was pertinent to the Liberty pit lake, a small acidic pit lake formed in a former Cu deposit in south-central Nevada where potentiometric and geochemical data potentially indicate pit-lake outflow. Potential discharge to groundwater from the pit lake was evaluated using a water-balance model, but uncertainty in hydraulic parameters led to ambiguity in the hydrologic character. Stable isotopes of water were then sampled from the pit lake and adjacent groundwater wells, which unambiguously indicated the lack of an evaporative signature in downgradient groundwater because the groundwater did not plot on a hypothetical mixing line between evaporated pit lake water and observed meteoric recharge. This methodology provided a more effective and more data-driven approach for understanding pit-lake hydrology. Although predictive models are required to quantify reasonable bounds on future conditions, many models contain substantial uncertainty and are not well suited in some environments. Datasets that provide more clear lines of evidence could be collected from existing pit lakes whenever possible to inform water-rock interaction, limnological behavior, and connectivity to adjacent groundwater.

Nevada↗
USGS2026

Characterizing precariously balanced rocks (PBRs) in the eastern United States for estimating maximum past earthquake ground motions

Fragile geologic features (FGFs) can be used to estimate an upper bound to the strength of earthquake ground shaking in the past thousands to tens of thousands of years. FGFs could be especially useful in the eastern United States, where few active faults are available for paleoseismic studies and where low seismic attenuation means that an earthquake could damage FGFs at large distances. This article describes a pilot study of one form of FGF, precariously balanced rocks (PBRs), in the Blue Ridge of Virginia, ∼160 km from the central Virginia seismic zone (CVSZ) and ∼100 km from Virginia’s Giles County seismic zone (GCSZ). The main goal is to develop and refine methods for studying PBRs in the eastern United States. Photogrammetric 3D models of the rocks were used to locate their centers of mass and delineate their basal contacts, with the latter corroborated by wrapping a wire around the bases to verify the circumferences. Relations between hill geometry and amplification were used to estimate topographic amplification. Seismometers on the PBRs provided the frequencies, amplitudes, and directions of rocking after giving the rocks gentle pushes. The levels of ground shaking likely to topple the rocks were compared with the median hazard curves from the U.S. Geological Survey’s National Seismic Hazard Model and were used to map magnitudes of earthquakes in the nearby region likely to topple the rocks. Results indicate that an M w ≥ 7.5 earthquake is needed in the CVSZ or near the 2020 Sparta, North Carolina, earthquake to topple the two most fragile PBRs and that an earthquake of ~ M w 7–7.5 in the GCSZ is needed. The fragilities and estimated ages of the PBRs (20–40 ka) are thus consistent with the absence of such events in the historical record in the nearby region.

eastern United States↗
USGS2000

Natural attenuation of chlorinated volatile organic compounds in ground water at Area 6, Naval Air Station Whidbey Island, Washington

Natural attenuation is a viable alternative to pump and treat for meeting remediation objectives in the vicinity of the southern contaminant plume. The combination of historically low contaminant concentrations in ground water, a landfill cap that limits source area contributions, favorable conditions for degradation of VC, and a relatively long downgradient distance to potential receptors are all favorable for natural attenuation as a remediation alternative. Natural attenuation could effectively meet all but one remediation goal that extraction wells PW- 2, PW-4, PW-6, PW-7, PW-8 and PW-9 are currently being employed to meet. The goal of preventing migration of all VC across the site boundary could not be met by natural attenuation. Some VC would migrate south of the Navy boundary, but the potential for subsequent VC mineralization downgradient of the base and the existing institutional controls would result in minimal additional risk from using natural attenuation. In the western contaminant plume natural attenuation is not currently a viable alternative to pump and treat for meeting remediation objectives. There is a possibility that rates for reductive dechlorination of TCE and TCA could increase substantially if the plume was allowed to migrate beneath the Oak Harbor landfill, but there are not enough data to be certain of such an increase. TCA and TCE concentrations at the leading edge of the western plume need to be reduced to at least 25 and 40 µg/L, respectively, to be protective of potential downgradient receptors. Source area TCA and TCE concentrations have decreased substantially over the past ten years, and the extraction wells PW-3 and PW-5 in particular are removing a significant mass of contaminants from ground water, so natural attenuation may be a viable alternative for the western plume in the future. There would be some possible side benefits of using natural attenuation as an alternative to pump and treat in the southern contaminant plume. The first would be that the resulting decrease in the amount of treated water that would need to be recharged in the swale north of the landfill would result in less off-base migration of contamination across the western site boundary. The second benefit would be a substantial reduction in the amount of dissolved iron and manganese being extracted from the shallow aquifer and run through the treatment system. Removing that source of operation and maintenance problems would result in more effective containment and removal of contamination in the western contaminant plume. The most critical data gap identified in this evaluation is the paucity of contaminant chemistry information downgradient of the Navy boundary in the vicinity of the southern contaminant plume. Without such data, the behavior of the plume and the protectiveness of natural attenuation to downgradient receptors cannot be verified, and field attenuation rates for VC cannot be determined directly. The existing long-term monitoring plan would need to be reviewed and revised if natural attenuation is selected as a remedy for the contamination in the southern plume. In particular, additional performance monitoring wells may be required downgradient of the property.

Washington↗
USGS2024

Twentieth century extreme precipitation detected in a high-resolution, coastal lake-sediment record from California

California faces increasing economic and societal risks from extreme precipitation and flooding associated with atmospheric rivers (ARs) under projected twenty-first century climate warming. Lake sediments can retain signals of past extreme precipitation events, allowing reconstructions beyond the period of instrumental records. Here, we calibrate AR-related extreme precipitation from the last century to proxy data from lake sediments collected in the latitudinal zone of the highest frequency landfall for modern ARs in California. Excursions in erosional proxy data (Ti/Al) are positively and significantly correlated (r median = 0.45, p median = 0.04) with modern records of integrated vapor transport (IVT, kg m −1 s −1 ), a key metric of AR intensity, using correlations that incorporate age-model uncertainty. Despite the land-use change near the study site, the data suggest intense and long-lasting AR storms are identifiable in this sedimentary record. These results allow conservative inferences concerning past extreme hydrology at this site.

California↗
USGS1976

Effect of deicing chemicals on the hydrologic environment in Massachusetts; evaluation of surface-water data collection

The objective of this study is to develop predictive relationships that can be used to describe the impact of highway deicing salts on the hydrologic environment. A method, presented in an earlier report, for estimating yearly mean chloride concentrations from estimated or actual runoff and salt-application data was applied to the period 1971 through1974 and the results compared with yearly mean chloride concentrations computed from records of specific conductance. Estimates were found to be from minus 63 percent to plus 56 percent in error when the error was expressed as the difference between estimated and computed values as a percentage of the computed values. Chloride concentrations used in this study are computed from specific conductance/chloride concentration relationships and records of specific conductance. Preliminary results from graphic analyses and least-squares regression analysis show that relationships between measured values of specific conductance and chloride concentration have correlation coefficients ranging from 0.30 to 0.97. Analysis of streamflow for all major dissolved constituents is recommended with the purpose of attempting to describe the variations in the specific conductance/chloride concentration relationships.

Massachusetts↗
USGS2024

A new water temperature modeling approach to predict thermal habitat suitability for nonnative cichlids in Florida rivers

As global temperatures increase, the spatiotemporal arrangement of thermal habitats in Florida rivers may shift, creating the potential for greater dispersal and establishment of nonnative tropical freshwater fishes. To understand how water temperature changes may affect the spatial distribution of these nonnative species, more effective water temperature prediction models are necessary. Currently, most models employ either a generalized air–water temperature relationship or require expensive and complicated tools to measure hydrometeorological factors (e.g. groundwater input). Thus, we developed a novel modeling approach that is accurate, accessible, and cost-effective in allowing fisheries managers to project water temperatures in rivers across Central and North Florida. To characterize the potential for nonnative fishes to spread northward, we evaluated two hardy and abundant species currently found primarily in South Florida: Mayan Cichlid ( Mayaheros urophthalmus ) and Oscar ( Astronotus ocellatus ). Our results show an increase in thermally suitable winter days for both species in 10 of 11 rivers studied, consistent with predicted water temperature warming under 16 climate-change scenarios spanning different levels of air temperature warming (+1 °C, +2 °C, +3 °C, +4 °C) and precipitation/groundwater thermal sensitivity (0, 0.33, 0.66, 1). Considering resource limitations, fisheries managers can use our water temperature modeling approach to predict effects of climate change on Mayan Cichlid and Oscar survival, growth, and dispersal and take actions to manage potential northward movement of these species.

Florida↗
USGS2025

Tracing metal sources and groundwater flow paths in the Upper Animas River watershed using rare earth elements and stable isotopes

Groundwater flow paths and processes that govern metal mobility and transport are difficult to characterize in mountainous bedrock watersheds. Despite the difficulty in holistic characterization, conceptual understanding of subsurface hydrologic and geochemical processes is key to developing remediation plans for locations affected by acid mine drainage, such as the Upper Animas River watershed in southwestern Colorado, USA. Stable isotopes of water and rare earth elements were utilized to evaluate groundwater flow and metal sources within this complex catchment. Stable isotope samples collected from draining mine adits and springs display systematic spatial variation wherein sample sites at higher elevations have greater seasonal variability than sites at lower elevations. The Upper Cement Creek watershed, where multiple draining mines are present, displays the lowest seasonal variation in stable isotopic signatures, potentially indicating the presence of a large, well-mixed volume of groundwater storage or interbasin groundwater flow. Rare earth elements display statistically significant variation between different alteration styles in the catchment. Overprinting of regional propylitic alteration is evident based on enrichment of middle rare earth elements in acidic springs and mines that are not spatially associated with surficial exposures of acid generating alteration styles. Europium anomaly and middle rare earth enrichment signatures from two flooded mine tunnels on opposite sides of a watershed divide indicate connections to the same subsurface flooded mine workings.

Colorado↗
USGS2025

Linking permafrost to the abundance, biomass, and energy density of fish in Arctic headwater streams

Permafrost thaw alters groundwater flow, river hydrology, stream-catchment interactions, and the availability of carbon and nutrients in headwater streams. The impact of permafrost on watershed hydrology and biogeochemistry of headwater streams has been demonstrated, but there is little understanding of how permafrost influences fish in these ecosystems. We examined relations among permafrost characteristics, the resulting changes in water temperature, stream hydrology (e.g., discharge flashiness), and macroinvertebrates, with the abundance, biomass, and energy density of juvenile Dolly Varden ( Salvelinus malma ) and Arctic Grayling ( Thymallus arcticus ) across 10 headwater streams in northwestern Alaska. Macroinvertebrate density was driven by concentrations of dissolved carbon and nutrients supporting stream food webs. Dolly Varden abundance was primarily related to water temperature with fewer fish in warmer streams, whereas Dolly Varden energy density decreased with the flashiness of the headwater streams. Dolly Varden biomass was related to both temperature and bottom-up food web effects. The energy density of Arctic Grayling decreased with warmer temperatures and discharge flashiness. These relations demonstrate the importance of terrestrial–aquatic connections in permafrost landscapes and indicate the complexity of landscape effects on fish. Because permafrost thaw is one of the most impactful changes occurring as the Arctic warms, an improved understanding of how stream temperature, hydrology, and bottom-up food web processes influence fish populations can aid forecasting of future conditions across the Arctic.

Alaska↗
USGS2024

Predicted exposure of communities in southeastern United States to climate-related coastal hazards

A rigorous analysis of 21st Century multi-hazard exposure for U.S. Southeast Atlantic coastal communities indicates that up to 70% of residents will be exposed daily to shallow and emerging groundwater by ~2100, 15 times higher than from surficial flooding alone. This threat further exacerbates other coastal stressors, such as flooding, subsidence, and beach erosion, that impact these vulnerable coastal communities and ecosystems under expected future climate change scenarios.

Florida, Georgia, North Carolina, South Carolina, ↗
USGS2025

High-pass corner frequency selection and review tool for use in ground-motion processing

Raw seismological waveform data contain noise from the instrument’s surroundings and the instrument itself that can dominate recordings at low and high frequencies. To use these data in ground‐motion modeling, the effects of noise on the signals must be reduced and the signals’ usable frequency range identified. We present automated procedures to efficiently reduce low‐frequency noise that are implemented in the software package gmprocess. These procedures check for, and as needed remove, low‐frequency artifacts in the displacement record using polynomial fits, which can be used in combination with existing signal‐to‐noise ratio (SNR)‐based corner‐frequency selection procedures. The automated selections are then efficiently verified and refined using a graphical user interface (GUI) that plots relevant ground‐motion time series and spectra and tracks modifications to signal processing parameters. We demonstrate these procedures using recordings from the 2020 M 5.1 Sparta, North Carolina, and the 2013 M 4.7 southern Ontario earthquakes. Data processed with the SNR‐only and polynomial criteria for these events contain displacement artifacts in 37% and 23% of processed traces, respectively. Records with remaining artifacts are corrected manually using the GUI. These processing steps illustrate the workflow for efficient data processing with quality control.

Seismological Research Letters↗
USGS2026

Multi-objective optimization of a hydro-economic model in an over-allocated agricultural basin

Groundwater depletion for agricultural irrigation poses significant environmental and economic challenges. This study introduces a proof-of-concept that combines hydro-economic modeling, scenario-based modeling, and multi-objective optimization to manage pumping curtailment in an over-allocated basin in the western United States. Three optimization scenarios were evaluated, each offering different degrees of management flexibility. Results reveal that scenarios with finer spatial resolution achieved greater environmental benefits per unit profit loss. Additionally, strategies allowing fractional reductions in curtailed wells–rather than complete shutdowns based on water rights seniority–substantially improved efficiency, highlighting the value of increased decision-making flexibility. Although scenario testing can aid stakeholder engagement and strategy exploration, multi-objective optimization provides a systematic framework to quantify tradeoffs between competing objectives. This combined approach demonstrates promise for building consensus and supporting the design of sustainable water management strategies that balance agricultural livelihoods with ecosystem preservation.

Oregon↗
USGS1999

Water-quality assessment of the New England coastal basins in Maine, Massachusetts, New Hampshire, and Rhode Island: Environmental settings and implications for water quality and aquatic biota

The New England Coastal Basins in Maine, Massachusetts, New Hampshire, and Rhode Island constitute one of 59 study units selected for water-quality assessment as part of the U.S. Geological Survey's National Water-Quality Assessment (NAWQA) program. England Coastal Basins study unit encompasses the fresh surface waters and ground waters in a 23,000 square-mile area that drains to the Atlantic Ocean. Major basins include those of the Kennebec, Androscoggin, Saco, Merrimack, Charles, Blackstone, Taunton, and Pawcatuck Rivers. Defining the environmental setting of the study unit is the first step in designing and conducting a multi-disciplinary regional water-quality assessment. The report describes the natural and human factors that affect water quality in the basins and includes descriptions of the physiography, climate, geology, soils, surface- and ground-water hydrology, land use, and the aquatic ecosystem. Although surface-water quality has greatly improved over the past 30 years as a result of improved wastewater treatment at municipal and industrial wastewater facilities, a number of water-quality problems remain. Industrial and municipal wastewater discharges, combined sewer overflows, hydrologic modifications from dams and water diversions, and runoff from urban land use are the major causes of water-quality degradation in 1998. The most frequently detected contaminants in ground water in the study area are volatile organic compounds, petroleum-related products, nitrates, and chloride and sodium. Sources of these contaminants include leaking storage tanks, accidental spills, landfills, road salting, and septic systems and lagoons. Elevated concentrations of mercury are found in fish tissue from streams and lakes throughout the study area.

Maine, Massachusetts, New Hampshire, Rhode Island↗