Evaluation of C and P Factors in Universal Soil Loss Equation on Trapping Sediment: Case Study of Santubong River ()
1. Introduction
Buffer zone is the vegetation area including trees, grasses and bushes adjacent to streams, rivers, creeks or wetlands [1]. The aim of Buffer zone is to remove sediment and other pollutants from surface water runoff through filtration, deposition and infiltration processes [1]. Buffer zone is usually introduced to protect water bodies by slowing and reducing the surface runoff, and allows it to be absorbed by the ground to prevent flooding, and provides habitat for wildlife and enhances the aesthetics of the surroundings.
An appropriate size of buffer zone is important because under-sized buffer zone is unable to provide adequate protection for water bodies. In contrast, over-sized buffer zone might result in economic losses. In the past, Best Management Practices (BMPs) such as grass buffer strips, vegetative filter strips, riparian buffer zones, grass waterways, have been suggested as potential controls to help reduce erosion and offsite transport of sediments. Generally, vegetated buffers show a positive effect on reducing the transfer of sediments, nutrients and pesticides to surface waters [2,3].
Sediment yield over an area is governed by the Universal Soil Loss Equation (USLE). USLE was firstly developed by Walter H. Wischmeier in 1958, United States. It is the most comprehensive technique available to predict the long term average annual rate of erosion on a field slope. This erosion model was created for use in selected cropping and management systems, but is also applicable to nonagricultural conditions such as construction sites.
In the case of riparian zone, USLE is applied as sediment yield before reaching a vegetated buffer. The equation for USLE is presented in Equation (1).
(1)
whereA = Average annual soil loss per unit area;
K = Soil erodibility factor;
R = Rainfall and runoff erodibility index;
C = Crop/vegetation and management factor;
P = support practice factor;
LS = slope length-gradient factor.
In previous studies, the parameters of soil erodibility factor (K), rainfall and runoff erodibility index (R), and slope length-gradient factor (LS) are extensively studied, for instances Meester and Jungerius (1978), Loch and Rosewell, (1992), had investigated on K factor [4,5], Yu (1999), Mitasova (2002), Ryan (2012) and Mcusburger et al. (2012) had explored on R factor [6-9] and LS factor was studied byJose and Martin (2010) and Liu et al. (2011) [10,11].
However, two other parameters named ascrop/vegetation and management factor (C) and support practice factor (P), are not fully studied yet. Due to the rising awareness of riparian conservation, it is the initiation of this paper to explore on vegetative cover and its conservation or support practice factors to trap the sediment. This study employs Buffer Zone Calculator as a tool to determine the sediment removal efficiency for different C and P.
2. Buffer Zone Calculator
Buffer zone calculator was built up from a series of formulas. Basically, the formulas can be divided into two parts, the determination of buffer width and the outcome assessment of percentage of sediment reduction.
2.1. Determination of Buffer Width
Determination of buffer width is calculated using three following equations:
1) Soil loss or sediment produced on site, calculated by Universal Soil Loss Equation (USLE) [12] as presented in Equation (1).
2) Annual Average total solid concentration is the concentration of sediment flows into the river, calculated using Equation (2).
(2)
where AATSC = Annual Average Total Solid Concentration;
AASLD = Average Annual Solid Loss from Development;
AAR = Average Annual Runoff;
(3)
(4)
3) The required buffer width takes into account the factor of slope [2] and the equation is:
(5)
2.2. Outcome Assessment
The outcome assessment including Total Suspended Solid (TSS) retained, Outfall TSS and percentage of sediment reduction, are calculated using following equations.
1) Sediment retained by buffer is a function of concentration of TSS retained by vegetation divided by slope of buffer as presented in Equation (6). When the slope is steep, water flows at higher speed through the buffer. Thus, it required wider vegetative buffer width and acquired more time to slowdown the water flow for vegetation in buffer zone to trap sediment [13].
(6)
where CRBV = Concentration Retained by Buffer Vegetation (mg/l);
RRBW = Required Riparian Buffer Width (m);
2) Sediment released from buffer is the difference between TSS concentration inflow and TSS retained by buffer, as presented in Equation (7).
(7)
where AATSSC = Annual Average TSS Concentration (mg/l);
TSSR = TSS Retained (mg/l).
3) Percentage of sediment reduction is the removal rate of the sediment concentration in the buffer zone [14], as represented in Equations (8) and (9) respectively.
(8)
(9)
3. Study Area
The chosen study area is located at Santubong, which is well known for its fascinating wildlife in Sarawak. As Santubong river flows into South China Sea, Irrawaddy dolph in is always spotted at Santubong river mouth. Along Santubong river, rare proboscis monkeys are always spotted within the mangrove swamps. During night fall, fireflies are flying around the branches of the mangroves and crocodiles are often seen on the mud banks.
Santubong area is categorized as rural catchment. The riparian buffer zone along the Santubong river is presented in Figure 1. The red lines indicate part of the horizontal alignment of Santubong river and Santubong Bridge is circled in yellow.
A site investigation was carried out at the early stage of this study. It was found that the length of buffer zone is about 2500 m and the hill slope length is approximately 200 - 250 m. This area currently is covered by forest, and it is expected to turn into agricultural hub in nearest future. For both sides of Santubong river, it was covered with the buffer zone with the average width of 50 m each site. Both hill slope and buffer gradients are averagely 5%.
4. Methodology
USLE consists of five major factors namely R, K, LS, C and P for calculating the soil loss. In recent years, the parameters of K, R, and LS are extensively studied. However, parameters of C and P are not fully studied yet. Therefore, there is an urgent need to explore and investigate the impact of P and C factors to trap the sediment. In this study, the impact of C and P factors are investigated using Buffer Zone Calculator.
Factor C is used to determine the relative effectiveness of soil and crop management systems in term of preventing soil loss. C factor is a ratio comparing the soil loss from land under a specific crop and management system to the corresponding loss from continuously fallow and tilled land. The C Factor can be determined by selecting the crop type and tillage method that corresponds to the field and then multiplying these factors together. This generalized C factor provides relative numbers for different cropping and tillage systems. Thereby it is able to weigh the merits of each system.
There are 17 alternatives of future land uses available in Buffer Zone Calculator. These 17 different future land uses are forest with undergrowth, forest with no undergrowth, logging 10^6, logging 20, logging 30, mixed dipt. forest, new shifting agriculture, old shifting agriculture, settlement/cleared land, cultivated grass, oil palm plantation, rubber plantation, beans plantation, cocoa plantation, coffee plantation, tea plantation and paddy plantation. The C values for different types of crop investigated are presented in Table 1.
For investigating the impact of C values to the outfall and % reduction of TSS, the selection of C values for different land uses in Buffer Zone Calculator is changing, while 9 other parameters obtained from site investigation are remained constant. The values for these 9 parameters are:
1) soil types = Rjn − Rajang;
Figure 1. Santubong River with buffer vegetation grown alongside the river. (a) Plan view of Santubong area; (b) View of Santubong river from bridge.
Table 1. C values for different types of crop investigated.
2) mean annual precipitation = 3000 mm;
3) length of hill slope = 215 m;
4) length of buffer = 2300 m;
5) types of buffer vegetation = Forest;
6) gradient of hill slope = 7.5%;
7) gradient of buffer future soil conservation practices = 5%;
8) buffer zone width = 50 m;
9) Soil support practice = none.
Meanwhile, P factor is the support practice factor. It reflects the effects of practices that will reduce the amount and rate of the water runoff and thus reduce the amount of erosion. The P factor represents the ratio of soil loss by a support practice to that of straight-row farming up and down the slope. The most commonly used supporting cropland practices are cross slope cultivation, contour farming and strip cropping.
Table 2 presents P values for different investigated support practice factors. The conservation and support practices available in Buffer Zone Calculator comprised of none, contouring, contouring strip-cropping and terracing. The “none” means that the existing buffer zone at the study area will not undergo any future development and conservation practice. The impact of P factor to the sediment trapping capacity of buffer zone was conducted by selecting one of the practices, while the other 9 factors are remained constant as listed below.