The Applicability of Remote Sensing Models of Soil Salinization Based on Feature Space

Soil salinization is a major challenge for the sustainable use of land resources. An optimal remote sensing inversion model could monitor regional soil salinity across diverse geographical areas. In this study, the feature space method was used to study the applicability of the inversion model for typical salt-affected soils in China (Yanqi Basin (arid area) and Kenli County (coastal area)), and to obtain soil salinity grade distribution maps. The salinity index (SI) surface albedo (Albedo)model was the most accurate in both arid and coastal regions with overall accuracy reaching 93.3% and 88.8%, respectively. The sensitivity factors for the inversion of salinity in both regions were the same, indicating that the SI-Albedo model is applicable for monitoring salinity in arid and coastal areas of China. We combined Landsat 8 Operational Land Imager image data and field data to obtain the distribution pattern of soil salinity using the SI-Albedo model and proposed corresponding countermeasures for soil salinity in the Yanqi Basin and Kenli County according to the degree of salinity. This study on soil salinity in arid and coastal areas of China will provide a useful reference for future research on soil salinity both in China and globally.


Introduction
Land degradation is a major global environmental problem [1][2][3]. Soil salinization, an important element of land degradation, significantly impacts the global ecological environment. Soil salinization is disturbed by the natural environment and human activities when capillary action or evapotranspiration leads to the rise of groundwater and the accumulation of soluble salts at the surface, resulting in different degrees of salinization [4]. Soil salinization not only leads to a decrease in soil fertility, acid-base imbalance, and an increasingly poor land environment but also intensifies the conflict between socioeconomic development and environmental protection [5]. Therefore, it is important to extract soil salinization information dynamically, quickly, and accurately on a regional scale, monitor the spatial distribution characteristics of saline soils, and quantify the degree of soil salinization to manage soil salinization and prevent further deterioration of the ecological environment. Managing soil salinization provides security for the development of efficient agriculture, working towards achieving the goal of sustainable development of agricultural production and the social economy. Current salinization research is mostly focused on semi-arid and arid regions, and coastal areas. As a typical area of coastal soil salinization in China [6], the Yellow River Delta has a shallow groundwater depth and uneven spatial distribution of mineralization, complex regional water-salt dynamic processes, and serious re-salinization phenomena. The difficulty in using this moderately to severely salinized soil seriously restricts the development and utilization of land resources and the sustainable socioeconomic development of the Yellow River Delta.
The primary and secondary salinization of soils in arid and semi-arid regions is an increasingly prominent problem. With the large-scale farmland reclamation and the promotion of dropper technology, ecological problems such as reduced irrigation water and increased soil salinization have emerged, which seriously restrict oasis agriculture and socioeconomic development [7]. However, few studies are available on the universality of soil salinity inversion models using feature space for coastal versus inland drylands.
In this study, we used Kenli County of the Yellow River Delta, a coastal area, and the Yanqi Basin, an arid inland area, as the research areas. We will use the feature space method combined with field survey data and Landsat 8 OLI images to quantitatively invert the distribution pattern of salinization. We will then determine the most suitable two-region model by analyzing and comparing the applicability of the three models to different geographic environments and different causes of salinization. Our study hopes to provide a scientific basis for future sustainable development and utilization of land resources, ecological environmental protection, and dynamic monitoring of salinization.
This study answers the following research questions (RQs): RQ 1: Which is the best performing feature space model for the inversion of salinization in the Yanqi basin and the Yellow River delta region? RQ 2: Are the models that are suitable for the dynamic monitoring of salinization in arid and coastal areas consistent? RQ 3: Which factors are the most influential in the development of salinization in the Yellow River Delta and the Yanqi Basin?
RQ4: What are the main recommendations for the management of salinized areas in both regions?
This study includes five more sections. The second section presents literature related to the application of remote sensing technology for salinization monitoring and establishes the objective of the present research. The third section explains the research methodology, including information about the data used, salinization evaluation index, and the feature space method that has been applied to define salinization classification models. Section 4 presents the obtained results. Finally, Sections 5 and 6 focus on explaining the achieved results, comparing them with other studies, highlighting the limitations and main conclusions of this study, and providing suggestions for future lines of work.

Related Work
Numerous studies in the existing literature have sought to build tools whose purpose is to monitor soil salinization dynamically and accurately. variables. The global salinization area is vast [8,9], and the traditional point survey method for monitoring soil salinization is costly and time-consuming. Therefore, monitoring soil salinization using remote sensing is of great scientific and practical significance for achieving global sustainable development [10,11]. Since the 1970s, analysis of soil salinization using remote sensing data from multiple sources has been getting increasing attention. In 2007, Brunner et al. [12] studied the Bosten Lake area in Xinjiang to evaluate soil salinity using the NDVI index and obtained a soil salinity distribution map.In 2013, Madani et al. [13] studied the Siwa oasis in the western desert of Egypt using Landsat data, wherein they converted the Digital Number (DN) values of the image into reflectance and obtained a strong correlation between soil salinity, short-wave, and thermal infrared bands. In 2016, Azabdaftari et al. [14] extracted the Normalized Difference Vegetation Index (NDVI), Ratio Vegetation Index (RVI), Soil-Adjusted Vegetation Index (SAVI), Normalized Difference Snow Index (NDSI), and other indexes using Landsat in Adana, Turkey. In order to create a soil salinity map of the area, field-measured conductivity data were used to analyze the relationship between the original image band values and those of different indices In 2019, Hoa et al. [15] selected the Gaussian Process (GP) through a multi-machine learning model comparison using Sentinel-1A data to map the distribution of soil salinity in the Mekong Delta.
In summary, most scholars have used supervised and unsupervised classification, automatic neural network extraction, and other research methods to study the dynamic monitoring and spatiotemporal characteristics of salinization areas [16,17]. These semiautomatic classification methods are time-consuming, labor-intensive, inefficient, have long decoding periods, and are subject to a high degree of subjective influence; therefore, it is difficult to use them to conduct long-term series salinization dynamic monitoring. Using various types of indicators in multispectral remote sensing data to construct a feature space for extracting and monitoring salinization information, is an advanced method for salinity remote sensing monitoring research. Many scholars are now using indicators closely related to soil salinity to establish feature spaces to extract soil salinization information and have achieved positive results. Wang et al. [18] proposed the NDVI-SI feature space concept by considering both vegetation and soil information; they used the normalized difference vegetation index (NDVI) and salinity index (SI) to quantitatively analyze and monitor salinization in Xinjiang Yutian Oasis. Ha et al. [19] selected SI and surface reflectance (Albedo) to build a feature space and construct a model to extract information on the distribution of soil salinization in the Kriya Oasis. Ding et al. [20] used Landsat -TM data to analyze the relationship between saline soils, the Modified Soil-adjusted Vegetation Index (MSAVI), and Wetness Index (WI). Subsequently, they proposed the MSAVI-WI feature space and a remote sensing monitoring index model for soil salinity.
Spectral feature space is an effective modern method for the inversion of surface parameters using multi-band or multi-indicator synergistic changes with simple and feasible model indicators that are easy to extract. This method is widely used in research on surface evapotranspiration, crop water content, and for monitoring soil moisture, droughts, and desertification [21]. Recently the feature space model method has been applied for the inversion of soil salinization parameters with promising results [22,23].
The above literature review demonstrates the successful application of remote sensing techniques in monitoring salinity-related problems. However, the following research gaps were observed: (1) most studies used a single evaluation index or a particular feature space model for monitoring salinity and did not compare multiple special spatial models; (2) in a previous study, the authors built the optimal inversion model in only one area, the Yellow River Delta, and they did not conduct a model applicability analysis in different geographical locations.
In response to these gaps in the literature, the objectives of this study are as follows: (1) to establish a quantitative model that combines soil spectral and vegetation biophysical characteristics, and soil salinity to determine the parameters suitable for salinity inversion in the two regions, and analyze the general applicability of the feature space model to achieve a quantitative inversion and dynamic monitoring of soil; (2) and to determine the different drivers of salinity in the Yellow River Delta and the Yanqi Basin and propose recommendations for the prevention and control of salinity in key areas for protecting ecosystems and maintaining sustainable agricultural development.

Overview of the Study Area
In this study, two study areas were selected: Kenli County in the Yellow River Delta, a coastal saline area, and the Yanqi Basin, an arid saline area.
The Yellow River Delta, located in northern Shandong Province, is an alluvial plain formed by sediment carried by the Yellow River Estuary in the coastal zone. Kenli County is in the Yellow River Estuary. The coordinate range is between 37 • 24 -38 • 10 N and 118 • 15 -119 • 19 E (Figure 1). According to the Chinese climate zones, it has an area of approximately 233,100 ha and is categorized as a humid zone. The landform type is influenced by the repeated diversion of the Yellow River and sediment deposition into the sea and is thus very complex; the terrain is slightly lowered in a fan shape from the southwest to the northeast [24]. Because of the influence of the continental monsoon, the distribution of precipitation in this area is not uniform and is mostly concentrated in the summer. The average annual precipitation is 530-630 mm, the average annual evaporation is 1750-2430 mm, and the ratio of evaporation to rainfall of approximately 3.5:1. Most of the soil types are fragile saline and tidal soils. Saline soil is primarily distributed in the area with a phreatic water level of 1-2 m and a salinity of 10-30 g/L, and the high content of soluble salts in the soil is not conducive to crop growth. The main crops are wheat, corn, cotton, and rice, with low yields. There is a large distribution of salt wasteland and serious confusion in the soil spectrum.
The Yanqi Basin, in the northeast of the Tarim Basin, is a semi-enclosed mountain basin. The coordinate range is between 41 • 45 -42 • 15 N and 85 • 30 -88 • 15 E (Figure 1). According to the Chinese climate zone, it has an area of approximately 1,300,400 ha and is categorized as an arid zone, with typical arid zone oasis characteristics, such as large daily differences in temperature, low precipitation, high evaporation, and dry air. The northern part of the area is a flooded micro-slope plain, the middle contains Bosten Lake, and the east is flowing sand dunes. The terrain is high in the west, low in the east, high in the north, and low in the south, and the Kaidu River runs from west to east into Bosten Lake. The region is affected by a continental desert climate, with an average annual precipitation of 50-80 mm, average annual evaporation of 2000.5-2449.7 mm, and an evapotranspiration ratio greater than 30:1. Under this specific environment, underground runoff and salinity lack outlets, and soil salinization is more common [25]. The soil types are mainly brown desert, meadow, tidal, saline, and weathered soil [26] and are nutrient deficient, which adversely affects local agricultural production. The main crops are corn, tomato, and pepper.

Remote Sensing Image Data
Landsat 8 Operational Land Imager (OLI) satellite data launched by the National Aeronautics and Space Administration (NASA) was selected as the data source. The imager has nine bands, a coving width of 185 × 185 km, and a resolution of 30 m. The Landsat 8 OLI satellite remote sensing image on 30 April 2020, was selected for the Yanqi Basin, and 29 April 2019, was selected for Kenli County and the Yellow River Delta because these two areas experience abundant evaporation and annual soil moisture lows during springtime, and salt in the soil gathered on the surface. Additionally, farmlands have not been cultivated in the spring, and the soil is mostly covered by sparse ground vegetation.

Field Sampling Data
The data of the measured sample points of the two areas included soil salinity and soil conductivity data. The sampling points of the field-measured verification data were evenly distributed among the farmlands in Kenli County and the Yanqi Basin. The sampling data for Kenli County (47) were acquired on-site by the research team in April 2019, and sampling data for the Yanqi Basin (221) were acquired on-site by the research team in May 2020. Sample points were evenly distributed throughout the Kenli and the Yanqi basin areas ( Figure 2). The coordinates of the sample points were measured by GPS in the field with centimeter-level accuracy. Three soil samples were randomly selected within 25 m 2 of each sampling point, mixed, and placed in a sampling bag for subsequent physical and chemical analyses. The soil sampling depth was 0-20 cm, and approximately 1000 g of soil was taken at each sampling point. The salinity content was averaged from multiple measurements taken at the same site using a LeiCi DDS-307 (ShengKe, ShangHai, China) to ensure the reliability of the results.

Methods
The technical flowchart of this study shows the main ideas of this research ( Figure 2).

Feature Parameter Extraction
In our study, we used different parameters that can be extracted to measure salinity in the two regions.

Feature Parameter Extraction
In our study, we used different parameters that can be extracted to measure salinity in the two regions.

Modified Soil-Adjusted Vegetation Index (MSAVI)
The modified soil-adjusted vegetation index (MSAVI) considers the bare soil line, which can better eliminate or reduce the influence of soil and vegetation canopy background and thus better indicate the vegetation information in the study area. High soil salinity inhibits crop growth and results in low vegetation cover, which means that MSAVI decreases with increasing soil salinity, and therefore, it can be used as a biophysical parameter to reflect vegetation growth. The MSAVI was used to extract vegetation information, as shown in Formula (1) [22]: where NIR and R represent the reflectivity values of the near-infrared and red bands in the OLI image, respectively.

Albedo
Surface albedo is an indicator of the absorption and reflection ability of the ground to solar radiation and refers to the ratio of surface-reflected radiation to incident radiation. A higher albedo indicates a higher reflectivity and, thus, a smaller absorption capacity of the ground to solar radiation. Surface albedo can reveal information about salinized soil. In this study, we used the inversion formula of albedo established by Liang [19], as shown in Formula (2): where B, R, NIR, SR1, and SR2 correspond to the reflectivity values of blue, red, nearinfrared, short-wave infrared band 1, and short-wave infrared band 2, respectively, in the OLI image.

Normalized Difference Vegetation Index (NDVI)
The normalized difference vegetation index (NDVI) is often used to extract vegetation information to indirectly characterize soil salinity [18]. The NDVI is very sensitive to the performance of green vegetation and can eliminate most of the changes in irradiance related to instrument calibration, solar angle, terrain, cloud shadow, and atmospheric conditions. Higher NDVI values indicate more green vegetation. We used NDVI as an indicator to monitor soil salinization, as shown in Formula (3): where NIR and R represent the reflectivity values of the near-infrared and red bands in the OLI image, respectively.
Salinity Index (SI) The salinity index (SI) is an important factor for describing soil salinization [27]. We used SI to extract salt information, as shown in Formula (4): where B and R represent the reflectivity values of the blue and red bands in the OLI image, respectively.

Data Normalization Processing
The maximum and minimum values of MSAVI, Albedo, SI, and NDVI were calculated using The Environment for Visualizing Images (ENVI) statistical tools. Then, data normalization was performed using the ENVI data, as shown in Formulas (5)-(8): where M is the normalized value of the MSAVI, A is the normalized Albedo, N is the normalized NDVI, and S is the normalized SI.

Feature-Space Principle
The feature space is a spatial system composed of two or more bands or spectral indicators used to measure the spectral characteristics of ground parameters in remote sensing images [28]. For extracting thematic information of salinization and constructing a feature space model, having a reasonable selection of characteristic parameters of the spectral response of salinized soil is vital. We constructed three feature spaces using different combinations of the parameters described above.

Albedo-MSAVI Feature Space
There is a significant relationship between the occurrence of saline soil and the characteristic parameters of the soil spectrum. The Albedo-MSAVI feature space [29] shown in Figure 3a shows that there is a non-linear relationship between albedo and MSAVI. As albedo increased, the degree of soil salinization increased, and MSAVI decreased. Point A in Figure 3a represents heavily salinized land with low vegetation coverage and high surface albedo (low MSAVI and high albedo); point B represents bare land, and point C represents a high vegetation coverage area with relatively low albedo and sufficient soil moisture content (high MSAVI and low albedo). As can be seen, the surface evapotranspiration resistance, soil salt content, and albedo decreased with an increase in surface vegetation coverage.

SI-Surface Albedo (SI-Albedo) Feature Space
As shown in Figure 3b, the SI-albedo feature space demonstrates that SI and surface albedo exhibit significant linear characteristics [19]; where an increase in SI increases surface albedo and the severity of salinization.

SI-NDVI Feature Space
The SI-NDVI feature space [18], shown in Figure 3c, demonstrates a non-linear relationship between SI and NDVI. As salinity content increased, the vegetation coverage declined, and the degree of salinity gradually increased. Point A in Figure 3c represents heavily salinized land with high soil salinity and low vegetation coverage, point B represents bare land, and point C represents a high vegetation coverage area. The SI-NDVI feature space [18], shown in Figure 3c, demonstrates a non-linear re tionship between SI and NDVI. As salinity content increased, the vegetation coverage clined, and the degree of salinity gradually increased. Point A in Figure 3c represe heavily salinized land with high soil salinity and low vegetation coverage, point B rep sents bare land, and point C represents a high vegetation coverage area. To establish a two-dimensional scatter plot, the Albedo-MSAVI feature space inv sion model took normalized MSAVI as the horizontal axis and normalized albedo as vertical axis. Figure 4a and Figure 5a show that there is a non-linear relationship betw surface albedo and MSAVI, similar to that shown in Figure 3a. Taking the MSAVI a albedo coordinates (0, 1) as reference points in Figure 4a, the distance from any poin the scatter diagram to the reference point reflects the degree of salinization. Shorter d tances signified a greater degree of salinization and vice versa. In the two-dimensio feature space, the distance from any point to the (0, 1) point was D, and a remote sens monitoring model was constructed on this basis, as shown in the following equation ( To establish a two-dimensional scatter plot, the Albedo-MSAVI feature space inversion model took normalized MSAVI as the horizontal axis and normalized albedo as the vertical axis. Figures 4a and 5a show that there is a non-linear relationship between surface albedo and MSAVI, similar to that shown in Figure 3a. Taking the MSAVI and albedo coordinates (0, 1) as reference points in Figure 4a, the distance from any point in the scatter diagram to the reference point reflects the degree of salinization. Shorter distances signified a greater degree of salinization and vice versa. In the two-dimensional feature space, the distance from any point to the (0, 1) point was D, and a remote sensing monitoring model was constructed on this basis, as shown in the following Equation (9):

SI-Albedo Feature Space Inversion Model
The SI-Albedo feature space inversion model used the normalized SI as the horizontal axis and normalized albedo as the vertical axis to establish a two-dimensional scatter plot. As shown in Figures 4b and 5b, the scatter diagram shows a band-like distribution like that in Figure 3b. With increasing SI, the degree of salinization increased. The reference point was the point with an SI and surface albedo of 0. The closer any point was to this point in the scatter diagram, the lower the degree of soil salinization and vice versa. The distance from any point in the two-dimensional feature space to the (0, 0) point was D; a remote sensing monitoring model was built on this basis, as shown in Equation (10):

SI-NDVI Feature Space Inversion Model
The SI-NDVI feature space inversion model considers NDVI as the horizontal axis and normalized SI as the vertical axis to establish a two-dimensional scatter plot, as shown in Figures 4c and 5c. Considering the NDVI and SI coordinates (1, 0) as the reference point, greater distances from this point in the scatter diagram indicate greater soil salinization (similar to Figure 3c). In the two-dimensional feature space, the distance from any point to the (1, 0) point was D, and a remote sensing monitoring model was constructed based on this, as shown in Equation (11): greater distances from this point in the scatter diagram indicate greater soil salinization (similar to Figure 3c). In the two-dimensional feature space, the distance from any point to the (1, 0) point was D, and a remote sensing monitoring model was constructed based on this, as shown in equation (11): (a) (b) (c)  greater distances from this point in the scatter diagram indicate greater soil salinization (similar to Figure 3c). In the two-dimensional feature space, the distance from any point to the (1, 0) point was D, and a remote sensing monitoring model was constructed based on this, as shown in equation (11): (a) (b) (c)

Inversion Results and Classifications
In this study, we classified the degree of soil salinity according to the Jenks natural break method. Compared to other classification methods, such as goodness-of-fit, the Jenks natural break classification method has the advantage of having the least variance within classes and the greatest variance between classes and was thus able to minimize the sum of variance within each salinity class [30]. The modeled Albedo-MSAVI, SI-Albedo, and SI-NDVI values were divided into five classes, resulting in five value range intervals. Salinization  Table 1).   Table 2). ArcGIS was used to extract salinity information and generate salinity grade distribution maps for the three monitoring models of Kenli County ( Figure 6) and the Yanqi Basin ( Figure 7). All three models of Kenli County showed that the region has a high level of salinization, among which the Albedo-MSAVI and SI-NDVI inversion distribution maps showed larger areas of moderate and severe salinization, mainly concentrated in the southwest and northeast; mild salinization was less common and mainly distributed in the central region. Overall, the degree of salinization gradually increased from east to west.
In the SI-Albedo inversion distribution map, moderate and severe salinization were mainly distributed in the central region and some coastal areas; mild salinization was mainly distributed on both sides of the Yellow River and near rivers and lakes; and nonsalinized areas were mainly distributed on both sides of the Yellow River, which is consistent with the results of Zhang et al. [31,32]. Moving west from the coast, the degree of salinization increases from mild salinity to moderate salinity and then to severe salinity, and the overall degree of salinization shows a high-low-high trend, meaning, the coastal and inland marginal areas were the most heavily salinized whereas the central areas were mildly salinized.
The distribution of salinization in the Yanqi Basin differs greatly among the three models. In the SI-Albedo inversion distribution map, moderate salinization and severe salinization were mainly distributed in the eastern part of Heshuo County, the Yellow Water Ditch irrigation area, the inlet of Bosten Lake, and the southern shore of Bosten Lake. In contrast, mild salinization was mainly distributed on the banks of the Kaidu River and around the dense drainage system of Bosten Lake. The non-salinization areas were mainly located in the southern part of Yanqi County, the middle and lower reaches of the Nangan canal, and the southeast edge of the basin. From east to west, the salinization degree is distributed along the water system echelon, and the degree of salinization is generally in a high-low-high direction; meaning, the northern and southern marginal areas of the basin were the most salinized, and the central areas were mildly salinized. Lake. In contrast, mild salinization was mainly distributed on the banks of the Kaidu River and around the dense drainage system of Bosten Lake. The non-salinization areas were mainly located in the southern part of Yanqi County, the middle and lower reaches of the Nangan canal, and the southeast edge of the basin. From east to west, the salinization degree is distributed along the water system echelon, and the degree of salinization is generally in a high-low-high direction; meaning, the northern and southern marginal areas of the basin were the most salinized, and the central areas were mildly salinized.
We classified the salinization of the field measurement sample data according to these standards. The accuracy of the Albedo-MSAVI, SI-Albedo, and SI-NDVI feature space monitoring models was determined by comparing the salinity of the measured values with the salinity of the predicted values of each model (Table 3). Table 3. Accuracy verification of the models
We classified the salinization of the field measurement sample data according to these standards. The accuracy of the Albedo-MSAVI, SI-Albedo, and SI-NDVI feature space monitoring models was determined by comparing the salinity of the measured values with the salinity of the predicted values of each model (Table 3). The overall classification accuracy of the Albedo-MSAVI feature space monitoring model of Kenli County was 82.9%, compared to 87.2% for the SI-Albedo feature space monitoring model and 80.9% for the SI-NDVI feature space monitoring model. In the Yanqi Basin, the overall classification accuracy of the Albedo-MSAVI feature space monitoring model was 87.6%, compared to 93.3% for the SI-Albedo feature space monitoring model, and 80.1% for the SI-NDVI feature space monitoring model. Comparing the accuracies of the two study areas, it is clear that the SI-Albedo feature space monitoring model is superior, meaning that SI and albedo are sensitive parameters for inversion of salinization and are applicable to both the coastal and inland arid areas of China.
According to the monitoring model with the highest precision, Si-Albedo, the ratio of salinized soil to total soil and the degree of salinization are more severe in Kenli County than the Yanqi Basin, with approximately 70% of the total area impacted by soil salinization. Among the salinized area, the area of moderate salinization was the largest at 54,234 ha, accounting for 29.5% of the total area of Kenli County, followed by the area of light salinization at 48,371 ha (26.3%), heavy salinization at 34,282 ha (18.7%), and saline soil at 6979 ha (4%). In the Si-Albedo model of the Yanqi Basin, the area of moderate salinization was also the largest at 345,055 ha, accounting for 27.2% of the total area, followed by heavy salinization at 313,877 ha (24.7%), light salinization at 306,924 ha (24.1%), and saline soil at 202,934 ha (15.9%).

Feature Space Model Construction and Analysis of Regional Applicability of Remote Sensing Models
This study shows that using feature parameters with clear biophysical significance to construct a feature space can effectively distinguish between different degrees of salinized land. The remote sensing monitoring model can lessen the problem of insufficient monitoring ability by using soil spectral characteristics. In addition, compared to other methods, the model index of this method is simple, feasible, and easy to extract. The feature spatial model method was accurate for modeling soil salinity in China's typical salinization areas (arid and coastal areas), which is conducive to dynamic monitoring of soil salinization. The SI-Albedo feature space monitoring model was more accurate than the Albedo-MSAVI and SI-NDVI models for both regions, indicating that the SI-Albedo feature space model is suitable for coastal and inland dynamic monitoring of salinization and can quickly and accurately extract salinization information. Based on the feature space approach study, the SI-Albedo inversion model has a higher accuracy compared with the Albedo-MSAVI feature space, which was used in the soil salinization study in arid regions by Feng et al. [30], and the NDVI-SI feature space, which was used in the soil salinization study by Wang et al. [18].
The inversion accuracy of this feature space model may be higher because SI and albedo, two important biophysical parameters, directly reflect the degree of salinization and have higher inversion accuracy than the indirect feature parameters (MSAVI and NDVI) they are combined with in the other models. The proportions of moderate salinity, heavy salinity, and saline soil in the Albedo-MSAVI and SI-NDVI inversion models were greater than the measured data, probably because the remote sensing images were selected in spring, and the vegetation index is less sensitive to the extraction of information from vegetation with weaker growth momentum, thus affecting the SI. In addition, the bands in which albedo and SI are located are more sensitive to soil salinity, which is not susceptible to the interference of vegetation information, reducing precision [35]. Khan et al. [36] found that the red and blue bands of Landsat images are sensitive to the degree of soil salinity. At the same time, it was concluded from experimental spectral analyses that the SI determined using the red and blue bands of remote sensing images can better reflect the degree of soil salinity.

Analysis of Salinization Distribution and Driving Factors and Suggestions for the Governance of Key Areas
The central and southwestern regions of Kenli County in the Yellow River Delta have dense populations, frequent human construction activities, complex land use structures, a large proportion of agricultural land, and a small proportion of green land and forest land; the exploitation of groundwater due to land use is intensive, which increases the salinization level in these areas. The northeastern region contains the mouth of the Yellow River, where the terrain is low, flat, and prone to seawater backflow, and therefore, salinity tends to accumulate here. In addition, groundwater is generally buried at a shallow depth and has a high degree of mineralization. Salt in groundwater accumulates continuously to the surface through capillary ascension, which results in a high degree of salinization [37]. Non-salinization and mild salinization areas are mainly located near rivers, which may be because the fresh water from rivers (especially the Yellow River) influences the soil, crop irrigation, and drainage, making the soil accumulate less salt with a lesser degree of salinization. This is consistent with the findings of Wang et al. [38]. There is an obvious spatial differentiation pattern of soil salinization in the Yanqi Basin; generally, from west to east and from north to south, soil water content gradually decreases, and soil salt content increases. There is a lot of newly reclaimed arable land in the upper reaches of the Kaidu River, the Yellow Water Ditch irrigation area, and the south shore of Bosten Lake, which poses the problem of leaking canal systems. This leaking increases irrigation and drainage on farmland and raises the groundwater level, which drives evaporation. The salt content of the parent material of the Piedmont proluvial-alluvial fan increases the salinization of the area.
Due to their lower topographic parts and generally less than 3 m burial depth of groundwater [39], the central alluvial and lakeshore plains are prone to backwater formation, dry climate, strong evaporation, and saline groundwater accumulation to the surface through capillary upwelling, which makes salinization more severe. The non-salinized and lightly salinized areas were mainly located near rivers, which may be because these areas are well equipped with irrigation and drainage facilities, and the land has been improved for many years, making the soil less saline. This is consistent with the results of Wang et al. [40] and Jiang et al. [41] on farmland soil salt distribution patterns and remote sensing inversion in the Yanqi Basin.
Next, we will propose recommendations for the prevention and control of salinity for protecting ecosystems and maintaining sustainable agricultural development. For Kenli County, areas with slightly salinized soil with salt content less than 3‰ and good forest site conditions, we recommend implementing simple plowing and leveling and increasing the application of organic fertilizer. For low-lying, poorly drained, and large square connected moderate saline-alkali land, we recommend changing the deep ditch and strip field to a strip field, upper forest, and lower fishing mode. For heavily saline or low-lying land with low terrain and a high water table, we recommend constructing raised fields according to the requirements of planning and design. Where the water table is not too high, the terrain is flat, the soil salt content is high, and the area is large, we recommend implementing open ditch plus concealed pipe alkali drainage projects. For the greening of important nodes such as roads and intersections in the suburbs of cities, the blind ditch drainage model should be used to improve the soil. The southwestern region, the central region, and other areas with lighter soil salinization should pay attention to field preparation, deep pine, and leveling of the soil over time. The eastern region with more serious soil salinization, low-lying terrain, and poor conditions for agricultural development should take advantage of the proximity to the sea to develop farming. To build a complete irrigation and drainage system, the middle region should level the land, set up salt drainage ditches to lower the groundwater level, and improve the saline-alkali soil using techniques such as tree pit cushion.
For the Yanqi Basin, in the areas with lightly saline soils in the Kaidu River alluvial plain, agricultural and biological measures should be used to improve saline soils, regulate soil water and salt dynamics, apply more organic fertilizers, loosen the soil tillage layer, ensure the growth and development of productive crops, enhance the effect of salt drenching, reduce evaporation, and inhibit salt reversion. Planting green manure is a key method for improving saline-alkali soils. It can improve the physical and chemical properties of soil, consolidate, enhance the desalination effect, and fertilize the soil. The treatment and improvement of moderate saline-alkali land should use chemical amendments to improve soil, and chemical amendments such as gypsum and zeolite which can effectively reduce the salinity of soil saline-alkali land [42]. In terms of the development and utilization of saline-alkali land and halophytes, the Xinjiang Institute of Ecology and Geography, Chinese Academy of Sciences, sampled many research results and determined that saline-alkali land and halophytes are valuable natural resources. The key to their use is rational development and effective utilization. Physical improvement measures such as land leveling and soil moving improvement can be used to improve severely saline-alkali land, especially land that is difficult to control, such as that at the estuary of Kaidu River and the salt collecting area of Yellow Water Ditch. In addition, using irrigation water to wash the salt to the subsoil layer and draining dissolved salts with a drainage ditch can improve saline-alkali soil. While channel infiltration should be implemented and irrigation management strengthened to wash nutrients from the soil into the groundwater, the drainage system should be built to control the groundwater level and ensure the stable desalination of the soil [43].
Based on the above study and in view of the problems of different types of soil salinization and industrial ecologicalization, four aspects for improving salinized soil should be considered and adopted according to local conditions: water conservancy improvement measures, agricultural improvement measures, biological improvement measures, and chemical improvement measures. In addition, saline land management should be elevated to the level of priority as national food security. The government can set up special public welfare funds to form a complete system for saline land restoration.

Conclusions
In this study, we combined Landsat 8 OLI remote sensing and field sampling data to address the persistent soil salinization problem in coastal (humid zone) and arid areas (dry zone) of China and constructed Albedo-MSAVI, SI-Albedo, and SI-NDVI feature space inversion models. Our study was the first comparative analysis of the generalizability of feature space inversion models of typical salinization regions of China. The inversion accuracy of the SI-Albedo monitoring model was higher than that of the Albedo-MSAVI and SI-NDVI models in both regions, showing applicability for dynamic monitoring of salinity in coastal and arid regions. We used SI-Albedo to generate a 30-m resolution map of soil salinization in the Yellow River Delta in 2019 and Yanqi Basin in 2020 and analyzed the patterns and drivers of salt distribution. Based on these results and experience, we proposed corresponding sustainable development countermeasures for different types and severities of soil salinization.
However, this technology must be improved. The relationships between remote sensing indices and soil salinity within fields are highly complex; they require further investigation with additional samples and various soil salinity classifications should be used. In this study, the number of study areas was low. Therefore, areas with diverse geographic locations within the same climate zone should be studied to verify the model's general admissibility. Additionally, in this study, sampling was conducted in spring; therefore, seasonal differences in the salinization phenomena could not be considered.
Thus, the salinity grade map obtained by inversion only represented the current salinity distribution in the region in spring.
In our study, only two climate zones (arid and humid) were investigated. In the follow-up work, we will evaluate the accuracy of salinity inversion for different remote sensing models in other climatic zones transition zone and Qinghai-Tibet Plateau area) in China to more comprehensively assess the applicability of the SI-Albedo remote sensing model across the country.
Author Contributions: Y.F. and H.W. were responsible for the research design, analysis, manuscript design and review; J.L. drafted the manuscript and was responsible for data preparation; experiments and analyses were conducted by J.W.; L.Z. was responsible for funding acquisition and resources; Q.G. was responsible for the research design field survey; manuscript review was done by T.D.; and data curation was completed by Z.Z. and Q.Y. All authors contributed to the editing and reviewing of the manuscript. All authors have read and agreed to the published version of the manuscript.