AUTOMATED METHODS FOR ESTIMATING BASEFLOW AND GROUND WATER RECHARGE FROM STREAMFLOW RECORDS1

ABSTRACT: To quantify and model the natural ground water recharge process, six sites located in the midwest and eastern United States where previous water balance observations had been made were compared to computerized techniques to estimate: (1) base flow and (2) ground water recharge. Results from an existing automated digital filter technique for separating baseflow from daily streamflow records were compared to baseflow estimates made in the six water balance studies. Previous validation of automated baseflow separation techniques consisted only of comparisons with manual techniques. In this study, the automated digital filter technique was found to compare well with measured field estimates yielding a monthly coefficient of determination of 0.86. The recharge algorithm developed in this study is an automated derivation of the Rorabaugh hydrograph recession curve displacement method that utilizes daily streamflow. Comparison of annual recharge from field water balance measurements to those computed with the automated recession curve displacement method had coefficients of determination of 0.76 and predictive efficiencies of 71 percent. Monthly estimates showed more variation and are not advocated for use with this method. These techniques appear to be fast, reproducible methods for estimating baseflow and annual recharge and should be useful in regional modeling efforts and as a quick check on mass balance techniques for shallow water table aquifers.     

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.     

GROUND WATER/SURFACE WATER RESPONSES TO GLOBAL CLIMATE SIMULATIONS,SANTA CLARA‐CALLEGUAS BASIN,VENTURA, CALIFORNIA1

ABSTRACT: Climate variations can play an important, if not always crucial, role in successful conjunctive management of ground water and surface water resources. This will require accurate accounting of the links between variations in climate, recharge, and withdrawal from the resource systems, accurate projection or predictions of the climate variations, and accurate simulation of the responses of the resource systems. To assess linkages and predictability of climate influences on conjunctive management, global climate model (GCM) simulated precipitation rates were used to estimate inflows and outflows from a regional ground water model (RGWM) of the coastal aquifers of the Santa Clara‐Calleguas Basin at Ventura, California, for 1950 to 1993. Interannual to interdecadal time scales of the El Niño Southern Oscillation (ENSO) and Pacific Decadal Oscillation (PDO) climate variations are imparted to simulated precipitation variations in the Southern California area and are realistically imparted to the simulated ground water level variations through the climate‐driven recharge (and discharge) variations. For example, the simulated average ground water level response at a key observation well in the basin to ENSO variations of tropical Pacific sea surface temperatures is 1.2 m/°C, compared to 0.9 m/°C in observations. This close agreement shows that the GCM‐RGWM combination can translate global scale climate variations into realistic local ground water responses. Probability distributions of simulated ground water level excursions above a local water level threshold for potential seawater intrusion compare well to the corresponding distributions from observations and historical RGWM simulations, demonstrating the combination's potential usefulness for water management and planning. Thus the GCM‐RGWM combination could be used for planning purposes and — when the GCM forecast skills are adequate — for near term predictions.     

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.     

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 ECONOMIC VALUE OF GROUND WATER RECHARGE FOR IRRIGATION USE1

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

FACTORS AFFECTING SPRING RUNOFF ON TWO FORESTED WATERSHEDS1

ABSTRACT Spring runoff from two forested watersheds in northern Minnesota is a function of annual snowfall, soil water recharge, and water supply rates. A drainage basin with a clay soil and a hardwood overstory had greater snowmelt and water supply rates than another drainage basin with a sandy soil and conifer overstory. The average soil water recharge rate for the clay soil was 28 percent less than for the sandy soil. The lower recharge rate of the clay soil resulted in spring runoff which averaged 40 percent of water supplied during the three year study while an average of two percent was produced on the sandy soil. Soil frost which affected soil water recharge varied between soil types and was influenced by amount of soil water storage and snow cover.     

Assessing ground water vulnerability with the type transfer function model in the San Joaquin Valley, California

The recently developed type transfer function (TTF) simulation approach was applied to generate a regional-scale nonpoint-source ground water vulnerability assessment for the San Joaquin Valley, California. The computationally comparatively inexpensive TTF approach produces quantitative estimates of contaminant concentrations for large regional scales through characteristic functions based on different soil textures and their leaching properties. The TTF simulations employed an extensive soil and recharge database to estimate atrazine (1-chloro-3-ethylamino-5-isopropylamino-2,4,6-triazine) concentrations at a compliance depth of 3 m resulting from a surface application. Two different sets of TTFs with two different levels of upscaling were used for spatially uniform and distributed recharge estimates. Results show that estimated atrazine concentrations can be related to soil survey information. Areas with high potential vulnerability to atrazine leaching were found for soils with low organic carbon content and sandy loam and loam textures. Travel times for atrazine peak concentrations to the compliance depth ranged from 350 to 730 d. The extent of areas with estimated atrazine concentrations above the maximum contaminant level was less extensive when uniform annual recharge values were used. Simulated TTF concentrations were highest for eastern Fresno County, a vulnerability pattern that is also supported by field observations. The TTF modeling approach is shown to be a useful tool for quantitative pesticide leaching estimates at regional scales significantly larger than those of previous studies.     

EVALUATION OF AN INDUCED INFILTRATION MODEL AS APPLIED TO GLACIAL AQUIFER SYSTEMS1

ABSTRACT: Numerical models were used to examine the limitations of the assumptions used in an analytical induced infiltration model. The assumptions tested included negligible streambed effects, negligible areal recharge, two-dimensional ground water flow, fully penetrating rivers and wells, and constant surface water stage. For situations that deviate from the underlying assumptions, the analytical model becomes a less reliable predictor of induced infiltration. The numerical experiments show that streambed effects cannot be neglected if the streambed conductivity is more than one order of magnitude lower than the aquifer hydraulic conductivity. Areal recharge cannot be neglected if the ground water basin receives more than 5 in/yr of areal recharge. Three-dimensional flow effects due to well partial penetration cannot be neglected if the ratio of horizontal hydraulic conductivity to vertical hydraulic conductivity (K_h/K_u) is greater than 10. Surface water elevation fluctuations can significantly influence the induced infiltration rate, depending on the degree of fluctuations and the ground water hydraulic gradient.     

DYNAMIC RESPONSE OF AQUIFER SYSTEMS TO LOCALIZED RECHARGE1

ABSTRACT. The response of stream-unconfined aquifer systems to localized recharge is investigated by means of a two-dimensional finite element model. A variational approach is used in conjunction with the finite element method to solve the ground water flow equation. Linear approximated triangular elements are used to calculate the hydraulic head distribution in the flow region. The Crank-Nicholson centered scheme of numerical integration is employed to approximate the time derivative in the flow equation. A computer program is developed to calculate the hydraulic head distribution in the flow region. Solutions provided by the finite element model should prove useful in the evaluation of quantitative and qualitative changes in aquifer systems due to natural or artificial recharge. In addition, they should prove useful in the study of irrigation and drainage problems.     

MODEL DEVELOPMENT FOR CONJUNCTIVE USE STUDY OF THE SAN JACINTO BASIN,CALIFORNIA1

ABSTRACT: A two-layered confined-unconfined numerical model for flow and mass transport is developed for the San Jacinto Basin. The model structure is determined by the geological structure of the Basin and model parameters are calibrated using 20 years of historical records. The total number of historical head observations used for the flow model calibration is 1,117 and the total number of the estimated parameters is 91. The two-layered transport model is also calibrated using historical water quality records. Sensitivity analysis of the flow model shows that only 68 parameters (out of a total of 91) are relatively sensitive and reliable. However, the unreliable parameters (23 of them) are found to be insensitive and thus not significant to the prediction and management of conjunctive use of surface water and ground water. The developed flow model has been used to study the two proposed artificial recharge scenarios for the San Jacinto Basin. We have found that during a relatively dry condition, an artificial recharge rate of 80 acre-ft/day can be achieved during the recharge period October through January. However, for a relatively wet condition, only 80 percent of the proposed rate can be effectively stored in the Basin during these months.     

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.     

Soil Moisture Data Assimilation Using Support Vector Machines and Ensemble Kalman Filter1

Abstract: A hybrid data assimilation (DA) methodology that combines two state‐of‐the‐art techniques, support vector machines (SVMs) and ensemble Kalman filter (EnKF), is applied for soil moisture DA in this work. The SVM methodology provides a statistically sound and robust approach to solving the inverse problem, and thus to building statistical models. EnKF is an extension of the Kalman Filter (KF), a well‐known tool in prediction updating. In the present research, ground measurements were used to build a SVM‐type soil moisture predictor. Subsequent observations and their statistics were assimilated to update predictions from the SVM model by coupling it with EnKF. In this way, both model predictions and ground data, as well as their statistics, are fused thus minimizing the prediction error and making the predictions and observations statistically consistent. The results are shown for two approaches; one in which update is done at every time step and the other which assumes that data is only available at alternate time steps (in window of 10 time steps) and hence update is performed at those occasions. The SVM‐EnKF coupling is shown to improve soil moisture forecasts in an example using data from the Soil Climate Analysis Network site at Ames, Iowa.     

GROUNDWATER CONTAMINATION BY NEMATICIDES: INFLUENCE OF RECHARGE TIMING UNDER PINEAPPLE CROP1

ABSTRACT: Toxic organic compounds, such as DBCP, EDB, and c TCP, that are associated with pineapple cultivation in Hawaii have been discovered in drinking water wells on Oahu. In order to reach and contaminate the Pearl Harbor aquifer, pesticides must be transported quickly downward away from the soil surface prior to complete volatilization, degradation, or adsorption of residuals. This paper assesses the role of pesticide application timing relative to subsequent rainfall-induced recharge events in determining the amount and extent of chemical leaching from the soil. A water balance model for a pineapple crop is developed to estimate the time series of recharge from two fields for which soil contamination profiles are available. In general, the amounts of DBCP, EDB, and TCP found in the soil profiles of the two fields are consistent with expectations of leaching based on an analysis of the recharge time series. The results indicate that recharge during and immediately following the application of pesticides is important in determining whether groundwater contamination will result.     

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.     

Field-scale preferential transport of water and chloride tracer by depression-focused recharge

A tracer study was initiated in November 1993 to investigate depression-focused recharge and to monitor solute movement through the vadose zone into the shallow ground water in southeastern North Dakota. Granular potassium chloride (KCl) was surface-applied to two areas overlying subsurface drains and to one area instrumented with soil solution samplers, ground water monitoring wells, time domain reflectometry (TDR) probes, and temperature probes. One of the subsurface drain tracer plots was located on level ground while the other two sites were in small topographic depressions. Formation of ground water mounds beneath the depressions indicated that these areas are recharge sites. The applied Cl- tracer was found to move rapidly to the shallow ground water under the depressional areas after infiltration of spring snowmelt in 1994. Excessive rainfall events were also responsible for focused recharge and the rapid transport of the applied Cl- tracer. Water flow through partially frozen soil at the bottom of the depressions during thaw enhanced preferential movement of the tracer.     

Ion partitioning among soil and plant components under drip,furrow, and sprinkler irrigation regimes: field and modeling assessments

Soil and water resources can be severely degraded by salinity when total salt input exceeds output in irrigated agriculture. This study was conducted to examine partitioning of Ca2+, Na+, and Cl- between soil and soybean [Glycine max (L.) Merr.] plants under different irrigation regimes with both field and modeling assessments. In drip and sprinkler treatments, the irrigation water was salinized with NaCl and CaCl2 salts to simulate a Cl- and Na+ dominant saline drainage water. In the furrow irrigation treatment, the soil was salinized, prior to planting, with NaCl and CaCl2 salts to simulate a Cl- and Na+ dominant saline soil. A total of 756 soil and 864 plant samples were collected and analyzed for the salt ions to obtain ion partitioning and mass balance assessments. Modeling of salt ion uptake by plants and distribution in the soil profile was performed with a two-dimensional solute transport model for the three irrigation regimes. Results indicated that about 20% of the applied Ca2+ was recovered in harvested soybean biomass in all treatments. Plant uptake of either Na+ or Cl- was less than 0.5% in the drip and furrow, and about 2% in the sprinkler irrigation treatment. Significant increases in soil salinity were found in the sprinkler plot that received the highest cumulative amount of salts. Simulated ion distributions in the soil were comparable with the measurements. Compared with the total seasonal salt input, mass balances between 65 and 108% were obtained. Most salt inputs accumulate in the soil, and need to be removed periodically to prevent soil salinization.     

MACROSCALE HYDROLOGIC MODELING FOR REGIONAL CLIMATE ASSESSMENT STUDIES IN THE SOUTHEASTERN UMTED STATES1

ABSTRACT: A macroscale hydrologic model is developed for regional climate assessment studies under way in the southeastern United States. The hydrologic modeling strategy is developed to optimize spatial representation of basin characteristics while maximizing computational efficiency. The model employs the “grouped response unit” methodology, which follows the natural drainage pattern of the area. First order streams are delineated and their surface characteristics are tested so that areas with statistically similar characteristics can be combined into larger computational zones for modeling purposes. Hydrologic response units (HRU) are identified within the modeling units and a simple three‐layer water balance model, Soil and Water Assessment Tool (SWAT), is executed for each HRU. The runoff values are then convoluted using a triangular unit hydrograph and routed by Muskingum‐Cunge method. The methodology is shown to produce accurate results relative to other studies, when compared to observations. The model is used to evaluate the potential error in hydrologic assessments when using GCM predictions as climatic input in a rainfall‐runoff dominated environment. In such areas, the results from this study, although limited in temporal and spatial scope, appear to imply that use of GCM climate predictions in short term quantitative analyses studies in rainfall‐runoff dominated environments should proceed with caution.     

Modeling the Spatially Varying Water Balance Processes in a Semiarid Mountainous Watershed of Idaho1

Stratton, Benjamin T., Venakataramana Sridhar, Molly M. Gribb, James P. McNamara, and Balaji Narasimhan, 2009. Modeling the Spatially Varying Water Balance Processes in a Semiarid Mountainous Watershed of Idaho. Journal of the American Water Resources Association (JAWRA) 45(6):1390‐1408. Abstract: The distributed Soil Water Assessment Tool (SWAT) hydrologic model was applied to a research watershed, the Dry Creek Experimental Watershed, near Boise Idaho to investigate its water balance components both temporally and spatially. Calibrating and validating SWAT is necessary to enable our understanding of the water balance components in this semiarid watershed. Daily streamflow data from four streamflow gages were used for calibration and validation of the model. Monthly estimates of streamflow during the calibration phase by SWAT produced satisfactory results with a Nash Sutcliffe coefficient of model efficiency 0.79. Since it is a continuous simulation model, as opposed to an event‐based model, it demonstrated the limited ability in capturing both streamflow and soil moisture for selected rain‐on‐snow (ROS) events during the validation period between 2005 and 2007. Especially, soil moisture was generally underestimated compared with observations from two monitoring pits. However, our implementation of SWAT showed that seasonal and annual water balance partitioning of precipitation into evapotranspiration, streamflow, soil moisture, and drainage was not only possible but closely followed the trends of a typical semiarid watershed in the intermountain west. This study highlights the necessity for better techniques to precisely identify and drive the model with commonly observed climatic inversion‐related snowmelt or ROS weather events. Estimation of key parameters pertaining to soil (e.g., available water content and saturated hydraulic conductivity), snow (e.g., lapse rates, melting), and vegetation (e.g., leaf area index and maximum canopy index) using additional field observations in the watershed is critical for better prediction.     

A REEVALUATION OF THE GROUND WATER BUDGET FOR LAS VEGAS VALLEY,NEVADA, WITH EMPHASIS ON GROUND WATER DISCHARGE1

ABSTRACT: An essential component to the ground water budget for the Las Vegas Valley (LVV) in southern Nevada is discharge from the ground water system. Discharge for the LW has been based on estimates made more than 50 years ago of 35,524,224 m3 per year as evapotranspiration (ET) and 0 m³ per year as subsurface outflow. Newly published values for recharge based on a more robust data set (70,308,360 m³) indicate a large imbalance associated with the earlier discharge estimates, providing the basis for the reevaluation conducted in this study. ET estimates in this study, as opposed to previous studies, were assigned a range in values that included an approach that assigned higher weight to the unique soil, plant, water, and climatic conditions that existed in predevelopment (1905) LW. The earlier discharge estimates also assumed that the basin was hydrologically closed; however, based on our evaluation, a range in yearly discharge by subsurface outflow from 1,480,176 m³ to 19,735,680 m³ could be assigned. Likewise, a range in yearly ET from 20,475,768 m³ to 78,819,372 m³ could be assigned. Based on newly published recharge values, closure can only occur if higher values are assigned to both the subsurface outflow and/or ET components of ground water discharge. We cannot provide a complete water balance closure with our ground water discharge estimate of 64,140,960 m³. However our reevaluation gives support to the higher recharge estimates and provides the rationale for future studies to be conducted based on a more rigorous scientific assessment.