A Comparative Study of Specific Yield Determinations for a Shallow Sand Aquifer

Using the type-curve methods of Boulton (1963) and Neuman (1972), and comparisons, at various times, of the cumulative volume of water pumped to the volume of the water-table drawdown cone (volume-balance method), values of specific yield were obtained from pumping test data from numerous piezometers in an unconfined sand aquifer. The long-term value of specific yield for the aquifer was determined from measurements of the laboratory drainage curve of the aquifer material. The volume-balance method gave specific yield values of 0.02, 0.05, 0.12, 0.20, 0.23, and 0.25 at times of 0.25, 0.66, 10, 26, 45, and 65 hours, respectively, indicating a gradual increase in specific yield and an asymptotic approach to the long-term value of 0.30 determined from the laboratory method. The type-curve methods provided values of 0.07 and 0.08, which correspond to the volume-balance values at early times, but which are less than one-third of the value obtained from the laboratory method and from the volume-balance method applied at the end of the pumping test (2.7 days). The type-curve procedures therefore provide unrealistically low values of specific yield for application to problems concerning the long-term yield characteristics of the aquifer. The observed trend towards increasing values of specific yield with increasing duration of pumping, and the vertical hydraulic head profiles that were measured during the pumping test indicate that both delayed drainage from above the water table and downward hydraulic gradients in the saturated zone can be important hydraulic effects contributing to the delayed-drawdown segment that is characteristic of time-drawdown graphs for unconfined aquifers.     

Application Procedures for the Electromagnetic Borehole Flowmeter in Shallow Unconfined Aquifers

It is increasingly common for the electromagnetic borehole flowmeter (EBF) to he used to measure hydraulic conductivity (K) distributions in subsurface flow systems. Past applications involving the EBF have been made mostly in confined aquifers (Kabala 1994; Boman et al. 1997; Podgorney and Ritzi 1997; Ruud and Kabala 1997a, 1997b; Flach et al. 2000), and it has been common to set up a flow field around a test well using a small pump that is located near the top of the well screen (Mob, and Young 1993). In thin, unconfined aquifers that exhibit ground water tables near the ground surface and that undergo drawdown during pumping, such a configuration can be problematical because pumping and associated drawdown may effectively isolate the upper portion of the aquifer from the flowmeter. In these instances, a steady-state flow field in the vicinity of the test well may be created using injection rather than pumping, allowing for testing in the otherwise isolated upper portion of the aquifer located near the initial water table position. Using procedures developed by Molz and Young (1993), which were modified for an injection mode application, testing was conducted to determine whether or not the injection mode would provide useful information in a shallow, unconfined aquifer that required the collection of data near the initial water table position. Results indicated that the injection mode for the EBF was well suited for this objective.     

Influence of aquifer heterogeneity and return flow on pumping test data interpretation

Analytical solutions of drawdown in unconfined aquifers are widely applied for determining the specific yield, S_y, and the horizontal and the vertical hydraulic conductivity K_r and K_z, respectively. In many previous studies, estimates of S_y and K_z were observed to be highly variable and physically unrealistic. This has been attributed to the conceptualization of flow above the declining water table and aquifer heterogeneity in the applied models. We present the analysis of time-drawdown data from a pumping test instrumented with depth-differentiated observation piezometers arranged in clusters. Applying homogeneous anisotropic aquifer models in combination with nonlinear least squares parameter identification techniques, the data were analyzed in different groups: analysis of data from individual piezometer clusters and simultaneous analysis of the entire data set from all piezometer clusters (global analysis). From the cluster analyses, estimates of S_y and K_z exhibit large variances and depart from a priori estimates inferred from the hydrostratigraphy. Parameter estimates from the global analysis do not fall within the parameter bounds (minimum and maximum values) defined by the cluster analyses. While heterogeneity appears to be the important reason for large parameter variances, we discuss the influence of rarely considered aquifer return flow on drawdown and the inconsistent results from the cluster and global analyses. We corroborate our findings with data on hydraulic gradients, slug test data, and results from the application of a more realistic numerical flow model.     

Importance of the vadose zone in analyses of unconfined aquifer tests

Analytical models commonly used to interpret unconfined aquifer tests have been based on upper-boundary (water table) conditions that do not adequately address effects of time-varying drainage from the vadose zone. As a result, measured and simulated drawdown data may not agree and hydraulic parameters may be inaccurately estimated. A 72-hour aquifer test conducted in Cape Cod, Massachusetts, in a slightly heterogeneous, coarse-grained, glacial outwash deposit was found to be a good candidate for testing models with different upper-boundary conditions. In general, under the commonly invoked assumption of instantaneous drainage, measured and simulated drawdowns were found to agree with one another only at late time and early time. In the intermediate-time range, because of delayed drainage, measured drawdowns always exceeded simulated values, most noticeably in piezometers located near the water table. To reduce these discrepancies, an analytical model was developed that can fully account for time-varying drainage given that the aquifer is not strongly heterogeneous. The approach is flexible as the model, which makes use of empirical relations, does not constrain drainage to follow any particular functional relation. By this approach, measured and simulated drawdowns agree over the complete time range, and the estimated parameters are consistent with prior studies and with what is known about the aquifer geometry, stratigraphy, and composition. By properly accounting for vadose zone drainage, it was found that realistic estimates of all hydraulic parameters, including specific yield, could be obtained with or without the use of late-time data.     

Aquifer Properties Determined from Two Analytical Solutions

In the analysis of pumping test data, the quality of the determined aquifer parameters can be greatly improved by using a proper model of the aquifer system. Moench (1995) provided an analytical solution for flow to a well partially penetrating an unconfined aquifer. His solution, in contrast to the Neuman solution (1974), accounts for the noninstantaneous decline of the water table (delayed yield). Consequently, the calculated drawdown in these two solutions is different under certain circumstances, and this difference may therefore affect the computation of aquifer properties from pumping test data. This paper uses an inverse computational method to calculate four aquifer parameters as well as a delayed yield parameter, α₁ from pumping test data using both the Neuman (1974) and Moench (1995) solutions. Time-drawdown data sets from a pumping test in an unconfined alluvial aquifer near Grand Island, Nebraska, were analyzed. In single-well analyses, horizontal hydraulic conductivity values derived from the Moench solution are lower, but vertical hydraulic conductivity values are higher than those calculated from the Neuman solution. However, the hydraulic conductivity values in composite-well analyses from both solutions become very close. Furthermore, the Neuman solution produces similar hydraulic conductivity values in the single-well and composite-well analyses, but the Moench solution does not. While variable α₁, seems to play a role in affecting the computation of aquifer parameters in the single-well analysis, a much smaller effect was observed in the composite-well analysis. In general, specific yield determined using the Moench solution could be slightly higher than the values from the Neuman solution; however, they are still lower than the realistic values for sand and gravel aquifers.     

Influence of calibration methodology on ground water flow predictions

We constructed a numerical model of transient ground water flow and solute transport for a portion of the Biscayne Aquifer in Florida, and calibrated the model with three different combinations of data from a 193-day period: head (h) data alone, data on h and ground water discharge to a canal (q), and data on h, q, and ground water chloride concentration (C). We used each of the three calibrated models to predict h and q during a 182-day test period separate from the calibration period. All three calibrated models predicted h equally well during the test period (r = 0.95, where r = 1 indicates perfect agreement between measured and simulated values), though the model calibrated on h alone had significantly different parameter values than the other two models. Predictions of q during the test period depended on calibration methodology; models calibrated with multiple targets simulated q more accurately than the model calibrated on h alone (r = 0.79 compared to r = 0.49). Based on the results of these simulations, we conclude: (1) Post-calibration prediction is important in assessing the value of different data types in automated calibration; (2) inverse-solution uniqueness is not a requirement for accurate h predictions; (3) relatively simple models can predict with reasonable accuracy transient ground water flow in a complex aquifer, and parameters governing this prediction can be estimated by nonlinear regression methods that incorporate both h and q data; (4) addition of C data to the calibration did not improve model predictive capacity because the information in the C data was similar to that in the q data, from the perspective of model calibration (the subsurface chemical signal in question was controlled mainly by seepage of high-chloride canal water into the low-chloride ground water system).     

An Analysis of Low-Flow Ground Water Sampling Methodology

Low-flow ground water sampling methodology can minimize well disturbance and aggravated colloid transport into samples obtained from monitoring wells. However, in low hydraulic conductivity formations, low-flow sampling methodology can cause excessive drawdown that can result in screen desaturation and high ground water velocities in the vicinity of the well, causing unwanted colloid and soil transport into ground water samples taken from the well. Ground water velocities may increase several fold above that of the natural setting. To examine the drawdown behavior of a monitoring well, mathematical relationships can be developed that allow prediction of the steady-state drawdown for constant low-flow pumping rates based on well geometry and aquifer properties. The equations also estimate the time necessary to reach drawdown equilibrium. These same equations can be used to estimate the relative contribution of water entering a sampling device from either the well standpipe or the aquifer. Such equations can be useful in planning a low-flow sampling program and may suggest when to collect a water sample. In low hydraulic conductivity formations, the equations suggest that drawdown may not stabilize for well depths, violating the minimal drawdown requirement of the low-flow technique. In such cases, it may be more appropriate to collect a slug or passive sample from the well screen, under the assumption that the water in the well screen is in equilibrium with the surrounding aquifer.     

Parameter Estimation/Sensitivity Analysis for an Aquifer Test with Skin Effect

Aquifer information carried by aquifer test data may be affected by the presence of a finite thickness skin around the wellbore. The mathematical treatment for an aquifer accounting for the skin zone can be characterized by five parameters, that is, the outer radius of the skin zone and the transmissivity and storativity for each of the skin and aquifer zones. Sensitivity analysis was performed to examine the ground water flow behavior in the skin and aquifer zones in terms of the constant-head test (CHT) data. The simulated annealing procedure was applied to simultaneously determine the skin and aquifer parameters from the analysis of CHT data. Toward the previously mentioned goals, four suites of CHT data were analyzed in this article. The analyses of wellbore flow rate at the test well and the specific drawdown at the observation well gave accurate estimates for the skin and aquifer parameters, respectively. Only the skin thickness and both the skin and the aquifer diffusivities could be accurately estimated from the analysis of drawdown data in the observation well. The estimates for all skin and aquifer parameters from the composite analysis of flow rate and drawdown data were the most accurate. The results of sensitivity analyses and parameter estimations provide instructive references in the analysis of the skin-affected CHT data.     

Estimating aquifer properties using groundwater hydrograph modelling

Estimation of aquifer hydraulic properties is essential for predicting the response of an aquifer to extractions and hence estimating the availability of the groundwater resources. Aquifer tests are commonly used for the estimation of aquifer properties; however, they can be expensive and often only characterize the short‐term response of the aquifer. This paper presents a time series modelling approach to estimating aquifer hydraulic properties. It is applied to 42 bores monitoring an unconfined aquifer within an irrigation region of south‐eastern Australia, and the resulting probabilistic estimate of hydraulic properties are evaluated against pumping test estimates. It is demonstrated that the time series modelling can provide a reliable estimate of the hydraulic properties that are typical of a very long‐term pumping test. Furthermore, the application of the time series modelling to 42 bores provided novel insights into the aquifer heterogeneity. We encourage others to further test the approach and the source code is available from: http://www.mathworks.com/matlabcentral/fileexchange/48546‐peterson‐tim‐j‐groundwater‐statistics‐toolbox Copyright © 2015 John Wiley & Sons, Ltd.     

Analytical methods for transient flow to a well in a confined-unconfined aquifer

Concurrent existence of confined and unconfined zones of an aquifer can arise owing to ground water withdrawal by pumping. Using Girinskii's potential function, Chen (1974, 1983) developed an approximate analytical solution to analyze transient ground water flow to a pumping well in an aquifer that changes from an initially confined system to a system with both unconfined and confined regimes. This article presents the details of the Chen model and then compares it with the analytical model developed by Moench and Prickett (1972) for the same problem. Hypothetical pumping test examples in which the aquifer undergoes conversion from confined to water table conditions are solved by the two analytical models and also a numerical model based on MODFLOW. Comparison of the results suggests that the solutions of the Chen model give better results than the Moench and Prickett model except when the radial distance is very large or aquifer thickness is large compared with drawdown.     

Quantifying the Effects of Well Development in Unconsolidated Material

A fully instrumented physical model was designed and built to reproduce development by surging and monitor its effects during surging and after development. The model simulates a horizontal layer in a confined aquifer with control of vertical overburden pressure. An automatic apparatus produced development by surging in successive phases up to 24 hours. Aquifer tests in steady-state conditions were performed between successive phases. The paper reports the main results of three experiments performed with Johnson screens 200 mm in diameter; they had slot sizes between the D₅₄ and D₇₀ of the aquifer soil. This soil was placed under controlled conditions, and initial homogeneity was obtained as confirmed by initial control tests. Pore pressures (and thus hydraulic heads) were continuously monitored during development phases and aquifer tests by 22 electronic piezometers at distances between 0 and 1 m from the screen. These piezometers measured water pressures every 0.1 s when required. Solid particles passing through the screen were recovered to study the solid yield and the gradation of particles. Positive and negative values of local gradients reached values up to 400 close to the screen at the beginning of development and decreased with time of development. These high values produced high seepage forces displacing particles in the aquifer. The well yield was increased by a factor of 6 after development. These model test results confirmed empirical criteria on entrance velocity, internal stability criteria, and field values of "sand" production by development. In addition, they enabled a quantification of skin effects to be considered in interpreting an aquifer test.     

A method for evaluating horizontal well pumping tests

Predicting the future performance of horizontal wells under varying pumping conditions requires estimates of basic aquifer parameters, notably transmissivity and storativity. For vertical wells, there are well-established methods for estimating these parameters, typically based on either the recovery from induced head changes in a well or from the head response in observation wells to pumping in a test well. Comparable aquifer parameter estimation methods for horizontal wells have not been presented in the ground water literature. Formation parameter estimation methods based on measurements of pressure in horizontal wells have been presented in the petroleum industry literature, but these methods have limited applicability for ground water evaluation and are based on pressure measurements in only the horizontal well borehole, rather than in observation wells. This paper presents a simple and versatile method by which pumping test procedures developed for vertical wells can be applied to horizontal well pumping tests. The method presented here uses the principle of superposition to represent the horizontal well as a series of partially penetrating vertical wells. This concept is used to estimate a distance from an observation well at which a vertical well that has the same total pumping rate as the horizontal well will produce the same drawdown as the horizontal well. This equivalent distance may then be associated with an observation well for use in pumping test algorithms and type curves developed for vertical wells. The method is shown to produce good results for confined aquifers and unconfined aquifers in the absence of delayed yield response. For unconfined aquifers, the presence of delayed yield response increases the method error.     

Monitoring of drawdown pattern during pumping in an unconfined physical aquifer model

In this study, we attempted to analyse a drawdown pattern around a pumping well in an unconfined sandy gravelly aquifer constructed in a laboratory tank by means of both experimental and numerical modelling of groundwater flow. The physical model consisted of recharge, aquifer and discharge zones. Permeability and specific yield of the aquifer material were determined by Dupuit approximation under steady‐state flow and stepwise gravitational drainage of groundwater, respectively. The drawdown of water table in pumping and neighbouring observation wells was monitored to investigate the effect of no‐flow boundary on the drawdown pattern during pumping for three different boundary conditions: (i) no recharge and no discharge with four no‐flow boundaries (Case 1); (ii) no recharge and reservoir with three no‐flow boundaries (Case 2); (iii) recharge and discharge with two no‐flow boundaries (Case 3). Based on the aquifer parameters, numerical modelling was also performed to compare the simulated drawdown with that observed. Results showed that a large difference existed between the simulated drawdown and that observed in wells for all cases. The reason for the difference could be explained by the formation of a curvilinear type water table between wells rather than a linear one due to a delayed response of water table in the capillary fringe. This phenomenon was also investigated from a mass balance study on the pumping volume. The curvilinear type of water table was further evidenced by measurement of water contents at several positions in the aquifer between wells using time domain reflectometry (TDR). This indicates that the existing groundwater flow model applicable to an unconfined aquifer lacks the capacity to describe a slow response of water table in the aquifer and care should be taken in the interpretation of water table formation in the aquifer during pumping. Copyright © 2001 John Wiley & Sons, Ltd.     

Electrical Resistivity Characterization of Anisotropy in the Biscayne Aquifer

Electrical anisotropy occurs when electric current flow varies with azimuth. In porous media, this may correspond to anisotropy in the hydraulic conductivity resulting from sedimentary fabric, fractures, or dissolution. In this study, a 28‐electrode resistivity imaging system was used to investigate electrical anisotropy at 13 sites in the Biscayne Aquifer of SE Florida using the rotated square array method. The measured coefficient of electrical anisotropy generally ranged from 1.01 to 1.12 with values as high as 1.36 found at one site. The observed electrical anisotropy was used to estimate hydraulic anisotropy (ratio of maximum to minimum hydraulic conductivity) which ranged from 1.18 to 2.83. The largest values generally were located on the Atlantic Coastal Ridge while the lowest values were in low elevation areas on the margin of the Everglades to the west. The higher values of anisotropy found on the ridge may be due to increased dissolution rates of the oolitic facies of the Miami formation limestone compared with the bryozoan facies to the west. The predominate trend of minimum resistivity and maximum hydraulic conductivity was E‐W/SE‐NW beneath the ridge and E‐W/SW‐NE farther west. The anisotropy directions are similar to the predevelopment groundwater flow direction as indicated in published studies. This suggests that the observed anisotropy is related to the paleo‐groundwater flow in the Biscayne Aquifer.     

Improving Groundwater Model Calibration with Repeat Microgravity Measurements

Groundwater-flow models depend on hydraulic head and flux observations for evaluation and calibration. A different type of observation—change in storage measured using repeat microgravity—can also be used for parameter estimation by simulating the expected change in gravity from a groundwater model and including the observation misfit in the objective function. The method is demonstrated using new software linked to MODFLOW input and output files and field data from the vicinity of the All American Canal in southeast California, USA. Over a 10-year period following lining of the previously highly permeable canal with concrete, gravity decreased by over 100 μGal (equivalent to about 2.5 m of free-standing water) at some locations as seepage decreased and the remnant groundwater mound dissipated into the aquifer or was removed by groundwater pumping. Simulated gravity from a MODFLOW model closely matched observations, and repeat microgravity data proved useful for constraining both hydraulic conductivity and specific yield estimates. Specific yield estimated using the infinite-horizontal slab approximation agreed well with model-derived values, and the departure from the linear, flat-water-table approximation was small, less than 2%, despite relatively large and dynamic water-table slope. First-order second-moment parameter uncertainty analysis shows reduction in uncertainty for all hydraulic conductivity and specific yield parameter estimates with the addition of repeat microgravity data, as compared to drawdown data alone.     

An inverse procedure to estimate transmissivity from heads and SP signals

Experimental hydraulic heads and electrical (self-potential) signals associated with a pumping test were used in an inverse model to estimate the transmissivity distribution of a real aquifer. Several works reported in the literature show that there is a relatively good linear relationship between the hydraulic heads in the aquifer and electrical signals measured at the ground surface. In this experimental test field, first, the current coupling coefficient was determined by the best fit between experimental and modeled self-potential signals at the end of the pumping phase. Soon afterward, with the hydraulic heads obtained from the self-potential signals, the transmissivity distribution of the aquifer was conditioned by means an inverse model based on the successive linear estimator (SLE). To further substantiate the estimated T field from the SLE analysis, we analyzed the drawdown rate, the derivative of the drawdown with respect to the ln(t), because the drawdown rate is highly sensitive to the variability in the transmissivity field. In our opinion, these results show that self-potential signals allow the monitoring of subsurface flow in the course of pumping experiments, and that electrical potentials serve as a good complement to piezometric observations to condition and characterize the transmissivity distribution of an aquifer.     

A new method for aquifer system identification and parameter estimation

The standard practice for assessing aquifer parameters is to match groundwater drawdown data obtained during pumping tests against theoretical well function curves specific to the aquifer system being tested. The shape of the curve derived from the logarithmic time derivative of the drawdown data is also very frequently used as a diagnostic tool to identify the aquifer system in which the pumping test is being conducted. The present study investigates the incremental area method (IAM) to serve as an alternative diagnostic tool for the aquifer system identification as well as a supplement to the aquifer parameter estimation procedure. The IAM based diagnostic curves for ideal confined, leaky, bounded and unconfined aquifers have been derived as part of this study, and individual features of the plots have been identified. These features were noted to be unique to each aquifer setting, which could be used for rapid evaluation of the aquifer system. The effectiveness of the IAM methodology was investigated by analyzing field data for various aquifer settings including leaky, unconfined, bounded and heterogeneous conditions. The results showed that the proposed approach is a viable method for use as a diagnostic tool to identify the aquifer system characteristics as well as to support the estimation of the hydraulic parameters obtained from standard curve matching procedures. Copyright © 2012 John Wiley & Sons, Ltd.     

NMR Imaging of Fluid Exchange between Macropores and Matrix in Eogenetic Karst

Sequential time-step images acquired using nuclear magnetic resonance (NMR) show the displacement of deuterated water (D₂O) by fresh water within two limestone samples characterized by a porous and permeable limestone matrix of peloids and ooids. These samples were selected because they have a macropore system representative of some parts of the eogenetic karst limestone of the Biscayne Aquifer in southeastern Florida. The macroporosity, created by the trace fossil Ophiomorpha , is principally well connected and of centimeter scale. These macropores occur in broadly continuous stratiform zones that create preferential flow layers within the hydrogeologic units of the Biscayne. This arrangement of porosity is important because in coastal areas, it could produce a preferential pathway for salt water intrusion. Two experiments were conducted in which samples saturated with D₂O were placed in acrylic chambers filled with fresh water and examined with NMR. Results reveal a substantial flux of fresh water into the matrix porosity with a simultaneous loss of D₂O. Specifically, we measured rates upward of 0.001 mL/h/g of sample in static conditions, and perhaps as great as 0.07 mL/h/g of sample when fresh water continuously flows past a sample at velocities less than those found within stressed areas of the Biscayne. These experiments illustrate how fresh water and D₂O, with different chemical properties, migrate within one type of matrix porosity found in the Biscayne. Furthermore, these experiments are a comparative exercise in the displacement of sea water by fresh water in the matrix of a coastal, karst aquifer since D₂O has a greater density than fresh water.     

Identifying connections in a fractured rock aquifer using ADFTs

Fractured rock aquifers are difficult to characterize because of their extremely heterogeneous nature. Developing an understanding of fracture network hydraulic properties in these aquifers is difficult and time consuming, and field testing techniques for determining the location and connectivity of fractures in these aquifers are limited. In the Clare Valley, South Australia, well interference is an important issue for a major viticultural area that uses a fractured aquifer. Five fracture sets exist in the aquifer, all dipping > 25 degrees . In this setting, we evaluate the ability of steady-state asymmetric dipole-flow tests (ADFTs) to determine the connections between a test well and a set of piezometers. The procedure involves dividing a test well into two chambers using a single packer and pumping fluid from the upper chamber to the lower chamber. By conducting a series of tests at different packer elevations, an "input" signal is generated in fracture zones connected to the test well. By monitoring the "output" response of the hydraulic dipole field at piezometers, the connectivity of the fractures between the test well and piezometers can be determined. Results indicate the test well used in this study is connected in a complex three-dimensional geometry, with drawdown occurring above and below areas of potentiometric buildup. The ADFT method demonstrates that the aquifer evaluated in this study cannot be modeled effectively on the well scale using continuum flow models.