GIS: Some Models for Preparing Erosion Features Map
Total Page:16
File Type:pdf, Size:1020Kb
GIS: Some Models for Preparing Erosion Features Map 1 2 Mohammadi Torkashvand, A. – Haghighat, N. 1Assistant Professor, Department of Soil Science, Islamic Azad University-Rasht Branch, Rasht, Iran; 2Soil Laboratory Expert, Islamic Azad University-Rasht Branch, Rasht, Iran; Tel.: 0098-131-4224307; Fax: 0098-131-4223621; E-mail: [email protected] 1. Abstract Soil erosion is a serious geo-environmental issue causing land degradation in sub-humid to arid mediteranean countries including parts of Iran. This study in Jajrood sub-basin, north-east Tehran, was conducted to investigate some methods for water-soil erosion types mapping by GIS. Four models were used for providing working units' maps by integration of data layers including A. plant cover, geology and slope B. land use, geology and slope C. land use, rocks sensitivity to erosion and slope and D. land use, rocks sensitivity to erosion and land units. The surface, rill, gully erosion Intensity in the 314 spot was controlled to provide ground truth map from each of these erosion features. Soil erosion type's map obtained from integration of these truth maps, and then this map was crossed by the maps A to D. Results indicated that the integration of land use, rocks erodibility and land units layers is a better model for providing erosion types map than other models from an economic and executive regards. The cross of the map D with the ground truth maps of surface, rill and gully erosions showed the greatest and least accuracy are related to providing gully erosion and erosion features maps, respectively. The greatest precision of model was related to providing gully erosion map (with coefficient of variation 17.8%), and the least precision (with 40.5% coefficient of variations) was related to providing erosion features map. In conclusion, model D is suitable for preparation of gully erosion map no erosion features map. 2. Introduction Soil erosion as a serious geo-environmental has both direct and indirect negative impacts: loss of soil, loss of green cover, deterioration of agriculture, desertification and, of course, economic reverberations (Khwlie et al., 2002). The implementation of effective soil conservation measures has to be preceded by a spatially distributed erosion risk assessment (Moussa et al., 2002; Souchere et al., 2005), and erosion features assessment and its intensities (Mohammadi Torkashvand, 2005). Most of erosion and sediment studies have been carried out to provide a quantitative erosion map (Singh et al., 1992; Martinez-Casanovas, 2003; and Ygarden, 2003) and less to preparing an erosion features map. Qualitative erosion mapping approaches are adapted to regional characteristics and data availability. Resulting maps usually depict classes ranging from very low to very high erosion or erosion risk. There is no standard method for qualitative data integration, and consequently many different methods exist (Vrieling, 2006). Khawlie et al. (2002) for providing a risk map of soil erosion in eastern Mediterranean rugged mountainous areas, Lebanon, applied remote sensing and GIS. With steep slopes, torrential rain, improper human interference, run-off is high and water-soil erosion is continuously deteriorating the land cover. Remote sensing can facilitate studying the factors enhancing the process, such as soil type, slope gradient, drainage, geology and land cover. Digital elevation models created from SAR imagery contribute significantly to assessing vulnerability of hydric-soil erosion over such a difficult terrain. GIS layers of the above factors are integrated with erosional criteria to produce a risk map of soil erosion. Results indicated that 36% of the Lebanese terrain is under threat of high-level erosion, and 52% of that is concentrated in the rugged mountainous regions. Bou Kheir et al. (2006) integrated two data layers including erodibility of rocks and soil and potential sensitivity to erosion as a model for providing risk map of soil erosion. The risk map corresponds well to field observations on the occurrence of rills and gullies. In recent researches, integration of data layers has been used in erosion and sediment different studies (Feoli et al. 2002; Navas & Machin, 1997; Bayramin et al. 2003). Shrimali et al. (2001) prioritized erosion-prone areas in hills using remote sensing and GIS as a case study of the sukhna lake catchment, northern India. The study indicated that 1. IRS ID LISSIII data can be used for land use/land cover mapping with a reasonably good (83 percent) classification accuracy for hydrological and erosion assessment applications, and 2. That a simple index-based approach using three main causative factors, i.e., slope, soil and land use/land cover, can give fairly good delineation of erosion-prone areas for prioritizing. In Isfahan province, Rahnama (2003) also concluded similar results in providing erosion types map by aerial photograph interpretation. He recommended satellite imagery and GIS for this aim. The extension of the use of modern spatial information technologies, such as geographical information systems (GIS), digital elevation modeling (DEM) and remote sensing, have created new possibilities for research as a key for erosion types mapping (Martinez Casasnovas, 2003) that is economical due to low costs as well as quickness (Raoofi et al. 2004). But it should be regarded that the field survey and photo-interpretation for erosion types mapping at the national scale is time consuming and expensive (Raoofi et al. 2004). The aim of this paper is to develop a methodology based on data layers integration in GIS for providing water-soil erosion types map and its intensities at the national scale (1:250000). 3. Methods The Jajrood sub-basin with 162,558 ha located between 51°34′ E and 52°6′E, 35°13′ N and 35 °48′ N was considered for the investigation of erosion features. It extends from northeast to southeast Tehran Province, Iran. The highest and the lowest height of basin are 3000 and 867 m, respectively. Land covers were rangeland, badland, sand borrow, agriculture land and urban regions. Within the basin, different lithic units include pyroclastic stones, tuffs, andesite, shale, conglomerate, gypsum and limestone. Also, Quaternary deposits have covered in the major part of the southern basin particularly in the Varamin plain (47.8% of area basin). Necessary maps such as topographic, geology, plant cover type and land units were scanned and georeferenced. Digital Elevation Model (DEM) was prepared by 1:50,000 topographic digital data, classified slope map-the DEM- derived slope map was classified into eight slope (percent) classes 0-2, 2-5, 5-8, 8-12, 12-20, 20-40, 40-70 and >70 based on Mahler (1979) classification, land use was derived using ETM+ satellite image and rocks erodibility layer based on Feiznia (1995). According to their sensitivity to erosion, the rocks were categorized in to the five classes. Four models applied for preparation of working unit's maps that derived from integration of data layers including: A. plants cover, lithology and slope layers, B. land use, lithology and slope layers, C. land use, rocks erodibility and slope layers, and D. land use, rocks erodibility and land units layers. Selection of data layers was regarding exploratory studies in Kan sub-basin (Mohammadi Torkashvand, et al. 2005). After this, these models are called models "a", "b", "c" and "d". Different methods were incorporated for classification of soil erosion type's severities (Flugel et al., 1999; Refahi, 2000; Boardman et al., 2003 and Sirvio et al., 2004) and experiences and expertise considerations (Mohammadi Torkashvand et al. 2005). A total of 314 points has been considered on color composite images (for field investigation) by classified randomized sampling. A primary polygon was determined for each control point regarding image characteristics. The magnitude of erosion in each erosion feature was investigated in these ground control points and then frontiers of each primary polygon were corrected with due attention to the field views for every one of the surface, rill and gully erosions. Modified polygons with regard to the intensity of each erosion features in the field, were marked. Polygons with same the intensity were combined together and ground truth maps of surface, rill and gully erosions were prepared. The erosion features map obtained from integration of the surface, rill and gully erosions maps. Erosion features maps were crossed with working units' maps to investigate the ability of each method on separating erosion features. Equation 1 was used for investigating method's accuracy. n Z * c (1) ∑ ( x i ) i i = 1 A = n Z * ∑ ( x i ) i = 1 That A is map accuracy or map conformity with actual conditions (percent), Z * is working units' area (ha) ( xi ) and ci is maximum area of each working unit that is uniform in compared to actual conditions (percent). Root mean squared error of working units’ accuracy was computed by equation (2). n 2 []Z − Z * ∑ ( xi ) ( xi ) (2) RMSE = i =1 n That Z is working units' area (ha) that is uniform in actual conditions. The precision each of method by (xi ) S equation 3 was obtained. CV = .100 . That CV and S are coefficient of variation and standard deviation of X working units' accuracy, respectively, and X is method accuracy ( A in equation 1). 4. Results Table 1 indicates the numbers of units' in the working units' maps. In layers integration models, the most and least numbers of working units are related to models "a" and "d", respectively. Most polygons of models "a", "b" and "c" have small area which is not possible to be presented in the maps 1:250,000 maps due to cartographic limitations. Table 1 presents the accuracy and error of data layers integration models in distinguishing soil erosion type's intensities. According to this table, the highest and the least accuracy belong to models "a" and "c" that is 68.3 and 53.4 percent, respectively.