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.     

GROUND WATER MANAGEMENT OPTIMIZATION USING GENETIC ALGORITHMS AND SIMULATED ANNEALING: FORMULATION AND COMPARISON1

ABSTRACT: Genetic algorithms (GA) and simulated annealing (SA), two global search techniques, are coupled with MODFLOW, a commonly used groundwater flow simulation code, for optimal management of ground water resources under general conditions. The coupled simulation-optimization models allow for multiple management periods in which optimal pumping rates vary with time to reflect the changing flow conditions. The objective functions of the management models are of a very general nature, incorporating multiple cost terms such as the drilling cost, the installation cost, and the pumping cost. The models are first applied to two-dimensional maximum yield and minimum cost water supply problems with a single management period, and then to a multiple management period problem. The strengths and limitations of the GA and SA based models are evaluated by comparing the results with those obtained using linear programming, nonlinear programming, and differential dynamic programming. For the three example problems examined in this study, the GA and SA based models yield nearly identical or better solutions than the various programming methods. While SA tends to outperform GA in terms of the number of forward simulations needed, it uses more empirical control parameters which have significant impact on solution efficiency but are difficult to determine.     

MODELING GROUND WATER FLOW USING FLUX BOUNDARY CONDITIONS1

ABSTRACT: Determination of the boundary conditions for modeling ground water flow is a critical point especially in regional models. Normally the regional models require model areas that are greater than the given area of interest. This work focuses on the prediction of hydraulic heads in regional models using flux boundary conditions. The model uses flux boundary conditions that were estimated using a radial flow analog and Darcy's law. The regional model that is presented uses no parameter identification (inverse estimation) procedures. In the present work, the Houston area was used. The simulation of the hydrological conditions of the Chicot and Evangeline Aquifers that underlie the Houston area were made using the available information about the geological profile in the Houston region and the current information about the existing production wells. The regional model works as a forward problem. The system parameters such as hydraulic conductivity, specific storage, and hydrological stresses were specified, and the model predicts the hydraulic head. Actual data from piezometers operated by the U.S. Geological Survey (USGS) in many places throughout Houston were used as initial conditions. Some piezometric head data were generated using the regional variable theory called kriging to supply head estimates in areas where data were unavailable. The Modular Three Dimensional Finite Difference Groundwater Flow Model developed by the USGS was used to predict the hydraulic heads. The predicted ground water heads are compared to the actual data. The results show that the model performs well for locations where data were available.     

USING UNCENSORED TRACE-LEVEL MEASUREMENTS TO DETECT TRENDS IN GROUND WATER CONTAMINATION1

ABSTRACT: Environmental decision making involving trace-levels of contaminants can be complicated by censoring, the practice of reporting concentrations either as less than the limit of detection (LOD) or as not detected (ND) when a test result is less than the LOD. Censoring can result in data series that are difficult to meaningfully summarize, graph, and analyze through traditional statistical methods. In spite of the relatively large measurement errors associated with test results below the LOD, simple and meaningful analyses can be carried out that provide valuable information not available if data are censored. For example, an indication of increasing levels of contamination at the fringe of a plume can act as an early warning signal to trigger further study, an increased sampling frequency, or a higher level of remediation at the source. This paper involves the application of nonparametric trend analyses to uncensored trace-level groundwater monitoring data collected between March 1991 and August 1994 on dissolved arsenic and chromium for seven wells at an industrial site in New York.     

A DIGITAL MODEL APPLIED TO GROUND WATER RECHARGE AND MANAGEMENT1

ABSTRACT: Under Colorado's appropriative water right system, withdrawals by junior ground water rights must be curtailed to protect senior surface water appropriators sharing the same river system unless the ground water users replace the amount of their injury to the river under an approved plan for augmentation. Compensation of such injury with surface water may not only be expensive but unreliable in dry years. As an alternative, the curtailment of pumping may be obviated by recharging unused surface water into the aquifer when available and withdrawing it when needed. In order to manage such an operation, a practical tool is required to accurately determine that portion of the recharge water that does not return to the river before pumping for irrigation. A digital model was used for this purpose in a demonstration recharge project located in the South Platte River basin in northeastern Colorado. This paper summarizes the experiences gained from this project, the results of the digital model, the economic value of recharge, and the feasibility of the operation. It was determined through the use of the digital model that, with the given conditions in the area, 77 percent of the recharged water remained available for pumping. Economic analyses showed that water could be recharged inexpensively averaging about two dollars per acre foot.     

GROUND WATER MONITORING NETWORK DESIGN USING MULTIPLE CRITERIA DECISION MAKING ANT) GEOSTATISTICS1

ABSTRACT: A multi-criteria approach to ground water quality monitoring network design is developed. The methodology combines multi-criteria decision making (MCDM) and modifications of geostatistical variance reduction analysis. Composite programming, a distance based optimization algorithm that employs a hierarchial structure, is used for the MCDM component of the design methodology. MCDM allows the consideration of numerous, often conflicting, design criteria. The methodology is useful for identifying the preferred combination of direct borehole and indirect geo-electric data. It also permits the use of prior information during initial stages of network development. Multi-variate kriging is employed to evaluate network performance using the combination of direct borehole data and indirect geoelectric data. Weighted measures of estimation variance are used as primary measures of performance, with the reduction in estimation variance being computed by the fictitious point method. Case study results demonstrate that the network design methodology can be used in both early and late phases of network development. It also leads to selection of the preferred combination and spatial orientation of direct and indirect data sources while considering cost-effectiveness and performance of alternative designs.     

UNSATURATED FLOW THROUGH SOLID WASTE LANDFILLS: MODEL AND SENSITIVITY ANALYSIS1

ABSTRACT: The movement of precipitation water infiltrating through the material (refuse) of solid waste landfills is examined via numerical solution of the equations of continuity, and motion (Darcy's Law). The solution of the equations is obtained by a fully implicit, finite-difference scheme. Both unsaturated and saturated surface conditions are considered, making the scheme suitable for real-time simulation of net precipitation and moisture redistribution events. A sensitivity analysis showed that for unsaturated surface conditions the solution is primarily affected by hydraulic conductivity and capillary diffusivity, and is relatively independent of the space and time steps. In addition, the precipitation averaging process is shown to be critical in the correct computation of moisture transport during the time period where the transition from unsaturated to saturated conditions occurs. The model presented herein is suitable for analysis of water movement through landfills, and the design of bottom collection systems.     

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.     

SIMULATION OF REGIONAL GROUND WATER FLOW IN BEDROCK,SOUTHWESTERN NEW YORK-NORTHWESTERN PENNSYLVANIA1

ABSTRACT: This paper presents the results of steady-state three-dimensional computer simulations to determine the hydrogeologic setting of formation water in the hydrocarbon producing formations of southwestern New York and northwestern Pennsylvania. Recharge areas for the regional ground water flow systems in the study area are the Valley Heads Moraine and Allegheny uplands; discharge areas are Lakes Erie and Ontario to the north and the northern margin of the Appalachian basin to the south. Simulated ground water flow in all model layers moves north from the ground water divide on the Valley Heads Moraine towards Lake Erie at a rate from 10^?-⁶ to 10^?-³ ft/day. South of the divide intermediate-scale and local-scale flow systems occur in the upper 4000 feet of the stratigraphic section and the directions of ground water flow diverge towards major rivers and other topographically low areas.     

USING COMPUTED STREAM FLOW IN WATERSHED ANALYSIS1

Modifications in the computed climatic water budget have made it possible to achieve good agreement between computed and measured stream flow on both a monthly and annual basis in basins without appreciable winter snow cover. Comparisons of computed and measured stream flow in 28 basins on the Delmarva peninsula show that for some basins the agreement is excellent (regression line essentially equals unity), for other basins the regression line has a slope of one but it is displaced above or below the y=x line, while for other basins, the slope of the regression line differs appreciably from unity. Study of the basins where agreement between computed and measured values is only fair to poor reveals that the patterns of disagreement can be used to provide information on the water holding capacity in the root zone of the soil, on the quantity of deep aquifer recharge within the basin, or on the effect of human modifications within the basin. The technique should also reveal the quantity of interbasin transfers or other consumptive uses within the basin. The water budget, thus, becomes a useful tool to study hydrologic characteristics or their changes over time within a basin.     

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.     

SIMULATION OF INTEGRATED SURFACE WATER AND GROUND WATER SYSTEMS - MODEL FORMULATION1

ABSTRACT: The unique characteristics of the hydrogeologic system of south Florida (flat topography, sandy soils, high water table, and highly developed canal system) cause significant interactions between ground water and surface water systems. Interaction processes involve infiltration, evapotranspiration (ET), runoff, and exchange of flow (seepage) between streams and aquifers. These interaction processes cannot be accurately simulated by either a surface water model or a ground water model alone because surface water models generally oversimplify ground water movement and ground water models generally oversimplify surface water movement. Estimates of the many components of flow between surface water and ground water (such as recharge and ET) made by the two types of models are often inconsistent. The inconsistencies are the result of differences in the calibration components and the model structures, and can affect the confidence level of the model application. In order to improve model results, a framework for developing a model which integrates a surface water model and a ground water model is presented. Dade County, Florida, is used as an example in developing the concepts of the integrated model. The conceptual model is based on the need to evaluate water supply management options involving the conjunctive use of surface water and groundwater, as well as the evaluation of the impacts of proposed wellfields. The mathematical structure of the integrated model is based on the South Florida Water Management Model (SFWMM) (MacVicar et al., 1984) and A Modular Three-Dimensional Finite-Difference Groundwater Flow Model (MODFLOW) (McDonald and Harbaugh, 1988).     

PHYTOPLANKTON RESPONSES TO GROUND WATER ADDITION IN CENTRAL FLORIDA LAKES1

ABSTRACT: Phytoplankton were studied in lakes augmented with water from the Floridan aquifer and in control lakes. Augmented lakes exhibited enhanced phytoplankton diversity which was believed to be related to chemical changes in the lakes brought about as a result of ground water addition. In particular, elevated concentrations of inorganic carbon appeared to influence phytoplankton populations.     

GROUND WATER QUALITY PREDICTION USING CLIMATIC INDICES1

ABSTRACT: Drought affects the quality of ground water in certain aquifers used by municipalities in Kansas. Water quality changes occur as a function of the amount of water available for recharge and hence to dilute more mineralized ground waters. Several measures of meteorological drought, including the Palmer Index and Eagleman Aridity Index, were correlated with water quality data to determine the degree of association. Several locations showed sharp delinces in water quality as the drought progressed. These relationships can be used to predict possible variations in present and future well-water supplies in locations subject to drought induced water quality deterioration.     

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.     

GROUND WATER FLOW AND RUNOFF IN A COASTAL PLAIN STREAM1

ABSTRACT: The quantity, seasonality, and sources of flow were analyzed for two segments of Four Mile Branch, a small stream on the Coastal Plain of South Carolina using data obtained from USGS gauging stations. Flows in the “upstream segment,” a 12.6-km² watershed comprising the head waters of Four Mile Branch, averaged 0.129 m³ s^?1 and showed a distinctly seasonal pattern, with maximum flows in February and March and minimum flows in September and October. Inflow to the “downstream segment,” a 2.2-km² watershed associated with the main channel, averaged 0.059 m³ s^?1 and showed no seasonal patterns. Discharges per unit area of watershed were greater for the downstream segment, 0.83 m³ per year per m² of land surface, than for the upstream segment, 0.32 m³ per year per m². The differences in discharge rates and seasonalities between the two segments reflect differences in aquifers supplying the different segments. Analyses of Streamflow by hydrograph separation and Streamflow partitioning methods indicated that greater than 90 percent of the flows in the upstream and downstream segments were due to ground water-driven base flows.     

CHARACTERIZING GROUND WATER FLOW IN THE MUNICIPAL WELL FIELDS OF CEDAR RAPIDS,IOWA, WITH SELECTED ENVIRONMENTAL TRACERS1

ABSTRACT: Cedar Rapids obtains its municipal water supply from a shallow alluvial aquifer along the Cedar River in east-central Iowa. Water samples were collected and analyzed for selected isotopes and chlorofluorocarbons to characterize the ground-water flow system near the municipal well fields. Analyses of deuterium and oxygen-18 indicate that water in the alluvial aquifer and in the underlying carbonate bedrock aquifer was recharged from precipitation during modern climatic conditions. Analyses of tritium indicate modern, post-1952, water in the alluvial aquifer and older, pre-1952, water in the bedrock aquifer. Mixing of the modern and older waters occurs in areas where (1) the confining layer between the two aquifers is discontinuous, (2) the bedrock aquifer is fractured, or (3) pumping of supply wells induces the flow of water between aquifers. Analyses of chlorofluorocarbons were used to determine the date of recharge of water samples. Water in the bedrock aquifer likely was recharged prior to the 1950s. Water in the alluvial aquifer likely was recharged from the 1960s to 1990s. Biodegradation or sorption probably affected some of the ground water analyzed for chlorofluorocarbons. These processes reduce the concentrations of CFCs, which results in older than actual calculated dates of recharge.     

MODEL ACCURACY IN SNOWMELT-RUNOFF FORECASTS EXTENDING FROM 1 TO 20 DAYS1

ABSTRACT: This paper examines the performance of snowmelt-runoff models in conditions approximating real-time forecast situations. These tests are one part of an intercomparison of models recently conducted by the World Meteorological Organization (WMO). Daily runoff from the Canadian snowmelt basin Illecille. waet (1155 km², 509–3150 m a.s.l.) was forecast for 1 to 20 days ahead. The performance of models was better than in a previous WMO project, which dealt with runoff simulations from historical data, for the following reasons: (1) conditions for models were more favorable than a real-time forecast situation because measured input data and not meteorological forecast inputs were distributed to the modelers; (2) the selected test basin was relatively easy to handle and familiar from the previous WMO project; and (3) all kinds of updating were allowed so that some models even improved their accuracy towards longer forecast times. Based on this experience, a more realistic follow-up project can be imagined which would include temperature forecasts and quantitative precipitation forecasts instead of measured data.     

A CRITICAL APPROACH TO THE CALIBRATION OF A WATERSHED MODEL1

ABSTRACT: A complex watershed-scale water quality simulation model, the Hydrological Simulation Program-FORTRAN (HSPF) model, was calibrated for a 16 km² catchment. The simulation step size was 0.33 hours with predicted and recorded hydrologic flows compared on an annual and monthly basis during a total calibration period of four years. Unguided numerical optimization when applied alone did not yield a model parameter set with acceptable predictive capability; instead, it was necessary to apply a critical process that included sensitivity analysis, numerical optimization, and testing of derived model parameter sets to evaluate their performance for periods other than those for which they were determined. Using this critical calibration process, the model was proven to have significant predictive capability. Numerical optimization is an aid for model calibration, but it must not be used blindly.     

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.