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.     

A COMPARISON OF MODELS PREDICTING GROUND WATER LEVELS ON HILLSIDE SLOPES1

ABSTRACT: A comparative study of ground water level predictions on hillside slopes using two models is presented. The models are a simplified mass balance model that has components for evapotran-spiration, recharge, and drainage; and a two-dimensional finite difference model that employs kriging to estimate soil parameters and accounts for non-uniform thickness of the soil layer. These models are representative of a wide range of modeling capabilities and are used to illustrate the sensitivity of ground water level predictions to the sophistication of the modeling techniques. The drainage and recharge components of the two models are evaluated and the importance of unsaturated flow in recharge computations is underscored. Piezometric observations in a small drainage depression on the slope of Kennel Creek Valley in Tongass National Forest, Alaska, were used to evaluate the two models. The results show that, although the predictions differ from the field observations, the simple physically-based mass balance model predicts the ground water levels as well as the two-dimensional model. It is suggested that caution should be exercised in using complex models to validate simpler models.     

AUTOMATED WEB GIS BASED HYDROGRAPH ANALYSIS TOOL,WHAT1

The separation of the base flow component from a varying streamflow hydrograph is called “hydrograph analysis.” In this study, two digital filter based separation modules, the BFLOW and Eckhardt filters, were incorporated into the Web based Hydrograph Analysis Tool (WHAT) system. A statistical component was also developed to provide fundamental information for flow frequency analysis and time series analysis. The Web Geographic Information System (GIS) version of the WHAT system accesses and uses U.S. Geological Survey (USGS) daily streamflow data from the USGS web server. The results from the Eckhardt filter method were compared with the results from the BFLOW filter method that was previously validated, since measured base flow data were not available for this study. Following validation, the two digital filter methods in the WHAT system were run for 50 Indiana gaging stations. The Nash‐Sutcliffe coefficient values comparing the results of the two digital filter methods were over 0.91 for all 50 gaging stations, suggesting the filtered base flow using the Eckhardt filter method will typically match measured base flow. Manual separation of base flow from streamflow can lead to inconsistency in the results, while the WHAT system provides consistent results in less than a minute. Although base flow separation algorithms in the WHAT system cannot consider reservoir release and snowmelt that can affect stream hydrographs, the Web based WHAT system provides an efficient tool for hydrologic model calibration and validation. The base flow information from the WHAT system can also play an important role for sustainable ground water and surface water exploitation, including irrigation and industrial uses, and estimation of pollutant loading from both base flow and direct runoff. Thus, best management practices can be appropriately applied to reduce and intercept pollutant leaching if base flow contributes significant amounts of pollutants to the stream. This Web GIS based system also demonstrates how remote, distributed resources can be shared through the Internet using Web programming.     

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).     

ESTIMATING TRANSMISSIVITY AND HYDRAULIC CONDUCTIVITY OF CHICOT AQUIFER FROM SPECIFIC CAPACITY DATA1

ABSTRACT: Specific capacity data obtained from Well Construction reports which are available from USGS offices, can provide useful estimates of tranamissivity (T), and hydraulic conductivity (K), of an aquifer. The Chicot Aquifer in Louisiana is one of the largest sources of fresh ground water in North America. Hydrologic data collected for the Chicot Aquifer indicate that specific capacity tests can be used in estimating local and regional values for T and K, if the Cooper-Jacob equation for transient flow is used with proper corrections for well loss and partial penetration. Where full scale pumping test data are scarce, specific capacity test data that are adequately distributed spatially can be used to map changes in T and K values and can be summarized statistically to indicate applicable regional values. A computer program called “TGUESS” which is available from International Ground Water Modeling Center, Holcomb Research Institute, was used in this study. The contour maps for T and K values are prepared for different well depth intervals to avoid wide variation of values.     

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.     

PARAMETER ESTIMATES FOR A GROUND WATER MODEL1

ABSTRACT: This paper presents a parameter sensitivity study of a two-dimensional flow and transport model of a contaminated site. Hydrogeological and site data from previous investigations were used for calibration. The USGS contaminant transport model (MOC) was used. After flow calibration to establish a reference model, parameters were varied to examine the effect each had on predictions of a contaminant plume. Hydrogeological parameters and a step size parameter were incrementally varied individually. Each result was compared to the reference model output to evaluate changes in concentration values and contaminant plume configuration. The study indicated that a generally predictable trend can be established for some parameters not affected by pumping or similar high stresses. Ranges were identified to relate concentration error or plume change to the amount of parameter error. Some parameter perturbations produced distorted model responses at high stress locations. Porosity and anisotropy were found to be the most influential of the model parameters studied on the plume predictions. (KEY TERMS: ground water hydrology; hydrogeology; pollution modeling; water quality; model calibration.)     

STORMFLOW SIMULATION USING A GEOGRAPHICAL INFORMATION SYSTEM WITH A DISTRIBUTED APPROACH1

ABSTRACT: With the increasing availability of digital and remotely sensed data such as land use, soil texture, and digital elevation models (DEMs), geographic information systems (GIS) have become an indispensable tool in preprocessing data sets for watershed hydrologic modeling and post processing simulation results. However, model inputs and outputs must be transferred between the model and the GIS. These transfers can be greatly simplified by incorporating the model itself into the GIS environment. To this end, a simple hydrologic model, which incorporates the curve number method of rainfall‐runoff partitioning, the ground‐water base‐flow routine, and the Muskingum flow routing procedure, was implemented on the GIS. The model interfaces directly with stream network, flow direction, and watershed boundary data generated using standard GIS terrain analysis tools; and while the model is running, various data layers may be viewed at each time step using the full display capabilities. The terrain analysis tools were first used to delineate the drainage basins and stream networks for the Susquehanna River. Then the model was used to simulate the hydrologic response of the Upper West Branch of the Susquehanna to two different storms. The simulated streamflow hydrographs compare well with the observed hydrographs at the basin outlet.     

MODELING METALS TRANSPORT AND SEDIMENT/WATER INTERACTIONS IN A MINING IMPACTED MOUNTAIN STREAM1

ABSTRACT: The U.S. Environmental Protection Agency (USEPA) Water Quality Analysis Simulation Program (WASP5) was used to model the transport and sediment/water interactions of metals under low flow, steady state conditions in Tenmile Creek, a mountain stream supplying drinking water to the City of Helena, Montana, impacted by numerous abandoned hard rock mines. The model was calibrated for base flow using data collected by USEPA and validated using data from the U.S. Geological Survey (USGS) for higher flows. It was used to assess metals loadings and losses, exceedances of Montana State water quality standards, metals interactions in stream water and bed sediment, uncertainty in fate and transport processes and model parameters, and effectiveness of remedial alternatives that include leaving contaminated sediment in the stream. Results indicated that during base flow, adits and point sources contribute significant metals loadings to the stream, but that shallow ground water and bed sediment also contribute metals in some key locations. Losses from the water column occur in some areas, primarily due to adsorption and precipitation onto bed sediments. Some uncertainty exists in the metal partition coefficients associated with sediment, significance of precipitation reactions, and in the specific locations of unidentified sources and losses of metals. Standards exceedances are widespread throughout the stream, but the model showed that remediation of point sources and mine waste near water courses can help improve water quality. Model results also indicate, however, that alteration of the water supply scheme and increasing base flow will probably be required to meet all water quality standards.     

EVAPOTRANSPIRATION CONCEPTUALIZATION IN THE HSPF‐MODFLOW INTEGRATED MODELS1

In 1988, the Florida Institute of Phosphate Research (FIPR) funded project to develop an advanced hydrologic model for shallow water table systems. The FIPR hydrologic model (FHM) was developed to provide an improved predictive capability of the interactions of surface water and ground water using its component models, HSPF and MODFLOW. The Integrated Surface and Ground Water (ISGW) model was developed from an early version of FHM and the two models were developed relatively independently in the late 1990s. Hydrologic processes including precipitation, interception, evapotranspiration, runoff, recharge, streamflow, and base flow are explicitly accounted for in both models. Considerable review of FHM and ISGW and their applications occurred through a series of projects. One model evolved, known as the Integrated Hydrological Model IHM. This model more appropriately describes hydrologic processes, including evapotranspiration fluxes within small distributed land‐based discretization. There is a significant departure of many IHM algorithms from FHM and ISGW, especially for soil water and evapotranspiration (ET). In this paper, the ET concepts in FHM, ISGW, and IHM will be presented. The paper also identifies the advantages and data costs of the improved methods. In FHM and IHM, ground water ET algorithms of the MODFLOW ET package replace those of HSPF (ISGW used a different model for ground water ET). However, IHM builds on an improved understanding and characterization of ET partitioning between surface storages, vadose zone storage, and saturated ground water storage. The IHM considers evaporative flux from surface sources, proximity of the water table to land surface, relative moisture condition of the unsaturated zone, thickness of the capillary zone, thickness of the root zone, and relative plant cover density. The improvements provide a smooth transition to satisfy ET demand between the vadose zone and deeper saturated ground water. While the IHM approach provides a more sound representation of the actual soil profile than FHM, and has shown promise at reproducing soil moisture and water table fluctuations as well as field measured ET rates, more rigorous testing is necessary to understand the robustness and/or limitations of this methodology.     

A LAND-USE POLICY BASED ON WATER SUPPLY1

ABSTRACT: Santa Fe Country, New Mexico, has adopted a land-use policy in which zoning densities provide a balance between the water use on a parcel of land and the water supply available bencath that land. In two of four mapped hydrologic zones, ground water in storage will be allowed to be mined to exhaustion in 100 years (40 years in urban areas). Elsewhere, the policy is for a steady state with use balanced by recharge. Equations to determine storage or recharge can be solved using site specific data or regional estimates of hydrologic conditions. Substantial reductions in the lot size requirements are allowed if water conservation convenants are adopted. Public acceptance indicates that the policy successfully integrates technical and political concerns. It is simple to administer, yet reflects widely expressed public goals and values.     

WATER LOSSES ALONG A REACH OF THE PECOS RIVER IN NEW MEXICO1

ABSTRACT: A reach of the Pecos River, located in eastern New Mexico, was examined to evaluate losses of river flows due to evaporation, seepage, and transpiration. An accurate assessment of the water losses along this reach is critical for determining how water rights are adjudicated for water users in the Pecos basin and interstate compact accounting. Water losses significantly impact flows through critical habitat for species protected under the Endangered Species Act. Daily losses of river flows were analyzed for the study reach that extends from immediately below the Pecos River confluence with Taiban Creek to the United States Geological Survey (USGS) gage near Acme. The analysis was completed with consideration for other processes including flood wave travel times and attenuation along with stream bank storage and returns. The analysis was completed using daily stream flow data from USGS gages located along the study reach. Empirical seasonal functions were developed to relate flow loss to the flow rate in the river. The functions were ultimately developed to provide a method for comparing the effects of different river flows on the available water supply.     

A PARAMETRIC MODEL CALIBRATED WITH A PHYSICALLY BASED MODEL FOR RUNOFF PREDICTION FROM UNGAGED STREAMS1

ABSTRACT: Recent developments in the numerical solution of the governing partial differential equations for overland and channel flow should make possible physically based models which predict runoff from ungaged streams. However, these models, which represent the watershed by sets of intersecting planes, are complex and require much computer time. Parametric models exist that have the advantage of being relatively simple, and once calibrated are inexpensive to use and require limited data input. In this study, a procedure was developed for calibrating a parametric model against a physically based model, utilizing base areas of one acre and one square mile, with the expectation that base areas can be combined to model real watersheds. Simulation experiments with the physically based model showed that, for the one-acre base area, the dominant parameter (cell storage ratio, K) related to the slope and friction of the planes, whereas for one square-mile areas, the dominant parameters (K plus a lag factor, L) relate to channel properties. These parameters decreased exponentially as rainfall intensity increased.     

MODELING THE FOREST HYDROLOGY OF WETLAND-UPLAND ECOSYSTEMS IN FLORIDA1

ABSTRACT: Few hydrological models are applicable to pine flat-woods which are a mosaic of pine plantations and cypress swamps. Unique features of this system include ephemeral sheet flow, shallow dynamic ground water table, high rainfall and evapotranspiration, and high infiltration rates. A FLATWOODS model has been developed specifically for the cypress wetland-pine upland landscape by integrating a 2-D ground water model, a Variable-Source-Area (VAS)-based surface flow model, an evapotranspiration (ET) model, and an unsaturated water flow model. The FLATWOODS model utilizes a distributed approach by dividing the entire simulation domain into regular cells. It has the capability to continuously simulate the daily values of ground water table depth, ET, and soil moisture content distributions in a watershed. The model has been calibrated and validated with a 15-year runoff and a four-year ground water table data set from two different pine flat woods research watersheds in northern Florida. This model may be used for predicting hydrologic impacts of different forest management practices in the coastal regions.     

ON THE CAPACITY-EXPANSION OF DESALINATION PLANTS AND STORAGE TANKS1

ABSTRACT A mathematical programming model is structured to find the optimal time and capacity expansion path of desalination plants and storage tanks for a hypothetical community which depends on desalination as its sole, or major, water supply source. The objective is to determine the least costly combination of sues and times of installation (of both desalting plants and storage tanks) which can meet a rising water demand over a finite planning horizon. The optimality criterion used in the model is based on two major economic elements: the economies of scale inherent in such facilities and the time-value of money represented by the interest rate, the former favoring large capacities and the latter small capacities. The model is applied using three population growth patterns and two interest rates. The water demand components for every year in the planning period are computed using empirical formulas which are based on population and other basic data. The model is solved for each of the above cases with the aid of a computer program based on the method of feasible conjugate directions. The results clearly reflect the balance between the economies of scale and the time-value of money under every demand growth function.     

EMERGENCY WATER SUPPLIES FROM GROUND WATER IN HUMID REGIONS1

ABSTRACT: This paper focuses on the development and testing of a mathematical model of an emergency ground water supply operated principally during periods of low streamflow. The process of ground water withdrawal and recharge is simulated taking account of streamflow, water demand, evapotranspiration, natural and artificial recharge and increased evapotranspiration due to artificial recharge, ground water pumpage, and streamflow contribution to pumped water. The model determines whether natural recharge is possible in less time than the return period of drought and also whether artificial recharge is needed. By simulating operation over a long period of time, the model can examine different droughts of short and long duration and can test the operating rules for ground water storage development in an area. Submodels analyze the components of the operating process including ground water flow into the stream, seepage losses, stream portion of well discharge due to induced infiltration and recharge from rainfall or water spreading. The model has been tested for areas in the humid northeastern United States.     

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.     

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.     

REGULATORY REQUIREMENTS AFFECTING AQUIFER THERMAL ENERGY STORAGE1

ABSTRACT: Thermal energy storage involves the capture and storage of thermal energy (either heat or chill) during one time period for use at a later period. Storage of thermal energy in aquifers on a seasonal basis is one promising application of the technology that has been implemented in several foreign countries and is currently undergoing field testing in the U.S. Potential developers of aquifer thermal energy storage projects will face a number of regulatory requirements at the federal, state, and local level of government. These can include meeting: (1) surface land and ground water use restrictions, (2) regulations relating to withdrawal of ground water, and (3) requirements for reinjecting thermally altered ground waters. Separate permits for ground water withdrawal and reinjection may be required. The permit process is likely to involve opportunities for public comment and may involve contested proceedings.     

SIMULATION OF SALTWATER INTRUSION IN VOLUSIA COUNTY,FLORIDA1

ABSTRACT: Volusia County, in east central Florida, comprises approximately 1,200 square miles situated between the St. Johns River and the Atlantic Ocean. Most of the County is underlain by a three-aquifer system. Population centers in Volusia County, which create a large water demand, are located near the coast. Saltwater intrusion into the ground water near these population centers has led to relocation of public water supply wells further inland. Regional management of the county's water resources commissioned construction of a three-dimensional computer model of the county. Predevelopment simulation results were used as initial conditions for the development simulations, which included well discharge data. The predevelopment model calibration consisted of reproducing field-determined potentiometric surfaces. As part of the calibration process, sensitivity analyses were performed on boundary conditions, recharge rates, permeability, and leakage properties. Results of the model study indicate the utility of computer models as a management tool for the complex ground-water system in Volusia County.