Erosion is a type of soil and rock weathering in which superficial layer is worn away through the action of glaciers, wind and water  . Although both soil and rock erosion are important, cultivation and life system soil erosion is more important, and makes social and economic problem and ecosystem changes    .
Erosion monitoring is more effective than control of sediment using different methods. By detecting different erosion problems and classifying the area into different levels of susceptibility to erosion, it will be possible to use cost effective erosion control measurements  .
Different erosion types like rill, gully, stream and channel are used for erosion severity and sediment production intensity from which gully erosion is more important, as it is occurred in the last stages of erosion progress, indicating high sediment production, but its study has been limited recently  . Whilst, study on gully erosion based on historical evidence, photology or remotely sensed images, is very practical. The result is also very applicable, as there is a strong correlation between the length and volume of the gully  .
Drainage area of any basin is assumed as its sediment source, so the area is mapped and classified into erosion susceptive aprons to make remedial programs in order to reduce the erosion intensity, decrease soil erosion and overcome to sediment production and transform to dam reservoirs. This study uses different factors such as lithology, drainage pattern, steep slopes, erosion types and land use, to map and classify the drainage area using different models, and make comparison for the results.
2. Study Area
The 78.2 Km2 study area is located in the northern Hashtrood City, between 37˚28'51''N - 37˚36'39'N and 46˚59'59''E - 47˚05'36''E (Figure 1) with 2110 and 1668 m above sea level maximum and minimum elevations respectively. Average precipitation is 325 mm, and could be assumed as cold semiarid. The main land use, 82 percent of the study area, is dry farm.
Geological studies illustrate the area has rarely been affected by tectonic forces, and sedimentary rocks are the predominant lithology. Unconsolidated conglomerates, sandy marls and alluvial sediments cover 81% of study area, which cause high potential for erosion.
Two empirical methods, Fargas and BLM were used in this study to identify the critical sediment resources. The first one as a qualitative method which was developed by Fargas et al.  was initially used in 1500 Km2 area watershed of the Joaquin Costa reservoir (NE Spain). This method uses reliable but limited parameters of lithology and drainage density which is rated according to erosion class. It is supposed that too less factors are easy to obtain, but could be the weak point.
The second one, BLM methodology is also applicable method for sediment risks evaluation. This technique uses seven factors, and the results are classified according to five major erosion classes. If these factors study carefully, the result will be reliable.
Figure 1. Location of study area.
Two main models, Fargas and BLM, which are generally used for erosion assessment and soil conservation planning, were used to categorize the erosion intensity, erosion risk assessment, sediment production and critical sediment resources detection. Finally the results of above mentioned methods were appraised in comparison with each other and the similarities were recognized.
4.1. Fargas Model
The Fargas model is based on two main factors, lithology and drainage density. So, these factors are prepared as geological maps, and used in GIS software (Figure 2) and each unit area was calculated as percentage of total area (Table 1). This study uses geological map which was prepared based on detail field study as well as laboratory analysis, considering its importance. According to the map and Table 1, Qt2 geological unit, composing of terrace alluviums which have high resistance index (erosion susceptible) covers about 46.5% of the area, whilst the hard rocks of gr, Or, Eva and Ea have frequency of 14.26, 0.68, 0.68 and 3.67 percent, respectively. Qplc, Qplm and Qal are other soft lithological units with 4.7, 27.56 and 2.21 percentage, correspondingly.
The drainage density, the second main factor of this model, was prepared based on the detailed topography maps (Figure 3); then, erosion class which is attached to an
Figure 2. Geological map of the study area.
Figure 3. Drainage density on geological map.
Table 1. Geological units resistance index  .
elementary catchment depending on its drainage density  , were prepared (Table 2).
The coefficient value for each lithological unit from which the emission risk is obtained, is determined by multiplying the erosion class according the drainage densityand the weighting factor  . Table 3 shows the calculated values and the erosion risk for different geological units of study area. The critical sediment source could be calculated based on the risk data. This table shows the erosion risk of slight, moderate, high, severe and very severe as 0.68%, 0.33%, 3.67%, 14.26% and 81.06% of total area, respectively. The data is presented as map of erosion risk and critical sediment source (Figure 4).
4.2. BLM Model
BLM methodology, initiated by Bureau of Land Management (BLM) of the United State of America, is mainly base on the erosion types. This model is integrated seven factors of surface erosion variable from (0 - 14), Litter covering the surface (0 - 14), rock fragment (0 - 14), superficial degradation (0 - 14), rill erosion (0 - 14), runoff sediments precipitation (0 - 15) and gully erosion (0 - 15); to calculate the Soil Surface Factor (SSF), which is the base for assessment of erosion intensity and erosion factor determination (Table 4). Every factor of this model was studied and mapped in GIS software to get the results as Table 5, Figure 5 and Figure 6. According to Table 5, the study area could be dissection in 3 major erosion risk, including moderate, high and severe with area of 32.1%, 42.96% and 24.94% respectively.
Fargas and BIM models are mostly used to determine erosion risk, erosion intensity and critical sediment source. Fargas model generally is more common, as it uses less factors, which are easy to produce the data. BLM model covers more components and needs data to be studied carefully and integrated to one aspect as soil surface factor. This study uses both model to evaluate erosion severity in the same area, in order to find the similarities and differences. The results show, generally there are similar classes between both model results, but deep investigation on the calculated classes area, re-
Table 2. Erosion class based on drainage density.
Table 3. Classification of sediment emission risk.
Table 4. The erosion risk based on the 7 factors (bureau of land management).
Table 5. Erosion risk in the study area based on SSF calculation.
Figure 4. Risk of erosion based on Fargas method.
Figure 5. Erosion types in the study area.
Figure 6. Erosion risk based on BLM method.
veals that there are no enough coverage between the final outcomes of these models. Fargas model predicts 3.67%, 14.26% and 81.06% of the area susceptible for high, severe and very severe erosion respectively; whilst, referring to BLM model outcome, 42.96% area has high sensitivity; 42.96% and 24.94% of the area have high and severe sensitivity for erosion, respectively. Furthermore, both models show the same severity for around 18% of the study area. This means, by applying two models in the same study and making comparison between the products, we find advantages and disadvantages of the methods, and understand the shortage of any model, or any faults during the study. Finally it is recommended to never rely on one of these models, or maybe any others, even very common methods.
This study is based on the research plan of “Evaluation of Sediment Production intensity in Different phases of basic studies (Case Study, North of Hashtrood City)” which is supported by Islamic Azad University (IAU), Ahar Branch. We are thankful to the chancellor, research deputy and all of the staffs of the university which gave support us in this issue.
 Moran Lu, H. and Prosser, C.J. (2005) Modeling Sediment Delivery Ratio over the Murray Darling Basin. Environmental Modelling & Software, 21, 1297-1308.
 Santhi, C., Srinivasan, R., Arnold, J.G. and Williams, J.R. (2006) A Modeling Approach to Evaluate the Impacts of Water Quality Management Plans Implemented in a Watershed in Texas. Environmental Modelling & Software, 21, 1141-1157.
 Miller, S.N., Semmens, D.J., Goodrich, D.C., Hernandez, M., Miller, R.C., Kepner, W.G. and Guertin, D.P. (2007) The Automated Geospatial Watershed Assessment Tool. Environmental Modelling & Software, 22, 365-377.
 Poesen, J., Nachtergaele, J., Verstraeten, G. and Valentin, C. (2003) Gully Erosion and Environmental Change: Importance and Research Needs. Catena, 50, 91-133.
 Cheng, H., Zou, X., Wu, C., Zhang, C., Zheng, Q. and Jiang, Z. (2007) Morphology Parameters of Ephemeral Gully in Characteristics Hillslopes on the Loess Plateau of China. Soil and Tillage Research, 94, 4-14.
 Fargas, D., Martínez-Casanovas, J.A. and Poch, R. (1997) Identification of Critical Sediment Source Areas at Regional Level. Journal of Physics & Chemistry of the Earth, 22, 355-359.