Field Study of Hydraulic Conductivity in a Heterogeneous Aquifer: Comparison of Single-Borehole Measurements Using Different Instruments

This field study compares three techniques for estimating the vertical distribution of horizontal hydraulic conductivity Kr in a heterogeneous aquifer and evaluates possible support volume effects. The dipole flow test (DFT), multilevel slug test (MLST), and borehole flowmeter test (BFT) are based on different kinematic flow structures and the shape and the size of the support volumes. The experiment design employed an identical characteristic linear scale for all tests. Vertical profiles of Kr ranging up to 260 m/day from tested wells in an alluvial aquifer exhibit a strong correlation in spite of the differences between test hydraulics. Results suggest that tested screen length is an important indicator of the averaging mechanism for hydraulic tests. Correlation between the DFT and MLST is especially strong. Correlation between data from the BFT and other tests is not as strong due to the absence of a distinct physical vertical scale, among other factors. The differences between the tests are discussed using the concept of a weighting function associated with the magnitude of instantaneous local velocity

Reply to comment by H. Lough, Department of Civil Engineering, University of Canterbury, Christchurch, New Zealand, on the paper “Stream depletion predictions using pumping test data from a heterogeneous stream–aquifer system (a case study from the Great Plains, USA)”

1. General remark 2. The study by Kollet and Zlotnik (2003) 3. Remark on the explanation of the drawdown behavior 4. Remark on the re-analysis of the data from piezometer C2d 5. Summar

An approach to hydrogeological modeling of a large system of groundwater-fed lakes and wetlands in the Nebraska Sand Hills, USA

The feasibility of a hydrogeological modeling approach to simulate several thousand shallow groundwater-fed lakes and wetlands without explicitly considering their connection with groundwater is investigated at the regional scale (~40,000 km2) through an application in the semi-arid Nebraska Sand Hills (NSH), USA. Hydraulic heads are compared to local land-surface elevations from a digital elevation model (DEM) within a geographic information system to assess locations of lakes and wetlands. The water bodies are inferred where hydraulic heads exceed, or are above a certain depth below, the land surface. Numbers of lakes and/or wetlands are determined via image cluster analysis applied to the same 30-m grid as the DEM after interpolating both simulated and estimated heads. The regional water-table map was used for groundwater model calibration, considering MODIS-based net groundwater recharge data. Resulting values of simulated total baseflow to interior streams are within 1% of observed values. Locations, areas, and numbers of simulated lakes and wetlands are compared with Landsat 2005 survey data and with areas of lakes from a 1979–1980 Landsat survey and the National Hydrography Dataset. This simplified process-based modeling approach avoids the need for field-based morphology or water-budget data from individual lakes or wetlands, or determination of lake-groundwater exchanges, yet it reproduces observed lake-wetland characteristics at regional groundwater management scales. A better understanding of the NSH hydrogeology is attained, and the approach shows promise for use in simulations of groundwater-fed lake and wetland characteristics in other large groundwater systems. Includes supplementary materials

Analysis of Groundwater Recharge in Mongolian Drylands Using Composite Vadose Zone Modeling

Knowledge of groundwater recharge (GR) is important for the effective management of water resources under semi-arid continental climates. Unfortunately, studies and data in Mongolia are limited due to the constraints in funding and lack of research infrastructures. Currently, the wide accessibility of freely available global-scale digital datasets of physical and chemical soil properties, weather data, vegetation characteristics, and depths to the water table offers new tools and basic information that can support low-cost physically based and process-oriented models. Estimates of GR over 41 study sites in Mongolia were obtained using HYDRUS-1D in a 2-m-thick soil profile with root depths of either 0.30 or 0.97m by exploiting the daily precipitation and biome-specific potential evapotranspiration values. The GR simulated by HYDRUS-1D arrives at the water table and becomes the actual GR with a lag time that has been calculated using a simplified form of the Richards equation and a traveling wave model. The mean annual precipitation ranges from 57 to 316mm year−1, and on average about 95% of it is lost by mean annual actual evapotranspiration. In the steppe region, the vegetation cover induces higher-than-normal actual transpiration losses and consequently lower GR. The mean annual GR rates span between 0.3 and 12.0mm year−1, while travel times range between 4 and 558 years. Model prediction uncertainty was quantified by comparing actual evapotranspiration and GR with available maps and by a sensitivity assessment of lag time to the soil moisture in the deep vadose zone. The partial least squares regression (PLSR) was used to evaluate the impact of available environmental properties in explaining the 47.1 and 59.1%variability of the spatially averaged mean annual GR and travel time, respectively. The most relevant contributors are clay content, aridity index, and leaf area index for GR, and depth to the water table and silt content for the lag time. In data-poor, arid, and semi-arid regions such as Mongolia, where the mean annual GR rates are low and poorly correlated to precipitation, the ever-increasing availability of world databases and remote sensing products offers promise in estimating GR

Recognition of Regional Water Table Patterns for Estimating Recharge Rates in Shallow Aquifers

We propose a new method for groundwater recharge rate estimation in regions with stream-aquifer interactions, at a linear scale on the order of 10 km and more. The method is based on visual identification and quantification of classically recognized water table contour patterns. Simple quantitative analysis of these patterns can be done manually from measurements on a map, or from more complex GIS data extraction and curve fitting. Recharge rate is then estimated from the groundwater table contour parameters, streambed gradients, and aquifer transmissivity using an analytical model for groundwater flow between parallel perennial streams. Recharge estimates were obtained in three regions (areas of 1500, 2200, and 3300 km2) using available water table maps produced by different methods at different times in the area of High Plains Aquifer in Nebraska. One region is located in the largely undeveloped Nebraska Sand Hills area, while the other two regions are located at a transition zone from Sand Hills to loess-covered area and include areas where groundwater is used for irrigation. Obtained recharge rates are consistent with other independent estimates. The approach is useful and robust diagnostic tool for preliminary estimates of recharge rates, evaluation of the quality of groundwater table maps, identification of priority areas for further aquifer characterization and expansion of groundwater monitoring networks prior to using more detailed methods. Includes supplemental materials

Pumping-Induced Drawdown and Stream Depletion in a Leaky Aquifer System

The impact of ground water pumping on nearby streams is often estimated using analytic models of the interconnected stream-aquifer system. A common assumption of these models is that the pumped aquifer is underlain by an impermeable formation. A new semianalytic solution for drawdown and stream depletion has been developed that does not require this assumption. This solution shows that pumping-induced flow (leakage) through an underlying aquitard can be an important recharge mechanism in many stream-aquifer systems. The relative importance of this source of recharge increases with the distance between the pumping well and the stream. The distance at which leakage becomes the primary component of the pumping-induced recharge depends on the specific properties of the aquifer, aquitard, and streambed. Even when the aquitard is orders of magnitude less transmissive than the aquifer, leakage can be an important recharge mechanism because of the large surface area over which it occurs. Failure to consider aquitard leakage can lead to large overestimations of both the drawdown produced by pumping and the contribution of stream depletion to the pumping-induced recharge. The ramifications for water resources management and water rights adjudication can be significant. A hypothetical example helps illustrate these points and demonstrates that more attention should be given to estimating the properties of aquitards underlying stream-aquifer systems. The solution presented here should serve as a relatively simple but versatile tool for practical assessments of pumping-induced stream-aquifer interactions. However, this solution should not be used for such assessments without site-specific data that indicate pumping has induced leakage through the aquitard

Dynamic interpretation of slug tests in highly permeable aquifers

This is the published version. Copyright American Geophysical Union[1] Considerable progress has been made in developing a theoretical framework for modeling slug test responses in formations with high hydraulic conductivity K. However, several questions of practical significance remain unresolved. Given the rapid and often oscillatory nature of test responses, the traditional hydrostatic relationship between the water level and the transducer-measured head in the water column may not be appropriate. A general dynamic interpretation is proposed that describes the relationship between water level response and transducer-measured head. This theory is utilized to develop a procedure for transforming model-generated water level responses to transducer readings. The magnitude of the difference between the actual water level position and the apparent position based on the transducer measurement is a function of the acceleration and velocity of the water column, test geometry, and depth of the transducer. The dynamic approach explains the entire slug test response, including the often-noted discrepancy between the actual initial water level displacement and that measured by a transducer in the water column. Failure to use this approach can lead to a significant underestimation of K when the transducer is a considerable distance below the static water level. Previous investigators have noted a dependence of test responses on the magnitude of the initial water level displacement and have developed various approximate methods for analyzing such data. These methods are re-examined and their limitations clarified. Practical field guidelines are proposed on the basis of findings of this work. The soundness of the dynamic approach is demonstrated through a comparison of K profiles from a series of multilevel slug tests with those from dipole-flow tests performed in the same wells

Prediction of Biome-Specific Potential Evapotranspiration in Mongolia under a Scarcity of Weather Data

We propose practical guidelines to predict biome-specific potential evapotranspiration (ETp) from the knowledge of grass-reference evapotranspiration (ET0) and a crop coefficient (Kc) in Mongolia. A paucity of land-based weather data hampers use of the Penman–Monteith equation (FAO-56 PM) based on the Food and Agriculture Organization (FAO) guidelines to predict daily ET0. We found that the application of the Hargreaves equation provides ET0 estimates very similar to those from the FAO-56 PM approach. The Kc value is tabulated only for crops in the FAO-56 guidelines but is unavailable for steppe grasslands. Therefore, we proposed a new crop coefficient, Kc adj defined by (a) net solar radiation in the Gobi Desert (Kc adjD) or (b) leaf area index in the steppe region (Kc adjS) in Mongolia. The mean annual ETp obtained using our approach was compared to that obtained by FAO-56 guidelines for forages (not steppe) based on tabulated Kc values in 41 locations in Mongolia. We found the differences are acceptable (RMSE of 0.40 mm d-1) in northern Mongolia under high vegetation cover but rather high (RMSE of 1.69 and 2.65 mm d-1) in central and southern Mongolia. The FAO aridity index (AI) is empirically related to the ETp/ET0 ratio. Approximately 80% and 54% reduction of ET0 was reported in the Gobi Desert and in the steppe locations, respectively. Our proposed Kc adj can be further improved by considering local weather data and plant phenological characteristics

A Model of IceWedge Polygon Drainage in Changing Arctic Terrain

As ice wedge degradation and the inundation of polygonal troughs become increasingly common processes across the Arctic, lateral export of water from polygonal soils may represent an important mechanism for the mobilization of dissolved organic carbon and other solutes. However, drainage from ice wedge polygons is poorly understood. We constructed a model which uses cross-sectional flow nets to define flow paths of meltwater through the active layer of an inundated low-centered polygon towards the trough. The model includes the eects of evaporation and simulates the depletion of ponded water in the polygon center during the thaw season. In most simulations, we discovered a strong hydrodynamic edge eect: only a small fraction of the polygon volume near the rim area is flushed by the drainage at relatively high velocities, suggesting that nearly all advective transport of solutes, heat, and soil particles is confined to this zone. Estimates of characteristic drainage times from the polygon center are consistent with published field observations

Impact of grassland conversion to forest on groundwater recharge in the Nebraska Sand Hills

Study region: Nebraska National Forest in the High Plains Aquifer, Nebraska Sand Hills, U.S.A. Study focus: This research aimed to investigate the effects of grassland conversions to forest on recharge rates in a century-old experimental forest. The DiffeRential Evolution Adaptive Metropolis (DREAMZS) global optimization algorithm was used to calibrate the effective soil hydraulic parameters from observed soil moisture contents for 220 cm deep uniform soil profiles. The historical recharge rates were then estimated by applying the numerical model HYDRUS 1-D for simulation of two plots representing grasslands and dense pine forest conditions. New hydrological insights: The results indicate that conversion from grasslands to dense pine forests led to vegetation induced changes in soil hydraulic properties, increased rooting depth, and greater leaf area index, which together altered the water budget considerably. The impacts of land use change, expressed in percent of gross precipitation, include a 7% increase in interception associated with an increase in leaf area index, a nearly 10% increase in actual evapotranspiration, and an overall reduction of groundwater recharge by nearly 17%. Simulated average annual recharge rates decreased from 9.65 cm yr−1 in the grassland to 0.07 cm yr−1 in the pine plot. These outcomes highlight the significance of the grassland ecology for water resources, particularly groundwater recharge, in the Nebraska Sand Hills and the overall sustainability and vitality of the High Plains Aquifer