Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

The Occurrence, Sources and Spatial Characteristics of Soil Salt and Assessment of Soil Salinization Risk in Yanqi Basin, Northwest China

  • Zhang Zhaoyong,

    Affiliations State Key Laboratory of Desert and Oasis Ecology, Xinjiang Institute of Ecology and Geography, Chinese Academy of Sciences, Urumqi, China, University of the Chinese Academy of Sciences, Beijing, China

  • Jilili Abuduwaili ,

    [email protected]

    Affiliation State Key Laboratory of Desert and Oasis Ecology, Xinjiang Institute of Ecology and Geography, Chinese Academy of Sciences, Urumqi, China

  • Hamid Yimit

    Affiliation Key Laboratory of Xingjiang Arid Land Lake Environment and Resource, Xinjiang Normal University, Urumqi, China

Abstract

In order to evaluate the soil salinization risk of the oases in arid land of northwest China, we chose a typical oasis-the Yanqi basin as the research area. Then, we collected soil samples from the area and made comprehensive assessment for soil salinization risk in this area. The result showed that: (1) In all soil samples, high variation was found for the amount of Ca2+ and K+, while the other soil salt properties had moderate levels of variation. (2) The land use types and the soil parent material had a significant influence on the amount of salt ions within the soil. (3) Principle component (PC) analysis determined that all the salt ion values, potential of hydrogen (pHs) and ECs fell into four PCs. Among them, PC1 (C1-, Na+, SO42-, EC, and pH) and PC2 (Ca2+, K+, Mg2+and total amount of salts) are considered to be mainly influenced by artificial sources, while PC3 and PC4 (CO3- and HCO32-) are mainly influenced by natural sources. (4) From a geo-statistical point of view, it was ascertained that the pH and soil salt ions, such as Ca2+, Mg2+ and HCO3-, had a strong spatial dependency. Meanwhile, Na+ and Cl- had only a weak spatial dependency in the soil. (5) Soil salinization indicators suggested that the entire area had a low risk of soil salinization, where the risk was mainly due to anthropogenic activities and climate variation. This study can be considered an early warning of soil salinization and alkalization in the Yanqi basin. It can also provide a reference for environmental protection policies and rational utilization of land resources in the arid region of Xinjiang, northwest China, as well as for other oases of arid regions in the world.

Introduction

Soil salinization is a global problem and it is a potential environmental problem in all continents with the exception of unassessed Antarctica. Soil saline levels are found within a wide range, and soil salinization occurs in much of the waterfront, arid and semi-arid zones of more than 100 countries and regions [1][3]. According to statistics done by the United Nations Educational, Scientific and Cultural Organization (UNESCO), and the Food and Agriculture Organization (FAO), salinized soil covers an area of about 9.543×106 km2 on Earth [4][6]. In China alone, the area of salinized soil is about 3.693×105 km2, which accounts for about a third of the total arable land [7][8]. The area of salinized soil in the oasis basin of Xinjiang in northwest China is about 1.05×104 km2, which accounts for 33.4% of the total land in this area and research has found that the salinity of this area is trending upwards [9], [10]. Soil salinization restricts agricultural development, especially when sustainable agricultural development and environmental quality improvement strategies are being considered. Studies have found that when the salt ions in the soil attain 8 g.kg−1, they can greatly harm and even kill crops in farmland [11], [12].

In oases of arid regions of northwest china, the environment is so weak, including a lack of precipitation and high envapotion, that economic activities such as fishing, agriculture, forestry and grassland farming of the oases have been strongly limited, especially for agriculture [13][15]. Therefore, it is necessary to identify the distribution characteristics, sources of the soil properties, such as salt ions, potential of hydrogen (pH) and electrical conductivity (EC), and also the status and causes of salinization of the land of the oasis, in order to provide a scientific basis for protection of the soil that sustains land plants.

Multivariate analyses and other statistical methods have been widely applied in studies to determine the sources of elements found in soil, such as total soil salt content and heavy metals [16][18]. The spatial variation model and spatial distribution are used to make a hazard risk map of soil salt properties in regions of interest. Correlation analysis, principal component (PC) analysis and cluster analysis are classic methods used to identify the natural and man-made sources of salt ions and to simplify data. Additionally, use of the comprehensive index results in a class of data with high correlations that better reflect the associations between the data.

The Geostatistical Analyst is based on GIS technology [19], [20]. Among these, the ordinary kriging is the most widely used one in the study of soil salt distributions [21], [22], [29]. In recent years, the Geostatistical Analyst method has been used in the field of hydrology and water resources, including studies of groundwater pollution risk, water potential research and spatial distribution of soil salinization in arid land [23][25], [28], [32]. Since the 1990s, geostatistical methods have been widely used to study spatial variability characteristics of soil salt properties (salt ions, EC and pH). Sylla et al. [26] studied the spatial variation characteristics of soil salt content of an agricultural ecosystem under different scales in West Africa. Ammari et al. [27] studied the soil salinity changes in the Jordan Valley and the potential threat against the sustainable irrigated agriculture. In China, Bai et al. [30] researched the spatial variation characteristics and composition of soil salt content in Huang Huai Hai plain, northeast China and found that the average influential range of the soil salt content was higher than 200 km, indicating the salt content of the soil is mainly gathered in a large area of the regions. In Xinjiang in northwest China, Lin et al. [33] researched the spatial variation characteristics of soil salt in the Wei Gan He irrigated area and found the agricultural irrigation has resulted in serious soil salinization in this area and these deserve serious attention.

In arid regions, oases in basins are the main places where humans live and life can survive [31]. Therefore, it is important to understand the spatial distribution characteristics of the soil salt properties, including total salt content, salt ions, pH and EC. A quantitative grasp of soil salinity levels could serve as a reference and a basis for maintaining soil quality, which would help to effectively control the human pollution and develop the regional economy in a reasonable and orderly fashion [34], [35]. However, previous research has focused on rapidly developing areas, such as coastal plains and large irrigation areas in eastern china and elsewhere of the world with the purpose of assessing land usability, environmental effects and soil salinization risk [36][38]. Since the 1990s, implementation of the “western development policy of China” has led to prodigious economic development in many oases in Xinjiang, and the agriculture in these regions has undergone rapid progress. However, the rational irrigation of the agriculture, lack of precipitation and high envapotion of these regions have negatively influenced soil salt properties, resulting in increased soil salt contents, ECs and pHs, which can result in serious soil salinization [39]. Unfortunately, research on the soil salt property distribution characteristics and soil salinization risk assessment in the oases of arid regions of northwest of China is lacking.

The Yanqi basin is a typical oasis in a basin in the southern Tianshan Mountains, Xinjiang in northwest China. Since the 1990s, both the implementation of the “western development policy of China” and the development policy made by the Xinjiang Province, China, have led to prodigious economic development in the Yanqi basin, but regional economic development and associated human activity have left the current ecological environment fragile [39], [40]. Together with economic development, the blind expansion of farmland and unrestrained surface water irrigation led to a rise in groundwater and an increase in soil salinization in the basin oasis. 64.12% of the area experienced mild soil salinization, 8.25% had moderate salinization, and 27.07% had severe salinization. Research has shown that excessive use of water resources by agriculture has made the soil salinization status severe and decreased agricultural production [41].

After a basic analysis of land use and soil parent materials types in the area, we created land use and soil type geological maps, and, using ArcGIS 10.0 software and combining the grid sampling method with 3S technology, we made sampling points to get soil samples across the whole area. We evaluated the soil salt properties in different land use types in the laboratory, and assessed the soil salinization status and the cause in the Yanqi basin. Then ordinary kriging of Geostatistical Analyst method was used to reveal the spatial distribution characteristics of the soil salt properties in this region. Then by combining these properties with the climate, precipitation, evaporation and temperature, we assessed the soil salinization risk of this area. From this we can provide helpful proposals to prevent the environmental risks that could lead to soil salinization in this area. This research can serve as a helpful reference for environmental protection in this region and for soil salinization prevention in arid regions of northwest China.

Materials and Methods

Study area

The area studied in this work is a desert basin oasis in the arid region of northwest China including four counties in the Yanqi basin: Yanqi County, Hejing County, Bohu County and Heshuo County. This region lies within the geographical coordinates of 85°50′–87°50′E and 41°40′–42°20′N with a length of about 85 km from north to south and width of 130 km from east to west, totaling an area of about 723100 km2. The terrain slopes up from the northwest down to the southeast. The northwest is mountainous and the south is low-lying desert that is 1050–2000 m above sea level. The western area has extensive intrusive rock and metamorphic rock from the Proterozoic era, Neoproterozoic and Cenozoic. Weathering of this rock results in brown earth soil, acidic rocky soil and an acidic soil skeleton in the west. The east is primarily made up of quaternary sediments, which form Takyic (Calcisols), Chemic (Phaeozems), Stagnic (Gleysols), Irragric (Anthrosols), Fragic (Arenosols), Eutric (Gleysols) and Yemic (Solonchaks) [42]. The area researched is in a continental desert climate temperate zone with an annual mean temperature of 14.6°C, 186 frost-free days per year, and 50.7–79.9 mm of annual precipitation. By calculating the potential evaporation (ET0) by the method of Hargreaves (1985) [43], we then got the annual average potential evaporation of this area as 2438.9 mm, the ≥10°C active accumulated temperature 3414.4–3694.1°C and an annual average relative humidity of 72%.

Soil sampling and analyses

In order to perform a basic analysis of the land use and soil type, geological maps were made of the study area using ArcGIS 10.0 software to lay out a grid of soil sampling points on a digital map of the Yanqi basin. All samples were acquired in July 2012 or July 2013 from a collection area. In order to best assess the ecological risk in the Yanqi basin, diverse land use types were encompassed in our study of salt ion distribution. Soil samples were collected at depths of 0–20 cm, where a hard plastic shovel was used to dig a vertical 20×20 cm soil profile. 1 kg uniform samples were collected, and then they were put into a clean cloth, numbered and sealed. The collection position, date, sample vegetation types and surrounding vegetation conditions of each sampling area were recorded. After the soil samples were taken back to the laboratory, they were air dried and impurities, such as plant residues and rocks, were removed. The samples were then pushed through a 20 mesh nylon sieve (0.84 mm) to eliminate the plant residue and stones. We then used agate to grind the soil samples through 100 mesh nylon sieves (0.25 mm) to prevent contamination and then stored the samples in plastic bottles [46].

Total soil salt content, soil salt ions, pH and EC tested are as follows: 50 g of ground sample were removed from the plastic bottles, dissolved in 250 ml of deionized water (CO2 has been removed) (1∶5, soil∶water) for 2 hours to fully dissolve salt ions contained in the soil. The samples were then put in a centrifuge tube, vibrated for 3 min with an oscillator and then centrifuged at a speed of 4,500–5,000 r·min−1. To get the supernatant prepared for analysis of total soil salt content, salt ions content, pH and EC, the method described by Lu (2000) was followed [44].

The total salt content of the soil was determined by gravimetry of the evaporation residue. First, the supernatant was absorbed in a porcelain dish, and hydrogen peroxide (H2O2) was used to ozidize organic matter. Then, the samples were boiled in a water bath at 105–110°C until it dried, and weighed. The drying quality of the residue is expressed as total salt content of the soil. The pH of the extracted supernatant was tested using a Potentiometric Titrimeter (G20, METTLER, and TOLEDO). Burette drive resolution was 1/20000. Mv/pH electrode measurement range was ±2000 mv. The ECs were tested using a Conductivity Meter (DDSJ-308A, Shanghai, China) with a measurement range of 0–1.999×105 µs/cm and a test error of ±0.5% (FS) ±1.

The supernatant was run through a 0.45 µm drainage cellulose acetate membrane. Then, the cation content (K+, Na+, Ca2+, Mg2+) of the solutions was determined using an inductively coupled plasma atomic emission spectrometer (Vista MPX, Varian, USA). The anion content (Cl-, SO42-, CO32-, and HCO3-) of the solutions was determined using an Ion Chromatograph (ICS-90, Dionex, USA). All tests were conducted using the following protocol: a standard solution was prepared for Na+, K+, Ca2+, Mg2+, Cl-, SO42-, CO32- and HCO3-. The salt ion content was determined by comparing each sample to the standard solution of known concentration. The standard solutions used for the salt ions in this study were national level standard material (Gss series, China). The coefficient of the best fitting curve was determined by the testing equipment based on the standard material and then the amount of the salt ions (Na+, K+, Ca2+, Mg2+, Cl-, SO42-, CO32-, and HCO3-) in the liquid supernatant was tested. After all the samples had been tested for their salt ion content (Na+, K+, Ca2+, Mg2+, Cl-, SO42-, CO32-, and HCO3-), we chose approximately 20% for retesting and found that 97.3% of the results were repeatable, inspiring confidence in the original data. After all the total salt ion contents (Na+, K+, Ca2+, Mg2+, Cl-, SO42-, CO32-, and HCO3-) of the solution were determined, we recalculated them from unit of µg/ml into mg/g (g/kg) using the method described by Bao (2005) [45]. To prevent contamination during the testing process, all glassware was soaked in 5% HNO3 for 24 hours, rinsed and then dried.

Statistical analyses

Descriptive and multivariate statistical analysis.

Descriptive statistical methods were used to analyze the range, mean, median, standard deviation, coefficient of variation, kurtosis and skewness of the total salt content, each salt ion, SAR, pH and EC of the soil samples. Correlation analysis, PC analysis and cluster analysis of the classic multivariate statistical method were used to process data and identify the soil salinity. Single factor analysis of variance (ANOVA) was used to analyze the differences in the amount of salt ions between different land use types. These analyses were all processed using the software SPSS 19.0.

Ordinary kriging method.

Ordinary kriging (OK) is a commonly used linear spatial interpolation method that estimates variables at unsampled locations by using information from neighboring points and assigning weights to these points based on their distance from the point and the spatial variability structure. The OK method can be formulated as(1)where is the OK estimation at an unsampled location , n is the number of samples in a search neighborhood, and are the weights assigned to the ith observation . Weights are assigned to each sample such that the estimation or kriging variance is minimized and the estimates are unbiased [47]. Weights are determined after computing a semivariogram that models spatial correlation and covariance structure between data points for each variable using Eq. 1 [48].(2)where is the semivariance between two observation points and separated by a distance h, and N is number of observation pairs at the distance h.

Soil salinization evaluation criteria used in this research as under below (Table 1).

thumbnail
Table 1. Indicators used for risk assessment of soil salinization.

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

Results and Discussion

Descriptive statistical analysis of soil salinity

The descriptive statistics concerning the soil properties in Yanqi basin in Table 2 show that the maximum and average values of HCO3-, Ca2+, CO3-, Na+, Mg2+, K+, SO42-, Cl-, SAR, total amount of salts, EC, and pH were 0.98(0.19) g.kg−1, 1.95(0.75) g.kg−1, 0.85(0.46) g.kg−1, 2.43(1.19) g.kg−1, 1.89(0.58) g.kg−1, 2.13(0.69) g.kg−1, 1.58(1.14) g.kg−1, 2.36(0.98) g.kg−1, 33.241(10.51), 14.77(9.73) g.kg−1, 1.39(0.95) dS.m−1, and 8.55(8.15), respectively. Within the analysis of the soil samples, a large amount of variation occurred, suggesting that the sources and influencing factors of the salt properties in Yanqi basin are complex. This work found that the main salt ions accounted for 84.41% of the total salt content and were K+, Ca2+, Na+, Cl-, Mg2+ and SO42. Meanwhile the amount of HCO3- and CO3- was very low. The pH of the soil of the study areas ranged from 7.85 to 8.55, which had only a small variation. However, there was a dramatic change within the total salt content of the soil, ranging from 1.16 to 14.77 g.kg−1.

thumbnail
Table 2. Descriptive statistics of the soil salt properties from Yanqi basin.

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

The coefficient of variation is the ratio between the standard deviation and average, and it can be used to compare different dimensions of indicators. The coefficients of variation of HCO3-, CO3-, Na+, Mg2+, SO42-, Cl-, EC, and pH were 35.37%, 32.08%, 23.01%, 20.07%, 19.63%, 12.94%, 21.27%, and 30.32%, respectively, and were of medium variation (10%<CV<100%). However, the coefficients of variation of SAR, the total amount of salts, Ca2+ and K+ were 148.56%, 126.73%, 191.67% and 226.25%, and, therefore, had high levels of variation (CV>100%)[52]. In particular, Ca2+ and K+ had higher coefficients of variation as compared to the other elements. From the perspective of skewness, the values of these ten soil properties are ordered as K+>Ca2+>total amount of salts>Na+>SAR>Cl->EC>Mg2+>CO3->SO42->pH>HCO3-.

The differences in soil salt properties from different land use types and soil parent materials

Land utilization types and soil parent materials are the main examples of human activity and geological background that influence the salt ion content of soil. We analyzed the relationship between soil salt properties, human activity and geological background to further explore the distribution characteristics and sources of the soil salt properties of the Yanqi basin. The salt properties of the soil from each land use type and soil parent material of Yanqi basin are in Table 3. The land use types and soil parent materials of Yanqi basin are shown in Fig. 1. The analysis of these data suggests that the manner in which land is used has a significant influence on the amount of Ca2+ and K+. In farmland, the average content of Ca2+ was 1.83 g.kg−1, K+ was 1.71 g.kg−1, and the total amount of salts was 10.88 g.kg−1. This was significantly higher than in the other land use types including areas of urban construction, forest, grassland, urban construction areas and desert. Meanwhile the maximum average values of Na+, Mg2+, SO42-, Ca2+, the total amount of salts, and K+ found in farmland were higher than in grassland, desert, and areas of urban construction. The variance test attained a significant level of 0.05, indicating these elements and their distribution are mainly controlled by human activity. The activities that had a significant influence on these soil salt properties include agricultural activities, such as irrigating, fertilizing and farming. The CV calculated indicates that the pH, EC, total amount of salts and all the soil salt ions measured belong to medium variability categories (10–100%). Among the five land use types, the differences between these groups were small and, therefore, the classes of element enrichment were not obvious.

thumbnail
Figure 1. Land use types and parental material pattern in Yanqi basin.

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

thumbnail
Table 3. Statistical parameters of the soil salt ions found at 0–20 cm depth within the investigated land use and land cover categories of the study area.

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

This research also found differences in the amount of organic soil salt ions. For example, HCO3- reached a maximum in lacustrine deposits of 0.74 mg.kg−1. The maximum average values of Na+, Mg2+, K+, SO42-, total amount of salts, and Cl- were found in coarse crystalline rock weathered material, sandy shale of weathered material and lacustrine deposits, while the maximum average values of CO3-, Ca2+, and EC were found to be significantly higher in diluvial material than in sandy shale or coarse crystalline rock from weathered material or lacustrine deposits. For the five land use types, the class of element enrichment was not obvious as the differences between these groups were small.

R2 represents the ratio of the sum of squares in groups and the total error of the sum of squares. It reflects the contribution of different factors on the soil salt properties [53]. For this study, the land use types explained the variances in Ca2+ (37.85%), K+ (71.45%) and total amount of salts (56.71%), which were higher than that of the soil parent material. This demonstrates that the way land was used played a major role in the accumulation of Ca2+, K+ and total amount of salts in Yanqi basin. This analysis also found that the variances of HCO3-, CO3-, Mg2+, Na+, SO42-, and Cl- in the soil parent material are higher than those of the land use factors (Table 4), indicating that the soil parent material played a major role in the accumulation of these elements. However, there was little difference in the variances of EC and pH, indicating that these soil salt properties were mainly influenced by land use and soil parent material. Overall, the R2 analysis fits well with the results of the multivariate statistical analysis.

thumbnail
Table 4. The correction matrix of soil salt properties in Yanqi basin.

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

Multivariable statistics

Correlation analysis.

Table 4 shows the Pearson correlation coefficients between the soil salinity variables. There is a significant correlation of 0.96 (P<0.01) between the total amount of salt cations and Na+, as well as the total amount of salt cations and Ca2+ at 0.51 (P<0.05) in the soil in the Yanqi basin, but there is no significant correlation with other salt cations. Furthermore, we found a significant correlation between the amount of salt anions and SO42- of 0.94 (P<0.01), and the correlation coefficients between the total amount of salt anions and Cl- or CO32- are 0.55 and 0.49 (P<0.05), respectively. This indicates that SO42- was the primary salt anion, Cl- was the secondary and CO32- was the tertiary. Together, the correlation coefficients between the total amount of salts and Cl-, and the total amount of salts and SO42- are 0.57 and 0.92, respectively (P<0.01), indicating that the main types of salt in the soil were sulfate and chloride. This study also found that there are close correlation coefficients between the EC and Na+, the EC and Cl-, and the EC and SO42- content, where the coefficients are 0.98, 0.66, and 0.87, respectively (P<0.01). Additionally, the correlation coefficient between the soil EC and CO32- is 0.48, which is also significant (P<0.05). Overall, a significant correlation was found between pH and EC, and total amount of salts (TS) and salt anions (Na+, Ca2+, Ca2+, Ca2+, Cl-, SO42-, HCO3-, and CO32-) (P<0.01). Further analysis shows that the correlation coefficient between the pH and SO42- is 0.67 and between pH and Cl- is 0.58 (P<0.01), indicating the pH of the soil is primarily influenced by the SO42- and Cl- content.

Principal component analysis and clustering analysis.

All of the elements studied were found to fall into four PCs (Table 5) with a cumulative variance of 95.72%. This is a reflection of their sources and main influences, particularly for the ten indicators of Cl-, Na+, EC, pH, K+, Mg2+, Ca2+, HCO3-, CO32-, total amount of salts, and SO42-. Among these, the variance contribution rate of the first PC (Cl-, Na+, SO42-, EC and pH) was 33.75% and the second PC (Ca2+, K+, Mg2+, and total amount of salts) was 28.54%. The primary contribution was by agricultural development, in particular from herbicide application and acid salt fertilizer [54], [55]. The variance contribution rates of the third (CO3-) and fourth PCs (HCO32-) were at 18.25% and at 15.18%, respectively. Upon combining the sampling sites within their respective land use types, it was found that the samples that had a high content of salts were often taken from the desert, grassland and forest, and appear to have originated from natural sources. This analysis also showed there were larger loads of Ca2+ in the first PC, SO42- in the second PC and K+, total salt content in the fourth PC, indicating these elements were influenced by both artificial and natural sources. We used clustering analysis to see if the results are consistent with the results from the PC analysis (Fig. 2). All the elements were classified into four categories: the first category consisted of Cl-, Na+, SO42-, pH, and EC; the second of Mg2+, Ca2+, K+, and the total amount of salts; the third of HCO3; and the fourth of CO32-.

thumbnail
Figure 2. Clustering tree of soil salt properties of Yanqi basin.

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

thumbnail
Table 5. Factors matrix of soil salt properties from Yanqi basin.

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

Spatial distribution of soil salt ions in Yanqi basin.

The spatial dependence of the soil salt properties was determined by semivariance analysis in order to quantitate the spatial variability. The parameters of the semivariogram included model type, nugget, sill and effective range. The nugget value (C0) represents the random variation derived from measurement inaccuracy or variations in properties that cannot be detected in the sample range [54]. The sill value is the upper limit of the fitted semivariogram model [47]. The ratio of nugget to sill was a criterion to classify the spatial dependence of soil properties and it reflects the influence of regional factors (nature) and the role of the non-regional factors (human factors). The range of the semivariogram (A0) represents the average distance through which the variable semivariance reaches its peak value. A small effective range implies a distribution pattern composed of small patches. The cross-validation value is the coefficient of determination (R2) of the correlation between the measured values and the cross-validation values, which were predicted based on the semivariogram and neighbor values [53]. Previous studies have shown that the semivariogram often differs considerably from its regional counterpart [33], [53]. Fig. 3 and Table 6 show the tested variables over the study area that were modeled using spherical, gaussian and exponential semivariograms with the lower nugget effect based on the coefficients of determination (R2) and residual sum of square (RSS). In this respect, a low ratio (less than 25% as was found for pH, Ca2+, Mg2+, and HCO3-) means that a large part of the variance is introduced spatially. This implies a strong spatial dependency of the variable, most likely due to intrinsic factors, including soil formation factors, such as soil parent materials, topography, and/or climate. A high ratio (more than 75% was found with Na+ and Cl-) often indicates a weak spatial dependency in the present sampling that was most likely due to extrinsic factors, including contamination, irrigation, and soil management practices, where the fertilization and the soil chemical properties were continuously shifting. In this research the Nug/Sill ratios of CO32-, and total salt content are 67.27 and 67.53% in the range between 25% and 75%, indicating the total salt content of the soil in Yanqi basin was influenced by the above factors. Additionally, all the other variables have a moderate spatial dependency for both intrinsic and extrinsic factors. The effective range calculated for variograms of different soil salt properties was 50–280 m, indicating that the sample distance was adequate for the characterization of the spatial variability of these properties.

thumbnail
Figure 3. The semivariance function diagram of the soil salt properties of Yanqi basin.

https://doi.org/10.1371/journal.pone.0106079.g003

thumbnail
Table 6. The spatial variation parameters of the soil salt properties of Yanqi basin.

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

The main application of geostatistics in soil science has been estimating and mapping the chemical properties in the soil of unsampled areas. Maps for each of the soil properties can be obtained using an ordinary kriging interpolation based on the best-fit semivariogram model. The skewness, defined as more than +1 or less than -1, indicates that some soil properties were not normally distributed. For these properties, it is difficult to estimate the semivariogram and doing so would result in a high value of kriging standard deviations [53]. Lognormal kriging with non-linear transformations is an alternative method for dealing with a data set with outliers or a non-normal distribution [31], [33]. Since the concentration of the variables with the non-normal distribution had a lognormal distribution, their concentrations were log-transformed, resulting in more regular variograms. The kriging interpolations were performed on the log-concentrations and the estimated values were back-transformed by an exponential function. As seen in the results, total amounts of salt, pH, EC, and the concentrations of Na+, K+, Mg2+, CO32-, and Cl- were relatively higher when from farmland and grassland. The spatial patterns of these variables had a significant geographical distribution with their primary occurrence in the western, north and central areas, which were higher than in other areas. This further proves the effect of the intrinsic factors of topography, soil forming factors and soil type, and extrinsic factors of soil management practices, such as fertilization, use of organic fertilizer and land management in farms, on the spatial distribution of soil chemical properties [33], [53]. The soil salinity tended to increase from the margin to the center across the study area, where agricultural wells are denser. It is clear that substantial soil salinization has taken place in these areas due to the effects of land management, farming and climate conditions, such as rainfall and high evaporation. This shows that more attention should be paid to these areas to prevent future problems. Soil pH in the study area ranged from 7.95 to 8.55 in most parts of the Yanqi basin. This pH range falls into the middle level of values that meets the fundamental conditions for plant growth and fertility. The ECs of the soil in all sections were within the acceptable value of <4 dS/m as based on the soil quality standard given by Bao et al. [45]. Areas presently not affected by salinity, but near saline areas, are potential areas for the development of salinity in the future, especially if they are also low-lying. In this respect, it is important to take necessary precautions and implement proper land use plans and cultivation practices.

Assessment of soil salinization risk in Yanqi basin.

Soil salinization occurs mainly in arid and semi-arid regions. It may arise due to climate, but is more likely to occur when irrigation practices alter the natural salt balance. Irrigation promotes soil salinization by raising the water table of the underlying aquifer, thus carrying salts upwards. Salts, unlike water, remain in the soil as evaporation and plant transpiration take place, thereby amplifying the salinity. The accumulation of salts in the surface and near-surface zones of soil is a major issue of environmental degeneration and is one of the main causes of low crop yields, loss of land and decreased production. The risk of soil salinization in the Yanqi basin is presented in Table 7. We observed that the whole area has a low salinization risk, as previously mentioned, mainly due to anthropogenic activities and climatic variation [39], [40]. The overall salinization risks for farmland, grassland, forest, urban construction areas and desert were none, none, none, low, and moderate, respectively. The mean values of the salinization risk indicators, except the SAR in farmland and desert, are lower than the corresponding maximum limitations that are predicted for the risk (Table 7). This indicates that farmland and desert in the study area have a low salinization risk, but they have mean values of 15.2 and 32.7 for SAR, which showed a moderate grade of salinization (Table 1, Table 7). This indicates a potential risk for environment declination. This research showed the grassland and forest have no risk of salinization, while the areas of urban construction have a moderate risk of salinization. In terms of spatial distribution, the soil salinization risk and the soil salt ions present at higher concentrations in the study area were similar. Environment declination due to soil salinity within the study area almost all took place on the edge of a lake, pond or river. These areas are heavily affected by climate warming, which in turn results in an increase in the amount of evapotranspiration that exceeds 2438.9 mm/a and in the average annual rainfall <100 mm (the drought and waterlog) (Fig. 4). Apart from natural factors, the main driving factors that jointly determined how local dwellers changed the landscape pattern were land use policies, economic systems and population growth. Human activity increases salinization through excessive application of irrigation water without adequate drainage. In addition, the cultivation of grassland is another major cause of inland salinization [40], [41]. On one side of the study area, plenty of grass landscape suffered damage, and, in some regions, the land salinization and desertification problem were serious enough to destroy the harmony between the material cycle and energy flow of the ecosystems [42]. On the other side of the study area, the habitat for wildlife was deteriorating, thus seriously threatening the biological diversity. Therefore, establishing and modifying policy, adjusting the irrigation system, improving drainage, dredging the surface water system, promoting circulation between surface water and underground water, and setting up wetland resource monitoring systems are necessary in order to restore damaged wetland and grassland [49], [50], [55].

thumbnail
Figure 4. The monthly evaporation, precipitation and temperature of Yanqi basin in 2011.

https://doi.org/10.1371/journal.pone.0106079.g004

thumbnail
Table 7. Classification of soil salinization risk taken at 0–20 cm depth within the investigated land use and land cover categories in Yanqi basin.

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

Conclusions

(1) From all the soil samples taken in this study, it can be gathered that there are numerous salt properties that vary largely between different samples. This analysis shows that the main salt ions in the soil were K+, Ca2+, Na+, Cl-, Mg2+, and SO42-, which accounted for 84.41% of the total salt content of the soil samples. Conversely, the concentrations of HCO3- and CO3- were very low. Except for the high variation found for the amounts of Ca2+ and K+ (191.67% and 226.25%, respectively), the other soil salt properties of Yanqi basin had moderate levels of variation (10%<CV<100%).

(2) From the analysis shown that the average values of Na+, Mg2+, SO42-, Ca2+, total amounts of salt and K+ being higher in farmland than grassland, desert, and urban construction areas. This work also determined that the maximum average values of Na+, Mg2+, K+, SO42-, total amount of salts, and Cl- were found in coarse crystalline rock weathered material, sandy shale of weathered material and lacustrine deposits. Within the five land use types examined, they were not obvious which classes of elements were enriched and the differences between these groups were small.

(3) PC analysis determined that PC1 (Cl-, Na+, SO42-, EC, and pH) and PC2 (Ca2+, K+, Mg2+, and total amount of salts) originated from artificial sources, while PC3 and PC4 (CO3- and HCO32-) originated from natural sources. Together, this research shows that Ca2+, K+, SO42 and the total amount of salts were influenced by both artificial and natural sources. Clustering analysis is consistent with the results from the PC analysis.

(4) From the geo-statistical point of view, it can be speculated that pH and soil salt ions, such as Ca2+, Mg2+ and HCO3-, had a strong spatial dependency. Meanwhile, Na+ and Cl- had only a weak spatial dependency, which was probably due to extrinsic factors, such as contamination, irrigation, and current soil management practices. We evaluated the EC, SAR and total salt content standard to reveal the risk of soil surface salinization. Soil salinization indicators suggest that the entire area had a low risk of salinization as mentioned previously, and this risk was mainly due to anthropogenic activities and climate variation. It is recommended that management of salinized land be preceded by an assessment of local factors and processes that may affect land composition.

Although the overall soil environment was healthy in the Yanqi basin, human activity, such as excessive groundwater pumping, have negatively impacted conditions by inducing soil salinization in the oasis. This matter deserves increased attention. This study can be considered an early warning of soil salinization and alkalization in the Yanqi basin. It can also provide a reference for environmental protection policies and for rational utilization of land resources in the arid region of Xinjiang, northwest China, as well as for other oases of arid regions in the world.

Acknowledgments

Many thanks to the members from the Institute of Geographical Science and Tourism of Xinjiang Normal University, Urumqi, China, and the members form College of Resources and Environment Sciences of Xinjiang University, Urumqi, China for the data collection and assistance in processing.

Author Contributions

Conceived and designed the experiments: ZZ JA. Performed the experiments: ZZ JA HY. Analyzed the data: ZZ JA. Contributed reagents/materials/analysis tools: ZZ JA HY. Contributed to the writing of the manuscript: ZZ.

References

  1. 1. Yang YJ, Yang SJ, Liu GM, Yang XY (2005) Space-Time Variability and Prognosis of Soil Salinization. Pedosphet 15(6): 797–804.
  2. 2. Amezketa E (2006) An integrated methodology for assessing soil salinization, a pre-condition for land desertification. Journal of Arid Environments 67(4): 594–606.
  3. 3. Masoud AA, Koike K (2006) Arid land salinization detected by remotely-sensed landcover changes: A case study in the Siwa region, NW Egypt. Journal of Arid Environments 66(1): 151–167.
  4. 4. Yang JS (2008) Development and prospect of the research on salt-affected soils in China. Acta Pedologica Sinica 45(5): 837–845
  5. 5. Stirzaker RJ, Cook FJ, Knight JH (1999) Where to plant trees on cropping land for control of dryland salinity: some approximate solutions. Agricultural Water Management 39(2): 115–133.
  6. 6. Ghassemi F, Jakeman AJ, Nix HA (1995) Salinisation of land and water resources: human causes, extent, management and case studies. CAB international.
  7. 7. Chai S X, YangBZ, Wang XY, Wei L, Wang P, et al. (2008) Analysis of salinization of saline soil in west coast area of Bohai gulf. Rock and Soil Mechanics 29(5): 1217–1221 (In Chinese).
  8. 8. Wang XL, Zhang FR, Wang YP, Feng T, Lian XJ, et al. (2013) Effect of irrigation and drainage engineering control on improvement of soil salinity in Tianjin. Transactions of the Chinese Society of Agricultural Engineering (20): 82–88 (In Chinese).
  9. 9. Ren JG, Heng XL, Xi JM, Li JL (2005) Salinization characteristics of the soil in Yeerqiang river valley, Xinjiang. Soils 37 (6): 635–639 (In Chinese).
  10. 10. Chen XB, Yang JS, Liu CQ, Hu SJ (2007) Soil Salinization Under Integrated Agriculture and Its Countermeasures in Xinjiang. Soils 39(3): 347–353 (In Chinese).
  11. 11. Wang YG, Xiao DN, Li Y (2008) Spatial and Temporal Dynamics of Oasis Soil Salinization in Upper and Middle Reaches of Sangonghe River, Northwest China. Journal of Desert Research 28(3): 478–484 (In Chinese).
  12. 12. Jiang L, Li PC, Hu AY, Yi X (2009) Analysis and evaluation of soil salinization in oasis of arid region. Arid Land Geography 32(2): 234–239 (In Chinese).
  13. 13. Sawut M, Eziz M, Tiyip T (2013) The effects of land-use change on ecosystem service value of desert oasis: a case study in Ugan-Kuqa River Delta Oasis, China. Canadian Journal of Soil Science 93(1): 99–108.
  14. 14. Ling H, Xu H, Fu J, Fan Z, Xu X (2013) Suitable oasis scale in a typical continental river basin in an arid region of China: A case study of the Manas River Basin. Quaternary International 286: 116–125.
  15. 15. Eziz M, Yimit H, Mamat Z, Li JT (2013) Driving forces of farmland dynamics and its ecological effects in Keriya Oasis in recent 60 years. Agricultural Research in the Arid Areas 3: 033 (In Chinese).
  16. 16. Triki I, Trabelsi N, Zairi M, Dhia HB (2014) Multivariate statistical and geostatistical techniques for assessing groundwater salinization in Sfax, a coastal region of eastern Tunisia. Desalination and Water Treatment 52(10-12): 1980–1989.
  17. 17. Fu S, Wei CY (2013) Multivariate and spatial analysis of heavy metal sources and variations in a large old antimony mine, China. Journal of Soils and Sediments 13(1): 106–116.
  18. 18. Gil PM, Saavedra J, Schaffer B, Navarro R, Fuentealba C, et al. (2014) Quantifying effects of irrigation and soil water content on electrical potentials in grapevines (Vitis vinifera) using multivariate statistical methods.Scientia. Horticulturae173: 71–78.
  19. 19. Hu W, Shao MA, Wan L, Si BC (2014) Spatial variability of soil electrical conductivity in a small watershed on the Loess Plateau of China. Geoderma 230: 212–220.
  20. 20. Oliver MA, Webster R (2014) A tutorial guide to geostatistics: Computing and modelling variograms and kriging.Catena. 113: 56–69.
  21. 21. Emadi M, Baghernejad M (2014) Comparison of spatial interpolation techniques for mapping soil pH and salinity in agricultural coastal areas, northern Iran. Archives of Agronomy and Soil Science 60(9): 1315–1327.
  22. 22. Elbasiouny H, Abowaly M, Abu-Alkheir A, Gad A (2014) Spatial variation of soil carbon and nitrogen pools by using ordinary Kriging method in an area of north Nile Delta, Egypt. Catena 113: 70–78.
  23. 23. Wu Y, Wang Y, Xie X (2014) Spatial occurrence and geochemistry of soil salinity in Datong basin, northern China. Journal of Soils and Sediments 1–11.
  24. 24. Li SJ, Sun YN, Wang HB, Chen ZW (2013) Impact of Soil Nutrient Contents and Spatial Variability in Different Sampling Schemes under the Conservation Tillage.Advanced Materials Research. 718: 316–320.
  25. 25. Bilgili AV (2013) Spatial assessment of soil salinity in the Harran Plain using multiple kriging techniques. Environmental monitoring and assessment 185(1): 777–795.
  26. 26. Sylla M, Stein A, Van Breemen N, Fresco LO (1995) Spatial variability of soil salinity at different scales in the mangrove rice agro-ecosystem in West Africa. Agriculture, ecosystems & environment 54(1): , 1–15.
  27. 27. Ammari TG, Tahhan R, Abubaker S, Al-Zu'Bi Y, Tahboub A, et al. (2013) Soil salinity changes in the Jordan Valley potentially threaten sustainable irrigated agriculture.Pedosphere. 23(3): 376–384.
  28. 28. Walter C, McBratney AB, Douaoui A, Minasny B (2001) Spatial prediction of topsoil salinity in the Chelif Valley, Algeria, using local ordinary kriging with local variograms versus whole-area variogram. Soil Research 39(2): , 259–272.
  29. 29. Jordán MM, Navarro-Pedreno J, García-Sánchez E, Mateu J, Juan P (2004) Spatial dynamics of soil salinity under arid and semi-arid conditions: geological and environmental implications. Environmental Geology 45(4): 448–456.
  30. 30. Bai YL, Li BG, Hu KL (1999) Spatial variability of soil salt and its composing ions in salt-affected soil in Huang-Huai-Hai plain. Soil and fertilizer (3): 22–26 (In Chinese).
  31. 31. Hu Kl, Li BG, Chen DL (2001) Spatial variability of soil water and salt in field and estimating soil salt using CoKriging. Advances in water science 12(4): 460–466 (In Chinese).
  32. 32. Xu Y, Chen YX, Shi HB (2004) Scale effect of spatial variability of soil water- salt. Transactions of the Chinese Society of Agricultural Engineering 20(2): 1–5 (In Chinese).
  33. 33. Lin J, Anwar M, Dilbar S (2007) Investigation of the spatial variability of soil salts in saline soil in Xinjiang. Research of Soil and Water Conservation 14(6): 189–192.
  34. 34. Aragüés R, Medina ET, Clavería I, Martínez-Cob A, Faci J (2014) Regulated deficit irrigation, soil salinization and soil sodification in a table grape vineyard drip-irrigated with moderately saline waters. Agricultural Water Management 134: , 84–93.
  35. 35. Bouksila F, Bahri A, Berndtsson R, Persson M, Rozema J, et al.. (2013) Assessment of soil salinization risks under irrigation with brackish water in semiarid Tunisia. Environmental and Experimental Botany 92: , 176–185.
  36. 36. Liang SY, Wu TN, Wu YS, Chou YC, Lee CH (2013) Assessment of Aquifer Salinization Beneath an Offshore Industrial Park Based on Solute Transport Calculation. Advanced Materials Research 779: , 1285–1288.
  37. 37. Chen L, Feng Q (2013) Geostatistical analysis of temporal and spatial variations in groundwater levels and quality in the Minqin oasis, Northwest China.Environmental Earth Sciences. 70(3): 1367–1378.
  38. 38. Wang LC, Wu RW, Gao J (2014) Spatial coupling relationship between settlement and land and water resources–based on irrigation scale–A case study of Zhangye Oasis.Advanced Engineering and Technology 225.
  39. 39. Mamat Z, Yimit H, Eziz M, Ablimit A (2013) Analysis of the Ecology-Economy Coordination Degree in Yanqi Basin, Xinjiang, China. Asian Journal of Chemistry 25(16): 9034–9040.
  40. 40. Mamat Z, Yimit H, Eziz A (2014) Oasis land-use change and its effects on the eco-environment in Yanqi Basin, Xinjiang, China. Environment Monitoring and Assessment 186(1): 335–348.
  41. 41. Wu JS, Zhang YQ, Liu ZH, Peng J, He JF (2010) Land salnization monitorng with remote sensing on Yanqi County, Xin jiang. Arid Land Geography 32(7): 251–257 (In Chinese).
  42. 42. IUSS Working Group, WRB (2006). World reference base for soil resources. World Soil Resources Report,103..
  43. 43. Hargreaves G H, Allen R G (2003) History and evaluation of hargreaves evapotranspiration equation. Journal of Irrigation and Drainage Engineering 129(1): 53–63..
  44. 44. Lu RK (2000) Soil and agricultural chemistry analysis. Beijing/China Agricultural.
  45. 45. Bao SD (2005) Soil Agricultural Chemistry Analysis. Beijing/China Agriculture Press. pp: 17–200.
  46. 46. Carter M R (1993) Soil sampling and methods of analysis. CRC Press.
  47. 47. Lark RM (2001) Geostatistics for environmental scientists.European Journal of Soil Science 52(3): , 526–526.
  48. 48. Cressie N (1988) Spatial prediction and ordinary kriging. Mathematical Geology 20(4): 405–421.
  49. 49. Fan LQ, Yang JG, Xu X, Sun ZJ (2012) Salinity characteristics and correlation analysis of saline soil in irrigation area of Ningxia. Soil and Fertilizer Sciences in China 6: 003 (In Chinese).
  50. 50. Metternicht G, Zinck JA (1997) Spatial discrimination of salt-and sodium-affected soil surfaces. International Journal Remote Sensing 18(12): 2571–2586.
  51. 51. Agriculture Department of Xinjiang (1996) Soil Survey Office in Xinjiang. Xinjiang Soil. Beijing/Science press 458–464.
  52. 52. Wilding LP (1984) Spatial variability: Its documentation, accommodation and implication to soil surveys//MNielson D R, Bouma J. Soil Spatial Variability. Purdoc, Wageningen: 166–193.
  53. 53. Fan XM, Liu GH, Liu HG (2014) Evaluating the spatial distribution of soil salinity in the Yellow river delta based on Kriging and Cokriging Methods. Resources Science 36(2): 0321–0327 (In Chinese).
  54. 54. Wang SX, Dong XG, Liu YF (2009) Spatio-Temporal variation of subsurface hydrology and groundwater and salt evolution of the oasis area of Yanqi Basin in 50 Years recently. Geological Science and Technology Information 28(5): 101–108 (In Chinese).
  55. 55. Li XG, Lai N, Chen SJ, Mamattursun E (2014) Spatial variability of soil salt based on geostatistics and GIS in the oasis of the lower reaches of Kaidu River: a case study on Yanqi County. Gco-graphy and Gco-information Science 30(1): 105–109 (In Chinese).