This research is aimed at predicting erosion and sedimentation in the Ghareh Aghach Basin in the central part of Iran using Erosion Potential Model (EPM) models incorporated into Geographic Information System (GIS) software. This basin has an area of about 8955 hectares. The region has a range of vegetation, geological, soil texture and land use types. The basin was subdivided into 5 sub-basins. Data required for this study were collected in part through published reports, whilst the remaining was derived by field surveys. Necessary maps in EPM models were prepared in Autocad-2006 medium and were transported to IILWIS, after some revision. After constructing topologies for all polygons, we entered weightings for all layers within the Arc-View software. Combinations of all layers were managed thereafter. Coefficient of each factor was determined, and erosion intensity coefficient (Z) was calculated. Four layers for EPM model were combined to develop the final layer of erosion and sedimentation. The results of the EPM model for homogenous and uniform sampling units showed that 0.19% (16.7 ha) of the total watershed area were classified as class I of erosion category with very low sedimentation and 15.1% (1352 ha) were classified as class II of erosion category with low sedimentation and 41.3% (3699 ha) were classified as class III of erosion category with medium sedimentation and 13.2% (1175 ha) were classified as class IV of erosion category with high sedimentation and finally 30.2% (2711 ha) were classified at class V of erosion category with very high sedimentation, respectively. The result of comparing erosion and sediment values using an EPM model with measured values showed that no significant difference was observed between the estimated and measured values (P < 0.05).
Introduction
The gradual soil salinization along with deterioration of rangeland vegetation covers has caused accelerating soil erosion and sedimentation in water reservoirs.1,2 The United Nations in its development plan3 has reported that at present rate of soil erosion in Iran is about 20 tones/ha, which has increased by 10 tones/ha compared to the last decade.3 There are not enough sediment measurement stations in most watersheds in the country which makes it more difficult to provide specific models with details of local watershed characteristics. One of the most important problems with empirical models of soil erosion is their lack of accuracy in processing the huge amount of data which should be digitalized by GIS systems and analyzed by mathematical models. EPM is an empirical model for estimating the quantity and quality of sediment. In fact, quantifying and digitalizing the sediment data is an important breakthrough in developing sediment assessment models.4 This problem could be partially solved by estimating models.5 Since soil erosion is a product of few different interacting factors, there is no simple model that can assess all the contributing elements simultaneously.6
Lack of information for preparing erosion maps for quantitative and qualitative evaluation of sedimentation rates is a major need in the watershed management in Iran. Finding a suitable model for evaluation of erosion and sedimentation processes is an important area of research in soil conservation. The goal of this research was to apply an EPM model and a GIS system (to reduce the error level to the minimum) to have the most possible accurate assessment of the soil erosion and sedimentation in the study area and make the best applicable suggestions to control it.
Material and Methods
Modeled area
The area of Ghareh Aghach watershed is about 8954 hectares which is located between 51° 45′ 53″ and 51° 34′ 54″ eastern latitudes and 31° 30′ 28″ and 31° 26′ 19″ northern altitudes. The maximum and minimum heights of the area are 3810 and 2630 m from the sea level. The mean altitude of the watershed is 2936 meters from the sea level. The watershed area is divided in five major hydrologic sections according to its topographic features.15 Long-term mean annual precipitation and temperature are 358 mm and 10.5 °C respectively, which classifies the site in semi-arid climatic conditions according to Ambergie (Q = 40.7) and Demartin (Ia = 17.4) categories. There are about 17 different vegetation types. The natural vegetation is mostly Astragalus spp., Agropyron trichophoum, Bromus tomentalus, forbs, and annual perennial grasses. Four major landscapes of mountainous lands, hilly lands, plateau and gravel lands consisting of nine minor land units.7
Model description
EMP model is the other package which has been used to estimate soil erosion in Iran. However, since fewer environmental factors are considered in this model, its accuracy in sediment estimation is less than in other models.8,9 This model was created based on erosion measurements during 40 years in previous Yugoslavia and for the first time was introduced in River Stream International Conference.10 In the use of this model, the first step was to develop a topographic map in 1:25000 scale of the study area. This was done by digitizing the data using the GIS ILWIS program11 to provide a contour map. Using the contour map, the slope and aspect maps as well as the digital evaluation model (DEM) of the site were prepared. An Erosion feature map was obtained using aerial photos of 1:40000 scale. In the second step a land use map was obtained from supervised classification of the satellite images (ETM+ bands 1-7) and data obtained from the field visits that were used to provide information on the local topography, climatic factors, geology, pedology and land use data of the study area. Precipitation and temperature maps were obtained from the climatic data measured by gages placed in the watershed. By integration of geology map, erosion faces map and slope maps, 21 homogenous uniform sampling units were identified to apply EPM model (Fig. 1).
Model information layers comprised four effective factors in erosion and sedimentation in the watershed site. These were obtained by digitalizing and classifying the basic information data in the GIS program. In frequent field visits all the necessary data for EPM model were collected, checked and improved in 21 homogenous and uniform sampling units. Following of finalizing the information layers, all of them were integrated to create the final erosion map (Fig. 2).
Therefore, sediment estimation in this model is based on four factors consisting: ϕ: Erosion coefficient of watershed, X: Land use Coefficient, Y: susceptibility of rock and soil to erosion, I: mean watershed slope each one of these factors will receive a proper value according to its contribution to erosion process (Table 1).
Explination of Point Procedure by Quadruplet Factors
Erosion coefficient
The accurate way for determinding of erosion coefficient is providing geomorphologic map of watershed up to facieses or sampling units, there for the amount of erosion coefficient can be obtained base on Table 2, for each sampling units.
Table 1
The contributing factors in EMP model to estimate the soil erosion.
Mean Watershed slope
By using topography map of watershed slope map was provide, then the average weight slope for each sampling units were calculate.
Coefficient Susceptibility Erosion of Rock and soil
The soil erosion sensitivity coefficient map was obtained based on point to geomorphologic and peddological map for each sampling units according to Table 3.
Table 2
Erossion coefficient point in facieses or sampling units.
Table 3
The soil erosion sensitivity coefficient.
The Relationship between Erosion and Sedimentation in Emp Model
This model is able to measure erosion, sediment carrying capacity as well as a primary as estimation of sediments behind the reservoirs.
The coefficient of erosion intensity (Z) is calculated by the following equation in this model:
Where:
Y: rock and soil susceptibility coefficient, Xa: land use coefficient, ϕ: Erosion coefficient of watershed, I: average slope (percent). Tabular values are available for each coefficient. Table 5 shows classification of Z values. According to Z values the erosion intensity map can be prepared. Table 5. Classification of Z coefficient values
Table 4
The land use coefficient.
Specific erosion, sediment yield coefficient, specific sediment yield, and total erosion and sediment yield can be computed using following equation. The volume of soil erosion is calculated by the following equation in this method.
Where:
H: mean annual rainfall (mm), π: 3.14, WSP: the volume of soil erosion (m3/km2/yr), T: coefficient of temperature which is calculated by equation No. 3.
Where: t = mean annual temperature
The sediment production rate in this model is calculated based on the ratio of eroded materials in each section of the stream to the total erosion in the whole watershed area (Equation No. 4).
Where:
P = circumference of the watershed, L = watershed length (km), D = height difference in watershed area (km), After calculation of the Ru value the special sediment rate is estimated by equations No. 5 and 6 (9).
Where:
GSP: special sediment rate, WSP: volume of special erosion, Ru: coefficient of sedimentation, GS: total sediment rate (m3/yr), F: total watershed area (km2).8,12
Table 5
Classification of Z coefficient values.
Study of the Validity of the Model
The model was checked against the ground data to study its validity. For this purpose, during the visits from the region at the data, the validity of the model was determined and with regard to the data from 5 control points in each subunits. T-Test was employed to compare the estimated erosion and sediment values by EPM model with measured values using SPSS statistical package.
Results
When the proper maps are prepared in EPM model, different contributing factors in erosion and sedimentation will receive their values in GIS system.
Epm model
The calculated values of four contributing factors for Ghareh Aghach watershed in EPM model is presented in Table 6. The maps of susceptibility values for each contributing factor created by ILWIS system are presented in figures 3 to 10. After the information layers were valuated to produce the erosion map, the map of erosion intensity (Z) was prepared and the susceptible sites to erosion were identified (Fig. 8). The volume of erosion (WSP) was then calculated by equation No. 2 (Table 7) Figure No. 9 shows the erosion in the watershed area. The sedimentation yield in the Ghareh Aghach watershed was classified in four erosion categories eroded (low, medium, high and very high) using EPM model.
Discussion
The Efficiency of Epm Model to Estimate Erosion and Sedimentation in Ghareh Aghach watershed
In order to determine the accuracy of EPM model, five sample points were randomly assigned to each subunits. Amounts of sediment yield can be estimated from available sediment concentration data collected at sedimentary gauge in the sub-watershed outlet. The t-student test was employed to compare the estimated erosion and sediment values by EPM model with measured values using SPSS statistical package.7 The results in Table 8 showed that there were no significant differences (P < 0.05) between the estimated and measured values. These results support the pervious achievements by the other researchers.10,13,14
Table 6
Coefficients of four contributing factors in erosion by EPM model.
Table 7
The calculated erosion and sediment for Gharehach watershed by the EPM model.
Soil erosion is a major environmental threat to the sustainability and productivity of agriculture in Iran. The average soil loss in Iran is estimated to be 20 to 30 ton/ha per year, which totals to about 5 billion tons per year, and is mainly due to inadequate agricultural land use. The agent of erosion is water and wind, but the main causes are improper land use and inadequate cultural practices.
Quantification of the actual rate and pattern of soil erosion and sedimentation is necessary for designing degradation control strategies.
Table 8
The statistical analysis and mean comparisons in EPM model.
The EPM model estimates potential amount of erosion in a specific region and it is possible that the results of this model concerns with several factors in each homogeneous unit that could determine the potential erosion more reliably. But the model has some limitation such as getting the effect of slope in resistant rocks and formations that have very low erosion. Furthermore, the scores given for different factors were developed for Yugoslavia condition. These need to be calibrated for other countries and circumstances.
The GIS system minimizes the personal effects in integrating the different information layers to identify the working units and detects and classifies the different erosion sites in the watershed.
Disclosure
This manuscript has been read and approved by the author. This paper is unique and is not under consideration by any other publication and has not been published elsewhere. The author and peer reviewers of this paper report no conflicts of interest. The author confirms they have permission to reproduce any copyrighted material.
