SIMULATING FLOW IN REGIONAL WETLANDS WITH THE MODFLOW WETLANDS PACKAGE1

ABSTRACT: As part of the Comprehensive Everglades Restoration Plan (CERP), various water supply projects have been proposed in a region located between the Miami metropolitan area and the extensive regional wetland systems that are part of the Everglades or remnant Everglades. A ground water flow model of the surficial aquifer within northern Miami‐Dade County was constructed using MODFLOW to evaluate the effects of these projects on water levels in the wetlands and the underlying surficial aquifer. The new Wetlands package was used to conjunctively simulate overland flow through these wetlands and the shallow ground water system. Comparisons of simulated to measured ground water levels and wetland stages were very satisfactory, where computed and measured water levels agreed within 0.5 ft over most of the period of record at nearly all of the monitoring sites. Temporal trends in water levels were also replicated. It was concluded that the assumptions and methodologies inherent to the Wetlands package were suitable for simulating regional wetland hydrology within the Everglades area.     

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.     

GEOSTATISTICAL ESTIMATION OF HORIZONTAL HYDRAULIC CONDUCTIVITY FOR THE KIRKWOOD‐COHANSEY AQUIFER1

ABSTRACT: The Kirkwood‐Cohansey aquifer has been identified as a critical source for meeting existing and expected water supply needs for southern New Jersey. Several contaminated sites exist in the region; their impact on the aquifer has to be evaluated using ground water flow and transport models. Ground water modeling depends on availability of measured hydrogeologic data (e.g., hydraulic conductivity, for parameterization of the modeling runs). However, field measurements of such critical data have inadequate spatial density, and their locations are often clustered. The goal of this study was to research, compile, and geocode existing data, then use geostatistics and advanced mapping methods to develop a map of horizontal hydraulic conductivity for the Kirkwood‐Cohansey aquifer. Spatial interpolation of horizontal hydraulic conductivity measurements was performed using the Bayesian Maximum Entropy (BME) Method implemented in the BMELib code library. This involved the integration of actual measurements with soft information on likely ranges of hydraulic conductivity at a given location to obtain estimate maps. The estimation error variance maps provide an insight into the uncertainty associated with the estimates, and indicate areas where more information on hydraulic conductivity is required.     

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.     

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.     

SHALLOW AQUIFERS RELATIVE TO SURFACE WATER,LOWER NORTH PLATTE RIVER VALLEY,WYOMING1

ABSTRACT. The occurrence of ground water in the lower North Platte Valley, Goshen County, Wyoming, was studied to determine safe yield within the alluvial aquifer under varying discharge and recharge conditions. The alluvium of the North Platte is extensively developed for irrigation purposes and the effects of large-scale pumpage are of major concern. Actual withdrawals are estimated to be 46,000 acre-feet. Should pumping reach potentially higher levels an overdraft is expected. Effect of ground water withdrawals are established from projections of the flow regime within the alluvial aquifer. A time dependent, numerical model was employed to predict aquifer response to increased withdrawals. The results suggest that more efficient use of surface waters and/or increased use of ground water will reduce the annual subsurface return flow to the North Platte River and its tributaries by an amount equal to the reduced ground water recharge increment. Alternatives are available for management of the lower North Platte alluvial aquifer. The preferred course is to correlate surface and subsurface water rights, in light of convenience, economics, and best means of storage for maximum utilization of the single water resource.     

Modeling Hydrology in a Small Rocky Mountain Watershed Serving Large Urban Populations1

Abstract: This article describes the development of a calibrated hydrologic model for the Blue River watershed (867 km²) in Summit County, Colorado. This watershed provides drinking water to over a third of Colorado’s population. However, more research on model calibration and development for small mountain watersheds is needed. This work required integration of subsurface and surface hydrology using GIS data, and included aspects unique to mountain watersheds such as snow hydrology, high ground‐water gradients, and large differences in climate between the headwaters and outlet. Given the importance of this particular watershed as a major urban drinking‐water source, the rapid development occurring in small mountain watersheds, and the importance of Rocky Mountain water in the arid and semiarid West, it is useful to describe calibrated watershed modeling efforts in this watershed. The model used was Soil and Water Assessment Tool (SWAT). An accurate model of the hydrologic cycle required incorporation of mountain hydrology‐specific processes. Snowmelt and snow formation parameters, as well as several ground‐water parameters, were the most important calibration factors. Comparison of simulated and observed streamflow hydrographs at two U.S. Geological Survey gaging stations resulted in good fits to average monthly values (0.71 Nash‐Sutcliffe coefficient). With this capability, future assessments of point‐source and nonpoint‐source pollutant transport are possible.     

DEVELOPMENT OF A FRACTURED MEDIUM FLOW MODEL FOR A WASTE SITE1

ABSTRACT: A three‐dimensional fractured medium flow model was developed for the Bear Creek Valley (BCV) S‐3 site of the Oak Ridge Reservation (ORR) using SWIFT III. The numerical modeling for this site focused on a conceptual model established through the analysis of heterogeneous geologic units and matrix fracture properties of the subsurface in the BCV area. The SWIFT III modeling analysis was based on the previous modeling studies that used MODFLOW and MODPATH. A rigorous calibration was obtained first by comparing simulated results with the existing data on ground water levels and then by comparing pumping test results with the simulated ground water levels. A satisfactory agreement between observed and simulated results was obtained. The calibrated model was used to determine sustained yield from a ground water interceptor trench. Different withdrawal rates were used to simulate the performance of the trench for the sustained withdrawal of ground water.     

Evaluation of Rio Grande Management Alternatives Using a Surface‐Water/Ground‐Water Model1

Abstract: Previous investigations observed significant seepage losses from the Rio Grande to the shallow aquifer between Socorro and San Antonio, New Mexico. High‐resolution telescopic modeling was used along a 10‐km reach of the Rio Grande and associated drains and canals to evaluate several management alternatives aimed at improving river conveyance efficiency. Observed data consisted of ground‐water and surface‐water elevations, seepage rates along the Rio Grande and associated canals and drains, and borehole geology. Model calibration was achieved by adjusting hydraulic conductivity and specific storage until the output matched observed data. Sensitivity analyses indicated that the system was responsive to changes in hydrogeologic properties, especially when such alterations increased vertical connectivity between layers. The calibrated model predicted that removal of the low flow conveyance channel, a major channel draining the valley, would not only decrease river seepage by 67%, but also decrease total flow through the reach by 75%. The decreased flow through the reach would result in increased water logging and an average increase in ground‐water elevations of 1.21 meter. Simulations of the system with reduced riparian evapotranspiration rates or a relocated river channel also predicted decreased river seepage, but to a much lesser degree.     

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

INTEGRATING GEOGRAPHIC INFORMATION SYSTEMS IN GROUND-WATER APPLICATIONS USING NUMERICAL MODELING TECHNIQUES1

ABSTRACT: The Sand and Gravel Aquifer is the sole source of potable water in Escambia County, Florida. In order to better understand the hydraulics of the aquifer, a numerical computer model of the aquifer was developed. The model applied a finite element technique which allowed for density-dependent transport and flow in three dimensions. The modeling technique was integrated with GIS to develop a system for optimal management of the resource. The GIS was the primary tool in the development of the model grid, as well as being the integral component in the modeling procedure. Multiple model grids were developed for simulating regional flow and local flow/transport phenomena. The model grids were generated by the GIS where nodal and element sequencing were recorded. The grid topology was stored in the GIS with the element numbers, node numbers, and the related hydrogeologic attributes. The Triangulated Irregular Network (TIN) module was used for transferring interpolated value between GIS coverages. TIN allowed a fit of the model grid to the physical dimensions of the aquifer and for interpolating boundary values for telescopically refined grids. Calculations between TIN surfaces provided the residuals of the dependent variable from observed TIN surfaces. Model calibration was conducted within the GIS environment through a combination of visual and relational querying. The GIS provided an integrated environment which facilitated model analyses and data storage and retrieval.     

EDWARDS AQUIFER EVALUATION: KINNEY COUNTY,TEXAS1

ABSTRACT: The Edwards Aquifer is one of the most studied and most prolific aquifers in the United States. The aquifer is a heavily fractured and faulted carbonate aquifer with transmissivities in excess of 100 ft²/s. The City of San Antonio relies upon the Edwards Aquifer as its sole source for water. Much work has been done on quantifying recharge to the aquifer and discharge from wells and acquiring aquifer characteristics from pumping tests, specific capacity tests, and geophysical logs. Although the aquifer has been well studied in Bexar County, much less is known about the Edwards Aquifer in Kinney County. This is partly due to the lower population within the county (approximately 3,500 people) relative to the eastern counties (Uvalde, Medina, Bexar, Comal, and Hays) and the great distance of Kinney County from high profile discharge areas such as the City of San Antonio and Comal and San Marcos Springs. Three key products resulted from this study: (1) exploratory well drilling and the largest aquifer test in the county that were conducted to evaluate the well yields within a 10,000 acre study area in which a drawdown of 2.5 ft approximately 1.2 miles away was observed while pumping at approximately 4,600 gpm; (2) a recharge estimate for the Edwards Aquifer within Kinney County of approximately 71,382 ac‐ft/yr; and (3) locating the Brackettville Groundwater Divide from an evaluation of ground water flow direction and hydrograph analysis. These results help evaluate the complex hydraulics occurring within Kinney County and aid in development of ground water modeling that will be used in managing the Edwards Aquifer.     

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.     

A Retrospective Look at the Water Resource Management Policies in Nassau County,Long Island,New York1

Abstract: The residents of Nassau County Long Island, New York receive all of their potable drinking water from the Upper Glacial, Jameco/Magothy (Magothy), North Shore, and Lloyd aquifers. As the population of Nassau County grew from 1930 to 1970, the demand on the ground‐water resources also grew. However, no one was looking at the potential impact of withdrawing up to 180 mgd (7.9 m³/s) by over 50 independent water purveyors. Some coastal community wells on the north and south shores of Nassau County were being impacted by saltwater intrusion. The New York State Legislature formed a commission to look into the water resources in 1972. The commission projected extensive population growth and a corresponding increase in pumping resulting in a projected 93.5 to 123 mgd (4.1 to 5.5 m³/s) deficit by 2000. In 1986, the New York Legislature passed legislation to strengthen the well permit program and also establish a moratorium on new withdrawals from the Lloyd aquifer to protect the coastal community’s only remaining supply of drinking water. Over 30 years has passed since the New York Legislature made these population and pumping projections and it is time to take a look at the accuracy of the projections that led to the moratorium. United States Census data shows that the population of Nassau County did not increase but decreased from 1970 to 2000. Records show that pumping in Nassau County was relatively stable fluctuating between 170 and 200 mgd (7.5 to 8.8 m³/s) from 1970 to 2004, well below the projection of 242 to 321 mgd (10.6 to 14.1 m³/s). Therefore, the population and water demand never grew to projected values and the projected threat to the coastal communities has diminished. With a stable population and water demand, its time to take a fresh look at proactive ground‐water resource management in Nassau County. One example of proactive ground‐water management that is being considered in New Jersey where conditions are similar uses a ground‐water flow model to balance ground water withdrawals, an interconnection model to match supply with demand using available interconnections, and a hydraulic model to balance flow in water mains. New Jersey also conducted an interconnection study to look into how systems with excess capacity could be used to balance withdrawals in stressed aquifer areas with withdrawals in unstressed areas. Using these proactive ground‐water management tools, ground‐water extraction could be balanced across Nassau County to mitigate potential impacts from saltwater intrusion and provide most water purveyors with a redundant supply that could be used during water emergencies.     

Transport and fate of nitrate and pesticides: hydrogeology and riparian zone processes

There is continuing concern over potential impacts of widespread application of nutrients and pesticides on ground- and surface-water quality. Transport and fate of nitrate and pesticides were investigated in a shallow aquifer and adjacent stream, Cow Castle Creek, in Orangeburg County, South Carolina. Pesticide and pesticide degradate concentrations were detected in ground water with greatest frequency and largest concentrations directly beneath and downgradient from the corn (Zea mays L.) field where they were applied. In almost all samples in which they were detected, concentrations of pesticide degradates greatly exceeded those of parent compounds, and were still present in ground waters that were recharged during the previous 18 yr. The absence of both parent and degradate compounds in samples collected from deeper in the aquifer suggests that this persistence is limited or that the ground water had recharged before use of the pesticide. Concentrations of NO(-)(3) in ground water decreased with increasing depth and age, but denitrification was not a dominant controlling factor. Hydrologic and chemical data indicated that ground water discharges to the creek and chemical exchange takes place within the upper 0.7 m of the streambed. Ground water had its greatest influence on surface-water chemistry during low-flow periods, causing a decrease in concentrations of Cl(-), NO(-)(3), pesticides, and pesticide degradates. Conversely, shallow subsurface drainage dominates stream chemistry during high-flow periods, increasing stream concentrations of Cl(-), NO(-)(3), pesticides, and pesticide degradates. These results point out the importance of understanding the hydrogeologic setting when investigating transport and fate of contaminants in ground water and surface water.     

HYDROLOGIC ANALYSIS OF DISCHARGE SUSTAINABILITY FROM AN ABANDONED UNDERGROUND COAL MINE1

ABSTRACT: Discharge from flooded abandoned subsurface coal mines is considered a potential source for water supplies where other acceptable water sources are not available. The objective of this study was to develop procedures for determining sustainability of mine‐water discharge using rainfall and discharge data for a case study site. The study site is located in southwest Virginia where Late Paleozoic sequences of sandstone, coal, and shale predominate. A rain gauge and a flow rate monitoring system were installed at the site and data were collected for a period of 100 days. The recording period corresponded with one of the driest periods in recent years and, therefore, provided valuable information regarding the flow sustainability during baseflow conditions. From available data on underground mining patterns, geology, and ground water flow regimes, it was determined that a coal mine aquifer exhibits hydraulic characteristics very similar to the extremely heterogeneous systems observed in karst aquifers, and the mine discharge is analogous to springflow. Thus, techniques commonly used in karst‐water systems and springflow analysis were used to develop rainfall/mine‐discharge relationships. Springflow recession analysis was performed on five rainfall recessions and the coefficient for each recession was compared and interpreted in light of known geologic information. It was found that the recession coefficients described the mine discharge adequately and the mine aquifer response to a rainfall pulse was very similar to the response from certain types of karst aquifers. A cross‐correlation analysis was performed to verify the results of the recession analysis and to develop a “black box” statistical model for discharge data. The correlation analysis proved the validity of springflow recession analysis for mine discharge. The recorded data length was not adequate to create a statistical model, however, but a procedure was proposed for a statistical model that could be used with large flow records. For the study site, the mine discharge was found to be sustainable for a prolonged period of time.     

INTEGRATING GEOGRAPHIC INFORMATION SYSTEMS AND MODFLOW FOR GROUND WATER RESOURCE ASSESSMENTS1

ABSTRACT: The Floridan Aquifer is the primary source of water in the coastal area of Santa Rosa County, Florida. In order to optimize well field design and analyze aquifer stress problems, the USGS MODFLOW code (McDonald and Harbaugh, 1988) is applied to develop a numerical computer model of the aquifer. The Geographical Information System (GIS) is the primary tool used in the development of the model grid, performance of the modeling procedure, and model analysis. The GIS is used in generating multiple grids in which to simulate both regional scale and local scale flow. The grid topology is recorded in geographic coordinates which facilitates geo-referencing and orientation of the grid to base maps and data coyerages. The GIS allows data transfer from various coverages to the nodes of the block centered grid where hydrogeologic information is stored as attributes to the grid coverage. From this grid coverage, pertinent information is queried within the GIS environment and used to generate the input files for the MODFLOW simulation. After MODFLOW execution, simulated heads and drawdown are imported into the grid coverage where residual error and recharge rates can be calculated. Contoured surfaces are then created for selected data sets including simulated heads, drawdown, residual error, and recharge rates. Model calibration is conducted utilizing the GIS to generate and process data sets associated with model simulations.     

AQUIFER SALINIZATION IN THE YUN‐LIN COASTAL AREA,TAIWAN1

ABSTRACT: This study analyzes possible causes of shallow ground water salinization in the coastal area of Yun‐Lin. The local hydro‐geologic setting is determined from geological drilling data and sea floor topography. Three possible causes (sea water intrusion, salt water percolation through wells, and infiltration of salty water from fish ponds) are evaluated. Chloride concentration is used as an index to measure ground water salinization. Sea water intrusion is modeled by the advective/dispersive equation, and salt water infiltration from wells and fish ponds is calculated by estimating the amount of water percolated. The determined local hydrogeologic setting suggests that the shallow aquifer may be connected to the sea water, resulting in salt water intrusion as a large amount of shallow ground water is withdrawn. The percent contributions of sea water intrusion, percolation through wells, and infiltration of water from fish ponds, to the salinization of the shallow aquifer at Ko‐Hu in the Yun‐Lin coastal area are approximately 27 percent, less than 1 percent and 73 percent, respectively. The results suggest that the vertical infiltration of salt water from fish ponds is the major cause of shallow ground water salinization in the coastal area of Yun‐Lin.     

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.     

Variable Source Area Hydrology Modeling with the Water Erosion Prediction Project Model

In nondegraded watersheds of humid climates, subsurface flow patterns determine where the soil saturates and where surface runoff is occurring. Most models necessarily use infiltration‐excess (i.e., Hortonian) runoff for predicting runoff and associated constituents because subsurface flow algorithms are not included in the model. In this article, we modify the Water Erosion Prediction Project (WEPP) model to simulate subsurface flow correctly and to predict the spatial and temporal location of saturation, the associated lateral flow and surface runoff, and the location where the water can re‐infiltrate. The modified model, called WEPP‐UI, correctly simulated the hillslope drainage data from the Coweeta Hydrologic Laboratory hillslope plot. We applied WEPP‐UI to convex, concave, and S‐shaped hillslope profiles, and found that multiple overland flow elements are needed to simulate distributed lateral flow and runoff well. Concave slopes had the greatest runoff, while convex slopes had the least. Our findings concur with observations in watersheds with saturation‐excess overland flow that most surface runoff is generated on lower concave slopes, whereas on convex slopes runoff infiltrates before reaching the stream. Since the WEPP model is capable of simulating both saturation‐excess and infiltration‐excess runoff, we expect that this model will be a powerful tool in the future for managing water quality.