USGS Science⌕ Search

Geology topics

Stephen P. Garabedian

Publications and source records attributed to Stephen P. Garabedian.

10 recordsLinked to original sources
USGS2010

The Massachusetts Sustainable-Yield Estimator: A decision-support tool to assess water availability at ungaged stream locations in Massachusetts

Federal, State and local water-resource managers require a variety of data and modeling tools to better understand water resources. The U.S. Geological Survey, in cooperation with the Massachusetts Department of Environmental Protection, has developed a statewide, interactive decision-support tool to meet this need. The decision-support tool, referred to as the Massachusetts Sustainable-Yield Estimator (MA SYE) provides screening-level estimates of the sustainable yield of a basin, defined as the difference between the unregulated streamflow and some user-specified quantity of water that must remain in the stream to support such functions as recreational activities or aquatic habitat. The MA SYE tool was designed, in part, because the quantity of surface water available in a basin is a time-varying quantity subject to competing demands for water. To compute sustainable yield, the MA SYE tool estimates a daily time series of unregulated, daily mean streamflow for a 44-year period of record spanning October 1, 1960, through September 30, 2004. Selected streamflow quantiles from an unregulated, daily flow-duration curve are estimated by solving six regression equations that are a function of physical and climate basin characteristics at an ungaged site on a stream of interest. Streamflow is then interpolated between the estimated quantiles to obtain a continuous daily flow-duration curve. A time series of unregulated daily streamflow subsequently is created by transferring the timing of the daily streamflow at a reference streamgage to the ungaged site by equating exceedence probabilities of contemporaneous flow at the two locations. One of 66 reference streamgages is selected by kriging, a geostatistical method, which is used to map the spatial relation among correlations between the time series of the logarithm of daily streamflows at each reference streamgage and the ungaged site. Estimated unregulated, daily mean streamflows show good agreement with observed unregulated, daily mean streamflow at 18 streamgages located across southern New England. Nash-Sutcliffe efficiency goodness-of-fit values are between 0.69 and 0.98, and percent root-mean-square-error values are between 19 and 283 percent. The MA SYE tool provides an estimate of streamflow adjusted for current (2000-04) water withdrawals and discharges using a spatially referenced database of permitted groundwater and surface-water withdrawal and discharge volumes. For a user-selected basin, the database is queried to obtain the locations of water withdrawal or discharge volumes within the basin. Groundwater and surface-water withdrawals and discharges are subtracted and added, respectively, from the unregulated, daily streamflow at an ungaged site to obtain a streamflow time series that includes the effects of these withdrawals and discharges. Users also have the option of applying an analytical solution to the time-varying, groundwater withdrawal and discharge volumes that take into account the effects of the aquifer properties on the timing and magnitude of streamflow alteration. For the MA SYE tool, it is assumed that groundwater and surface-water divides are coincident. For areas of southeastern Massachusetts and Cape Cod where this assumption is known to be violated, groundwater-flow models are used to estimate average monthly streamflows at fixed locations. There are several limitations to the quality and quantity of the spatially referenced database of groundwater and surface-water withdrawals and discharges. The adjusted streamflow values do not account for the effects on streamflow of climate change, septic-system discharge, impervious area, non-public water-supply withdrawals less than 100,000 gallons per day, and impounded surface-water bodies.

Massachusetts↗
USGS2008

Sampling design for groundwater solute transport: Tests of methods and analysis of Cape Cod tracer test data

Tests of a one-dimensional sampling design methodology on measurements of bromide concentration collected during the natural gradient tracer test conducted by the U.S. Geological Survey on Cape Cod, Massachusetts, demonstrate its efficacy for field studies of solute transport in groundwater and the utility of one-dimensional analysis. The methodology was applied to design of sparse two-dimensional networks of fully screened wells typical of those often used in engineering practice. In one-dimensional analysis, designs consist of the downstream distances to rows of wells oriented perpendicular to the groundwater flow direction and the timing of sampling to be carried out on each row. The power of a sampling design is measured by its effectiveness in simultaneously meeting objectives of model discrimination, parameter estimation, and cost minimization. One-dimensional models of solute transport, differing in processes affecting the solute and assumptions about the structure of the flow field, were considered for description of tracer cloud migration. When fitting each model using nonlinear regression, additive and multiplicative error forms were allowed for the residuals which consist of both random and model errors. The one-dimensional single-layer model of a nonreactive solute with multiplicative error was judged to be the best of those tested. Results show the efficacy of the methodology in designing sparse but powerful sampling networks. Designs that sample five rows of wells at five or fewer times in any given row performed as well for model discrimination as the full set of samples taken up to eight times in a given row from as many as 89 rows. Also, designs for parameter estimation judged to be good by the methodology were as effective in reducing the variance of parameter estimates as arbitrary designs with many more samples. Results further showed that estimates of velocity and longitudinal dispersivity in one-dimensional models based on data from only five rows of fully screened wells each sampled five or fewer times were practically equivalent to values determined from moments analysis of the complete three-dimensional set of 29,285 samples taken during 16 sampling times.

Water Resources Research↗
USGS2008

Large-scale natural gradient tracer test in sand and gravel, Cape Cod, Massachusetts: 1. Experimental design and observed tracer movement

A large-scale natural gradient tracer experiment was conducted on Cape Cod, Massachusetts, to examine the transport and dispersion of solutes in a sand and gravel aquifer. The nonreactive tracer, bromide, and the reactive tracers, lithium and molybdate, were injected as a pulse in July 1985 and monitored in three dimensions as they moved as far as 280 m down-gradient through an array of multilevel samplers. The bromide cloud moved horizontally at a rate of 0.42 m per day. It also moved downward about 4 m because of density-induced sinking early in the test and accretion of areal recharge from precipitation. After 200 m of transport, the bromide cloud had spread more than 80 m in the direction of flow, but was only 14 m wide and 4–6 m thick. The lithium and molybdate clouds followed the same path as the bromide cloud, but their rates of movement were retarded about 50% relative to bromide movement because of sorption onto the sediments.

Massachusetts↗
USGS2008

Large-scale natural gradient tracer test in sand and gravel, Cape Cod, Massachusetts: 2. Analysis of spatial moments for a nonreactive tracer

A large-scale natural gradient tracer test was conducted to examine the transport of reactive and nonreactive tracers in a sand and gravel aquifer on Cape Cod, Massachusetts. As part of this test the transport of bromide, a nonreactive tracer, was monitored for about 280 m and quantified using spatial moments. The calculated mass of bromide for each sampling date varied between 85% and 105% of the injected mass using an estimated porosity of 0.39, and the center of mass moved at a nearly constant horizontal velocity of 0.42 m per day. A nonlinear change in the bromide longitudinal variance was observed during the first 26 m of travel distance, but afterward the variance followed a linear trend, indicating the longitudinal dispersivity had reached a constant value of 0.96 m. The transverse dispersivities were much smaller; transverse horizontal dispersivity was 1.8 cm, and transverse vertical dispersivity was about 1.5 mm.

Massachusetts↗
USGS2004

Delineation of Areas Contributing Water to the Dry Brook Public-Supply Well, South Hadley, Massachusetts

Areas contributing water to the Dry Brook public-supply well in South Hadley, Massachusetts, were delineated with a numerical ground-water-flow model that is based on geologic and hydrologic information for the confined sand and gravel aquifer pumped by the supply well. The study area is along the Connecticut River in central Massachusetts, about 12 miles north of Springfield, Massachusetts. Geologic units in the study area consist of Mesozoic-aged sedimentary and igneous bedrock, late-Pleistocene glaciolacustrine sediments, and recent alluvial deposits of the Connecticut River flood plain. Dry Brook Hill, immediately south of the supply well, is a large subaqueous lacustrine fan and delta formed during the last glacial retreat by sediment deposition into glacial Lake Hitchcock from a meltwater tunnel that was likely near where the Connecticut River cuts through the Holyoke Range. The sediments that compose the aquifer grade from very coarse sand and gravel along the northern flank of the hill, to medium sands in the body of the hill, and to finer-grained sediments along the southern flank of the hill. The interbedded and overlapping fine-grained lacustrine sediments associated with Dry Brook Hill include varved silt and clay deposits. These fine-grained sediments form a confining bed above the coarse-grained aquifer at the supply well and partially extend under the Connecticut River adjacent to the supply well. Ground-water flow in the aquifer supplying water to Dry Brook well was simulated with the U.S. Geological Survey ground-water-flow modeling code MODFLOW. The Dry Brook aquifer model was calibrated to drawdown data collected from 8 observation wells during an aquifer test conducted by pumping the supply well for 10 days at a rate of 122.2 cubic feet per minute (ft3/min; 914 gallons per minute) and to water levels collected from observation wells across the study area. Generally, the largest hydraulic conductivity values used in the model were in the sand and gravel aquifer near the Dry Brook well, which is consistent with the geologic information. Results of aquifer-test simulation indicated that spatially variable aquifer hydraulic properties and boundary conditions affected heads and ground-water flow near the well. A comparison and analysis of water-level fluctuations in study area wells to fluctuations in the Connecticut River indicated a hydraulic connection of the aquifer with the river, which is also consistent with geologic information. Simulated ground-water levels indicated that most ground water in the study area flowed toward and discharged to the Connecticut River and the Dry Brook well. Small amounts of ground water also discharged to smaller streams (Dry Brook and Bachelor Brook) in the study area. Areas contributing water to the well were delineated with the MODPATH particle-tracking routine. Results of the contributing-area analysis indicated that the greatest sources of water to the well were recharge in the Dry Brook Hill area and infiltration of Connecticut River water in an area beyond the extent of the confining bed where the aquifer is in hydraulic connection with the river. The amount of water entering the Dry Brook well from recharge dominated at a lower pumping rate (40.0 ft3/min); about 90 percent of the pumped water originated from recharge and boundary flow, and infiltration from the Connecticut River supplied the remaining 10 percent. At a high pumping rate (122.2 ft3/min), however, about half of the water pumped from the Dry Brook well originated from recharge and boundary flow (49 percent), and half originated from infiltration of water from the Connecticut River (51 percent). Results of a sensitivity analysis of the extent of areas contributing water to the Dry Brook well when pumped at 122.2 ft3/min indicated that the size of these areas did not substantially change when aquifer properties were varied. In contrast, however, the size of these areas changed most when the recharge

Water-Resources Investigations Report↗
USGS2000

Estimation of hydraulic parameters from an unconfined aquifer test conducted in a glacial outwash deposit, Cape Cod, Massachusetts

An aquifer test conducted in a sand and gravel, glacial outwash deposit on Cape Cod, Massachusetts was analyzed by means of a model for flow to a partially penetrating well in a homogeneous, anisotropic unconfined aquifer. The model is designed to account for all significant mechanisms expected to influence drawdown in observation piezometers and in the pumped well. In addition to the usual fluid-flow and storage processes, additional processes include effects of storage in the pumped well, storage in observation piezometers, effects of skin at the pumped-well screen, and effects of drainage from the zone above the water table. The aquifer was pumped at a rate of 320 gallons per minute for 72-hours and drawdown measurements were made in the pumped well and in 20 piezometers located at various distances from the pumped well and depths below the land surface. To facilitate the analysis, an automatic parameter estimation algorithm was used to obtain relevant unconfined aquifer parameters, including the saturated thickness and a set of empirical parameters that relate to gradual drainage from the unsaturated zone. Drainage from the unsaturated zone is treated in this paper as a finite series of exponential terms, each of which contains one empirical parameter that is to be determined. It was necessary to account for effects of gradual drainage from the unsaturated zone to obtain satisfactory agreement between measured and simulated drawdown, particularly in piezometers located near the water table. The commonly used assumption of instantaneous drainage from the unsaturated zone gives rise to large discrepancies between measured and predicted drawdown in the intermediate-time range and can result in inaccurate estimates of aquifer parameters when automatic parameter estimation procedures are used. The values of the estimated hydraulic parameters are consistent with estimates from prior studies and from what is known about the aquifer at the site. Effects of heterogeneity at the site were small as measured drawdowns in all piezometers and wells were very close to the simulated values for a homogeneous porous medium. The estimated values are: specific yield, 0.26; saturated thickness, 170 feet; horizontal hydraulic conductivity, 0.23 feet per minute; vertical hydraulic conductivity, 0.14 feet per minute; and specific storage, 1.3x10-5 per foot. It was found that drawdown in only a few piezometers strategically located at depth near the pumped well yielded parameter estimates close to the estimates obtained for the entire data set analyzed simultaneously. If the influence of gradual drainage from the unsaturated zone is not taken into account, specific yield is significantly underestimated even in these deep-seated piezometers. This helps to explain the low values of specific yield often reported for granular aquifers in the literature. If either the entire data set or only the drawdown in selected deep-seated piezometers was used, it was found unnecessary to conduct the test for the full 72-hours to obtain accurate estimates of the hydraulic parameters. For some piezometer groups, practically identical results would be obtained for an aquifer test conducted for only 8-hours. Drawdowns measured in the pumped well and piezometers at distant locations were diagnostic only of aquifer transmissivity.

Massachusetts↗
USGS1999

Estimated hydraulic properties for the surficial- and bedrock-aquifer system, Meddybemps, Maine

Analytical and numerical-modeling methods were used to estimate hydraulic properties of the aquifer system underlying the Eastern Surplus Company Superfund Site in Meddybemps, Maine. Estimates of hydraulic properties are needed to evaluate pathways for contaminants in ground water and to support evaluation and selection of remediation measures for contaminated ground water at this site. The hydraulic conductivity of surficial materials, determined from specific-capacity tests, ranges from 17 to 78 feet per day for wells completed in coarse-grained glaciomarine sediments, and from about 0.1 to 1.Ofoot per day for wells completed in till. The transmissivity of fractured bedrock determined from specific-capacity tests and aquifer tests in wells completed in less than 200 feet of bedrock ranges from about 0.09 to 130 feet squared per day. Relatively high values of transmissivity at the south end of the study area appear to be associated with a high-angle fracture or fracture zone that hydraulically connects two wells completed in bedrock. Transmissivities at six low-yielding (less than 0.5 gallon per minute) wells, which appear to lie within a poorly transmissive block of the bedrock, are consistently in a range of about 0.09 to 0.5 foot squared per day. The estimates of hydraulic conductivity and transmissivity in the southern half of the study area are supported by results of steady-state calibration of a numerical model and simulation of a 24-hour pumping test at a well completed in bedrock. Hydraulic conductivity values for the surficial aquifer used in the model were 30 feet per day for coarse-grained glaciomarine sediments, 0.001 to 0.01 foot per day for fine-grained glaciomarine sediments, and 0.1 to 0.5 foot per day for till. As part of model calibration, a relatively transmissive zone in the surficial aquifer was extended beyond the hypothesized extent of coarse-grained sediments eastward to the Dennys River. Hydraulic conductivity values used for bedrock in the model ranged from 3x10 -4 to 1.5 feet per day. The highest values were in the fracture zone that hydraulically connects two wells and apparently extends to the Dennys River. The transmissivity of bedrock used in the model ranged from 0.03 to 150 feet squared per day, with the majority of the bedrock transmissivities set at 0.3 foot squared per day. Numerical modeling results indicated that a very low ratio of vertical hydraulic conductivity to thickness (1x10 -9 days -l ) was required to simulate a persistent cone of depression near a residential well that lies in the previously identified poorly transmissive block of bedrock.

Maine↗
USGS1998

Using transport model interpretations of tracer tests to study microbial processes in groundwater

It has long been known that microorganisms affect the geochemistry of groundwater. But despite this recognition, little detailed information is available regarding the rates and the factors controlling microbial processes in groundwater. Part of the reason stems from the relatively inaccessible nature of most groundwater and the difficulties encountered in obtaining representative samples of groundwater and subsurface sediments. At the same time, most groundwater systems are nutrient poor or oligotrophic environments in which the resident microorganisms are severely stressed and often nearly inactive. These populations are functioning so slowly that many types of activity measurements designed to assess microbial processes in more productive environments are ineffective for groundwater. However, because groundwater is by far the largest reservoir of freshwater in the world ( Freeze and Cherry 1979 ), our lack of knowledge about groundwater microorganisms and their processes represents a significant void in the study of microbial ecology and in our ability to predict the outcome when these reserves are compromised by contamination.

Book chapter↗