Next Article in Journal
Positioning Locality Using Cognitive Directions Based on Indoor Landmark Reference System
Previous Article in Journal
A Monolithic Electrochemical Micro Seismic Sensor Capable of Monitoring Three-Dimensional Vibrations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Mapping Soil Salinity/Sodicity by using Landsat OLI Imagery and PLSR Algorithm over Semiarid West Jilin Province, China

1
Northeast Institute of Geography and Agroecology, Chinese Academy of Sciences, Key Laboratory of Wetland Ecology and Environment, Changchun 130102, China
2
University of Chinese Academy of Sciences, Beijing 100049, China
3
Northeast Normal University, Key Laboratory for Vegetation Ecology Science of Ministry of Education, Changchun 130021, China
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Sensors 2018, 18(4), 1048; https://doi.org/10.3390/s18041048
Submission received: 6 February 2018 / Revised: 14 March 2018 / Accepted: 29 March 2018 / Published: 31 March 2018
(This article belongs to the Section Remote Sensors)

Abstract

:
Soil salinity and sodicity can significantly reduce the value and the productivity of affected lands, posing degradation, and threats to sustainable development of natural resources on earth. This research attempted to map soil salinity/sodicity via disentangling the relationships between Landsat 8 Operational Land Imager (OLI) imagery and in-situ measurements (EC, pH) over the west Jilin of China. We established the retrieval models for soil salinity and sodicity using Partial Least Square Regression (PLSR). Spatial distribution of the soils that were subjected to hybridized salinity and sodicity (HSS) was obtained by overlay analysis using maps of soil salinity and sodicity in geographical information system (GIS) environment. We analyzed the severity and occurring sizes of soil salinity, sodicity, and HSS with regard to specified soil types and land cover. Results indicated that the models’ accuracy was improved by combining the reflectance bands and spectral indices that were mathematically transformed. Therefore, our results stipulated that the OLI imagery and PLSR method applied to mapping soil salinity and sodicity in the region. The mapping results revealed that the areas of soil salinity, sodicity, and HSS were 1.61 × 106 hm2, 1.46 × 106 hm2, and 1.36 × 106 hm2, respectively. Also, the occurring area of moderate and intensive sodicity was larger than that of salinity. This research may underpin efficiently mapping regional salinity/sodicity occurrences, understanding the linkages between spectral reflectance and ground measurements of soil salinity and sodicity, and provide tools for soil salinity monitoring and the sustainable utilization of land resources.

1. Introduction

Basically, mapping and assessing the soil salinity and sodicity is a first step for salinity monitoring, salt-soil agriculture development, and environmental sustainability maintenance throughout the world [1,2]. Such work is also valuable to the scientific understanding of soil salinity and sodicity occurrences and pursuing effective measures for evaluating the affected soils [3].
The west Jilin Province of China, typically a semiarid area, is among the world’s three largest sodic saline-alkali regions. Due to overgrazing and grassland conversion to farmland, its area of salt-affecting lands increased drastically during the 1950s and the 1990s [4,5]. Consequently, soil salinity and sodicity decreased the productivity of the croplands and affected diversity of the grassland plant species, which thwarted ecological sustainability and societal and economic development of the west Jilin Province. For better understanding the occurring extent and severity of soil salinity and sodicity, depicting accurate projection of future salinity dynamics and making decisions about land management and utilization in the region, mapping the local salinity, and sodicity is much needed [6].
Soil salinity refers to the amount of soluble salt conserved in the soils, whereas sodicity points to the sodium as often means poor soil structures that impede soil hydraulic conductance and water uptake of plants [1,7]. As soil salinity and sodicity aggravates, more salt ions appear at the soil surface, which favors the use of conventional remote sensing (RS) technology for surveillance [8]. In the last decades, multispectral and microwave RS data with varied spatio-temporal resolutions (e.g., Landsat, SPOT XS, IKONOS, ASTER, MODIS, IRS, and Radar, etc.) have been used for monitoring [9,10,11] and mapping [12,13] salt-affected soils and halophytic vegetation [3,14]. In most cases, multiple bands were transformed into one index that was more sensitive to soil salinity than using a single band [15,16]. As for the single band, reflectance at the visible (VIS) bands (400~750 nm) was found to be more sensitive to soil pH than in other spectral ranges; Landsat 8 OLI-band1 (coastal) was also found closely correlated with soil salinity [17]. Nonetheless, substantial elements in linking RS spectral information and field salinity measurements are still to be elucidated, particularly across diverse regions. In particular, soil salinity was more often mentioned than sodicity in previous studies [18,19]. Yet, discerning the disparity effects between salinity and sodicity could underpin full-view recognition of soil health and making decisions about land restoration, and improve the comprehensive utilization of the lands. But, this critical issue should be fully examined especially in reference to the soils with complicated coexistence of soil salinity and sodicity in the west Jilin Province.
The Partial Least Squares Regression (PLSR) has been widely used in soil parameter estimation [3]. This technique is a robust multivariate regression method appropriate particularly when the predictors exhibit multicollinearity [20]. Accordingly, many studies reported success of using PLSR for soil salinity assessment [21,22,23]. By contrast, Janik et al. [24] and Pang et al. [25] implied that the ANN (Artificial Neutral Network) method could outperform the basic PLSR for prediction of some soil properties and non-linear regression may be more accurate than linear regression. However, ANN is difficult to accurately analyze the various performance indices of a biblical network. Therefore, we chose PLSR to establish the inversion models of soil salinity and soil sodicity. Another consideration to improve the accuracy of modelling was that we tried to transform spectral reflectivity through non-linear functions before modelling. Here, the advantages of using the PLSR approach rely on its inference capabilities for modelling the probable linear relationships that are existing between spectral reflectance and salt content in soils, and power to include multiple response variables simultaneously for addressing the likely collinearity among the noise-engaged predictor variables.
Accordingly, the objectives of this research were as follows: (1) to build up the relationships between reflectance bands, spectral indices (vegetation indices and salinity/sodicity indices), and field measurements of soil salinity/sodicity; (2) to improve soil salinity/sodicity estimating accuracy by optimizing PLSR; (3) to map soil salinity/sodicity of the study area, particularly the HSS; and, (4) to unravel the occurring severity and distribution of soil salinity, soil sodicity, and HSS with regard to land cover and soil types. This research may provide a clue for pursuing tools about soil salinity monitoring, understanding correlations between spectral reflectance and field measurements of salinity, as well as tackling regional land degradation and development issues.

2. Materials and Methods

2.1. Study Area Descriptions

The west Jilin Province, as the largest sodic saline-alkali soil distribution in Northeast China, is located in the southern Songnen Plains (43°59′~46°18′ N, 121°38′~126°11′ E. Figure 1), with low-lying terrains and elevations between 100~200 m. The groundwater depth is about 2 m. The study area is under a temperate, semi-arid continental monsoon climate, with precipitation dwindling from 420~460 mm in the east to 350~420 mm in the west, and evaporation increasing from 1200~1600 mm in the east to 1500~2000 mm in the west [26,27]. About 80% of the total precipitation occurs between mid-June and mid-August. The soils are Chernozem (Haplic Chernozem, FAO), meadow soil (Eutric Vertisol, FAO), aeolian soil (Arenosol, FAO), Solonetz (Solonetz, FAO), and Chestnut soil (Haplic Kastanozem, FAO). Pristine vegetation of this area is temperate steppes grown on the salted soils, with widespread fragmented scar-like saline lands that were spotted on the landscapes. However, grassland conversion to cropland has been common since the 1980s in this region, leading to less than one-quarter of the natural steppes are merely remained and cultivated croplands are ubiquitous with salinity to some extent.

2.2. Data Sources and Pre-Processing

Six scenes of satellite images created by Landsat 8 OLI covering the study region were used, of which five scenes were acquired on August 2016, and one image at the end of July. The images were geo-rectified to the Universal Transverse Mercator (UTM) coordinate system using World Geodetic System (WGS) 1984 datum assigned to north UTM Zone 51. Atmospheric correction was performed using FLAASH Model (Fast Line-of-sight Atmospheric Analysis of Spectral Hypercubes) that can correct both additive and multiplicative atmospheric effects [28]. Using ENVI 5.3, the study region was clipped from a mosaic of all the processed images. The spectral indices that link vegetation performance and soil salinity were calculated by the reflectance of corresponding bands, including NDVI, EVI, SIs [2,15,29], DVI [30], SAVI, and SRSI [31]. The spectral indices are listed in Table 1.
For enhancing the remotely sensed data that built the linkages between spectral signatures and soil salinity and verifying the inversion models on the retrieval of soil salinity and sodicity, we conducted field measurements and soil sampling in October 2016. As a result, all of the soil samples were collected in triplicate from the 0–5 cm layer. The number of soil samples obtained from each land cover type is shown in Figure 1. All of the samples were air-dried, sieved through a 2 mm sieve. Soil pH and EC were measured in 1:5 soil:water extracts using LEICI pH5-3E and LEICI DDS-307A meters, respectively [17]. Wetland, water body, and built-up land were not involved in the inversion because we focused on the predominantly terrestrial salinity, and thus excluded samples from the lands. By excluding four data points with aberrant values following experimental analysis, we finally accepted 157 sampling points for analysis.
Auxiliary data included Digital Elevation Model (DEM, resolution 30 m) data; a soil type map (1:1000,000), provided by Northeast Institute of Geography and Agroecology, Chinese Academy of Sciences; and, reference/validation site for classification of land cover by field survey or from Google Earth in 2016.

2.3. Random Forest to Classifying Land Cover

The RS images were classified by Random Forest (RF) to draw a land cover map for analyzing soil salinity and sodicity within each land cover type. RF is an ensemble of classification trees, where each tree contributes with a single vote for the assignment of the most frequent class to the input data. Key advantages of RF include non-parametric nature, high classifying accuracy, and capability to weigh variable importance [32,33]. RF uses a random subset of inputs or predictive variables in a node division instead of using best variables, which can reduce the generalizing error and finally lead to highest accuracy. Here, the classification system consisted of seven basic land cover classes: woodland, grassland, wetland, water body, barren land, cropland, and built-up [34,35]. The overall land cover classification accuracy was 95.76%, with a Kappa value 0.91 for RF.

2.4. The PLSR

The PLSR is an extension of multiple regression analysis, in which the effects of linear combinations of several predictors on a response variable (or multiple response variables) are analyzed. The method combines the traits of principal component analysis and multiple linear regressions [20]. Therefore, to eliminate the influence of collinearity among variables on multiple linear regression models, PLSR is usually used to improve the stability of estimators [36]. For developing the PLSR algorithms in this research, soil EC and pH were used each as an indicator of soil salinity and sodicity. Based on field sampling, the PLSR models can be established for soil salinity and sodicity retrieval [3,16,36]. The variables were transformed into six formulas to better retrieve soil EC and pH before three times of PLSR analyses. For the first time, soil EC and pH were retrieved merely by the reflectance of bands. Subsequently, we only used spectral indices to build up the retrieval models. Finally, we combined the bands and spectral indices to retrieve soil EC and pH. By using R2 and RMSE, the models were assessed to determine best models. Consequently, the optimized latent variables for regression were obtained from the models for soil salinity and sodicity. Model performance can be evaluated with R2, Bias, Standard Deviation (SD), and RMSE. R2 indicates the strength of statistical correlation between measured and predicted values. The model can be accurate for R2 > 0.91, good for 0.82 < R2 < 0.90, moderate for 0.66 < R2 < 0.81, and poor for 0.5 < R2 < 0.65, according to Farifteh et al. [37]. Bias measures the mean difference between prediction and measurements, and SD represents the random component of total uncertainties [38]. The metrics are described, as follows:
R 2 = 1 ( y ^ i y i ) 2 ( y i y ¯ ) 2
Bias = ( y ^ i y i ) N
SD = [ ( y ^ i y i Bias ) 2 N 1 ] 1 / 2
RMSE = [ ( y ^ i y i ) 2 N ] 1 / 2
where y denotes the measures with a mean value of y ¯ , y ^ denotes the predicted values, and N the number of samples. The 0.05 probability level was used to test the significance of correlation.

3. Results

3.1. Sensitivity of RS Variables to Soil EC and pH

In order to distil the variables sensitive to soil pH and EC, spectral analyses of soil samples were implemented. Regression coefficients were used to indicate the impact of independent variables on the dependent variable, for instance, the higher the regression coefficient value, the greater the impact on the dependent variable [3]. Before modelling, the 157 soil samples were divided into two subsets. The first set (110 samples) was used for calibration and the second set (47 samples) for validation. Laboratory analysis showed that the pH values of the sampled soils ranged from 5.90 to 10.45, with an average 8.70 that generally pointed to sodicity [1]. EC measurements of the samples ranged from 0.03 to 3.99 mS/cm, with an average 0.54 mS/cm that indicated lower soluble salt concentration. Thus, most soils in the study area were more affected by sodicity rather than by salinity.
Table 2 shows the correlation coefficients between EC, pH, and reflectance of bands, and between bands. EC showed the lowest coefficient positive with NIR (R = 0.348), moderately correlated with SWIR1, SWIR2, and PAN (R = 0.704, 0.738, and 0.763), and higher with Coastal, Red, Green, and Blue (R = 0.821, 0.810, 0.826, and 0.818). In particular, the correlation coefficient between Red vs. pH (R = 0.805), or between EC vs. Green (R = 0.826) was the highest each among the group, implying that the band Red had priority for soil salinity retrieval and the band Green for sodicity.
Table 3 indicates the correlation coefficients between soil EC, pH, and all of the spectral indices one another. Vegetation indices (NDVI, EVI, DVI, and SAVI) were negatively correlated with EC and pH, with NDVI/SAVI the most negatively correlated. By contrast, salinity indices (SI, SI2, SI3, and SI4) were positively correlated with EC and pH. Since SI (R = 0.803) and SAVI (R = 0.814) exhibited more sensitive to soil sodicity than other spectral indices, thus we stipulated that SAVI was the priority spectral index to reflect soil sodicity. For EC, SI3 was the best estimator of soil salinity (R = 0.826), yet EVI had the poorest (R = −0.449). Overall, the RS-derived vegetation indices were more sensitive to soil sodicity than to soil salinity.

3.2. The PLSR Models for Soil Salinity/Sodicity Retrieval

We conducted mathematical transformations on all of the normalized variables, which included the bands and spectral indices before building regression models. The results implied that mathematical transformations upon the normalized variables that integrated bands and spectral indices should be carried out before building regression models, which improved the retrieval accuracy of soil salinity and sodicity than before.
For better comparing accuracy of the inversion models after different forms of mathematical transformation, all of the variables were normalized, as follows [39]:
x = ( x x m i n ) / ( x m a x x m i n )
where x is an original value, x is the normalized value.
The variables transformed into “1/log(r)” or “er” exhibited better correlations with soil EC (Table 4). For retrieval of soil pH, the “1/r” was best transform (Table 5). At the same time, we found that some nonlinear transformations can enhance the correlations between the variables and EC or pH, and the accuracy of the regression models that were established by combining bands and spectral indices was improved than that by bands or spectral indices per se. The prediction results showed that the R2 for EC increased from 0.704 to 0.735, and from 0.666 to 0.694 for pH. The optimized models were highlighted in Table 4 and Table 5. For the EC regression model, the most influential variables were SI3 and the band Green, which were transformed through “er”. For the pH regression model, the most influential variables were SAVI and the band Coastal, which were transformed through “1/r”. Obviously, the optimized variables and the mathematical transformations were pivotal in well predicting soil EC and pH.
The prediction accuracy was assessed with 47 validating samples using RMSE and R2, which revealed the RMSE of the pH reaching 0.427 (R2 = 0.747) and that of the EC reaching 0.532 mS/cm (R2 = 0.698), respectively. These indicated that the two models were acceptable for estimating soil salinity and sodicity.

3.3. The Spatial Distributions of Soil Salinity and Sodicity

For demonstration purpose, the retrieval results of soil EC and pH were categorized into four soil salinity levels (EC) and four sodicity levels (pH). In this research, EC values were relatively lower than some previous research [1,23,40]. It was likely due to the precipitation prior to sampling time because rainfall had a certain effect on the salinity of surface soils (0~5 cm) [41]. Thus, the salinity levels were ranked as non-salt affected (EC < 0.2 mS/cm), slightly salted (0.2 mS/cm < EC < 0.4 mS/cm), moderately salted (0.4 mS/cm < EC < 0.8 mS/cm), and intensively salted (EC > 0.8 mS/cm), following Taylor [40] and Zhang et al. [23]; and, the sodicity levels were non-sodium affected (pH < 8.5), slightly sodic (8.5 < pH < 9.0), moderately sodic (9.0 < pH < 9.5), and intensively sodic (pH > 9.5), according to Ali [1] and Farifteh et al. [42].
The mapping results showed that the non-salt affecting area was located in northeast and southwest of the study area, and the intensively salted area was in the middle, whilst the moderately salted area was in the south and middle of the area (Figure 2). The salt-affecting area was 1.61 × 106 hm2, accounting for 34.28% of the study area; slightly salted was common in the salt-affecting lands (1.17 × 106 hm2), accounting for 72.90% of the area. The areas of moderately salted and intensively salted were 0.23 × 106 hm2 and 0.21 × 106 hm2, accounting for 14.31% and 12.79% of the salt-affecting area, respectively.
The sodium-affecting area was 1.46 × 106 hm2, accounting for 31.15% of the study area; slightly sodic was the dominant level (0.79 × 106 hm2), accounting for 54.51% of the area. The areas of moderately sodic and intensively sodic were 0.41 × 106 hm2 and 0.25 × 106 hm2, accounting for 28.14% and 17.35% in the sodium-affecting area, respectively. The intensively sodic area was distributed in the middle and the north of the study area, and moderately and slightly sodic areas were in the south and middle. By contrast, the non-sodium affecting areas were close to the non-salt affecting area, which was mainly located in the northeast and southwest of the study area (Figure 3).
Spatial distribution of the soils that were affected by HSS was obtained by overlay analysis with the maps of soil salinity and sodicity. The HSS levels are shown in Table 6. The spatial distribution of HSS indicated that the occurrence of HSS appeared roughly at the same place of soil salinity and soil sodicity (Figure 4). The HSS affecting area was 1.36 × 106 hm2, accounting for 28.91% of the study area. The areas of moderate and intensive HSS were 0.24 × 106 hm2 and 0.19 × 106 hm2; accounting for 18.07% and 14.14% in the HSS affecting area, respectively. The intensive HSS area was distributed in the middle of the study area, as revealed manifest spatial overlap between the salinity and sodicity occurrences.

3.4. Land Cover Classifications

The spatial distribution of land cover in the west Jilin Province is shown in Figure 5. The results indicated that cropland was the largest land cover, with an area of 3.05 × 106 hm2, accounting for 65.00% of the study area, as followed by grassland (12.93%), woodland (5.03%), water body (4.77%), wetland (4.64%), built-up land (4.34%), and barren land (3.29%). Grassland and barren land were mainly distributed in the middle of the study area, where many small lakes and ponds were available; woodland scattered across the area, with many as shelter belts on the fields. Accordingly, these classification results can be a basis for assessing soil salinity and sodicity in relation to land cover in the study area.

3.5. The Distribution of Soil Salinity/Sodicity with Land Cover

In this research, we merely analyzed soil salinity, sodicity, and HSS for woodland, cropland, grassland, and barren land in the study area. Figure 6 illustrates the areas of soil salinity, sodicity, and HSS with regard to land cover. The results indicated that cropland was the largest in area among all of the land cover types, with the areas of soil salinity, sodicity, and HSS to be 0.91 × 106 hm2, 0.76 × 106 hm2, and 0.69 × 106 hm2, respectively. Grassland, barren land, and woodland followed then. However, the slightly salted level resided the prevalent degree in the affected area among cropland, grassland, and woodland, which emerged the same for the slightly HSS level. Further, the barren land was mostly intensively salted and/or sodic (with about 60.00% of the barren land). As for grassland, moderately sodic level was dominant in the sodium-affecting area (32.24%), whereas woodland exhibited the least affected by salinity and sodicity.

3.6. The Distribution of Soil Salinity/Sodicity with Soil Types

Figure 7 presents the areas of soil salinity, sodicity and HSS with regard to soil types. In this research, soil types encompassed eight classes in the west Jilin Province. Chernozem soil and Meadow soil were dominant in this area. The areas of soil salinity, sodicity, and HSS, with regard to soil types, followed as: Chernozem > Meadow soil > Aeolian soil > Solonetz > Chestnut soil > Bog soil > Solonchak > Black soil. The results stipulated that Chernozem was the largest in area among the soil types, in which the areas of soil salinity, sodicity, and HSS were 0.76 × 106 hm2, 0.71 × 106 hm2, and 0.65 × 106 hm2, respectively. Meadow soil and Aeolian soil followed as the second and the third. Solonetz soil was the most seriously affected by salinity, sodicity and HSS, in which the influenced percentages in area were 72.14%, 68.34%, and 66.27%, respectively. Solonchak soil then followed, with an affected percentage of 57.72% in area. In addition, the soils with affected percentages above 35.00% in area included Aeolian soil, Meadow soil, and Bog soil. The soils with affected percentages below 35.00% in area pointed to Chestnut soil, Chernozem, and Black soil (Mollisol). Naturally, the Black soils were least affected among the soils.

4. Discussion

This research endeavored to map soil salinity and sodicity via building the relationships between Landsat 8 OLI imagery and in-situ measurements of soil EC and pH over the west Jilin of China. We pursued establishing and improving the models for soil salinity and sodicity retrieval using PLSR and mathematical transformation of selected variables. As a result, we corroborated that combining Landsat 8 OLI imagery and improved PLSR algorithm can efficiently optimize the retrieval models of soil salinity (measured with soil EC) and sodicity (measured with soil pH).
In the past decades, soil salinity and/or sodicity retrieval models have been proposed often by relating soil physi-chemical properties and spectral performance with RS data and ground measurements [29,43]. By using three dielectric retrieval algorithms, small perturbation, optical physics, and Dubois models to implement analysis, Bell et al. [14] retrieved an improved estimation of a dielectric constant for soil salinity discrimination (optimized R2 = 0.649). Fan et al. [16] established a soil salinity retrieval model (R2 = 0.749) with advanced Multi-Spectral Sensor and PLSR. Janik et al. [24] documented a method by combining partial least square (PLS) and ANN (R2 = 0.94). Nawar et al. [36] suggested the multivariate adaptive regression splines (MARS) model (R2 = 0.73) for soil salinity prediction, which performed better than PLSR. Sidike et al. [3] found that PLSR was better than Stepwise Multiple Regressions (SMR) in retrieving soil salinity. However, we found that the spectral reflectance increased with soil EC, which became dramatic when EC exceeded 1 mS/cm; while pH exceeded 8, the spectral reflectance increased dramatically too (Figure 8). This implied nonlinear relationship might exist between soil EC, pH and the spectral data. Therefore, unlike previous research, we improved the PLSR models’ accuracy via combining reflectance bands and spectral indices through nonlinear mathematical transformation before modelling. Consequently, we confirmed that mathematical transformation indeed could improve the correlation between spectral reflectivity and surface sampling measurements, which further improved the accuracy of the PLSR models. Here, “er”and “1/r” were the optimal mathematical transforms for building the inversion models of soil EC and pH, respectively.
Then, we verified that Landsat 8 OLI imagery was appropriate and economic for extracting the information of soil salinity and sodicity over a large area. Generally, the advantages of RS technology involve quick access and saving time, wider coverage, and constant time series provision for long term monitoring [6,16,44]. By contrast, OLI imagery owns broader views, comparable spatial resolution relative to satellite hyperspectral data (e.g., Hyperion [45]), and abundant band/spectral information competent with higher resolution imagery like QuickBird, IKONOS, and GeoEye-1 [46]. Most importantly, Landsat 8 OLI is free to obtain [47].
Further, we witnessed that linking Landsat OLI-derived variables and field measurements of salinity and sodicity indeed made sense for the mapping and assessment in our case. We found that the band Green best correlated with soil EC, and the band Red best correlated with soil pH (Table 2). Yet, Bai et al. [17] found the band Coastal was best correlated with soil EC and pH. This may lie in the imagery acquisition time difference between spring season in their case and summer growing season in ours. Undoubtedly, it should be better to represent soil EC by reading vegetation performance during summer, as can be readily described by the band Green [29]. As for spectral indices, we found the vegetation indices revealed negatively, and the salinity indices were positively correlated with soil EC and pH. Moreover, the vegetation indices were more sensitive to soil pH than to soil EC, NDVI/SAVI was best correlated spectral index, and SI3 was the most sensitive salinity index to soil EC (Table 2). Nonetheless, Shrestha [11] found that NDVI was poorly negatively correlated with soil EC, and NIR was best correlated with soil EC.
Finally, our mapping results indicated that about one third of the area was affected by salt and sodium in the west Jilin Province. The salt-affecting area was larger than that by sodicity, however, the severity of soil salinity was less than that of sodicity, because the moderate and intensive sodicity accounted for 45.49% of the sodic area, and the moderate and intensive salinity accounted for 29.03% of the saline area. Bai et al. [17] have ever estimated the soil sodicity and salinity for the northern Songnen Plains. But soil sodicity in our study was much more serious (mean pH = 8.70) than that in Bai’s study area (mean pH = 7.92). The causes may lie in less precipitation and higher temperatures in the west Jilin [26]. The soil EC values were lower when compared with other research, which may be due to the regional variations of the study areas [12,15,16,48] and field sampling uncertainties [15,38]. Despite that the soil samples were collected during a desalting period, we classified the levels of soil salinity and sodicity by expert knowledge and characteristics of the study area [14,23]. This ensures that our results can represent the spatial distributions of soil salinity and sodicity in the region well. Cropland (68.51%) was the major contributor to slight salinity, while in the moderate salinity area, grassland was dominant (50.23%). Among all of the land cover types, grassland and barren land were the major contributors to intensive salinity, with percentages 41.48% and 44.31% in area. Cropland and grassland were also dominant contributors to slight and moderate sodicity in the area. However, some differences existed between soil salinity and sodicity, e.g., the contribution rate of cropland to intensive sodicity was twice of that to intensive salinity in area.
However, several aspects should be consolidated in our research. First, we did not sample the soils from built-up land, wetland, and waters for the sake of targeting terrestrial off-town salinity issues. However, the excluded areas were also likely affected by salt problems. Second, human activities played critical roles in affecting the soil salinity and sodicity, e.g., worsening salinity and sodicity of the top soils by bringing salt to soil surface through altering natural water cycles, or allowing for excessive recharge of groundwater and salt accumulation by concentrating [1]. Built-up land is the central of human activities, however, soil salinity and sodicity in built-up land was poorly assessed [49]. Khaledian et al. [50] found that the soil quality of built-up land was poorer than that of other land types. Also, as wetland and waters are important natural resources, their salinity and sodicity will directly affect vegetation and water quality [48,51]. Therefore, to accurately analyze soil degradation by salinity and sodicity, sampling sites should be much representative, better covering all sorts of landscapes and human activities. In future research, soil samples from wetlands and urban areas are needed. Moreover, it is promising to retrospect soil salinity and sodicity during historical stages in accordance with the employed inversion models, and to unveil the areal trending of salt affection through inversion results of soil salinity and sodicity throughout multiple periods.
Land degradation, featured with land salinity mapping in our case, can be interpreted as a decline in land quality and quantity with spatio-temporal variations. However, the soils and RS data used in this study were for only one year. Hence, the use of richer time series of soils and RS data would help to boost an understanding the land degradation processes. We may focus more on the spatio-temporal dynamics of soil sodicity and salinity at diverse scales and the underlying drivers in our future investigations.

5. Conclusions

We successfully mapped multiple forms of soil salinity and sodicity in the west Jilin Province of China using Landsat 8 OLI imagery and improved PLSR algorithm through the nonlinear transformation to retrieve soil salinity and sodicity. Relationships between field measurements of salinity and sodicity, reflectance bands, and spectral indices were built up; and, spatial distributions and occurring severity of salinity and sodicity were also analyzed. The results showed that the most influential variables were SI3 and the band Green were transformed through “er”. For the pH regression model, the most influential variables were SAVI and the band Coastal were transformed through “1/r”. RMSE of pH reached 0.427 (R2 = 0.747) and that of EC reached 0.532 mS/cm (R2 = 0.698).
By statistics from the mapping, about one-third of the study area was influenced by salinity and sodicity. The intensive salinity, sodicity, and HSS lands were mostly distributed in the low-elevation areas near lakes. Soil degradation due to salinity and sodicity was common with regard to land cover in the area. Barren land was the most affected land cover that was subjected to salinity and sodicity. Grassland was the most notable vegetated land cover affected by soil salinity and sodicity, in which exceeding 75.00% in area was salted. By contrast, cropland was the largest affected area among all the land cover, in which the areas that were subjected to salinity, sodicity and HSS were 0.91 × 106 hm2, 0.76 × 106 hm2, and 0.69 × 106 hm2, respectively. With regard to soil types, Chernozem had the largest area subjected to salinity and sodicity, while Solonetz was the most serious type that was affected by salinity. In short, our results may provide realistic guidance for land degradation monitoring, ecological reclamation, salt-soil agriculture management, and comprehensive utilization of natural resources in the west Jilin Province of China.

Acknowledgments

This study was supported by the Strategic Priority Research Program of the Chinese Academy of Sciences (XDA19040500), China Ministry of S&T Special Funds for S&T Infrastructure Establishment (2015FY110500), the National Natural Science Foundation of China (No. 41671219) and the Outstanding Young Scientist Foundation of Northeast Institute of Geography and Agroecology (IGA), Chinese Academy of Sciences (No. Y5H1061001). We are very grateful to those who participated in the field surveys.

Author Contributions

Wang, Hu and Zhang conceived and designed the experiments; Yu and Liu performed the experiments; Yu and Liu analyzed the data; Zhang, Du contributed reagents/ materials/analysis tools; Yu and Liu wrote the paper; Yu prepared all figures and tables. Wang and Hu drafted the paper and revised it critically. All authors reviewed the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ali, M.H. Practices of Irrigation & on-Farm Water Management; Ali, M.H., Ed.; Springer Science + Business Media, LLC: New York, NY, USA, 2011; Volume 2, pp. 271–325. [Google Scholar]
  2. Bouaziz, M.; Matschullat, J.; Gloaguen, R. Improved remote sensing detection of soil salinity from a semi-arid climate in northeast Brazil. C. R. Geosci. 2011, 343, 795–803. [Google Scholar] [CrossRef]
  3. Sidike, A.; Zhao, S.; Wen, Y. Estimating soil salinity in Pingluo county of China using Quickbird data and soil reflectance spectra. Int. J. Appl. Earth Obs. Geoinf. 2014, 26, 156–175. [Google Scholar] [CrossRef]
  4. Liu, X.T. Discussion on ecological conservation strategy in western Songnen Plain. Syst. Sci. Compr. Stud. Agric. 2003, 19, 282–289. (In Chinese) [Google Scholar]
  5. Zhang, G.X.; Yang, J.F.; Liu, Q. Discussion on the problems of water environment ecologization in ecological housing district. Ecol. Environ. 2004, 13, 290–292. (In Chinese) [Google Scholar]
  6. Setia, R.; Lewis, M.; Marschner, P.; Raja Segaran, R.; Summers, D.; Chittleborough, D. Severity of salinity accurately detected and classified on a paddock scale with high resolution multispectral satellite imagery. Land Degrad. Dev. 2013, 24, 375–384. [Google Scholar] [CrossRef]
  7. Thimmappa, K.; Sharma, D.K.; Dagar, J.C.; Raju, R. Reclamation of salt-affected soils: Socioeconomic impact assessment. In Innovative Saline Agriculture; Dagar, J.C., Sharma, P.C., Sharma, D.K., Singh, A.K., Eds.; Springer Nature: New Delhi, India, 2016; pp. 489–505. [Google Scholar]
  8. Metternicht, G.I.; Zinck, J.A. Remote sensing of soil salinity: Potentials and constraints. Remote Sens. Environ. 2003, 85, 1–20. [Google Scholar] [CrossRef]
  9. Periasamy, S.; Shanmugam, R.S. Multispectral and microwave remote sensing models to survey soil moisture and salinity. Land Degrad. Dev. 2017, 28, 1412–1425. [Google Scholar] [CrossRef]
  10. Wu, J.; Vincent, B.; Yang, J.; Bouarfa, S.; Vidal, A. Remote sensing monitoring of changes in soil salinity: A case study in Inner Mongolia, China. Sensors 2008, 8, 7035–7049. [Google Scholar] [CrossRef] [PubMed]
  11. Shrestha, R.P. Relating soil electrical conductivity to remote sensing and other soil properties for assessing soil salinity in northeast Thailand. Land Degrad. Dev. 2006, 17, 677–689. [Google Scholar] [CrossRef]
  12. Taghizadeh-Mehrjardi, R.; Minasny, B.; Sarmadian, F.; Malone, B.P. Digital mapping of soil salinity in Ardakan region, central Iran. Geoderma 2014, 213, 15–28. [Google Scholar] [CrossRef]
  13. Yu, R.; Liu, T.; Xu, Y.; Zhu, C.; Zhang, Q.; Qu, Z.; Liu, X.; Li, C. Analysis of salinization dynamics by remote sensing in Hetao Irrigation District of North China. Agric. Water Manag. 2010, 97, 1952–1960. [Google Scholar] [CrossRef]
  14. Bell, D.; Menges, C.; Ahmad, W.; Zyl, J.J.V. The application of dielectric retrieval algorithms for maping soil salinity in tropical coastal environment using airborne polarimetric SAR. Remote Sens. Environ. 2001, 75, 375–384. [Google Scholar] [CrossRef]
  15. Allbed, A.; Kumar, L.; Aldakheel, Y.Y. Assessing soil salinity using soil salinity and vegetation indices derived from IKONS high-spatial resolution imageries: Applications in a date palm dominated region. Geoderma 2014, 230–231, 1–8. [Google Scholar] [CrossRef]
  16. Fan, X.; Liu, Y.; Tao, J.; Weng, Y. Soil salinity retrieval from advanced multi-spectral sensor with partial least square regression. Remote Sens. 2015, 7, 488–511. [Google Scholar] [CrossRef]
  17. Bai, L.; Wang, C.Z.; Zang, S.Y.; Zhang, Y.H.; Hao, Q.N.; Wu, Y.X. Remote sensing of soil alkalinity and salinity in the Wuyu’er-shuangyang River Basin, Northeast China. Remote Sens. 2016, 8, 163. [Google Scholar] [CrossRef]
  18. Bromham, L.; Saslis-Lagoudakis, C.H.; Bennett, T.H.; Flowers, T.J. Soil alkalinity and salt tolerance: Adapting to multiple stresses. Biol. Lett. 2013, 9, 20130566. [Google Scholar] [CrossRef] [PubMed]
  19. Bui, E. Possible role of soil alkalinity in plant breeding for salt-tolerance. Biol. Lett. 2013, 9, 20130566. [Google Scholar] [CrossRef] [PubMed]
  20. Yan, B.; Fang, N.F.; Zhang, P.C.; Shi, Z.H. Impacts of land use change on watershed streamflow and sediment yield: An assessment using hydrologic modelling and partial least squares regression. J. Hydrol. 2013, 484, 26–37. [Google Scholar] [CrossRef]
  21. Sawut, M.; Ghulam, A.; Tiyip, T.; Zhang, Y.-j.; Ding, J.-l.; Zhang, F.; Maimaitiyiming, M. Estimating soil sand content using thermal infrared spectra in arid lands. Int. J. Appl. Earth Obs. Geoinf. 2014, 33, 203–210. [Google Scholar] [CrossRef]
  22. Udelhoven, T.; Emmerling, C.; Jarmer, T. Quantitative analysis of soil chemical properties with diffuse reflectance spectrometry and partial least-square regression: A feasibility study. Plant Soil 2003, 319–329. [Google Scholar] [CrossRef]
  23. Zhang, T.T.; Zeng, S.L.; Gao, Y.; Ouyang, Z.T.; Li, B.; Fang, C.M.; Zhao, B. Using hyperspectral vegetation indices as a proxy to monitor soil salinity. Ecol. Indic. 2011, 11, 1552–1562. [Google Scholar] [CrossRef]
  24. Janik, L.J.; Forrester, S.T.; Rawson, A. The prediction of soil chemical and physical properties from mid-infrared spectroscopy and combined partial least-squares regression and neural networks (PLS-NN) analysis. Chemom. Intell. Lab. Syst. 2009, 97, 179–188. [Google Scholar] [CrossRef]
  25. Pang, G.; Wang, T.; Liao, J.; Li, S. Quantitative model based on field-derived spectral characteristics to estimate soil salinity in Minqin county, China. Soil Sci. Soc. Am. J. 2014, 78, 546–554. [Google Scholar] [CrossRef]
  26. Fang, L.; Bai, Z.; Wei, S.; Yanfen, H.; Zongming, W.; Kaishan, S.; Dianwei, L.; Zhiming, L. Sandy desertification change and its driving forces in western Jilin Province, North China. Environ. Monit. Assess. 2008, 136, 379–390. [Google Scholar] [CrossRef] [PubMed]
  27. Wang, Z.; Huang, N.; Luo, L.; Li, X.; Ren, C.; Song, K.; Chen, J.M. Shrinkage and fragmentation of marshes in the west Songnen Plain, China, from 1954 to 2008 and its possible causes. Int. J. Appl. Earth Obs. Geoinf. 2011, 13, 477–486. [Google Scholar] [CrossRef]
  28. Wu, W.; Mhaimeed, A.S.; Al-Shafie, W.M.; Ziadat, F.; Dhehibi, B.; Nangia, V.; De Pauw, E. Mapping soil salinity changes using remote sensing in central Iraq. Geoderma Reg. 2014, 2–3, 21–31. [Google Scholar] [CrossRef]
  29. Douaoui, A.E.K.; Nicolas, H.; Walter, C. Detecting salinity hazards within a semiarid context by means of combining soil and remote-sensing data. Geoderma 2006, 134, 217–230. [Google Scholar] [CrossRef]
  30. Huete, A.R.; Liu, H.Q.; Batchily, K.; van Leeuwen, W. A comparison of vegetation indices global set of TM images for EOS-MODIS. Remote Sens. Environ. 1997, 59, 440–451. [Google Scholar] [CrossRef]
  31. Alhammadi, M.S.; Glenn, E.P. Detecting date palm trees health and vegetation greenness change on the eastern coast of the United Arab Emirates using SAVI. Int. J. Remote Sens. 2008, 29, 1745–1765. [Google Scholar] [CrossRef]
  32. Rodriguez-Galiano, V.F.; Chica-Olmo, M.; Abarca-Hernandez, F.; Atkinson, P.M.; Jeganathan, C. Random forest classification of mediterranean land cover using multi-seasonal imagery and multi-seasonal texture. Remote Sens. Environ. 2012, 121, 93–107. [Google Scholar] [CrossRef]
  33. Guo, L.; Chehata, N.; Mallet, C.; Boukir, S. Relevance of airborne lidar and multispectral image data for urban scene classification using Random Forests. ISPRS J. Photogramm. 2011, 66, 56–66. [Google Scholar] [CrossRef]
  34. Ouyang, Z.; Zhang, L.; Wu, B.; Li, X.; Xiao, Y.; Zheng, H. An ecosysytem classification system based on reomte sensor information in China. Acta Ecol. Sin. 2015, 35, 219–226. (In Chinese) [Google Scholar]
  35. Zhang, L.; Wu, B.; Li, X.; Xing, Q. Classification system of China land cover for carbon budget. Acta Ecol. Sin. 2014, 34, 7158–7166. (In Chinese) [Google Scholar] [CrossRef]
  36. Nawar, S.; Buddenbaum, H.; Hill, J.; Kozak, J. Modeling and mapping of soil salinity with reflectance spectroscopy and landsat data using two quantitative methods (PLSR and ANN). Remote Sens. 2014, 6, 10813–10834. [Google Scholar] [CrossRef]
  37. Farifteh, J.; Van der Meer, F.; Atzberger, C.; Carranza, E.J.M. Quantitative analysis of salt-affected soil reflectance spectra: A comparison of two adaptive methods (PLSR and ANN). Remote Sens. Environ. 2007, 110, 59–78. [Google Scholar] [CrossRef]
  38. Weng, Y.; Gong, P.; Zhu, Z. A spectral index for estimating soil salinity in the yellow river delta region of China using EO-1 hyperion data. Pedosphere 2010, 20, 378–388. [Google Scholar] [CrossRef]
  39. Mirlas, V.; Benyamini, Y.; Marish, S.; Gotesman, M.; Fizik, E.; Agassi, M. Method for normalization of soil salinity data. J. Irrig. Drain. Eng. ASCE 2003, 129, 64–66. [Google Scholar] [CrossRef]
  40. Taylor, S. Dryland Salinity-Introductory Extension Notes, 2nd ed.; Department of Conservation and Land Management: Perth, Australia, 1993. [Google Scholar]
  41. Fang, S.; Chen, X. Influence of atmospheric precipitation on soil leaching and desalinization in the North China plain. Acta Pedol. Sin. 2005, 42, 730–736. (In Chinese) [Google Scholar]
  42. Farifteh, J.; van der Meer, F.; van der Meijde, M.; Atzberger, C. Spectral characteristics of salt-affected soils: A laboratory experiment. Geoderma 2008, 145, 196–206. [Google Scholar] [CrossRef]
  43. Elhag, M. Evaluation of different soil salinity mapping using remote sensing techniques in arid ecosystems, Saudi Arabia. J. Sens. 2016, 2016, 1–8. [Google Scholar] [CrossRef]
  44. Navarro-Pedreño, J.; Jordan, M.M.; Meléndez-Pastor, I.; Gómez, I.; Juan, P.; Mateu, J. Estimation of soil salinity in semi-arid land using a geostatistical model. Land Degrad. Dev. 2007, 18, 339–353. [Google Scholar] [CrossRef]
  45. Ghosh, G.; Kumar, S.; Saha, S.K. Hyperspectral satellite data in mapping salt-affected soils using linear spectral unmixing analysis. J. Indian Soc. Remote Sens. 2012, 40, 129–136. [Google Scholar] [CrossRef]
  46. De Benedetto, D.; Castrignano, A.; Diacono, M.; Rinaldi, M.; Ruggieri, S.; Tamborrino, R. Field partition by proximal and remote sensing data fusion. Biosyst. Eng. 2013, 114, 372–383. [Google Scholar] [CrossRef]
  47. Aldabaa, A.A.A.; Weindorf, D.C.; Chakraborty, S.; Sharma, A.; Li, B. Combination of proximal and remote sensing methods for rapid soil salinity quantification. Geoderma 2015, 239–240, 34–46. [Google Scholar] [CrossRef]
  48. Cemek, B.; Guler, M.; Kilic, K.; Demir, 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]
  49. Zornoza, R.; Acosta, J.A.; Bastida, F.; Domínguez, S.G.; Toledo, D.M.; Faz, A. Identification of sensitive indicators to assess the interrelationship between soil quality, management practices and human health. Soil 2015, 1, 173–185. [Google Scholar] [CrossRef]
  50. Khaledian, Y.; Kiani, F.; Ebrahimi, S.; Brevik, E.C.; Aitkenhead-Peterson, J. Assessment and monitoring of soil degradation during land use change using multivariate analysis. Land Degrad. Dev. 2017, 28, 128–141. [Google Scholar] [CrossRef]
  51. An, D.; Zhao, G.; Chang, C.; Wang, Z.; Li, P.; Zhang, T.; Jia, J. Hyperspectral field estimation and remote-sensing inversion of salt content in coastal saline soils of the Yellow River Delta. Int. J. Remote Sens. 2016, 37, 455–470. [Google Scholar] [CrossRef]
Figure 1. Location of the study area, digital elevation model (DEM) and soil sampling locations.
Figure 1. Location of the study area, digital elevation model (DEM) and soil sampling locations.
Sensors 18 01048 g001
Figure 2. Distributions of soil salinity in the west Jilin.
Figure 2. Distributions of soil salinity in the west Jilin.
Sensors 18 01048 g002
Figure 3. Distributions of soil sodicity in the west Jilin.
Figure 3. Distributions of soil sodicity in the west Jilin.
Sensors 18 01048 g003
Figure 4. Distributions of hybridized salinity and sodicity (HSS) in the west Jilin.
Figure 4. Distributions of hybridized salinity and sodicity (HSS) in the west Jilin.
Sensors 18 01048 g004
Figure 5. Spatial distribution of land cover in the west Jilin.
Figure 5. Spatial distribution of land cover in the west Jilin.
Sensors 18 01048 g005
Figure 6. The areas of soil salinity, sodicity, and HSS at different levels with regard to land cover (a) Salinity (b) Sodicity (c) HSS (d) Affected soils.
Figure 6. The areas of soil salinity, sodicity, and HSS at different levels with regard to land cover (a) Salinity (b) Sodicity (c) HSS (d) Affected soils.
Sensors 18 01048 g006
Figure 7. The areas affected by soil salinity, sodicity, and HSS at different levels with regard to soil types. (a) Salinity (b) Sodicity (c) HSS (d) Affected soil types (CN: Chestnut soil, MS: Meadow soil, AS: Aeolian soil, CZ: Chernozem, BL: Black soil ST: Solonetz, BO: Bog soil SC: Solonchak).
Figure 7. The areas affected by soil salinity, sodicity, and HSS at different levels with regard to soil types. (a) Salinity (b) Sodicity (c) HSS (d) Affected soil types (CN: Chestnut soil, MS: Meadow soil, AS: Aeolian soil, CZ: Chernozem, BL: Black soil ST: Solonetz, BO: Bog soil SC: Solonchak).
Sensors 18 01048 g007
Figure 8. Correlations between optimal variables and soil EC and pH measurements.
Figure 8. Correlations between optimal variables and soil EC and pH measurements.
Sensors 18 01048 g008
Table 1. The expressions of spectral indices.
Table 1. The expressions of spectral indices.
Spectral IndexExpressionFull NameReference
NDVI(NIR − Red)/(NIR + Red)Normalized Differential Vegetation IndexShrestha (2006)
DVINIR − RedDifferential Vegetation IndexClevers et al. (1988)
EVINIR/RedEnhanced Vegetation IndexHuete et al. (1997)
SAVI(NIR − Red)/(NIR+Red + 0.5)   ×   1.5Soil Adjusted Vegetation IndexBouaziz et al. (2011)
SI Green   ×   Red   Salinity IndexBouaziz et al. (2011)
SI2 Green 2 + Red 2 + NIR 2   Salinity Index2Douaoui et al. (2006)
SI3 Green 2 + Red 2   Salinity Index3Douaoui et al. (2006)
SI4SWIR1/NIRSalinity Index4Douaoui et al. (2006)
SRSI ( NDVI 1 ) 2 + SI 2 Salinization Remote Sensing IndexAlhammadi et al. (2008)
Table 2. Correlation coefficients between band reflectance, soil pH and EC.
Table 2. Correlation coefficients between band reflectance, soil pH and EC.
ECpHCoastalRedGreenBlueNIRSWIR1SWIR2PAN
EC1---------
pH0.703 **1--------
Coastal0.821 **0.791 **1-------
Red0.810 **0.805 **0.988 **1------
Green0.826 **0.793 **0.991 **0.992 **1-----
Blue0.818 **0.795 **0.998 **0.992 **0.992 **1----
NIR0.348 **0.068 **0.294 **0.245 *0.334 **0.283 **1---
SWIR10.704 **0.788 **0.912 **0.944 **0.925 **0.923 **0.160 *1--
SWIR20.738 **0.791 **0.939 **0.967 **0.948 **0.948 **0.143 *0.989 **1-
PAN0.763 **0.773 **0.958 **0.964 **0.966 **0.963 **0.289 **0.906 **0.931 **1
Note: * represents significance at the 0.05 level; ** represents significance at the 0.01 level; bold number means the maximum value in a column.
Table 3. Correlation coefficients between spectral indices, soil pH and EC.
Table 3. Correlation coefficients between spectral indices, soil pH and EC.
ECpHSISI2SI3SI4SRSISAVINDVIEVIDVI
EC1----------
pH0.703 **1---------
SI0.818 **0.803 **1--------
SI20.722 **0.507 **0.770 **1-------
SI30.826 **0.793 **1.000 **0.804 **1------
SI40.472 **0.704 **0.731 **0.209 *0.695 **1-----
SRSI0.818 **0.803 **1.000 **0.771 **0.997 **0.731 **1----
SAVI−0.665 **−0.814 **−0.931 **−0.467 **−0.888 **−0.920 **−0.912 **1---
NDVI−0.665 **−0.814 **−0.931 **−0.467 **−0.888 **−0.920 **−0.912 **1.000 **1--
EVI−0.449 **−0.694 **−0.708 **−0.159 *−0.677 **−0.877 **−0.708 **0.889 **0.889 **1-
DVI−0.555 **−0.725 **−0.775 **−0.198 *−0.736 **−0.912 **−0.775 **0.935 **0.935 **0.909 **1
Note: * represents significance at the 0.05 level; ** represents significance at the 0.01 level; bold number means the maximum value in a column.
Table 4. The regression models for soil EC in relation to bands, spectral indices and hybrid variables.
Table 4. The regression models for soil EC in relation to bands, spectral indices and hybrid variables.
TransformImaging FeatureR2ExpressionSDRMSE
rBand0.706EC = 1.601Green − 3.093SWIR1 − 1.1210.4560.451
Index0.704EC = 21.114NDVI + 5.411SI3 − 12.0240.4570.453
Hybrid0.706EC = 5.854SI3 − 3.093SWIR1 − 0.0290.4550.451
1/rBand0.646EC = − 45.758Blue − 1.031NIR + 2.6700.5010.527
Index0.707EC = 8.506EVI − 11.121NDVI + 19.7170.4540.451
Hybrid0.712EC = 7.915EVI + 0.161Red − 2.2360.450.446
erBand0.717EC = 8.553Green − 1.969SWIR1 − 7.5230.4460.442
Index0.711EC = 2.783SI3 − 3.1650.450.448
Hybrid0.735EC = 14.090SI3 − 24.833Green + 13.1740.4320.429
log (r)Band0.635EC = 9.244Green − 3.032Red + 4.5570.5070.503
Index0.601EC = 1.692SI3 + 2.764SI2 + 2.5490.530.525
Hybrid0.709EC = 20.220Green − 9.125SI + 8.0660.4520.449
1/log(r)Band0.726EC = − 2.052Green + 0.363SWIR1 − 1.9690.4320.435
Index0.727EC = − 0.486SI3 − 0.4500.4370.435
Hybrid0.727EC = − 0.486SI3 − 0.4500.4370.435
1/erBand0.691EC = 4.634SWIR1 − 15.431Green + 9.4780.4670.462
Index0.691EC = − 6.391SI3 − 7.701EVI + 3.169DVI + 8.6180.4690.463
Hybrid0.734EC = − 90.819Green + 43.953SI3 + 38.9160.4320.429
Sqrt(r)Band0.675EC = 11.977Green − 3.844SWIR1 − 2.9510.4790.474
Index0.682EC = 5.083SI3 + 7.866EVI − 3.338DVI − 4.7270.4740.469
Hybrid0.724EC = 58.351Green − 27.823SI3 − 13.2140.4410.437
Note: r represents the normalized variables (bands and spectral indices).
Table 5. The regression models for soil pH in relation to bands, spectral indices, and hybrid parameters.
Table 5. The regression models for soil pH in relation to bands, spectral indices, and hybrid parameters.
TransformImaging FeatureR2ExpressionSD RMSE
rBand0.666pH = 6.082Red − 1.743NIR + 8.1560.5730.568
Index0.663pH = − 62.762NDVI + 41.7010.5750.571
Hybrid0.686pH = 3.753Green − 40.192NDVI + 28.9700.5570.551
1/rBand0.689pH = − 63.919Blue − 0.499NIR + 43.3300.5540.548
Index0.691pH = 3.541SAVI − 0.346Red − 18.8630.5520.547
Hybrid0.694pH = 1.705SAVI − 1.348Coastal − 1.8980.5510.546
erBand0.651pH = 4.442Red − 1.306NIR + 5.2000.5860.581
Index0.682pH = − 155.862SAVI + 1.049SI3 + 183.8620.560.56
Hybrid0.684pH = − 147.848SAVI + 2.551Green + 173.0260.5580.553
log (r)Band0.689pH = 3.466Blue + 11.9680.5520.548
Index0.682pH = 2.761SI + 10.5820.5580.554
Hybrid0.689pH = 3.466Blue + 11.9680.5520.548
1/log(r)Band0.653pH = − 2.287Blue + 0.229NIR + 6.8660.5840.579
Index0.677pH = 493.206SAVI − 44.705NDVI + 394.7190.5640.559
Hybrid0.691pH = 41.952SAVI − 0.559Green + 54.2730.5600.555
1/erBand0.678pH = − 8.159Red + 2.137NIR + 13.9220.5640.558
Index0.688pH = 129.434SAVI − 3.412SRSI − 102.8210.5540.549
Hybrid0.688pH = 129.434SAVI − 3.412SRSI − 102.8210.5540.549
Sqrt(r)Band0.678pH = 6.139Red + 5.8510.5610.558
Index0.681pH = 4.617SRSI + 6.5290.5580.555
Hybrid0.689pH = 4.718SRSI − 1.647NIR + 7.6430.5530.548
Note: r represents the normalized variables (bands and spectral indices).
Table 6. The standard for soil HSS severity ranking.
Table 6. The standard for soil HSS severity ranking.
LevelRange of EC and pH
UnaffectedpH < 8.5 or EC < 0.2 mS/cm
Slightly affected0.2 mS/cm < EC < 0.4 mS/cm and 8.5 < pH < 9.0;
0.4 mS/cm < EC < 0.8 mS/cm and 8.5 < pH < 9.0;
0.2 mS/cm < EC < 0.4 mS/cm and 9.0 < pH < 9.5;
Moderately affected0.4 mS/cm < EC < 0.8 mS/cm and9.0 < pH < 9.5;
EC > 0.8 mS/cm and 8.5 < pH < 9.0;
0.2 mS/cm < EC < 0.4 mS/cm and pH > 9.5;
Intensively affected0.8 mS/cm < EC and pH > 9.5;
0.4 mS/cm < EC < 0.8 mS/cm and pH > 9.5;
0.8 mS/cm < EC and 9.0 < pH < 9.5;

Share and Cite

MDPI and ACS Style

Yu, H.; Liu, M.; Du, B.; Wang, Z.; Hu, L.; Zhang, B. Mapping Soil Salinity/Sodicity by using Landsat OLI Imagery and PLSR Algorithm over Semiarid West Jilin Province, China. Sensors 2018, 18, 1048. https://doi.org/10.3390/s18041048

AMA Style

Yu H, Liu M, Du B, Wang Z, Hu L, Zhang B. Mapping Soil Salinity/Sodicity by using Landsat OLI Imagery and PLSR Algorithm over Semiarid West Jilin Province, China. Sensors. 2018; 18(4):1048. https://doi.org/10.3390/s18041048

Chicago/Turabian Style

Yu, Hao, Mingyue Liu, Baojia Du, Zongming Wang, Liangjun Hu, and Bai Zhang. 2018. "Mapping Soil Salinity/Sodicity by using Landsat OLI Imagery and PLSR Algorithm over Semiarid West Jilin Province, China" Sensors 18, no. 4: 1048. https://doi.org/10.3390/s18041048

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