Multivariate Analysis of Soil Salinity in a Semi-Humid Irrigated District of China: Concern about a Recent Water Project

: The Chaobai River (CBR) basin in northern China is experiencing an unprecedented continuous inﬂow of external water via the South–North Water Diversion Project, which has channeled water from the southern part of the country to the north. Consequently, the steady rise of groundwater table in recent years is threatening soil salinity regulation. The purpose of this study was to describe the status of salinity of the surface soil in the CBR basin and to evaluate the impact of environmental factors including groundwater table on the spatial distribution of soil salinity using multivariate analysis (MVA) technique. In this study, 10 chemical variables of soil samples collected in 204 sites along CBR were analyzed, considering their interaction with three environmental factors: the density of irrigation canals, groundwater depth and topography. Statistical analysis mainly consisted of principal component analysis (PCA), redundancy analysis (RDA) and clustering analysis (CA). The results allow deﬁning the surface soil in the CBR basin as a slightly saline and moderately alkaline media. The ﬁrst two axes of multivariate model approximately explains 51% of the observed variability and allows distinguishing two main domains: the saline and the alkaline. The variability of the saline domain, deﬁned by major cations and anions, is obviously controlled by macro environmental factors, of which the density of irrigation canals and groundwater depth contributes 71% and 28%, respectively, while that of the alkaline domain, related to pH and bicarbonate, mainly manifests as singular behaviors of soil groups like rice cultivation or sewage irrigation. The results suggests that more attention should be paid to the ongoing water table rise to help inform future land management decisions and to prevent a double threat of both groundwater and surface water on soil salinization. Meanwhile, this study shows the enormous potential of MVA technique, speciﬁcally the complementary duo of RDA and CA, for integrating both global and local information of soil salinity and environmental factors.


Introduction
The mother river of the Chinese national capital Beijing and Tianjin municipality-the Chaobai River (CBR), has given birth to generations of civilization and now is the important grain and vegetable base of the two megacities [1,2]. Nevertheless, the CBR basin was once known as a saline-alkaline, waterlogged, flood-prone and low-yielding region due to high groundwater table in the 1950s, and at that time, soil salinity constituted a major factor limiting crop production [3,4] because it affects plant growth and survival [5]. The period from 1950s onward for a few decades witnessed successful The plowing layer mainly consists of silty clay, clay silt and silt, and was formed in the Holocene. The primary type of surface soil is Chao soil, which can be divided into four subclasses: alkaline soil; meadow cinnamon soil; saline Chao soil and loamy calcareous soil [20]. In these conditions, irrigation practices began at least 2000 years ago at the riverbanks and were extended to include most of the land in the 1980s [21]. At the same time in the 1980s, a history of sewage irrigation began and continued until today. Nowadays, Tianjin has the largest area of sewage irrigation in China, followed by Beijing [22]. The ionic concentrations of irrigation water are relatively high and are increasing year-by-year [23], so a significant amount of irrigated lands are salt-affected. The main crops are wheat, corn and soybeans, where rice is only grown in the southeast part of the area [24].

Sampling and Laboratory Analyses
Sampling was carried out in July and August 2019, under dry soil conditions. A total of 204 soil samples were taken in the topsoil at 20-30 cm deep. The sampling sites were uniformly distributed and had good coverage of the study area, except some individual areas conditioned by their limited accessibility. An Eijkelkamp auger was adopted to obtain the soil samples, then they were air dried and sieved through a 2-mm mesh before laboratory analyses.
The suspensions were made by 1:5 soil-water extracts and all samples were analyzed for cations such as potassium (K + ), sodium (Na + ), calcium (Ca 2+ ) and magnesium (Mg 2+ ) and for anions such as sulfate (SO4 2− ), chloride (Cl − ), nitrate (NO3 -), carbonate (CO3 2− ) and bicarbonate (HCO3 − ). Soil pH were measured in 1:2.5 suspension in water at the same time. Standard methods were used for laboratory testing: Na + and K + (atomic emission spectrophotometry), Ca 2+ , Mg 2+ , SO4 2− , Cl − , NO3 − , CO3 2− and HCO3 − (absorption spectrophotometry with a continuum-flux autoanalyzer), pH (potentiometric determination). The plowing layer mainly consists of silty clay, clay silt and silt, and was formed in the Holocene. The primary type of surface soil is Chao soil, which can be divided into four subclasses: alkaline soil; meadow cinnamon soil; saline Chao soil and loamy calcareous soil [20]. In these conditions, irrigation practices began at least 2000 years ago at the riverbanks and were extended to include most of the land in the 1980s [21]. At the same time in the 1980s, a history of sewage irrigation began and continued until today. Nowadays, Tianjin has the largest area of sewage irrigation in China, followed by Beijing [22]. The ionic concentrations of irrigation water are relatively high and are increasing year-by-year [23], so a significant amount of irrigated lands are salt-affected. The main crops are wheat, corn and soybeans, where rice is only grown in the southeast part of the area [24].

Sampling and Laboratory Analyses
Sampling was carried out in July and August 2019, under dry soil conditions. A total of 204 soil samples were taken in the topsoil at 20-30 cm deep. The sampling sites were uniformly distributed and had good coverage of the study area, except some individual areas conditioned by their limited accessibility. An Eijkelkamp auger was adopted to obtain the soil samples, then they were air dried and sieved through a 2-mm mesh before laboratory analyses.
The suspensions were made by 1:5 soil-water extracts and all samples were analyzed for cations such as potassium (K + ), sodium (Na + ), calcium (Ca 2+ ) and magnesium (Mg 2+ ) and for anions such as sulfate (SO 4  The considered three environmental factors include topography, groundwater depth and the density of irrigation canals because they are intimately related to soil chemical properties. The topography and groundwater depth are, respectively in quantitative terms of elevation (meters above sea level) and the distance between the surface and groundwater table (negative value). Both of them were measured in the field survey, which conducted simultaneously with soil sampling. The depth of the wells selected for measurement is less than 80 m, so as to ensure that the data can represent the state of shallow groundwater. The groundwater depth at the sampling sites can be obtained by kriging interpolation [25]. The density of irrigation canals is defined as the area proportion of canals in the circle of radius one kilometer around sampling sites. The information of these canals were interpreted from the Worldview-4 remote sensing data (a resolution of 0.31 m) [26].

Statistical Analyses
A descriptive statistical analysis was conducted first, in which mean, median, maximum, minimum, standard deviation, coefficients of variation, skewness and kurtosis were calculated. Then Pearson correlation analysis was performed to assess the relationships and associations among variables.
A multivariate gradient analysis, combining ordination and multiple regression techniques, was carried out to integrate the information of spatial variability and relationships between soil salinity and environmental factors. The ordination is aimed at finding canonical axes that explain the maximum variability in the chemical properties of the samples. Unconstrained ordination (i.e., indirect gradient analysis) seeks the variables that could best explain the saline ions composition and takes them as ordination axes. On the other hand, in constrained ordination (i.e., direct gradient analysis), variability is confined to be only explained in terms of the environmental factors. The two approaches complement each other. The main gradient of variability is obtained via the former indirect method, while with the latter direct one only the variability associated with the environmental factors is involved.
According to the analysis of the gradient length [27], which antecedently conducted on the response dataset, a linear analysis model was selected and thus determining the two methods of PCA and RDA, respectively used for indirect gradient analysis and direct gradient analysis. Before analysis, the data were standardized by standard deviation for PCA and by error variance for RDA. Monte Carlo permutation tests were adopted to evaluate the signification of ordination axes.
CA was used in this study to explain the principal component axes that can't be interpreted effectively in RDA. To identify the type of spatial correlation, the Local Moran's I [28] is used as an indicator of local spatial association to provide a measure of association for each unit. The local measure of Moran's I is defined as: where x i is the normalized score of the principal component at a particular sample, x j is the normalized score of the principal component at another neighboring sample, x is the average normalized score, w ij is a weight that denotes the proximity between samples i and j, n is the number of samples, and S 2 is the variance of the normalized scores of the principal components. The Moran scatter plot helps identify areas where the point values are similar or dissimilar and categorizes them into the following clusters: low values surrounded by low values (LL), high values surrounded by high values (HH) and outliers, either a high value surrounded by low values (HL) or a low value surrounded by high values (LH) and the others without significant abnormality (O). Data standardization was completed in SPSS 24 (IBM, Chicago, IL, USA). The PCA and RDA were conducted with Canoco 5.0 (Microcomputer Power, New York, NY, USA) [29]. The Local Moran's I was performed with Arcgis 10.1 (ESRI, Redlands, CA, USA) [30].

General Features of the Soil Ionic Composition
Descriptive statistics for soil ionic composition are shown in Table 1. The average concentration soil salt is 1.1 g/kg, meeting the standard of slight salinization [31]. The average concentrations of major anions and cations decrease as follows: bicarbonate > sulfate > chloride > nitrate and sodium > calcium > magnesium > potassium. The first three cations are the most abundant in salt affected soils [32,33]. The concentration of bicarbonate, which has the highest level among the anions, is approximately in accordance with the records in the Luan River delta 60 km east of the study area [34] and the Manas River watershed in northwestern China [35] and is larger than many other regions in the world [15,17,36]. The accumulation of bicarbonate leads to a trend of soil alkalization [37], so the measured soil pH ranges mainly between 7.5 and 8.5, which is a typical range of soil pH in the North China plain [38], up to 9.85. Therefore, according to the threshold values of soil chemical variables set by Chinese soil society [31], the soil can be described as slightly saline and moderately alkaline. Generally high coefficients of variation (CV) were observed in the data with the exception of pH. Agreeing with the results from previous studies, CV is particularly high in the case of NO 3 − and K + due to the uneven application of chemical fertilizer [39,40], although the average concentrations of them are lower than in the adjacent area [34]. From a global perspective, the high variability of NO 3 − seems to be widespread, which not only existed in irrigation districts, but also in natural wetland with little anthropogenic activities [15]. All anions and cations show positive skewness because of the presence of large values. This suggests an either natural or anthropogenic enrichment of main ions in the surface soil [41]. The observed asymmetry and the high kurtosis coefficients indicate the non-normal distribution of the data, which is verified by the Kolmogorov-Smirnov test. Only pH appears to be approximately normal distribution. The non-normality of the distribution of soil chemical properties is more of a rule than an exception due to the heterogeneity of the soil system and its complex interactions with the environmental factors, as has been demonstrated by other researchers [42,43]. Logarithmic transformation can partially overcome the general absence of symmetry in main ions [44].
The results of bivariate correlation analysis are shown in Table 2. Strong correlations generally exist among major anions and cations. The highest correlation is found between sodium and chloride (r = 0.86), comparable to that found in other works [32,45,46]. Potassium is the cation least correlated with the rest of major anions and cations (r < 0.2), indicating that in highly saline points dominated by Water 2020, 12, 2104 6 of 14 sodium and calcium, the concentration of potassium is low because the two of sodium and calcium have high correlations with the major anions. pH was slightly negatively correlated with the Ca 2+, but strongly and positively correlated with bicarbonate and carbonate, following the law of chemical reactions among the three ions when pH ranges between 8 and 10 [47].

Main Gradients in the Soil Ionic Composition
The contribution of PCA can be measured by the eigenvalues (λ i ) of principal component axes. The first four axes explain 75% (λ 1 = 0.34, λ 2 = 0.17, λ 3 = 0.14, λ 4 = 0.10) of the total variance. Table 3 shows the loadings of soil variables on the four principal components axes. Considering small contributions of axis 3 and axis 4 to the total variance, as well as the absence of actual significance of axis 3, the following analysis would focus only on the first two axes.  Figure 2 shows the ordination results of the first two axes. Correlations can be evaluated by perpendicularly projecting the arrows of the other variables onto those of the targeted variables [48]. In this case, two different chemical environments can be identified in the soils of the CBR basin: the saline domain and the alkaline domain. The former is defined by major cations and anions, which assigned high positive value on axis I. The latter is related to pH, bicarbonate and carbonate, which take main positive loadings of axis II. The gathering of the arrows within the two domains agrees with the already commented correlation between the chemical variables. Moreover, the separation of potassium from the other ions is also observed, further confirming their low correlation. Therefore, the biplot synthesizes the results of the correlation analysis. The unconstrained variability explained by the first two principal component axes (51%) can be comparable with the 54% observed by Vestin et al. (2006), considering that their work was carried out in a homogenous forest growing in cold climates, but a bit poorer than the 59% in a wetland [15], the 64% in an irrigated lowland [32] and the 69% in an irrigated hilly region [17]. It is not difficult to find that good homogeneity, little artificial interference and closeness of the soil system can be beneficial to the gathering of variability of soil chemical properties toward the principal component axes [15]. The CBR basin is more than 15 times the study area of any other researches mentioned above and has active exchange of surface water and groundwater with not only the adjacent basin, but also remote ecological system (Figure 1). Moreover, the two megacities of Beijing and Tianjin bring strong anthropogenic activities and intensively changed the pedologic characteristic in topsoil. All of these led to the relative low percentage of explained variability. Even so, the significances that the first two axes in this and previous studies stand for were similar. The first axis were highly correlated with soil salinity (electrical conductivity and main ions), while their second axis were linked to alkalinity and nutrients (pH, bicarbonate and dissolved organic carbon).
In general, the samples are distributed evenly around the origin of the coordinate, except for some individual points showing high scores on both two axes. The distance between the samples is a measure of the degree of their dissimilarity in chemical composition, and the farther the distance, the greater the difference. Viewed this way, the observed commonly proximity between samples can be attributed to the fact that they have been taken in surface soil in the same basin. The unconstrained variability explained by the first two principal component axes (51%) can be comparable with the 54% observed by Vestin et al. (2006), considering that their work was carried out in a homogenous forest growing in cold climates, but a bit poorer than the 59% in a wetland [15], the 64% in an irrigated lowland [32] and the 69% in an irrigated hilly region [17]. It is not difficult to find that good homogeneity, little artificial interference and closeness of the soil system can be beneficial to the gathering of variability of soil chemical properties toward the principal component axes [15]. The CBR basin is more than 15 times the study area of any other researches mentioned above and has active exchange of surface water and groundwater with not only the adjacent basin, but also remote ecological system (Figure 1). Moreover, the two megacities of Beijing and Tianjin bring strong anthropogenic activities and intensively changed the pedologic characteristic in topsoil. All of these led to the relative low percentage of explained variability. Even so, the significances that the first two axes in this and previous studies stand for were similar. The first axis were highly correlated with soil salinity (electrical conductivity and main ions), while their second axis were linked to alkalinity and nutrients (pH, bicarbonate and dissolved organic carbon).
In general, the samples are distributed evenly around the origin of the coordinate, except for some individual points showing high scores on both two axes. The distance between the samples is a measure of the degree of their dissimilarity in chemical composition, and the farther the distance, the greater the difference. Viewed this way, the observed commonly proximity between samples can be attributed to the fact that they have been taken in surface soil in the same basin.

Saline Domain
In the RDA, ordination was constrained to be in terms of the three environmental factors: the density of irrigation canals, groundwater depth and topography (the information of them is shown in Figure 1 and Table 1). The first two ordination axes in the RDA account for 21% of the variance in the data, dominated by the first axis (λ 1 = 0.20, λ 2 = 0.01, Figure 3). It means that the environmental factors account for 59% and 6% of the unconstrained variability explained by PCA respective in axes I and II, the significance of which is proved by the Monte Carlo tests based on 999 permutations. According to these results, it can be inferred that the saline domain in PCA could be well explained by the environmental factors, whereas the alkaline domain can't.
Water 2020, 12, x FOR PEER REVIEW 8 of 14

Saline Domain
In the RDA, ordination was constrained to be in terms of the three environmental factors: the density of irrigation canals, groundwater depth and topography (the information of them is shown in Figure 1 and Table 1). The first two ordination axes in the RDA account for 21% of the variance in the data, dominated by the first axis (λ1 = 0.20, λ2 = 0.01, Figure 3). It means that the environmental factors account for 59% and 6% of the unconstrained variability explained by PCA respective in axes Ⅰ and Ⅱ, the significance of which is proved by the Monte Carlo tests based on 999 permutations. According to these results, it can be inferred that the saline domain in PCA could be well explained by the environmental factors, whereas the alkaline domain can't. Saline ions tend to accumulate in the areas where irrigation canals densely distributed ( Figure  3). The contribution of canals to axis I is 71% (p = 0.002), the most of the three factors. Given the history of sewage irrigation derived from 1980s in the CBR basin [22], a close relationship between the chemical characters of irrigation water and topsoil can be inferred. A long-term record of chemical characteristic of irrigation water in a canal in the middle reach of CBR shows sustained high value of bicarbonate content [49], which in accordance with soil tests. Continuous increase of chloride, sulfate and sodium between 1980-2007 was also observed. Scholars attributed this increase to domestic and industrial sewage water, during which the population exploded, and industry rapidly developed. It is worth noting that sulfate saw the largest increase, approximate 50%, among all of the ions, which may explains the accumulation of sulfate in the soil. The reason for its large rise is probably related to the absence of regulation for limiting sulfate in national standard (Chinese quality of farmland irrigation water, GB20922-2007). In fact, even though the emission standard of TDS was specified Saline ions tend to accumulate in the areas where irrigation canals densely distributed (Figure 3). The contribution of canals to axis I is 71% (p = 0.002), the most of the three factors. Given the history of sewage irrigation derived from 1980s in the CBR basin [22], a close relationship between the chemical characters of irrigation water and topsoil can be inferred. A long-term record of chemical characteristic of irrigation water in a canal in the middle reach of CBR shows sustained high value of bicarbonate content [49], which in accordance with soil tests. Continuous increase of chloride, sulfate and sodium between 1980-2007 was also observed. Scholars attributed this increase to domestic and industrial sewage water, during which the population exploded, and industry rapidly developed. It is worth noting that sulfate saw the largest increase, approximate 50%, among all of the ions, which may explains the accumulation of sulfate in the soil. The reason for its large rise is probably related to the absence of regulation for limiting sulfate in national standard (Chinese quality of farmland irrigation water, GB20922-2007). In fact, even though the emission standard of TDS was specified (less than one gram per liter), substandard salt was often detected in the processed sewage water [22]. Considering the significance for water conservation and regulation, sewage irrigation will no doubt continue. In this condition, it is essential to improve the supervision of wastewater treatment plant. Some positive changes were predicted that the water with low salinity in Danjiangkou Reservoir would dilute the surface water in CBR basin after South-North Water Diversion Project launched [50]. However, the dilution has not been observed so far.
Topography has the opposite effect of groundwater depth on soil salinity, so their arrows are drawn in a straight line in opposite directions (Figure 3). The groundwater depth contributes 28% to axis I (p = 0.002) and the topography contributes 1% (p = 0.01). The little percentage of the topography can be attributed to its strong correlation with the groundwater table [51]. These results have much in common with previous studies carried out in other different landforms [52,53]. These researches demonstrated that, in most cases, the distribution of soil salinity under natural conditions is in connection with the location of the soil in the topographic unit, ranging from the high and sloping fields having deep groundwater table without salt-affected soils to the moderately to strongly saline soils in the bottoms having shallow groundwater table [17].
According to recent monitoring, TDS of the shallow groundwater in the study area ranges from 0.5 to 1 g/L [2], which seems safe compared to surface water. However, there are threat of groundwater to soil salinization, more depending on evaporation intensity. Once the groundwater table reach salt-returning line, such groundwater is capable to provide much salt to topsoil just like it did in the 1950s [50]. In the past few decades, the influence of groundwater table on soil salinity in CBR basin has been of little concern, because most areas have a safe groundwater level which benefited from long-term overexploitation of groundwater. However, the present situation seems somewhat different. According to previous studies, the critical depth of groundwater in salinization of soil in western China is believed to range between 2.5-2.7 m [54,55], while about 3.5 m [56] in eastern China where the CBR basin belongs to. In our investigation, of groundwater depth, 25% of the data are less than 3.5 m and the ones less than 6 m account for about half. Moreover, the results would have been more severe if the survey were conducted in winter, the peak season of groundwater table in a year [57]. What seems certain is that the groundwater level will keep rising in the coming years due to the continuous inflow of water by South-North Water Diversion Project, and there will be more soils at the risk of salinization. Therefore, managers must find the right balance between their stand of groundwater and salinization of soil.

Alkaline Domain
As the RDA had shown, the three environmental factors, characterized by macroscopy, were incapable to explain the variation of the soil alkaline domain. Therefore, CA was conducted to help distinguish groups of soils having similar behaviors. On the basis of the scores of axis II in PCA, CA divided all the samples into three categories including 4 clusters: LL (Cluster 1 and Cluster 2 with relatively weak alkalinity), HH (Cluster 3 with relatively strong alkalinity) and O (the others without significant abnormality) (Figure 4a, Figure 5).
As the RDA had shown, the three environmental factors, characterized by macroscopy, were incapable to explain the variation of the soil alkaline domain. Therefore, CA was conducted to help distinguish groups of soils having similar behaviors. On the basis of the scores of axis II in PCA, CA divided all the samples into three categories including 4 clusters: LL (Cluster 1 and Cluster 2 with relatively weak alkalinity), HH (Cluster 3 with relatively strong alkalinity) and O (the others without significant abnormality) (Figure 4a, Figure 5).   Cluster 1 is distributed in the lowest part of the entire basin (Figure 4a,b). The average pH in this group was well below its overall level ( Figure 5). It is worth mentioning that the upper reaches of this cluster are adjacent to a perennial rice-growing area. Since the pH value suitable for the growth of rice is less than 7 [58], particularly 4.5-5.5 for rice seedling [59], the application of acidulent is essential to prevent some diseases against the background of alkalinity. Obviously, the acidulent reduced soil pH and this influence expanded downstream with the acid flowing in canals and groundwater.
Rice growing and alkalinity of soil affect each other and that varies by location. In China, the dominant areas of this crop are mainly distributed to the south of the Yangtze River and in the Northeast Plain, where the soils are naturally acidic and suitable for rice cultivation [60]. By contrast, in the regions with moderate alkaline soil like north china plain, manual intervention on soil pH is necessary to ensure good health of seedlings [59]. Some scholars have already paid attention to the influence of application of acidulants on soil pH but find some limits during the acidization [58] due to the presence of pH buffer system existed in paddy soils under the persistent flooding condition. Consequently, after long-term planting of rice, the soil pH generally tends to be neutral [58], which is consistent with our results (the average pH value of Cluster 1 is 7.04, Figure 5). In some cases, rice also grown in highly saline and alkaline soils, for example, Flumen in Spain [17] where rice is the only profitable crop. Nonetheless, rice cultivation here seems powerless to reduce pH because of inherent strong alkalinity. Given the above, although it matters, rice-growing cannot simply be an indicator of acidity. Cluster 2 is characterized by linear distribution along the Huangsha canal near the western border of lower reaches (Figure 4c). Therefore, the abnormality of this cluster, high concentrations of major ions such as sodium, potassium, chloride and sulfate ( Figure 5), may have something to do with this canal. On one hand, Huangsha canal is a branch of the Longfeng River (Figure 4c), one of Cluster 1 is distributed in the lowest part of the entire basin (Figure 4a,b). The average pH in this group was well below its overall level ( Figure 5). It is worth mentioning that the upper reaches of this cluster are adjacent to a perennial rice-growing area. Since the pH value suitable for the growth of rice is less than 7 [58], particularly 4.5-5.5 for rice seedling [59], the application of acidulent is essential to prevent some diseases against the background of alkalinity. Obviously, the acidulent reduced soil pH and this influence expanded downstream with the acid flowing in canals and groundwater.
Rice growing and alkalinity of soil affect each other and that varies by location. In China, the dominant areas of this crop are mainly distributed to the south of the Yangtze River and in the Northeast Plain, where the soils are naturally acidic and suitable for rice cultivation [60]. By contrast, in the regions with moderate alkaline soil like north china plain, manual intervention on soil pH is necessary to ensure good health of seedlings [59]. Some scholars have already paid attention to the influence of application of acidulants on soil pH but find some limits during the acidization [58] due to the presence of pH buffer system existed in paddy soils under the persistent flooding condition. Consequently, after long-term planting of rice, the soil pH generally tends to be neutral [58], which is consistent with our results (the average pH value of Cluster 1 is 7.04, Figure 5). In some cases, rice also grown in highly saline and alkaline soils, for example, Flumen in Spain [17] where rice is the only profitable crop. Nonetheless, rice cultivation here seems powerless to reduce pH because of inherent strong alkalinity. Given the above, although it matters, rice-growing cannot simply be an indicator of acidity. Cluster 2 is characterized by linear distribution along the Huangsha canal near the western border of lower reaches (Figure 4c). Therefore, the abnormality of this cluster, high concentrations of major ions such as sodium, potassium, chloride and sulfate ( Figure 5), may have something to do with this canal. On one hand, Huangsha canal is a branch of the Longfeng River (Figure 4c), one of the three most rivers loaded by sewage in the basin [22]. The industrial and domestic wastewater in this river showed typical features of high salinity and acidity, which is highly consistent with the chemical statistics of Cluster 2. On the other hand, a brick factory in operation was found to be located upstream of Cluster 2 through the interpretation of remote sensing, confirmed by field survey. In order to obtain materials for brick making, several gigantic pits were excavated to carry an exposure of groundwater so that industrial pollutants can easily enter into groundwater and spread downstream. However, further investigation is suggested to verify the above views.
Compared with the first two groups of soil, Cluster 3 is the one with stronger alkalinity, distributed in the lower reaches. It is difficult to determine whether its identification is due to human activities, because a slightly natural enrichment of alkalinity originally tends to appear in the bottom of the basin with the accumulation of alkalis like sodium ( Figure 3). In this condition, the spatial relationships that Cluster 3 is surrounded by the other two clusters with opposite characteristics would amplify its abnormality. Statistically, the pH value of category HH is not so much higher than category O, while the concentrations of bicarbonate almost equal ( Figure 5). Therefore, Cluster 3 is suggested to be a joint result of natural enrichment of alkalinity and spatial relationships. Given the above, the results of the CA stress the power of the MVA method for detecting singular behaviors of small groups of observations.

Conclusions
The present article reveals the complexity of spatial distribution of the soil salinity in the CBR basin, China. These soil chemical variables differ depending on the irrigation, groundwater table, terrain, irrigation regulation and the characteristics of the crops. The first two axes of PCA approximately explains 51% of the observed variability and allows distinguishing two main domains: the saline and the alkaline. In the RDA, the variability of the saline domain is obviously controlled by macro environmental factors, of which the density of canals and groundwater depth, respectively contributes 71% and 28%. In contrast, the variability of the alkaline domain mainly manifests as singular behaviors of soil groups like rice-growing or sewage irrigation, which is detected by CA.
The last decades witnessed the transformation of the source of soil salt from groundwater to surface water. A double threat of them on soil salinization is coming in next few years due to the ongoing rise of groundwater table. Therefore, managers must find the right balance between their stand of groundwater and salinization of soil.
Sustainable agriculture in irrigated districts requires information about the salinity of the soils, incorporating many influence factors to achieve sound interpretations. The moderate percent of the variances explained in this study suggests that other unconsidered environmental and management factors should be considered in future research. Even so, MVA, specifically the complementary duo of RDA and CA, no doubt has enormous potential for integrating both global and local information to improve our understanding of linkages among soil salinity, environmental factors and management practices.