1. Introduction
Soil salinization represents a widespread global phenomenon associated with soil degradation, posing significant challenges to both agricultural development and ecological stability. China, hosting the world’s third-largest area of saline soils spanning 3.69 × 10
7 hm
2, equivalent to approximately 4.88% of its available land area, is particularly impacted [
1]. Among the salt-affected regions in China, the Songnen Plain stands out with a saline–alkali area of 3.73 × 10
6 hm
2, primarily consisting of sodium bicarbonate and sodium carbonate as the main salt minerals. Due to the relative high clay content and adsorbable cations, the salt-affected soils in the Songnen Plain suffer from a serious desiccation cracking process during water evaporation, suggesting that exploring the relationship between crack extents and salinity levels can therefore help to expeditiously and accurately obtain soil salinity data, assess the degree of soil salinization, and delineate the extent of salinized soil, which is also pivotal for guiding soil amelioration efforts, optimizing soil land utilization, and safeguarding the local food security.
Field sampling followed by traditional laboratory analysis is considered the most direct method for detecting soil salinity [
2,
3]. However, these conventional methods are not efficient in providing real-time soil salinity parameters due to their time-consuming testing processes, long measurement periods, and high labor costs. To overcome these limitations, the electromagnetic induction method is commonly employed to determine salt content without physical contact. This method involves measuring the correlation between the primary magnetic field and the induced secondary magnetic field to quickly assess the apparent soil conductivity [
4,
5]. While electromagnetic induction sensors like the EM-38 are widely used for monitoring soil salinity, soil texture, and other soil properties in the field [
6,
7], their accuracy is affected by the physical and chemical properties of different soil types, which are also sensitive to environmental influences [
8,
9]. In contrast, hyperspectral remote-sensing technology offers a powerful and effective approach for monitoring soil salinity. This technique relies on the diagnostic spectroscopy of different salt minerals and utilizes the wide band range and high spectral resolution of hyperspectral sensors [
10,
11]. Numerous studies have emphasized the significance of hyperspectral remote sensing in extracting and predicting soil salinity levels. These studies have provided valuable insights into the spectral characteristics and sensitive bands associated with different levels of salinity [
12,
13,
14,
15]. However, it should be noted that previous ground-based spectral measurements have mainly focused on soil samples obtained through sieving (particles with sizes of 2 mm or less), and spectral data derived from large-scale remote-sensing images often encounter mixed pixel issues, limiting their ability to accurately represent the actual surface conditions of salt-affected soils.
Under natural conditions, it is a common occurrence for clayey soil to undergo shrinkage and develop desiccation cracks during the process of water evaporation. This phenomenon is influenced by several soil properties, including soil salinity, clay content, and mineralogy. Specifically, soils with a higher clay content are more susceptible to severe shrinkage deformation, resulting in more prominent desiccation cracks [
16,
17,
18]. The presence of different clay minerals can also impact swell–shrinkage cracking, as indicated by their respective plasticity indexes. Soils with lower smectite content and higher amounts of kaolinite, illite, and quartz tend to exhibit smaller shrinkage strains [
19,
20,
21]. Regarding salt-affected soils, soil salinity also has a significant influence on the cracking process of cohesive salt-affected soils. Extensive research has demonstrated that changes in soil physicochemical parameters, specifically the type of salt minerals and the exchangeable sodium percentage, play a crucial role in the cracking process. In a study conducted by Ren et al. [
22], the authors utilized a large number of soil samples collected under natural conditions to establish a relationship between the texture features of desiccation cracks and the salinity of soda saline–alkali soils. They employed a four-level orthogonal wavelet decomposition and found that high-frequency texture features, such as energy and the L1 norm, exhibited an exponential correlation with EC
1:5 values of salt-affected soils. Another study by Zhang et al. [
23] focused on predicting the main salt parameters of salt-affected soils in the Songnen Plain of China. They compared three different methods and discovered that total salinity promotes the process of desiccation cracking in soda saline–alkali soils. They also demonstrated that the artificial neural network (ANN) method can accurately predict Na
+ and EC
1:5 levels. Fujimaki et al. [
24] investigated the relationship between the crack characteristics and physicochemical parameters of saline soils. They evaluated EC
1:2 values in agricultural fields with desiccation cracks in the Nile Delta and found that cracks facilitate the rapid movement of soil salts through water, thereby influencing the distribution of soil salts. Additionally, Zhang et al. [
25] conducted a study on pure clay of bentonite and observed that the salt concentration significantly induced and altered the morphology and patterns of desiccation cracks. As the concentration of the sodium chloride (NaCl) solution increased, larger crack lengths and lower crack densities were obtained. In order to gain a comprehensive understanding of the mechanism underlying desiccation cracking in saline–alkali soil and provide valuable insights for engineering practices, Xing et al. [
26] conducted a study wherein they treated soil samples with saturated solutions of NaCl at concentrations of 5 g/L, 50 g/L, and 100 g/L. The results of their experiments demonstrated a gradual decrease in both soil-crack area and crack length with increasing Na
+ concentration. However, an opposite trend was observed for crack number, which increased with increasing Na
+ concentration. Wang et al. [
27] suggested different results in their research, wherein NaCl treatments were found to increase both the total area and average width of desiccation cracks in shrink–swell light clays. This discrepancy highlights the complexity of the desiccation cracking process and emphasizes the need for further investigation. Although previous studies have examined the desiccation cracking process under varying salinity conditions, the specific mechanism through which salinity influences the soil-cracking process remains unresolved. This uncertainty may be attributed to the complex interactions between different cations and soil particles, which are influenced by factors such as ion concentration, valence state, topography, climate, and soil type.
Although numerous studies have investigated the effects of soil salinity and salt mineral composition on the cracking process of different salt-affected soils, most of these studies have focused on determining the control of the soil-cracking process and its impact on the extent of soil surface cracking through experimentation with a limited number of sample variations under controlled desiccation cracking conditions. However, there is a lack of research that aims to predict soil salinity in real-world scenarios based on crack characteristics. Furthermore, the hyperspectral characteristics of salt-affected soil, which provide important information about the type and content of soil salt minerals, have not been extensively studied in relation to cracked saline–alkali soil. Additionally, the complexity of cracks can further influence the spectral response of the soil surface, while the crack area can also affect the volume scattering process of electromagnetic radiation. Therefore, integrating crack characteristics and the spectral response through ground-based remote-sensing methods can greatly improve the accuracy of soil salinity content prediction, thereby facilitating soil salinity monitoring. This study attempts to perform desiccation cracking experiments under controlled indoor conditions using soil samples with different salinity levels. Three types of crack parameters, including crack length, crack area, and texture feature, were quantitatively extracted from the crack patterns of soil samples under completely dry conditions. Meanwhile, spectral data of cracked soil samples were also measured in this study. After that, structural equation modeling (SEM) was explored to realize single-parameter and multi-parameter prediction methods for soil salt parameters including EC1:5, pH, Na+, HCO3−, CO32−, and the total salinity by combing the crack characteristics and spectral responses of cracked saline–alkali soils. This not only enables the effective prediction of salt information in cracked soil from a ground-based remote-sensing perspective but also enhances the experimental process and scheme design for the hyperspectral measurement of saline–alkali soil surfaces. As a result, this study provides valuable data support for local soil remediation efforts and strengthens the remediation of soda saline–alkali soil in the Songnen Plain, China. Furthermore, this approach propels the field of hyperspectral remote sensing for soil salinization monitoring and extends the application of soil hyperspectral remote-sensing technology. Additionally, this research also attempts to improve the experimental process and scheme design for hyperspectral measurement of saline alkali soil surfaces, providing crucial data for local soil remediation efforts and advancing the application of hyperspectral remote sensing for soil salinization monitoring.
4. Discussion
Soil cracking is a complex mechanical process that is influenced by various factors, including soil texture, soil chemical properties, and environmental conditions. In this study, the environmental factors, such as temperature and relative humidity during the drying process, were maintained at a constant level, suggesting that the impact of experimental conditions and sample preparation can be disregarded. Moreover, the distribution of soil samples in this study was relatively concentrated, indicating that the samples had similar clay content and mineral types. From the findings presented in
Table 4, it can be observed that soil texture had a minimal effect on the desiccation cracking process. The correlation coefficients between clay content and crack characteristics ranged from only 0.13 to 0.26, suggesting a weak relationship. Similar results were reported by Zhang et al. [
51] in their study on soda saline–alkali soil in the Songnen Plain. They found that the clay content in the soil ranged from 20% to 29%, which aligns closely with the measurements obtained in this study. They also discovered that the activity index of salt-affected soil in the Songnen Plain ranged from 0.33 to 0.48, indicating that clay particles have little influence on the expansion and contraction of the salt-affected soil in this region. This can be attributed to the predominant presence of illite and montmorillonite mixed layers, with a mixed-layer ratio greater than 50%, along with smaller amounts of illite and kaolinite. Notably, the montmorillonite content, which exhibits the highest activity, is evidently low.
From
Table 4, it can be seen that soil salinity plays a very important role in the process of desiccation cracking, which may be attributed to several factors. In soda saline–alkali soils, the presence of exchangeable cations in the salt solution is a crucial factor that determines the physical properties of the soil during the process of dehydration and cracking. This, in turn, influences the cracking of the soil surface due to exposure to salt. From
Table 1, it is evident that Na
+ is the predominant cation in the Songnen Plain of China, with a significantly higher content than other types of soil cations. The surface of cohesive saline soil samples generally carries negative charges, which allows for the exchange with exchangeable cations [
52]. The alteration in charge distribution further affects the adsorption mechanism of water molecules and consequently leads to changes in their adsorption behavior [
53]. Specifically, the arrangement and behavior of water molecules are influenced by cations, particularly Na
+, resulting in the formation of a specific bound water film on the soil surface [
54]. This kind of water film acts as a dispersant, increasing the distance between soil particles and reducing their contact and bonding. As a result, the cohesion and tensile strength of the soil are reduced. With the further thickening of the bound water film, soil particles become more susceptible to sliding or separating under stress, leading to a decrease in the overall structural stability and strength of the soil [
55,
56,
57]. Additionally, the lubricating effect of the bound water film between soil particles decreases the internal friction angle and shear strength of the soil, promoting the desiccation cracking process of salt-affected soils [
58]. Moreover, the content of monovalent Na
+ in the soda saline–alkali soil solution in the Songnen Plain is significantly higher than that of the divalent Ca
2+ and Mg
2+, which can negatively impact the stability of soil aggregates. This weakens the bonding force between soil particles and makes them more easily separated by water, thereby promoting the process of soil surface cracking [
59]. In addition, salinity reduces the water absorption capacity of the soil. This may be because clay minerals swell when the concentration of sodium ions in the soil is too high, leading to a reduction in soil pore space. This reduces the rate of water infiltration and permeability and disperses soil particles, further decreasing pores and impeding the flow of water and air. Consequently, the soil becomes more prone to drying and shrinkage, thereby facilitating the process of desiccation cracking [
60]. Notably, studies have suggested that the diffuse double layer (DDL) also plays a critical role in the process of desiccation cracking of saline–alkali soils during water evaporation. Specifically, as water evaporates, the DDL thins, reducing the distance between soil particles and facilitating the propagation of desiccation cracks on the soil surface [
61]. Therefore, a higher salt content leads to stronger cohesion between soil particles, resulting in increased soil volume shrinkage and more complex soil-cracking phenomena.
It has been shown that the reflectance of salinized soil surface decreases with increasing soil salinity [
62,
63]. It is also affected by factors such as soil moisture content, surface roughness and organic carbon [
64]. Soil cracking increases the roughness condition of the surface, leading to an increase in its scattering at the surface. At the same time, electromagnetic radiation will also cause multiple body scattering in the cracked region, leading to a significant attenuation of electromagnetic radiation energy, suggesting that the presence of cracks significantly reduces the reflectivity of the salinized soil surface. In addition, for samples with different salinity levels, the area of the cracked region increases due to the increasing salinity. The attenuation of electromagnetic radiation is more intense compared to the low-salinity region, which further reduces the spectral reflectance of the samples and also increases the reflectance difference between samples. In the present study, a desiccation cracking experiment was conducted under controlled laboratory conditions. The measurement environment had the same effects on the different samples during the whole drying process, and the moisture content of the samples was so low that its effect on the spectral reflectance can thus be ignored. It was shown that the main factor affecting the spectral reflectance was the salt content. From the correlation curves of spectral reflectance and physicochemical parameters of soil samples, it can be found that the 1990 nm band can be selected as the characteristic band, which corresponds to the similar findings reported by Ren et.al. [
22].
Figure 5 indicates that when the salinity increased, the reflectance showed a decreasing trend, which was especially obvious at wavelengths greater than 1400 nm; this finding is consistent with the findings of Wang et al. [
65] and Howari et al. [
66]. It is because NaHCO
3 and Na
2CO
3 are the main components of salt minerals in saline soils. In addition, according to Wang et al. [
65], Na
2CO
3 saline soil has the lowest albedo shape compared to Na
2SO
4 saline soil, NaCl saline soil and non-saline soil. Their results also showed that at wavelengths above 1400 nm, reflectance decreased significantly with increasing Na
2CO
3 content. Howari et al. also [
66] showed that the sharp decrease in the reflectance of NaHCO
3 and Na
2CO
3 with increasing wavelength due to unique resonances induced by stretching and bending can be considered diagnostic. In addition, according to Tedeschi and Dell’Aquila [
67] and Huang et al. [
68], the aggregation capacity of the soil decreases with soil salinity, probably due to the salt-induced anti-flocculation of the mucilaginous colloids, leading to the inverse ratio of salinity to spectral reflectance.
Figure 5 also shows that the soil reflectance tended to increase at 350–1400 nm, while the increase slowed down at 1400–2200 nm and a decreasing trend was observed at 2200–2500 nm. The reason for this is that the reflectance of NaHCO
3 in soil shows a significant decreasing trend at wavelengths greater than 1400 nm, and the trend becomes more obvious with the increase in soil salinity.
The use of SEMs has gained significant popularity and has been widely applied in various fields, including the social sciences, economics, and psychology. In recent years, there has been an increasing interest in using SEMs in the field of soil science [
69,
70]. In this study, all the regression models developed in
Table 5 and
Table 6 demonstrated high accuracy in modeling the total soil salinity, EC
1:5, and Na
+ levels. This can be attributed to the fact that Na
+, being the predominant cation in saline–alkali soils of the Songnen Plain, is closely associated with soil salinity. The EC
1:5 of the soil is mainly influenced by the presence of soluble ions, making it a commonly used indicator for assessing soil salinity in salt-affected soils. These findings highlight the strong predictive capabilities of SEMs for salinity, EC
1:5, and Na
+. With regards to the prediction of other ions, the overall accuracy of the SEMs displayed a trend, with CO
32− > HCO
3− > pH. This can be attributed to the hydrolysis of CO
32− and HCO
3− ions in solution, which is often incomplete and reversible. As a result, the content of CO
32− and HCO
3− ions remains highly unstable. Additionally, the pH level is directly determined by the OH
− content generated through the hydrolysis of CO
32− and HCO
3− ions.
From
Table 1 and
Table 7, it can be seen that the soil samples are very representative, and the model predictions are very close to the actual measurements. In addition,
Table 7 also indicates that the SEM is able to combine different crack parameters and spectral responses to quickly and effectively achieve the synchronous prediction of salt parameters, including Na
+, EC
1:5, and salinity, with desirable prediction accuracy. However, several complex statistical assumptions still need to be clarified and satisfied when an SEM is employed, and multicollinearity and data bias may also affect the precision and explanatory nature of SEMs. In addition, the complexity of SEMs may increase the difficulty and resource requirements for model interpretation, validation, and computation. Therefore, in future studies, both the types of crack parameters and the number of diagnostic spectral bands can be further considered in an attempt to improve the application potential and accuracy of the method in online field measurements of soil salt parameters in the Songnen Plain. Meanwhile, multi-source hyperspectral remote-sensing data from UAV measurements and satellite images can also be added in future studies to realize the high-precision, large-scale and fast-rate inversion of regional soil salinity, thus promoting the application development of soil salinity monitoring using hyperspectral remote sensing.