SIMULATING FLOW IN REGIONAL WETLANDS WITH THE MODFLOW WETLANDS PACKAGE1

ABSTRACT: As part of the Comprehensive Everglades Restoration Plan (CERP), various water supply projects have been proposed in a region located between the Miami metropolitan area and the extensive regional wetland systems that are part of the Everglades or remnant Everglades. A ground water flow model of the surficial aquifer within northern Miami‐Dade County was constructed using MODFLOW to evaluate the effects of these projects on water levels in the wetlands and the underlying surficial aquifer. The new Wetlands package was used to conjunctively simulate overland flow through these wetlands and the shallow ground water system. Comparisons of simulated to measured ground water levels and wetland stages were very satisfactory, where computed and measured water levels agreed within 0.5 ft over most of the period of record at nearly all of the monitoring sites. Temporal trends in water levels were also replicated. It was concluded that the assumptions and methodologies inherent to the Wetlands package were suitable for simulating regional wetland hydrology within the Everglades area.     

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.     

INTERACTIONS BETWEEN GROUND WATER AND SURFACE WATER IN THE SUWANNEE RIVER BASIN,FLORIDA1

ABSTRACT: Ground water and surface water constitute a single dynamic system in most parts of the Suwannee River basin due to the presence of karst features that facilitate the interaction between the surface and subsurface. Low radon-222 concentrations (below background levels) and enriched amounts of oxygen-18 and deuterium in ground water indicate mixing with surface water in parts of the basin. Comparison of surface water and regional ground water flow patterns indicate that boundaries for ground water basins typically do not coincide with surface water drainage subbasins. There are several areas in the basin where ground water flow that originates outside of the Suwannee River basin crosses surface water basin boundaries during both low-flow and high-flow conditions. In a study area adjacent to the Suwannee River that consists predominantly of agricultural land use, 18 wells tapping the Upper Floridan aquifer and 7 springs were sampled three times during 1990 through 1994 for major dissolved inorganic constituents, trace elements, and nutrients. During a period of above normal rainfall that resulted in high river stage and high ground water levels in 1991, the combination of increased amounts of dissolved organic carbon and decreased levels of dissolved oxygen in ground water created conditions favorable for the natural reduction of nitrate by denitrification reactions in the aquifer. As a result, less nitrate was discharged by ground water to the Suwannee River.     

A Stochastic Modeling Approach to Address Hydrogeologic Uncertainties in Modeling Wellhead Protection Boundaries in Karst Aquifers1

Abstract: Development of any numerical ground‐water model is dependent on hydrogeologic data describing the subsurface. These data are obtained from geologic core analyses, stratigraphic analyses, aquifer performance tests, and geophysical studies. But typically in remote areas, these types of data are very sparse and site‐specific in terms of the aerial extent of the resource to be modeled. Uncertainties exist as to how well the available data from a few locations defines a heterogeneous surficial aquifer such as the Biscayne Aquifer in Miami‐Dade County, Florida. This is particularly the case when an exceptionally conductive horizontal flow zone is detected at one site due to specialized testing that was not historically conducted at the other at sites that provided data for the model. Not adequately accounting for the potential effect of the high flow zone in the aquifer within a ground‐water numerical model, even though the zone may be of very limited thickness, might underpredict the well field protection capture boundaries. Applied Stochastic ground‐water modeling in determining well field protection zones is steadily becoming important in addressing the uncertainty of the hydrogeologic subsurface parameters, specifically in karstic heterogeneous aquifers. This is particularly important in addressing the uncertainty of a 60‐day travel time capture zone in the Northwest Well Field, Miami‐Dade County, where a predominantly high flow zone controls much of the flow in the production wells. A stochastic ground‐water modeling application along with combination of pilot points and regularization technique is presented to further consolidate the uncertainty of the subsurface.     

POTENTIAL EFFECTS OF CLIMATE CHANGE ON GROUND WATER IN LANSING,MICHIGAN1

ABSTRACT: Computer simulations involving general circulation models, a hydrologic modeling system, and a ground water flow model indicate potential impacts of selected climate change projections on ground water levels in the Lansing, Michigan, area. General circulation models developed by the Canadian Climate Centre and the Hadley Centre generated meteorology estimates for 1961 through 1990 (as a reference condition) and for the 20 years centered on 2030 (as a changed climate condition). Using these meteorology estimates, the Great Lakes Environmental Research Laboratory's hydrologic modeling system produced corresponding period streamflow simulations. Ground water recharge was estimated from the streamflow simulations and from variables derived from the general circulation models. The U.S. Geological Survey developed a numerical ground water flow model of the Saginaw and glacial aquifers in the Tri‐County region surrounding Lansing, Michigan. Model simulations, using the ground water recharge estimates, indicate changes in ground water levels. Within the Lansing area, simulated ground water levels in the Saginaw aquifer declined under the Canadian predictions and increased under the Hadley.     

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.     

DEVELOPMENT OF A FRACTURED MEDIUM FLOW MODEL FOR A WASTE SITE1

ABSTRACT: A three‐dimensional fractured medium flow model was developed for the Bear Creek Valley (BCV) S‐3 site of the Oak Ridge Reservation (ORR) using SWIFT III. The numerical modeling for this site focused on a conceptual model established through the analysis of heterogeneous geologic units and matrix fracture properties of the subsurface in the BCV area. The SWIFT III modeling analysis was based on the previous modeling studies that used MODFLOW and MODPATH. A rigorous calibration was obtained first by comparing simulated results with the existing data on ground water levels and then by comparing pumping test results with the simulated ground water levels. A satisfactory agreement between observed and simulated results was obtained. The calibrated model was used to determine sustained yield from a ground water interceptor trench. Different withdrawal rates were used to simulate the performance of the trench for the sustained withdrawal of ground water.     

SIMULATING THREE-DIMENSIONAL GROUND WATER RESPONSE IN A SMALL MOUNTAINOUS WATERSHED1

ABSTRACT: Snowmelt from deep mountainous snowpacks is seldom rapid enough to exceed infiltration rates; thus, the source of streamflow in many mountainous watersheds is snowmelt recharge through shallow ground water systems. The hydrologic response and interaction between surface and sub-surface flow processes in these watersheds, which is controlled by basin structure, the spatial distribution of snowmelt, and the hydrogeology of the subsurface, are not well understood. The purpose of this study was to test a three-dimensional ground water model using simulated snowmelt input to simulate ground water response to spatially distributed snowmelt on the Upper Sheep Creek Watershed located within the Reynolds Creek Experimental Watershed in Southwestern Idaho. The model was used to characterize the mountainous aquifer and to delineate the subsurface flow mechanisms. Difficulty in finding a reasonable combination of grid spacing and time stepping within the model was encountered due to convergence problems with the Picard solution to the non-linear variably saturated ground water flow equations. Simulation results indicated that flow may be either unconfined or confined depending on inflow rate and hydrogeologic conditions in the watershed. The flow mechanism had a much faster response time when confined flow occurred. Response to snowmelt from a snow drift approximately 90 m away took only a few hours when flow was confined. Simulated results showed good agreement with piezometer measurements both in magnitude and timing; however, convergence problems with the Picard solution limited applicability of the model.     

Optimal Location and Management of Waste Disposal Facilities Affecting Ground Water Quality1

ABSTRACT: Numerical simulation of ground water solute transport is combined with linear programming to optimize waste disposal. A discretized form of the equation governing solute transport is included as a set of constraints in a linear program. Two problems are described. First, the management model is used to maximize ground water waste disposal. The model constrains disposal activities so that the quality of local ground water supplies is protected. Parametric programming is shown to be important in evaluating waste disposal tradeoffs at the various facilities. Changes in the velocity field induced by waste water injection cause a nonlinearity in the solute transport equation which is dealt with by employing an iterative procedure. The second problem is aimed at identifying all sites which are suitable for waste disposal in the subsurface. The management model is manipulated so that the optimal value of the dual variables are “unit source impact indicators.” This physical interpretation is valuable in identifying feasible disposal sites. The joint simulation and optimization approach permits the management of complex ground water systems where the aquifer is used simultaneously for waste disposal and water supply.     

SUBSURFACE RETURN FLOW AND GROUND WATER RECHARGE OF TERRACE FIELDS IN NORTHERN TAIWAN1

ABSTRACT: This study estimates subsurface return flow and effective ground water recharge in terraced fields in northern Taiwan. Specifically, a three dimensional model, FEMWATER, was applied to simulate percolation and lateral seepage in the terraced fields under various conditions. In the terraced paddy fields, percolation mainly moves vertically downward in the central area, while lateral seepage is mainly focused around the bund. Although the simulated lateral seepage rate through the bund exceeded the percolation rate in the central area of the paddy field, annual subsurface return flow at Pei‐Chi and Shin‐Men was 0.17 × 10⁶ m³ and 0.37 × 10⁶ m³, representing only 0.17 percent and 0.21 percent of the total irrigation water required for rice growth at Pei‐Chi and Shin‐Men, respectively. For upland fields, the effective ground water recharge rate during the second crop period (July to November) exceeded that during the first crop period (January to May) because of the wet season in the second crop period. Terraced paddy fields have the most efficient ground water recharge, with 21.2 to 23.4 percent of irrigation water recharging to ground water, whereas upland fields with a plow layer have the least efficient ground water recharge, with only 4.8 to 6.6 percent of irrigation water recharging to ground water. The simulation results clearly revealed that a substantial amount of irrigation water recharges to ground water in the terraced paddy, while only a small amount of subsurface return flow seeps from the upstream to the downstream terraced paddy. The amounts of subsurface flow and ground water recharge determined in the study are useful for the irrigation water planning and management and provide a scientific basis to reevaluate water resources management in the terrace region under irrigated rice.     

EDWARDS AQUIFER EVALUATION: KINNEY COUNTY,TEXAS1

ABSTRACT: The Edwards Aquifer is one of the most studied and most prolific aquifers in the United States. The aquifer is a heavily fractured and faulted carbonate aquifer with transmissivities in excess of 100 ft²/s. The City of San Antonio relies upon the Edwards Aquifer as its sole source for water. Much work has been done on quantifying recharge to the aquifer and discharge from wells and acquiring aquifer characteristics from pumping tests, specific capacity tests, and geophysical logs. Although the aquifer has been well studied in Bexar County, much less is known about the Edwards Aquifer in Kinney County. This is partly due to the lower population within the county (approximately 3,500 people) relative to the eastern counties (Uvalde, Medina, Bexar, Comal, and Hays) and the great distance of Kinney County from high profile discharge areas such as the City of San Antonio and Comal and San Marcos Springs. Three key products resulted from this study: (1) exploratory well drilling and the largest aquifer test in the county that were conducted to evaluate the well yields within a 10,000 acre study area in which a drawdown of 2.5 ft approximately 1.2 miles away was observed while pumping at approximately 4,600 gpm; (2) a recharge estimate for the Edwards Aquifer within Kinney County of approximately 71,382 ac‐ft/yr; and (3) locating the Brackettville Groundwater Divide from an evaluation of ground water flow direction and hydrograph analysis. These results help evaluate the complex hydraulics occurring within Kinney County and aid in development of ground water modeling that will be used in managing the Edwards Aquifer.     

SHALLOW AQUIFERS RELATIVE TO SURFACE WATER,LOWER NORTH PLATTE RIVER VALLEY,WYOMING1

ABSTRACT. The occurrence of ground water in the lower North Platte Valley, Goshen County, Wyoming, was studied to determine safe yield within the alluvial aquifer under varying discharge and recharge conditions. The alluvium of the North Platte is extensively developed for irrigation purposes and the effects of large-scale pumpage are of major concern. Actual withdrawals are estimated to be 46,000 acre-feet. Should pumping reach potentially higher levels an overdraft is expected. Effect of ground water withdrawals are established from projections of the flow regime within the alluvial aquifer. A time dependent, numerical model was employed to predict aquifer response to increased withdrawals. The results suggest that more efficient use of surface waters and/or increased use of ground water will reduce the annual subsurface return flow to the North Platte River and its tributaries by an amount equal to the reduced ground water recharge increment. Alternatives are available for management of the lower North Platte alluvial aquifer. The preferred course is to correlate surface and subsurface water rights, in light of convenience, economics, and best means of storage for maximum utilization of the single water resource.     

SIMULATING THE EFFECTS OF REDUCED PRECIPITATION ON GROUND WATER AND STREAMFLOW IN THE NEBRASKA SAND HILLS1

ABSTRACT: The Nebraska Sand Hills have a unique hydrologic system with very little runoff and thick aquifers that constantly supply water to rivers, lakes, and wetlands. A ground water flow model was developed to determine the interactions between ground water and streamflow and to simulate the changes in ground water systems by reduced precipitation. The numerical modeling method includes a water balance model for the vadose zone and MOD‐FLOW for the saturated zone. The modeling results indicated that, between 1979 and 1990, 13 percent of the annual precipitation recharged to the aquifer and annual ground water loss by evapotranspiration (ET) was only about one‐fourth of this recharge. Ground water discharge to rivers accounts for about 96 percent of the streamflow in the Dismal and Middle Loup rivers. When precipitation decreased by half the average amount of the 1979 to 1990 period, the average decline of water table over the study area was 0.89 m, and the streamflow was about 87 percent of the present rate. This decline of the water table results in significant reductions in ET directly from ground water and so a significant portion of the streamflow is maintained by capture of the salvaged ET.     

STABLE ISOTOPE CHARACTERIZATION OF THE IMPACTS OF ARTIFICIAL GROUND WATER RECHARGE1

ABSTRACT: Stable isotopes of deuterium and oxygen-18 of surface and ground water, together with anion concentrations and hydraulic gradients, were used to interpret mixing and flow in ground water impacted by artificial recharge. The surface water fraction (SWF), the percentage of surface water in the aquifer impacted via recharge, was estimated at different locations and depths using measured deuterium/hydrogen (DIH) ratios during the 1992, 1993, and 1994 recharge seasons. Recharged surface water completely displaced the ground water beneath the recharge basins from the regional water table at 7.60 m to 12.16 m below the land surface. Mixing occurred beneath the recharge structures in the lower portions of the aquifer (>12.16 m). Approximately 12 m down-gradient from the recharge basin, the deeper zone (19.15 m depth) of the primary aquifer was displaced completely by recharged surface water within 193, 45, and 55 days in 1992, 1993, and 1994, respectively. At the end of the third recharge season, recharged surface water represented ～50 percent of the water in the deeper zone of the primary aquifer ～1000 m downgradient from the recharge basin. A classic asymmetrical distribution of recharged surface water resulted from the recharge induced horizontal and vertical hydraulic gradients. The distribution and breakthrough times of recharged surface water obtained with stable isotopes concurred with those of major anions and bromide in a tracer test conducted during the 1995 recharge season. This stable isotope procedure effectively quantified mixing between surface and ground water.     

THE ECONOMIC VALUE OF GROUND WATER RECHARGE FOR IRRIGATION USE1

ABSTRACT: Artificial recharge as a means of augmenting water sup plies for irrigation is a management alternative which policy makers in ground water decline areas are beginning to consider seriously. A conceptual model is developed to evaluate the economic benefits from ground water recharge under conditions where the major water use is irrigation. The methodology presented separates recharge benefits into two components: pumping cost savings and aquifer extension benefits. This model is then applied to a Nebraska case to approximate the value of recharge benefits as a function of aquifer response. discount rate, and commodity prices. It was found that recharge benefits vary from less than $2 to over $6 an acre foot recharged.     

REGIONAL SCALE MODELING OF SURFACE AND GROUND WATER INTERACTION IN THE SNAKE RIVER BASIN1

ABSTRACT: Changes in irrigation and land use may impact discharge of the Snake River Plain aquifer, which is a major contributor to flow of the Snake River in southern Idaho. The Snake River Basin planning and management model (SRBM) has been expanded to include the spatial distribution and temporal attenuation that occurs as aquifer stresses propagate through the aquifer to the river. The SRBM is a network flow model in which aquifer characteristics have been introduced through a matrix of response functions. The response functions were determined by independently simulating the effect of a unit stress in each cell of a finite difference groundwater flow model on six reaches of the Snake River. Cells were aggregated into 20 aquifer zones and average response functions for each river reach were included in the SRBM. This approach links many of the capabilities of surface and ground water flow models. Evaluation of an artificial recharge scenario approximately reproduced estimates made by direct simulation in a ground water flow model. The example demonstrated that the method can produce reasonable results but interpretation of the results can be biased if the simulation period is not of adequate duration.     

GROUND WATER REPLENISHMENT FROM RIVERFLOW1

ABSTRACT: Ground water contours and gradients were determined within and adjacent to a gravelly riverbed with about 100 m active channel width and about 800 m flood plain width during low flow, making use of about 60 shallow dug wells and one deeper cased well Ground water slopes away from the riverbed to the south and slopes to the riverbed from the north. Replenishment of ground water to the south is past a “shell,” a zone of low permeability at the margin of and below the riverbed, (1) into a shallow aquifer from “sideway-spill-over” and (2) into an underlying aquifer that is seemingly confined by the “shell.” Seepage from river channels to underflow takes place all along channel lengths except below riffles where the flow direction is reversed. Half the seepage to underflow was lost to ground water replenishment Underflow seems to constitute merely a few percent of total riverflow but makes a disproportionate contribution to ground water replenishment.     

GROUND WATER SIMULATION USING A ONE DIMENSIONAL FLOW MODEL1

ABSTRACT: In projects involving ground water problems, dependence on the mathematical modeling of the ground water flow phenomena is inescapable. At present, two dimensional flow models, which require tremendous amounts of computer time and storage, are generally used. When such bulky models are used for planning purposes, the two requirements (computer time and storage) can severely limit the number of alternatives that can be considered. A simple quantity and quality simulation model is developed here which requires considerably less computer time and storage and gives reasonably accurate results. The model was applied to simulate a ground water basin in San Luis Rey River in Southern California. The results were compared with those obtained by a USGS model. It was found that the simple model gave results which were consistentaly within five percent of the USGS model results, while the requirements on computer time and storage were drastically reduced.     

EFFECT OF SEMI‐PERCHED GROUND WATER ON MONITORING OF GROUND WATER LEVELS IN A DEVELOPED AQUIFER1

ABSTRACT: Interpretation of ground water level changes in a developed aquifer usually relies on reference to some benchmark such as “predevelopment” ground water levels, changes from fall to fall and/or spring to spring, or to determination of maximum stress during the pumping season. The assumption is that ground water levels measured in the monitoring well accurately reflect the state of the ground water resource in terms of quantity in storage and the effects of local pumping. This assumption is questionable based on the patterns shown in continuous hydrographs of water levels in monitoring wells in Nebraska, and wells installed to determine vertical gradients. These hydrographs show clear evidence for vertical ground water gradients and recharge from overlying parts of the aquifer system to deeper zones in which production wells are screened. The classical concept of semi‐perched ground water, as described by Meinzer, is demonstrated by these hydrographs. The presence of semi‐perched ground water (Meinzer definition, there is no intervening unsaturated zone) invalidates the use of measured ground water levels in regional observation programs for detailed numerical management of the resource. Failure to recognize the Meinzer effect has led to faulty management. The best use of data from the observation well network would be for detection of trends and education unless it is clearly understood what is being measured.     

OPTIMAL GROUND WATER MANAGEMENT IN TWO-AQUIFER SYSTEMS1

ABSTRACT; This paper presents a numerical model for the prediction of optimal ground water withdrawal from a two-aquifer system by observing a set of constraints determined by the ecological conditions of the ground water basin. The aquifer system consists of an upper unconfined and a lower confined aquifer with a leaky stratum between them. It is assumed that water is withdrawn from the confined aquifer only, but the unconfined aquifer will also be affected due to the leakiness of the layer separating the upper and lower aquifers. Simulation and linear programming are employed for developing a computer model for the optimal management of such systems, with the objectives of determining withdrawal rates for predetermined ground water levels.