Assessing Environmental Sensitivity Areas to Desertification in North of Iran

The sensitivity map to desertification of Chehel-chai basin (North of Iran) has been elaborated by the crossing of fives criteria including climate, water erosion, soil, vegetation, and management, which have the main impact on the evolutionary process of desertification. The methodology is resulted from the desertification model of MEDALUS, which use qualitative index to define the land sensitive area to desertification. According to the factorial scaling technique, score-ranging from 1(good condition) to 2(deteriorated condition) is assigned to each indicator. Each index was assessed based on selected indicators, which resulted in qualitative mapping of each index based on geometric average of the indicators. Thematic databases, with a 1:50000 scale resolution, were integrated and elaborated in GIS software based on arc view3.2, Ilwis3.1 and ArcGIS9.3. The obtained results indicate that 39.39% of Chehel-chai basin is highly sensitive to desertification, 2.13 % has moderate sensitivity, only 2.43% has low sensitivity and 56.05% is non-sensitive. Studying the mean weight of numerical value, it is distinguished that the intensity of desertification for the total area is 1.32 which is classified as fragile (2).

Defined as the change of fertile and productive lands into non-productive desert areas, desertification remains as a major environmental concern in most arid, semiarid, and sub-humid areas of the world (Dregne, 1977;UNCCD, 1994;Herrmann & Hutchinson, 2005).Desertification may results from various natural and humanassociated factors such as wind or/and water erosion, destruction of vegetation covers and water resources, water logging, soil salinization and alkalization, etc. (Kassas, 1977;Mainguet, 1994;Blum, 1998;Krishan et al., 2009).Among these factors, unsustainable human activities, such as urbanization, overgrazing, over-cultivation, poorly drained irrigation systems, deforestation, etc. are currently considered as the most important factors accelerating the process of desertification (Barbero-Sierra et al., 2013; Bo et al., 2013).In addition to these direct effects, human activities may also indirectly contribute to desertification, for example through affecting environmental factors such as climate (Kassas, 1977;Khosravi, 2004).
The United Nations Convention to Combat Desertification (UNCCD), which is the only internationally legally organization, offers strategies to combat desertification and mitigate the effects of drought through national action programs that incorporate long-term strategies supported by international cooperation and partnership arrangements.Several models have been developed and a large number of studies have been carried out to assess the risk of desertification in different areas most of which are specialized to particular geographic areas (Santinia et al., 2010).To develop these models to local areas of interest, it is necessary to re-investigate and adjust them to local environments (Geeson et al., 2002).
During the last decades, many areas of Iran have been subjected to intense desertification because of drought, unsustainable land use, and increasing pressures on land and water resources (Amiraslani & Dragovich, 2011).In this study, we used the Mediterranean Desertification and Land Use (MEDALUS) project to assess and map desertification sensitivity in a watershed located at north of Iran.In this project, the focus is primarily on Mediterranean environments where physical loss of soil by water erosion, and the associated loss of soil nutrient status are identified as the dominant problem (Brandt & Thornes, 1996;Geeson et al., 2002;Arar et al., 2009).The MEDALUS project is a popular GISbased technique which computes an index, such as the Environmentally Sensitive Areas (ESAs) index containing several environmental (e.g.climate, soil vegetation) and anthropogenic (e.g.management) parameters (Kosmas et al., 1999).The indicators selected to evaluate the sensitivity to desertification can provide an overview for evolution of ecosystems and environments, which can be applied as efficient tools for decision-making and planning (Basso et al. 2000).Results of this study would provide a map of desertification sensitivity exhibiting the area of degraded land, process type, dominant indicator involved, and intensity classes of desertification.

Study area
The Chehel-chai River is one of the largest branches of the Gorganrood River located at Golestan province, northeastern part of Iran.This river drains an area of about 250 km 2 across geographic coordinates ranging from 55°23´ E to 55°38´ E and from 36°59´ N to 37°13´ N (Figure 1).The elevation of this watershed ranges from 190 to 2570 m (a.s.l) and the mean slope is 45.82%.Limestone formations, alluvial deposits near the streams and quaternary sedimentary formations, loess, dominate the underlying geology in the lowlands of the study area.The region has a Mediterranean climate with warm and dry summers and temperate and rainy winters.Falling occurs mainly in winter and spring with the mean annual precipitation is about 750 mm.Therefore, the river Fig. 1: the location of chehel-chai watershed Fig. 2: the view of the chehel-chai watershed has a relatively high discharge during winter and spring, but low during summer and autumn.Forests (60%) and croplands (40%) are predominant land uses in the watershed (Figure 2).

mEthoDology
We used the MEDALUS methodology, with slight modifications, to map vulnerable areas to desertification using means of environment-state and response indices (Geeson et al., 2002, Kosmas et al., 1999).This methodology has been presented as a manual of key indicators of desertification and mapping environmentally sensitive area to desertification (Kosmas et al., 1999).
Based on their applicability and susceptibility to environmental processes, the risk of desertification in Chehel-Chai watershed was evaluated on a regional scale by defining levels on the basis of some parameters or indices of five categories including climate, vegetation, soil, water erosion, and management trough field investigations and spatial data.Each parameter was weighted in relation to its impact and contribution to desertification process.
All data required for this method were elicited from previous studies or calculated or collected from field samplings (Table 1).These data were then introduced to GIS system to calculate required indices and visualize a desertification map for the study region (Rafiei Emam, 2003).
To establish the relative scale of severity, those threshold values with critical conditions were selected to normalize outputs, and the results were then weighted using lands in the study area (Al-Adamat et al., 2003).Finally, a single index obtained via integrating all of these indices was selected as a representative of the desertification risk in the study area.
Five general indices (Quality indices) were considered by the model each containing various parameters (indicators) (Table 1).To calculate the individual indices, detailed studies were conducted on hydrologic, geomorphologic, erosion, and soil and vegetation cover characteristics of the region.The detailed information of the individual indicators and their scores has been summarized in Tables 2,  3, 4, 5, and 6.
For each index, a score ranging between 1, as good condition, and 2, as deteriorated condition, was assigned based on factorial scaling technique.Additionally, A zero value was assigned when the measure was not appropriate for an area or where that area was not classified.In most cases, a linear function ranging from 1 to 2 (extreme values) represents the variation of the indicators (score).Although in some particular cases, a non-linear variation is possible.The value of quality index for each elementary unit within a layer was obtained as geometric average of scores of single indices according to the following formula: Where Index-X is a given index, Layer is the indicators of each index and N is the number of indicators for each index.
The five quality indices were computed using ArcGIS software v. 9.3 as the geometric mean of the indicators related to each single index.The importance of each index in desertification was determined using the following formula:

•
Water Erosion Index (WaEI) = (Erosion critical limit × Water erosion risk)  The maps of abovementioned indicators were obtained based on modified model and a final desertification map was provided by multiplying the prepared maps.Here, the value of each index was Where DM is desertification mapping and the other abbreviations are the same as the above five formula.
The intensity ranges of desertification sensitivity for each class have been shown in Table 8.
In order to better clarify the boundaries between each fragile and critical classes, they were divided by three sub-groups ranging from 1 as low sensitivity to 2 as high sensitivity.

rESultS
The maps of sensitivity to desertification of different indices, provided by the ESAs model, have been shown in Figure 3.
The geometric average and qualitative classes of indicators of each index have been also summarized in Tables 9 and10.Analysis of desertification indices in Chehel-chai basin showed that soil quality index was the major problem of desertification in the study area with a geometric average of 1.57, which shows very high desertification class.The vegetation quantity index with a weighted average of 1.53 was the second major factor underlying desertification in the studied area.
Based on five studied indices, the assessment of frequency distribution classes in terms of desertification risk for current statues showed that mean weight of quantitative values was 1.32 for all study area compared to ESAs table (Table 9) implying that this area has a desertification class of fragile (class 2) (Figure 4).
In addition, the result showed that 56.05% of the area fits into "absent class" of no sign of desertification.2.43% of the region was allocated in low class, which means that land degradation has some little effect on it.This area should be considered as a sensitive area.2.43% of the region was allocated in the medium class where desertification has obvious effects.Finally, 39.39% of the region was allocated to high classes of desertification.These areas, mainly located in the south and southwestern part of the region, are critically sensitive to desertification (Figure 5).

DIScuSSIoN
In this study, we used the MEDALUS method to assess risk of desertification in Chehel-Chai watershed.The ESAs model is very advantageous and user-friendly model for determination of different classes of vulnerable areas to desertification (Arar  This Chehel-Chai watershed is believed to be exposed to high rates of desertification, such that more than 57% of the total surface had an ESAs value of more than 1.5.This has been attributed to a combination of factors including land weakening and degradation as a result of unprincipled human activities, low qualities of vegetation surfaces, and low quality soils.The low-quality vegetation surfaces occur predominately in southern part of the region, Results of this study showed that the use of GIS is helpful in assessing areas at risk of desertification because of saving time, providing precision, and reliability in future studies.This system provided a general overview of desertification status as a map of desertification vulnerability (Chenchouni et al., 2010).This approach will be worth for decision makers for developing the best guidelines to combat against desertification in sensitive lands. of course, some other detailed studies are still required to recognize the most sensitive areas to desertification.Human-related factors seem to be currently the most important factors influencing desertification in Chehl-Chai basin, since environmental limitation such as rainfall inadequacy, prolonged dry periods, extreme temperature and evaporation, salinity and soil alkalinity do not exist in this area.The role of climate changes in desertification is more evident in arid and semi-arid regions.Deforestation and conversion of forests to rain-fed lands is one of the main causes for starting the process of land degradation and desertification in the study area. of course, other natural factors such as climatic factors like storms, early cold, drought, geological factors like landslides (mass erosion) also plays a role in the destruction of the forest of basin.
Improper management of land use in a watershed has undesirable effects on the available resources.optimization of land use is one of the most useful strategies for achieving sustainable development and reducing wasted resources (Jalili et al., 2007).Use of lands based on their potential in proper managerial planning can reduce destruction and losses of land (Chapi, 1997).

Fig. 3 : 3 •
Fig. 3: maps of water erosion index (a), vegetation quality index (b), soil quality index (c), climate quality index (d) and management quality index (e)

table 3 :
classes and corresponding weights assigned for calculation of the vegetation quality index (vQI)

Fig. 4 :
Fig. 4: characterization of the area sensitive desertification in chehel-chai with a whole set of supplementary indices

table 1 : Desertification indicators and quality indices used in ESAs model for visualization of the chehel-chai desertification status map
* Nitrogen (N), Phosphate (P), and Potassium (K) content ** Defined as the means of Transeau ratio (P-ETP) where P is the average annual rainfall (mm) and ETP is the average annual potential evapotranspiration (mm)