DYNAMIC RESPONSE OF AQUIFER SYSTEMS TO LOCALIZED RECHARGE1

ABSTRACT. The response of stream-unconfined aquifer systems to localized recharge is investigated by means of a two-dimensional finite element model. A variational approach is used in conjunction with the finite element method to solve the ground water flow equation. Linear approximated triangular elements are used to calculate the hydraulic head distribution in the flow region. The Crank-Nicholson centered scheme of numerical integration is employed to approximate the time derivative in the flow equation. A computer program is developed to calculate the hydraulic head distribution in the flow region. Solutions provided by the finite element model should prove useful in the evaluation of quantitative and qualitative changes in aquifer systems due to natural or artificial recharge. In addition, they should prove useful in the study of irrigation and drainage problems.     

Seasonal and Long‐Term Variations in Hydraulic Head in a Karstic Aquifer: Roswell Artesian Basin,New Mexico1

Abstract: Water levels in the karstic San Andres limestone aquifer of the Roswell Artesian Basin, New Mexico, display significant variations on a variety of time scales. Large seasonal fluctuations in hydraulic head are directly related to the irrigation cycle in the Artesian Basin, lower in summer months and higher in winter when less irrigation occurs. Longer‐term variations are the result of both human and climatic factors. Since the inception of irrigated farming more than a century ago, over‐appropriation of water resources has caused water levels in the artesian aquifer to fall by as much as 70 m. The general decline in hydraulic head began to reverse in the mid‐1980s due to a variety of conservation measures, combined with a period of elevated rainfall toward the end of the 20th Century.     

ANALYSIS OF PUMPING TESTS: SIGNIFICANCE OF WELL DIAMETER,PARTIAL PENETRATION,AND NOISE1

ABSTRACT: The nonlinear least squares (NLS) method was applied to pumping and recovery aquifer test data in confined and unconfined aquifers with finite diameter and partially penetrating pumping wells, and with partially penetrating piezometers or observation wells. It was demonstrated that noiseless and moderately noisy drawdown data from observation points located less than two saturated thicknesses of the aquifer from the pumping well produced an exact or acceptable set of parameters when the diameter of the pumping well was included in the analysis. The accuracy of the estimated parameters, particularly that of specific storage, decreased with increases in the noise level in the observed drawdown data. With consideration of the well radii, the noiseless drawdown data from the pumping well in an unconfined aquifer produced good estimates of horizontal and vertical hydraulic conductivities and specific yield, but the estimated specific storage was unacceptable. When noisy data from the pumping well were used, an acceptable set of parameters was not obtained. Further experiments with noisy drawdown data in an unconfined aquifer revealed that when the well diameter was included in the analysis, hydraulic conductivity, specific yield and vertical hydraulic conductivity may be estimated rather effectively from piezometers located over a range of distances from the pumping well. Estimation of specific storage became less reliable for piezometers located at distances greater than the initial saturated thickness of the aquifer. Application of the NLS to field pumping and recovery data from a confined aquifer showed that the estimated parameters from the two tests were in good agreement only when the well diameter was included in the analysis. Without consideration of well radii, the estimated values of hydraulic conductivity from the pumping and recovery tests were off by a factor of four.     

EFFECTS OF A RECREATIONAL IMPOUNDMENT ON NEARBY NATURAL LAKES: A CASE STUDY1

ABSTRACT: Completion of a 1270 acre recreational impoundment (Legend Lake) in the glacial sands of Menominee County, Wisconsin, produced geochemical and hydrologic alterations in some nearby natural lakes. The impoundment was produced by the construction of three dams, one of which proved to be temporary, connecting 9 natural lakes and ponds of 383 acres with 951 acres of flooded lands. Water levels were raised 3–15 feet within the impounded area. Much of the flooded area was peat rich wetland associated with the prior drainage. Water depths are less than 15 feet in 70% of the impoundment. Three seepage lakes, located less than 1/2 mile from the impoundment, experienced shoreline flooding, shoreline and soil erosion, some tree kills, and increased turbidity. These lakes also experienced concentration increases in several chemical constituents which indicate an influx of impoundment water through a regional alternation in the groundwater flow paths. The three lakes were connected by canals, and a 2.3 cfs gravity drain with an auxiliary pumping station was built to return excess water to the outflow of the impoundment. Future projects of this type would benefit from a more extensive hydrologic and geochemical analysis prior to initiation. Had environmental assessments been required at the time of this development, as they now are in Wisconsin for similar projects, some of the problems encountered might have been alleviated.     

Management of the Spring Snowmelt Recession in Regulated Systems

In an effort to restore predictable ecologically relevant spring snowmelt recession flow patterns in rivers regulated by dams, this study defined a methodology by which spring flow regimes can be modeled in regulated systems from the quantifiable characteristics of spring snowmelt recessions in unregulated rivers. An analysis of eight unregulated rivers across the Sierra Nevada mountain range in California found that unregulated systems behaved similarly with respect to seasonal spring patterns and recession limb curvature, and thus prescribed flows could be designed in a manner that mimics those predictable characteristics. Using the methodology to quantify spring recession flows in terms of a daily percent decrease in flow, a series of flow recession scenarios were created for application in an existing hydrodynamic model for the regulated Rubicon River. The modeling results showed that flow recessions with slow ramping rates similar to those observed in unregulated rivers (less than 10% per day) were likely to be protective of native aquatic species, such as the Foothill yellow‐legged frog, while flows that receded at greater rates would likely result in desiccation of egg masses and potential stranding of tadpoles and fry. Furthermore, recession rates of less than 10% per day provided the most spatially diverse hydraulic habitat in the modeled domain for an appropriate duration in spring to support all native species guilds and maximize aquatic biodiversity.     

RECYCLING INPUT DATA DURING ANALYTIC ELEMENT MODELING NEAR INDIANAPOLIS,INDIANA1

ABSTRACT: Overlapping and adjacent ground water investigations are common in areas where aquifers are threatened by industrial development. In the Indianapolis area in Marion County, Indiana, a patchwork of ground water flow models have been used during the past 20 years to evaluate ground water resources and to determine the effects of local contamination. In every case these ground water models were constructed from scratch. Site specific finite difference grids or finite element meshes inhibit the direct reuse of input data when the area of interest shifts. Because the aquifer is not discretized into a grid or mesh with analytic element models, there are unique opportunities for direct reuse of model input data. In two applications of this principle we illustrate how the newly emerging analytic element method allows a fairly straightforward reuse of model input data from previous models in the same general area. In analytic element models of Central Indiana, streams and their tributaries are represented in different resolutions. Input data items of several modeling studies are stored and cataloged on disk in such a manner that they can be selectively retrieved by a data management program PREPRO. In this manner, a new ground water model can be set up quickly with input data which have been previously defined and tested during model calibration.     

Effects of near-surface hydraulic gradients on nitrate and phosphorus losses in surface runoff

Phosphorous (P) and nitrogen (N) in runoff from agricultural fields are key components of nonpoint-source pollution and can accelerate eutrophication of surface waters. A laboratory study was designed to evaluate effects of near-surface hydraulic gradients on P and N losses in surface runoff from soil pans at 5% slope under simulated rainfall. Experimental treatments included three rates of fertilizer input (control [no fertilizer input], low [40 kg P ha(-1), 100 kg N ha(-1)], and high [80 kg P ha(-1), 200 kg N ha(-1)]) and four near-surface hydraulic gradients (free drainage [FD], saturation [Sa], artesian seepage without rain [Sp], and artesian seepage with rain [Sp + R]). Simulated rainfall of 50 mm h(-1) was applied for 90 min. The results showed that near-surface hydraulic gradients have dramatic effects on NO(3)-N and PO(4)-P losses and runoff water quality. Under the low fertilizer treatment, the average concentrations in surface runoff from FD, Sa, Sp, and Sp + R were 0.08, 2.20, 529.5, and 71.8 mg L(-1) for NO(3)-N and 0.11, 0.54, 0.91, and 0.72 mg L(-1) for PO(4)-P, respectively. Similar trends were observed for the concentrations of NO(3)-N and PO(4)-P under the high fertilizer treatment. The total NO(3)-N loss under the FD treatment was only 0.01% of the applied nitrogen, while under the Sp and Sp + R treatments, the total NO(3)-N loss was 11 to 16% of the applied nitrogen. These results show that artesian seepage could make a significant contribution to water quality problems.     

Dams in the Amazon: Belo Monte and Brazil’s Hydroelectric Development of the Xingu River Basin

Hydroelectric dams represent major investments and major sources of environmental and social impacts. Powerful forces surround the decision-making process on public investments in the various options for the generation and conservation of electricity. Brazil’s proposed Belo Monte Dam (formerly Kararaô) and its upstream counterpart, the Altamira Dam (better known by its former name of Babaquara) are at the center of controversies on the decision-making process for major infrastructure projects in Amazonia. The Belo Monte Dam by itself would have a small reservoir area (440 km²) and large installed capacity (11, 181.3 MW), but the Altamira/Babaquara Dam that would regulate the flow of the Xingu River (thereby increasing power generation at Belo Monte) would flood a vast area (6140 km²). The great impact of dams provides a powerful reason for Brazil to reassess its current policies that allocate large amounts of energy in the country’s national grid to subsidized aluminum smelting for export. The case of Belo Monte and the five additional dams planned upstream (including the Altamira/Babaquara Dam) indicate the need for Brazil to reform its environmental assessment and licensing system to include the impacts of multiple interdependent projects.     

OVERPUMPED ARTESIAN WELLS AMONG A WELL GROUP1

ABSTRACT. Overpumped artesian wells are defined and explained. Unified well equations which can be used for artesian, gravity and overpumped artesian wells are derived under both the steady and transient conditions of flow. The location of the boundary between artesian and gravity steady and transient flow systems in artesian aquifers, due to multiple wells among which overpumped artesian wells may exist, is explained.     

DETERMINATION OF AQUIFER PARAMETERS USING REGRESSION ANALYSIS1

ABSTRACT: Transmissivity and storativity of an aquifer are usually determined by analysis of steady or nonsteady pumping test data. The classical methods of nonsteady pumping test analysis are mostly graphical in nature and are, therefore, subject to errors of judgment in curve fitting, interpolating, and reading graphs and charts. A method is described here which does not require construction of graphs or use of charts and tables. The transmissivity and the storativity are calculated using regression analysis of the nonsteady time drawdown field data. The calculations can readily be performed on a hand held calculator. The procedure is described using four examples, and the results are compared with those obtained from graphical techniques. It is shown that the method is a viable alternative to the type curve solution of Theis or Straight line solution of Jacob for nonleaky artesian aquifers. However, the regression method poses problem in the cases of leaky artesian and water table aquifers.     

AN APPLICATION OF THE ANALYTIC ELEMENT METHOD IN MODELING FLORIDA EVERGLADES HYDROLOGY1

ABSTRACT: The large volumes of ground water that are discharged from the Everglades toward the Miami metropolitan area have historically posed a significant environmental water supply problem. In order to analyze the effects of seepage barriers on these subsurface outflows, the analytic element modeling code GFLOW was used to construct a ground water flow model of a region that includes a portion of the Everglades along with adjacent developed areas. The hydrology of this region can be characterized by a highly transmissive surficial aquifer in hydraulic contact with wetlands and canals. Calibration of the model to both wet and dry season conditions yielded satisfactory results, and it was concluded that the analytic element method is a suitable technique for modeling ground water flow in the Everglades environment. Finally, the model was used to evaluate the potential effectiveness of a subsurface barrier approximately two miles long for increasing water levels within the adjacent fringes of the Everglades National Park. It was found that the barrier had a negligible effect on water levels due to both its relatively short length and the high transmissivity of the surficial aquifer.     

Economics of reservoir sedimentation and sustainable management of dams

Accepted practice has been to design and operate reservoirs to fill with sediment, generating benefits from remaining storage over a finite period of time. The consequences of sedimentation and project abandonment are left to the future. This 'future' has already arrived for many existing reservoirs and most others will eventually experience a similar fate, thereby imposing substantial costs on society. Such costs could be avoided if sedimentation was minimized and dams were allowed to live forever. The fact that the world's inventory of suitable reservoir sites is limited provides an additional reason for encouraging the sustainable management of dams. This paper provides a framework for assessing the economic feasibility of sediment management strategies that would allow the life of dams to be prolonged indefinitely. Even if reduced accumulation or removal of sediment is technically possible, its economic viability is likely to depend on physical, hydrological and financial parameters. The model presented incorporates such factors and allows a characterization of conditions under which sustainable management would be desirable. The empirical implementation of the model draws upon the substantial amount of technical information available. We analyze the sustainability of reservoirs, with a focus on the trade-off between such sustainability and the short to medium term benefits which a reservoir is expected to produce. The results show that, for a very wide range of realistic parameter values, sustainable management of reservoirs is economically more desirable than the prevailing practice of forcing a finite reservoir life through excessive sediment accumulation.     

FACTORS INVOLVED IN EVALUATING GROUND WATER IMPACTS OF DEEP COAL MINE DRAINAGE1

ABSTRACT: The determination of probable ground water impacts of proposed deep coal mining is required as part of permit applications. Impact prediction generally involves well production test analysis and modeling of ground water systems associated with coal seams. Well production tests are often complicated due to the relatively low permeabilities of sandstones and shales of ground water systems. The effects of the release of water stored within finite diameter production wells must be considered. Well storage capacity appreciably affects early well production test time drawdown or time recovery data. Low pumping rates, limited cones of depression, and length of required pumping periods ate important well production test design factors. Coal seam ground water system models are usually multilayered and leaky artesian. Mine drafts partially penetrate the ground water system. Simulation of coal mine drainage often involves the horizontal permeability and storage coefficient of the coal seam zone, vertical permeabilities of sandstones and shales (aquitard) above and below the coal seam zone, and the hydrologic properties of the source bed above the aqultard overlying the coal seam zone. Ground water level declines in both the coal seam zone and source bed near land surface are necessary factors in impact analysis. An example of evaluation studies in southwest Indiana will illustrate factors involved in deep coal mine drainage modeling efforts.     

CRITICAL STATE OF SALT-WATER UPCONING BENEATH ARTESIAN DISCHARGE WELLS1

ABSTRACT: An approach, based on the realization of the vertical components of flow is presented to determine the free surface of gravity wells and the shape of salt-water upconing in artesian aquifers. The transitional stages from pumping fresh water to pumping salt water at the critical condition are discussed. Vertical hydraulic resistances, being determined and included in the flow domain, allow the use of the finite-difference approximation in a simple form. The suggested approach implies the use of either simple modeling or computing techniques.     

A DECISION THEORY APPROACH TO EARTH DAM CLASSIFICATION1

ABSTRACT: The fact that dams have failed indicates that there is still a risk involved, in spite of the major effort to ensure reliability. In seeking ways to reduce the risk, all aspects of the design, construction, monitoring, inspection, and rehabilitation of earth dams should be examined. Because dam classification is a central element, ways of minimizing uncertainty associated with classification need to be considered. A Bayesian approach to classification has numerous advantages over existing methods of earth dam classification because it directly evaluates the effects of uncertainties. A framework for incorporating Bayesian decision theory into the classification process is presented.     

MODELING PERCOLATION LOSSES FROM A PONDED FIELD UNDER VARIABLE WATER-TABLE CONDITIONS1

ABSTRACT: A numerical simulation model was developed to predict the vertical and lateral percolation losses from a ponded agricultural field. The two-dimensional steady-state unsaturated/ saturated flow equation was solved using the finite-difference technique. A constant ponding depth was maintained at the soil surface with different water table conditions in an application of the model for rice fields bordered by bunds. Field experiments were conducted for two different water table depths to collect data on the spatial distribution of volumetric soil-moisture content for model verification. The measured soil-moisture content values were found to be in close agreement with those predicted by the model. The sensitivity analysis of the model with selected hydrologic conditions shows that the model is most sensitive to the values of saturated hydraulic conductivity, but relatively less sensitive to water table depth, ponding depth, and evaporation rate from the soil surface. It implies that, in a ponded rice field condition, the lateral and vertical percolation losses are mostly governed by the hydraulic conductivity of the soil. The vertical percolation losses were almost equal to the saturated hydraulic conductivity values and, in most cases, these losses increased with deeper water table depths. The lateral percolation losses also increased with deeper water table depths; however, these losses were relatively small in comparison to the vertical percolation losses. The vertical and lateral percolation losses increased with the increase in ponding depths. The lateral percolation losses through the bund decreased when the evaporation losses increased from the soil surface. The results of this study indicate that the percolation losses from a ponded field may be predicted accurately for a wide range of soil and hydrological conditions when the values of hydraulic conductivity, evaporation rate, depth of ponding, and water table depth are accurately known.     

A FINITE ELEMENT MODEL OF CONTAMINANT MOVEMENT IN GROUNDWATER1

ABSTRACT. Transient, two-dimensional solutions are developed which describe the movement and distribution of a conservative substance in a stream-aquifer system. The solutions are obtained by solving sequentially the groundwater flow and mass transport equations. A variational approach in conjunction with the finite element method is used to solve the groundwater flow equation. Galerkin's approach coupled with the finite element method is used to solve the mass transport equation. Linear approximated triangular elements and a centered scheme of numerical integration are employed to calculate the hydraulic head distribution and the concentration of solute in the flow region. The linear approximation used to define the concentration function within each element is not appropriate for cases involving steep concentration gradients. For such cases, higher order approximations are necessary to assure the continuity of gradients across interelemental boundaries. Numerical examples that illustrate the applicability of the model are presented.     

WATER BUDGETS FOR SMALL DIKED MARSHES1

ABSTRACT: This paper reports an analysis of the water budgets of 10 small (5–6 ha) diked areas (cells) within the Delta Marsh in southcentral Manitoba, Canada. The important terms in the water budget equation in this study were precipitation (P), water pumped in (SWI), evapotranspiration (ET), seepage in (GWI) and out (GWO), and change in storage (ΔS). P, SWI, and S were measured directly, and the sum of ET and GWO determined by difference. Estimating ET as 0.7 pan evaporation gave a seepage loss of 2.9 mm/day from the most intensively studied cell. Other methods of estimating ET produced estimates of GWO ranging from 2.4 to 3.8 mm/day. Water budgets for less intensively studied cells indicated seepage loss increased as perimeter available for seepage increased, but not proportionately. Efforts to measure seepage directly or estimate it from measured hydraulic gradients and hydraulic conductivity produced estimates much lower than the estimates from the water budget equation. Hydraulic conductivities were very heterogeneous, reflecting the sorting of water deposited sediments. Comparison of the hydraulic conductivities with seepage estimates from the water budget strongly suggests water movement downward as well as laterally from these diked areas.     

Upstream Sediment‐Control Dams: Five Decades of Experience in the Rapidly Eroding Dahan River Basin,Taiwan

Sedimentation is emerging as a key issue in sustainable reservoir management. One approach to controlling reservoir sedimentation is to trap sediment in hydraulic structures upstream of the reservoir. In the 1,163‐km² catchment of the Dahan River (Taiwan) over 120 “sabo” dams were built to reduce sediment yield to Shihmen Reservoir. Built in 1963 for water supply, Shihmen has lost over 40% of its 290‐Mm³ storage capacity to sedimentation. Most of these upstream structures were small, but three had capacities >9 Mm³. Field measurements and historical data from the Water Resources Agency show most smaller dams had filled with sediment by 1976. The three largest were full or nearly so by 2007, when one (Barlin Dam) failed, releasing a pulse of 7.5 Mm³, most of its 10.4 Mm³ stored sediment downstream. The Central Range of Taiwan is rapidly eroding (denudation rates 3‐6 mm/yr), so geologically high loads make sediment problems manifest sooner. Even in other environments, however, eventually small dams built upstream of large reservoirs are likely to fill themselves, creating multiple small sediment‐filled reservoirs, some located in sites inaccessible to mechanical removal. Our analysis suggests sabo dams do not offer a long‐term basis for controlling reservoir sedimentation in such a high‐sediment yield environment. Sustainable solutions must somehow pass sediment downstream, as would be accomplished by a sediment bypass around Shihmen Reservoir, as now being studied.     

Can Volunteers Collect Data that are Comparable to Professional Scientists? A Study of Variables Used in Monitoring the Outcomes of Ecosystem Rehabilitation

Having volunteers collect data can be a cost-effective strategy to complement or replace those collected by scientists. The quality of these data is essential where field-collected data are used to monitor progress against predetermined standards because they provide decision makers with confidence that choices they make will not cause more harm than good. The integrity of volunteer-collected data is often doubted. In this study, we made estimates of seven vegetation attributes and a composite measure of six of those seven, to simulate benchmark values. These attributes are routinely recorded as part of rehabilitation projects in Australia and elsewhere in the world. The degree of agreement in data collected by volunteers was compared with those recorded by professional scientists. Combined results showed that scientists collected data that was in closer agreement with benchmarks than those of volunteers, but when data collected by individuals were analyzed, some volunteers collected data that were in similar or closer agreement, than scientists. Both groups’ estimates were in closer agreement for particular attributes than others, suggesting that some attributes are more difficult to estimate than others, or that some are more subjective than others. There are a number of ways in which higher degrees of agreement could be achieved and introducing these will no doubt result in better, more effective programs, to monitor rehabilitation activities. Alternatively, less subjective measures should be sought when developing monitoring protocols. Quality assurance should be part of developing monitoring methods and explicitly budgeted for in project planning to prevent misleading declarations of rehabilitation success.