Search USGSSearch

SEARCH · Search USGS

Results for “Contribution”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 recordsLinked to original sources

Elevated CO 2 enhances biological contributions to elevation change in coastal wetlands by offsetting stressors associated with sea-level rise

1. Sea-level rise, one indirect consequence of increasing atmospheric CO 2 , poses a major challenge to long-term stability of coastal wetlands. An important question is whether direct effects of elevated CO 2 on the capacity of marsh plants to accrete organic material and to maintain surface elevations outweigh indirect negative effects of stressors associated with sea-level rise (salinity and flooding). 2. In this study, we used a mesocosm approach to examine potential direct and indirect effects of atmospheric CO 2 concentration, salinity and flooding on elevation change in a brackish marsh community dominated by a C 3 species, Schoenoplectus americanus, and a C 4 grass, Spartina patens. This experimental design permitted identification of mechanisms and their role in controlling elevation change, and the development of models that can be tested in the field. 3. To test hypotheses related to CO 2 and sea-level rise, we used conventional anova procedures in conjunction with structural equation modelling (SEM). SEM explained 78% of the variability in elevation change and showed the direct, positive effect of S. americanus production on elevation. The SEM indicated that C 3 plant response was influenced by interactive effects between CO 2 and salinity on plant growth, not a direct CO 2 fertilization effect. Elevated CO 2 ameliorated negative effects of salinity on S. americanus and enhanced biomass contribution to elevation. 4. The positive relationship between S. americanus production and elevation change can be explained by shoot-base expansion under elevated CO 2 conditions, which led to vertical soil displacement. While the response of this species may differ under other environmental conditions, shoot-base expansion and the general contribution of C 3 plant production to elevation change may be an important mechanism contributing to soil expansion and elevation gain in other coastal wetlands. 5. Synthesis. Our results revealed previously unrecognized interactions and mechanisms contributing to marsh elevation change, including amelioration of salt stress by elevated CO 2 and the importance of plant production and shoot-base expansion for elevation gain. Identification of biological processes contributing to elevation change is an important first step in developing comprehensive models that permit more accurate predictions of whether coastal marshes will persist with continued sea-level rise or become submerged. ?? 2008 The Authors.

Journal of Ecology

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

Contributions of Great Salt Lake playa- and industrially-sourced priority pollutant metals in dust contribute to possible health hazards in the communities of northern Utah

Communities and ecosystems of northern Utah, USA receive particulate pollution from anthropogenic activity and dust emissions from sources including the Great Salt Lake (“the Lake”) playa. In addition to affecting communities, anthropogenic pollution is delivered to the Lake's playa sediments, which are eroded during dust events. Yet, spatial variability in dust flux and composition and their risks to human health are poorly understood. We analyzed dust in 17 passive samplers proximal to the Lake during fall 2022 for dust flux, the dust fraction of particulate matter, 87 Sr/ 86 Sr, and elemental geochemistry. We evaluated spatial patterns of 11 priority pollutant metals and estimated the hypothetical non-cancer dust and soil ingestion health hazard for six age cohorts. We observed the highest dust fluxes proximal to the Lake's playa. The highest concentrations of and greatest number of metals occurred in and south of Ogden, UT. Sites to the northeast of Farmington Bay had the highest fluxes. Metal concentrations and 87 Sr/ 86 Sr suggest that the dust composition near Bountiful represents contributions from anthropogenic sources, whereas the dust composition to the northeast of Farmington Bay reflects the Lake's playa emissions. Evaluations of potential health hazards from dust ingestion suggest that children between birth and 6 years are vulnerable at higher ingestion rates. Thallium, As, Pb, Co and Cr contributed most to the estimated hazard. Among these, As and sometimes Pb are likely derived from the Lake's playa emissions. Thus, suppression of dust emissions from the Lake's playa may decrease possible health risks for children in northern Utah.

Utah

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

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

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

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

Estimation of the recharge area contributing water to a pumped well in a glacial-drift, river-valley aquifer

The highly permeable, unconfined, glacial-drift aquifers that occupy most New England river valleys constitute the principal source of drinking water for many of the communities that obtain part or all of their public water supply from ground water. Recent events have shown that these aquifers are highly susceptible to contamination that results from a number of sources, such as seepage from wastewater lagoons, leaking petroleum-product storage tanks, and road salting. To protect the quality of water pumped from supply wells in these aquifers, it is necessary to ensure that potentially harmful contaminants do not enter the ground in the area that contributes water to the well. A high degree of protection can be achieved through the application of appropriate land-use controls within the contributing area. However, the contributing areas for most supply wells are not known. This report describes the factors that affect the size and shape of contributing areas to public supply wells and evaluates several methods that may be used to delineate contributing areas of wells in glacial-drift, river-valley aquifers. Analytical, two-dimensional numerical, and three-dimensional numerical models were used to delineate contributing areas. These methods of analysis were compared by applying them to a hypothetical aquifer having the dimensions and geometry of a typical glacial-drift, river-valley aquifer. In the model analyses, factors that control the size and shape of a contributing area were varied over ranges of values common to glacial-drift aquifers in New England. The controlling factors include the rate of well discharge, rate of recharge to the aquifer from precipitation and from adjacent till and bedrock uplands, distance of a pumping well from a stream or other potential source of induced recharge, degree of hydraulic connection of the aquifer with a stream, horizontal hydraulic conductivity of the aquifer, ratio of horizontal to vertical hydraulic conductivity, and degree of well penetration. Analytical methods proved easiest to apply but gave results that are considered to be less accurate than those obtainable by means of numerical-model analysis. Numerical models have the capability to more closely reflect the variable geohydrologic conditions typical of glacial-drift valley aquifers. For average conditions in the hypothetical aquifer, the analytical method predicts a contributing area limited to the well side of the river because a constant-head boundary simulated by image wells is used in the analytical model. For typical glacial-drift, river-valley aquifers, this simulation is unrealistic because drawdowns, caused by a pumping well, and the contributing area of the well can extend beneath and beyond a river or stream. A wide range of hydrologic conditions was simulated by using the two-dimensional numerical model. The resulting contributing area for a well pumped at 1.0 million gallons per day--a common pumping rate--ranged from about 0.9 to 1.8 square miles. Model analyses also show that the contributing area of pumped wells may be expected to extend to the opposite side of the river and to include significant areas of till uplands adjacent to the aquifer on both sides of the valley. Simulations done with the three-dimensional model allow a full three-dimensional delineation of the zone of contribution for a pumped well. For the relatively thin (100 feet or less) unconfined aquifers considered in this analysis, the three-dimensional model showed that the zone of contribution extended throughout the entire saturated thickness of aquifer; therefore, the two-dimensional simulations were considered adequate for delineating contributing areas in this particular hydrologic setting. For thicker aquifers, especially those having partially penetrating wells, three-dimensional models are preferable. Values for several of the factors that affect the size and shape of contributing recharge areas cannot be det

Water Supply Paper

Descriptions of anisotropy and heterogeneity and their effect on ground-water flow and areas of contribution to public supply wells in a karst carbonate aquifer system

Delineation of areas of contribution to wells tapping a karst carbonate aquifer system can be extremely difficult using conventional approaches designed for isotropic and homogeneous aquifers, because ground-water flow tends to be through solution-enhanced conduits. Nonradial flow along preferential zones can result in inaccurate estimates of flow paths and traveltimes. Because of the large variability in factors affecting contributing areas and an imperfect understanding of how these factors can vary, the estimation of contributing areas is an approximation at best. To better understand the effects of aquifer anisotropy and heterogeneity on areas of contribution, an exploratory modeling approach was used. MODFLOW, a numerical flow model, and MODPATH, a particle tracking program, were used to generate time-related areas of contribution for six hypothetical carbonate aquifer system types. The six types were conceptualized to approximate different types of aquifer anisotropy and heterogeneity. These include: (1) an isotropic and homogeneous single-layer system; (2) an anisotropic in a horizontal plane single layer system; (3) a discrete vertically fractured single-layer system; (4) a multilayered system; (5) a doubly porous single-layer system; and (6) a vertically and horizontally interconnected heterogeneous system. The simulated aquifer anisotropy was 5:1 (Kxx/Kyy) determined from TENSOR2D results. The simulated discrete vertical fracture network represents locations inferred from mapped photolineaments. The simulated enhanced flow zones were determined from borehole video and geophysical logs. Areas of contribution were simulated for two prototype regions. The two prototypes were selected to be representative of the hydrologic diversity within the study area and were designated the Central Swamp and Lake Terrace regions. Localized conditions in pumping, production well distribution, and aquifer transmissivity affect the size, shape, and orientation of areas of contribution to public supply wells. The simulated areas of contribution are 60 percent larger in the Central Swamp region where pumpage is more than double and transmissivity is about half that of the Lake Terrace region. Although these factors are important, this study focused on the effects from hydrogeologic factors common to karst carbonate aquifer systems. This study indicates that the distribution and type of aquifer anisotropy and heterogeneity will affect the size, shape, and orientation of areas of contribution in a karst carbonate aquifer system. The size of the 50-year time-related areas of contribution ranged from 8.2 to 39.1 square miles in the Central Swamp region and from 4.0 to 18.3 square miles in the Lake Terrace region. Simulations showed that the size of areas of contribution is primarily affected by simulated withdrawal rates, effective porosity of the carbonate rock, and transmissivity. The shape and orientation of the simulated areas of contribution primarily result from aquifer anisotropy, well distribution, flow along solution-enhanced zones, and short-circuiting of flow through fracture networks. Comparisons also were made between protection zones delineated using analytical models and areas of contribution delineated using numerical models. The size of the 5-year time-related protection zone in the Central Swamp region using an analytical model was almost twice as large as the numerically simulated area of contribution, and more than eight times larger than the numerically simulated area of contribution in the Lake Terrace region. The differences in size are primarily the result of how the flow field is approximated. The analytical method assumes only lateral flow to wells but numerical methods allow particles to move laterally and vertically. Additionally, multiple-well-interference effects resulting from the close proximity of several pumping wells cause individual capture zones to converge or diverge, depending on the difference in pumping rates and orientation among the wells. Such an interpretation is not available from analytical methods. The simulated distributions of aquifer anisotropy and heterogeneity, in this study, were highly conceptualized, but were based on plausible occurrences of anisotropy and heterogeneity inherent in carbonate aquifer systems.

Florida

Dissolved-solids contribution to the Colorado River from public lands in southeastern Nevada, through September 1993

The Bureau of Land Management administers about 9,300 square miles of public lands in southeastern Nevada that are part of the Colorado River Basin. The U.S. Geological Survey, in cooperation with the Bureau of Land Management, began a 5-year program in October 1988 to assess the contribution of dissolved solids to the fiver from those lands. About 6,200 square miles of public lands are in the Muddy River subbasin in Nevada. The estimated average dissolved-solids load contributed to the Colorado River from those lands was 28,000 tons per year from October 1988 through September 1993. Subsurface flow contributed about 86 percent (24,000 tons per year) of that load. About 730 square miles of public lands in the Las Vegas Wash subbasin contribute dissolved-solids load to the Colorado River. (About 120 square miles of public lands do not contribute to the river.) The estimated average dissolved-solids load contributed to the river from those lands was about 1,300 tons per year from October 1988 through September 1993. Subsurface flow contributed almost all of that load. About 1,100 square miles of public lands are in the Virgin River subbasin in Nevada. The estimated average dissolved- solids load contributed to the Colorado River from Nevada public lands in the subbasin was 8,700 tons per year. Subsurface flow contributed almost the entire load. About 1,200 square miles of Nevada public lands are in ephemeral tributaries that drain direcfly to the Colorado River or its impoundments (Lake Mead and Lake Mobave). The estimated average dissolved-solids load contributed to the river from those lands was 50 tons per year from surface runoff; however, the dissolved-solids load contributed by subsurface flow was not estimated. From October 1992 to September 1993, the Colorado River carried about 6,600,000 tons of dissolved solids past a streamflow gaging station 0.3 mile downstream from Hoover Dam. In contrast, surface runoff and subsurface flow contribute an estimated average dissolved-solids load of 38,000 tons per year from public lands in southeastern Nevada to the Colorado River. Land-management practices probably would not substantially reduce this contribution.

Water-Resources Investigations Report

Highly specialized recreationists contribute the most to the citizen science project eBird

Contributory citizen science projects (hereafter “contributory projects”) are a powerful tool for avian conservation science. Large-scale projects such as eBird have produced data that have advanced science and contributed to many conservation applications. These projects also provide a means to engage the public in scientific data collection. A common challenge across contributory projects like eBird is to maintain participation, as some volunteers contribute just a few times before disengaging. To maximize contributions and manage an effective program that has broad appeal, it is useful to better understand factors that influence contribution rates. For projects capitalizing on recreation activities (e.g., birding), differences in contribution levels might be explained by the recreation specialization framework, which describes how recreationists vary in skill, behavior, and motives. We paired data from a survey of birders across the United States and Canada with data on their eBird contributions (n = 28,926) to test whether those who contributed most are more specialized birders. We assigned participants to 4 contribution groups based on eBird checklist submissions and compared groups’ specialization levels and motivations. More active contribution groups had higher specialization, yet some specialized birders were not active participants. The most distinguishing feature among groups was the behavioral dimension of specialization, with active eBird participants owning specialized equipment and taking frequent trips away from home to bird. Active participants had the strongest achievement motivations for birding (e.g., keeping a life list), whereas all groups had strong appreciation motivations (e.g., enjoying the sights and sounds of birding). Using recreation specialization to characterize eBird participants can help explain why some do not regularly contribute data. Project managers may be able to promote participation, particularly by those who are specialized but not contributing, by appealing to a broader suite of motivations that includes both appreciation and achievement motivations, and thereby increase data for conservation.

Ornithological Applications

Use of a ground-water flow model to delineate contributing areas to the Puchack Well Field, Pennsauken township and vicinity, Camden county, New Jersey

The New Jersey Department of Environmental Protection (NJDEP) Well Head Protection Program, developed in response to the 1986 Federal Safe Drinking Water Act Amendments, requires delineation of Well Head Protection Areas (WHPA's), commonly called contributing areas, for all public and non-community water-supply wells in New Jersey. Typically, WHPA's for public community water-supply wells in New Jersey are delineated using a two-dimensional ground-water flow model incorporating the regional hydraulic gradient; however, NJDEP guidelines allow for the use of a three-dimensional flow model to delineate contributing areas to wells in complex hydrogeologic settings. The Puchack well field in Pennsauken Township, Camden County, N.J., is an area of strong hydraulic connection between the Lower aquifer of the Potomac-Raritan-Magothy aquifer system and the Delaware River. Interactions among and within the public-supply well fields in the area are complex. To delineate the contributing area to the Puchack well field, the U.S. Geological Survey, in cooperation with the NJDEP, developed an 11-layer ground-water flow model of the Potomac-Raritan-Magothy aquifer system in the Pennsauken Township area to simulate flow in the vicinity of the well field. The model incorporates the interaction between the aquifer system and the Delaware River, and includes boundary flows from an existing regional model of the Camden area. Recharge used in the model ranged from 4.5 to 14 inches per year, and horizontal hydraulic conductivity ranged from 50 to 250 feet per day. Values of vertical hydraulic conductivity ranging from 0.001 to 0.5 feet per day were assigned to zones created on the basis of variations in hydrogeologic conditions observed in geophysical logs from wells. A steady-state simulation was used to calibrate the model to synoptic water-level data collected in March 1998. Near the Puchack well field, simulated heads generally were within 1 foot of the measured heads in both the Middle and Lower aquifers. Simulated water-level differences across the confining units at most of the nested wells were within ? 0.5 feet of the differences calculated from measured water levels. The existing flow model was modified to meet NJDEP guidelines for delineating contributing areas in complex hydrogeologic settings. These modifications included rediscretizing the model grid to a finer grid and preparing the water-use data set for use in the rediscretized model. The contributing area to the Puchack well field was delineated by means of particle tracking. An uncertainty analysis was conducted in which 36 model-input parameters were both increased and decreased until the resulting change in simulated heads exceeded the model-calibration criterion of ? 5 feet at any model cell. Porosity most affected the size and shape of the contributing area. The distribution of withdrawals at the Morris/Delair well field and variations in recharge affected both the size and shape of contributing area to the Puchack well field and the source of water to the Puchack wells. The results of the uncertainty analysis were combined to determine the 'aggregate' contributing area to the Puchack well field--a composite of areas on the land surface that contributed flow to the Puchack well field in less than 12 years in any uncertainty simulation. The shape of the aggregate contributing area was most similar to that associated with a reduction in porosity, which indirectly affected the size and shape of the contributing areas by changing travel time.

Scientific Investigations Report

Case study for delineating a contributing area to a well in a fractured siliciclastic-bedrock aquifer near Lansdale, Pennsylvania

A supply well used by the North Penn Water Authority near Lansdale, Pa., was selected as a case study for delineating a contributing area in a fractured siliciclastic-bedrock aquifer. The study emphasized the importance of refining the understanding of factors that control ground-water movement to the well by conducting (1) geophysical logging and flow measurements, (2) ground-water level monitoring, (3) aquifer testing, and (4) geochemical sampling. This approach could be applicable for other wells in siliciclastic-bedrock terranes, especially those of Triassic age in southeastern Pennsylvania. The principal methods for refining the understanding of hydrology at supply well MG-1125 were aquifer testing, water-level measurements, and geophysical logging. Results of two constant-discharge aquifer tests helped estimate the transmissivity of water-producing units and evaluate the anisotropy caused by dipping beds. Results from slug tests provided estimates of transmissivity that were used to evaluate the results from the constant-discharge aquifer tests. Slug tests also showed the wide distribution of transmissivity, indicating that ground-water velocities must vary considerably in the well field. Water-level monitoring in observation wells allowed maps of the potentiometric surface near the well field to be drawn. The measurements also showed that the hydraulic gradient can change abruptly in response to pumping from nearby supply wells. Water levels measured at a broader regional scale in an earlier study also provided a useful view of the potentiometric surface for purposes of delineating the contributing area. Geophysical logging and measurements of flow within wells showed that about 60 percent of water from supply well MG-1125 probably is contributed from relatively shallow water-producing fractures from 60 to 125 feet below land surface, but measurable amounts of water are contributed by fractures to a depth of 311 feet below land surface. Chemical samples supported the evidence that shallow fractures probably contribute significant amounts of water to well MG-1125. The large contribution of water from shallow fractures indicates that the area providing part of the recharge to the well is not far removed from the wellhead. Preliminary delineations of the contributing area and the 100-day time-of travel area were computed from a water budget and time-of-travel equation. These delineations provided insight into the size (but not the shape) of the contributing areas. Three other approaches were used and results compared: (1) uniform-flow equation, (2) hydrogeologic mapping, and (3) numerical modeling. The uniform-flow equation predicted a contributing area that seemed unrealistic—extending far across the ground-water divide into an adjacent watershed. Hydrogeologic mapping, if used with the potentiometric surface and constrained by the water budget, produced contributing area that was similar to that from numerical modeling. Numerical modeling allowed the incorporation of anisotropy caused by dipping water-producing units, differing transmissivity values of geologic units, and ground-water withdrawals from nearby supply wells. The numerical modeling showed that groundwater withdrawals from nearby supply wells affected the contributing area to supply well MG-1125 but had less effect on the 100-day time-of-travel area.

Water-Resources Investigations Report

Hydrogeologic investigation and simulation of ground-water flow in the Upper Floridan Aquifer of north-central Florida and southwestern Georgia and delineation of contributing areas for selected city of Tallahassee, Florida, water-supply wells

A 4-year investigation of the Upper Floridan aquifer and ground-water flow system in Leon County, Florida, and surrounding counties of north-central Florida and southwestern Georgia began in 1990. The purpose of the investigation was to describe the ground-water flow system and to delineate the contributing areas to selected City of Tallahassee, Florida, water-supply wells. The investigation was prompted by the detection of low levels of tetrachloroethylene in ground-water samples collected from several of the city's water-supply wells. Hydrologic data and previous studies indicate that; ground-water flow within the Upper Floridan aquifer can be considered steady-state; the Upper Floridan aquifer is a single water-bearing unit; recharge is from precipitation; and that discharge occurs as spring flow, leakage to rivers, leakage to the Gulf of Mexico, and pumpage. Measured transmissivities of the aquifer ranged from 1,300 ft 2 /d (feet squared per day) to 1,300,000 ft2/d. Steady-state ground-water flow in the Upper Floridan aquifer was simulated using a three-dimensional ground- water flow model. Transmissivities ranging from less than 5,000 ft2/d to greater than 11,000,000 ft 2 /d were required to calibrate to observed conditions. Recharge rates used in the model ranged from 18.0 inches per year in areas where the aquifer was unconfined to less than 2 inches per year in broad areas where the aquifer was confined. Contributing areas to five Tallahassee water-supply wells were simulated by particle- tracking techniques. Particles were seeded in model cells containing pumping wells then tracked backwards in time toward recharge areas. The contributing area for each well was simulated twice, once assuming a porosity of 25 percent and once assuming a porosity of 5 percent. A porosity of 25 percent is considered a reasonable average value for the Upper Floridan aquifer; the 5 percent porosity simulated the movement of ground-water through only solution-enhanced bedding plains and fractures. The contributing areas were generally elliptical in shape, reflecting the influence of the sloping potentiometric surface. The contributing areas delineated for a 5 percent porosity were always much larger than those determined using a 25 percent porosity. The lowest average ground-water velocity computed within a contributing area, using a 25 percent porosity, was 1.0 ft/d (foot per day) and the highest velocity was 1.6 ft/d. The lowest average ground-water velocity, determined using a 5 percent porosity, was 2.4 ft/d and the highest was 7.4 ft/d. The contributing areas for each of the five wells was also determined analytically and compared to the model-derived areas. The upgradient width of the simulated contributing areas were larger than the upgradient width of the analytically determined contributing areas for four of the five wells. The model could more accurately delineate contributing areas because of the ability to simulate wells as partially penetrating and by incorporating complex, three-dimensional aquifer characteristics, which the analytical method could not.

Florida, Georgia

Population-specific life histories contribute to metapopulation viability

Restoration efforts can be improved by understanding how variations in life-history traits occur within populations of the same species living in different environments. This can be done by first understanding the demographic responses of natural occurring populations. Population viability analysis continues to be useful to species management and conservation with sensitivity analysis aiding in the understanding of population dynamics. In this study, using life-table response experiments and elasticity analyses, we investigated how population-specific life-history demographic responses contributed to the metapopulation viability of the Federally threatened Pitcher's thistle ( Cirsium pitcheri ). Specifically, we tested the following hypotheses: (1) Subpopulations occupying different environments within a metapopulation have independent demographic responses and (2) advancing succession results in a shift from a demographic response focused on growth and fecundity to one dominated by stasis. Our results showed that reintroductions had a positive contribution to the metapopulation growth rate as compared to native populations which had a negative contribution. We found no difference in succession on the contribution to metapopulation viability. In addition, we identified distinct population-specific contributions to metapopulation viability and were able to associate specific life-history demographic responses. For example, the positive impact of Miller High Dunes population on the metapopulation growth rate resulted from high growth contributions, whereas increased time of plant in stasis for the State Park Big Blowout population resulted in negative contributions. A greater understanding of how separate populations respond in their corresponding environment may ultimately lead to more effective management strategies aimed at reducing extinction risk. We propose the continued use of sensitivity analyses to evaluate population-specific demographic influences on metapopulation viability. In understanding the underlying causes of the projected extinction probabilities of each population and identifying broad-scale contributions of different populations to the metapopulation, the process of pinpointing target populations is simplified. More detailed analyses can then be applied to the target populations to increase population viability and consequently metapopulation viability. Based on our research, we suggest that the best approach to improve the overall metapopulation viability is to manage the contributions to population growth for each population separately.

Ecosphere

Estimating the per-capita contribution of habitats and pathways in a migratory network: A modelling approach

Every year, migratory species undertake seasonal movements along different pathways between discrete regions and habitats. The ability to assess the relative demographic contributions of these different habitats and pathways to the species’ overall population dynamics is critical for understanding the ecology of migratory species, and also has practical applications for management and conservation. Metrics for assessing habitat contributions have been well-developed for metapopulations, but an equivalent metric is not currently available for migratory populations. Here, we develop a framework for estimating the demographic contributions of the discrete habitats and pathways used by migratory species throughout the annual cycle by estimating the per capita contribution of cohorts using these locations. Our framework accounts for seasonal movements between multiple breeding and non-breeding habitats and for both resident and migratory cohorts. We illustrate our framework using a hypothetical migratory network of four habitats, which allows us to better understand how variations in habitat quality affect per capita contributions. Results indicate that per capita contributions for any habitat or pathway are dependent on habitat-specific survival probabilities in all other areas used as part of the migratory circuit, and that contribution metrics are spatially linked (e.g. reduced survival in one habitat also decreases the contribution metric for other habitats). Our framework expands existing theory on the dynamics of spatiotemporally structured populations by developing a generalized approach to estimate the habitat- and pathway-specific contributions of species migrating between multiple breeding and multiple non-breeding habitats for a range of life histories or migratory strategies. Most importantly, it provides a means of prioritizing conservation efforts towards those migratory pathways and habitats that are most critical for the population viability of migratory species.

Ecography

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

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