Search USGS⌕ Search

Geology topics

S. P. Garabedian

Publications and source records attributed to S. P. Garabedian.

14 recordsLinked to original sources

Effects of sea-level rise on ground water flow in a coastal aquifer system

The effects of sea-level rise on the depth to the fresh water/salt water interface were simulated by using a density-dependent, three-dimensional numerical ground water flow model for a simplified hypothetical fresh water lens that is similar to shallow, coastal aquifers found along the Atlantic coast of the United States. Simulations of sea-level rise of 2.65 mm/year from 1929 to 2050 resulted in an increase in water levels relative to a fixed datum, yet a net decrease in water levels relative to the increased sea-level position. The net decrease in water levels was much greater near a gaining stream than farther from the stream. The difference in the change in water levels is attributed to the dampening effect of the stream on water level changes in response to sea-level rise. In response to the decreased water level altitudes relative to local sea level, the depth to the fresh water/salt water interface decreased. This reduction in the thickness of the fresh water lens varied throughout the aquifer and was greatly affected by proximity to a ground water fed stream and whether the stream was tidally influenced. Away from the stream, the thickness of the fresh water lens decreased by about 2% from 1929 to 2050, whereas the fresh water lens thickness decreased by about 22% to 31% for the same period near the stream, depending on whether the stream was tidally influenced. The difference in the change in the fresh water/salt water interface position is controlled by the difference in the net decline in water levels relative to local sea level. ?? 2007 National Ground Water Association.

Ground Water↗

Use of colloid filtration theory in modeling movement of bacteria through a contaminated sandy aquifer

\A filtration model commonly used to describe removal of colloids during packed-bed filtration in water treatment applications was modified for describing downgradient transport of bacteria in sandy, aquifer sediments. The modified model was applied to the results of a small-scale (7 m), natural-gradient tracer test and to observations of an indigenous bacterial population moving downgradient within a plume of organically contaminated groundwater in Cape Cod, MA. The model reasonably accounted for concentration histories of labeled bacteria appearing at samplers downgradient from the injection well in the tracer experiment and for the observed 0.25-μm increase in average cell length for an unlabeled, indigenous bacterial population, 0.6 km downgradient from the source of the plume. Several uncertainties were apparent in applying filtration theory to problems involving transport of bacteria in groundwater. However, adsorption (attachment) appeared to be a major control of the extent of bacterial movement downgradient, which could be described, in part, by filtration theory. Estimates of the collision efficiency factor, which represents the physicochemical factors that determine adsorption of the bacteria onto the grain surfaces, ranged from 5.4 x 10 -3 to 9.7 x 10 -3 .

Environmental Science & Technology↗

Simulation of variable-density flow and transport of reactive and nonreactive solutes during a tracer test at Cape Cod, Massachusetts

A multispecies numerical code was developed to simulate flow and mass transport with kinetic adsorption in variable-density flow systems. The two-dimensional code simulated the transport of bromide (Br − ), a nonreactive tracer, and lithium (Li + ), a reactive tracer, in a large-scale tracer test performed in a sand-and-gravel aquifer at Cape Cod, Massachusetts. A two-fraction kinetic adsorption model was implemented to simulate the interaction of Li + with the aquifer solids. Initial estimates for some of the transport parameters were obtained from a nonlinear least squares curve-fitting procedure, where the breakthrough curves from column experiments were matched with one-dimensional theoretical models. The numerical code successfully simulated the basic characteristics of the two plumes in the tracer test. At early times the centers of mass of Br − and Li + sank because the two plumes were closely coupled to the density-driven velocity field. At later times the rate of downward movement in the Br − plume due to gravity slowed significantly because of dilution by dispersion. The downward movement of the Li + plume was negligible because the two plumes moved in locally different velocity regimes, where Li + transport was retarded relative to Br − . The maximum extent of downward transport of the Li + plume was less than that of the Br − plume. This study also found that at early times the downward movement of a plume created by a three-dimensional source could be much more extensive than the case with a two-dimensional source having the same cross-sectional area. The observed shape of the Br − plume at Cape Cod was simulated by adding two layers with different hydraulic conductivities at shallow depth across the region. The large dispersion and asymmetrical shape of the Li + plume were simulated by including kinetic adsorption-desorption reactions.

Massachusetts↗

Comparison of denitrification activity measurements in groundwater using cores and natural-gradient tracer tests

The transport of many solutes in groundwater is dependent upon the relative rates of physical flow and microbial metabolism. Quantifying rates of microbial processes under subsurface conditions is difficult and is most commonly approximated using laboratory studies with aquifer materials. In this study, we measured in situ rates of denitrification in a nitrate-contaminated aquifer using small-scale, natural-gradient tracer tests and compared the results with rates obtained from laboratory incubations with aquifer core material. Activity was measured using the acetylene block technique. For the tracer tests, co-injection of acetylene and bromide into the aquifer produced a 30 μM increase in nitrous oxide after 10 m of transport (23−30 days). An advection−dispersion transport model was modified to include an acetylene-dependent nitrous oxide production term and used to simulate the tracer breakthrough curves. The model required a 4-day lag period and a relatively low sensitivity to acetylene to match the narrow nitrous oxide breakthrough curves. Estimates of in situ denitrification rates were 0.60 and 1.51 nmol of N 2 O produced cm - 3 aquifer day - 1 for two successive tests. Aquifer core material collected from the tracer test site and incubated as mixed slurries in flasks and as intact cores yielded rates that were 1.2−26 times higher than the tracer test rate estimates. Results with the coring-dependent techniques were variable and subject to the small-scale heterogeneity within the aquifer, while the tracer tests integrated the heterogeneity along a flow path, giving a rate estimate that is more applicable to transport at the scale of the aquifer.

Environmental Science & Technology↗

Hydrology and digital simulation of the regional aquifer system, eastern Snake River Plain, Idaho

The occurrence and movement of water in the regional aquifer system that underlies the eastern Snake River Plain, Idaho, de- pend on the transmissivity and storage capacity of rocks that compose the geologic framework and on the distribution and amount of recharge and discharge of water within that frame- work. On a regional scale, most water moves horizontally through interflow zones in Quaternary basalt of the Snake River Group. In recharge and discharge areas, water also moves vertically along joints and interfingering edges of basalt flows. Aquifer thickness is largely unknown, but geophysical studies suggest that locally the Quaternary basalt may exceed several thousand feet. Along the margins of the plain, sand and gravel several hundred feet thick transmit large volumes of water. Regional ground-water movement is generally from northeast to southwest, from areas of recharge to areas of discharge. Re- charge is from seepage of surface water used for irrigation, stream and canal losses, underflow from tributary drainage ba- sins, and infiltration of precipitation. Aquifer discharge is largely spring flow to the Snake River and water pumped for irrigation. Major springs are near American Falls Reservoir and along the Snake River from Milner Dam to King Hill. Regional ground-water flow was simulated with numerical models. Initially, a two-dimensional steady-state model that in- cluded a nonlinear, least-squares regression technique was used to estimate aquifer properties. Later, a three-dimensional steady- state and transient model was used to replace the two-dimen- sional model. Three-dimensional model results indicated that average total transmissivity ranged from about 0.05 to 120 feet squared per second and vertical leakance ranged from about 3 x 10 -10 to 5 x 10 -6 feet per second per foot of aquifer thickness. The three-dimensional transient model was used to compare measured and estimated long-term changes in ground-water dis- charge and water levels with simulated values. Initial head con- ditions used in transient simulations were derived from a steady-state solution of estimated preirrigation hydrologic condi- tions. Transient simulations were 5-year stress periods beginning in 1891 and ending in 1980. Recharge for each stress period from 1926 to 1980 was estimated from surface-water irrigation, pre- cipitation, and streamflow records. Recharge for stress periods from 1891 to 1925 was based on the average value for stress peri- ods from 1926 to 1980 and was indexed to estimated irrigated acreages. Average annual tributary drainage-basin underflow for stress periods from 1891 to 1910 was calculated by using basin- yield equations. Underflow for stress periods from 1911 to 1980 was varied by use of streamflow records. Transient simulations reasonably approximated measured changes in aquifer head and ground-water discharge that re- sulted from use of surface water for irrigation. Irrigation with surface water peaked in about 1950; subsequent increases in irri- gation have been supplied largely by ground water. The three- dimensional model simulated water-level declines and reduced ground-water discharge caused in part by increases in ground- water pumping. The transient model was used to simulate aquifer changes from 1981 to 2010 in response to three hypothetical development alternatives: (1) Continuation of 1980 hydrologic conditions, (2) increased pumpage, and (3) increased recharge. Simulation of continued 1980 hydrologic conditions for 30 years indicated that head declines of 2 to 8 feet might be expected in the central part of the plain. The magnitude of simulated head declines was con- sistent with head declines measured during the 1980 water year. Larger declines were calculated along model boundaries, but these changes may have resulted from underestimation of tribu- tary drainage-basin underflow and inadequate aquifer definition. Simulation of increased ground-water pumpage (an additional 2,400 cubic feet per second) for 30 years indicated head declines of 10 to 50 feet in the central part of the plain. These relatively large head declines were accompanied by increased simulated river leakage of 50 percent and decreased spring discharge of 20 percent. The effect of increased recharge (800 cubic feet per sec- ond) for 30 years was a rise in simulated heads of 0 to 5 feet in the central part of the plain.

Idaho;Oregon↗

In situ measurement of methane oxidation in groundwater by using natural-gradient tracer tests

Methane oxidation was measured in an unconfined sand and gravel aquifer (Cape Cod, Mass.) by using in situ natural-gradient tracer tests at both a pristine, oxygenated site and an anoxic, sewage-contaminated site. The tracer sites were equipped with multilevel sampling devices to create target grids of sampling points; the injectate was prepared with groundwater from the tracer site to maintain the same geochemical conditions. Methane oxidation was calculated from breakthrough curves of methane relative to halide and inert gas (hexafluoroethane) tracers and was confirmed by the appearance of 13 C-enriched carbon dioxide in experiments in which 13 C-enriched methane was used as the tracer. A V max for methane oxidation could be calculated when the methane concentration was sufficiently high to result in zero-order kinetics throughout the entire transport interval. Methane breakthrough curves could be simulated by modifying a one-dimensional advection-dispersion transport model to include a Michaelis-Menten-based consumption term for methane oxidation. The K m values for methane oxidation that gave the best match for the breakthrough curve peaks were 6.0 and 9.0 μ M for the uncontaminated and contaminated sites, respectively. Natural-gradient tracer tests are a promising approach for assessing microbial processes and for testing in situ bioremediation potential in groundwater systems.

Massachusetts↗

Application of a parameter-estimation technique to modeling the regional aquifer underlying the eastern Snake River Plain, Idaho

A nonlinear, least-squares regression technique for the estimation of ground-water flow model parameters was applied to the regional aquifer underlying the eastern Snake River Plain, Idaho. The computer program simulates two-dimensional, steady-state ground-water flow. Hydrologic data for the 1980 water year were used to calculate recharge rates, boundary fluxes, and spring discharges. Groundwater use was estimated from irrigated land maps and crop consumptive-use figures. These estimates of ground-water withdrawal, recharge rates, and boundary flux, along with leakance, were used as known values in the model calibration of transmissivity. Leakance values were adjusted between regression solutions by comparing model-calculated to measured spring discharges. In other simulations, recharge and leakance also were calibrated as prior information regression parameters, which limits the variation of these parameters using a normalized standard error of estimate. Results from a best-fit model indicate a wide areal range in transmissivity from about 0.05 to 44 feet squared per second and in leakance from about 2.2 x 10 -9 to 6.0 x 10 -8 feet per second per foot. Along with parameter values, model statistics also were calculated, including the coefficient of correlation between computed and observed head (0.996), standard error of the estimates of head (40 feet), and parameter coefficients of variation (about 10-40 per-cent). Additional boundary flux was added in some areas during calibration to achieve proper fit to ground-water flow directions. Model fit improved significantly when areas that violated model assumptions were removed. Model fit also improved when y-direction (northwest-southeast) transmissivity values were larger than x-direction (northeast-southwest) transmissivity values. The model was most sensitive to changes in recharge, and in some areas, to changes in transmissivity, particularly near the spring discharge area from Milner to King Hill.

Idaho↗

Front-tracking model for convective transport in flowing ground water

A finite-difference numerical model that simulates the convective transport of water or tracer particles through porous media is described. It can be applied to one- or two-dimensional problems involving either steady-state or transient flow. The model tracks representative water or tracer particles, initially located along specified lines, as they move in response to the groundwater velocity filed. Aquifer properties may be both anisotropic and nonhomogeneous. Included in the report is a listing of the program along with input formats and test problem results. The front-tracking model provides a useful first approximation determining the movement of solutes in an aquifer, particularly in cases when dispersion and dilution is of minor consideration. (USGS)

Water-Resources Investigations Report↗