Land Suitability Assessment on a Watershed of Loess Plateau Using the Analytic Hierarchy Process

Abstract

In order to reduce soil erosion and desertification, the Sloping Land Conversion Program has been conducted in China for more than 15 years, and large areas of farmland have been converted to forest and grassland. However, this large-scale vegetation-restoration project has faced some key problems (e.g. soil drying) that have limited the successful development of the current ecological-recovery policy. Therefore, it is necessary to know about the land use, vegetation, and soil, and their inter-relationships in order to identify the suitability of vegetation restoration. This study was conducted at the watershed level in the ecologically vulnerable region of the Loess Plateau, to evaluate the land suitability using the analytic hierarchy process (AHP). The results showed that (1) the area unsuitable for crops accounted for 73.3% of the watershed, and the main factors restricting cropland development were soil physical properties and soil nutrients; (2) the area suitable for grassland was about 86.7% of the watershed, with the remaining 13.3% being unsuitable; (3) an area of 3.95 km2, accounting for 66.7% of the watershed, was unsuitable for forest. Overall, the grassland was found to be the most suitable land-use to support the aims of the Sloping Land Conversion Program in the Liudaogou watershed. Under the constraints of soil water shortage and nutrient deficits, crops and forests were considered to be inappropriate land uses in the study area, especially on sloping land. When selecting species for re-vegetation, non-native grass species with high water requirements should be avoided so as to guarantee the sustainable development of grassland and effective ecological functioning. Our study provides local land managers and farmers with valuable information about the inappropriateness of growing trees in the study area along with some information on species selection for planting in the semi-arid area of the Loess Plateau.

About the Authors:

Xiaobo Yi

Affiliations College of Resources and Environment, Northwest A&F University, Yangling, Shaanxi, China, State Key Laboratory of Soil Erosion and Dryland Farming on the Loess Plateau, Northwest A&F University, Yangling, Shannxi, China

Li Wang

* E-mail: [email protected]

Affiliations College of Resources and Environment, Northwest A&F University, Yangling, Shaanxi, China, State Key Laboratory of Soil Erosion and Dryland Farming on the Loess Plateau, Northwest A&F University, Yangling, Shannxi, China

Introduction

Sandy desertification and soil erosion are two of the most serious problems affecting China's water and land resources. The World Bank has suggested that potentially more than 331 million hectares of land are susceptible to desertification (about one third of the area of China) while about 262 million hectares are actually affected [1]. Previous studies have shown that soil erosion affects about 360 million hectares of land in China, which is about 38% of its total area, and this proportion is more than three times the global average [2]–[3]; land in China is considered to be amongst the most severely eroded in the world. For example, on the Loess Plateau, erosion rates are about 8000–25,000 tonnes km−2 year−1 in the gully areas of Shanxi and Shaanxi Provinces. Intensive cultivation of the steep hillsides has resulted in the loss of an estimated 1.6 billion tonnes of soil annually to the Yellow River [4]. In addition, dry conditions combined with the fine texture of the loess soil make the area very susceptible to dust storms. Due to the severe soil erosion and sandy desertification, the eco-environment of the Loess Plateau has been severely degraded, which has seriously affected sustainable development in the region [5]. With the environmental goals of reducing soil erosion and desertification, the Chinese Central Government launched the Sloping Land Conversion Program (SLCP, also known as Grain for Green or Grain to Green) in the late 1990s, with the intension of increasing the country's forest and grassland cover by retiring steeply sloping and marginal land from agricultural production [6]–[9]. This program aimed to restore degraded ecosystems, reduce poverty and assist rural households to move towards more sustainable economic activities [2], [10].

However, after more than 10 years, this large-scale ecological-recovery project has faced some key problems that have limited the successful development of the current SLCP policy [11]. Surveys and case studies of SLCP have consistently identified insufficient technical support and arid conditions as being the key constraints to achieving program goals [3]. For example, in the north part of the Loess Plateau, the potential for converting cropland into forests has been over estimated: tree species have been planted in areas that are better suited to growing shrubs or grass, and species with high water requirements have been planted in areas where drought-tolerance is required. The result has been soil desiccation and the development of dry soil layers because the planted trees exploited water stored in the deeper soil layers and prevented them from being recharged by rainwater [12], [7]. Consequently, once the stored water was exhausted, the limited precipitation was insufficient to maintain normal growth in the re-vegetated areas [13], leading to the vegetation dying or the production of stunted trees (colloquially referred to as “little old man trees”, which are only about 20% of the normal height for their age) [14]–[17]. Thus, implementing the SLCP policy on the Loess Plateau requires careful consideration of several factors that have a significant effect on land use change and land suitability associated with carrying out the re-vegetation program. Currently, China is facing increased environmental pressures with shortages of water potentially limiting development, especially in its dryer northern and western regions, including the Loess Plateau [12]. In order to reduce soil water depletion and implement a successful re-vegetation program for environmental improvement on the Loess Plateau it is essential to identify appropriate species to plant. It is important to evaluate the land suitability and to develop land use plans that maximize reduction in soil erosion while at the same time minimizing water yield reduction, thus ensuring the sustainable growth and succession of vegetation. Only the proper implementation of the SLCP and successful re-vegetation can make a real contribution to efforts being made to combat the urgent environmental problems of soil erosion and desertification, as well as of climate change and loss of biodiversity, currently confronting the Loess Plateau, in particular, and China as a whole. In this study, the analytic hierarchy process (AHP) was used to assess the suitability of cropland, grassland and forestland in the Liudaogou watershed, and the main limiting factors for different land use options were quantitatively analyzed. The AHP has been widely used by decision-makers and researchers. It is a mathematical method, developed by Saaty in 1977 [18] and improved by the same author in 1980 [19], for analyzing complex decisions involving many criteria [20]. It has been widely used in site selection, suitability analysis, regional planning, and land consolidation analysis [21]–[22].The success of the AHP as a practical and reliable method is highlighted by its extensive application in the past two decades [23]–[24]. Furthermore, its simplicity in relation to its power was a significant factor in the choice for its use in the presented study, and all the mentioned factors ensured that the study objectives would be successfully achieved. The objectives of the study are to determine which land use is best suited to the re-vegetation program by assessing the land suitability, to identify the constraints for future land conservation, and to provide a scientific basis for decision-making in the successful implementation of the Sloping Land Conversion Program, not only for the study's watershed but across the whole Loess Plateau.

Materials and Methods

Study area

This study was conducted in the Liudaogou watershed, located in Shenmu County, Shaanxi Province, China (110°21′–110°23′E, 38°46′–38°51′N; Fig. 1). The Liudaogou watershed covers 6.9 km2 and is located at the center of the wind-water erosion crisscross region in the north part of the Loess Plateau. This area suffers its most serious water erosion in summer and autumn and its most serious wind erosion in winter and spring [25]–[26]. The watershed is representative of the ecotone between the grass–pastoral and the agricultural areas, in the transitional zone between the desert aeolian deflation zone (the Mu Us Desert) and the loess hilly area (Loess Plateau), as well as being between the arid and the semi-arid regions. For these reasons, the Shenmu Erosion and Environmental Research Station (SEERS) of the Chinese Academy of Sciences was built within the watershed [27]. The study area is located at altitudes ranging from 1080 to 1270 m above mean sea level and has a semiarid continental monsoon climate, with a mean annual temperature of 8.4°C. The monthly mean temperature ranges from 9.7°C in January to 23.7°C in July. The mean annual precipitation is 437 mm, 77% of which occurs from June to September. Mean annual potential evapotranspiration can be as high as 1800 mm, which would result in a water deficit of 1350 mm year−1. The area has a deep (up to 100 m) loess layer, which originated during the Quaternary period. The dominant soil type (cultivated loessial soil), is a Ust-Sandic Entisol, which is loess-derived and is consequently easily eroded by both water and wind.

[Figure omitted. See PDF.]

Figure 1. Location of the study site.

(a) China; (b) the Liudaogou watershed.

https://doi.org/10.1371/journal.pone.0069498.g001

Analytical hierarchy process (AHP)

The application of AHP for making a decision about land suitability in this study involved four main steps, as follows:

The first step was to decompose the decision problem into a hierarchical structure where the attributes and plans were present as inter-related elements. Based on a qualitative analysis of the environment in the study area, the final hierarchical structures were separated into four levels (Fig. 2). The first level was the overall land suitability. The second level was composed of subsystems: geological and topographical conditions; nutrient status; and soil physical properties. The third level consisted of the specific factors that affected the land suitability. The fourth level comprised each assessment unit (cell).

[Figure omitted. See PDF.]

Figure 2. Hierarchical structure of land suitability.

https://doi.org/10.1371/journal.pone.0069498.g002

Elevation (C1) and slope gradient (C2) were selected to represent geological and topographical conditions (B1). Slope is a crucial factor affecting vegetation structure and soil erosion on the Loess Plateau. Variation in elevation has an impact on soils, microclimatic effects, and other processes that could affect land suitability. Soil organic matter (C3), total nitrogen (C4), available phosphorous (C5), and available potassium (C6) were selected as factors that represented nutrient status (B2). These four nutrition factors are closely related to land use and cover change, and understanding the effects of land use change on soil organic matter and nitrogen is important to sustainable management of land resources and associated watershed processes, as well as regional responses to global climatic change [28]. Three factors were selected to represent the soil physical properties. Soil bulk density (C7), soil texture (C8) and soil water content (C9) strongly influence plant growth and land use, and soil properties and plant recovery processes have distinct characteristics in the typical desertified sandy land of the Loess Plateau [17]. Socio-economic factors were not considered in this study because it was limited to a small watershed where such factors are generally less important than the physical ones.

The second step involved a pair-wise comparison of the elements based on a nine point weighting scale; this generated the input data (Table 1). The comparison was carried out for each decision element at 1–(–1) levels, where was the matrix size. A matrix was generated as a result of the pair-wise comparisons and weights for the criteria were derived from these calculations.

[Figure omitted. See PDF.]

Table 1. AHP pair-wise comparison scale for variables i and j.

https://doi.org/10.1371/journal.pone.0069498.t001

A matrix of scores could be developed from the comparisons, given bywhere indicated how much more important the objective was than the objective, while making a suitable material handling/equipment selection decision. For all and , it was necessary that and . The possible assessment values of in the pair-wise comparison matrix, along with their corresponding interpretations, are shown in Table 1.

The scores were normalized by dividing each element within the matrix by the sum of the column , to create a normalized matrix, :

The weight associated with each objective () could be estimated as the mean of the normalized scores in row of the matrix. Hence, was calculated to give a matrix of weights, C:where is the number of objectives being compared, and the value indicates the relative degree of importance (weight) of the objective.

The third step was to check for the consistency of the weight values () underlying the theoretical validity of the comparison matrix. In order to determine consistency, the consistency vector ( matrix) and were calculated as:The eigenvalue of the pair-wise comparison matrix, , was then estimated using the following equation:An approximation to the consistency index () was calculated and the consistency judgment was checked for the appropriate value of by . and were calculated as [18]:where is the random consistency index. The values for different numbers of are listed in Table 2.

[Figure omitted. See PDF.]

Table 2. Random consistency index for AHP.

https://doi.org/10.1371/journal.pone.0069498.t002

The rule commonly used in the AHP was applied, whereby if was less than 0.10 (i.e. 10%), the degree of consistency was acceptable; and if it was greater than 0.10 it was considered that there were significant inconsistencies and, in such a case, the AHP would not produce meaningful results [29], and it would be necessary to review and improve the judgments.

The fourth step was a complete evaluation based on the rating of the final weights in the decision plan. Matrixes for layers B and C and for C and D were generated using the same method as that used for layers A and B. By comparing the final values by simple rankings, the weights of all elements in each level of the hierarchy relative to the entire level could be obtained. These in turn were all ranked, and were carried from the upper layer to the lower layer. The combined weight () of each assessment factor was determined for the integrated assessment of land suitability in the Liudaogou watershed.

Membership value standardization for assessment factors

In the process of land suitability assessment, a primary step was to ensure a standardized measurement system for all the factors considered. Since those factors have different standards of measurements, they had to be standardized to a uniform rating scale; in this study the scale was between 1 and 4 for ease of analysis. Assigning values to specific factors required specific decision rules in the form of thresholds for each factor. Various statistical and empirical guidelines from the related national codes and literature were used to determine the boundary values. As a general guideline, a positive correlation between the value awarded and suitability was employed. The class boundaries and standardized measurements employed for each factor are given for arable land in Section 3.1.1 and for forestland and grassland in Section 3.2.1. The integer numbers ranging from 1 to 4 were assigned to high, moderate, marginal and unsuitable classes, respectively. The next step involved assigning a new value for the degree of membership of every attribute to each index at each level. All these new values ranged from 0 to 1 where 0 indicated a poor fit and 1 indicated a perfect fit for membership of an attribute to a given index.. The membership degree function (Fi) of each assessment factor is given as follows:

The membership degree of soil organic matter, total nitrogen, available phosphorus and available potassium to the soil nutrient status criterion were calculated using the following equations:where and are the lower and upper limits of the value range of a factor.

The membership degree of slope gradient was given by:

The membership degree of physical clay content was given by:where means the best physical clay content in the study area.

The membership degree of altitude was given by:

The degree of membership of soil water content and bulk density was given by:

Assessment of land suitability

The integrated classification of land suitability for the study area was obtained by integrating the combined weight () and the total degree of membership for different factors () as follows:

The maximum membership principle was adopted to define the comprehensive classification of land suitability. If , the comprehensive classification was .

Data

Field survey and soil sampling were conducted during August and September 2011. The land use status in 1993 and 2011 is presented in Table 3. When compared with 1993, there was significantly less arable land in 2011, whilst there had been increases in the extent of grassland, land that had been build on and land that was not in use. The main reason was the construction of four collieries during this period; this resulted in the transfer of the work force (local farmers) from traditional cultivation to the coal industry [30]. Based on land use and topography, we selected 15 assessment units: three for representing arable land, five for forestland and seven for grassland. Details of the assessment units are presented in Table 4.

[Figure omitted. See PDF.]

Table 3. Land use in the Liudaogou watershed (%) in 1993 and 2011.

https://doi.org/10.1371/journal.pone.0069498.t003

[Figure omitted. See PDF.]

Table 4. Detailed information about the assessment units.

https://doi.org/10.1371/journal.pone.0069498.t004

Altitude was measured in the field by double frequency RTK-GPS. Slope gradient was determined using a compass incorporating an inclinometer. Sampling transects were placed along the slope in each assessment unit. For each transect, 3 plots (15×15 m) were established for the field investigation and sampling. In each sampling plot, a 1.0 m long×0.7 m wide×0.5 m deep pit was dug to allow measurements of soil bulk density. Samples were collected using a 5.0 cm diameter by 5.0 cm long stainless cutting ring with which samples were collected from the top 20 cm of the soil profile. Soil samples in the ring were extracted and the roots were carefully removed by hand. The soils were dried at 105°C to constant weight and bulk density was calculated by dividing the dry mass by the known volume. Soils were also collected with a soil auger (4 cm-diameter) to estimate profile water content in 10 cm layers to a depth of 40 cm. Soil water content was determined by oven-drying samples at 105°C.

Composite samples of about 1 kg from each assessment unit were collected and then air-dried and ground to pass through 1.00 mm and 0.25 mm nylon screens prior to laboratory analysis. The sieved soil was used to measure soil particle composition (physical clay content) and soil organic matter (SOM) , respectively. Particle composition was measured by the laser diffraction technique using a MasterSizer 2000 (Malvern Instruments, Malvern, England), equipped with a low-power (2 mW) Helium-Neon laser with a wavelength of 633 nm as the light source. Soil organic carbon (SOC) was measured by an Elementar Vario EL element analyzer and soil organic matter (SOM) was obtained by multiplying SOC content values by 1.723.

Soil total nitrogen (TN) was measured following the Kjeldahl digestion method [31]. Total phosphorus (TP) was determined colorimetrically after wet digestion with sulfuric acid and perchloric acid, and available phosphorus (AP) was extracted following the Olsen bicarbonate extractable P method [32]. Available K was determined by the ammonium acetate extraction method [33].

Results and Discussion

Suitability assessment for arable land

Standardized evaluation factors for arable land.

Arable land suitability was divided into 4 classes, i.e. highly suitable (S1), moderately suitable (S2), marginally suitable (S3) and unsuitable (N). All factors were standardized to a uniform rating scale for ease of analysis. Standardized values for specific factors required specific decision rules in the form of thresholds for each factor. Various statistical and empirical guidelines from the related national–provincial codes and literature were used to determine the boundary values. In addition, twelve experts were invited to act as the decision makers to guarantee the reliability of boundary values. The experts involved were ecologists, pedologists, geographers, forest and grass managers, land resource experts and environmental protection experts; all have undertaken related studies of vegetation restoration and ecological– environmental issues in the study area and all are familiar with the Loess Plateau. Thus, they were able to suggest reasonable values for the thresholds. As a general guideline, there was a positive correlation between the value assigned and suitability. The class boundaries and standardized measurements applied for each factor are shown in Table 5.

[Figure omitted. See PDF.]

Table 5. Land characteristics, thresholds and degree of suitability for arable land.

https://doi.org/10.1371/journal.pone.0069498.t005

Distribution of suitable cropland and associated statistics.

Based on the method described above, suitable cropland classes were calculated and the main limiting factors were identified for each assessment unit (Table 6). Assessment units 1 and 3 in arable land and unit 14 in grassland were considered to be of marginal suitability, assessment unit 2 in arable land was classified as being of moderate suitability, and the other units from forest and grasslands were classified as being unsuitable for crops. Generally, an area of 0.40 km2, accounting for 6.7% of the total area of the watershed (not including the built-up and unused lands), was classified as being moderately suitable for crops, and 1.19 km2 (20.0%) as marginally suitable for crops (Table 7). In total, the area suitable for crops comprised about 26.7% of the watershed, so the area unsuitable for crops accounted for 73.3% of the watershed. The main factors restricting cropland development were soil physical properties and soil nutrients. First, the soil organic matter content and total nitrogen were generally low. This was particularly true in forest and grasslands, where the total nitrogen ranged between 0.05 and 0.35 g/kg (on average 0.17±0.09 g/kg), and was significantly lower than the generally reported values for soils (0.48–0.91 g/kg; on average 0.71) in the center of the Loess Plateau [34]. Soil organic matter was in the range 1.99–13.84 g/kg (on average 6.99±4.05 g/kg), and lower than the generally reported values for soils (7.0–15.0 g/kg; on average 11.0) in the center of the Loess Plateau [35]. Secondly, soil texture and soil water content obviously affected the land's suitability to grow crops in the watershed. In forestland in particular, the clay content was significantly lower than that of crop and grasslands, suggesting that planting trees adversely affected the soil structure in the study area. Our data were consistent with those presented in previous studies [26], [30] indicating that species such as Populus simonii can cause soil degradation. The watershed is located in the ecotone between arid and the semi-arid regions, and soil water is always a limiting factor for plant growth because of the low levels of precipitation (about 437 mm) and high potential evapotranspiration (1800 mm). In forestland, soil water content ranged between 2.00% and 4.46%, which was close to the wilting point. Grasslands also exhibited soil water deficits compared with the croplands, although soil water content was higher than that in forestlands. Previous studies have shown that planting P. simonii (a tree species) and Medicago sativa (alfalfa, a forage legume) could cause large reductions in soil water content due to their high water consumption, and this could result in the formation of dry soil layers on land where they have been planted [36]–[37]. Usually, soil drying due to high plant water use adversely affects soil physical properties, for example, reducing aggregate stability and soil surface roughness [38]–[39]. Thus, soil water deficit and soil physical properties may have a compounding effect.

[Figure omitted. See PDF.]

Table 6. Cropland suitability and limiting factors for each assessment unit.

https://doi.org/10.1371/journal.pone.0069498.t006

[Figure omitted. See PDF.]

Table 7. Areas of land suitable for growing crops in the Liudaogou watershed.

https://doi.org/10.1371/journal.pone.0069498.t007

It is noteworthy that available phosphorus was in the range 2.46–5.87 mg/kg (on average 4.29±1.39 mg/kg) in forest and grass lands, and was, therefore, very much higher than generally reported for soils (0.91–1.76 g/kg; on average 1.10) in the center of the Loess Plateau. Available potassium was in the range 85.04–297.40 mg/kg (on average 165.60±78.49 mg/kg), which is also significantly higher than generally reported for soils (77.2–170.0 mg/kg) in the center of the Loess Plateau [35]. Therefore, these two nutrients were not the limiting factors for cropland development in the watershed.

Suitability assessment for forest and grassland

Standardized evaluation factors for forest and grassland.

As for the arable land suitability assessment, forest and grassland suitability was also divided into 4 categories: highly suitable (S1), moderately suitable (S2), marginally suitable (S3) and unsuitable (N). The class boundaries and standardized measurements applied for each factor are shown in Table 8. Altitude is not a limiting factor for forest and grassland development.

[Figure omitted. See PDF.]

Table 8. Land characteristics, thresholds and degree of suitability for forest and grasslands.

https://doi.org/10.1371/journal.pone.0069498.t008

Distribution of forest and grassland suitability.

Only assessment unit 6 in the forestland and unit 11 in the grassland were classified as being unsuitable for grassland (Table 9). The limiting factor for unit 6 was the low soil water content resulting from the high water consumption by tree species, and the limiting factor for unit 11 was soil texture because of its very low clay content (Table 4). Overall, an area of 0.79 km2, accounting for 13.3% of the total area of the watershed, was classified as highly suitable for grassland, 1.18 km2 (20.0%) was classified as being moderately suitable for grassland, and 3.16 km2 (53.4%) as being marginally suitable for grassland (Table 10). In total, the area suitable for grassland was about 86.7% of the watershed, with the remaining 13.3% being unsuitable. This result is consistent with a previous study suggesting that development of vegetation in the watershed should focus on grassland [30]. Hou et al. [27] reported that the watershed is located within the boundary of the 400 mm rainfall isoline, representing the demarcation between cropping and pastoral regions. From the perspective of the climate, the watershed is more suitable for the growth and development of herbaceous plants than for trees. It is also notable that the area is not entirely appropriate for grassland, with 53.4% being classified as only marginally suitable. Besides soil water content and soil texture, total nitrogen may be a subsidiary factor limiting the suitability of this land for use as grassland – it is present at low levels in forest and grasslands (Table 11).

[Figure omitted. See PDF.]

Table 9. Grassland suitability and limiting factors for each assessment unit.

https://doi.org/10.1371/journal.pone.0069498.t009

[Figure omitted. See PDF.]

Table 10. Areas suitable for use as grassland in the Liudaogou watershed.

https://doi.org/10.1371/journal.pone.0069498.t010

[Figure omitted. See PDF.]

Table 11. Forestland suitability and limiting factors for each assessment unit.

https://doi.org/10.1371/journal.pone.0069498.t011

With respect to forestland suitability, assessment units 2 and 3 in the arable land were classified as being moderately suitable, assessment unit 1 in the arable land and units 13 and 14 in the grassland were considered marginally suitable, and other units in the forest and grasslands were unsuitable for use as forest. Overall, an area of 0.79 km2, accounting for 13.3% of the total area of the watershed, was classified as moderately suitable for forestland and 1.19 km2 (20.0%) as marginally suitable for forestland; the remaining 3.95 km2, accounting for 66.7% of the watershed was unsuitable for forest (Table 12).

Like cropland, the main factors restricting the establishment of forest were soil water content, organic matter and total nitrogen. This means soil water and nutrition conditions were not suitable for widespread afforestation. The slope gradients of the arable and grassland assessment units that were suitable for forest were all below 10u, indicating that slope gradient may be a subsidiary limiting factor with respect to forests. Usually, steep slopes are associated with low soil water content, with increasing soil erosion and decreasing infiltration [40]–[42].

[Figure omitted. See PDF.]

Table 12. Areas of land suitable for forest growth in the Liudaogou watershed.

https://doi.org/10.1371/journal.pone.0069498.t012

Implications for the Sloping Land Conversion Program

On the Loess Plateau, rapid population growth and an economic boom, coupled with severe soil erosion and desertification, have led to deterioration of the natural environment and a reduction in biodiversity [43], [4]. The Sloping Land Conversion Program is, therefore, considered to be a necessary step for increasing vegetation cover and restoring ecosystem service functions, thus promoting sustainable environmental and economic development. However, during implementation of the SLCP, tree planting was overemphasized, resulting in inappropriate planting schemes [34]ü. L et al. [8]; reported that the large areas converted from farmland to woodland have resulted in decreased regional water yield as the climate warms and dries on the Loess Plateau. Successful ecological rehabilitation programs have, thus, been largely dependent on innovative ecosystem management systems and technical support. Land suitability evaluation indicates that the Liudaogou watershed is suited to grass growth and is unsuitable for forest growth on the sloping land this is due mainly to water and nitrogen deficits. In fact, there is still 1.79 km2% (26 of the total area) of forestland in the watershed, with Populus simonii and Pinus tabulaeformis that were mainly planted in the late 1970s and after 1999 [27], [30]. Hou [27] and Yang et al. [44] reported that the forest in the Liudaogou watershed could not function as a stable forest ecosystem. The height of the 20-year-old Populus simonii–– trees averaged about 46 m (the smallest just 2 m), the diameter at breast height (DBH) was about 56 cm, and the volume of timber was just 0.0031 m3 per individual. In general, the growth of Populus simonii stopped after about 15 years. Under normal circumstances, the height of a 20-year-old Populus simonii≥ tree would be more than 20 m with DBH 15 cm. Due to limited cover and poor growth, the Populus simonii“” forest is very restricted in its ecological function. It was originally planted as a part of a shelterbelt in the Three North Protective Forest Program, aiming to combat desertification and soil erosion in the Northwest of China [45]–[46]. However, according to a survey conducted by Hou et al. [47], the wind erosion depth in the Populus simonii–× forestland amounts to 1.914.68 cm/a, corresponding to 1.9104–×4.7104 m3;– soil loss per year. Compared with the natural grass vegetation, it does not prevent wind erosion, but greatly accelerates it. The main reasons for the stunted tree growth are: (1) there is insufficient soil water available to maintain normal growth rates (2) there is insufficient soil fertility (e.g. besides low total nitrogen, the Soil Organic Carbon Density is about 1.182.81 kg/m2–, significantly lower than the national mean values for soils of 11.5212.04 kg/m2) [48];– and (3) there is insufficient management because of the very limited economic value of forest local farmers are more willing to devote themselves to coal-mining which generates a much higher income.

%Although our forestland suitability assessment indicated that about 33.3 of total area is suitable for forest growth (Table 12), such areas are currently arable land and some grassland with gentle gradients (Table 4). Compared to tree species, crop and grass species have shallow roots and consume relatively little water [30];, so do not deplete soil water to any great depth however, the current arable lands are dammed, terraced fields that collect valuable rainfall because of the low soil erosion rates. Therefore, those arable lands and grasslands with gentle gradients were classified as being moderately or marginally suitable for forest development. We predict that planting trees in these areas would cause soil water depletion because of high water consumption combined with the low precipitation rates in the watershed.

%The grassland suitability assessment indicated that about 86.7 of the area of the watershed is suitable for grass growth, thus the SLCP in the watershed should focus on conversion of cropland to grassland. Wang et al. [17]>– reported that soil physical properties such as bulk density, hydraulic conductivity, mean weight diameter, and the stability of 1 mm macro-aggregates have been significantly ameliorated in the 020 cm soil layer under secondary natural grasslands, implying that natural grass (Stipa bungeana– Trin.) restoration is an appropriate and sensible approach to re-vegetation in the windwater erosion region of the northern Loess Plateau of China. Planting alfalfa (Medicago sativa) and korshinsk peashrub (Caragana korshinskii) has no effect on the soil physical conditions, and may even reduce bulk density and soil permeability because of their very high water consumption with deep and vigorous root systems. Therefore, care is required when selecting species for re-vegetation to promote shifts from arable land to grasslands in a more environmentally compatible manner. Non-native grassland species, such as the korshinsk peashrub and alfalfa, have high water requirements and should be avoided when choosing species to plant in the watershed.

Conclusion

Land use management involves complex decision-making that requires an understanding of many factors. This paper presents the AHP as a decision support tool for use when selecting an appropriate land use, which is an important issue for the sustainable development of the Liudaogou watershed. Based on the AHP, grassland was found to be the most suitable land-use type to support the aims of the Sloping Land Conversion Program in the watershed. Under the constraints of soil water shortage and nutrient deficit, crops and forests were considered to be inappropriate land uses in the study area, especially on sloping land. All forest should be converted to grassland because continued growth of trees will damage the soil water environment and increase desertification problems. When selecting species for re-vegetation, non-native grass species with high water requirements should be avoided so as to guarantee the sustainable development of grassland and effective ecological functioning. In the future, the area of abandoned cropland is likely to increase rapidly due to the government policy embodied in the Sloping Land Conversion Program. Our study provides local land managers and farmers with valuable information about the inappropriateness of growing trees in the study area along with some information on species selection for planting in the semi-arid area of the Loess Plateau.

Acknowledgments

We wish to thank graduate students Wang Jianguo, Shi Zhanfei and Wang Mei for their help in field survey and data analysis in the Laboratory. We also wish to thank Dr. David Warrington for his invaluable comments and suggestions on this paper.

Author Contributions

Conceived and designed the experiments: LW XBY. Performed the experiments: XBY LW. Analyzed the data: XBY. Contributed reagents/materials/analysis tools: XBY LW. Wrote the paper: LW XBY.

Citation: Yi X, Wang L (2013) Land Suitability Assessment on a Watershed of Loess Plateau Using the Analytic Hierarchy Process. PLoS ONE 8(7): e69498. https://doi.org/10.1371/journal.pone.0069498

References

1. World Bank (2001) China: air, land and water, environmental priorities for a New Millennium. World Bank, Washington, DC.

2. –SFA (2003) Sloping land conversion program plan. (20012010) (in Chinese).

3. Michael T, Bennett MT (2008) '?Chinas sloping land conversion program: Institutional innovation or business as usual Ecological economics 65: 699–711.

4. Wang L, Shao MA, Wang QJ, Gale WJ (2006) Historical changes in the environment of the Chinese Loess Plateau. &Environmental science policy 9: 675–684.

5. Ding CR (2003) Land policy reform in China: assessment and prospects. Land Use Policy 20: 109–120.

6. Uchida E, Xu JT, Rozelle S (2005) 'Grain for Green: cost-effectiveness and sustainability of Chinas conservation set-aside programme. Land Economics 81: 247–264.

7. Wang L, Wang QJ, Wei SP, Shao M, Li Y (2008) Soil desiccation for Loess soils on natural and regrown areas. Forest Ecology and Management 255 (7) 2467–2477.

8. üL Y, Fu B, Feng X, Zeng Y, Liu Y, et al. (2012) A policy-driven large scale ecological restoration: quantifying ecosystem services changes in the Loess Plateau of China. PLoS ONE 7 (2) e31782

9. Wei J, Cheng J, Li W, Liu W (2012) Comparing the effect of naturally restored forest and grassland on carbon sequestration and its vertical distribution in the Chinese Loess Plateau. PLoS ONE 7 (7) e40123

10. Yin RS, Zhao MJ (2012) —'Ecological restoration programs and payments for ecosystem services as integrated biophysical and socioeconomic processesChinas experience as an example. Ecological Economics 73: 56–65.

11. Wang XH, Lu CH, Fang JF, Shen YC (2007) Implications for development of grain-for-green policy based on cropland suitability evaluation in desertification-affected north China. Land Use Policy 24: 417–424.

12. McVicar TR, Li LT, Van Nie TG, Zhang L, Li R, et al. (2007) 'Developing a decision support tool for Chinas revegetation program: Simulating regional impacts of afforestation on average annual streamflow in the Loess Plateau. Forest Ecology and Management 251: 65–81.

13. Yang RJ, Fu BJ, Liu GH, Ma KM (2004) Research on the relationship between water and eco-environment construction in Loess Hilly and Gully Region. Chinese Journal of Environmental Science 25 (2) 37–42 (in Chinese with an English abstract)..

14. Yang WZ, Tian JL (2004) Essential exploration of soil aridization in the Loess Plateau. Acta Pedologica Sinica 41 (1) 1–6 (in Chinese with English abstract)..

15. McVicar TR, Van Niel TG, Li LT, Wen ZM, Yang QK, et al. (2010) 'Parsimoniously modelling perennial vegetation suitability and identifying priority areas to support Chinas re-vegetation program in the Loess Plateau: Matching model complexity to data availability. Forest Ecology and Management 259: 1277–1290.

16. Wang L, Wei SP, Horton R, Shao M A (2011) Effects of vegetation and slope aspect on water budget in the hill and gully region of the Loess Plateau of China. Catena 87: 90–100.

17. Wang L, Wei SP, Shao HB, Wu YJ, Wang QJ (2012) Simulated water balance of forest and farmland in the hill and gully region of the Loess Plateau of China. Plant Biosystems dx.doi.org/10.1080/11263504.2012.709198.

18. Saaty TL (1977) A scaling method for priorities in hierarchical structures. Journal of Mathematical Psychology 15: 234–281.

19. Saaty TL (1980) The analytical hierarchy process. McGraw Hill, New York.

20. Kurttila M, Pesonen M, Kangas J, Kajanus M (2000) Utilizing the analytic hierarchy process (AHP) in SWOT analysis-a hybrid method and its application to a forest-certification case. Forest Policy and Economics 1 (1) 41–52.

21. Ayalew L, Yamagishi H, Marui H, Kanno T (2005) Landslides in Sado Island of Japan: Part II. GIS-based susceptibility mapping with comparisons of results from two methods and verifications. Engineering Geology 81: 432–445.

22. Cay T, Uyan M (2012) Evaluation of reallocation criteria in land consolidation studies using the Analytic Hierarchy Process (AHP). Land Use Policy 30: 541–548.

23. Mardle S, Pascoe S (2004) Management objective importance in fisheries: an evaluation using the Analytic Hierarchy Process (AHP). Environmental Management 33 (1) 1–11.

24. Pascoe S, Dichmont CM, Brooks K, Pears R, Jebreen E (2013) Management objectives of Queensland fisheries: Putting the horse before the cart. Marine Policy 37: 115–122.

25. Li M, Li ZB, Liu PL, Yao WY (2005) –Using cesium-137 technique to study the characteristics of different aspect of soil erosion in the WindWater Erosion Crisscross Region on Loess Plateau of China. Applied Radiation and Isotopes 62: 109–113.

26. Wang L, Mu Y, Zhang QF, Jia Z K (2012) Effects of vegetation restoration on soil physical properties in the wind-water erosion region of the northern Loess Plateau of China. CLEAN - Soil, Air, Water 40 (1) 7–15.

27. Hou QC (1994) Comprehensive analysis on natural conditions and environmental harnessing in the experimental Area. Memoir of NISWC, Academia Sinica and Ministry of Water Resources 18: 136–143 (in Chinese with English abstract)..

28. Fahey B, Jackson R (1997) Hydrological impacts of converting native forests and grasslands to pine plantations, South Island, New Zealand. Agricultural and Forest Meteorology 84: 69–82.

29. Chakraborty S, Banik D (2006) Design of a material handling equipment selection model using analytic hierarchy process. The International Journal of Advanced Manufacturing Technology 28: 1237–1245.

30. Wang L, Zhang QF, Wei SP, Wang QJ (2009) Vegetation restoration model in a watershed of a coal mining area in the water and wind erosion crossing zone of the Loess Plateau. Journal of Beijing Forestry University 31 (2) 36–43.

31. ——–Bremmer JM, Mulvaney CS (1982) Nitrogentotal. In: Page AL, Miller RH, Keeney DR (Eds.), Methods of soil analysis, Part 2Chemical and microbiological properties. ASA-SSSA, Madison, WI, pp. 595624.

32. –Olsen SR, Sommers LE (1982) Phosphorous. In: Page AL, Miller RH, Keeney DR (eds) Methods of soil analysis Part 2, Chemical and microbial properties. Agronomy Society of America, Agronomy Monograph 9, Madison, Wisconsin, pp 403430.

33. –Pratt PF (1965) Potassium. In: Black, C.A. (Ed.), Methods of Soil Analysis. Part 2. Chemical and Microbiological Properties. Am. Soc. of Agron, Inc., Madison, p. 10221030.

34. Wang XL, Guo SL, Ma YH, Huang DY, Wu JS (2007) Effects of land use type on soil organic C and total N in a small watershed in loess hilly-gully region. Chinese Journal of Applied Ecology 18 (6) 1281–1285 (in Chinese with English abstract)..

35. Bai WJ, Jiao JY, Ma XH, Jiao F (2005) Soil environmental effects of artificial woods in abandoned croplands in the Loess Hilly-gullied region. Journal of Arid Land Resources and Environment 19 (7) 135–141 (in Chinese with English abstract)..

36. Jiang N, Shao MA, Lei TW (2007) –Soil water characteristics of different typical land use patterns in WaterWind Erosion Interlaced Region. Journal of Beijing Forestry University 29: 134–137 (in Chinese with English abstract)..

37. Huo Z, Shao MA, Horton R (2008) –Impact of gully on soil moisture of shrubland in WindWater Erosion Crisscross Region of the Loess Plateau. Pedosphere 18: 674–680.

38. Reid JB, Goss MJ (1982) Suppression of decomposition of C-labeled plant roots in the presence of living roots of maize and perennial ryegrass. Journal of Soil Science 33: 387–395.

39. Grant CD, Dexter AR (1987) Generation of microcracks in moulded soils by rapid wetting. Australian Journal of Soil Research 27 (1) 169–182.

40. Tang KL, Zhang KL, Lei AL (1998) “”Research on up-limit slope of grain for green in the plowlands of Loess Hilly and Gully Region. Chinese Science Bulletin 43 (2) 200–203 (in Chinese)..

41. Luk SH, Cai Q, Wang GP (1993) Effects of surface crusting and slope gradient on soil and water losses in the hilly Loess region, North China. Catena 24: 29–45.

42. Ng CWW, Shi Q (1997) A numerical investigation of the stability of unsaturated soil slopes subjected to transient seepage. Computers and Geotechnics 22 (1) 1–28.

43. Li WH (2004) Degradation and restoration of forest ecosystems in China. Forest Ecology and Management 201: 33–41.

44. Yang G, Song YX (1993) Analysis on the current situation of forestry and forestry developing strategy in Liudaogou Watershed of Shenmu Experimental Area. Memoir of Northwestern Institute of Soil and Water Conservation 18: 106–112.

45. Fang JY, Chen AP, Peng CH, Zhao SQ, Ci LJ (2001) Changes in forest biomass carbon storage in China between 1949 and 1998. Science 292: 2320–2322.

46. Wang XM, Zhang CX, Hasi E, Dong ZB (2010) ?Has the Three Norths Forest Shelterbelt Program solved the desertification and dust storm problems in arid and semiarid China Journal of Arid Environments 74 (1) 13–22.

47. Hou QC, Wang YK, Yang G (1996) Several problems of vegetation construction in criss2cross belt of bind2water erosion. Bulletin of Soil and Water Conservation 16 (5) 36–40.

48. Li YY, Shao MA, Zheng JY, Li QF (2007) Impact of grassland recovery and reconstruction on soil organic carbon in the northern Loess Plateau. Acta Ecologica Sinica 27 (6) 1–9.