USGS ScienceSearch

Geology topics

Carl S. Carlson

Publications and source records attributed to Carl S. Carlson.

At least 19 recordsLinked to original sources

An automated geographic information system-based hydraulic modeling tool for developing preliminary culvert designs for stream crossings in Massachusetts

Introduction Currently (2026), many of the about 25,000 roadway crossing structures over rivers and streams in Massachusetts are undersized. Undersized culverts and bridges can be detrimental to fish and wildlife movement, habitat continuity, and the health of aquatic organisms. Undersized culverts also can lack the resiliency needed to withstand large floods, which could be worsened by potential increases in flood magnitude and frequency due to climate change. Improving culvert and bridge designs for stream crossing projects may improve aquatic organism passage, stream continuity, and resiliency during future floods by decreasing upstream overbank flooding, road flooding and erosion, and degradation of aquatic habitat. The U.S. Geological Survey (USGS), Massachusetts Department of Environmental Protection (MassDEP), and University of Massachusetts Amherst began a series of cooperative studies in July 2019 to develop an automated geographic information system (GIS) hydraulic modeling tool for preliminary culvert designs for stream crossings. The USGS plans to provide preliminary culvert designs in the web-based StreamStats application, which enables municipalities and engineers to view potential designs and related information for stream crossing replacement projects in Massachusetts. This application can (a) provide information on hydrology, hydraulics, and ecological conditions at stream crossing sites, (b) provide users with potential culvert designs to improve aquatic organism passage and flood resiliency, and (c) assist MassDEP in implementing the Massachusetts Wetlands Protection Act regulations for stream crossing projects.

Massachusetts

Three-dimensional distribution of residence time metrics in the glaciated United States using metamodels trained on general numerical models

Residence time distribution (RTD) is a critically important characteristic of groundwater flow systems; however, it cannot be measured directly. RTD can be inferred from tracer data with analytical models (few parameters) or with numerical models (many parameters). The second approach permits more variation in system properties but is used less frequently than the first because large‐scale numerical models can be resource intensive. Using a novel automated approach, a set of 115 inexpensive general simulation models (GSMs) was used to create RTD metrics (fraction of young groundwater, defined as < 65 years old; mean travel time of young fraction; median travel time of old fraction; and mean path length). GSMs captured the general trends in measured tritium concentrations in 431 wells. Boosted Regression Tree metamodels were trained to predict these RTD metrics using available wall‐to‐wall hydrogeographic digital sets as explanatory features. The metamodels produced a three‐dimensional distribution of predictions throughout the glacial system that generally matched with the numerical model RTD metrics. In addition to the expected importance of aquifer thickness and recharge rate in predicting RTD metrics, two new data sets, Multi‐Order Hydrologic Position (MOHP) and hydrogeologic terrane were important predictors. These variables by themselves produced metamodels with Nash‐Sutcliffe efficiency close to the full metamodel. Metamodel predictions showed that the volume of young groundwater stored in the glaciated U.S. is about 6,000 km 3 , or about 0.5% of globally stored young groundwater.

Water Resources Research

Development of simulated groundwater-contributing areas to selected streams, ponds, coastal water bodies, and production wells in the Plymouth-Carver region and Cape Cod, Massachusetts

Introduction The U.S. Geological Survey (USGS), in support of the Massachusetts Estuaries Project (MEP), delineated groundwater-contributing areas to various hydrologic receptors including ponds, streams, and coastal water bodies throughout southeastern Massachusetts, including portions of the Plymouth-Carver aquifer system and all of Cape Cod. These contributing areas were delineated over a 6-year period from 2003 through 2008 by using previously published regional USGS groundwater-flow models for the Plymouth-Carver region (Masterson and others, 2009), the Sagamore (western) and Monomoy (eastern) flow lenses of Cape Cod (Walter and Whealan, 2005), and lower Cape Cod (Masterson, 2004). The original USGS groundwater-contributing areas were subsequently revised in some locations by the MEP to remove modeling artifacts or to make the contributing areas more consistent with site-specific hydrologic conditions without further USGS review. This report describes the process used to create the USGS groundwater-contributing areas and provides these model results in their original format in a single, publicly accessible publication.

Massachusetts

Groundwater and surface-water interaction and effects of pumping in a complex glacial-sediment aquifer, phase 2, east-central Massachusetts

The U.S. Geological Survey, in cooperation with the Town of Framingham, Massachusetts, has investigated the potential of proposed groundwater withdrawals at the Birch Road well site to affect nearby surface water bodies and wetlands, including Lake Cochituate, the Sudbury River, and the Great Meadows National Wildlife Refuge in east-central Massachusetts. In 2012, the U.S. Geological Survey developed a Phase 1 numerical groundwater model of a complex glacial-sediment aquifer to synthesize hydrogeologic information and simulate potential future pumping scenarios. The model was developed with MODFLOW-NWT, an updated version of a standard USGS numerical groundwater flow modeling program that improves solution of unconfined groundwater flow problems. The groundwater model and investigations of the aquifer improved understanding of groundwater&ndash;surface-water interaction and the effects of groundwater withdrawals on surface-water bodies and wetlands in the study area. The initial work also revealed a need for additional information and model refinements to better understand this complex aquifer system. In this second phase of the study, the original groundwater flow model was revised to improve representation of groundwater and surface-water hydrology, stabilize the model, and reduce model error. The model was simplified by reducing the number of layers from 5 to 3 and adding the MODFLOW lake package (LAK) to simulate Lake Cochituate and Pod Meadow Pond and better represent interaction between the lakes and the aquifer. Model revisions improved stability and shortened run times, allowing use of automated parameter estimation software (PEST) to further refine the model hydraulic parameters and reduce simulation errors. Model simulations indicate that under average base-flow conditions, the Birch Road wells have a small effect on flow in the Sudbury River during most months, even at the maximum pumping rate of 4.9 ft 3 /s (3.17 Mgal/d). Maximum percent streamflow depletion in the Sudbury River caused by simulated pumping takes place during simulated drought conditions, when streamflow decreased by as much as 21 percent under maximum continuous pumping. Simulations also indicate that groundwater withdrawals at the Birch Road site could be managed so that adverse streamflow impacts are substantially ameliorated. Under the most ecologically conservative simulated drought conditions, simulated streamflow depletion was reduced from 21 percent to 3 percent by pumping at the maximum rate for 6 months rather than for 12 months. Simulations that return 10 percent of the Birch Road well withdrawals to Pod Meadow Pond indicate a modest reduction in the Sudbury River streamflow depletion and provide a larger percentage increase to streamflow just downstream of the pond. The groundwater model also indicates that well locations can have a large effect on the sustainable pumping rate and so should be chosen carefully. The model provides a tool for evaluating alternative pumping rates and schedules not included in this analysis.

Massachusetts

Simulated responses of streams and ponds to groundwater withdrawals and wastewater return flows in southeastern Massachusetts

Water use, such as withdrawals, wastewater return flows, and interbasin transfers, can alter streamflow regimes, water quality, and the integrity of aquatic habitat and affect the availability of water for human and ecosystem needs. To provide the information needed to determine alteration of streamflows and pond water levels in southeastern Massachusetts, existing groundwater models of the Plymouth-Carver region and western (Sagamore flow lens) and eastern (Monomoy flow lens) Cape Cod were used to delineate subbasins and simulate long-term average and average monthly streamflows and pond levels for a series of water-use conditions. Model simulations were used to determine the extent to which streamflows and pond levels were altered by comparing simulated streamflows and pond levels under predevelopment conditions with streamflows and pond levels under pumping only and pumping with wastewater return flow conditions. The pumping and wastewater return flow rates used in this study are the same as those used in previously published U.S. Geological Survey studies in southeastern Massachusetts and represent the period from 2000 to 2005. Streamflow alteration for the nontidal portions of streams in southeastern Massachusetts was evaluated within and at the downstream outlets of 78 groundwater subbasins delineated for this study. Evaluation of streamflow alteration at subbasin outlets is consistent with the approach used by the U.S. Geological Survey for the topographically derived subbasins in the rest of Massachusetts. The net effect of pumping and wastewater return flows on streamflows and pond levels varied by location and included no change in areas minimally affected by water use, decreases in areas affected more by pumping than by wastewater return flows, or increases in areas affected more by wastewater return flows than by pumping. Simulated alterations to long-term average streamflows at subbasin outlets in response to pumping with wastewater return flows were within about 10 percent of predevelopment streamflows for most of the subbasins in the study area. Alterations ranged from a decrease (depletion) of 43.9 percent at an unnamed tributary to Salt Pond in the Plymouth-Carver region to an increase (surcharge) of 18.2 percent at an unnamed tributary to the Centerville River on western Cape Cod. In general, the relative effects of pumping and wastewater return flows typically were larger in the subbasins with low streamflows than in the subbasins with high streamflows, and there were more depleted streamflows than surcharged streamflows. Increases in streamflows in response to wastewater return flows were generally largest in subbasins with a high density of septic systems or a centralized wastewater treatment facility. For average monthly conditions, streamflow alteration results were similar spatially to results for long-term average conditions. However, differences in the extent of alteration by month were observed; percentage streamflow depletions in most subbasins typically were greatest during the low-streamflow months of August and October. The percentages of the total number of ponds affected by pumping with wastewater return flows under long-term average conditions in the modeled areas were 28 percent for the Plymouth-Carver region, 67 percent for western Cape Cod, and 75 percent for eastern Cape Cod. Pond-level alterations ranged from a decrease of 4.6 feet at Great South Pond in the Plymouth Carver region to an increase of 0.9 feet at Wequaquet Lake in western Cape Cod. The magnitudes of monthly alterations to pond water levels were fairly consistent throughout the year.

Massachusetts

Simulation of nitrogen attenuation in a subterranean estuary, representative of the southern coast of Cape Cod, Massachusetts

A two-dimensional model was developed by the U.S. Geological Survey, in cooperation with the U.S. Environmental Protection Agency, to assess flow and chemical reaction associated with groundwater discharge through the subterranean estuary representative of coastal salt ponds of southern Cape Cod. The model simulated both the freshwater and saltwater flow systems and accounted for density-dependent flow, tidal fluctuation, and chemical reactivity among oxygen, dissolved organic carbon, nitrate, and ammonia. Not previously incorporated into one model, the interaction of these effects can now be simulated in the subterranean estuary context. An analysis of the flow system under mean-tide conditions was conducted first to provide the initial conditions for a subsequent analysis that included the effects of tidal fluctuations. Tidal fluctuations were simulated with a repeated couplet that represented a high tide-low tide sequence and alternating locations of head-dependent flux boundaries placed along the simulated seabed, above and below the levels of the respective high and low tides. Boundary conditions for chemical species included nitrate in recharge, and oxygen and organic matter (including organic nitrogen) in infiltrating solutions of head-dependent boundaries. Reaction chemistry was limited to oxidative degradation of organic matter (including remineralization of ammonia) with oxygen or nitrate as electron acceptors and nitrification of ammonia in the presence of oxygen. Simulations using the SEAWAT-2000 computer program resulted in two mixing zones-between freshwater and saltwater in a deep saltwater wedge and in an intertidal salt zone, which results from tidal fluctuation. The mixing zones are the principal locations where nitrogen attenuation reactions occurred-between organic matter in the saltwater zones of the aquifer and nitrate in the freshwater zone. In mean-tide PHT3D model simulations, 15 percent of nitrogen that is recharged was attenuated because of reaction with dissolved organic matter, a denitrification reaction that reduces nitrate to nitrogen gas. When a fluctuating tide was simulated, the amount of recharged nitrogen that was denitrified increased to 20 percent. Chemical reaction was controlled by the rate of mixing of freshwater and saltwater, which contained the reactants nitrate and dissolved organic matter, respectively, necessary for nitrogen attenuation reactions to take place. Reaction occurred in both the deep saltwater wedge and in an increased denitrification. However, mixing may also have been enhanced partly by numerical dispersion.

Massachusetts

Development of a numerical model to simulate groundwater flow in the shallow aquifer system of Assateague Island, Maryland and Virginia

A three-dimensional groundwater-flow model was developed for Assateague Island in eastern Maryland and Virginia to simulate both groundwater flow and solute (salt) transport to evaluate the groundwater system response to sea-level rise. The model was constructed using geologic and spatial information to represent the island geometry, boundaries, and physical properties and was calibrated using an inverse modeling parameter-estimation technique. An initial transient solute-transport simulation was used to establish the freshwater-saltwater boundary for a final calibrated steady-state model of groundwater flow. This model was developed as part of an ongoing investigation by the U.S. Geological Survey Climate and Land Use Change Research and Development Program to improve capabilities for predicting potential climate-change effects and provide the necessary tools for adaptation and mitigation of potentially adverse impacts.

Open-File Report

Geohydrology, geochemistry, and groundwater simulation (1992-2011) and analysis of potential water-supply management options, 2010-60, of the Langford Basin, California

Groundwater withdrawals began in 1992 from the Langford Basin within the Fort Irwin National Training Center (NTC), California. From April 1992 to December 2010, approximately 12,300 acre-feet of water (averaging about 650 acre-feet per year) has been withdrawn from the basin and transported to the adjacent Irwin Basin. Since withdrawals began, water levels in the basin have declined by as much as 40 feet, and the quality of the groundwater withdrawn from the basin has deteriorated. The U.S. Geological Survey collected geohydrologic data from Langford Basin during 1992–2011 to determine the quantity and quality of groundwater available in the basin. Geophysical surveys, including gravity, seismic refraction, and time-domain electromagnetic induction surveys, were conducted to determine the depth and shape of the basin, to delineate depths to the Quaternary-Tertiary interface, and to map the depth to the water table and changes in water quality. Data were collected from existing wells and test holes, as well as 11 monitor wells that were installed at 5 sites as part of this study. Water-quality samples collected from wells in the basin were used to determine the groundwater chemistry within the basin and to delineate potential sources of poor-quality groundwater. Analysis of stable isotopes of oxygen and hydrogen in groundwater indicates that present-day precipitation is not a major source of recharge to the basin. Tritium and carbon-14 data indicate that most of the basin was recharged prior to 1952, and the groundwater in the basin has an apparent age of 12,500 to 30,000 years. Recharge to the basin, estimated to be less than 50 acre-feet per year, has not been sufficient to replenish the water that is being withdrawn from the basin. A numerical groundwater-flow model was developed for the Langford Basin to better understand the aquifer system used by the Fort Irwin NTC as part of its water supply, and to provide a tool to help manage groundwater resources at the NTC. Measured groundwater-level declines since the initiation of withdrawals (1992–2011) were used to calibrate the groundwater-flow model. The simulated recharge was about 46 acre-feet per year, including approximately 6 acre-feet per year of natural recharge derived from precipitation runoff and as much as 40 acre-feet per year of underflow from the Irwin Basin. Between April 1992 and December 2010, an average of about 650 acre-feet per year of water was withdrawn from the Langford Basin. Groundwater withdrawals in excess of natural recharge resulted in a net loss of 11,670 acre-feet of groundwater storage within the basin for the simulation period. The Fort Irwin NTC is considering various groundwater-management options to address the limited water resources in the Langford Basin. The calibrated Langford Basin groundwater-flow model was used to evaluate the hydrologic effects of four groundwater-withdrawal scenarios being considered by the Fort Irwin NTC over the next 50 years (January 2011 through December 2060). Continuation of the 2010 withdrawal rate in the three existing production wells will result in 70 feet of additional drawdown in the central part of the basin. Redistributing the 2010 withdrawal rate equally to the three existing wells and two proposed new wells in the northern and southern parts of the basin would result in about 10 feet less drawdown in the central part of the basin but about 100 feet of additional drawdown in the new well in the northern part of the basin and about 50 feet of additional drawdown in the new well in the southern part of the basin. Reducing the withdrawals from the three existing production wells in the central part of the basin from about 45,000 acre-feet to about 32,720 acre-feet would result in about 40 feet of additional drawdown in the central basin near the pumping wells, about 25 feet less than if withdrawals were not reduced. The combination of reducing and redistributing the cumulative withdrawals to the three existing and two proposed new wells results in about 40 feet of additional drawdown in the central and southern parts of the basin and about 70 feet in the northern part of the basin. These results show that reducing and redistributing the groundwater withdrawals would maintain the upper aquifer at greater than 50 percent of its predevelopment saturated thickness throughout the groundwater basin. The scenarios simulated for this study demonstrate how the calibrated model can be utilized to evaluate the hydrologic effects of different water-management strategies.

California

Simulation of groundwater and surface-water interaction and effects of pumping in a complex glacial-sediment aquifer, east central Massachusetts

The effects of groundwater pumping on surface-water features were evaluated by use of a numerical groundwater model developed for a complex glacial-sediment aquifer in northeastern Framingham, Massachusetts, and parts of surrounding towns. The aquifer is composed of sand, gravel, silt, and clay glacial-fill sediments up to 270 feet thick over an irregular fractured bedrock surface. Surface-water bodies, including Cochituate Brook, the Sudbury River, Lake Cochituate, Dudley Pond, and adjoining wetlands, are in hydraulic connection with the aquifer and can be affected by groundwater withdrawals. Groundwater and surface-water interaction was simulated with MODFLOW-NWT under current conditions and a variety of hypothetical pumping conditions. Simulations of hypothetical pumping at reactivated water supply wells indicate that captured groundwater would decrease baseflow to the Sudbury River and induce recharge from Lake Cochituate. Under constant (steady-state) pumping, induced groundwater recharge from Lake Cochituate was equal to about 32 percent of the simulated pumping rate, and flow downstream in the Sudbury River decreased at the same rate as pumping. However, surface water responded quickly to pumping stresses. When pumping was simulated for 1 month and then stopped, streamflow depletions decreased by about 80 percent within 2 months and by about 90 percent within about 4 months. The fast surface water response to groundwater pumping offers the potential to substantially reduce streamflow depletions during periods of low flow, which are of greatest concern to the ecological integrity of the river. Results indicate that streamflow depletion during September, typically the month of lowest flow, can be reduced by 29 percent by lowering the maximum pumping rates to near zero during September. Lowering pumping rates for 3 months (July through September) reduces streamflow depletion during September by 79 percent as compared to constant pumping. These results demonstrate that a seasonal or streamflow-based groundwater pumping schedule can reduce the effects of pumping during periods of low flow.

Massachusetts

Effects of water use and land use on streamflow and aquatic habitat in the Sudbury and Assabet River Basins, Massachusetts

Water withdrawals from surface-water reservoirs and groundwater have affected streamflow in the Sudbury and Assabet River Basins. These effects are particularly evident in the upper Sudbury River Basin, which prompted the need to improve the understanding of water resources and aquatic habitat in these basins. In 2004, the U.S. Geological Survey, in cooperation with the Massachusetts Department of Conservation and Recreation, developed a precipitation-runoff model that uses Hydrologic Simulation Program-FORTRAN (HSPF) to evaluate the effects of water use and projected future water-use and land-use change on streamflow. As part of this study, the aquatic habitat in the basins and the effects of streamflow alteration also were evaluated. Chapter 1 of the report covers the development of the HSPF model that focuses on the upper Sudbury River Basin (106 square miles) but covers the entire Sudbury and Assabet River Basins (339 square miles). The model was calibrated to an 11-year period (1993-2003) using observed or estimated streamflow at four streamgages. The model was then used to simulate long-term (1960-2004) streamflows to evaluate the effects of average 1993-2003 water use and projected 2030 water-use and land-use change over long-term climatic conditions. Simulations indicate that the average 1993-2003 withdrawals most altered streamflow relative to no withdrawals in small headwater subbasins where the ratios of mean annual withdrawals to mean annual streamflow are the highest. The effects of withdrawals are also appreciable in other parts of the upper Sudbury River Basin as a result of the perpetuation of the effects of large withdrawals in upstream reaches or in subbasins that also have a high ratio of withdrawal to streamflow. The simulated effects of potential 2030 water-use and land-use change indicate small decreases in flows as a result of increased water demands, but these flow alterations were offset as a result of decreased evapotranspiration associated with the loss of deep-rooted vegetation. Simulations of reactivating production wells near the north end of Lake Cochituate indicate pumping could substantially affect lake levels and flows at the lake outlet or in nearby reaches in the Sudbury River during periods of low flow, but the effects vary depending on the source of the water to the wells, which is largely unknown. Chapter 2 of the report covers the fish-community assessment and comparison of streamflow-setting standards for protecting aquatic habitat. The fish-community assessment indicates the main stems of the Sudbury and Assabet Rivers are dominated by macrohabitat generalists. Water temperatures recorded in seven free-flowing reaches in the upper Sudbury River Basin at three sites unaffected by withdrawals or impoundments are generally suitable for cold-water fish; however, summer temperatures often rose to a level considered critical to long-term survival of brook trout. At four sites downstream from withdrawals or reservoirs, or both, summer water temperatures were often in the upper critical range for brook trout survival. Physically and statistically based methods for determining streamflows for protecting aquatic habitat were applied at 10 selected riffle sites in the Sudbury and Assabet River Basins. Physically based methods, R2Cross and Wetted-Perimeter, use site-specific physical and hydraulic information and a one-dimensional hydraulics model, HEC-RAS, to determine flows that meet the criteria set forth by the method. The median flow that meets 2-of-3 of the R2Cross hydraulic criteria (percentage of bankfull wetted perimeter, average velocity, and mean depth) ranged from about 0.07 to 0.72 cubic feet per second per square mile (ft3/s/mi2) with an overall median of about 0.24 ft3/s/mi2; the median Wetted-Perimeter target flow ranged from about 0.10 to 0.51 ft3/s/mi2 with an overall median of about 0.25 ft3/s/mi2. Statistically based methods?Tennant, New England Aquatic Base Flow (ABF)

Scientific Investigations Report

Hydrogeology and numerical simulation of the unconsolidated glacial aquifer in the Pootatuck River Basin, Newtown, Connecticut

A study of the groundwater and stream-aquifer interaction in the Pootatuck River Basin, Newtown, Connecticut, was conducted to analyze the effect of production wells on the groundwater levels and streamflow in the Pootatuck River as part of a cooperative program between the U.S. Geological Survey and Newtown, Connecticut. This study will help address concerns about the increasing competition for water for human uses and protection of aquatic habitat. The groundwater-flow model developed in the study was designed for use as a tool to assist planners in assessing the effects of potential future development, which will change the amount and distribution of recharge available to the groundwater system. Several different techniques were used to investigate the interconnection between the stream and the aquifer. Temperature, groundwater levels, stream stage, and stable-isotope data collected during aquifer tests at the principal production wells in the Pootatuck River Basin, as well as groundwater-flow simulations of the system, indicate that more than half of the water pumped from the wells comes from the Pootatuck River. This finding potentially has a large effect on approaches for protecting the water quality of the pumped water. Increases in the amount of impervious surface from future development will reduce and redistribute recharge to the groundwater system. The simulation of future development scenarios showed a decrease in the simulated base flow in the main stem of the Pootatuck River and in all of the 26 simulated subbasins, with some of the subbasins showing a decrease of more than 20 percent when new development had 85 percent impervious area. The groundwater-flow model and particle tracking were used to determine areas that contribute recharge to the five production wells available for use in the Pootatuck River Basin. These areas included narrow portions of the aquifer that extended beyond the immediate upgradient areas, probably because of deeper groundwater-flow paths.

Scientific Investigations Report

Hydrologic conditions and a firm-yield assessment for J.B. Converse Lake, Mobile County, Alabama, 1991-2006

J.B. Converse (Converse) Lake is the primary source of drinking water for the city of Mobile, Alabama. Concerns regarding the ability of the reservoir to meet current and future water demands during drought conditions have prompted this study. The 1991 through 2006 water years included a drought that occurred during 2000, and drought conditions currently (2007) are affecting the area. To assist officials of the Mobile Area Water and Sewer System in planning for future demands for drinking water in the Mobile metropolitan area, the firm yield for Converse Lake was estimated by the U.S. Geological Survey. The firm yield of Converse Lake was estimated using the Massachusetts Department of Environmental Protection’s firm-yield-estimator (FYE) model, which recently was refined by the U.S. Geological Survey. The model uses a mass-balance approach to determine the maximum average daily withdrawal rate that can be sustained during a period of record that includes a drought of record. If the reservoir is in contact with an aquifer, the FYE also includes routines that estimate the volume of ground-water and surface-water exchange between the aquifer and the reservoir. The average daily firm yield for Converse Lake was estimated to be 79 million gallons per day using the FYE routine that does not include ground-water exchange between the reservoir and the adjacent aquifer. Observed lake levels and withdrawals during the drought of 2000 indicate that more than 74 million gallons per day of water were withdrawn without complete depletion of reservoir storage. Therefore, it is likely that ground-water exchange with the reservoir may supplement available reservoir storage. If water exchange occurs between the aquifer and the reservoir, an increase in the volume of water available to the reservoir may occur during a drought. To quantify the potential ground-water contribution to reservoir storage, an analytical solution was applied to the FYE simulation of Converse Lake to estimate ground-water exchange between the reservoir and the aquifer. Aquifer properties required by the FYE were estimated by model calibration to observed water levels that occurred during the drought of 2000. When ground-water exchange between the reservoir and adjacent aquifer is included, the average daily firm yield increased to 83 million gallons per day. The estimate of 83 million gallons per day incorporates both total surface-water flow and ground-water exchange components. This analysis indicated that direct ground-water interaction contributes about 5 percent of the firm yield of Converse Lake. However, the average daily firm yield of 83 million gallons per day, based in part on calibrated values for aquifer transmissivity and storage, can be used only as a guideline until these aquifer properties can be defined better by field investigation in the Converse Lake watershed.

Alabama

Hydrogeology and Simulation of Groundwater Flow in the Plymouth-Carver-Kingston-Duxbury Aquifer System, Southeastern Massachusetts

The glacial sediments that underlie the Plymouth-Carver-Kingston-Duxbury area of southeastern Massachusetts compose an important aquifer system that is the primary source of water for a region undergoing rapid development. Population increases and land-use changes in this area has led to two primary environmental effects that relate directly to groundwater resources: (1) increases in pumping that can adversely affect environmentally sensitive groundwater-fed surface waters, such as ponds, streams, and wetlands; and (2) adverse effects of land use on the quality of water in the aquifer. In response to these concerns, the U.S. Geological Survey, in cooperation with the Massachusetts Department of Environmental Protection, began an investigation in 2005 to improve the understanding of the hydrogeology in the area and to assess the effects of changing pumping and recharge conditions on groundwater flow in the Plymouth-Carver-Kingston-Duxbury aquifer system. A numerical flow model was developed based on the USGS computer program MODFLOW-2000 to assist in the analysis of groundwater flow. Model simulations were used to determine water budgets, flow directions, and the sources of water to pumping wells, ponds, streams, and coastal areas. Model-calculated water budgets indicate that approximately 298 million gallons per day (Mgal/d) of water recharges the Plymouth-Carver-Kingston-Duxbury aquifer system. Most of this water (about 70 percent) moves through the aquifer, discharges to streams, and then reaches the coast as surface-water discharge. Of the remaining 30 percent of flow, about 25 percent of the water that enters the aquifer as recharge discharges directly to coastal areas and 5 percent discharges to pumping wells. Groundwater withdrawals are anticipated to increase from the current (2005) rate of about 14 Mgal/d to about 21 Mgal/d by 2030. Pumping from large-capacity production wells decreases water levels and increases the potential for effects on surface-water bodies, which are affected by pumping and wastewater disposal locations and rates. Pumping wells that are upgradient of surface-water bodies potentially capture water that would otherwise discharge to these surface-water bodies, thereby reducing streamflow and pond levels. The areas most affected by proposed increases in groundwater withdrawals are in the Towns of Plymouth and Wareham where more than half of the proposed increase in pumping will occur. In response to an increase of about 7 Mgal/d of pumping, groundwater discharge to streams is reduced by about 6 cubic feet per second (ft3/s) (about 4 Mgal/d) from a total of about 325 ft3/s. Reduction in streamflow is moderated by an increase of artificial recharge from wastewater returned to the aquifer by onsite domestic septic systems and centralized wastewater treatment facilities. It is anticipated that about 3 Mgal/d of the 7 Mgal/d of increase in pumped water will be returned to the aquifer as wastewater by 2030. Currently (2005) about 3 percent of groundwater discharge to streams is from wastewater return flow to the aquifer during average conditions. During drought conditions, the component of streamflow augmented by wastewater return flow doubles as wastewater recharge remains constant and aquifer recharge rates decrease. Wastewater return flow, whether as direct groundwater discharge to streams or as an additional source of aquifer recharge, increases the height of the water table near streams, thereby moderating the effects of increased groundwater withdrawals on streamflow. An analysis of a simulated drought similar to the 1960s drought of record indicates that the presence of streams moderates the effects on water levels of reduced aquifer recharge. The area where water-table altitudes were least affected by drought was in the Weweantic River watershed in the Town of Carver. Water levels decreased by less than 2 feet from current average conditions compared to decreases of greater than 5

Scientific Investigations Report

Characteristics of the April 2007 Flood at 10 Streamflow-Gaging Stations in Massachusetts

A large 'nor'easter' storm on April 15-18, 2007, brought heavy rains to the southern New England region that, coupled with normal seasonal high flows and associated wet soil-moisture conditions, caused extensive flooding in many parts of Massachusetts and neighboring states. To characterize the magnitude of the April 2007 flood, a peak-flow frequency analysis was undertaken at 10 selected streamflow-gaging stations in Massachusetts to determine the magnitude of flood flows at 5-, 10-, 25-, 50-, 100-, 200-, and 500-year return intervals. The magnitude of flood flows at various return intervals were determined from the logarithms of the annual peaks fit to a Pearson Type III probability distribution. Analysis included augmenting the station record with longer-term records from one or more nearby stations to provide a common period of comparison that includes notable floods in 1936, 1938, and 1955. The April 2007 peak flow was among the highest recorded or estimated since 1936, often ranking between the 3d and 5th highest peak for that period. In general, the peak-flow frequency analysis indicates the April 2007 peak flow has an estimated return interval between 25 and 50 years; at stations in the northeastern and central areas of the state, the storm was less severe resulting in flows with return intervals of about 5 and 10 years, respectively. At Merrimack River at Lowell, the April 2007 peak flow approached a 100-year return interval that was computed from post-flood control records and the 1936 and 1938 peak flows adjusted for flood control. In general, the magnitude of flood flow for a given return interval computed from the streamflow-gaging station period-of-record was greater than those used to calculate flood profiles in various community flood-insurance studies. In addition, the magnitude of the updated flood flow and current (2008) stage-discharge relation at a given streamflow-gaging station often produced a flood stage that was considerably different than the flood stage indicated in the flood-insurance study flood profile at that station. Equations for estimating the flow magnitudes for 5-, 10-, 25-, 50-, 100-, 200-, and 500-year floods were developed from the relation of the magnitude of flood flows to drainage area calculated from the six streamflow-gaging stations with the longest unaltered record. These equations produced a more conservative estimate of flood flows (higher discharges) than the existing regional equations for estimating flood flows at ungaged rivers in Massachusetts. Large differences in the magnitude of flood flows for various return intervals determined in this study compared to results from existing regional equations and flood insurance studies indicate a need for updating regional analyses and equations for estimating the expected magnitude of flood flows in Massachusetts.

Scientific Investigations Report

Simulated Effects of Year 2030 Water-Use and Land-Use Changes on Streamflow near the Interstate-495 Corridor, Assabet and Upper Charles River Basins, Eastern Massachusetts

Continued population growth and land development for commercial, industrial, and residential uses have created concerns regarding the future supply of potable water and the quantity of ground water discharging to streams in the area of Interstate 495 in eastern Massachusetts. Two ground-water models developed in 2002-2004 for the Assabet and Upper Charles River Basins were used to simulate water supply and land-use scenarios relevant for the entire Interstate-495 corridor. Future population growth, water demands, and commercial and residential growth were projected for year 2030 by the Metropolitan Area Planning Council. To assess the effects of future development on subbasin streamflows, seven scenarios were simulated by using existing computer-based ground-water-flow models with the data projected for year 2030. The scenarios incorporate three categories of projected 2030 water- and land-use data: (1) 2030 water use, (2) 2030 land use, and (3) a combination of 2030 water use and 2030 land use. Hydrologic, land-use, and water-use data from 1997 through 2001 for the Assabet River Basin study and 1989 through 1998 for the Upper Charles River Basin study were used to represent current conditions - referred to as 'basecase' conditions - in each basin to which each 2030 scenario was compared. The effects of projected 2030 land- and water-use change on streamflows in the Assabet River Basin depended upon the time of year, the hydrologic position of the subbasin in the larger basin, and the relative areas of new commercial and residential development projected for a subbasin. Effects of water use and land use on streamflow were evaluated by comparing average monthly nonstorm streamflow (base flow) for March and September simulated by using the models. The greatest decreases in streamflow (up to 76 percent in one subbasin), compared to the basecase, occurred in September, when streamflows are naturally at their lowest level. By contrast, simulated March streamflows decreased less than 6.5 percent from basecase streamflows in all subbasins for all scenarios. The simulations showed similar effects in the Upper Charles River Basin, but increased water use contributed to decreased simulated streamflow in most subbasins. Simulated changes in March streamflows for 2030 in the Upper Charles River Basin were within +- 6 percent of the basecase for all scenarios and subbasins. Percentage decreases in simulated September streamflows for 2030 were greater than in March but less than the September decreases that resulted for some subbasins in the Assabet River Basin. Only two subbasins of the Upper Charles River Basin had projected decreases greater than 5 percent. In the Mill River subbasin, the decrease was 11 percent, and in the Mine Brook subbasin, 6.6 percent. Changes in water use and wastewater return flow generally were found to have the greatest effect in the summer months when streamflow and aquifer recharge rates are low and water use is high. September increases in main-stem streamflow of both basins were due mainly to increased discharge of treated effluent from wastewater-treatment facilities on the main-stem rivers. In the Assabet River Basin, wastewater-treatment-facility discharge became a smaller proportion of total streamflow with distance downstream. In contrast, wastewater-treatment facility discharge in the Upper Charles River Basin became a greater proportion of streamflow with distance downstream. The effects of sewer-line extension and low-impact development on streamflows in two different subbasins of the Assabet River Basin also were simulated. The result of extending sewer lines with a corresponding decrease in septic-system return flow caused September streamflows to decrease as much as 15 percent in the Fort Pond Brook subbasin. The effect of low-impact development was simulated in the Hop Brook subbasin in areas projected for commercial development. In this simulation, the greater the area where low-i

Scientific Investigations Report

Ground-Water Contributions to Reservoir Storage and the Effect on Estimates of Firm Yield for Reservoirs in Massachusetts

Potential ground-water contributions to reservoir storage were determined for nine reservoirs in Massachusetts that had shorelines in contact with sand and gravel aquifers. The effect of ground water on firm yield was not only substantial, but furthermore, the firm yield of a reservoir in contact with a sand and gravel aquifer was always greater when the ground-water contribution was included in the water balance. Increases in firm yield ranged from 2 to 113 percent, with a median increase in firm yield of 10 percent. Additionally, the increase in firm yield in two reservoirs was greater than 85 percent. This study identified a set of equations that are based on an analytical solution to the ground-water-flow equation for the case of one-dimensional flow in a finite-width aquifer bounded by a linear surface-water feature such as a stream. These equations, which require only five input variables, were incorporated into an existing firm-yield-estimator (FYE) model, and the potential effect of ground water on firm yield was evaluated. To apply the FYE model to a reservoir in Massachusetts, the model requires that the drainage area to the reservoir be clearly defined and that some surface water flows into the reservoir. For surface-water-body shapes having a more realistic representation of a reservoir shoreline than a stream, a comparison of ground-water-flow rates simulated by the ground-water equations with flow rates simulated by a two-dimensional, finite-difference ground-water-flow model indicate that the agreement between the simulated flow rates is within ?10 percent when the ratio of the distance from the reservoir shoreline to the aquifer boundary to the length of shoreline in contact with the aquifer is between values of 0.5 and 3.5. Idealized reservoir-aquifer systems were assumed to verify that the ground-water-flow equations were implemented correctly into the existing FYE model; however, the modified FYE model has not been validated through a comparison of simulated and observed data. A comparison of simulated and observed reservoir water levels would further define limitations to the applicability of the ground-water-flow equations to reservoirs in Massachusetts whose shorelines are in contact with a sand and gravel aquifer.

Scientific Investigations Report

Ground-water levels near the top of the water-table mound, western Cape Cod, Massachusetts, 2002-04

In January 2002 the U.S. Geological Survey began continuous water-level monitoring in three wells in the vicinity of the Southeast Ranges of Camp Edwards, near the Impact Area of the Massachusetts Military Reservation on Cape Cod. The purpose of this effort was to examine how water levels at sites with different unsaturated-zone thicknesses near the top of the water-table mound beneath western Cape Cod are affected by temporally variable recharge from precipitation, which is the sole source of water to the sand and gravel aquifer. The depths to water at the well sites are about 18, 30, and 101 feet below land surface. This report presents the first 3 years of water-level records and an estimate of aquifer recharge calculated from climatological measurements by the Jensen and Haise method and the Thornthwaite method. The water levels in the three wells varied temporally by about 4.5 feet during the study period. A comparison of the water levels with those measured in a nearby monitoring well with about 42 years of monthly measurements indicates that the 3-year monitoring period included the lowest water levels on western Cape Cod since the drought of the 1960's. The response of water levels to recharge was related to the depth to water. Water levels in the two wells with shallow depths to water responded quickly (within hours or days) to recharge, whereas the water-level response in the well with the greatest depth to water often lagged the recharge event by a month or more. The variations in the water levels among the wells changed as the location of the top of the water-table mound moved with the changing water-table altitude.

Massachusetts

Simulated ground-water flow for a pond-dominated aquifer system near Great Sandy Bottom Pond, Pembroke, Massachusetts

A ground-water flow simulation for a 66.4-square-mile area around Great Sandy Bottom (GSB) Pond (105 acres) near Pembroke, Massachusetts, was developed for use by local and State water managers to assess the yields for public water supply of local ponds and wells for average climatic and drought conditions and the effects of water withdrawals on nearby water levels and streamflows. Wetlands and ponds cover about 30 percent of the study area and the aquifer system is dominated by interactions between ground water and the ponds. The three largest surface-water bodies in the study area are Silver Lake (640 acres), Monponsett Pond (590 acres), and Oldham Pond (236 acres). The study area is drained by tributaries of the Taunton River to the southwest, the South and North Rivers to the northeast, and the Jones River to the southeast. In 2002, 10.8 million gallons per day of water was exported from ponds and 3.5 million gallons per day from wells was used locally for public supply. A transient ground-water-flow model with 69 monthly stress periods spanning the period from January 1998 through September 2003 was calibrated to stage at GSB Pond and nearby Silver Lake and streamflow and water levels collected from September 2002 through September 2003. The calibrated model was used to assess hydrologic responses to a variety of water-use and climatic conditions. Simulation of predevelopment (no pumping or export) average monthly (1949-2002) water-level conditions caused the GSB Pond level to increase by 6.3 feet from the results of a simulation using average 2002 pumping for all wells, withdrawals, and exports. Most of this decline can be attributed to pumping, withdrawals, and exports of water from sites away from GSB Pond. The effects of increasing the export rate from GSB Pond by 1.25 and 1.5 times the 2002 rate were a lowering of pond levels by a maximum of 1.6 and 2.8 feet, respectively. Simulated results for two different drought conditions, one mild drought similar to that of 1979-82 and a more severe drought similar to that of 1963-66, but with current (2002) pumping, were compared to results for average monthly recharge conditions (1949-2002). Simulated mild drought conditions showed a reduction of GSB Pond level of about 1.3 feet and a lower streamflow of about 1.7 percent in the nearby stream. Simulated severe drought conditions reduced the pond level at GSB Pond by almost 7 feet and lowered streamflow by about 37 percent. Varying cranberry-irrigation practices had little effect on simulated GSB Pond water levels, but may be important in other ponds. The model was most sensitive to changes in areal recharge. An increase and decrease of 22 percent in recharge produced changes in the GSB Pond water level of +1.4 feet and -2.4 feet, respectively. The accuracy of simulation results was best in the central portion of the study area in the immediate location of GSB Pond. The model was developed with the study-area boundary far enough away from the GSB Pond area that the boundary would have minimal effect on the water levels in GSB Pond, nearby ponds, and the underlying aquifer system. The model is best suited for use by local and State water managers to assess the effects of different withdrawal scenarios for wells and ponds near GSB Pond and for general delineation of areas contributing recharge to wells and ponds in the vicinity of GSB Pond. The model in its current form may not be well suited to detailed analyses of water budgets and flow patterns for parts of the study area farther from GSB Pond without further investigation, calibration, and data collection.

Scientific Investigations Report