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.     

MODELING RUNOFF FROM VARIABLE SOURCE AREAS IN HUMID,SHALLOW WATER TABLE ENVIRONMENTS1

ABSTRACT: Variable Source Areas (VSAs) are zones with water saturated soils in forested wetlands fringing streams and creeks. Runoff from these areas is generated by saturation excess after a shallow water table rises and inundates the ground surface. In humid regions, like Florida and the Southeast, VSAs are believed to produce most of the runoff in shallow water table environments. Modeling the spatial extent and temporal fluctuation of a VSA is difficult because the formation of a VSA depends on a number of hydrological and morphological factors like rainfall intensity, soil texture, water table depth, and topographic attributes of the terrain. In this paper, we couple a digital elevation model with a two‐dimensional variable saturation model to illustrate the formation of a VSA at the hillside scale. The topography derived from the digital elevation model forms the upper domain geometry for the two‐dimensional finite element simulations of variable saturated flow. The objectives are: (1) to model the spatial and dynamic fluctuation of a VSA, and (2) to understand the roles of rainfall variability and terrain attributes on the formation of a VSA. Results show that hillsides with shallow water table depths, low saturated hydraulic conductivity, mild slopes, and concave slope curvature were more susceptible to runoff from a variable source. Runoff from a variable source showed little sensitivity to rainfall intensity. In general, landscapes with steep slopes generated a small VSA and a seepage face that vanished rapidly with time. In contrast, flat terrains are more amenable to VSA and retain ground surface inundation for longer periods of time.     

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.     

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.     

Physical and Temporal Isolation of Mountain Headwater Streams in the Western Mojave Desert,Southern California1

Abstract: Streams draining mountain headwater areas of the western Mojave Desert are commonly physically isolated from downstream hydrologic systems such as springs, playa lakes, wetlands, or larger streams and rivers by stream reaches that are dry much of the time. The physical isolation of surface flow in these streams may be broken for brief periods after rainfall or snowmelt when runoff is sufficient to allow flow along the entire stream reach. Despite the physical isolation of surface flow in these streams, they are an integral part of the hydrologic cycle. Water infiltrated from headwater streams moves through the unsaturated zone to recharge the underlying ground‐water system and eventually discharges to support springs, streamflow, isolated wetlands, or native vegetation. Water movement through thick unsaturated zones may require several hundred years and subsequent movement through the underlying ground‐water systems may require many thousands of years – contributing to the temporal isolation of mountain headwater streams.     

GENERALIZED WATERSHED LOADING FUNCTIONS FOR STREAM FLOW NUTRIENTS1

ABSTRACT: Loading functions are proposed as a general model for estimating monthly nitrogen and phosphorus fluxes in stream flow. The functions have a simple mathematical structure, describe a wide range of rural and urban nonpoint sources, and couple surface runoff and ground water discharge. Rural runoff loads are computed from daily runoff and erosion and monthly sediment yield calculations. Urban runoff loads are based on daily nutrient accumulation rates and exponential wash off functions. Ground water discharge is determined by lumped parameter unsaturated and saturated zone soil moisture balances. Default values for model chemical parameters were estimated from literature values. Validation studies over a three-year period for an 850 km² watershed showed that the loading functions explained at least 90 percent of the observed monthly variation in dissolved and total nitrogen and phosphorus fluxes in stream flow. Errors in model predictions of mean monthly fluxes were: dissolved phosphorus - 4 percent; total phosphorus - 2 percent; dissolved nitrogen - 18 percent; and total nitrogen - 28 percent. These results were obtained without model calibration.     

EQUIVALENCE BETWEEN TOPMODEL AND THE NRCS CURVE NUMBER METHOD IN PREDICTING VARIABLE RUNOFF SOURCE AREAS1

An equivalence is proposed between two rainfall‐runoff methods with a long history of use in the United States and Europe. In watersheds where variable source areas dominate runoff, the two methods can have comparable probability distribution functions of moisture deficit, and therefore predict similar saturated runoff source areas. A novel approach is introduced to determine the S parameter in the Natural Resources Conservation Service (NRCS) method. This approach constrains S by the physical soil and topography characteristics of the catchment and depth to water table. The NRCS curve number method is at the core of many rainfall‐runoff models in hydrology. As a simple lumped parameter method, it is often scrutinized because it is not obvious how to derive S from catchment hydromorphological characteristics. The novel approach provides a clear physical meaning for S, allowing better estimation of this parameter in humid shallow water table environments where the variable source area can be the dominant runoff mechanism.     

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.     

A LANDSCAPE‐BASED APPROACH TO ESTIMATE RIPARIAN HYDROLOGICAL AND NITRATE REMOVAL FUNCTIONS1

ABSTRACT: This study evaluates a conceptual model developed for riparian zones in Ontario, Canada, that links landscape hydrogeological characteristics to riparian ground water hydrology and nitrate removal efficiency. Data from a range of riparian sites in the United States and Europe suggest that the riparian zone types identified in the model are consistent with patterns of riparian hydrology and nitrate flux and removal in many humid temperate landscapes. These data also support the view that a riparian width of less than 20 m is often sufficient for effective nitrate removal unless riparian sediments are coarse grained or nitrate transport occurs mainly in surface‐fed ground water seeps. This study assesses the possibility of using topographic, soil, surficial geology, and vegetation maps to determine landscape attributes linked by the model to riparian zone hydrological functioning and nitrate removal efficiency. Although mappable data can help in determining broad classes of riparian zones, field visits are necessary to determine non‐mappable riparian attributes such as seeps, organic horizons, and permeable sediment depth in the riparian zone. This research suggests that the conceptual model could be used for landscape management purposes in most temperate landscapes with minor modifications and that the hydrological component of the model could be adapted for contaminants other than nitrate.     

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

The Role of Ground Water in Generating Streamflow in Headwater Areas and in Maintaining Base Flow1

Abstract: The volume and sustainability of streamflow from headwaters to downstream reaches commonly depend on contributions from ground water. Streams that begin in extensive aquifers generally have a stable point of origin and substantial discharge in their headwaters. In contrast, streams that begin as discharge from rocks or sediments having low permeability have a point of origin that moves up and down the channel seasonally, have small incipient discharge, and commonly go dry. Nearly all streams need to have some contribution from ground water in order to provide reliable habitat for aquatic organisms. Natural processes and human activities can have a substantial effect on the flow of streams between their headwaters and downstream reaches. Streams lose water to ground water when and where their head is higher than the contiguous water table. Although very common in arid regions, loss of stream water to ground water also is relatively common in humid regions. Evaporation, as well as transpiration from riparian vegetation, causing ground‐water levels to decline also can cause loss of stream water. Human withdrawal of ground water commonly causes streamflow to decline, and in some regions has caused streams to cease flowing.     

A CONCEPTUAL FRAMEWORK FOR EVALUATING AGRICULTURAL WATER CONSERVATION POTENTIALS1

ABSTRACT: This paper defines types of water losses in irrigated agriculture and outlines potentials for water conservation. Recoverable water “losses” (seepage, leakage, and spillage during storage and conveyance, and surface runoff and deep percolation during irrigation) and irrecoverable losses (evaporation from water and soil surfaces and transpiration from plants) are described and illustrated. Some conservation terms are defined, particularly the distinction between on-farm irrigation efficiency and areawide efficiency. Briefly reviewed are agricultural water conservation technologies and their applicability. The biggest untapped potential for water conservation may be a reduction in irrecoverable losses, especially evapotranspiration. The advantages and disadvantages of reducing recoverable and irrecoverable water losses are described, including possible effects on ground water, energy, salinity, crops, wildlife, and in-stream uses. Such information may be useful in several policy and management issues, e.g., ground water overdraft and possible constraints on crops and sites to be irrigated.     

USING INFILTRATION ENHANCEMENT AND SOIL WATER MANAGEMENT TO REDUCE DIAZINON IN RUNOFF1

ABSTRACT: Pesticide runoff from dormant sprayed orchards is a major water quality problem in California's Central Valley. During the past several years, diazinon levels in the Sacramento and San Joaquin Rivers have exceeded water quality criteria for aquatic organisms. Orchard water management, via post‐application irrigation, and infiltration enhancement, through the use of a vegetative ground cover, are management practices that are believed to reduce pesticide loading to surface waters. Field experiments were conducted in Davis, California, to measure the effectiveness of these management practices in reducing the toxicity of storm water runoff. Treatments using a vegetative ground cover significantly reduced peak concentrations and cumulative pesticide mass in runoff for first flush experiments compared with bare soil treatments. Post‐application irrigation was found to be an effective means of reducing peak concentrations and cumulative mass in runoff from bare soil treatments, but showed no significant effect on vegetated treatments.     

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.     

Upland Controls on the Hydrological Functioning of Riparian Zones in Glacial Till Valleys of the Midwest1

Abstract: Identifying relationships between landscape hydrogeological setting, riparian hydrological functioning and riparian zone sensitivity to climate and water quality changes is critical in order to best use riparian zones as best management practices in the future. In this study, we investigate water table dynamics, water flow path and the relative importance of precipitation, deep ground water (DG) and seep water as sources of water to a riparian zone in a deeply incised glacial till valley of the Midwest. Data indicate that water table fluctuations are strongly influenced by soil texture and to a lesser extent by upland sediment stratigraphy producing seeps near the slope bottom. The occurrence of till in the upland and at 1.7‐2 m in the riparian zone contributes to maintaining flow parallel to the ground surface at this site. Lateral ground‐water fluxes at this site with a steep topography in the upland (16%) and loam soil near the slope bottom are small (<10 l/d/m stream length) and intermittent. A shift in flow path from a lateral direction to a down valley direction is observed in the summer despite the steep concave topography and the occurrence of seeps at the slope bottom. Principal component and discriminant analysis indicate that riparian water is most similar to seep water throughout the year and that DG originating from imbedded sand and gravel layers in the lower till unit is not a major source of water to riparian zones in this setting. Water quality data and the dependence of the riparian zone for recharge on seep water suggest that sites in this setting may be highly sensitive to changes in precipitation and water quality in the upland in the future. A conceptual framework describing the hydrological functioning of riparian zones on this setting is presented to generalize the finding of this study.     

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

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

Ground water stratification and delivery of nitrate to an incised stream under varying flow conditions

Ground water processes affecting seasonal variations of surface water nitrate concentrations were investigated in an incised first-order stream in an agricultural watershed with a riparian forest in the coastal plain of Maryland. Aquifer characteristics including sediment stratigraphy, geochemistry, and hydraulic properties were examined in combination with chemical and isotopic analyses of ground water, macropore discharge, and stream water. The ground water flow system exhibits vertical stratification of hydraulic properties and redox conditions, with sub-horizontal boundaries that extend beneath the field and adjacent riparian forest. Below the minimum water table position, ground water age gradients indicate low recharge rates (2-5 cm yr(-1)) and long residence times (years to decades), whereas the transient ground water wedge between the maximum and minimum water table positions has a relatively short residence time (months to years), partly because of an upward increase in hydraulic conductivity. Oxygen reduction and denitrification in recharging ground waters are coupled with pyrite oxidation near the minimum water table elevation in a mottled weathering zone in Tertiary marine glauconitic sediments. The incised stream had high nitrate concentrations during high flow conditions when much of the ground water was transmitted rapidly across the riparian zone in a shallow oxic aquifer wedge with abundant outflow macropores, and low nitrate concentrations during low flow conditions when the oxic wedge was smaller and stream discharge was dominated by upwelling from the deeper denitrified parts of the aquifer. Results from this and similar studies illustrate the importance of near-stream geomorphology and subsurface geology as controls of riparian zone function and delivery of nitrate to streams in agricultural watersheds.     

Comparative study of transport processes of nitrogen, phosphorus, and herbicides to streams in five agricultural basins, USA

Agricultural chemical transport to surface water and the linkage to other hydrological compartments, principally ground water, was investigated at five watersheds in semiarid to humid climatic settings. Chemical transport was affected by storm water runoff, soil drainage, irrigation, and how streams were linked to shallow ground water systems. Irrigation practices and timing of chemical use greatly affected nutrient and pesticide transport in the semiarid basins. Irrigation with imported water tended to increase ground water and chemical transport, whereas the use of locally pumped irrigation water may eliminate connections between streams and ground water, resulting in lower annual loads. Drainage pathways in humid environments are important because the loads may be transported in tile drains, or through varying combinations of ground water discharge, and overland flow. In most cases, overland flow contributed the greatest loads, but a significant portion of the annual load of nitrate and some pesticide degradates can be transported under base-flow conditions. The highest basin yields for nitrate were measured in a semiarid irrigated system that used imported water and in a stream dominated by tile drainage in a humid environment. Pesticide loads, as a percent of actual use (LAPU), showed the effects of climate and geohydrologic conditions. The LAPU values in the semiarid study basin in Washington were generally low because most of the load was transported in ground water discharge to the stream. When herbicides are applied during the rainy season in a semiarid setting, such as simazine in the California basin, LAPU values are similar to those in the Midwest basins.     

AN INTEGRATED APPROACH TO WATER RESOURCES DEVELOPMENT OF THE TEHRAN REGION IN IRAN1

ABSTRACT: In this paper, a system approach to water resources development in Tehran Metropolitan Area, with its complex system of water supply and demands, is discussed. Water resources in this region include water storage in the Lar, Latyan, and Karaj reservoirs, the Tehran aquifer, as well as water discharge in local rivers and in drainage channels (mainly supplied by urban runoff and wastewater). This study consists of three phases of long‐term water resources planning and management in the Tehran metropolitan area. In each phase, a different level of details among different components of the system is considered. In the first phase, optimal operating policies for Tehran reservoirs and a decision support system are developed. In the second phase, interactions between surface and ground water resources as well as surface runoffs and wastewater disposal in different subareas are investigated. The water table fluctuations as a result of implementing sewerage collection project was also simulated. In the last phase, long‐term scenarios for water resources and agricultural development in the Southern part of Tehran are defined, and the effects of each scenario on the quality and quantity of surface and ground water resources are studied.     

HERBICIDES AND TRANSFORMATION PRODUCTS IN SURFACE WATERS OF THE MIDWESTERN UNITED STATES1

ABSTRACT: Most herbicides applied to crops are adsorbed by plants or transformed (degraded) in the soil, but small fractions are lost from fields and either move to streams in overland runoff, near surface flow, or subsurface drains, or they infiltrate slowly to ground water. Herbicide transformation products (TPs) can be more or less mobile and more or less toxic in the environment than their source herbicides. To obtain information on the concentrations of selected herbicides and TPs in surface waters of the Midwestern United States, 151 water samples were collected from 71 streams and five reservoir outflows in 1998. These samples were analyzed for 13 herbicides and 10 herbicide TPs. Herbicide TPs were found to occur as frequently or more frequently than source herbicides and at concentrations that were often larger than their source herbicides. Most samples contained a mixture of more than 10 different herbicides or TPs. The ratios of TPs to herbicide concentrations can be used to determine the source of herbicides in streams. Results of a two‐component mixing model suggest that on average 90 percent or more of the herbicide mass in Midwestern streams during early summer runoff events originates from the runoff and 10 percent or less comes from increased ground water discharge.