Shuhui Dun

Effects of DEM source and resolution on WEPP hydrologic and erosion simulation: a case study of two forest watersheds in Northern Idaho.

The recent modification of the Water Erosion Prediction Project (WEPP) model has improved its applicability to hydrology and erosion modeling in forest watersheds. To generate reliable topographic and hydrologic inputs for the WEPP model, carefully selecting digital elevation models (DEMs) with appropriate resolution and accuracy is essential because topography is a major factor controlling water erosion. Light detection and ranging (LIDAR) provides an alternative technology to photogrammetry for generating fine-resolution and high-quality DEMs. In this study, WEPP (v2006.201) was applied to hydrological and erosion simulation for two small forest watersheds in northern Idaho. Data on stream flow and total suspended solids (TSS) in these watersheds were collected and processed. A total of six DEMs from the National Elevation Dataset (NED), Shuttle Radar Topography Mission (SRTM), and LIDAR at three resolutions (30 m, 10 m, and 4 m) were obtained and used to calculate topographic parameters as inputs to the WEPP model. WEPP-simulated hydrologic and erosion results using the six DEMs were contrasted and then compared with field observations. For the study watersheds, DEMs with different resolutions and sources generated varied topographic and hydrologic attributes, which in turn led to significantly different erosion simulations. WEPP v2006.201 using the 10 m LIDAR DEM (vs. using other DEMs) produced a total amount of as well as seasonal patterns of watershed discharge and sediment yield that were closest to field observations.

Development of a GIS Interface for WEPP Model Application to Great Lakes Forested Watersheds

This presentation will highlight efforts on development of a new online WEPP GIS interface, targeted toward application in forested regions bordering the Great Lakes. The key components and algorithms of the online GIS system will be outlined. The general procedures used to provide input to the WEPP model and to display model output will be demonstrated.

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.

Modifying WEPP to improve streamflow simulation in a Pacific Northwest watershed

Abstract. The assessment of water yield from hillslopes into streams is critical in managing water supply and aquatic habitat. Streamflow is typically composed of surface runoff, subsurface lateral flow, and groundwater baseflow; baseflow sustains the stream during the dry season. The Water Erosion Prediction Project (WEPP) model simulates surface runoff, subsurface lateral flow, soil water, and deep percolation. However, to adequately simulate hydrologic conditions with significant quantities of groundwater flow into streams, a baseflow component for WEPP is needed. The objectives of this study were (1) to simulate streamflow in the Priest River Experimental Forest in the U.S. Pacific Northwest using the WEPP model and a baseflow routine, and (2) to compare the performance of the WEPP model with and without including the baseflow using observed streamflow data. The baseflow was determined using a linear reservoir model. The WEPP-simulated and observed streamflows were in reasonable agreement when baseflow was considered, with an overall Nash-Sutcliffe efficiency (NSE) of 0.67 and deviation of runoff volume (D v ) of 7%. In contrast, the WEPP simulations without including baseflow resulted in an overall NSE of 0.57 and D v of 47%. On average, the simulated baseflow accounted for 43% of the streamflow and 12% of precipitation annually. Integration of WEPP with a baseflow routine improved the model’s applicability to watersheds where groundwater contributes to streamflow.

Residue Management Impacts on Field-Scale Snow Distribution and Soil Water Storage

Spatial variation of soil water affects crop performance, fertilizer use efficiencies, and other important economic and environmental factors. Soil water storage could be increased and field variability reduced by residue management practices such as no tillage (NT), as surface residues can retain more snow, enhance water infiltration, and reduce evaporation as compared to conventional tillage (CT). Our objectives were to evaluate the residue effects on snow distribution and the spatial variation of soil water storage for two adjacent fields near Pullman, Washington: one under NT, and the other under CT. The fields were surveyed during the winter and spring of 2007-2008 to assess topographic variations in snow depth, snow water equivalent (SWE), and soil water storage. Standing stubble under NT retained 10 to 20 cm more snow on ridge tops and steeply sloped ground than CT, and the snowpack was distributed more evenly with less spatial variation. SWE followed the same pattern of larger spatial variation in CT than in NT. Soil water (0 to 1.5 m) in the spring was lowest for ridge tops and highest in valleys for NT and CT. Under NT, however, soil water varied less across different field locations than under CT, and overall water storage was 60, 29, and 13 mm more for NT than CT on ridge top, south slope, and valley locations, respectively. Although many factors can contribute to the spatial variation of soil water, standing wheat residue under NT retained more snow on ridge tops and steeply sloped areas, reduced soil water spatial variation, and increased soil water recharge.

IMPROVING FROST-SIMULATION SUBROUTINES OF THE WATER EROSION PREDICTION PROJECT (WEPP) MODEL

Erosion models play an important role in assessing the influence of human activities on the environment. For cold areas, adequate frost simulation is crucial for predicting surface runoff and water erosion. The Water Erosion Prediction Project (WEPP) model is a physically based erosion prediction software program developed by the USDA. One of the major components of WEPP is the simulation of winter processes, which include snow accumulation and melt as well as soil freeze and thaw. WEPP is successfully used in the evaluation of important natural resource issues throughout the U.S. and in a number of other countries. However, previous studies revealed problems in the winter component of the WEPP model, especially the routine for frost simulation. The main purpose of this study was to improve the WEPP model (v2006.5) by changing the soil profile discretization and computation of key thermal and hydraulic parameters in the frost simulation routines so that the model can adequately simulate soil freeze-thaw and winter runoff and erosion. WEPP v2006.5 and the modified version (v2010.1) were applied to experimental plots in Pullman, Washington, and Morris, Minnesota. The simulated snow and frost depths as well as runoff and sediment yield were contrasted and compared with field observations; the results from v2010.1 showed substantial improvement compared to those from v2006.5.

WEPP simulations of dryland cropping systems in small drainages of northeastern Oregon

Computer simulation models are essential tools for evaluating soil erosion potential over large areas of cropland. Small-plot and field-scale evaluations are commonly conducted for federal farm program compliance, but producers are now faced with off-farm water quality concerns. Predicting the potential contribution of upland sediment is of interest to producers and state and federal agencies. The purpose of this study was to evaluate the applicability of the Water Erosion Prediction Project (WEPP) model for quantifying hydrological and erosion processes in the semiarid croplands of the Columbia Plateau. Two headwater drainages managed using conventional inversion tillage (CT) or no-tillage (NT) management techniques were monitored from 2001 through 2007 in the dryland cropping region of northeastern Oregon. The WEPP model was parameterized primarily from field data, including management and weather data. Crop parameters (above-ground biomass and crop yield), water balance components (volumetric soil water, evapotranspiration [ET], and surface runoff), and soil loss were observed and subsequently used to evaluate WEPP simulations. This detailed dataset allowed for a unique opportunity to evaluate not only the WEPP routines for runoff and erosion but also the routine for crop growth, which greatly influences the erodibility and hydraulic conductivity of top soil layers. Graphical and goodness-of-fit analyses indicate that WEPP generated satisfactory estimates for volumetric soil water and crop yields in NT and CT and above-ground biomass production in NT. Gross patterns of ET simulated by WEPP were compatible with those determined using observed precipitation and soil water data. Observed annual runoff and erosion values from both drainages were low (NT: 0.1 mm [0.004 in], 2.5 kg ha−1 [0.001 tn ac−1]; CT: 0.9 mm [0.04 in], 72.0 kg ha−1 [0.03 tn ac−1]). On average only 0.3% and 0.03% of total precipitation left the catchment as runoff during the six-year study period for CT and NT, respectively. No runoff was predicted by WEPP when default input values for a Walla Walla silt loam soil were used in the model. Simulated runoff and erosion agreed well with field observations after the effective surface hydraulic conductivity Keff and rill erodibility Kr were calibrated. With minimal calibration, the WEPP model was able to successfully represent the hydrology, sediment transport, and crop growth for CT and NT cropping systems in northeastern Oregon during years of below normal precipitation, mild weather, and little runoff.

Seasonal Change of WEPP Erodibility Parameters for Two Fallow Plots on a Palouse Silt Loam

Abstract. In cold regions, frozen soil has a significant influence on runoff and water erosion. In the U.S. Inland Pacific Northwest, major erosion events typically occur during winter as frozen soil thaws and exhibits low cohesion. Previous applications of the WEPP (Water Erosion Prediction Project) model to a continuous bare tilled fallow (CBF) runoff plot at the Palouse Conservation Field Station (PCFS) in southeastern Washington State showed that WEPP reproduced the occurrence of the large erosion events, but the amount of sediment yield was either under- or overestimated. The inability of WEPP to reproduce the magnitude of field-observed erosion events at the PCFS suggests the need for an examination of the dynamic changes in soil erosion properties and for improving the representation of such dynamics. The objective of this study was to evaluate the seasonal changes of rill erosion parameters for two CBF runoff plots at the PCFS. Field-observed runoff and erosion events during 1984-1990 were used to estimate WEPP hydraulic and erosion parameters, including soil effective hydraulic conductivity, critical shear stress (I„ c ), and rill erodibility (K r ). The parameters for each event were best-fitted using WEPP single-event simulations to reproduce the observed runoff and sediment yield on both plots. The results suggest that the adjustment factors for I„ c and K r of frozen and thawing soils in the WEPP model could be modified to better describe the changes of these parameters in winter.

Effects of DEM resolution on forest hydrologic and erosion prediction using WEPP

The recent modification of WEPP (Water Erosion Prediction Project) has improved the original model’s applicability to hydrology and erosion modeling in forest watersheds. To generate reliable topographic and hydrologic inputs for the WEPP model, carefully selecting Digital Elevation Models (DEMs) with appropriate resolution and accuracy is essential because topography is a major factor controlling water erosion. LIght Detection and Ranging (LIDAR), a new remote sensing technology, provides an alternative for generating fine and high-quality DEMs. This study applies WEPP (v2006.201) for hydrological and erosion simulation under forest conditions and evaluates the effects of DEM resolution and accuracy on watershed hydrology and water erosion prediction at a watershed scale. Stream flow and total suspended solids (TSS) in two small forest watersheds located in northern Idaho were collected and processed. A total of six DEMs from three sources (NED, SRTM, and LIDAR) at three resolutions (30 m, 10 m, and 4 m) were obtained and used to calculate topographic parameters as inputs to the WEPP model. WEPP-simulated hydrologic and erosion results using the six DEMs were compared with the field-observed data. For both study watersheds, DEMs with different resolutions and sources generated varied topographic and hydrologic attributes, which in turn led to significantly different erosion predictions by WEPP.

Applying Online WEPP to Assess Forest Watershed Hydrology

Abstract. A new version of the online Water Erosion Prediction Project (WEPP) GIS interface has been developed to assist in evaluating sediment sources associated with forests and forest management within the Great Lakes basin. WEPP watershed structure and topographical inputs for each watershed element are generated from the USGS 30 m National Elevation Dataset (NED), soil inputs are automatically retrieved from the USDA-NRCS SSURGO database, and land use and management inputs are selected from the WEPP database based on the USGS National Land Cover Database 2001 (NLCD2001). Additionally, ground cover and soil properties of the WEPP management and soil input files can be customized to represent site-specific conditions. Daily climate inputs are generated from long-term climate parameters using CLIGEN, a stochastic climate generator embedded in the online interface. Alternatively, a registered user can upload and use observed daily climate data for online WEPP simulation. Long-term observational data, including runoff and water chemistry, from two mature forest watersheds of the Fernow Experimental Forest in West Virginia were used to assess the online WEPP GIS interface. Online WEPP simulations were carried out using both observed and CLIGEN-generated climate inputs, and model performance was examined by comparing simulated and observed runoff and simulated and estimated (from measured water chemistry data) sediment yield. The online WEPP reasonably simulated average annual runoff and the annual maximum runoff series for both watersheds, but overpredicted sediment yield for the annual average and annual maximums. The online WEPP simulation results accurately reflected the differences between the two watersheds in their hydrological characteristics. The online WEPP GIS interface is a user-friendly, web-based computer package that can be used by scientists, researchers, and practitioners as a cost-effective simulation tool for watershed management.