Translator Disclaimer
1 January 2020 Estimate of Erosion and Sedimentation in Semi-arid Basin using Empirical Models of Erosion Potential within a Geographic Information System
Author Affiliations +
Abstract

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

Figure 1

Scheme for developing homogeneous unit map.

10.4137_ASWR.S3427-fig1.tif

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.

Figure 2

EPM model for soil erosion.

10.4137_ASWR.S3427-fig2.tif

Table 1

The contributing factors in EMP model to estimate the soil erosion.

10.4137_ASWR.S3427-table1.tif

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.

10.4137_ASWR.S3427-table2.tif

Table 3

The soil erosion sensitivity coefficient.

10.4137_ASWR.S3427-table3.tif

Land Use coefficient

The amount of land use coefficient was determined according to Table 4.

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:

10.4137_ASWR.S3427-eq1.tif
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.

10.4137_ASWR.S3427-table4.tif

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.

10.4137_ASWR.S3427-eq2.tif

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.

10.4137_ASWR.S3427-eq3.tif
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).

10.4137_ASWR.S3427-eq4.tif

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

10.4137_ASWR.S3427-eq5.tif
10.4137_ASWR.S3427-eq6.tif
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.

10.4137_ASWR.S3427-table5.tif

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.

10.4137_ASWR.S3427-table6.tif

Figure 3

The Slope map of Ghareh Aghach watershed.

10.4137_ASWR.S3427-fig3.tif

Figure 4

Uniform Sampling Units map of Ghareh Aghach watershed.

10.4137_ASWR.S3427-fig4.tif

Figure 5

The Geomorphology Facieses map of Ghareh Aghach watershed.

10.4137_ASWR.S3427-fig5.tif

Figure 6

The Geology erosion map of Ghareh Aghach watershed Based on erosion sensitivity.

10.4137_ASWR.S3427-fig6.tif

Figure 7

The Land use map of Ghareh Aghach watershed Based on erosion sensitivity.

10.4137_ASWR.S3427-fig7.tif

Figure 8

The Erosion classes map (Z) of Ghareh Aghach watershed.

10.4137_ASWR.S3427-fig8.tif

Figure 9

The Intensity Erosion map (WSP) of Ghareh Aghach watershed.

10.4137_ASWR.S3427-fig9.tif

Figure 10

The Vegetation map of Ghareh Aghach watershed.

10.4137_ASWR.S3427-fig10.tif

Table 7

The calculated erosion and sediment for Gharehach watershed by the EPM model.

10.4137_ASWR.S3427-table7.tif

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.

10.4137_ASWR.S3427-table8.tif

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.

REFERENCES

1.

Meijering, A.M.J.Erosion and sediment yield in catchments, ITC, Iecture notes. The Netherlands. 1995. Google Scholar

2.

Lu, H.J., Geollant prosser, L., Priestly, Moran GPrediction of sheet and rill erosion over the Australian continent: Incorporating monthly soil loss distribution. National Land and Water Resource. Technical Report; 2001. No. 31, Australia. Google Scholar

3.

UNDP. Human Development Report of the Islamic Republic of IRAN. 1999; Chapter 8, pp. 109–21. Google Scholar

4.

Nearing, M.A., Govers, G., Norton, L.D.Variability in soil erosion data from replicated plots, soil. Sci Soc Am J. 1999: 1829–35. Google Scholar

5.

Lufafa, A., Tenywa, M.M., Isabirye, M., Majaliwa, M.J.G., woomer, P.L.Prediction of soil erosion in alake victoria basin catchment using a GIS- based universal soil loss model. Agri Systems. 2003; 76: 883–94. Google Scholar

6.

Daroussin, J., King, D.Mapping erosion risk for cultivated soil in France. Catena. 2001: 46; 207–20. Google Scholar

7.

Isfahan Regional Water Organization. Study of Physiology, Climatic, Hydrology, and Geomorphology in Ghareh Aghach watershed. Isfahan University of Technology. 2006. Google Scholar

8.

Ahmadi, H. Application Geomorphology, Vol.1 (Water Erosion). Publish by Tehran University. 2006. Google Scholar

9.

Rafahi, H., Namati, M.Using EPM for study erodibility and sediment yield in Alamout watershed. Iranian Journal of Agriculture Science. 1994; 26(2): 32–45. Google Scholar

10.

Gaverlovic, Z.The use of an empirical method (Erosion Potential Method) for calculating sediment production and transportation on in unstudied or torrential streams. Proc Int Conf on River Regime. 1988: 411–22. Google Scholar

11.

ITC. ILWIS 5.4 Users Manual. Computer Department, International Institute for Aerospace Survey and Earth Science, Enscheda.The Netherlands.2007. Google Scholar

12.

Rafahi, H Gh Water Erosion and Conservation, Vol. 1. Publish by Tehran University. 2004. Google Scholar

13.

Shakeri, Sh, Blalipour, F.Comparing MPSIAC and EPM model for estimating erosion and sediment yield in Dou Abe watershed. Conference of water resource management. Isfahan University of Technology. 1996: 272–81. Google Scholar

14.

Shahbazi kia, S., Restami, S.Comparing MPSIAC and EPM model for estimate erosion and sedimentation in in Azarbayejan Province. International Conference of Erosion and Sedimentation. 2006: 598–602. Google Scholar

15.

Amiri, F.Modeling multiple use of rangeland by using GIS. PhD thesis, Islamic Azad University Research and Science Branch. Tehran. 2008: pp. 560. Google Scholar
© 2010 SAGE Publications. This article is distributed under the terms of the Creative Commons Attribution-NonCommercial 3.0 License (http://www.creativecommons.org/licenses/by-nc/3.0/) which permits non-commercial use, reproduction and distribution of the work without further permission provided the original work is attributed as specified on the SAGE and Open Access page (https://us.sagepub.com/en-us/nam/open-access-at-sage).
Fazel Amiri "Estimate of Erosion and Sedimentation in Semi-arid Basin using Empirical Models of Erosion Potential within a Geographic Information System," Air, Soil and Water Research 3(1), (1 January 2020). https://doi.org/10.1177/ASWR.S3427
Published: 1 January 2020
JOURNAL ARTICLE
PAGES


SHARE
ARTICLE IMPACT
Back to Top