Next Article in Journal
Climate and Land Use Changes on Streamflow and Subsurface Recharge in the Fluvià Basin, Spain
Next Article in Special Issue
Optimal Allocation of Water Resources Based on Water Supply Security
Previous Article in Journal
European Water Footprint Scenarios for 2050
Previous Article in Special Issue
Variations in the Runoff-Sediment Relationship of the Weihe River Basin Based on the Copula Function
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Variations on Soil Salinity and Sodicity and Its Driving Factors Analysis under Microtopography in Different Hydrological Conditions

1
Northeast Institute of Geography and Agroecology, Chinese Academy of Sciences, 4888 Shengbei Street, Changchun 130102, China
2
Department of Agricultural and Resource Economics, The University of Tennessee, Knoxville, TN 37996, USA
3
Dynamic Solutions, LLC, Knoxville, TN 37919, USA
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Water 2016, 8(6), 227; https://doi.org/10.3390/w8060227
Submission received: 7 April 2016 / Revised: 15 May 2016 / Accepted: 23 May 2016 / Published: 27 May 2016
(This article belongs to the Special Issue Tackling Complex Water Problems in China under Changing Environment)

Abstract

:
Over three million hectares of salt-affected soils characterized with high salinity and sodicity caused serious land degradation in Songnen Plain, northeast China. Soil salinity–sodicity heterogeneous distribution under microtopography is usually influenced by several environmental factors. The side direction movement of soil water driven by water from depression is the key factor that aggravates the soil salinization under microtopography in dry condition. In this study, the differences in surface soil salinity–sodicity (0–10 cm) between dry year and wet year were compared, and the relationship between soil salinity–sodicity and environment factors such as ground elevation, surface ponding time, surface ponding depth, and soil moisture at four soil layers (0–10, 10–30, 30–60, and 60–100 cm) were analyzed using redundancy analysis (RDA) and simple correlation analysis (Pearson analysis) for two different hydrological years. Analyzed soil salinity–sodicity parameters include soluble ions (Na+, K+, Ca2+, Mg2+, CO32−, HCO3, Cl and SO42−), salt content (SC), electrical conductivity (EC), sodium adsorption ratio (SAR), and pH. Results showed that values of SAR, Cl, and SO42− were significantly higher in dry year than in wet year, while Ca2+, Mg2+, K+, and HCO3 showed the opposite results. Values of Na+, CO32−, and EC were significantly higher at higher ground elevation gradient (20–40 cm) in dry year than wet year. Redundancy analysis indicated that spatial distributions and variations of salinity and sodicity in surface soil layer were related with environmental factors of ponding depth, ponding time and ground elevation in wet year, and they were related with ground elevation, ponding depth, ponding time, and soil moisture at 30–60 and 60–100 cm soil layer in dry year. Ponding depth and ground elevation rank first and second as the influential factors of the spatial distribution and variation of soil salinity–sodicity in wet year. However in dry year, primary and secondary influential factors are ground elevation and soil moisture at 60–100 cm soil layer.

1. Introduction

Salt-affected area covers about 1 × 109 ha in the world [1]. Soil salinization is an increasing environmental problem because of adverse effects of soil salinity and sodicity on soil properties and behavior [2,3]. Songnen Plain, northeast China, is one of three major saline-sodic soil distribution regions in the world. The area of saline-sodic soil is about 2.57 × 106 ha, approximately 22.62% of its total area [4,5]. The area of saline-sodic soil is increasing at an annual rate of 1.35 × 104 ha. In some counties, the salinization land accounts for more than 20% of the total area and has been changed into wasteland [6]. Soil salinization has caused serious soil degradation, which negatively influences the economic development and decreased the yield of local crops.
The mechanism of soil salinization and its development have been a great concern for researchers and growers in the arid or semiarid regions [5]. Saline-affected soil formation was correlated with critical depth of shallow groundwater (GW) level, except when controlled by arid climate [7]. Salinization occurred when GW depth was less than 2.0 m with salinity higher than 0.48 dS/m [8,9]. For the GW with lower salinity, the GW table should not be shallower than 1.5 m [10]. By research of the effect of shallow underground water on the formation process of salinity–sodicity soil in Songnen Plain, results indicate that the shallow ground water level has already exceeded the critical water level (2–3 m) [11,12]. However, when the groundwater table reached 5 m, the soil salinization still continued and aggravated [13,14], which is unexplainable from the viewpoint that soil salinization is solely caused by the rise of underground capillary water. One hypothesis about the evolution of saline-sodic soil in Songnen Plain is proposed that the saline sodic soil development is affected by human activities and climatic conditions. The influence of underground water on soil salinization may be weakened when soil salinization of Songnen Plain has developed to a certain degree. However, the existence of microtopography may cause the uneven movement of water and soil salt, which affects the spatial distribution of soil salinity–sodicity. In dry climate conditions, ponding water with soluble salts moves laterally from high soil moisture to low moisture soil along the moisture gradient because of evaporation. Hence, salt was left in surface soil layer at the middle and upper part of microtopography, which will aggravate the soil salinization in this area.
The scholars have noticed the heterogeneous distribution of soil salinity–sodicity under microtopography and carried out a lot of studies, such as soil moisture transform and distribution characteristics of soil salt under microtopography [10,15] and analysis of spatial variation of soil salinity–sodicity [16,17,18]. It is shown that the spatial variation of soil salinity is mainly affected by the environmental factors (the topography, hydrological process and climate conditions). However, very little research has been done on the spatial variation of soil salinity–sodicity and their relation with environment factors in Songnen Plain [12,19]. In addition, regular statistical analysis and regression analysis are mainly employed. In order to examine the relationships between surface soil salinity–sodicity and environmental factors, ordination techniques are preferred approaches. In the numerous studies, two multivariate ordination techniques were primarily used in the soil–environment research [20,21,22]. One of the approaches consisted of e linear ordination methods such as principal components analysis (PCA) and redundancy analysis (RDA). These analyses are performed to determine the variation in the objective data within a narrow range, which could be attributed to linear coupling. The other approach used unimodal based methods of correspondence analysis (CA), detrended correspondence analysis (DCA), and canonical correspondence analysis (CCA) to determine the variation in the objective data over a wide range. DCA should be carried out on the soil salinity–sodicity data to determine the most appropriate method, with the gradient length as the criterion [23].
The results of DCA showed that linear methods (PCA and RDA) were preferable for the dataset considered in this study. Redundancy analysis was chosen because of the relation analysis between the two sets of variables (soil salinity–sodicity and environment factor), while principal component analysis (PCA) involved a single set of soil data. In order to gain a better understanding of the relationship between soil salinity–sodicity and environmental factors. RDA was used to analyze the relationships among soil salinity–sodicity, microtopography and hydrological elements. The objectives of this research were to: (a) compare the differences of surface soil salinity–sodicity between dry year and wet year; (b) analyze the correlations between environmental attributes and surface soil salinity–sodicity changes in different hydrological years; and (c) identify the key driving factors of spatial variations in surface soil salinity–sodicity in different hydrological years.

2. Material and Methods

2.1. Study Area

The experiment is carried out in Da’an Sodic Land Ecological Experiment Station of Northeast Institute of Geography and Agroecology, Chinese Academy of Sciences (Figure 1). This region is located in 123°50′27″–123°51′31″ of east longitude and 45°35′58″–45°36′28″ of north latitude, and features temperate continental monsoon climate with an average annual precipitation of 413.7 mm and an average annual evaporation of 1696.9 mm. The experiment station is surrounded by low flood plain for Nenjiang ancient channel. The land surface is featureless and gently undulating, with the alternate distribution of small mounds and depressions. The complex coexistence of different level of soil salinity–sodicity is the common characteristics in Songnen Plain, northeast China [24].

2.2. Field and Laboratory Analysis

A representative saline-sodic area measuring 100 m × 100 m was selected with the saline-sodic soil complex of different soil salinity–sodicity levels. A ridge was constructed along the boundary of the area in case surface runoff entered the area after rainfall events. The ridge was 36 cm high and 40 cm wide. The experiment area was divided into 10 m × 10 m grid squares, and 40 grids were randomly selected as sampling points (Figure 2). The land surface in the area is gently undulating with microtopography of coexistence of mounds and depressions. The maximum difference in elevation is about 36 cm. The grids with the highest elevation (such as sampling points of No. 2, No. 3, No. 4 and No. 21) are located in the middle bottom (Figure 2). The grids having the lowest elevation (such as sampling points of No. 33, No. 34, No. 36, No. 37, No. 39 and No. 40) are located in the top right corner (Figure 2). Soil of experiment area was naturally saline-sodic soils without tillage and grazing disturbances. The plant communities in this area mainly include Phragmitesaustralis, Chlorisvirgata, Suaedaglauca and Puccinelliatenuiflora. Typically, Suaedaglauca is distributed in the highest elevation area (mounds), and Phragmitesaustralis is found in the lowest elevation area (depressions). The soil type is loam soil by USDA textural classification triangle. Dry bulk density ranged 1.49 to 1.67 g/cm3. Shallow groundwater level is lower than 5 m. The shallow groundwater has the high values of TDS (0.2%–0.4%), and pH ranges from 7.2 to 7.9. The main ions of the shallow groundwater are Na+, Mg2+, HCO3 and Cl.
Experimental data were obtained in 2005 (wet year) and 2007 (dry year). The experiment lasts for 5 months for each year, from 1 May to 1 October. The surface ponding time and surface ponding depth were observed every five days for no-rainfall period, and every day after the rainfall event. At the end of May, July and September, the soil moistures were measured, and the soil moistures is the mean value of three observations. Soil moistures were observed at four depths of soil layers: 0–10, 10–30, 30–60, and 60–100 cm. Ground elevation was measured by DSN232 balance level.
At the end of September, soil samples were collected at soil depth of 0–10 cm at these 40 sampling points for laboratory analysis of pH, EC, Na+, K+, Ca2+, Mg2+, CO32−, HCO3, Cl, and SO42−. All soil samples were air-dried and then passed a 1-mm round-hole sieve for chemical analyses. Soluble salt estimates were based on soil: water 1:5 extracts. The pH and EC of the extracts were determined using a pH meter and a conductivity meter, respectively. The concentration of Na+, K+, Ca2+, and Mg2+ were determined by atomic absorption spectrometry (GBC-906AAS). Anion concentrations (CO32−, HCO3, Cl, and SO42−) were determined by standard methods [25]. The concentration of CO32− and HCO3 were determined by neutral titration method. The concentration of Cl was determined by silver nitrate titration method, and the concentration of SO42− was determined by barium sulfate turbidimetric method. Sodium adsorption ratio (SAR) was calculated by the following equation where the concentrations of soluble cations are expressed in mmolc L−1 [26].
S A R = [ N a + ] ( [ C a 2 + ] + [ M g 2 + ] ) / 2

2.3. Redundancy Analysis

RDA computation was performed using CANOCO for Windows 4.5 [27]. In the RDA computations, the correlation matrix option was selected, and scaling was conducted on a correlation bi-plot. In each analysis, step-wise forward selection was used to reduce the environmental variables to those most correlated with the axes. The significance of each environmental variable was tested using 499 unrestricted Monte Carlo permutations. For each analysis, the variance explanation of each environmental variable (marginal effect), and their additional variance explained when included in the model (conditional effect) [27] were also provided. The redundancy analysis (RDA) method is capable of presenting the research objects ordination and environment factors ordination in one diagram, and the relationships can be directly observed between them. The research objects are represented by the dotted arrow, the solid arrow represents the environment factors, and the quadrant that the arrow lies in indicates positive and negative correlations between environment factors (research objects) and ordination axes. The length of solid arrow connecting line means correlation between an environment factor and research object distribution, i.e., the longer line indicates that the environment factor has more influences on the distribution of research object. The angle between two arrows means the correlation between two environment factors (research objects), i.e., the smaller angle determines greater relevance. The angle between the arrow and ordination axis means the correlation between an environment factor (research objects) and ordination axis [23]. In this study, research objectives are soil ions (Na+, Ca2+, Mg2+, K+, CO32−, HCO3, SO42− and Cl), soil salt content (SC), electrical conductivity (EC), sodium adsorption ratio (SAR), and pH. Considered environmental factors include surface ponding time (d), surface ponding depth (h), soil moistures in different soil layers (W0–10, W10–30, W30–60, and W60–100), and ground elevation (H).

3. Results and Analysis

3.1. Rainfall during Different Hydrological Years

In this experiment, two years data were selected to analyze correlation between surface soil salinity–sodicity and environment factors under dry and wet conditions. The 50-year (1951–2009) average rainfall during May to October in the study area is approximately 350.7 mm. Year 2005, with rainfall of 455.9 mm during May to October, was selected to represent the wet condition. Year 2007 was selected to represent the dry condition with the rainfall of 224.5 mm during May to October (Figure 3).

3.2. Comparison of Surface Soil Salinity–Sodicity Parameters between Wet and Dry Years

The descriptive statistics of soil ion contents and associated properties at the soil layer (0–10 cm) at these 40 measurement points are given in Table 1. The content of Na+ was significantly higher than other cations with a mean of 809 mg/kg in wet year (2005) and 1238 mg/kg in dry year (2007) (p < 0.05), and the content of HCO3 was significantly higher than other anions with a mean of 3369 mg/kg (2005) and 1583 mg/kg (2007), respectively (p < 0.05) (Table 1). The mean pH value in 2005 and 2007 were 9.51, and 9.48, respectively. Soil alkalinity is associated with the presence of sodium carbonates (Na2CO3) or sodium bicarbonate (NaHCO3) in the soil, so the extracts were alkaline or strongly alkaline.
Variation coefficient (CV) was the most important factor in describing the distribution and variability of soil property. Spatial variation can be divided into three categories: high spatial variation (CV > 1), moderate spatial variation (0.1 < CV < 1), and low spatial variation (CV < 0.1) [28]. Experimental results indicated soil soluble salts and their chemical properties displayed high and moderate spatial variations except the pH. The calculated coefficients of variation for soil ion ranged from 34% to 203%. CV value for EC, SAR, and SC are from 74% to 197%. CV values of Na+, CO32−, SO42−, SAR and EC were relatively higher under dry weather condition (2007) than wet weather condition (2005) (Table 1).
Parameter values of SAR, Cl, and SO42− were significantly higher in dry year (2007) than in wet year (2005) (p < 0.05), while Ca2+, Mg2+, K+, and HCO3 showed the opposite results (Table 1). There were no significant differences in parameter values in EC, SC, Na+, and CO32− between these two years. Because of high CV values for these parameters (CV > 1), soil salinity–sodicity in different ground elevation may be different under microtopography [10,15]. To further evaluate the impact of ground elevation gradient on these soil salinity–sodicity parameters (EC, SC, Na+, and CO32−) between dry and wet weather conditions, ground elevation was divided into three gradients: 0–10, 10–20, and 20–40 cm. Results of comparison of these parameter values between dry and wet years at each ground elevation gradient are given in Table 2. Ground elevation had significant impact on these soil salinity–sodicity parameters (Na+, CO32−, EC, and SC) (Table 2). Parameter values of Na+, CO32−, SC, and EC showed significant differences between dry and wet years only at higher ground elevation (20–40 cm): values of Na+, CO32−, SC, and EC were significantly higher in dry year (2007) than wet year (2005) (Table 2). While parameter values of SC were significantly lower in dry year (2007) than the wet year (2005) at middle and lower ground elevation gradients (10–20 cm and 0–10 cm). Thus, high evaporation and low rainfall can increase soil salinization in soil surface layer (0–10 cm), especially at higher ground elevation gradient in dry climate condition.

3.3. Relationship between Surface Soil Salinity–Sodicity and Environment Factors

The RDA was used to analyze the relationship between surface soil salinity–sodicity and environmental factors based on the soil salinity–sodicity data matrix and an environmental data matrix. The connection between soil salinity–sodicity of the sample plots and the measured environmental factors in the study area as revealed by RDA are shown in Table 3. The results showed that 51.6% of soil salinity–sodicity spatial distribution was explained by these seven environment variables (sum of all eigenvalues is 1.0; sum of all canonical eigenvalues is 0.516) during the wet year, whereas 68.8% was explained during the dry year (sum of all eigenvalues is 1.0; sum of all canonical eigenvalues is 0.688). In addition, the cumulative variance explained for the first two ordination axes accounts for 100% of total cumulative variance. The research objects, environmental correlations, are higher for the first canonical axes, explaining 94.6% and 99.8% of the total cumulative variance in the wet and dry year, respectively (Table 3). Hence, these seven environment factors selected may, to a large extent, explain the spatial distribution of soil salinity–sodicity, by the first two ordination axes during different hydrological years. Through the above analysis, it is evident that the two ordination axes can describe the spatial distribution of the soil salinity–sodicity. These seven environment factors are able to explain the distribution characteristics of the surface soil salinity–sodicity by the first ordination axes, especially under dry weather conditions.
Monte Carlo permutations (499 iterations) were performed to test for significant relationship between soil salinity–sodicity data and the canonical axes. Results indicated that all canonical axes were significant (p = 0.002). From the intra-set correlations of the environmental factors with the axes of the RDA shown in Table 3, it can be seen that the selected environmental factors were mainly correlated with the first axis. During the wet year, there was significant relationship between environmental factors (ponding depth, ponding time and ground elevation) and the first axis (r = −0.5281, −0.4864 and 0.4339, respectively) (p < 0.05) while soil moistures at all four soil layers showed no significant relationship with the first axis (Table 3). Under wet year, ground elevation showed significant positive relationship with both the first and second axes (Table 3). The results showed surface ponding and microtopography were key factors, which affect distribution and variation of surface soil salinity–sodicity in wet year. During dry year, all environmental parameters showed significant relationship with the first axis except soils moistures at the upper layer (W0–10 and W10–30) (Table 3), which has also been revealed by the RDA ordination biplots (Figure 3 and Figure 4). Therefore, microtopography, surface ponding and soil moisture in deep soil layer are the significant factors affecting the distribution of salinity and sodicity in the soil during the dry year.
The relationship between parameters of soil salinity–sodicity and environment variations were analyzed using RDA and Pearson correlation analysis method. During the wet year, soil salinity–sodicity parameters (Na+, CO32−, Cl, EC, SAR, pH, and SC) were all positively correlated with ground elevation (Table 4 and Figure 4). Soil salinity–sodicity parameters (Na+, Ca2+, CO32−, Cl, HCO3, EC, pH, and SC) were negatively correlated with ponding time and ponding depth (Table 4 and Figure 4). There is an increase in soil ions (Na+, CO32−, and Cl) at the upper part of microtopography (mound). For Ca2+, Mg2+ and HCO3, they are richer at the middle part of microtopography than upper and lower part. SAR and EC values increase with the increase of ground elevation. SAR and EC were the highest at upper part of microtopography (Table 2 and Figure 3). However, during the dry year, surface soil salinity–sodicity parameters (Na+, CO32−, Cl, HCO3, SO42−, EC, SAR, pH, and SC) were all positively correlated with ground elevation (Table 4 and Figure 4). Surface soil salinity–sodicity parameters (Na+, CO32−, Cl, HCO3, SO42−, pH, SC, and SAR) were negatively correlated with ponding time and ponding depth, and were negatively correlated with soil moistures (soil layer: 30–60 cm and 60–100 cm) (p < 0.05) (Table 4 and Figure 5). The results indicated that spatial distribution and variations of soil salinity–sodicity in surface soil layer were related with ponding and microtopography in wet year, and they were related with ponding, microtopography, and moisture at deep soil layer in dry year. High values of surface soil salinity–sodicity are found at the upper part of microtopography (mound). Moderate values of salinity and sodicity are found at the middle part of microtopography. Mild value of salinity and sodicity are observed at the lower part of microtopography (depression). Surface soil salinity–sodicity was severe in dry year than wet year, especially at the upper part of microtopography.

3.4. Key Factors Affecting Spatial Distribution of Surface Soil Salinity–Sodicity

It is often of great interest to rank environmental variables in their importance for determining the surface soil salinity and sodicity distribution and variation. In the RDA analysis, Monte Carlo was used to determine the significance of each environment factor. The marginal effects of all environmental variables (i.e., the independent effect of each environmental variable) and conditional effects (i.e., the effect that each variable brings in addition to all the variables already selected) are shown in Table 4 and Table 5. During the wet year, ponding depth (h) is the most important environmental factor, which explains 56% of the spatial variation in the total cumulative variance, and ground elevation (H) is the second key environmental factor, which explains 16.7% of the spatial variation in the total cumulative variance. Results indicated that ponding depth and ground elevation, are the most important environmental factors affecting the spatial distribution and variation of soil salinity–sodicity in the wet year. Under dry year (2007), ground elevation (H) and soil water content at 60–100 cm soil layer (W60–100) are the key environmental factors (Table 6), which explain 83.8% of the spatial variation in the total cumulative variance together. During dry year, ground elevation is the most important environmental factor, which explains 67.4% of the spatial variation in the total cumulative variance, and W60–100 is the second key environmental factor, which explains 16.4% of the spatial variation in the total cumulative variance. Hence, in wet year, primary and secondary influence factors on the spatial distribution of surface soil salinity–sodicity are ponding and microtopography, separately. But in dry year, primary and secondary influence factors are microtopography and soil moisture in deep soil layer, separately. Soil moisture in deep soil layer is mainly supplied by shallow groundwater. Therefore, distribution of soil salinity and sodicity is affected by both microtopography and shallow groundwater, but the effects of the shallow groundwater are relatively lower than that of microtopography. Results from this research are consistent with other studies. Yang et al. [29] found influence of capillary rise from shallow groundwater was less important due to aggravation of soil salinization. Luo et al. [10] reported that the ground elevation has more effects on soil salinization than that of underground water.
Microtopography affects hydrological processes in local scale, and affects the distribution and variations of surface soil salinity and sodicity, which in turn affects the movement of soil water and salts. Low level of soil salinity–sodicity was found to be located in the depressions. Similar results have been documented by Kilic et al. [30]. Level of soil salinity–sodicity was observed to increase with the increase of ground elevation. Soil salinity–sodicity is aggravated due to the increase of Na+ and SAR. Soil tended to disperse due to the existence of excessive sodium. High sodicity degraded the soil physical properties, and caused soil compaction and slow internal drainage and the reduction of soil water storage and poor aeration [1,2,31]. The permeability of the upper soil is extremely low with saturated hydraulic conductivity less than 1 mm/day at the ground elevation of 10–40 cm [32]. Effective soil depth affected by rainfall and evapotranspiration was found not to exceed at 40 cm soil layer, and the shallow groundwater had slight effect on soil water contents above the depth of 80 cm soil layer [29]. Hence, the vertical movement of water (capillary rise and infiltration) is comparatively slow due to the aggravation of soil salts. Because of mild salinity and sodicity located in the depressions, soil physical properties were improved compared to the mound. In addition, owing to the difference in plant distribution, the infiltration capacity is higher in the depressions than in the mounds. Depressions were covered by Phragmitesaustralis due to seasonal waterlogging. Roots are intricately crossed and extended to 60 cm in the soil profile [33]. The advanced root system of Phragmitesaustralis can create more active soil pores, which in turn will increase the infiltration capacity [34].
Salts, derived from ponding and soil weathering reactions, accumulate in the surface of zones in the soil profile. This salinity fluctuation also changes with reason and rainfall [35]. In wet year, runoff with salt from mound flows into the depression due to of the excessive rainfall in rainy season. Because of higher infiltration capacity in the depressions, which can infiltrate a lot of salts into the deep soil or shallow groundwater to form desalination water layer. Thus, depressions were salt water exchange places. In a dry year, because the upper part and middle part of microtopography has higher temperature than that at lower part [10], the evaporation at the upper part and middle part of microtopography is higher than that at lower part. In the dry year, the evaporation shows relatively higher impact than the wet year. Based on the principles that the soil water flows from high to low gradient, soil moisture at the upper part and middle part of microtopography is supplied with water from the depression by side direction moving of moisture from low elevation to high elevation. During the process of side movement of moisture from low elevation to high elevation, the salts were also transported and accumulated through evaporation. The higher the evaporation, the faster the moisture is supplied, and the more the salt is accumulated. Low rainfall and high evaporation can increase the salt contents at the surface soil, especially at the upper part and middle part of microtopography [15]. The microtopography could influence the pattern and magnitude of spatial variability in salinity and alkalinity [36]. Therefore, under arid climate, the side direction movement of soil water driven by water with soluble salts from depression is the key factor, which aggravates the soil salinization under microtopography in dry condition.
Based on research of mechanism of soil salinization formation and its development, it is found that the side direction movement of soil water driven by ponding water under microtopography in dry conditions is the key factor which aggravates the soil salinization, which will lead to co-distribute high saline and sodic areas and non or low saline and sodic areas at different part of microtopography. Most of the traditional saline-sodic soil reclamation approaches are only effective for saline-sodic soils with good uniformity in ground elevations. Therefore, these approaches may not work in areas with spatial variations in ground elevations, where there are also high spatial variations in soil salinity–sodicity.
Continuous cropping was suggested together with good surface water drainage and improvement of soil structure prior to cropping was found to be beneficial [1]. The best way to ameliorate saline-sodic soil for the Patch with salinization involves increase in Ca2+ on the cation exchange site at the expense of Na+. The replaced Na+ together with excess soluble salts is removed from the root zone through infiltrating water as a result of excessive irrigations to achieve soil improvement.

4. Conclusions

The hydrological elements and microtopography are important factors affecting the spatial distribution of surface soil salinity and sodicity in Songnen Plain, northeast China. The values of SAR, Cl, and SO42− were significantly higher in dry year than in wet year, while Ca2+, Mg2+, K+, and HCO3 showed the opposite results. Values of Na+, CO32−, and EC were significantly higher at higher ground elevation gradient (20–40 cm) in dry year than wet year. The spatial distribution and variations of soil salinity–sodicity in surface soil layer were related with ponding and microtopography in wet year, and they were related with ponding, microtopography, and moisture at deep soil layer in dry year. The results of this research can provide some useful information for explaining mechanism of salinization process and utilization of saline-sodic soils in Songnen Plain, northeast China.

Acknowledgments

This research is supported by the National Basic Research Program of China (2015CB150802), the National Natural Science Foundation of China (41571210, 41371260,) and the National Science and Technology Basic Work of China (2015FY110500). We also wish to acknowledge the support of Da’an Sodic Land Experimental Station.

Author Contributions

The research presented here was carried out in collaboration between all authors. Fan Yang and Zhichun Wang had the original idea for the study. Fenghua An, Fan Yang and Hongyuan Ma designed/conducted the research methods. Xia Zhou and Zhijun Liu contributed to the writing of the paper. All authors discussed the structure and commented on the manuscript at all stages.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Qadir, M.; Schubert, S.; Ghafoor, A.; Murtaza, G. Amelioration strategies for sodic soils: A review. Land Degrad. Dev. 2001, 12, 357–386. [Google Scholar] [CrossRef]
  2. Shaw, R.; Brebber, L.; Ahern, C.; Wein, M. A review of sodicity and sodic soil behavior in Queensland. Aust. J. Soil Res. 1994, 32, 143–172. [Google Scholar] [CrossRef]
  3. Ammari, T.G.; Tahhan, R.; Abubaker, S.; Al-Zu’Bi, Y.; Tahboub, A.; Ta’any, R.; Abu-Romman, S.; Al-Manaseer, N.; Stietiya, M.H. Soil salinity changes in the Jordan Valley potentially threaten sustainable irrigated agriculture. Pedosphere 2013, 23, 376–384. [Google Scholar] [CrossRef]
  4. Liu, X.T. Management on Degraded Land and Agricultural Development in the Songnen Plain; Science Press: Beijing, China, 2001. (In Chinese) [Google Scholar]
  5. Wang, L.; Seki, K.; Miyazaki, T.; Ishihama, Y. The causes of soil alkalinization in the Songnen Plain of Northeast China. Paddy Water Environ. 2009, 7, 259–270. [Google Scholar] [CrossRef]
  6. Li, X.Y.; Wang, Z.M.; Song, K.S.; Zhang, B.; Liu, D.W.; Guo, Z.X. Assessment for salinized wasteland expansion and land use change using GIS and remote sensing in the west part of Northeast China. Environ. Monit. Assess. 2007, 131, 421–437. [Google Scholar] [CrossRef] [PubMed]
  7. Jordán, M.M.; Navarro-Pedreño, J.; Garcia-Sánchez, E.; Mateu, J.; Juan, P. Spatial dynamics of soil salinity under arid and semi-arid conditions geological and environmental implications. Environ. Geol. 2004, 45, 448–456. [Google Scholar] [CrossRef]
  8. Rakhimbaev, F.M.; Bezpalov, N.F.; Khamidov, M.K. Peculiarities of Crop Irrigation in Lower Amu Dary River Areas; Fan Iashkent: Tashkent, Uzbekistan, 1992. (In Russian) [Google Scholar]
  9. Ibrakhimov, M.; Khamzina, A.; Forkutsa, I.; Paluasheva, G.; Lamers, J.P.A.; Tischbein, B.; Vlek, P.L.G.; Martius, C. Groundwater table and salinity: Spatial and temporal distribution and influence on soil salinization in Khorezm region (Uzbekistan, Aral Sea Basin). Irrig. Drain. Syst. 2007, 21, 219–236. [Google Scholar] [CrossRef]
  10. Luo, J.M.; Yang, F.; Wang, Y.J.; Ya, Y.J.; Deng, W.; Zhang, X.P.; Liu, Z.J. Mechanism of soil sodification at the local scale in Songnen Plain, northeast China, as affected by shallow groundwater table. Arid Land Res. Manag. 2011, 25, 234–256. [Google Scholar] [CrossRef]
  11. Song, C.C.; Deng, W. Characters of groundwater and influence on the interior salt-affected soil in the west of Jilin Province. Sci. Geogr. Sin. 2000, 20, 246–250. (In Chinese) [Google Scholar]
  12. Yang, F.; Zhang, G.X.; Yin, X.R.; Liu, Z.J.; Huang, Z.G. Study on capillary rise from shallow groundwater and critical water table depth of a saline-sodic soil in western Songnen Plain of China. Environ. Earth Sci. 2011, 64, 2119–2126. [Google Scholar] [CrossRef]
  13. Zhang, G.X. Study on the Evolution of Groundwater Environment and Safety Management in the Taoer River Basin. Ph.D. Thesis, Chinese Academy of Sciences, Beijing, China, 2004. [Google Scholar]
  14. Zhou, Y.; Fu, Z.; Liu, D.; Wang, L.; Wang, L. A study on the processes of grassland degradation, land salinization and soil desertification in the west of Jilin province. J. Jilin Univ. Earth Sci. Ed. 2003, 33, 348–354. (In Chinese) [Google Scholar]
  15. Joeckel, R.M.; Clement, B.J.A. Soils, surficial geology, and geomicrobiology of saline-sodic wetlands, North Platte River Valley, Nebraska, USA. CATENA 2005, 61, 63–101. [Google Scholar] [CrossRef]
  16. Ardahanlioglu, O.; Oztas, T.; Evren, S.; Yilmaz, H.; Yildirim, Z.N. Spatial variability of exchangeable sodium, electrical conductivity, soil pH and boron content in salt- and sodium-affected areas of the Igdir plain (Turkey). J. Arid Environ. 2003, 54, 495–503. [Google Scholar] [CrossRef]
  17. Corwin, D.L.; Lesch, S.M.; Oster, J.D.; Kaffka, S.R. Monitoring management-induced spatio-temporal changes in soil quality through soil sampling directed by apparent electrical conductivity. Geoderma 2006, 131, 369–387. [Google Scholar] [CrossRef]
  18. Emadi, M.; Baghernejad, M. Comparison of spatial interpolation techniques for mapping soil pH and salinity in agricultural coastal areas, northern Iran. Arch. Agron. Soil Sci. 2014, 60, 1315–1327. [Google Scholar] [CrossRef]
  19. Zhao, C.W.; Xu, L.; Wang, Z.C.; Wang, Y.H.; Chi, C.M. Studying Soil Salinity and Its Relations with Microtopography and Vegetation at Field Scale. Pol. J. Environ. Study 2013, 22, 1227–1237. [Google Scholar]
  20. Odeh, I.O.A.; Chittleborough, D.J.; McBratney, A.B. Elucidation of soil-landform interrelationships by canonical ordination analysis. Geoderma 1991, 49, 1–32. [Google Scholar] [CrossRef]
  21. McKenzie, N.J.; Austin, M.P. A quantitative Australian approach to medium and small scale surveys based on soil stratigraphy and environmental correlation. Geoderma 1993, 57, 329–355. [Google Scholar] [CrossRef]
  22. Gui, D.W.; Lei, J.Q.; Zeng, F.J.; Runge, M.; Mu, G.J.; Yang, F.X.; Zhu, J.T. Ordination as a tool to characterize soil particle size distribution, applied to an elevation gradient at the north slope of the Middle Kunlun Mountains. Geoderma 2010, 158, 352–358. [Google Scholar] [CrossRef]
  23. TerBraak, C.J.F. CANOCO—A Fortran Program for Canonical Community Ordinatin by Partial, Detrended, Canonical Correspondence Analysis and Redundancy Analysis (Version 2.1); Technical Report LWA-88-2 TNO; Agricultural Mathematics Group: Wageningen, The Netherlands, 1988. [Google Scholar]
  24. Deng, W.; Qiu, S.W.; Liang, Z.W. Background of Regional Eco-Environment in Da’an Sodic Land Experiment Station of China; Science Press: Beijing, China, 2005. (In Chinese) [Google Scholar]
  25. Eaton, A.D.; Clesceri, L.S.; Rice, E.W.; Greenberg, A.E. (Eds.) Standard Methods for the Examination of Water and Wastewater, 21st ed.; APHA Publication: Washington, DC, USA, 2005.
  26. Mahmoodabadi, M.; Yazdanpanah, N.; Sinobas, L.R.; Pazira, E.; Neshat, A. Reclamation of calcareous saline sodic soil with different amendments (I): Redistribution of soluble cations within the soil profile. Agric. Water Manag. 2013, 120, 30–38. [Google Scholar] [CrossRef]
  27. TerBraak, C.J.F.; Šmilauer, P. CANOCO Reference Manual and Canocodraw for Windows User’s Guide: Software for Canonical Community Ordination (Version 4.5); Microcomputer Power: New York, NY, USA, 2002. [Google Scholar]
  28. Lei, Z.D.; Yang, S.X.; Xu, Z.R.; Vachaud, G. Preliminary investigation of the spatial variability of properties. J. Hydraul. Eng. 1988, 9, 10–21. (In Chinese) [Google Scholar]
  29. Yang, J.F.; Deng, W.; Zhang, G.X. Water balance and temporal variation in water content in alkaline sodic soils with shallow water table. Hydrogeol. Eng. Geol. 2007, 4, 61–69. (In Chinese) [Google Scholar]
  30. Kilic, K.; Kilic, S. Spatial variability of salinity and alkalinity of a field having salination risk in semi-arid climate in northern Turkey. Environ. Monit. Assess. 2007, 127, 55–65. [Google Scholar] [CrossRef] [PubMed]
  31. Rengasamy, P.; Olson, K.A. Sodicity and Soil Structure. Aust. J. Soil Res. 1991, 29, 935–952. [Google Scholar] [CrossRef]
  32. Yang, J.F.; Zhang, G.X. Effect of hydrological process change on soil salinization in western Songnen Plain. J. Arid Land Resour. Environ. 2010, 24, 168–172. (In Chinese) [Google Scholar]
  33. Yang, F.; Zhang, G.X.; Yin, X.R.; Liu, Z.J. Field-Scale Spatial Variation of Saline-Sodic Soil and Its Relation with Environmental Factors in Western Songnen Plain of China. Int. J. Environ. Res. Public Health 2011, 8, 374–387. [Google Scholar] [CrossRef] [PubMed]
  34. Chi, C.M.; Wang, Z.C. Soil salinity, sodicity and alkalinity characteristics of degraded grassland in the west of Jinlin province, China. Ecol. Environ. Sci. 2013, 22, 1534–1537. (In Chinese) [Google Scholar]
  35. Rengasamy, P. Transient salinity and subsoil constraints to dry land farming in Australian sodic soils: An overview. Aust. J. Exp. Agric. 2002, 42, 351–361. [Google Scholar] [CrossRef]
  36. Cemek, B.; Guler, M.; Kilic, K.; Demer, Y.; Arslan, H. Assessment of spatial variability in some soil properties as related to soil salinity and alkalinity in Bafra plain in northern Turkey. Environ. Monit. Assess. 2007, 124, 223–234. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Location of Da’an Sodic Land Experimental Station.
Figure 1. Location of Da’an Sodic Land Experimental Station.
Water 08 00227 g001
Figure 2. Distribution of measuring points and relative height in experiment Site.
Figure 2. Distribution of measuring points and relative height in experiment Site.
Water 08 00227 g002
Figure 3. Rainfall Regulations during Different Hydrological Years.
Figure 3. Rainfall Regulations during Different Hydrological Years.
Water 08 00227 g003
Figure 4. RDA 2D Ordination Diagram for Soil Salinization and Environment Factors in wet year.
Figure 4. RDA 2D Ordination Diagram for Soil Salinization and Environment Factors in wet year.
Water 08 00227 g004
Figure 5. RDA 2D Ordination Diagram for Soil Salinity and Environment Factors in dry year.
Figure 5. RDA 2D Ordination Diagram for Soil Salinity and Environment Factors in dry year.
Water 08 00227 g005
Table 1. Descriptive statistics on soil ions and salinity parameters in 2005 and 2007.
Table 1. Descriptive statistics on soil ions and salinity parameters in 2005 and 2007.
Chemical PropertyYearNumber of PointsMinimumMaximumMeanStandard DeviationCoefficient of Variation
EC (µs·cm−1)20054082.82710392 ± 63.2a4001.02
2007401234730660 ± 145.4a9201.39
pH2005407.8710.79.51 ± 0.13a0.810.09
2007408.2510.69.48 ± 0.11a0.690.07
SAR2005403.0986.8510.61 ± 2.0a12.591.19
2007405.55305.8935.92 ± 11.2b70.871.97
SC (mg/kg)200540587.713,6165212 ± 606a38331.25
200740873.822,8264069 ± 807a51050.74
Cation content (mg/kg)
Na+20054083.65062809 ± 123a7810.96
200740161.584991238 ± 312a19761.60
K+2005404.9465.627.5 ± 3.1a19.30.70
2007403.4617.88.89 ± 0.49b3.090.34
Mg2+2005407.76620108 ± 19.6a1241.14
2007400.5159.913.4 ± 1.78b11.30.84
Ca2+20054042.41515501 ± 71.4a4510.90
20074044571130 ± 16.9b1070.81
Anion content (mg/kg)
CO32−2005400.01386189 ± 62.9a2581.37
2007400.03240368 ± 118a7482.03
HCO320054036696993369 ± 408a25810.77
200740439.268811583 ± 237b15000.95
Cl200540711793198 ± 41.9a2651.33
2007401423186543 ± 127b8051.48
SO42−2005400.036.59.14 ± 1.07a6.780.74
20074034.76951183 ± 32.4b2051.12
Note: EC is electrical conductivity; SC is soil salinity content; SAR is sodium adsorption ratio; Means followed by the same letters within each column (lower case letters) are not significantly different at 0.05 level.
Table 2. Mean values of soil ions and index of salinity in different ground elevation.
Table 2. Mean values of soil ions and index of salinity in different ground elevation.
Ground Elevation (cm)YearNa+ mg/kgCO32− mg/kgEC µs/cmpHSC mg/kg
0–102005594 ± 70a103 ± 30a286 ± 23a9.18 ± 0.12a4444 ± 3675a
2007496 ± 275a90 ± 138a316 ± 35a9.18 ± 0.12a2044 ± 968b
10–202005844 ± 102a238 ± 178a397 ± 144a9.8 ± 0.22a6227 ± 3752a
20071013 ± 184a263 ± 215a617 ± 432a9.7 ± 0.14a3550 ± 1874b
20–4020052219 ± 1433a660 ± 379a1200 ± 462a10.4 ± 0.15a7324 ± 2952a
20077797 ± 520b2844 ± 310b3263 ± 766b10.5 ± 0.04a20,850 ± 1342b
Note: EC is electrical conductivity; SC is soil salinity content; Means followed by the same letters within each column (lower case letters) are not significantly different at 0.05 level.
Table 3. Intra-set correlations of the environmental variables, eigenvalue, and cumulative percentage variance and species-environment correlation coefficients for the first four axes of RDA.
Table 3. Intra-set correlations of the environmental variables, eigenvalue, and cumulative percentage variance and species-environment correlation coefficients for the first four axes of RDA.
Wet Year (2005)Dry Year (2007)
AX1AX2AX3AX4AX1AX2AX3AX4
Eigenvalue0.480.0360.0000.0000.6870.0010.0000.000
Species-environment correlation0.6950.7710.3240.1530.8350.4290.3100.249
Cumulative percentage variance (%)94.610010010099.8100100100
Ground elevation (H)0.4339 *0.5213 **−0.10010.00290.8092 **−0.0433−0.0349−0.0243
W0–10−0.0972−0.27980.1884−0.0888−0.29370.0064−0.11470.2058
W10–30−0.2478−0.05130.0897−0.0941−0.2904−0.0334−0.07570.0922
W30–60−0.2424−0.1314−0.0242−0.0808−0.3866 *−0.17470.16450.0853
W60–100−0.27060.17140.1879−0.0422−0.3564 *−0.01900.02910.0222
Ponding depth (h)−0.5281 **−0.06510.0825−0.0136−0.5032 **0.22170.11080.1025
Ponding time (d)−0.4864 **−0.11550.1106−0.0009−0.5532 **0.33120.08300.0852
Notes: * Correlation significant at the 0.05 level; ** Correlation significant at the 0.01 level. W0–10, W10–30, W30–60 and W60–100 are the soil moistures at 0–10, 10–30, 30–60, and 60–100 cm, respectively.
Table 4. Pearson correlation analysis between salinity and sodicity and environmental factors.
Table 4. Pearson correlation analysis between salinity and sodicity and environmental factors.
YearVariableNa+Ca2+CO32−HCO3ClSO42−SCpHECSAR
2005H0.728 **0.2030.744 **0.2890.688 **0.3100.471 **0.558 **0.733 **0.629 **
W30–60−0.286−0.140−0.306−0.200−0.266−0.057−0.252−0.256−0.269−0.214
W60–100−0.076−0.291−0.192−0.2920.0490.02−0.262−0.351 *−0.0310.122
h−0.452 **−0.429 **−0.530 **−0.481 **−0.324 *0.12−0.537 **−0.553 **−0.399 *−0.240
d−0.458 **−0.382 *−0.531 **−0.430 **−0.340 *0.06−0.498 **−0.560 **−0.413 **−0.265
2007H0.783 **0.1990.788 **0.805 **0.763 **0.768 **0.809 **0.736 **0.625 **0.784 **
W30–60−0.394 *0.064−0.352 *−0.337 *−0.429 **−0.407 **−0.386 *−0.370 *−0.225−0.396 *
W60–100−0.346 *−0.113−0.344 *−0.344 *−0.348 *−0.327 *−0.356 *−0.317 *−0.290−0.362 *
h−0.461 **−0.423 **−0.454 **−0.545 **−0.460 **−0.438 **−0.505 **−0.397 *−0.527 **−0.458 **
d−0.511 **−0.411 **−0.508 **−0.598 **−0.501 **−0.471 **−0.555 **−0.444 **−0.591 **−0.504 **
Notes: * Correlation significant at the 0.05 level; ** Correlation significant at the 0.01 level. W0–10, W10–30, W30–60 and W60–100 are the soil moistures at 0–10, 10–30, 30–60, and 60–100 cm, respectively.
Table 5. Ranking environmental variables in importance in 2005.
Table 5. Ranking environmental variables in importance in 2005.
Marginal EffectsConditional Effects
VariableLam.VariableLam.PF
h0.27h0.270.00213.79
d0.23H0.080.0244.58
H0.19W0100.040.0982.79
W601000.07W601000.050.1002.88
W10300.06d0.040.1182.71
W30600.06W10300.000.6660.22
W0100.01W30600.010.8260.10
Notes: P = significance level of the effect, as obtained with a Monte Carlo permutation test under the null model with 499 random permutations. W0–10, W10–30, W30–60 and W60–100 are the soil moistures at 0–10, 10–30, 30–60, and 60–100 cm, respectively.
Table 6. Ranking environmental variables in importance in 2007.
Table 6. Ranking environmental variables in importance in 2007.
Marginal EffectsConditional Effects
VariableLam.VariableLam.PF
H0.35H0.350.00219.27
d0.25W601000.100.0498.68
W30600.15h0.020.1062.69
W601000.13W30600.010.6500.21
W0100.09d0.000.7420.17
W10300.08W0100.000.9600.03
Note: W0–10, W10–30, W30–60 and W60–100 are the soil moistures at 0–10, 10–30, 30–60, and 60–100 cm, respectively.

Share and Cite

MDPI and ACS Style

Yang, F.; An, F.; Ma, H.; Wang, Z.; Zhou, X.; Liu, Z. Variations on Soil Salinity and Sodicity and Its Driving Factors Analysis under Microtopography in Different Hydrological Conditions. Water 2016, 8, 227. https://doi.org/10.3390/w8060227

AMA Style

Yang F, An F, Ma H, Wang Z, Zhou X, Liu Z. Variations on Soil Salinity and Sodicity and Its Driving Factors Analysis under Microtopography in Different Hydrological Conditions. Water. 2016; 8(6):227. https://doi.org/10.3390/w8060227

Chicago/Turabian Style

Yang, Fan, Fenghua An, Hongyuan Ma, Zhichun Wang, Xia Zhou, and Zhijun Liu. 2016. "Variations on Soil Salinity and Sodicity and Its Driving Factors Analysis under Microtopography in Different Hydrological Conditions" Water 8, no. 6: 227. https://doi.org/10.3390/w8060227

APA Style

Yang, F., An, F., Ma, H., Wang, Z., Zhou, X., & Liu, Z. (2016). Variations on Soil Salinity and Sodicity and Its Driving Factors Analysis under Microtopography in Different Hydrological Conditions. Water, 8(6), 227. https://doi.org/10.3390/w8060227

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop