ESTIMATING WATER BUDGET IN A REGIONAL AQUIFER USING HSPF‐MODFLOW INTEGRATED MODEL1

ABSTRACT: Integrated water resources management is important, especially in watersheds where substantial interactions exist between the ground and surface water sources. This management warrants the need for reliable estimates of both an overall basin water budget and hydrologic fluctuations between ground water and surface water sources. The objectives of this study were to estimate the total water budget and to simulate the effects of the management of water in the Big Lost River Basin in Idaho. The study used the FIPR Hydrological Model (FHM), a hydrological model developed by the University of South Florida for the Florida Institute of Phosphate Research (FIPR). The FHM is an integrated model that simulates the full water budget of the surface and ground water systems. It has two public domain components: Hydrological Simulation Program ‐ FORTRAN (HSPF) and Modular Three‐Dimensional Finite‐Difference Ground Water Flow Model (MODFLOW). This study quantified the hydrologic fluxes between ground water and surface water and determined a comprehensive and accurate water budget for the Big Lost River. The study showed an annual amount of 10.44 m3/sec leaves the basin and never to return to the system. The study is useful in developing and calculating the annual water budget in the Big Lost River, and this process should be applicable to estimating water budgets in other basins.     

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.     

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.     

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.     

IMPACTS OF SURFACE DRAINAGE ON GROUND WATER HYDRAULICS1

ABSTRACT: The current dredge and fill practices in locating canals along the periphery of wetlands in south Florida are transforming natural basins that originally had primarily slower subsurface drainage to ones that discharge larger quantities of water faster, via a surface drainage system. The objective of this paper is to develop an analytical technique and a numerical model in quantifying the difference of surface and subsurface runoff before and after the construction of drainage canals, and for delineating the effects of drains on channel level and regional water tables in adjacent areas in south Florida. The surface runoff model is formulated on the climatic water balance technique, and the ground water model is treated as a one dimensional transient phenomenon that forms a nonlinear flow problem. Analytical solutions are derived through problem linearization. These two models are coupled to estimate the impact of drainage canals on the adjacent water table drawdown.     

CHALLENGES IN GROUND WATER REMEDIATION: DENSE NONAQUEOUS PHASE LIQUIDS1

ABSTRACT: Dense nonaqueous phase liquid (DNAPL) contaminated ground water has proven to be exceptionally difficult to remediate for both physical and chemical reasons. Since DNAPL's are denser than water, their movement is not governed by the direction of ground water flow as is generally the case for other ground water contaminants. Additionally, DNAPLs' interactions with aquifer solids through processes such as sorption tend to make the pollution linger or sometimes apparently disappear, only to return later. Unfortunately, pump-and-treat systems, the traditional way ground water contamination is addressed, have not been effective in cleaning DNAPL contaminated water. Other remediation technologies continue to be developed to address these problems. Policy changes will also be necessary to effectively address the difficulties associated with ground water remediation.     

A System Dynamics Model for Conjunctive Management of Water Resources in the Snake River Basin

The Pacific Northwest is expected to witness changes in temperature and precipitation due to climate change. In this study, we enhance the Snake River Planning Model (SRPM) by modeling the feedback loop between incidental recharge and surface water supply resulting from surface water and groundwater extraction for irrigation and provide a case study involving climate change impacts and management scenarios. The new System Dynamics‐Snake River Planning Model (SD‐SRPM) is calibrated to flow at Box Canyon Springs located along a major outlet of the East Snake Plain Aquifer. A calibration of the model to flow at Box Canyon Springs, based on historic diversions (1950‐1995) resulted in an r² value of 0.74 and a validation (1996‐2005) r² value of 0.60. After adding irrigation entities to the model an r² value of 0.91, 0.88, and 0.87 were maintained for modeled vs. observed (1991‐2005) end‐of‐month reservoir content in Jackson Lake, Palisades, and American Falls, the three largest irrigation reservoirs in the system. The scenarios that compared the impacts of climate change were based on ensemble mean precipitation change scenarios and estimated changes to crop evapotranspiration (ET). Increased ET, despite increased precipitation, generally increased surface water shortages and discharge of springs. This study highlights the need to develop and implement models that integrate the human‐natural system to understand the impacts of climate change.     

GROUND WATER QUALITY IMPLICATIONS OF SOIL CONSERVATION MEASURES: AN ECONOMIC PERSPECTWE1

ABSTRACT: An evaluation of the intermedia movement of pesticides applied under various land management systems already in place, or to be implemented, under the Conservation Reserve and Conservation Compliance programs is presented. The simulation modeling approach followed in this analysis consists of a mathematical programming model and leaching/surface runoff, Pesticide Root Zone Model (PRZM) models. Special care was taken to ensure that the physical model was sensitive to the chemical characteristics of individual pesticides and the important physical changes brought about by different agricultural practices. Results show that, although these programs as now planned, increase farm income and achieve soil conservation goals, they may adversely affect ground water quality. Also, depending on soil and location characteristics, there are tradeoffs between surface and ground water quality implications. Hence, if these programs are to address water quality problems, the recommended practices must be evaluated for their impact on water quality, particularly in potentially vulnerable areas.     

Mason Valley Groundwater Model: Linking Surface Water and Groundwater in the Walker River Basin,Nevada1

Carroll, Rosemary W.H., Greg Pohll, David McGraw, Chris Garner, Anna Knust, Doug Boyle, Tim Minor, Scott Bassett, and Karl Pohlmann, 2010. Mason Valley Groundwater Model: Linking Surface Water and Groundwater in the Walker River Basin, Nevada. Journal of the American Water Resources Association (JAWRA) 46(3):554-573. DOI: 10.1111/j.1752-1688.2010.00434.x Abstract: An integrated surface water and groundwater model of Mason Valley, Nevada is constructed to replicate the movement of water throughout the different components of the demand side of water resources in the Walker River system. The Mason Valley groundwater surface water model (MVGSM) couples the river/drain network with agricultural demand areas and the groundwater system using MODFLOW, MODFLOW’s streamflow routing package, as well as a surface water linking algorithm developed for the project. The MVGSM is capable of simulating complex feedback mechanisms between the groundwater and surface water system that is not dependent on linearity among the related variables. The spatial scale captures important hydrologic components while the monthly stress periods allow for seasonal evaluation. A simulation spanning an 11-year record shows the methodology is robust under diverse climatic conditions. The basin-wide modeling approach predicts a river system generally gaining during the summer irrigation period but losing during winter months and extended periods of drought. River losses to the groundwater system approach 25% of the river’s annual budget. Reducing diversions to hydrologic response units will increase river flows exiting the model domain, but also has the potential to increase losses from the river to groundwater storage.     

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.     

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.     

A NONLINEAR MULTILEVEL MODEL FOR REGIONAL WATER RESOURCES PLANNING1

ABSTRACT: A nonlinear multilevel transportation model is developed to study large-scale allocations in a water resources system. The model uses a modified transportation matrix formulated with nonlinear cost functions as the basic subregional model and the goal coordination method for multilevel decomposition and optimization of the overall regional system. The model is applied to projected water requirements for Salt Lake County in 1985. Sources of water supply - surface water, ground water, import water, and reuse of reclaimed wastewater on a restricted basis - are available to satisfy water requirements for municipal, industrial, and agricultural sectors in four subregions. The conjugate gradient projection method is used to optimize the first level subregional models having cost functions of the form of C = aX^b, and the second level problem is solved using the conjugate gradient method.     

A WATER QUALITY MODEL FOR A CONJUNCTIVE SURFACE-GROUNDWATER SYSTEM: AN OVERVIEW1

ABSTRACT. A mathematical model to predict water quality in a surface-groundwater system is under development. This project is being sponsored by the Environmental Protection Agency. The ultimate goal of this study is to obtain cause and effect relationships between pollutant sources and the ensuing concentrations at different locations in a basin. Several programs are used to model the various hydrologic processes occurring in nature, namely: rainfall, runoff, flow in surface bodies of water, infiltration, and groundwater flow. At every time step in the simulation, the water quantity computations for the above hydrologic models are performed first. Subsequently, the results of these computations, typically in the form of flow velocities, are used as input to the water quality calculations. The water quality routines involve the modeling of the associated physical, chemical, and biological processes. In this study, emphasis is being placed on pollution in agricultural areas. Accordingly the Lake Apopka basin in Central Florida is being used as the application site.     

A Calibrated,High‐Resolution GOES Satellite Solar Insolation Product for a Climatology of Florida Evapotranspiration1

Paech, Simon J., John R. Mecikalski, David M. Sumner, Chandra S. Pathak, Quinlong Wu, Shafiqul Islam, and Taiye Sangoyomi, 2009. A Calibrated, High‐Resolution GOES Satellite Solar Insolation Product for a Climatology of Florida Evapotranspiration. Journal of the American Water Resources Association (JAWRA) 45(6):1328‐1342. Abstract: Estimates of incoming solar radiation (insolation) from Geostationary Operational Environmental Satellite observations have been produced for the state of Florida over a 10‐year period (1995‐2004). These insolation estimates were developed into well‐calibrated half‐hourly and daily integrated solar insolation fields over the state at 2 km resolution, in addition to a 2‐week running minimum surface albedo product. Model results of the daily integrated insolation were compared with ground‐based pyranometers, and as a result, the entire dataset was calibrated. This calibration was accomplished through a three‐step process: (1) comparison with ground‐based pyranometer measurements on clear (noncloudy) reference days, (2) correcting for a bias related to cloudiness, and (3) deriving a monthly bias correction factor. Precalibration results indicated good model performance, with a station‐averaged model error of 2.2 MJ m^?2/day (13%). Calibration reduced errors to 1.7 MJ m^?2/day (10%), and also removed temporal‐related, seasonal‐related, and satellite sensor‐related biases. The calibrated insolation dataset will subsequently be used by state of Florida Water Management Districts to produce statewide, 2‐km resolution maps of estimated daily reference and potential evapotranspiration for water management‐related activities.     

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.     

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.     

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.     

AGRICULTURAL WATER MODELING AND ECONOMIC INFORMATION NEEDS UNDER THE MODEL WATER CODE1

ABSTRACT: Economic information for efficient water allocation is difficult and costly to acquire under administrative water systems evolving under the Model Water Code. One approach to obtaining more information is to use a simulator like the Florida AGWATER model. The advantage of AGWATER is the potential for realistic prediction, because it operates at the field and day levels, using detailed information for each crop and tract. Unfortunately, such simulators are complex and require large amounts of costly input data. A better solution to the information problem may be to use markets for the marketable goods associated with water, because information is inherent in such a process. This approach will allow limited modeling and management resources to be put into using water models to generate information for the goods dependent on water that are difficult to market, like wildlife services.     

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.     

Riparian Ground‐Water Flow Patterns Using Flownet Analysis: Evapotranspiration‐Induced Upwelling and Implications for N Removal1

Abstract: Ground‐water flow paths constrain the extent of nitrogen (N) sinks in deep, stratified soils of riparian wetlands. We examined ground‐water flow paths at four forested riparian wetlands in deep, low gradient, stratified deposits subjected to Southern New England’s temperate, humid climate. Mid‐day piezometric heads were recorded during the high water table period in April/May and again in late November at one site. Coupling field data with a two‐dimensional steady‐state ground‐water flow model, flow paths and fluxes were derived to 3 m depths. April/May evapotranspiration (ET) dominated total outflux (44‐100%) while flux to the stream was <10% of total outflux. ET exerted upward ground‐water flux through shallow carbon‐rich soils, increasing opportunities for N transformations and diverting flow from the stream. Dormant season results showed a marked increase in flux to the stream (27% of the total flux). Riparian sites with deep water tables (naturally or because of increased urbanization or other hydrologic modifications) or shallow root zones may not generate ground‐water upwelling to meet evaporative demand, thereby increasing the risk of N movement to streams. As water managers balance issues of water quality with water quantity, they will be faced with decisions regarding riparian management. Further work towards refining our understanding of ET mediation of N and water flux at the catchment scale will serve to inform these decisions.