GEOSTATISTICAL ESTIMATION OF HORIZONTAL HYDRAULIC CONDUCTIVITY FOR THE KIRKWOOD‐COHANSEY AQUIFER1

ABSTRACT: The Kirkwood‐Cohansey aquifer has been identified as a critical source for meeting existing and expected water supply needs for southern New Jersey. Several contaminated sites exist in the region; their impact on the aquifer has to be evaluated using ground water flow and transport models. Ground water modeling depends on availability of measured hydrogeologic data (e.g., hydraulic conductivity, for parameterization of the modeling runs). However, field measurements of such critical data have inadequate spatial density, and their locations are often clustered. The goal of this study was to research, compile, and geocode existing data, then use geostatistics and advanced mapping methods to develop a map of horizontal hydraulic conductivity for the Kirkwood‐Cohansey aquifer. Spatial interpolation of horizontal hydraulic conductivity measurements was performed using the Bayesian Maximum Entropy (BME) Method implemented in the BMELib code library. This involved the integration of actual measurements with soft information on likely ranges of hydraulic conductivity at a given location to obtain estimate maps. The estimation error variance maps provide an insight into the uncertainty associated with the estimates, and indicate areas where more information on hydraulic conductivity is required.     

ESTIMATING TRANSMISSIVITY AND HYDRAULIC CONDUCTIVITY OF CHICOT AQUIFER FROM SPECIFIC CAPACITY DATA1

ABSTRACT: Specific capacity data obtained from Well Construction reports which are available from USGS offices, can provide useful estimates of tranamissivity (T), and hydraulic conductivity (K), of an aquifer. The Chicot Aquifer in Louisiana is one of the largest sources of fresh ground water in North America. Hydrologic data collected for the Chicot Aquifer indicate that specific capacity tests can be used in estimating local and regional values for T and K, if the Cooper-Jacob equation for transient flow is used with proper corrections for well loss and partial penetration. Where full scale pumping test data are scarce, specific capacity test data that are adequately distributed spatially can be used to map changes in T and K values and can be summarized statistically to indicate applicable regional values. A computer program called “TGUESS” which is available from International Ground Water Modeling Center, Holcomb Research Institute, was used in this study. The contour maps for T and K values are prepared for different well depth intervals to avoid wide variation of values.     

ASSESSMENT OF HYDRAULIC PROPERTIES IN AN UNCONFINED ALLUVIAL AQUIFER NEAR GRAND ISLAND,NEBRASKA1

ABSTRACT: A method is presented to analyze time-drawdown data from one or more observation wells for the calculation of four hydraulic parameters for unconfined aquifers: vertical hydraulic conductivity, horizontal hydraulic conductivity, storage coefficient, and specific yield. The hydraulic parameter results are further analyzed for reliability and the possible ranges of the actual parameter values. After verification using a theoretical example, the method was used to analyze pumping test data from 22 observation wells in an unconfined alluvial aquifer near Grand Island, Nebraska. Results indicate that this method can be used to efficiently calculate the four hydraulic parameters in this type of aquifer. The method can also identify the impact of measurement errors on the parameter estimates, and provide ranges of the actual parameter values. The parameter values calculated using this method were compared to those determined using other theories. It is found that this method is very useful for calculating the hydraulic properties from pumping test data and for analyzing the parameter reliability.     

PRESENCE AND DISTRIBUTION OF CHLORINATED ORGANIC COMPOUNDS IN STREAMBED SEDIMENTS,NEW JERSEY1

ABSTRACT: Concentrations of 18 hydrophobic chlorinated organic compounds in streambed sediments from 100 sites throughout New Jersey were examined to determine (1) which compounds were detected most frequently, (2) whether detection frequencies differed among selected drainage basins, and (3) whether concentrations differed significantly among selected drainage basins. Twelve drainage basins across New Jersey that contain a range of land-use patterns and population densities were selected to represent various types and degrees of development. To ensure an adequate number of samples for statistical comparison among drainage basins, the 12 selected basins were consolidated into seven drainage areas on the basis of similarities in land-use patterns and population densities. Additionally, data for three classes of chlorinated organic compounds in streambed sediments from 255 sites throughout New Jersey were examined to determine whether the presence of these compounds in streambed sediments is related to the type and degree of development within the drainage area of each sampling site. Chlorinated organic compounds detected most frequently within the seven representative drainage areas were DDT, DDE, DDD, chlordane, dieldrin, and PCBs. DDT, DDE, and DDD, which were the most widely distributed organic compounds, were detected in about 60 to 100 percent of the samples from all drainage areas but one (where the detection rate for these compounds was about 20 to 40 percent). Chiordane and dieldrin were detected in about 80 to 100 percent of samples from highly urbanized and populated drainage areas; detection frequencies for these compounds tended to be smaller in less developed and populated areas. PCBs were detected in about 40 to 85 percent of samples from all drainage areas; detection frequencies were highest in the most heavily developed and populated areas. Analysis of variance on rank-transformed organic compound concentrations normalized to sediment organic carbon content was used to evaluate differences in concentrations among the seven representative drainage areas. Chlordane and PCBs were the chlorinated organic compounds with the most highly elevated concentrations in streambed sediments across the State. Median normalized concentrations of all six of the most frequently detected chlorinated organic compounds were highest in the most heavily urbanized and populated drainage area and lowest in the less populated, predominantly agricultural or forested areas. Concentrations of DDT and DDE, however, did not differ significantly among most of the drainage areas. Concentrations of DDD, chlordane, dieldrin, and PCBs differed significantly among drainage areas. The highest median normalized concentrations were found in samples from the most heavily urbanized and populated areas, and the lowest were in samples from the least developed, most heavily forested area. Logistic regression was used to examine relations between the presence of hydrophobic chlorinated organic compounds in streambed sediments at specified concentrations and variables that characterize the type and degree of development within the drainage areas of 255 sites across New Jersey. The explanatory variables found most useful for predicting the presence of chlorinated organic compounds in streambed sediments include total population and amounts (in square kilometers) of various land-use categories. Logistic regression equations were developed to identify significant relations between population and amounts of specific land-use categories within drainage areas and the probability of detecting chlorinated organic contaminants in streambed sediments. These relations can be used to assist in the identification of geographic regions of primary concern for contamination of bed sediments by chlorinated organic compounds across the State.     

EVALUATION OF A STREAM AQUIFER ANALYSIS TEST USING ANALYTICAL SOLUTIONS AND FIELD DATA1

ABSTRACT: Considerable advancements have been made in the development of analytical solutions for predicting the effects of pumping wells on adjacent streams and rivers. However, these solutions have not been sufficiently evaluated against field data. The objective of this research is to evaluate the predictive performance of recently proposed analytical solutions for unsteady stream depletion using field data collected during a stream/aquifer analysis test at the Tamarack State Wildlife Area in eastern Colorado. Two primary stream/aquifer interactions exist at the Tamarack site: (1) between the South Platte River and the alluvial aquifer and (2) between a backwater stream and the alluvial aquifer. A pumping test is performed next to the backwater stream channel. Drawdown measured in observation wells is matched to predictions by recently proposed analytical solutions to derive estimates of aquifer and streambed parameters. These estimates are compared to documented aquifer properties and field measured streambed conductivity. The analytical solutions are capable of estimating reasonable values of both aquifer and streambed parameters with one solution capable of simultaneously estimating delayed aquifer yield and stream flow recharge. However, for long term water management, it is reasonable to use simplified analytical solutions not concerned with early‐time delayed yield effects. For this site, changes in the water level in the stream during the test and a varying water level profile at the beginning of the pumping test influence the application of the analytical solutions.     

SOLUTE TRANSPORT IN THREE‐DIMENSIONAL SATURATED POROUS MEDIA HAVING SELF‐SIMILAR HYDRAULIC CONDUCTIVITY DISTRIBUTION1

ABSTRACT: Customarily, it has been assumed that hydraulic conductivity is a stationary, homogeneous stochastic process with a finite variance for stochastic analysis of solute transport in the subsurface. That the distribution of hydraulic conductivity may have a fractal behavior with long range correlations was suggested from field data analyses. This motivates us to further investigate how the fractal behavior of permeability distribution impacts solute transport in porous media. This study provides longitudinal and transverse macrodispersivity coefficients and the variance of the solute concentration. Longitudinal and transverse macrodispersivity coefficients are found to depend strongly on the fractal dimension (D) of logarithmic hydraulic conductivity (logK). The longitudinal and transverse macrodispersivity coefficients are the highest when D = 1, and the values decrease monotonically to zero at D = 2. Both coefficients correspond to the characteristic length scale of the logK distribution, thus are scale dependent parameters. The ratio of the transverse to the longitudinal macrodispersivity coefficient is on the order of 10^‐1 to 10^‐4. Concentration variance also decreases with the fractal dimension of logK. There is no spatial spreading of solute for D = 2, and the concentration variance reaches zero for this case.     

VELOCITY DISTRIBUTIONS IN STREAMBED SIMULATION CULVERTS USED FOR FISH PASSAGE1

ABSTRACT: Stream crossings designed to simulate a natural streambed inside a culvert are commonly used to meet criteria where fish passage is required. Currently, there has been little research that quantitatively examines the flow patterns within these streambed simulation culverts. A model is presented that estimates the percent of a cross section that is within the swimming ability of juvenile fish developed from velocity measurements on first‐generation stream simulation culverts — those culverts with a continuous gravel bed. The model, developed with regression techniques, uses physical and hydraulic parameters including discharge, total cross sectional area, Froude number, and relative roughness, and was tested directly against velocity distributions computed from field measurements and in a culvert design mode. Results were favorable, although larger percentage errors exist, particularly at small flow depths. The model appears to underestimate the percent of channel cross section at or below a limiting velocity, hence it is generally conservative in design mode.     

IMPACT OF IRRIGATION WELLS ON BASEFLOW OF THE BIG BLUE RWER,NEBRASKA1

ABSTRACT: The Kansas-Nebraska Big Blue River compact requires that the state of Nebraska insure a minimum flow of the Big Blue River across the state line. There are two options that the state of Nebraska may use to ensure minimum flows. The obvious option is to limit surface-water irrigators along the river. However, under the terms of the compact, a second option may be to regulate irrigation wells that are within one mile of the river and were installed after November 1, 1968. The objective of this study is to quantify the effects of 17 irrigation wells that may be regulated on baseflow of the Big Blue River. A finite-element model is used to study the hydrogeologic system between DeWitt and Beatrice, Nebraska. The 17 wells that may be regulated are located between these towns and are developed in sediments deposited in a cross-cutting paleovalley anchor alluvium associated with the Big Blue River. While there wore considerable existing data, additional data were gathered by drilling an additional nine test holes, conducting several aquifer tests, stream-stage measurements, and baseflow calculation through extensive stream-discharge measurements, establishment of a ground water-level monitoring network, determining the amount of water pumped for irrigation and municipal use in the area, and a short-term precipitation network. The model was calibrated using observed baseflow and ground water level data. The model clearly shows that regulating the 17 wells to maintain baseflow would have a minimal effect on the overall water budget. This is reasonable, especially considering that there are over 250 irrigation wells in the project area. The 17 wells considered pumped only 6 percent of the total pumpage within the modeled area during the irrigation season of 1984. The computer model provides the documentation needed to demonstrate this fact. Although much of the resources spent and a significant amount of hydrogeologic data are being collected over a period of three years on a relatively small area, the simulation model could be improved through further field testing of the aquifer and stream-bed sediment characteristics and quantification of ground water recharge, discharge, and evapotranspiration rates.     

DETERMINATION OF UNCONFINED AQUIFER HYDRAULIC PROPERTIES FROM RECOVERY TEST DATA1

ABSTRACT: This paper analyzes the sensitivity of drawdown to four hydraulic parameters in unconfined aquifers: horizontal and vertical hydraulic conductivity K_r and K_z, storage coefficient S, and specific yield S_y. Sensitivity coefficients indicate that the sensitivity vanes with time for each aquifer parameter, and K_r, K_z, S, and S_y are identifiable from recovery test data. An inverse method was used to calculate the four parameters from residual drawdowns. Results of application examples demonstrate that residual data provide valid information in the determination of unconfined aquifer hydraulic parameters.     

Analyzing the Effects of Groundwater Pumping on an Urban Stream‐Aquifer System

When the cone of influence of a pumping well reaches a nearby river, the resulting hydraulic gradient can induce enhanced seepage of streamflow into the aquifer. The rate of seepage is often modeled using analytical solutions that are simple to apply but may not reproduce field data due to mathematical assumptions not being met in the field. Furthermore, the appropriateness of such models has not been investigated in detail due to difficulty in measuring streamflow loss in the field. In this study, a field experiment was conducted on a reach of the South Platte River near Denver, Colorado to estimate pumping‐induced streamflow loss. A network of stream gauges, monitoring wells, and in situ measurements was used to observe streamflow rates, groundwater levels, and temperature to assess if pumping wells have a significant impact on streamflow, and to compare observed streamflow depletion against analytical solutions. Data collected suggest that pumping wells have a noticeable impact on streamflow. The analytical solutions proved accurate if streamflow was low and constant but performed poorly if streamflow was high and variable. Therefore, for this reach, the use of analytical solutions to predict streamflow may only be appropriate under low‐flow, constant‐flow conditions. Methods and results can be used to guide other streamflow depletion studies and to inform cases of pumping‐induced streamflow depletion, particularly in regard to water rights.     

ASSESSING CHANGES IN WATERSHED FLOW REGIMES WITH SPATIALLY EXPLICIT HYDRAULIC MODELS1

ABSTRACT: Urbanization, farming, and other watershed activities can significantly alter storm hydrographs and sediment erosion rates within a watershed. These changes routinely cause severe economic and ecological problems manifested in the form of increased flooding and significant changes in channel morphology. As the activities within a watershed influence the hydrologic, hydraulic, and ecological conditions within a river, interdisciplinary approaches to predict and assess the impacts that different land uses have on streams need to be developed. An important component of this process is ascertaining how hydrologic changes induced by specific watershed activities will affect hydraulic conditions and the accompanying flood levels, sediment transport rates, and habitat conditions within a stream. A conceptual model for using spatially explicit (two‐dimensional) hydraulic models to help evaluate the impacts that changes in flow regime might have on a river is presented. This framework proposes that reproducing and quantifying flow complexity allows one to compare the hydraulic conditions within urban, urbanizing, and non‐urban streams in a more biologically and economically meaningful way. The justification, advantage, and need for such a method is argued through the results of one‐ and two‐dimensional hydraulic model studies. The implementation of this methodology in watershed urbanization studies is described.     

THE DEPENDENCE OF THE RESIDUAL GRAVITY ON HYDRAULIC CONSTANTS IN GLACIAL DEPOSITS1

ABSTRACT: Buried glacial stream channels contain large and easily accessible groundwater resources. Gravity surveys have been frequently applied for their location. A gravity survey in the geohydrologically explored Wood River Valley Area of southern Rhode Island shows extreme lows of -2 mgals over channel depths of maximal 300 feet. Three gravity profiles were observed in east-west direction across a north-south striking stream channel. The bedrock depth increases rapidly towards the south from 130 to 300 feet. The gravity lows observed across each profile are not related to the bedrock depth but rather to the saturated thickness of the main quifer and its hydraulic transmissivity. Well logs indicate that the large change of bedrock depth is solely due to an increase of till of low permeability. The volume of the glacial outwash, which is the major groundwater resource, changes little underneath the three profiles. The gravity lows apear to be directly related to the density contrast between glacial outwash and till. The response to the hydraulically more pertinent units renews the interest in the gravity method as it may have a potential to estimate yields of hydrologically complex aquifers     

SATURATED SOIL HYDRAULIC CONDUCTWITY IN LOBLOLLY PINE PLANTATIONS ON DRAINED SITES1

ABSTRACT: Variation of in situ measured saturated hydraulic conductivity (KS) with stand age was examined in drained and intensively managed loblolly pine (Pinus taeda L.) plantations on very poorly drained Bayboro loam soils. Stand ages studied were 1-year-old and 14-years old. No differences in Ks values were found between the stand ages. In addition, no differences in core measured soil properties were found between the stand ages, indicating that there were no differences in the pore structure of the soil matrix. There was large variation of Ks within stands and between stands within ages. The mean within stand Ks values ranged from 0.66 cm/hr to 4.85 cm/hr. The frequency of tests exhibiting pipe flow through large non-capillary voids was significantly greater in the older stands; however, the continuity of the voids in the soil, and whether or not non-Darcy type flow would occur in a saturated profile, could not be determined.     

TRANSIENT RESPONSE FUNCTIONS FOR CONJUNCTIVE WATER MANAGEMENT IN THE SNAKE RIVER PLAIN,IDAHO1

ABSTRACT: Increasing demands on western water are causing a mounting need for the conjunctive management of surface water and ground water resources. Under western water law, the senior water rights holder has priority over the junior water rights holder in times of water shortage. Water managers have been reluctant to conjunctively manage surface water and ground water resources because of the difficulty of quantification of the impacts to surface water resources from ground water stresses. Impacts from ground water use can take years to propagate through an aquifer system. Prediction of the degree of impact to surface water resources over time and the spatial distribution of impacts is very difficult. Response functions mathematically describe the relationship between a unit ground water stress applied at a specific location and stream depletion or aquifer water level change elsewhere in the system. Response functions can be used to help quantify the spatial and temporal impacts to surface water resources caused by ground water pumping. This paper describes the theory of response functions and presents an application of transient response functions in the Snake River Plain, Idaho. Transient response functions can be used to facilitate the conjunctive management of surface and ground water not only in the eastern Snake River Plain basin, but also in similar basins throughout the western United States.     

Estimating Hydraulic Properties of Volcanic Aquifers Using Constant‐Rate and Variable‐Rate Aquifer Tests1

Abstract: In recent years the ground‐water demand of the population of the island of Maui, Hawaii, has significantly increased. To ensure prudent management of the ground‐water resources, an improved understanding of ground‐water flow systems is needed. At present, large‐scale estimations of aquifer properties are lacking for Maui. Seven analytical methods using constant‐rate and variable‐rate withdrawals for single wells provide an estimate of hydraulic conductivity and transmissivity for 103 wells in central Maui. Methods based on constant‐rate tests, although not widely used on Maui, offer reasonable estimates. Step‐drawdown tests, which are more abundantly used than other tests, provide similar estimates as constant‐rate tests. A numerical model validates the suitability of analytical solutions for step‐drawdown tests and additionally provides an estimate of storage parameters. The results show that hydraulic conductivity is log‐normally distributed and that for dike‐free volcanic rocks it ranges over several orders of magnitude from 1 to 2,500 m/d. The arithmetic mean, geometric mean, and median values of hydraulic conductivity are respectively 520, 280, and 370 m/d for basalt and 80, 50, and 30 m/d for sediment. A geostatistical approach using ordinary kriging yields a prediction of hydraulic conductivity on a larger scale. Overall, the results are in agreement with values published for other Hawaiian islands.     

ESTIMATING STREAMFLOW AND ASSOCIATED HYDRAULIC GEOMETRY,THE MID‐ATLANTIC REGION,USA1

ABSTRACT: Methods to estimate streamflow and channel hydraulic geometry were developed for unpaged streams in the Mid‐Atlantic Region. Observed mean annual streamflow and associated hydraulic geometry data from 75 gaging stations in the Appalachian Plateau, the Ridge and Valley, and the Piedmont Physiographic Provinces of the Mid‐Atlantic Region were used to develop a set of power functions that relate streamflow to drainage area and hydraulic geometry to streamflow. For all three physiographic provinces, drainage area explained 95 to 98 percent of the variance in mean annual streamflow. Relationships between mean annual streamflow and water surface width and mean flow depth had coefficients of determination that ranged from R²= 0.55 to R²= 0.91, but the coefficient of determination between mean flow velocity and mean annual streamflow was lower (R²= 0.44 to R²= 0.54). The advantages of using the regional regression models to estimate streamflow over a conceptual model or a water balance model are its ease of application and reduced input data needs. The prediction of the regression equations were tested with data collected as part of the U.S. Environmental Protection Agency (USEPA) Environmental Monitoring and Assessment Program (EMAP). In addition, equations to transfer streamflow from gaged to ungaged streams are presented.     

Evaluation of Ground‐Water and Surface‐Water Exchanges Using Streamflow Difference Analyses1

Abstract: In the karstic lower Flint River Basin, limestone fracturing, jointing, and subsequent dissolution have resulted in the development of extensive secondary permeability and created a system of major conduits that facilitate the exchange of water between the Upper Floridan aquifer and Flint River. Historical streamflow data from U.S. Geological Survey gaging stations located in Albany and Newton, Georgia, were used to quantify ground‐water and surface‐water exchanges within a 55.3 km section of the Flint River. Using data from 2001, we compared estimates of ground‐water flux using a time adjustment method to a water balance equation and found that these independent approaches yielded similar results. The associated error was relatively large during high streamflow when unsteady conditions prevail, but much lower during droughts. Flow reversals were identified by negative streamflow differences and verified with in situ data from temperature sensors placed inside large spring conduits. Long‐term (13 years) analysis showed negative streamflow differentials (i.e., a losing stream condition) coincided with high river stages and indicated that streamflow intrusion into the aquifer could potentially exceed 150 m³/s. Although frequent negative flow differentials were evident, the Flint River was typically a gaining stream and showed a large net increase in flow between the two gages when examined over the period 1989‐2003. Ground‐water contributions to this stream section averaged 2‐42 m³/s with a mean of 13 m³/s. The highest rate of ground‐water discharge to the Flint River occurred during the spring when regional ground‐water levels peaked following heavy winter and spring rains and corresponding rates of evapotranspiration were low. During periods of extreme drought, ground‐water contributions to the Flint River declined.     

Evaluation of Rio Grande Management Alternatives Using a Surface‐Water/Ground‐Water Model1

Abstract: Previous investigations observed significant seepage losses from the Rio Grande to the shallow aquifer between Socorro and San Antonio, New Mexico. High‐resolution telescopic modeling was used along a 10‐km reach of the Rio Grande and associated drains and canals to evaluate several management alternatives aimed at improving river conveyance efficiency. Observed data consisted of ground‐water and surface‐water elevations, seepage rates along the Rio Grande and associated canals and drains, and borehole geology. Model calibration was achieved by adjusting hydraulic conductivity and specific storage until the output matched observed data. Sensitivity analyses indicated that the system was responsive to changes in hydrogeologic properties, especially when such alterations increased vertical connectivity between layers. The calibrated model predicted that removal of the low flow conveyance channel, a major channel draining the valley, would not only decrease river seepage by 67%, but also decrease total flow through the reach by 75%. The decreased flow through the reach would result in increased water logging and an average increase in ground‐water elevations of 1.21 meter. Simulations of the system with reduced riparian evapotranspiration rates or a relocated river channel also predicted decreased river seepage, but to a much lesser degree.     

OBSERVATION WELL NETWORK DESIGN FOR PUMPING TESTS IN UNCONFINED AQUIFERS1

ABSTRACT: This paper presents a method for determining locations of observation wells to be used in conducting pumping tests in unconfined aquifers. Sensitivity coefficients, the distribution of relative errors, and the correlation coefficients between four aquifer parameters (horizontal and vertical hydraulic conductivities K_r and K_z, storage coefficient S, and specific yield S_y) are used as the criteria for the design of observation well networks and the interpretation of pumping tests. The contours of the relative errors over a vertical profile are very useful in selecting the “best” location of an observation well. Results from theoretical analyses suggest that a wide range of locations is suitable for the determination of K_r and that good locations for the determination of K_z and S may be poorly suited for the determination of S_y. Consideration must be given to the position and lengths of the pumping well screen in the selection of observation well locations. For a given location, the quality of test data can be improved by using high pumping rates and frequent sampling of drawdowns. We found that a minimum of two and preferably three observation locations are needed along a given transect. Results of the four parameters from a single well analysis may contain higher uncertainties. However, composite analyses of multiple observation wells can reduce the correlation between the four aquifer parameters, particularly between K_r and S_y, thus, improving the quality of parameter estimation. Results from two pumping tests conducted at sites located in Nebraska were examined with regard to the proposed methodology.     

REMOTE MONITORING OF SELECTED GROUND‐WATER DOMINATED LAKES IN THE NEBRASKA SAND HILLS1

ABSTRACT: The Landsat‐Muitispectral Scanner (MSS) data were used to measure lake area fluctuations (1972–1989) for 130 ground‐water dominated lakes in the Western Lakes Region of the Nebraska Sand Hills. In general, the pattern shown in lake area hydrographs was similar to that for in‐situ lake elevations. In‐situ lake‐elevation data verify that remote monitoring of surface‐area fluctuations, even at relatively coarse spatial resolution, is not only practical and useful, but also it elucidates the hydrologic characteristics of groundwater‐dominated lakes of the Sand Hills. The apparent differences in behavior between lakes in the northern and southern portions of the study area may be related to both their location in the regional ground water system and the substantial local hydrologic complexity.