Search USGS⌕ Search

Geology topics

Joshua C. Woda

Publications and source records attributed to Joshua C. Woda.

16 recordsLinked to original sources

Hydrogeology and water quality of the unconsolidated aquifer system in the lower Fall Creek valley, Tompkins County, New York

The U.S. Geological Survey, in cooperation with Tompkins County and the Town of Dryden, conducted a hydrogeologic investigation of the lower Fall Creek unconsolidated aquifer system in Tompkins County, New York. This study is part of an ongoing effort to characterize major valley-fill aquifers in the region to support sustainable groundwater management. The lower Fall Creek valley, encompassing the Village of Freeville, the unincorporated communities of Etna, and the hamlets of McLean and Varna, has experienced moderate development and relies heavily on groundwater for residential, agricultural, and commercial use. The study area spans approximately 51 square miles along the boundary of the Allegheny Plateau and North Rim of the Allegheny Plateau physiographic provinces. The aquifer system, composed primarily of glacially derived sand and gravel deposits, includes both unconfined and confined aquifer units. For the approximately 2,300 people who live over the aquifers, groundwater is the sole source of drinking water. In some areas, such as the Village of Freeville, artesian conditions exist, with potentiometric heads exceeding 25 feet above land surface. This report describes the hydrogeologic framework, aquifer distribution, groundwater-flow dynamics, and water quality of the lower Fall Creek aquifer system. Groundwater is generally hard and contains elevated concentrations of redox-sensitive constituents, such as iron and manganese relative to U.S. Environmental Protection Agency secondary maximum contaminant levels. Arsenic was not detected above the U.S. Environmental Protection Agency maximum contaminant level during this study but remains a regional concern due to past detections in the study area and nearby aquifers. The findings provide data for local planners and resource managers to use in understanding and protecting this water resource.

New York↗

Hydrogeology and water quality of the unconsolidated aquifer system in the Lower Fall Creek Valley, Tompkins County, New York

The U.S. Geological Survey, in cooperation with Tompkins County and the Town of Dryden, conducted a hydrogeologic investigation of the lower Fall Creek unconsolidated aquifer system in Tompkins County, New York. This study is part of an ongoing effort to characterize major valley-fill aquifers in the region to support sustainable groundwater management. The lower Fall Creek valley, encompassing the Village of Freeville, the unincorporated communities of Etna, and the hamlets of McLean and Varna, has experienced moderate development and relies heavily on groundwater for residential, agricultural, and commercial use. The study area spans approximately 51 square miles along the boundary of the Allegheny Plateau and North Rim of the Allegheny Plateau physiographic provinces. The aquifer system, composed primarily of glacially derived sand and gravel deposits, includes both unconfined and confined aquifer units. For the approximately 2,300 people who live over the aquifers, groundwater is the sole source of drinking water. In some areas, such as the Village of Freeville, artesian conditions exist, with potentiometric heads exceeding 25 feet above land surface. This report describes the hydrogeologic framework, aquifer distribution, groundwater-flow dynamics, and water quality of the lower Fall Creek aquifer system. Groundwater is generally hard and contains elevated concentrations of redox-sensitive constituents, such as iron and manganese relative to U.S. Environmental Protection Agency secondary maximum contaminant levels. Arsenic was not detected above the U.S. Environmental Protection Agency maximum contaminant level during this study but remains a regional concern due to past detections in the study area and nearby aquifers. The findings provide data for local planners and resource managers to use in understanding and protecting this water resource.

EarthArXiv↗

Deep learning models for mapping surficial geology in selected physiographic regions of New York

Surficial geologic mapping is required for decisions including infrastructure and water-resource development and management. Deep learning, a type of machine learning that uses training data to self-learn and perform tasks, is explored as a tool to help increase efficiency of the labour- and time-intensive mapping process. Deep learning models were trained, and their potential to aid in mapping surficial geology was explored for two physiographic regions in New York: the high-relief Allegheny Plateau and the low-relief Erie-Ontario Lowlands. Key to the development of deep learning models is the availability of highly detailed surficial geologic maps in each of these regions from which the models can learn. Through exploring different groupings of surficial deposits and their corresponding spatial data, deep learning models were iteratively developed to reproduce published mapping for greater than 79% of the training areas, and with similar accuracy in test areas within the same physiographic region. These models were trained using only two inputs: high-resolution lidar data and previously published surficial geologic maps. This straightforward approach is intended to make the methods and models created easily reproducible using widely accessible datasets. This research describes the strengths and limitations of lidar-derived surficial geologic models and how broadly grouping surficial types, characteristics of different physiographic regions, and testing models in physiographic regions outside of their training areas affect model creation and performance. In addition, the models created could be used as a tool for mapping surficial geology in similar physiographic regions, in addition to establishing a framework for creating similar models elsewhere.

New York↗

Shallow hydrogeologic framework of the Tully Valley mudboil area, Onondaga County, New York

Mudboils have been documented in the Tully Valley in southern Onondaga County, New York, since the late 1890s. Sediment-laden water from the mudboils flows into Onondaga Creek, which empties into Onondaga Lake at Syracuse 15 miles to the north. Turbidity from the mudboils has degraded the water quality of Onondaga Creek despite a series of mitigation efforts that began in the early 1990s. Turbidity mitigation actions presently (2025) being considered include creek relocation and offline sediment settling. In support of these proposed actions during 2021–23, the U.S. Geological Survey, in cooperation with the New York State Department of Environmental Conservation, U.S. Environmental Protection Agency, Onondaga Nation, Onondaga Environmental Institute, and Central New York Regional Planning and Development Board, collected and analyzed geologic, hydrologic, geophysical, and geotechnical data to characterize the shallow hydrogeology along four proposed creek-relocation paths and in the proposed offline settling basin area. The investigation indicated that the four proposed creek-relocation paths, two east of Onondaga Creek and two west of Onondaga Creek, are underlain by sediments including muck, alluvium, mudboil deposits, alluvial-fan sand and gravel, and lacustrine fines. The proposed excavations would penetrate partially to fully saturated conditions: generally, the water table is shallow near the creek and deep on the alluvial fans. The shallowest excavation, about 5 feet below land surface, would be near the creek and primarily in alluvium, and the deepest excavation, as much as 30 feet below land surface, would be in the alluvial-fan deposits. Brackish waters would be penetrated by proposed channel excavations on the eastern side of Onondaga Creek in an area downgradient from a potentially leaking historical salt-exploration borehole and near the main mudboil area. Excavation in these areas likely would provide a continuous source of brackish groundwater to the relocated creek. Proposed channel excavations of muck, soft to very soft lacustrine fines, and mudboil-type sediments in mudboil and suspected mudboil areas would pose an excavation and slope stability challenge and would have the greatest potential to create new mudboils. Proposed channel excavations below the water table on the Rattlesnake Gulf and Rainbow Creek alluvial fans would intercept groundwater and make the constructed streambank susceptible to seepage-induced slope instability. The substantial water-level fluctuation in the sediments of both alluvial fans would aggravate the stability condition. In addition, excavation on the Rattlesnake Gulf alluvial fan would have the potential to affect water-supply springs at the toe of the fan. The proposed offline settling basin area is in the northern part of the Rattlesnake Gulf alluvial fan. Natural and man-made diversions of Rattlesnake Gulf have resulted in saturated conditions in the general area of the proposed basin. The proposed offline settling basin would be excavated in, and berms would be constructed on, alluvial-fan deposits and lacustrine fines. In the proposed basin area, the alluvial deposits overlying the lacustrine fines are less than 10 feet thick. Excavation, berm construction, and loading of the saturated, soft to very soft lacustrine fines may be problematic and require soil strengthening.

New York↗

A targeted approach for mapping groundwater discharge to surface water and fish thermal refuge in four Lake Ontario tributaries

The duration, magnitude, and frequency of heatwaves are predicted to increase in the coming decades, a combination that can reduce the survival of many fish species. Across the world, there is broad interest in identifying thermal refuge for heat-intolerant fish species and exploring opportunities to enhance or protect these areas. Because deeper groundwater maintains a relatively constant temperature, groundwater-influenced areas along streams can provide cool-water refuge for fish during periods of extreme heat. A targeted approach was developed for identifying existing cold-water zones and areas of substantial groundwater discharge in four high priority Lake Ontario tributaries. Our approach included: (1) predicting where groundwater discharge is most likely with a simple geospatial model and (2) using model predictions to select field sites for intensive high-resolution study, including ground-based mapping of groundwater features (springs, seeps, tributaries) as well as drone-based optical and thermal infrared surveys. Results from field sites were used to both verify model performance and map different types and aerial extents of thermal anomalies. Geospatial modelling successfully predicted regions of widespread groundwater upwelling, later verified and mapped by field and drone surveys. Comparison of model and field survey results further highlighted specific geospatial layers, such as soil/bedrock types and topographic wetness index, as being particularly useful for predicting groundwater influence on streams in the study area. In addition, a comparison of geospatial model results with a model of fish abundances along the studied streams showed significant positive correlations for many heat-intolerant fish species over a wide geographic area. The approach developed in this study can be applied to other watersheds to highlight areas of probable groundwater discharge and could be used by fishery and water resource managers to support cold-water fish habitat management decision-making and resource conservation.

New York↗

Detecting hidden sedimentary geothermal systems in the Upper Colorado River Basin

Geothermal resources exist in sedimentary rock where circulation of water for efficient extraction or storage of heat is possible. Except in rare instances where hot water is expressed at the land surface, sedimentary geothermal resources are hidden, so the identification of these systems is optimally accomplished using predictive subsurface modeling. An integrated approach using detailed paleogeographic interpretations, subsurface geologic mapping, and numerical modeling has produced regional geologic and temperature models for the Upper Colorado River Basin, a large watershed in central North America that contains many sedimentary basins. These models identify areas of hidden sedimentary geothermal resource potential in low temperature (<90°C), moderate temperature (90–150°C), and high temperature (>150°C) fairways across the study area. These models incorporate maps of key horizons in outcrop and the subsurface to create a robust structural framework that can be used to target favorable geology for natural or engineered permeability. This framework is populated with lithologies derived from detailed palaeogeographical maps and over 40,000 bottom hole temperature (BHT) values were used to create a calibrated three-dimensional (3D) temperature model across the region. The resulting maps serve as a regional sedimentary geothermal play fairway screening tool for evaluating different grades of sedimentary geothermal resources and for identifying areas of interest where more detailed, prospect-scale studies can be undertaken.

Arizona, Colorado, New Mexico, Utah, Wyoming↗

An evaluation of the effects of different deicing salt application rates on three watersheds in Essex County, New York

The U.S. Geological Survey, in cooperation with the New York State Department of Transportation, evaluated the effects of different deicing salt application rates on surface water, groundwater, and highway runoff quality near State highways in northern New York. Three reaches of State highways were tested with different deicing treatments between October 2019 and November 2022: a salt-sand mixture (Treatment A), a salt mixture applied at a lower rate (Treatment B), and a control mixture consistent with typical deicing salt amounts and application rates. Data on pavement conditions and the quality of surface water, highway runoff, and groundwater were collected. Surface electromagnetic data were also collected. Surface-water and groundwater quality downgradient from the State highways were compared with water quality at upgradient locations. The percentage of snow or ice coverage was used to evaluate the effectiveness of the salt applications. This report provides an overview of the transport of deicing salt. The Treatment B watershed had deicing mixture applied more frequently than other highway reaches, which caused it to have the highest annual total chloride application. Despite differences in chloride application, flow-weighted mean chloride concentrations in highway runoff were comparable across treatments. Chloride concentrations were elevated in surface water and groundwater downgradient from highways relative to chloride concentrations upgradient from highways. A chloride mass balance, calculated for one treatment watershed, indicated that groundwater affected by legacy deicing practices may be contributing additional chloride to surface water. Spatial patterns from electromagnetic surveys show a shallow saline plume alongside the highway in that area. Differences in winter severity and pavement-surface conditions drove deicing salt applications in the treatment areas. This study found that several factors affect chloride loads in the watersheds, including variable winter conditions, adaptive snow and ice management, legacy management practices, and area-specific aquifer and groundwater conditions.

New York↗

Hydrogeology of unconsolidated and bedrock aquifers along the Salmon River, including Malone, Franklin County, New York

The U.S. Geological Survey, in cooperation with the New York State Department of Environmental Conservation, investigated the hydrogeology of the unconsolidated and bedrock aquifers along the Salmon River corridor in northern Franklin County, New York. The study area covers roughly 147 square miles and includes the village of Malone and parts of the Towns of Malone, Bellmont, Burke, Constable, Westville, Bangor, Duane, and Franklin, New York. Groundwater is the primary source for water supply within the study area. Eighty-three percent of all residents use public water supplied from two production wells that draw water from a thick, highly productive sand-and-gravel aquifer likely receiving induced infiltration from the Salmon River. Twenty-four percent of the 187 verified domestic wells in the study area outside of the production well service area boundary are screened in typically discontinuous deposits of stratified sand and sand and gravel. Seventy-six percent of the wells are completed in bedrock aquifers including the Potsdam Sandstone, metamorphic rocks, Theresa Formation, and unknown bedrock. Characterizing and understanding potential groundwater resources is critical for protecting the quality of the groundwater. The information in this report may be used to guide delineation of groundwater contributing areas, assess potential threats to aquifers from both point and nonpoint sources, respond to contamination from spills or leaks from underground storage facilities or other sources, and to support assessments for future development of municipal water supplies.

New York↗

Random forest regression models for estimating low-streamflow statistics at ungaged locations in New York, excluding Long Island

This report is currently being updated with re‑calibrated random‑forest models for StreamStats implementation. The basin‑characteristic inputs were modified to align with StreamStats requirements. Abstract Models to estimate low-streamflow statistics at ungaged locations in New York, excluding Long Island and including hydrologically connected basins from bordering States, were developed for the first time by the U.S. Geological Survey, in cooperation with the New York State Department of Environmental Conservation. A total of 224 basin characteristics were developed for 213 unaltered streamgages (locations where the human effects on streamflow were limited), across the following categories: basin geometry, climate, land cover, soils, surficial geology, and other characteristics. The basins with unaltered streamgages were evaluated for potential redundancy, and streamgages in close proximity and with similar drainage areas were flagged and removed from the testing and cross-validation datasets to prevent data leaking from the training dataset to the testing dataset. Random forest regression models were created by using basin characteristics as predictor variables and by developing a workflow to train, tune, and test the model. Models were developed to estimate the ungaged lowest annual 7-day and 30-day average streamflow that occurs (on average) once every 10 years (7Q10 and 30Q10). The top four basin characteristics used for the 7Q10 and 30Q10 models were drainage area, total stream length, perimeter of the basin, and length of the longest flow path. Results for the 7Q10 and 30Q10 models had coefficients of determination ( R 2 ) of 0.796 and 0.853, respectively. The output model results were bias-corrected for ungaged locations across New York and are available within the interactive StreamStats tool.

New York↗

Hydrogeologic framework and considerations for drilling and grouting of closed-loop geothermal bores in the Erie-Ontario Lowlands and Allegheny Plateau of New York State

The hydrogeologic framework at closed-loop geothermal sites in the Erie-Ontario Lowlands and Allegheny Plateau of central and western New York is the result of the complex interaction of bedrock geology, glacial geology, and groundwater hydrology, and the occurrence of petroleum and gas. Considerations for closed-loop geothermal bore installation include the thickness and character of glacial deposits, bedrock solubility and depth to competent rock, karst development, the distribution of highly permeable zones and their hydraulic heads, and the presence of saline water, gas, and oil. The hydrogeology of the Erie-Ontario Lowlands and Allegheny Plateau poses challenges to closed-loop geothermal bore drilling and casing; managing drill cuttings, discharge water, and gas; and grouting. The potential to encounter severe challenges typically increases with bore depth. This report highlights hydrogeologic considerations for closed-loop geothermal bore installation in New York’s Erie-Ontario Lowlands and Allegheny Plateau to help guide the efficient and safe development of geothermal resources in the regions.

New York↗

A geospatial analysis of water-quality threats from orphan wells in principal and secondary aquifers of the United States

Throughout the history of oil and gas production in the United States, millions of wells have been drilled for exploration and energy production. Hundreds of thousands of unplugged wells are no longer actively producing and are currently under orphan status, with no responsible party obligated for plugging. Orphan wells can pose threats to water resources by providing pathways for contaminants such as hydrocarbons and brines to migrate into water-supply aquifers. In this study, we investigate the potential threats to groundwater resources posed by orphan wells at the national scale. Water-quality data is extremely sparse in relation to orphan wells nationally and may not be suitable for identifying contamination from oil and gas development. We used geospatial and statistical methods to evaluate which principal and secondary aquifer systems may be most susceptible to contamination from orphan wells. Analysis involved three sets of susceptibility factors including: 1) factors related to the number and density of orphan wells; 2) factors that can threaten well integrity and contribute to transport of contaminants; and 3) factors related to groundwater withdrawal rates and the affected populations/communities in the event of water quality disturbances. From a dataset of 117,672 documented orphan wells, 64,203 fall within a principal aquifer system, while the remainder fall within a secondary aquifer system. By assessing the combination of well integrity and hydrogeologic factors within these aquifer systems, five groupings of principal aquifers were identified, where groups ranged from aquifer systems with high numbers of orphan wells, multiple well integrity threats and high withdrawals, to aquifers with a relatively low number of orphan wells, limited well integrity threats and minimal water use. Three regions of the country emerge containing aquifers with higher susceptibility to contamination from orphan oil and gas wells. These regions include 1) The Appalachian Basin (including the Pennsylvanian Aquifer System), 2) The Gulf Coast Aquifers (including the Coastal Lowlands Aquifer system) and 3) The California Aquifers (including the California Coastal Basin Aquifer system). This work is the first multivariate geospatial investigation of orphan wells and groundwater resources on a national scale, and sheds light on which aquifers are most susceptible to groundwater contamination from orphan wells.

Science of the Total Environment↗

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↗

Low-flow statistics for selected streams in New York, excluding Long Island

The U.S. Geological Survey, in cooperation with the New York State Department of Environmental Conservation, updated low-streamflow statistics for New York, excluding Long Island and including hydrologically connected watersheds in bordering States, for the first time since 1972. Historical daily streamflow data for active and inactive gages were considered for this study with periods of record as recent as March 31, 2022, adding 50 years of data to the last published low-streamflow statistics report for New York and including 119 new gages where low-streamflow statistics are calculated for the first time. Gages were evaluated across several criteria to identify gages that were not suitable for low-streamflow frequency analysis. In addition, gages were evaluated for the presence of alteration within the streamflow period of record based on previous studies and U.S. Geological Survey National Water Information System site metadata including peak flow codes. A trend analysis was performed using the Wilcoxon rank-sum hypothesis test comparing data from the most recent 30 years of record to data from 30 years and earlier for each long-record gage (30 years or more of available data). Results from the trend analysis indicated that 45 unaltered and 32 altered long-record sites had a statistically significant trend for the annual minimum n -day time series; most gages showed increasing trends in the annual minimum n -day time series. Low-streamflow statistics were calculated using the most recent 30 years of record for gages with a statistically significant trend. Before and after 1972, the lowest annual 7-day and 30-day average streamflow that occurs (on average) once every 10 years (7Q10 and 30Q10 statistics respectively) increased significantly at 41 unaltered gages and decreased significantly at 3 unaltered gages where data were available.

New York↗

Areas contributing recharge to selected production wells in unconfined and confined glacial valley-fill aquifers in Chenango River Basin, New York

In the Chenango River Basin of central New York, unconfined and confined glacial valley-fill aquifers are an important source of drinking-water supplies. The risk of contaminating water withdrawn by wells that tap these aquifers might be reduced if the areas contributing recharge to the wells are delineated and these areas protected from land uses that might affect the water quality. The U.S. Geological Survey, in cooperation with the New York State Department of Environmental Conservation and the New York State Department of Health, began an investigation in 2019 to improve understanding of groundwater flow and delineate areas contributing recharge to 16 production wells clustered in three study areas in the basin as part of an effort to protect the source of water to these wells. Areas contributing recharge were delineated on the basis of numerical steady-state groundwater-flow models representing long-term average hydrologic conditions. In the Cortland study area, four water suppliers operate 10 production wells that withdraw a total average rate of 2,480 gallons per minute from an unconfined aquifer consisting of well-sorted sand and gravel deposits. Simulated areas contributing recharge to these wells at their average pumping rates covered a total area of 6.93 square miles. Simulated areas contributing recharge extend upgradient from the wells to upland till deposits and to groundwater divides. Some simulated areas contributing recharge include isolated areas remote from the wells. Short simulated groundwater traveltimes from recharging locations to discharging wells indicated that the wells are vulnerable to contamination from land-surface activities; 50 percent of the traveltimes were 10 years or less. Land cover in some of the areas contributing recharge included a substantial amount of urban and agriculture land use. The groundwater-flow model of the Cortland study area was calibrated to available hydrologic data by inverse modeling using nonlinear regression. The parameter variance-covariance matrix from model calibration was used to create parameter sets that reflect the uncertainty of the parameter estimates and the correlation among parameters to evaluate the uncertainty associated with the single, predicted contributing areas to the wells. This analysis led to contributing areas expressed as a probability distribution. Because of the effects of parameter uncertainty, the size of the probabilistic contributing areas was larger than the size of the single, predicted contributing area for the wells. Thus, some areas not in the single, predicted contributing area might actually be in the contributing area, including additional areas of urban and agriculture land use that have the potential to contaminate groundwater. Additional areas that might be in the contributing area included recharge originating near the pumping wells that have relatively short groundwater-flow paths and traveltimes. In each of the Greene and Cincinnatus study areas, one water supplier operates three wells that are screened near the top of the bedrock surface in a confined aquifer consisting of poorly to well-sorted sand and gravel deposits. This confined aquifer is overlain by a lacustrine confining unit of very fine sand, silt, and clay, which in turn is overlain by a thin unconfined aquifer of sand and gravel. The groundwater-flow models for these two areas were manually calibrated because of the limited hydrologic data. Simulated areas contributing recharge to the Greene study area wells covered a total area of 0.35 square mile for the average pumping rate of 170 gallons per minute. The contributing areas extended southeastward of the wells to the groundwater divide in the till uplands. The contributing areas also included remote, isolated areas on the opposite side of the Chenango River from the wells primarily in the till uplands. For the Cincinnatus study area wells, which have a low average pumping rate (34 gallons per minute), the simulated contributing areas totaled 0.06 square mile and were on the same side of the river as the wells, but they are isolated areas remote from the wells primarily in the till-covered bedrock uplands. Land cover in these contributing areas for both study areas is primarily agriculture and forested, with the contributing areas to the Greene study area wells also including some urban land uses. Because the Greene and Cincinnatus study area wells are screened relatively deep and some flow paths to the wells partly travel through the confining unit, which impedes the connection with surface sources of recharge, overall groundwater traveltimes are greater than for wells in the Cortland study area. Fifty percent of Cortland study area wells, but only 9 and 44 percent of Greene and Cincinnatus study area wells, respectively, have groundwater traveltimes of 10 years or less.

New York↗

Areas contributing recharge to priority wells in valley-fill aquifers in the Neversink River and Rondout Creek drainage basins, New York

In southeastern New York, the villages of Ellenville, Wurtsboro, Woodridge, the hamlet of Mountain Dale, and surrounding communities in the Neversink River and Rondout Creek drainage basins rely on wells that pump groundwater from valley-fill glacial aquifers for public water supply. Glacial aquifers are vulnerable to contamination because they are highly permeable and have a shallow depth to water table. To protect the quality of these water resources, water managers need accurate information about the areas that contribute recharge to production wells that pump from these aquifers. The New York State Department of Environmental Conservation and the New York State Department of Health designated eight priority wells in this region for which water supply protection is of primary concern. The U.S. Geological Survey, in cooperation with the New York State Department of Environmental Conservation and the New York State Department of Health, began an investigation in 2019 with the general objectives of (1) improving understanding of regional groundwater-flow system, (2) delineating areas contributing recharge to eight priority production wells, and (3) quantifying the uncertainty of these contributing areas in a probabilistic way that can be used to inform decision-making related to priority well source-water protection. To complete these objectives, a MODFLOW 6 groundwater model was created encompassing the eight priority wells and the surrounding flow system, which includes parts of the Neversink River and Rondout Creek Basins in Sullivan County and Ulster County, New York. The model was built using Python tools (such as flopy , modflow-setup , and sfrmaker ) that facilitate transparent and repeatable model development using existing datasets. The model parameters were estimated with a stepwise approach using an iterative ensemble smoother implementation of the Parameter ESTimation software PEST++ (version 5.0.0). We evaluated initial “best guess” parameter bounds with a prior Monte Carlo analysis. Results of the first prior Monte Carlo analysis were used to make informed adjustments to model parameter bounds (typically resulting in expanded bounds), and a second prior Monte Carlo analysis was run to identify improved ranges for model parameters during history matching. The history matching effort produced an ensemble of parameter values for the groundwater-flow model that spans the range of values within prior uncertainty bounds. The ensemble is informed by the historical observation data, within a reasonable range of uncertainty on those observations. This history-matched ensemble was used in a particle tracking Monte Carlo analysis to delineate the areas contributing recharge to priority wells. The groundwater-flow and particle tracking (MODPATH7) models were run once for each ensemble member. Deterministic contributing areas computed for each ensemble member were aggregated to produce maps showing the probability that a location contributes recharge to priority wells. Finally, the particle tracking Monte Carlo analysis was repeated for six pumping scenarios, representing a wide range of possible pumping levels, to incorporate uncertainty in future pumping rates related to population growth or other management decisions. Increasing pumping rates generally led to larger contributing recharge areas and larger areas of high probability that a location contributes recharge to priority wells. These maps show the overall uncertainty of the areas contributing recharge to priority wells in the study area and provide a tool for risk-based decision making for protection of well source water.

New York↗

Data sources and methods for digital mapping of eight valley-fill aquifer systems in upstate New York

Digital hydrogeologic maps were developed in eight study areas in upstate New York by the U.S. Geological Survey in cooperation with the New York State Department of Environmental Conservation. The digital maps define the hydrogeologic framework of the valley-fill aquifers and surrounding till-covered uplands in the vicinity of the villages of Ellenville and Wurtsboro and hamlets of Woodbourne and South Fallsburg in Sullivan and Ulster Counties, town of Greene in Chenango County, city of Cortland and town of Cincinnatus in Cortland County, city of Jamestown in Chautauqua County, city of Olean and village of Ellicottville in Cattaraugus County, and villages of Fishkill and Wappinger Falls in Dutchess County. The hydrogeologic framework provided the foundation for groundwater-flow models that were used in the delineation of areas contributing groundwater flow to production wells screened in four of the eight valley-fill aquifers considered in this study. The hydrogeologic framework for the other four study areas was developed for potential future use in groundwater contributing-area studies. Data used in the creation of all digital surfaces and thicknesses included published surficial geology; aquifer maps and hydrogeologic sections; light detection and ranging (lidar) datasets; the Soil Survey Geographic Database; and lithologic well logs from the National Water Information System, New York State Department of Environmental Conservation, New York State Department of Transportation, and Empire State Organized Geologic Information System databases. Digital maps of the surficial geology; thickness of the surficial sand and gravel aquifers; and tops of the confining lacustrine silt and clay units, confined sand and gravel aquifers, and bedrock surfaces were created by using ArcGIS (a geographic information system). All surfaces and thicknesses were generated by using one of the following ArcGIS interpolation tools: Topo to Raster, Natural Neighbors, Kriging, or Empirical Bayesian Kriging. The datasets developed in this study provide a greater understanding of the underlying hydrogeologic framework in glacial valley-fill aquifers and can be applied in the evaluation of groundwater-supply development and protection.

New York↗