Search

In 2016, the U.S. Geological Survey began a multiyear cooperative study with the New York State Department of Environmental Conservation to evaluate the sustainability of Long Island’s sole-source aquifer system through hydrogeologic mapping, compilation of groundwater chloride concentrations, and groundwater flow modeling. In the initial phase of the islandwide study, the hydrogeologic framework and extent of saltwater intrusion in aquifers in Kings, Queens, and Nassau Counties on western Long Island, N.Y., were investigated. The aquifer system underlying western Long Island has been under stress from pumping of public, irrigation (golf course), and industrial supply wells. Saltwater intrusion has occurred from surrounding embayments (East River, Long Island Sound, Jamaica Bay, and the Atlantic Ocean) due to pumping. Eighteen boreholes were drilled and cores taken during 2019–21 to collect hydrogeologic, geochemical, and geophysical data to delineate the complex subsurface hydrogeology and extent of saltwater intrusion within the study area. Evaluation of the new cores, reexamination of legacy core descriptions, and analysis of borehole geophysical logs was used to refine the previously published hydrogeologic framework of Pleistocene and Cretaceous unconsolidated sediments in the area, including delineation of a previously undefined hydrogeologic unit between the Magothy aquifer and the Raritan confining unit, herein named the “upper Raritan aquifer.” The upper Raritan aquifer was first recognized in southeastern Nassau County from an analysis of about 50 closely spaced boreholes with high-resolution core descriptions and gammaray (gamma) logs. Further analysis of borehole logs across the study area indicated that the upper Raritan aquifer was also present in Kings and Queens Counties. Nuclear magnetic resonance (NMR) logging was used for the first time on Long Island to provide estimates of the hydraulic properties of the major aquifer and confining units. Unlike other geophysical logs that record responses to the rock matrix and fluid properties and are strongly dependent on mineralogy, NMR logs record responses to the presence of hydrogen protons in the formation fluid to determine water fraction and pore-size distribution. NMR log analysis provided estimates of the clay-bound, capillary-bound, and mobile water fractions and hydraulic conductivity of aquifers and confining units penetrated by five wells in Nassau County. Pumpage for public-supply and industrial wells on Long Island began in the 1870s with small, localized suppliers of populated areas in Kings and Queens Counties. By 1904–16, pumpage for public water supply in Kings County averaged 21 million gallons per day, and averaged 37 million gallons per day in Queens County, mostly from the upper glacial aquifer. Saltwater intrusion was reported as early as the beginning of the 20th century and included the upper glacial- Jameco-Magothy and Lloyd-North Shore aquifer systems. By 1936, pumping in central Kings County created a major cone of depression in the water table extending to the south shore of much of Kings County and into southwestern Queens County. Saltwater intrusion has caused the shutdown of public-supply wells in Kings, Queens, and Nassau Counties. A large saltwater intrusion wedge in the Lloyd aquifer was indicated in southern Queens County in the early part of the 20th century, and the saltwater interface may have been onshore predevelopment. Most of Kings and Queens Counties are intruded with saltwater in both the upper glacial-Jameco-Magothy and Lloyd-North Shore aquifers systems. Saltwater increased during the 20th century and continues to increase to the present (2023) in the Lloyd-North Shore aquifer system in Great Neck and Manhasset Neck in northern Nassau County. A major wedge of saltwater intrusion in the upper glacial-Jameco- Magothy aquifer in southwestern Nassau County appears to be increasing. [ABSTRACT FROM AUTHOR]

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 decisionmaking 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 1U.S. Geological Survey 2INTERA Geoscience and Engineering Solutions 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. [ABSTRACT FROM AUTHOR]