|
X-DRAIN
|
---|
William J. Elliot, Project Leader David E. Hall, Computer Programmer/Analyst Susan R. Graves, Hydrologist Dayna L. Scheele, Civil Engineer U.S.D.A. Forest Service Rocky Mountain Research Station and San Dimas Technology and Development Center October 1999 |
---|
X-DRAIN is designed to estimate sediment yield from
X-DRAIN can be used to determine optimum cross drain spacing for existing or planned roads, and for developing and supporting recommendations concerning road construction, reconstruction, realignment, closure, obliteration, or mitigation efforts based on sediment yield. |
Roads have been identified as the major source of sediment in most forest watersheds due to surface erosion or mass failure. Practices to control sedimentation from roads are well known, and have been incorporated into road design and cross drain spacing guides for many years (Packer and Christensen 1977). Such guidelines, however, merely provide percentage estimates of sediment reduction at best, and their applications are limited to the specific soils and climates where they were developed. There have been numerous cases in recent years where forest planners needed sediment yield predictions from a given road, but had limited tools for estimating the amount.
One of the most common forest road conditions leading to sedimentation of streams is where a forest road experiences erosion between cross drains.
Figure 1. Relationship of road, fill slope, forest buffer, and stream for WEPP cross drain sediment yield study
All insloping, flat, rutted, or outsloping roads can be described by this model. |
Whether the cross drain is a culvert or an open drain will have minimal impact on the sediment delivery.
Outsloping roads generally have an equivalent cross drain spacing of about 10 m (Foltz 1996),a but this soon increases as wheel tracks begin to flatten the road surface, and water bars or cross drains are generally recommended on outsloping roads. Roads that cross streams or drain directly into streams have buffer width of zero.
Version 1 vs Version 2.000
Version 1 of X-DRAIN was released in June 1998 by the San Dimas Technology and Development Center as a stand-alone program as part of their Water/Road Interaction series. At the same time it was also made available to run over the Internet using Forest Service servers. X-DRAIN version 1 was based on WEPP version 95.7 results.
X-DRAIN version 2.000:
It was developed by a group of scientists from such agencies as the U.S. Department of Agriculture's Agricultural Research Service (ARS), Forest Service, and Natural Resources Conservation Service; and the U.S. Department of Interior's Bureau of Land Management and Geological Survey.
For every day being modeled, WEPP simulates the daily conditions that impact erosion-- such as the amount of vegetation canopy, the surface residue, and the soil water content. For each day that has a precipitation event, WEPP determines whether the precipitation is rain or snow, and calculates the infiltration and runoff. If there is runoff, WEPP routes the runoff over the surface, calculating erosion or deposition rates for at least 100 points on the hillslope. It then calculates the average sediment yield from the hillslope.
For a run described in figure 1, more than 400 input variables are required to describe soil, topography, and management (vegetation) in addition to a daily weather data. The model has been validated for numerous conditions including forest roads (Elliot and others 1995). Included with WEPP is a weather generator, CLIGEN, along with a database of climate statistics covering the U.S.
To determine the road surface erosion rate, a 1-m wide nonerodible element was inserted between the road and the fillslope. The fillslope was 4 m long and 50 percent steepness for all the runs.
For each run, both the erosion rate of the road and the sediment yield from the bottom of the buffer was recorded. The output files from all runs were inspected to ensure that there were no missing values, and that the results were reasonable. The results from all the runs were stored by climate in hybrid text/binary files for access by the X-DRAIN programs.
Users will generally select the station nearest to their site. In mountainous areas, however, a more distant station with a closer match in elevation may be preferred. Comparisons between two nearby stations may also be helpful to note the local sensitivity of erosion rates to climates.
Recently, there have been some errors identified in the CLIGEN database. These errors are currently in the database distributed by the Agricultural Research Service, but they have been corrected for all the climates in the X-DRAIN database.
When gravel is added to either silt loam or clay loam, there are sufficient fines in the gravel to change the soil textural category to a loam. To fully describe each set of soils for WEPP requires 48 soil parameter values. These values are listed in the Appendix. Further details describing these parameters are available in the WEPP Technical Documentation (Alberts and others 1995).
Table 1. Categories of Common Forest Soils in relation to cross-drain soils Soil type Soil Description Universal Soil Classification Clay loam Native-surface roads on shales and similar decomposing fine-grained sedimentary rock MH, CH Silt loam Ash cap native-surface road; alluvial loess native-surface road ML,CL Sandy loam Glacial outwash areas; finer-grained granitics and sand SW, SP, SM, SC Graveled loam Clay and silt loams that have been graveled GC Graveled sand Course-grained granitics, and fine-grained granitics and sandy loams that have been graveled GM
The addition of gravel has two major impacts on erosion rate. Research has shown that gravel alters the hydraulic conductivity of a soil, and it changes the water flow path length of the road (Foltz and Truebe 1995). Generally, the addition of gravel increases the porosity and increases the hydraulic conductivity of the road, decreasing the runoff (Flerchinger and Watts 1987). Within the WEPP model, however, gravel is treated as part of the rock content of the soil. WEPP reduces the conductivity and porosity of the soil with increasing rock content. To offset the impacts of gravel on conductivity, the hydraulic conductivity values shown in the appendix have been increased for gravel. In some drier climates, this increase may not have been enough to fully describe the effects of gravel. On roads where added gravel has been rolled into the base, the adjustment may be too much.
The impacts of gravel on soil properties, however, are not as important in reducing road erosion rates as the impact of gravel on flow path length. Gravel reduces the formation of ruts (Foltz and Truebe 1995). On a gravel road that is outsloping, or that has a nonerodible ditch, the flow path generally is 10 m long. On roads without gravel where flattening or rutting of wheel tracks is common, the flow path length will be the spacing between cross drains. The reduction in flow path length with the addition of gravel is likely the reason for the 80 percent reduction in erosion often credited to the addition of gravel. (Example 5 demonstrates how gravel reduces road erosion through alteration of topography.)
The topographic options for X-DRAIN are shown in Table 2. The effect of changes in topography tends to be consistent, so interpolation between values is acceptable. The effect is not, however, linear (figure 2), so extrapolation beyond these values is not recommended. The 10-m value for cross drain spacing is reasonable for outsloped roads, considering road gradient, outslope steepness, and effects of traffic (Elliot and Hall (1997), Foltz 1996). On roads with no cross drains, flow path lengths along the road surface seldom exceed 120 m unless the road is severely rutted, or the road is insloped with a long inside ditch.
Table 2. Topographic conditions in X-DRAIN database Variable Values Steepness of buffer 4, 10, 25, and 60 percent Length of buffer between road and stream 0, 10, 50, 100, and 200 m Road gradient 2, 4, 8, and 16 percent Spacing of cross drains 10, 30, 60, and 120 and 240 m
If a road is crowned, then the user will need to consider both sides of the crown and combine the results. If there are ditches which are either completely vegetated or lined with a non-erodible material, and the ditch is not eroding, then only the surface will be eroding, and the user should select the 10-m value for spacing of drains to describe the relatively short eroding path length.
The most complex input file for the WEPP model is the management file, which describes the vegetation. For the X-DRAIN program (figure 1), the road and border are described as having no vegetation, the fillslope as sufficient vegetation to give about 50 percent ground cover, and the forest buffer as a 20-year old forest with 100 percent ground cover. The values are similar to those described for the PLUME3 Template for WEPP in Elliot and Hall (1997). For other conditions, such as a burned buffer or a vegetated road, users will have to run the WEPP model.
If it is applied to skid trails, or similar temporary trails, X-DRAIN can be run to determine the erosion rate the first year. Generally, following revegetation, erosion rates rapidly decline to near zero within five years. The X-DRAIN value can be reduced by 20 percent each year following the disturbance for five years to obtain an estimate of total erosion associated with the skid trail (example 4).
On the output screen, the sediment yield results from 20 WEPP runs are presented for different road gradients and cross drain spacings. The average annual sediment yields are presented in kg or lbs from 30 years of simulation for each topographic condition. Figure 2 shows a graph of some typical results.
Figure 2. Annual sediment yield vs. cross drain spacing for different road gradients for the Deadwood Dam, Idaho climate, with a buffer length of 50 m, a sandy loam soil, a buffer slope of 25%, and a road width of 5 m.
Interpolation between values is reasonable, but because the relationships are not linear, the results should not be extrapolated. Although the values are presented to two decimal places, the user should not overly emphasize small differences. Any predicted erosion value--by any model--will be, at best, within plus or minus 50 percent of the true value. Erosion rates are highly variable, and most models predict only a single value. Replicated research has shown that observed values vary widely for seemingly identical plots, or the same plot from year to year (Elliot and others 1994; Elliot and others 1995; Tysdal and others 1999). Additional discussion of the results is presented in Morfin and others (1996).
In some cases, users may wish to estimate for environmental impact analysis the percentage of the eroded sediment that is in the 'sand' size class. WEPP predicts the size distribution of eroded sediment by dividing the sediment into sand, silt, and clay particles; small aggregates made up of clay and silt; and large aggregates made up of clay, silt, and sand. Generally, sand deposits first, and clay and small aggregates deposit last. A series of runs were made with the WEPP 97.3 model for the Eagle, CO climate, to determine the predicted sand content values as particles and in the large aggregates for a range of sediment yields (Table 3). The predicted values appear to be reasonable, but users are encouraged to compare these predictions with local observations, as size distribution predictions have not been validated for forest conditions.
Table 3. Predicted sand content as particles and in large aggregates for different sediment yield amounts for the Eagle, CO climate. Percent sand particles in large aggregates Soil Clay Loam Silt Loam Sandy Loam Graveled loam Graveled sand On site sand content 30 30 60 40 70 Sediment Yield (kg) Percent sand particles in delivered sediment 3 28 14 60 33 63 5 29 9 56 33 64 30 30 17 55 34 60 60 30 20 50 34 57 300 29 25 (52)* 33 (58)* * These values predicted for Wallace, ID climate, the Eagle, CO climate did not have sediment yields as large as 300 kg for these soils.
There are several applications for the X-DRAIN model.
Optimum cross drain spacing for existing or planned roads can be determined. Recommendations concerning road reconstruction, realignment, closure, obliteration, or mitigation can be developed and supported by X-DRAIN results.
Table 4. Adapting X-DRAIN inputs to model different road designs. Condition Cross drain application Road with flat travelled way Enter width of travelled way in width box and read output directly Insloping road with no ditch treatment and no ruts Enter width of travelled way plus inside ditch in width box Insloping road with rocked or gravel ditch and no ruts Enter width of travelled way in width box and select 10 m (30 ft) for spacing of cross drains, divide the answer by 10 m and multiply by the true road length Outsloping road without ruts Enter width of travelled way in width box and select 10 m (30 ft) for spacing of cross drains, divide the answer by 10 m and multiply by the true road length (example 2) Outsloping road with ruts Enter width of travelled way. Read the results for the observed spacing of cross drains Bladed and compacted skid trail Select appropriate native surface soil and appropriate topographic variables for first year erosion. Subsequent years will decline rapidly as vegetation is reestablished on the skid trail, to near zero by year 5 (example 4). Added gravel Select graveled soil type, select cross drain spacing of 10 m (example 5). More complex conditions Run the WEPP model for the specific conditions
There has been no direct validation of the X-DRAIN program, nor has anyone collected data on the amount of sediment that enters a stream that was detached from a road and transported across a forested buffer. There have been field observations on erosion rates from road surfaces and road prisms, and on lengths of sediment plumes. Some of these observations are presented in Table 5, along with the sediment yields and sediment plume lengths predicted by the X-DRAIN program.
Table 5. Erosion rates and sediment plume lengths below for roads observed and predicted by X-DRAIN for 60-m long, 4 percent gradient road, 4 m wide. Site and Ecoregion Erosion Rate
(t/ha/yr)Source Observed at Zena Creek, ID 18 Megahan and Kidd 1972a. Included entire new road prism. Predicted for Zena Creek, ID 13 Deadwood Dam, ID climate Observed bare and graveled road in Southern Appalachians 11 to 150 Swift 1984a and 1984b Predicted for Southern Appalachians 42 - 51 Cullowhee, NC climate. Graveled sand and sandy loam soil types. Observed for Alum Creek, AR 15 - 76 Beasley and others 1984. Vegetated fill slope and ditch. Predicted for Alum Creek, AR 53 - 57 Clarksville, AR climate. Graveled sand and sandy loam soil types. Observed bare and gravelled road for Fernow NF, WV 13.5 - 118 Kochenderfer and Helvey 1987 Predicted for Fernow NF, WV 26 - 69 Clarksburg, WV climate. Graveled sand and sandy loam soil types. Sediment Plume Length
(m)
Observed on Silver Creek Watershed, ID
Cross drain
Below fill11 - 183
(mean = 50)
0.4 - 66
(mean = 3.8)Ketcheson and Megahan 1996 Observed in Nez Perce NF, Central ID
below culverts80 percent less than a mean of 24 Wasniewski 1994 Predicted for Silvercreek
Cross drain
Below fill< 50 m
< 50 mDeadwood Dam, ID climate. Sediment carried past 10 m, but not past 50 m. Sandy loam soil. Predicted for Wine Springs Watershed, NC L = 5.1 + 0.00197 M
= 7.0 mMcNulty and others 1995. L is length in meters, M is predicted sediment yield from road. Length with a X-DRAIN predicted yield of 1013 kg. Filter strip widths in SE U.S. 13 + 0.426 x percent slope = 24 m Swift 1986. Assumed slope of 25 percent. Predicted for SE U.S. 50 m For Cullowhee, NC climate. Graveled sand soil. 700 of 1013 kg deposited within 50 m Observed in Tuskegee NF, AL 50 - 60 Grace 1998. Occurred Aug '97 - Jan '98 Predicted for Tuskegee NF, AL 50 - 100 Birmingham, AL climate. Graveled sand soil. 10 percent buffer gradient. 1200 out of 1400 kg deposited within 50 m.
The presence of a distinct sediment plume does not necessarily mean that sediment does not travel beyond the observed plume. Numerous runs with the WEPP model have shown that a deposition plume can develop in a buffer zone during small runoff events, but larger runoff events route sediment over the zone of deposition, across the buffer, delivering large amounts of sediment to the stream system.
If you have X-DRAIN installed on your isolated PC, start your web server and run X-DRAIN locally.
If you are running stand-alone, please register with us so we can inform you of updates to X-DRAIN. Also, check either of the two Forest Service servers above frequently to keep current.
Field research shows that the range of sedimentation amounts observed will vary by at least 30 percent from the mean. Minimum observed values are frequently less than half the maximum observed values (Elliot and others 1994; Elliot and others 1995; Tysdal and others 1999). Validation work has shown values predicted by WEPP generally fall within the range of observed values. Users are encouraged to avoid placing too much emphasis on small differences between values. For conditions not modeled, the output relationships appear to be continuous, and interpolation between results appears to be valid. It is not advisable to extrapolate beyond the values presented, as the relationships are not linear.
In a limited analysis of the sensitivity of sediment yield to the various input factors, Morfin and others (1996) found sediment yield was particularly sensitive to cross drain spacing, road gradient, and buffer length. It was less sensitive to buffer slopes above 25 percent. In their study, sediment yield was sensitive to both climate and soil type.
This study assumed runoff water followed the road from one cross drain to the next. This template can be applied to a variety of conditions and provides a reasonable estimate of sediment yield. If the site template is not adequate to describe the site conditions, then site-specific runs can be made with the WEPP model with the aid of the templates developed by Elliot and Hall (1997).
If the road drains directly into channels, then the cross drain template is not valid, and the WEPP watershed version may be the more appropriate modeling tool (Tysdal and others 1999).
To run the most current version of X-DRAIN and the other FS WEPP interfaces to WEPP, go to the Forest Service Intranet FS WEPP or the Forest Service Internet FS WEPP site.
If you have a local version, start the Forest Service server, fsws, then start your local FS WEPP interfaces.
If X-DRAIN is started from the FS WEPP screen, you can specify whether you want to use metric or English units in X-DRAIN.
The unit system is selected by clicking the hole next to the desired unit description on the FS WEPP screen.
Climate
A scroll bar located within the climate list allows the user to scroll among the various available climates.
Clicking on the "Climate" hyperlink located directly above the climate list produces a descriptive table for all of the climates. This table is designed to assist the user in selecting a climate most appropriate to the location being studied. In some cases, the most appropriate climate may not be the one geographically closest to the location being simulated. Select the most appropriate climate based on the latitude, longitude, elevation, and annual precipitation.
Soil Type
Five soil types are available within X-DRAIN.
Buffer Topography
The buffer is assumed to be a forested area between the outlet from a road cross drain or waterbar and an ephemeral or perennial channel. The buffer slope steepness and horizontal length nearest to those of the site conditions being simulated are selected. Users may wish to look at several combinations of buffer length and steepness, and then interpolate between values to determine a more exact topography. Extrapolation beyond the largest or smallest values is discouraged.
Road Width
You may enter any numeric value for road width between 1 and 30 meters (3 and 100 feet). See table 3 for guidance on entries for road width; it may be only the traveled-way width or the sum of the traveled-way and ditch width, depending upon what is being modeled.
The rule of thumb is: if it is eroding, include it as part of the road width.
Run
Once the required inputs have been selected by the user, click the "Run" button to display a table of sediment yield values on the X-DRAIN display screen.
Examples of interpreting the sediment yield values provided by X-DRAIN are given in the report body.
Saving and printing Results
Use your browser's "File�Save-As" menu option to save the results as a browser-readable file, or print the page with "File�Print...".
Copy to Clipboard
Sediment yield information may be copied to the Windows clipboard by selecting the text with the mouse, and using the standard CTRL-C or your browser's Edit�Copy command.
From the clipboard, sediment yield information may be pasted directly into other Windows applications such as word processors or spreadsheets for further processing.
Exiting X-DRAIN
If you are running a local copy of X-DRAIN, return to the FS WEPP screen and click on "Close down FS WEPP interfaces" to quit X-DRAIN and shut down the FS WEPP server.
If you are running X-DRAIN from a remote server, there is nothing particular to do to quit X-DRAIN. Simply shut down your browser, or travel else where with it.
Soil Type Element and Property Units Clay
LoamSilt
LoamSandy
LoamGraveled
LoamGraveled
SandRoad Surface Albedo of the bare soil Fraction 0.6 0.6 0.6 0.6 0.6 Initial saturation of the soil m/m 0.5 0.5 0.5 0.5 0.5 Baseline interrill erodibility (ki) kg*s/m4 1.50E+06 2.00E+06 2.00E+06 1.50E+06 2.00E+06 Baseline rill erodibility (kr) s/m 0.0002 0.0003 0.0004 0.0002 0.0004 Baseline critical shear N/m2 0.04 0.04 0.04 0.04 0.04 Hydraulic conductivity mm/h 0.1 0.27 3.8 8.4 10.2 Depth of soil layer mm 200 200 200 200 200 Percentage of sand % 30 30 70 35 70 Percentage of clay % 30 15 5 30 5 Organic matter (by volume) % 0.01 0.01 0.01 0.01 0.01 Cation exchange capacity Meq/100 g 24 12 4 24 4 Rock fragments (by volume) % 20 10 20 65 65 Border Albedo of the bare soil Fraction 0.6 0.6 0.6 0.6 0.6 Initial saturation of the soil m/m 0.5 0.5 0.5 0.5 0.5 Baseline interrill erodibility (ki) kg*s/m4 1 1 1 1 1 Baseline rill erodibility (kr) s/m 0.0001 0.0001 0.0001 0.0001 0.0001 Baseline critical shear N/m2 40 40 40 40 40 Hydraulic conductivity mm/h 7.3 9.9 15.3 8.1 15.3 Depth of soil layer mm 300 300 300 300 300 Percentage of sand % 30 30 60 30 60 Percentage of clay % 30 15 5 30 5 Organic matter (by volume) % 0.01 0.01 0.01 0.01 0.01 Cation exchange capacity Meq/100 g 24 12 4 24 4 Rock fragments (by volume) % 20 10 20 30 30 Fillslope Albedo of the bare soil Fraction 0.12 0.12 0.12 0.12 0.12 Initial saturation of the soil m/m 0.45 0.45 0.45 0.45 0.45 Baseline interrill erodibility (ki) kg*s/m4 1.50E+06 2.00E+06 2.00E+06 1.50E+06 2.00E+06 Baseline rill erodibility (kr) s/m 0.0002 0.0003 0.0004 0.0002 0.0004 Baseline critical shear N/m2 2 2 2 2 2 Hydraulic conductivity mm/h 6.3 8.9 12.5 7.1 14.3 Depth of soil layer mm 300 300 300 300 300 Percentage of sand % 30 30 60 30 60 Percentage of clay % 30 15 5 30 5 Organic matter (by volume) % 4 4 4 4 4 Cation exchange capacity Meq/100 g 26 13 4 26 4 Rock fragments (by volume) % 20 10 20 30 30 Forest Buffer Soil Albedo of the bare soil Fraction 0.02 0.02 0.02 0.02 0.02 Initial saturation of the soil m/m 0.4 0.4 0.4 0.4 0.4 Baseline interrill erodibility (ki) kg*s/m4 1.00E+04 1.20E+05 1.30E+05 1.00E+04 1.30E+05 Baseline rill erodibility (kr) s/m 0.0002 0.0003 0.0004 0.0002 0.0004 Baseline critical shear N/m2 2 2 2 2 2 Hydraulic conductivity mm/h 20 25 30 20 30 Depth of soil layer mm 300 300 300 300 300 Percentage of sand % 30 30 60 30 60 Percentage of clay % 30 15 5 30 5 Organic matter (by volume) % 8 8 8 8 8 Cation exchange capacity Meq/100 g 27 14 5 27 5 Rock fragments (by volume) % 20 10 20 20 20 For further details of soil properties see Alberts and others 1995.
Useful Conversions
Multiply by to get mm (millimeters) 0.0394 in. (inches) m (meters) 39.4 in. (inches) m (meters) 3.28 ft (feet) m2 (square meters) 10.8 ft2 (square feet) kg (kilograms) 2.2 lbs (pounds mass) t (metric tonnes) 1,000 kg (kilograms) t (metric tonnes) 1.1 short tons short tons 2,000 lbs (pounds mass) t/ha 0.445 short tons/acre