Geology topics

Paul J. Friesz

Publications and source records attributed to Paul J. Friesz.

13 recordsLinked to original sources

USGS2023

Evaluation of alternative groundwater-withdrawal scenarios on water levels in Kingsbury Pond, upper Charles River Basin, eastern Massachusetts

Kingsbury Pond is a glacial kettle pond in the town of Norfolk, Massachusetts, in the Mill River Basin, which is part of the Upper Charles River Basin in eastern Massachusetts. The pond is hydraulically connected to the surrounding groundwater-flow system, and water levels in the pond fluctuate in response to recharge to the aquifer from precipitation and wastewater return flows through septic systems, to withdrawals from the aquifer at nearby wells, and to precipitation directly on the pond surface. Concerns about the effects of withdrawals on water levels in the pond prompted an investigation to better understand the hydrology of Kingsbury Pond and its response to groundwater withdrawals and to determine if withdrawals from wells in Franklin, Mass., can be modified to simultaneously reduce the effect on water levels in the pond and yet meet the water-supply demands of the Town of Franklin. An existing, transient groundwater-flow model of the Upper Charles River Basin was modified for this study in the area near Kingsbury Pond to improve representation of the hydrologic system near the pond. The mean annual water-level altitude simulated for the pond for nonpumping conditions using the modified model is 136 feet (ft), which falls within the range of likely annual pond-altitude fluctuations of 135 to 140 ft estimated for average hydrologic conditions before the beginning of withdrawals at two nearby wells operated by the Town of Franklin (wells FR–04 and FR–05). The mean annual water-level altitude at the pond decreased by 3.8 ft to 132.2 ft for simulated mean monthly withdrawal rates at all wells within the Upper Charles River Basin from 2010 to 2019 (referred to as the baseline withdrawal condition). A groundwater management model that links the groundwater-flow model with a mathematical optimization method was developed to evaluate the effects of three alternative groundwater-withdrawal scenarios for the Franklin public-water system on water levels in Kingsbury Pond. In the first scenario, monthly withdrawal rates at wells FR–04 and FR–05 were increased from the baseline withdrawal rates to their maximum authorized rates for all months of the year; all other Franklin wells were specified at their baseline withdrawal rates. This scenario resulted in a mean annual water-level altitude at the pond of 129.3 ft, or a mean annual decline of 6.7 ft compared with nonpumping conditions and a decline of 2.9 ft compared with baseline conditions. The results of the second scenario showed that water levels in the pond can be increased relative to 2010–19 conditions while meeting Franklin’s 2010–19 monthly water-supply demands if withdrawals at wells FR–04 and FR–05 were shifted to other Franklin wells. In this scenario, monthly withdrawal rates at wells FR–04 and FR–05 were decreased from their baseline rates to one-third their maximum practical rates for all months of the year; increased withdrawal rates at other Franklin wells were determined by the management model. The decrease in withdrawal rates at wells FR–04 and FR–05 resulted in a mean water-level altitude at the pond of 134.1 ft, which was equivalent to a 51 percent increase (improvement) in the mean annual water level of the pond relative to the baseline condition. A third scenario was done to determine if Franklin’s existing water-supply system has the capacity to meet the mean annual maximum permitted withdrawal rate of the system of 3.45 million gallons per day while maintaining monthly withdrawal rates at wells FR–04 and FR–05 at their 2010–19 rates and water levels in Kingsbury Pond at baseline conditions. The analysis indicated that the capacity of the system cannot meet the increased demand during some months of the year with withdrawal rates at the two wells fixed at their monthly 2010–19 rates; however, the existing system is capable of meeting about 90 percent of the maximum permitted rate (3.10 million gallons per day) by increasing withdrawal rates at other Franklin wells above their 2010–19 rates.

Massachusetts↗

USGS2022

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↗

USGS2021

Approaches for assessing long-term annual yields of highway and urban runoff in selected areas of California with the Stochastic Empirical Loading and Dilution Model (SELDM)

The California Department of Transportation, commonly known as CalTrans, and other municipal separate storm sewer system permittees in California as well as other State departments of transportation nationwide need information about potential loads and yields (loads per unit area) of constituents of concern in stormwater runoff and discharges from stormwater best management practices (BMPs). Although its National Pollution Discharge Elimination System stormwater permit is focused on areas subject to total maximum daily load (TMDL) regulations, CalTrans builds and maintains BMPs to minimize the adverse effects of roadway runoff on receiving waters throughout the State. This report describes approaches used by the U.S. Geological Survey in cooperation with CalTrans for using the Stochastic Empirical Loading and Dilution Model (SELDM) to assess long-term annual yields of highway and urban runoff in selected areas of California. In this study, a series of regional and local yields were simulated to provide statewide planning-level estimates and more refined TMDL-specific yield values. SELDM was used to analyze 368 State roadway and urban runoff yields for 53 runoff quality constituents. The analyses included 222 random-seed analyses, 60 regional State roadway-runoff analyses, 24 regional urban roadway-runoff analyses, and 62 focused TMDL-area analyses. This report describes approaches and statistics used to analyze available hydrologic and runoff quality data in all analyses. Results for all analyses are provided in the model archive, but only a selected subset of results are presented as examples in this report. State roadway runoff, urban runoff, and BMP discharge yields for total suspended solids, total nitrogen, total phosphorus, and total zinc were selected as examples because they are widespread constituents of concern with substantial amounts of State roadway and urban runoff monitoring data. In this report, a hypothetical basin was specified by using available geographic information to demonstrate use of the State roadway and urban runoff yields to estimate long-term annual stormwater loads from developed areas. Application of these yields to the hypothetical basin indicates that although State-roadway yields may be higher than urban-runoff yields for some constituents, State-roadway loads may be a small proportion of total stormwater loads because State roadways themselves are a small fraction of the total impervious area in such basins. Although application of results from this study may have considerable uncertainty for any particular stormwater outfall, the study does provide robust estimates to support basin-scale runoff-load analyses in California. These analyses also provide estimates for the 12 U.S. Environmental Protection Agency level III ecoregions that are completely or partially within the boundaries of the State of California.

California↗

USGS2019

Hydrogeologic framework and delineation of transient areas contributing recharge and zones of contribution to selected wells in the upper Santa Fe Group aquifer, southeastern Albuquerque, New Mexico, 1900–2050

The Santa Fe Group aquifer is an important source of water to communities within the Middle Rio Grande Basin, including the Albuquerque-Rio Rancho metropolitan area and Kirtland Air Force Base, New Mexico. In November 1999, Kirtland Air Force Base personnel observed fuel-stained soils at the Bulk Fuels Facility on the base. Subsequent pressure tests identified pipeline leaks. Fuels stored at the Bulk Fuels Facility have included aviation gasoline, jet propellant 4, and jet propellant 8. The fuels migrated about 480 feet down to the water table. Ethylene dibromide, the constituent making up the most extensive part of the plume and a component of leaded aviation gasoline, has formed a plume that, in December 2016, was 400 to 1,300 feet wide, extended about 5,800 feet northeast from the Bulk Fuels Facility, and was about 3,700 feet from the nearest downgradient water-supply well. Prior to widespread development of groundwater resources in southeastern Albuquerque, groundwater near the present-day location of the Bulk Fuels Facility flowed to the southwest. Groundwater began flowing northeast in about 1980 towards a large area of lowered water levels caused by groundwater pumping. In 2013 and 2014 the Albuquerque Bernalillo County Water Utility Authority, the U.S. Air Force, and the U.S. Geological Survey began a cooperative study to characterize the geology and hydrology of the Santa Fe Group aquifer in the vicinity of the ethylene dibromide plume and to develop a local-scale groundwater flow model to delineate areas contributing recharge and zones of contribution to selected water-supply wells. For this study, a previously developed Middle Rio Grande Basin regional groundwater-flow model was updated, and a smaller local-scale model was developed. Advective groundwater-flow paths were delineated and visualized with the MODPATH particle-tracking program. Of 11 wells included in the historical pumping analysis of areas contributing recharge, only wells K-3, K-7, and RC-4 derived a portion of their water from simulated recharge sources within the local-scale model. None of the areas contributing recharge overlap the Bulk Fuels Facility area or the ethylene dibromide plume footprint as delineated using December 2016 ethylene dibromide data. For the historical pumping analysis of zones of contribution, particles for the 11 selected wells generally moved southwest from the north and east boundaries of the local-scale model, moved past their target well, but reversed direction and moved back towards their target well after 1980 when groundwater flow changed to the northeast. Of the 11 wells, only BR-5, RC-5, and VH-2 had 1980–2013 particle pathlines that overlap the December 2016 ethylene dibromide plume footprint, and wells BR-5 and VH-2 have 1980–2013 particle pathlines that overlap the Bulk Fuels Facility area. Particles that were north of the Bulk Fuels Facility when groundwater flow reversed direction would not have the opportunity to interact with the ethylene dibromide plume. Wells BR-5, K-15, and VH-2 did have particles southwest of the Bulk Fuels Facility in 1980. Particles traveling to BR-5 and K-15 passed under or very near the Bulk Fuels Facility area in the 1980–2013 period, but none of the pathlines were shallow enough to interact with ethylene dibromide at the Bulk Fuels Facility. A few particles traveling to VH-2 passed through the Bulk Fuels Facility area at shallow enough depths to interact with ethylene dibromide at the Bulk Fuels Facility in the 1980–2013 period. Ethylene dibromide has not been detected in water samples collected in 2012 through 2015 from the VH-2 well. Of 10 water-supply wells near the ethylene dibromide plume included in the future pumping analysis of areas contributing recharge, only wells K-3, RC-3, and RC-4 had areas contributing recharge within the local-scale model. The areas contributing recharge for wells RC-3 and RC-4 do not overlap the Bulk Fuels Facility area or the December 2016 ethylene dibromide plume footprint, but K-3 derives part of its recharge prior to 1980 and during 1980–2015 from within the area of the December 2016 plume footprint. The analysis of the future pumping scenarios indicated that wells BR-5, K-3, K-16, RC-5, and VH-2 have pathlines for 1980–2015 and wells K-16 and VH-2 have pathlines for 2015–50 that when projected in plan view pass through the December 2016 plume footprint. Of these five wells, only K-3 and RC-5 have pathlines for 1980–2015 that are above an elevation of 4,800 feet and could interact with the ethylene dibromide plume if ethylene dibromide was present when the particles were present.

New Mexico↗

USGS2015

Areas contributing recharge to production wells and effects of climate change on the groundwater system in the Chipuxet River and Chickasheen Brook Basins, Rhode Island

The Chipuxet River and Chickasheen Brook Basins in southern Rhode Island are an important water resource for public and domestic supply, irrigation, recreation, and aquatic habitat. The U.S. Geological Survey, in cooperation with the Rhode Island Department of Health, began a study in 2012 as part of an effort to protect the source of water to six large-capacity production wells that supply drinking water and to increase understanding of how climate change might affect the water resources in the basins. Soil-water-balance and groundwater-flow models were developed to delineate the areas contributing recharge to the wells and to quantify the hydrologic response to climate change. Surficial deposits of glacial origin ranging from a few feet to more than 200 feet thick overlie bedrock in the 24.4-square mile study area. These deposits comprise a complex and productive aquifer system. Simulated areas contributing recharge to the production wells covered a total area of 0.63 square miles for average well withdrawal rates from 2007 through 2011 (total rate of 583 gallons per minute). Simulated areas contributing recharge for the maximum well pumping capacities (total rate of 3,700 gallons per minute) covered a total area of 2.55 square miles. Most simulated areas contributing recharge extend upgradient of the wells to morainal and upland till deposits and to groundwater divides. Some simulated areas contributing recharge include small, isolated areas remote from the wells. Relatively short groundwater traveltimes from recharging locations to discharging wells indicated that the wells are vulnerable to contamination from land-surface activities; median traveltimes ranged from 3.5 to 8.6 years for the production wells examined, and 57 to 91 percent of the traveltimes were 10 years or less. Land cover in the areas contributing recharge includes a substantial amount of urban and agriculture land use for five wells adjacent to the Chipuxet River; for one well adjacent to a tributary stream, land use is less developed. The calibrated groundwater-flow model provided a single, best representation of the areas contributing recharge to a production well. 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 predicted contributing areas to the wells. A Monte Carlo analysis led to contributing areas expressed as a probability distribution that differed from a single deterministic contributing area. Because of the effects of parameter uncertainty, the size of the probabilistic contributing areas for both average and maximum pumping rates was larger than the size of the deterministic contributing areas for the wells. Thus, some areas not in the deterministic contributing area might actually be in the contributing area, including additional areas of urban and agricultural land use that has 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. At the maximum pumping rates, areas associated with low probabilities extended long distances along groundwater divides in the uplands remote from the wells. Climate projections for the Chipuxet River and Chickasheen Brook Basins from downscaled output from general circulation models indicate that mean annual temperature might increase by 4.7 degrees Fahrenheit and 8.0 degrees Fahrenheit by the late 21st century (2070–99) compared with the late 20th century (1970–99) under scenarios of lower and higher emissions of greenhouse gases, respectively. By the late 21st century, winter and spring precipitation is projected to increase by 12 to 17 percent, summer precipitation to increase by about the same as mean annual precipitation (8 percent), and fall precipitation to decrease by 5 percent for both emission scenarios compared with the late 20th century. Soil-water-balance simulations indicate that, although precipitation is expected to increase in three seasons, only in winter do precipitation increases exceed actual evapotranspiration increases. Recharge is projected to decrease in fall and generally change little in spring and summer. By the late 21st century, winter recharge is expected to increase by 13 percent for the lower emissions scenario and by 15 percent for the higher emissions scenario. In fall, recharge is projected to diminish by 13 percent for the lower emissions scenario and by 24 percent for the higher emissions scenario. Although recharge is projected to change seasonally in the 21st century, mean annual recharge changes minimally. Soil moisture is projected to decrease in the 21st century from spring through fall because of increases in potential evapotranspiration, and in fall because of decreases in precipitation in addition to increases in potential evapotranspiration. By the late 21st century, soil moisture for the lower emissions scenario is expected to decrease by 11 percent in summer and 15 percent in fall, and for the higher emissions scenario, decrease by 23 percent for both seasons. These decreases in soil moisture during the growing season might have implications for agriculture in the study area. Predicted changes in the magnitude and seasonal distribution of recharge in the 21st century increase simulated base flows and groundwater levels in the winter months for both emission scenarios, but because of less recharge in the fall and less or about the same recharge in the preceding months of spring and summer, base flows and groundwater levels in the fall months decrease for both emission scenarios. October has the largest base flow and groundwater level decreases. By the late 21st century, base flows at the Chipuxet River in October are projected to decrease by 9 percent for the lower emissions scenario and 18 percent for the higher emissions scenario. For a headwater stream in the upland till with shorter groundwater-flow paths and lower storage properties in its drainage area, base flows in October are projected to diminish by 28 percent and 42 percent for the lower and higher emissions scenarios by the late 21st century. Groundwater level changes in the uplands show substantial decreases in fall, but because of the large storage capacity of stratified deposits, water levels change minimally in the valley. By the late 21st century, water levels in large areas of upland till deposits in October are projected to decrease by up to 2 feet for the lower emissions scenario, whereas large areas decrease by up to 5 feet, with small areas with decreases of as much as 10 feet, for the higher emissions scenario. For both emission scenarios, additional areas of till go dry in fall compared with the late 20th century. Thus projected changes in recharge in the 21st century might extend low flows and low water levels for the year later in fall and there might be more intermittent headwater streams compared with the late 20th century with corresponding implications to aquatic habitat. Finally, the size and location of the simulated areas contributing recharge to the production wells are minimally affected by climate change because mean annual recharge, which is used to determine the contributing areas to the production wells, is projected to change little in the 21st century.

Rhode Island↗

USGS2012

Evaluating prediction uncertainty of areas contributing recharge to well fields of multiple water suppliers in the Hunt-Annaquatucket-Pettaquamscutt River Basins, Rhode Island

Three river basins in central Rhode Island-the Hunt River, the Annaquatucket River, and the Pettaquamscutt River-contain 15 production wells clustered in 4 pumping centers from which drinking water is withdrawn. These high-capacity production wells, operated by three water suppliers, are screened in coarse-grained deposits of glacial origin. The risk of contaminating water withdrawn by these well centers may be reduced if the areas contributing recharge to the well centers are delineated and these areas protected from land uses that may affect the water quality. The U.S. Geological Survey, in cooperation with the Rhode Island Department of Health, began an investigation in 2009 to improve the understanding of groundwater flow and delineate areas contributing recharge to the well centers as part of an effort to protect the source of water to these well centers. A groundwater-flow model was calibrated by inverse modeling using nonlinear regression to obtain the optimal set of parameter values, which provide a single, best representation of the area contributing recharge to a well center. Summary statistics from the calibrated model were used to evaluate the uncertainty associated with the predicted areas contributing recharge to the well centers. This uncertainty analysis was done so that the contributing areas to the well centers would not be underestimated, thereby leaving the well centers inadequately protected. The analysis led to contributing areas expressed as a probability distribution (probabilistic contributing areas) that differ from a single or deterministic contributing area. Groundwater flow was simulated in the surficial deposits and the underlying bedrock in the 47-square-mile study area. Observations (165 groundwater levels and 7 base flows) provided sufficient information to estimate parameters representing recharge and horizontal hydraulic conductivity of the glacial deposits and hydraulic conductance of streambeds. The calibrated value for recharge to valley-fill deposits was 27.3 inches per year (in/yr) and to upland till deposits was 18.7 in/yr. Calibrated values for horizontal hydraulic conductivity of the valley-fill deposits ranged from 20 to 480 feet per day (ft/d) and of the upland till deposits was 16.2 ft/d. Calibrated values of streambed hydraulic conductance ranged from 10,000 to 52,000 feet squared per day. Values of recharge and horizontal hydraulic conductivity of the valley-fill deposits were the most precisely estimated, whereas the horizontal hydraulic conductivity of till deposits was the least precisely estimated. Simulated areas contributing recharge to the well centers on the basis of the calibrated model ranged from 0.19 to 1.12 square miles (mi 2 ) and covered a total area of 2.79 mi 2 for average well center withdrawal rates during 2004-08 (235 to 1,858 gallons per minute (gal/min)). Simulated areas contributing recharge for the maximum well center pumping capacities (800 to 8,500 gal/min) ranged from 0.37 to 3.53 mi2 and covered a total area of 7.99 mi2 in the modeled area. Simulated areas contributing recharge extend upgradient of the well centers to upland till and to groundwater divides. Some areas contributing recharge include small, isolated areas remote from the well centers. Relatively short groundwater traveltimes from recharging locations to discharging wells indicated the wells are vulnerable to contamination from land-surface activities: median traveltimes ranged from 2.9 to 5.0 years for the well centers, and 78 to 93 percent of the traveltimes were 10 years or less for the well centers. Land cover in the areas contributing recharge includes a substantial amount of urban land use for the two well centers in the Hunt River Basin, agriculture and sand and gravel mining uses for the well center in the Annaquatucket River Basin, and, for the well center in the Pettaquamscutt River Basin, land use is primarily undeveloped. Model-prediction uncertainty was evaluated using a Monte Carlo analysis. The parameter variance-covariance matrix from nonlinear regression was used to create parameter sets that reflect the uncertainty of the parameter estimates and the correlation among parameters. The remaining parameters representing the glacial deposits (vertical anisotropy of valley-fill deposits and of till deposits, maximum groundwater evapotranspiration, and hydraulic conductance for headdependent cells representing a groundwater divide) that could not be estimated with nonlinear regression were incorporated into the variance-covariance matrix using prior information on parameters. Thus the uncertainty analysis was an outcome of calibrating the parameters to available observations and to information that the modeler provided. A water budget and model-fit statistical criteria were used to assess parameter sets so that prediction uncertainty was not overestimated. Because of the effects of parameter uncertainty, the size of the probabilistic contributing areas for each well center for both average and maximum pumping rates was larger than the size of the deterministic contributing areas for the well center. Thus, some areas not in the deterministic contributing area may actually be in the contributing area, including additional areas of urban and agricultural land use. Generally, areas closest to the well centers with short groundwater traveltimes are associated with higher probabilities, whereas areas distant from the well centers with long groundwater traveltimes are associated with lower probabilities. The deterministic contributing areas generally corresponded to areas associated with high probabilities (greater than 50 percent). Areas associated with low probabilities extended long distances along groundwater divides in the uplands remote from the well centers.

Rhode Island↗

USGS2010

Delineation and Prediction Uncertainty of Areas Contributing Recharge to Selected Well Fields in Wetland and Coastal Settings, Southern Rhode Island

Areas contributing recharge to four well fields in two study sites in southern Rhode Island were delineated on the basis of steady-state groundwater-flow models representing average hydrologic conditions. The wells are screened in sand and gravel deposits in wetland and coastal settings. The groundwater-flow models were calibrated by inverse modeling using nonlinear regression. Summary statistics from nonlinear regression were used to evaluate the uncertainty associated with the predicted areas contributing recharge to the well fields. In South Kingstown, two United Water Rhode Island well fields are in Mink Brook watershed and near Worden Pond and extensive wetlands. Wetland deposits of peat near the well fields generally range in thickness from 5 to 8 feet. Analysis of water-level drawdowns in a piezometer screened beneath the peat during a 20-day pumping period indicated vertical leakage and a vertical hydraulic conductivity for the peat of roughly 0.01 ft/d. The simulated area contributing recharge for average withdrawals of 2,138 gallons per minute during 2003-07 extended to groundwater divides in mostly till and morainal deposits, and it encompassed 2.30 square miles. Most of a sand and gravel mining operation between the well fields was in the simulated contributing area. For the maximum pumping capacity (5,100 gallons per minute), the simulated area contributing recharge expanded to 5.54 square miles. The well fields intercepted most of the precipitation recharge in Mink Brook watershed and in an adjacent small watershed, and simulated streams ceased to flow. The simulated contributing area to the well fields included an area beneath Worden Pond and a remote, isolated area in upland till on the opposite side of Worden Pond from the well fields. About 12 percent of the pumped water was derived from Worden Pond. In Charlestown, the Central Beach Fire District and the East Beach Water Association well fields are on a small (0.85 square mile) peninsula in a coastal setting. The wells are screened in a coarse-grained, ice-proximal part of a morphosequence with saturated thicknesses generally less than 30 feet on the peninsula. The simulated area contributing recharge for the average withdrawal (16 gallons per minute) during 2003-07 was 0.018 square mile. The contributing area extended southwestward from the well fields to a simulated groundwater mound; it underlay part of a small nearby wetland, and it included isolated areas on the side of the wetland opposite the well fields. For the maximum pumping rate (230 gallons per minute), the simulated area contributing recharge (0.26 square mile) expanded in all directions; it included a till area on the peninsula, and it underlay part of a nearby pond. Because the well fields are screened in a thin aquifer, simulated groundwater traveltimes from recharge locations to the discharging wells were short: 94 percent of the traveltimes were 10 years or less, and the median traveltime was 1.3 years. Model-prediction uncertainty was evaluated using a Monte Carlo analysis; the parameter variance-covariance matrix from nonlinear regression was used to create parameter sets for the analysis. Important parameters for model prediction that could not be estimated by nonlinear regression were incorporated into the variance-covariance matrix. For the South Kingstown study site, observations provided enough information to constrain the uncertainty of these parameters within realistic ranges, but for the Charlestown study site, prior information on parameters was required. Thus, the uncertainty analysis for the South Kingstown study site was an outcome of calibrating the model to available observations, but the Charlestown study site was also dependent on information provided by the modeler. A water budget and model-fit statistical criteria were used to assess parameter sets so that prediction uncertainty was not overestimated. For the scenarios using maximum pumping rates at both study

Scientific Investigations Report↗

USGS2007

Simulation of Ground-Water Flow and Areas Contributing Recharge to Production Wells in Contrasting Glacial Valley-Fill Settings, Rhode Island

Areas contributing recharge and sources of water to a production well field in the Village of Harrisville and to a production well field in the Town of Richmond were delineated on the basis of calibrated, steady-state ground-water-flow models representing average hydrologic conditions. The study sites represent contrasting glacial valley-fill settings. The area contributing recharge to a well is defined as the surface area where water recharges the ground water and then flows toward and discharges to the well. In Harrisville, the production well field is composed of three wells in a narrow, approximately 0.5-mile-wide, valley-fill setting on opposite sides of Batty Brook, a small intermittent stream that drains 0.64 square mile at its confluence with the Clear River. Glacial stratified deposits are generally less areally extensive than previously published. The production wells are screened in a thin (30 feet) but transmissive aquifer. Paired measurements of ground-water and surface-water levels indicated that the direction of flow between the brook and the aquifer was generally downward during pumping conditions. Long-term mean annual streamflow from two streams upgradient of the well field totaled 0.72 cubic feet per second. The simulated area contributing recharge for the 2005 average well-field withdrawal rate of 224 gallons per minute extended upgradient to ground-water divides in upland areas and encompassed 0.17 square mile. The well field derived 62 percent of pumped water from intercepted ground water and 38 percent from infiltrated stream water from the Batty Brook watershed. For the maximum simulated well-field withdrawal of 600 gallons per minute, the area contributing recharge expanded to 0.44 square mile to intercept additional ground water and infiltration of stream water; the percentage of water derived from surface water, however, was the same as for the average pumping rate. Because of the small size of Batty Brook watershed, most of the precipitation recharge in the watershed was withdrawn by the well field at the maximum rate either by intercepted ground water or indirectly by infiltrated stream water. Because the production wells are screened in a thin and transmissive aquifer in a small watershed, simulated ground-water traveltimes from recharge locations to the discharging wells were relatively short: 93 percent of the traveltimes were 10 years or less. In Richmond, the production well field is composed of two wells adjacent to and east of the Wood River in a moderately broad, approximately 1.2-mile-wide, valley-fill setting. The wells are screened in a transmissive aquifer with saturated thickness greater than 60 feet. Streamflow measurements in Baker Brook, a tributary to the Wood River 0.4 mile north of the well-field site, indicated that natural net loss of streamflow between the upland-valley contact and a downstream site was 0.12 cubic feet per second under average hydrologic conditions. Simulated areas contributing recharge for the maximum well-field pumping rate of 675 gallons per minute and for one-half the maximum rate extended northeastward from the well field to ground-water divides in upland areas. The area contributing recharge also included a remote, isolated area on the opposite side of the Wood River from the well field. The model simulation indicated that the well field did not derive any of its water from the Wood River because of the large watershed and associated quantity of ground water available for capture by the well field. The area contributing recharge for one-half the maximum rate was 0.31 square mile and the primary source of water to the well field was direct precipitation recharge. Fifteen percent of the water withdrawn from the production wells, however, was obtained from Baker Brook, indicating the importance of even small, distant tributary streams to the contributing area to a well. The area contributing recharge on the opposite side of the Wood River is

Scientific Investigations Report↗

USGS2004

Delineation of areas contributing recharge to selected public-supply wells in Glacial Valley-Fill and Wetland Settings, Rhode Island

Areas contributing recharge and sources of water to one proposed and seven present public-supply wells, screened in sand and gravel deposits and clustered in three study areas, were determined on the basis of calibrated, steady-state ground-water-flow models representing average hydrologic conditions. The area contributing recharge to a well is defined as the surface area where water recharges the ground water and then flows toward and discharges to the well. In Cumberland and Lincoln, public-supply well fields on opposite sides of the Blackstone River are in a narrow valley bordered by steep hillslopes. Ground-water-level and river-stage measurements indicated that river water was infiltrating the aquifer and flowing toward the wells during pumping conditions. Simulated areas contributing recharge to the Cumberland well field operating alone for both average (324 gallons per minute) and maximum (1,000 gallons per minute) pumping rates extend on both sides of the river to the lateral model boundaries, which is the contact between the valley and uplands. The area contributing recharge at the average pumping rate is about 0.05 square mile and the well field derives 72 percent of pumped water from upland runoff. At the maximum pumping rate, the area contributing recharge extends farther up and down the valley to 0.12 square mile and the primary source of water to the well field was infiltrated river water (53 percent). Upland areas draining toward the areas contributing recharge encompass 0.58 and 0.66 square mile for the average and maximum rates, respectively. By incorporating the backup Lincoln well-field withdrawals (2,083 gallons per minute) into the model, the area contributing recharge to the Cumberland well field operating at its maximum rate is reduced to 0.08 square mile; part of the simulated area which contributes recharge to the Cumberland well field when it is operating alone contributes instead to the Lincoln well field when both well fields are pumped. The Cumberland well field compensates by increasing the percentage of water it withdraws from the river by 11 percent. The upland area draining toward the Cumberland contributing area is 0.55 square mile. The area contributing recharge to the Lincoln well field is 0.08 square mile and infiltrated river water contributes 88 percent of the total water; the upland area draining toward the contributing area is 0.34 square mile. In North Smithfield, a public-supply well in a valley-fill setting is close to Trout Brook Pond, which is an extension of the Lower Slatersville Reservoir. A comparison of water levels from the pond and underlying sediments indicates that water is not infiltrated from Trout Brook Pond when the supply well is pumped at its maximum rate of 200 gallons per minute. Simulated areas contributing recharge for the maximum pumping rate and for the estimated maximum yield, 500 gallons per minute, of a proposed replacement well extend to the ground-water divides on both sides of Trout Brook Pond. For the 200 gallons-per-minute rate, the area contributing recharge is 0.23 square mile; the well derives almost all of its water from intercepted ground water that normally discharges to surface-water bodies. For the pumping rate of 500 gallons per minute, the area contributing recharge is 0.45 square mile. The increased pumping rate is balanced by additional intercepted ground water and by inducing 25 percent of the total withdrawn water from surface water. In Westerly, one public-supply well is in a watershed where the primarily hydrologic feature is a wetland. Water levels in piezometers surrounding the well site indicated a downward vertical gradient and the potential for water in the wetland to infiltrate the underlying aquifer. The simulated area contributing recharge for the average pumping rate (240 gallons per minute) and for the maximum pumping rate (700 gallons per minute) extends to the surrounding uplands (surficial materials not covered by t

Scientific Investigations Report↗

USGS2001

Pond-aquifer interaction at South Pond of Lake Cochituate, Natick, Massachusetts

A U.S. Army facility on a peninsula in South Pond of Lake Cochituate was designated a Superfund site by the U.S. Environmental Protection Agency in 1994 because contaminated ground water was detected at the facility, which is near the Natick Springvale public-supply wellfield. The interaction between South Pond and the underlying aquifer controls ground-water flow patterns near the pond and determines the source of water withdrawn from the wellfield.A map of the bathymetry and the thickness of fine-grained pond-bottom sediments was prepared on the basis of fathometer, ground-penetrating radar, and continuous seismic-reflection surveys. The geophysical data indicate that the bottom sediments are fine grained toward the middle of the pond but are coarse grained in shoreline areas. Natick Springvale wellfield, which consists of three active public-supply wells adjacent to South Pond, is 2,200 feet downgradient from the boundary of the Army facility. That part of South Pond between the Natick Springvale wellfield and the Army facility is 18 feet deep with at least 14 feet of fine-grained sediment beneath the pond-bottom. Water levels from the pond and underlying sediments indicate a downward vertical gradient and the potential for infiltration of pond water near the wellfield. Head differences between the pond and the wellfield ranged from 1.66 to 4.41 feet during this study. The velocity of downward flow from South Pond into the pond-bottom sediments, determined on the basis of temperature profiles measured over a diurnal cycle at two locations near the wellfield, was 0.5 and 1.0 feet per day. These downward velocities resulted in vertical hydraulic conductivities of 1.1 and 2.9 feet per day for the pond-bottom sediments.Naturally occurring stable isotopes of oxygen and hydrogen were used as tracers of pond water and ground water derived from recharge of precipitation, two potential sources of water to a well in a pond-aquifer setting. The isotopic composition of pond water varied seasonally and was distinctly different from the isotopic composition of ground water. The isotopic composition of shallow water beneath and adjacent to South Pond near the wellfield corresponds to the temporal variation of pond water, indicating that nearly all water at shallow depths was derived from pond water. A two-component mixing model based on the average stable isotope values of the source waters indicated that 64 ?15 percent at the 95-percent confidence interval of the water withdrawn at the public-supply wells was derived from the pond; pond water accounted for most of the uncertainty in the result. The rate of infiltration of pond water into the aquifer and discharging to the wellfield was 1.0 million gallons per day at the average pumping rate.

Water-Resources Investigations Report↗

USGS2001

Geohydrology and limnology of Walden Pond, Concord, Massachusetts

The trophic ecology and ground-water contributing area of Walden Pond, in Concord and Lincoln, Mass., were investigated by the U.S. Geological Survey in cooperation with the Massachusetts Department of Environmental Management from April 1997 to July 2000. Bathymetric investigation indicated that Walden Pond (24.88 hectares), a glacial kettle-hole lake with no surface inlet or outlet, has three deep areas. The maximum depth (30.5 meters) essentially was unchanged from measurements made by Henry David Thoreau in 1846. The groundwater contributing area (621,000 square meters) to Walden Pond was determined from water-table contours in areas of stratified glacial deposits and from land-surface contours in areas of bedrock highs. Walden Pond is a flow-through lake: Walden Pond gains water from the aquifer along its eastern perimeter and loses water to the aquifer along its western perimeter. Walden Pond contributing area also includes Goose Pond and its contributing area. A water budget calculated for Walden Pond, expressed as depth of water over the lake surface, indicated that 45 percent of the inflow to the lake was from precipitation (1.215 meters per year) and 55 percent from ground water (1.47 meters per year). The groundwater inflow estimate was based on the average of two different approaches including an isotope mass-balance approach. Evaporation accounted for 26 percent of the outflow from the lake (0.71 meters per year) whereas lake-water seepage to the groundwater system contributed 74 percent of the outflow (1.97 meters per year). The water-residence time of Walden Pond is approximately 5 years. Potential point sources of nutrients to ground water, the Concord municipal landfill and a trailer park, were determined to be outside the Walden Pond groundwater contributing area. A third source, the septic leach field for the Walden Pond State Reservation facilities, was within the groundwater contributing area. Nutrient budgets for the lake indicated that nitrogen inputs (858 kilograms per year) were dominated (30 percent) by plume water from the septic leach field and, possibly, by swimmers (34 percent). Phosphorus inputs (32 kilograms per year) were dominated by atmospheric dry deposition, background ground water, and estimated swimmer inputs. Swimmer inputs may represent more than 50 percent of the phosphorus load during the summer. The septic-system plume did not contribute phosphorus, but increased the nitrogen to phosphorus ratio for inputs from 41 to 59, on an atom-to-atom basis. The ratio of nitrogen to phosphorus in input loads and within the lake indicated algal growth would be strongly phosphorus limited. Nitrogen supply in excess of plant requirements may mitigate against nitrogen fixing organisms including undesirable blooms of cyanobacteria. Based on areal nutrient loading, Walden Pond is a mesotrophic lake. Hypolimnetic oxygen demand of Walden Pond has increased since a profile was measured in 1939. Currently (1999), the entire hypolimnion of Walden Pond becomes devoid of dissolved oxygen before fall turnover in late November; whereas historical data indicated dissolved oxygen likely remained in the hypolimnion during 1939. The complete depletion of dissolved oxygen likely causes release of phosphorus from the sediments. Walden Pond contains a large population of the deep-growing benthic macro alga Nitella, which has been hypothesized to promote water clarity in other clear-water lakes by sequestering nutrients and keeping large areas of the sediment surface oxygenated. Loss of Nitella populations in other lakes has correlated with a decline in water quality. Although the Nitella standing crop is large in Walden Pond, Nitella still appears to be controlled by nutrient availability. Decreasing phosphorus inputs to Walden Pond, by amounts under anthropogenic control would likely contribute to the stability of the Nitella population in the metalimnion, may reverse oxygen depletion in the hypolimnion, and decreas

Water-Resources Investigations Report↗

USGS2001

Hydrology and trophic ecology of Walden Pond, Concord, Massachusetts

Explore the source record for details and available documents.

Water-Resources Investigations Report↗

USGS2000

Methods for estimating low-flow statistics for Massachusetts streams

Methods and computer software are described in this report for determining flow duration, low-flow frequency statistics, and August median flows. These low-flow statistics can be estimated for unregulated streams in Massachusetts using different methods depending on whether the location of interest is at a streamgaging station, a low-flow partial-record station, or an ungaged site where no data are available. Low-flow statistics for streamgaging stations can be estimated using standard U.S. Geological Survey methods described in the report. The MOVE.1 mathematical method and a graphical correlation method can be used to estimate low-flow statistics for low-flow partial-record stations. The MOVE.1 method is recommended when the relation between measured flows at a partial-record station and daily mean flows at a nearby, hydrologically similar streamgaging station is linear, and the graphical method is recommended when the relation is curved. Equations are presented for computing the variance and equivalent years of record for estimates of low-flow statistics for low-flow partial-record stations when either a single or multiple index stations are used to determine the estimates. The drainage-area ratio method or regression equations can be used to estimate low-flow statistics for ungaged sites where no data are available. The drainage-area ratio method is generally as accurate as or more accurate than regression estimates when the drainage-area ratio for an ungaged site is between 0.3 and 1.5 times the drainage area of the index data-collection site. Regression equations were developed to estimate the natural, long-term 99-, 98-, 95-, 90-, 85-, 80-, 75-, 70-, 60-, and 50-percent duration flows; the 7-day, 2-year and the 7-day, 10-year low flows; and the August median flow for ungaged sites in Massachusetts. Streamflow statistics and basin characteristics for 87 to 133 streamgaging stations and low-flow partial-record stations were used to develop the equations. The streamgaging stations had from 2 to 81 years of record, with a mean record length of 37 years. The low-flow partial-record stations had from 8 to 36 streamflow measurements, with a median of 14 measurements. All basin characteristics were determined from digital map data. The basin characteristics that were statistically significant in most of the final regression equations were drainage area, the area of stratified-drift deposits per unit of stream length plus 0.1, mean basin slope, and an indicator variable that was 0 in the eastern region and 1 in the western region of Massachusetts. The equations were developed by use of weighted-least-squares regression analyses, with weights assigned proportional to the years of record and inversely proportional to the variances of the streamflow statistics for the stations. Standard errors of prediction ranged from 70.7 to 17.5 percent for the equations to predict the 7-day, 10-year low flow and 50-percent duration flow, respectively. The equations are not applicable for use in the Southeast Coastal region of the State, or where basin characteristics for the selected ungaged site are outside the ranges of those for the stations used in the regression analyses. A World Wide Web application was developed that provides streamflow statistics for data collection stations from a data base and for ungaged sites by measuring the necessary basin characteristics for the site and solving the regression equations. Output provided by the Web application for ungaged sites includes a map of the drainage-basin boundary determined for the site, the measured basin characteristics, the estimated streamflow statistics, and 90-percent prediction intervals for the estimates. An equation is provided for combining regression and correlation estimates to obtain improved estimates of the streamflow statistics for low-flow partial-record stations. An equation is also provided for combining regression and drainage-area ratio estimates to obtain improved e

Massachusetts↗