Quantitative Analysis of Spectral Response to Soda Saline-AlkaliSoil after Cracking Process: A Laboratory Procedure to Improve Soil Property Estimation

: Cracking on the surface of soda saline-alkali soil is very common. In most previous studies, spectral prediction models of soil salinity were less accurate since spectral measurements were usually performed on 2 mm soil samples which cannot represent true soil surface condition very well. The objective of our research is to provide a procedure to improve soil property estimation of soda saline-alkali soil based on spectral measurement considering the texture feature of the soil surface with cracks. To achieve this objective, a cracking test was performed with 57 soil samples from Songnen Plain of China, the contrast (CON) texture feature of crack images of soil samples was then extracted from grey level co-occurrence matrix (GLCM). The original reﬂectance was then measured and the mixed reﬂectance considering the CON texture feature was also calculated from both the block soil samples (soil blocks separated by crack regions) and the comparison soil samples (soil powders with 2 mm particle size). The results of analysis between spectra and the main soil properties indicate that surface cracks can reduce the overall reﬂectivity of the soda saline-alkali soil and thus increasing the spectral di ﬀ erence among the block soil samples with di ﬀ erent salinity levels. The results also show that both univariate and multivariate linear regression models considering the CON texture feature can greatly improve the prediction accuracy of main soil properties of soda saline-alkali soils, such as Na + , EC and salinity, which also can reduce the intensity of ﬁeld spectral measurements under natural condition.


Introduction
Soil salinization is one of the most common worldwide land degradation processes, which seriously affects major soil degradation phenomenon such as soil dispersion, soil erosion, and engineering problems, and thus reduces agricultural production, and impacts environmental health [1][2][3][4]. According to Food and Agriculture Organization of the United Nations, the global surface area affected by soil salinization was more than 831 million ha by 2000, covering more than 6% of the Earth's surface [5]. Therefore, identification of the soil salinity and the salt-affected areas is essential for sustainable agricultural management, especially in arid and semi-arid areas [6,7].
Due to its high temporal sampling frequency and rapid nondestructive measurement, reflectance spectroscopy has been widely used to estimate soil salinity [8][9][10][11][12][13]. Hyperspectral data has become a promising source of data in many studies since it contains a large number of narrow contiguous spectral bands with high spectral resolution and thus enables the quantitative measurement of saline soils with different salinity levels [14][15][16][17]. During the last several decades, many studies have been focused on the effect of salinity on the spectral characteristics of saline soils. Taylor et al. [18] found absorption valleys at 980 nm, 1170 nm, 1450 nm, and 1900 nm, and strong reflection peaks at 800 nm. Farifteh et al. [19] carried out a laboratory experiment involving soils with six salt minerals and three soil textures in order to study the quantitative relationship between spectral response and salt concentration of saline soil samples, their results showed that the observed absorption features were broadened at wavelength above 1300 nm and the overall reflectance changed proportionally as the soil salinity increased. Wang et al. [20] carried out a spectra measurement in of saline soil, the laboratory and found significant absorption features at 505 nm, 920 nm, 1415 nm, 1915 nm, and 2205 nm. Moreira et al. [21] evaluated the variations in spectral reflectance of soil samples treated with different salt crystals using both principal component analysis and the continuum-removed method, after that they found the obvious absorption depths of major features at 1450 nm, 1950 nm, 1750 nm, and 2200 nm, their results also indicated the best negative correlations between reflectance and EC in the 1500-2400 nm range for CaCl 2 and MgCl 2 .
In many studies, spectral indices were commonly used to process the hyperspectral data of saline soil after spectral measurements. Zovko et al. [22] selected 246 soil samples for visible and near-infrared spectra measurements and defined a new geostatistical spectral index as covariate in ordinary cokriging for predicting electrical conductivity, spatial variation of soil salinity was then mapped in the Neretva River Valley, Croatia to propose an effective approach based on VNIR spectroscopy and geostatistics for mapping soil salinity. In order to predict the salt abundance in slightly saline soils in Burg Al-Arab, Egypt, Masoud [23] carried out spectrometry measurements for 21 soil samples with different salinity levels to evaluate the simple wetness index (WI) and sophisticated spectral mixture analysis techniques including mixture tuned matched filtering (MTMF) and linear spectral unmixing (LSU).The prediction results of salinity levels showed significant high accuracy with R 2 of 0.88, 0.84 for MTMF, LSU and WI, respectively. After comparing three spectral indices including difference index (DI), normalized difference index (NDI) and ratio index (RI). Nawar et al. [24] developed regression models with the EC in the El-Tina Plain, Egypt, their results showed RI can provide a reliable estimations of the soil EC with produced R 2 , RMSE and RPD values of 0.8, 6.2, and 2.06.
Due to its superiority in dealing with high dimensional multi-collinearity, partial least square regression (PLSR) has been used as a very common multivariate statistical method in predicting the soil salinity. Weng et al. [25] measured the reflectance spectra (350-2500 nm) from EO-1 Hyperion and tested PLSR as prediction model to evaluate the potential of predicting salt content in the Yellow River Delta of China, their results indicated that spectral bands at 1487-1527 nm, 1971-1991 nm, 2032-2092 nm, and 2163-2355 nm possessed large absolute values of regression coefficients. After collecting in situ field spectra of saline soils in the Yellow River Delta of China, Fan et al. [26] established a PLSR model between soil salinity and field spectra extracted from Advanced Land Imager (ALI), the effectiveness of PLSR model was then demonstrated in retrieving soil salinity from new-generation sensors, their results indicated that the prediction accuracy was quite acceptable with R 2 , RPD of 0.749, 3.584, respectively. In addition, many other multivariate statistical techniques such as artificial neural networks (ANN), principal component analysis (PCA), and multiple adaptive regression spline (MARS) were used for quantitative determination of soil salinity and illustrating the relationships between spectral characteristics and soil properties [27][28][29]. Although the spectral characteristics of soil samples can reflect the difference of soil salt content to a certain degree, only soil samples after grinding over 2 mm sieve were used in most previous spectral measurements, which usually cannot represent the true surface state of saline soils. For many saline soils with complex surface conditions or Remote Sens. 2019, 11, 1406 3 of 18 high roughness, spectral measurement of 2mm soil samples often cannot give a good estimation of salt content.
For saline soil samples, especially for those soda saline-alkali soils with relatively high clay content, it is very common to shrink and crack for the soil surface during water evaporation. Volume scattering was generated in the crack regions and result in a significant mixed pixel spectral effect. However, studies on the correlation between salt content and the overall spectral reflectance of saline soils with desiccation cracks are still limited and qualitative [30,31]. Texture features always contain important information about the structural arrangement of the soil surface with desiccation cracks, which can also accurately reflect the hue information and the overall arrangement of crack patterns [32][33][34][35]. However, spectral measurements considering the texture features of soil surface cracks were rarely focused on. To address this gap, the quantitative relationship between the reflectance spectra and soil properties of soda saline-alkali soils were tried to be analyzed based on consideration of the grey level co-occurrence matrix (GLCM) texture feature. Optimal predictive models for four main soil properties of interest including Na + , EC, pH and total salinity were also developed based on the original reflectance and the mixed reflectance considering the contrast texture feature (CON) for the block soil samples (soil blocks separated by crack regions generated during the desiccation cracking process) and the comparison soil samples (soil powders with 2 mm particle size) respectively. This can be used to improve the prediction accuracy of soil salinity with less work intensity of soil spectral measurement in field conditions.

Study Area and Sampling
Songnen Plain contains one of the most three largest amounts of heavy soda salt-affected soils in the world. The main salt components of saline soils are NaHCO 3 and Na 2 CO 3 , which make the soils very sodic [36]. This study was conducted from 2013 to 2014 with all soil samples obtained in 12 towns of Da'an City, which is located on the southwestern Songnen Plain, Jilin Province, China with the geographic coordinates lie between 123 • 47'35"E~123 • 59'55"E and 45 • 29'30"N~45 • 36'26"N. The annual mean precipitation of this area is approximately 400~500 mm with 80% occurs between June and September, while the average annual evaporation is about 1400~1900 mm, reaching 3.5 to 4.75 times of precipitation. Moreover, the soil clay content in the study area is quite high, which makes the soil surface shrink and crack very commonly during water evaporation.
After consideration of heterogeneity, 57 composite soil samples of different salinity levels were collected at a depth of 20 cm in July 2013, April and October 2014. Figure 1a,b shows the study area location and the spatial coordinates of all sampling locations recorded using a hand-held GPS. The soil surface from Figure 1c-e indicates that it is very common to crack on the surface of soda saline-alkali soils in Songnen Plain. All soil samples were air-dried, ground, passed through a 2 mm sieve and then divided into three subsamples. The first subsamples were prepared for the physical and chemical properties. The second subsamples were prepared as the block soil samples the third subsamples were prepared as comparison soil samples for spectral collection.

Soil Property Measurement
The main ions, such as K + , Na + , Mg 2+ , Ca 2+ , Cl − , CO 3 2− and HCO 3 − were measured in a laboratory. Soil extracts with water/soil mass ratio of 1:5 were used to measure the contents of Na + and K + using a flame photometer. The contents of Ca 2+ and Mg 2+ were determined by the complexometric ethylenediaminetetraacetic acid (EDTA) titration method. The contents of CO 3 2− and HCO 3 − were obtained based on the double indicator neutralization method. Moreover, EC and pH were also measured by the potentiostatic method and conductometry, using soil suspensions with a soil/water mass ratio of 1:5 [37]. In addition, a Malvern MS-200 laser particle size analyzer was used in this study to determine the contents of clay, silt, and sand of the soil samples [38].
location and the spatial coordinates of all sampling locations recorded using a hand-held GPS. The soil surface from Figure1c-e indicates that it is very common to crack on the surface of soda saline-alkali soils in Songnen Plain. All soil samples were air-dried, ground, passed through a 2 mm sieve and then divided into three subsamples. The first subsamples were prepared for the physical and chemical properties. The second subsamples were prepared as the block soil samples the third subsamples were prepared as comparison soil samples for spectral collection.

Desiccation Cracking Test
In order to remove the effects of environmental factors, the desiccation cracking process of soda saline-alkali soils was carried out under controlled experimental conditions. Saturated slurry samples were firstly prepared with a mass water content of 80% using the subsamples treated previously. The obtained slurry samples were then well hand-mixed and poured into wooden sample box of 50 cm × 50 cm × 3 cm size with the top flattened by a spatula. Afterwards, all soil samples were subjected in the laboratory with temperature of 25 • C, humidity of 35% and pressure of 101 kPa, and the soil samples were then air-dried until masses no longer decreased, which can be considered as completely dried samples. The gravimetric water contents of the soil samples were then measured. In order to obtain the crack images of soil samples, a digital camera was installed on a fixed experimental platform. The lens of the digital camera was set 1 m above the floor. A 50 cm × 50 cm square area was determined with the projection of the camera lens on the floor as the center, which was used to ensure that the geometric distortion effects were the same for all the samples during the photographing process. Each cracked soil specimen was then placed in the determined square area to be photographed ( Figure 2

Soil Property Measurement
The main ions, such as K + , Na + , Mg 2+ , Ca 2+ , Cl − , CO3 2− and HCO3 − were measured in a laboratory. Soil extracts with water/soil mass ratio of 1:5 were used to measure the contents of Na + and K + using a flame photometer. The contents of Ca 2+ and Mg 2+ were determined by the complexometric ethylenediaminetetraacetic acid (EDTA) titration method. The contents of CO3 2-and HCO3 -were obtained based on the double indicator neutralization method. Moreover, EC and pH were also measured by the potentiostatic method and conductometry, using soil suspensions with a soil/water mass ratio of 1:5 [37]. In addition, a Malvern MS-200 laser particle size analyzer was used in this study to determine the contents of clay, silt, and sand of the soil samples [38].

Desiccation Cracking Test
In order to remove the effects of environmental factors, the desiccation cracking process of soda saline-alkali soils was carried out under controlled experimental conditions. Saturated slurry samples were firstly prepared with a mass water content of 80% using the subsamples treated previously. The obtained slurry samples were then well hand-mixed and poured into wooden sample box of 50 cm × 50 cm × 3 cm size with the top flattened by a spatula. Afterwards, all soil samples were subjected in the laboratory with temperature of 25 °C, humidity of 35% and pressure of 101 kPa, and the soil samples were then air-dried until masses no longer decreased, which can be considered as completely dried samples. The gravimetric water contents of the soil samples were then measured. In order to obtain the crack images of soil samples, a digital camera was installed on a fixed experimental platform. The lens of the digital camera was set 1 m above the floor. A 50 cm × 50 cm square area was determined with the projection of the camera lens on the floor as the center, which was used to ensure that the geometric distortion effects were the same for all the samples during the photographing process. Each cracked soil specimen was then placed in the determined square area to be photographed ( Figure 2).

Texture Feature Extraction
The CON is an important GLCM texture feature，which can well express the clarity and reflect the intensity of the texture of adjacent pixels of crack patterns, thus the CON is considered as an

Texture Feature Extraction
The CON is an important GLCM texture feature, which can well express the clarity and reflect the intensity of the texture of adjacent pixels of crack patterns, thus the CON is considered as an important measure of local variations present in crack patterns in this study. The CON can be calculated by the following equation: where p(i,j) represents the second-order combined conditional probability density of gray level i and j, corresponding to the element of GLCM at (i,j) position. For an Lr × Lc grayscale image with gray level of Ng and sampling interval of d pixels, p(i,j)at directions of 0 • , 45 • , 90 • , and 135 • can be computed by the following equations separately [39][40][41][42]: In this study, the crack images were transformed into binary images (Ng = 2) considering the effects of the salt crystals separated out on the surface of soil samples. In detail, to determine the CON of dried soil samples from crack patterns, digital image processing technique was carried out in this study ( Figure 3). Firstly, the crack pattern of the cracked soil samples was geometry corrected using polynomial methods and cropped ( Figure 3a).After that, the corrected color image was transformed into a gray image ( Figure 3b).Thirdly, the grayscale image was binarized using a threshold ( Figure 3c) and then inversed for further processing (Figure 3d), which indicated that the crack pattern was segmented into white crack areas and black block areas. Finally, the open operation was then performed with a given threshold to reduce the small spot noise in every binary image ( Figure 3e). After the above processes, a distance of 40 pixels was selected as sampling interval in order to fully preserve the GLCM texture information and reduce the computational complexity, the mean CON texture features of all cracked soil samples were then computed based on the CON computed from GLCMs in 0 • , 45 • , 90 • , and 135 • directions [43].

Spectral Features Acquisition
In order to obtain the spectral characteristics of both the block soil samples and the comparison soil samples, an ASD FieldSpec3 hand-held spectrometer produced by Analytical Spectral Devices Inc. was used in this experiment, which was widely used in the spectral measurements of vegetation and soil. This spectrometer covers wavelengths between 350-2500 nm with two separate detectors namely VNIR (visible and near-infrared, 350-1000 nm) and SWIR (short-wave infrared, 1000-2500 nm), it can savethe spectral data and display the spectral curve of the objects. The sampling intervals are 1.4 nm with a spectral resolution of 3 nm and 2 nm for VNIR, and with a spectral resolution of 10 nm SWIR, respectively. All spectral reflectance was measured in clear and cloudless weather conditions.

Spectral Features Acquisition
In order to obtain the spectral characteristics of both the block soil samples and the comparison soil samples, an ASD FieldSpec3 hand-held spectrometer produced by Analytical Spectral Devices Inc. was used in this experiment, which was widely used in the spectral measurements of vegetation and soil. This spectrometer covers wavelengths between 350-2500 nm with two separate detectors namely VNIR (visible and near-infrared, 350-1000 nm) and SWIR (short-wave infrared, 1000-2500 nm), it can savethe spectral data and display the spectral curve of the objects. The sampling intervals are 1.4 nm with a spectral resolution of 3 nm and 2 nm for VNIR, and with a spectral resolution of 10 nm SWIR, respectively. All spectral reflectance was measured in clear and cloudless weather conditions.
For both block and comparison soil samples, the spectrometer with an 8° field of view was held on a fixed platform 10 cm vertically above the surface of soil samples to ensure that the zenith angle and the azimuth angle were kept the same, the corresponding measurement area is then determined as a circle area with a very small radius of 0.7 cm. After that, each soil sample was put uniformly right on a 1m×1m blackboard on the ground in order to minimize the effects of the background and the adjacent objects. 10 randomly spots on the surface of soil samples were then chosen as the center, and the spectral measurements were carried out. Note that before spectral measurements, all comparison soil samples were firstly put in the same aluminum boxes with an inner diameter of 10 cm and a depth of 2 cm (1.5 cm of soil can be considered as infinite thickness in optical remote sensing) and then flattened with a spatula to remove the effect of compaction on the spectral reflectance. During spectral measurement of each soil sample, 10 readings of spectral reflectance were taken, and the mean of all spectral reflectance was calculated as the actual reflectance in order to average out the difference caused by instrument noise. A standard white panel was also measured as calibration before each spectral measurement. In order to diminish the noise of the equipment and further reduce the data dimensions, the spectral reflectance was resampled at 10 nm for all the block and comparison soil samples. For both block and comparison soil samples, the spectrometer with an 8 • field of view was held on a fixed platform 10 cm vertically above the surface of soil samples to ensure that the zenith angle and the azimuth angle were kept the same, the corresponding measurement area is then determined as a circle area with a very small radius of 0.7 cm. After that, each soil sample was put uniformly right on a 1m×1m blackboard on the ground in order to minimize the effects of the background and the adjacent objects. 10 randomly spots on the surface of soil samples were then chosen as the center, and the spectral measurements were carried out. Note that before spectral measurements, all comparison soil samples were firstly put in the same aluminum boxes with an inner diameter of 10 cm and a depth of 2 cm (1.5 cm of soil can be considered as infinite thickness in optical remote sensing) and then flattened with a spatula to remove the effect of compaction on the spectral reflectance. During spectral measurement of each soil sample, 10 readings of spectral reflectance were taken, and the mean of all spectral reflectance was calculated as the actual reflectance in order to average out the difference caused by instrument noise. A standard white panel was also measured as calibration before each spectral measurement. In order to diminish the noise of the equipment and further reduce the data dimensions, the spectral reflectance was resampled at 10 nm for all the block and comparison soil samples. Table 1 lists the statistical characteristics of the soil properties of soil samples. The total salinity is the sum of ion content. Na + is the main cation while HCO 3 − and CO 3 2− are the main anions. the sample point selection accurately reflects the distribution range of the salt parameters. Furthermore, the differences of particle size distributions of the soil samples are very small especially for those particles with clay content <0.005mm, the CV and standard are 5.49% and 1.54 mg/g, which means no significant differences were observed in the clay content for the soil samples.  Figure 4 shows the distributions of the CON texture feature calculated from the crack images of soil samples, which indicates a very large heterogeneity among different soil samples with CON ranges from 0.007 to 0.342.

Spectral Characteristics
In order to study the spectral response of soda saline-alkali soil to salt contents of both block soil samples and comparison soil samples, sample 29, sample 8, sample 35, sample 43, sample 28 and sample 21 were selected according to the distribution of total salinity at an interval of 5 mg/g, the corresponding salinities were 5.17 mg/g, 10.37mg/g, 15.47mg/g, 20.71mg/g, 25.62mg/g, 29.73 mg/g respectively. The reflectance curves of the above soil samples are shown in Figure5. From Figure  5a,b, it can be seen that the overall reflectance curves of both block soil samples and comparison soil samples are relatively low in the visible region. The spectral reflectance increases from 350 to 1400 nm and turns basically flat between 1400 and 2200 nm, and then the reflectance curve gradually

Spectral Characteristics
In order to study the spectral response of soda saline-alkali soil to salt contents of both block soil samples and comparison soil samples, sample 29, sample 8, sample 35, sample 43, sample 28 and sample 21 were selected according to the distribution of total salinity at an interval of 5 mg/g, the corresponding salinities were 5.17 mg/g, 10.37mg/g, 15.47mg/g, 20.71mg/g, 25.62mg/g, 29.73 mg/g respectively. The reflectance curves of the above soil samples are shown in Figure 5. From Figure 5a,b, it can be seen that the overall reflectance curves of both block soil samples and comparison soil samples are relatively low in the visible region. The spectral reflectance increases from 350 to 1400 nm and turns basically flat between 1400 and 2200 nm, and then the reflectance curve gradually decreases with the increasing salinity. Figure 5a,b also indicates that the block soil samples are more capable distinguishing the soil salinity compared with the comparison soil samples, the higher the salinity, the lower the reflectance curve. In the sensitive wavelengths above 1400 nm, the spectral curve difference of the block soil samples becomes smaller, and the reflectance interval is more uniform. In addition, Figure 5c,d illustrates the difference in spectra between the block soil sample and the comparison soil sample at the same salinity, it can also be seen that the reflectance curve of the block soil sample is significantly lower than that of the comparison soil sample for both the highest salinityand the lowest salinity.

Correlation Analysis
In order to extract the optimum reflectance and the corresponding diagnostic wavelengths for developing prediction models with high accuracy, correlograms were calculated in this study across seven soil properties of particular interest (four main ions including Na + , Cl − , HCO3 − , CO3 2− , total salinity, pH and EC) for reflectance of soil samples under different soil surface conditions. The correlograms can display absolute value curves of the correlation coefficient between a given soil property and spectral reflectance as a function of wavelength (Figure6).
Figure6a indicates that for the block soil samples, the shapes of correlation coefficient curves for EC, Na + , Cl − , and salinity are basically the same, with an upward trend from 350 to 1400 nm and a small peak near 500 nm. The correlation between soil parameters and reflectivity decreases firstly and then increases between 1400 nm and 1900 nm, but the change is not obvious. The curve remains basically flat with a small uplift around 1900 nm, and then slightly decreases. The maximum value of the correlation coefficient appears at 1990 nm. For CO3 2− , the shape of the correlation coefficient curve is nearly the same as EC, Na + , Cl − and salinity from 350 to 1400 nm, but appears opposite trend from 1400 to 2500 nm. The curve shapes of HCO3 − and pH are basically the same, show a downward trend below 1400 nm and a small absorption valley around 1000 nm. The correlation coefficient curves of HCO3 − and pH show a trend of decreasing first and then increasing from 1400 to 2500 nm. The correlation between pH and reflectance is significantly lower than that of HCO3 − . This is because

Correlation Analysis
In order to extract the optimum reflectance and the corresponding diagnostic wavelengths for developing prediction models with high accuracy, correlograms were calculated in this study across seven soil properties of particular interest (four main ions including Na + , Cl − , HCO 3 − , CO 3 2− , total salinity, pH and EC) for reflectance of soil samples under different soil surface conditions. The correlograms can display absolute value curves of the correlation coefficient between a given soil property and spectral reflectance as a function of wavelength ( Figure 6). Figure 6a indicates that for the block soil samples, the shapes of correlation coefficient curves for EC, Na + , Cl − , and salinity are basically the same, with an upward trend from 350 to 1400 nm and a small peak near 500 nm. The correlation between soil parameters and reflectivity decreases firstly and then increases between 1400 nm and 1900 nm, but the change is not obvious. The curve remains basically flat with a small uplift around 1900 nm, and then slightly decreases. The maximum value of the correlation coefficient appears at 1990 nm. For CO 3 2− , the shape of the correlation coefficient curve is nearly the same as EC, Na + , Cl − and salinity from 350 to 1400 nm, but appears opposite trend from 1400 to 2500 nm. The curve shapes of HCO 3 − and pH are basically the same, show a Remote Sens. 2019, 11, 1406 9 of 18 downward trend below 1400 nm and a small absorption valley around 1000 nm. The correlation coefficient curves of HCO 3 − and pH show a trend of decreasing first and then increasing from 1400 to 2500 nm. The correlation between pH and reflectance is significantly lower than that of HCO 3 − . This is because the pH is determined by the OH − content produced by the CO 3 2− and HCO 3 − hydrolysis reaction, which is often incomplete and reversible due to conditions including temperature and water content. Therefore, in the sensitive band above 1400 nm, the correlation between HCO 3 − and pH and the spectral reflectance is quite poor, but the overall trend of correlograms is still the same as that of salinity. Figure 6b shows that the correlograms of comparison soil samples are opposite to those of the block soil samples in the visible range from 350 to 650 nm. The correlation coefficient curves of comparison soil samples share the same trends with the block soil samples at the wavelength above 650 nm but with rougher shapes and greater magnitudes. In addition, Figure 6c,d illustrates the difference in coefficient between the block soil samples and the comparison soil samples for salinity and EC, it can be seen that the correlogram curves of the block soil samplesaresignificantly higher than those of the comparison soil samples for both salinity and EC.
Remote Sens. 2019, 11, x FOR PEER REVIEW 10 of 19 difference in coefficient between the block soil samples and the comparison soil samples for salinity and EC, it can be seen that the correlogram curves of the block soil samplesaresignificantly higher than those of the comparison soil samples for both salinity and EC.

Regression Models
According to Figure 6, the most relevant wavelength was selected at 1990 nm as the diagnostic band, the reflectance of all the block soil samples, and the comparison soil samples was then extracted respectively. Mixed reflectance at 1990 nm was also computed considering the CON texture feature, which was calculated for both the block soil samples and the comparison soil samples based on Equation (6).

Regression Models
According to Figure 6, the most relevant wavelength was selected at 1990 nm as the diagnostic band, the reflectance of all the block soil samples, and the comparison soil samples was then extracted respectively. Mixed reflectance at 1990 nm was also computed considering the CON texture feature, which was calculated for both the block soil samples and the comparison soil samples based on Equation (6).
where R(λ) represents the mixed reflectance, R 0 (λ) represents the reflectance of crack regions (ignored in this study), R 1 (λ) represents the original reflectance of pure soil surface corresponding to the block and the comparison soil samples. In order to effectively predict the soil properties of soda saline-alkali soils and take into account the sufficient reduction of the work intensity, four main soil parameters of interest including Na+, EC, pH, and salinity were selected, both the 57 block soil samples and the 57 comparison soil samples were further arranged according to salinity and evenly selected 38 of them as calibration samples for regression models and other 19 were used for validation. After that, univariate linear regression models were established using both the original reflectance and the mixed reflectance computed from Equation (6), and the results are shown in Table 2. It can be seen from Table 2 that the modeling accuracy using the original reflectance of comparison soil samples is lowest with R 2 are only 0.41, 0.39, 0.10 and 0.45 for Na + , EC, pH, and salinity, that the modeling accuracy using the original reflectance of block soil samples is higher with R 2 are 0.74, 0.68, 0.23, 0.72 for Na + , EC, pH and salinity. Moreover, the modeling accuracy of the univariate regression model is greatly improved for both the block and the comparison soil samples after considering the CON texture feature. In order to make full use of the spectral reflectance, 1990 nm, and another three sensitive bands of 990 nm, 1470 nm, and 1990 nm were selected as the characteristic spectrum bands. Multivariate regression models of soil properties were then developed based on both the original reflectance and the mixed reflectance computed from Equation (6) at above four sensitive wavelengths for the block soil samples and the comparison soil samples respectively (Table 3). From Table 3, it can be seen that the modeling accuracy of multivariate linear models for soil properties is higher than that of univariate linear models. Table 3 also indicates that the accuracy based on the original reflectance of comparison soil samples is the poorest with R 2 of Na + , EC, pH, and salinity are only 0.44, 0.42, 0.17 and 0.49, that the modeling accuracy using the original block soil samples is higher with R 2 are 0.75, 0.70, 0.29, 0.75 for Na + , EC, pH and salinity. In addition, the CON texture feature can also help to improve the modeling accuracy of multivariate regression models for both the comparison soil samples and the block soil samples.

Discussion
Many studies have shown that the cracking process of clayey soil is mainly affected by the mechanical properties such as the soil matrix suction, the tensile strength, and the surface energy, these factors generally depend on the mineral type and clay content [45][46][47]. Du et al. [48] and Zhang and Peng [49] have shown that soil clay is an important factor affecting soil shrinkage and cracking process due to its obvious plasticity, swelling, and cohesiveness. Gray and Allbrook [50] found that soil shrinkage and cracking ability is significantly positively correlated with soil clay content. Greene-Kelly [51] showed that soil shrinkage and crack ability is stronger when the expansive clay minerals (such as smectite) dominate in the soil. After measuring the relationship between internal friction angle and cohesion and smectite content, Chen et al. [52] illustrated that when the smectite content is above 5%, the internal friction angle and soil cohesion of soil particles begin to decrease while shrinkage cracking enhances. Zhang et al. [53] carried out X-ray diffraction complete analysis for soil minerals, their measurement results showed that the mineral composition of the saline soils in Songnen Plain is dominated by quartz and feldspar, which accounts for 87.6% to 91.6% of mineral contents, and that only a small amount of kaolinite and almost no highly active smectite is contained. Illite/smectite mixed layer formation is the secondary mineral with an interlaying ratio above 0.5 and the activity index only from 0.33 to 0.48, indicating that the soda saline-alkali soil in Songnen Plain has relatively high compressibility, poor water permeability, low shear strength and limited expansion [54]. From Table 1, it can be seen that although relatively high clay contents (26.28% to 29.37%) were measured in this study, the Std and CV are only 1.54% and 5.49%, this indicates that the clay contents of soil samples are almost the same in this study. On the other hand, the soil sample preparation process was kept the same, and the cracking test was carried out under the same controlled conditions (temperature of 25 • C, humidity of 35% and pressure of 101 kPa), indicating that laboratory conditions also had no effect on the cracking process in this study. Table 4 shows that the crack degree has a good correlation with the soil properties related to the salt content, the correlation coefficients between the soil properties and the CON are all above 0.68. This is because during the dehydration process of the soil samples, a bounded water film is forming between soil particles (especially Na+ with the large hydrolysis radius) due to the interaction of soil particles with exchangeable cations. The bounded water film reduces the soil cohesion and the tensile strength [55], which also has a great influence on the internal friction angle and the shear strength of the soil [56][57][58]. In addition, the well-known diffuse double layer (DDL) thickness can also be considered as an important role in soil surface cracking, which always becomes thinner as the salt contents of soil samples increase and directly enhances the shrinkage and desiccation cracking of soils during water evaporation [59][60][61]. It can be seen from Figure 5 that surface states of soil samples cannot affect the morphology of the reflectance curves but can largely change the amplitudes of the curves. The reflectance curves of the comparison soil samples are higher than those of the block soil samples. This is because after the desiccation cracking process, the area scattering on the surface of block soil samples strengthens with the complexity of surface morphology characterized by the CON texture feature, and the relative surface roughness of block soil samples increases, which reduces the energy reflected to the spectrometer, and thus enhances the spectral difference between different block soil samples. In addition, due to the desiccation cracking test was carried out under natural conditions, all block soil sampleswere air-dried with some residual water remainingin the block soil samples and, therefore, reduced the overall reflectivity compared with the comparison soil samplestotally driedinan oven.Moreover, the shadows appeared within crack regions and made the energy reflect back to the spectrometer close to zero since the incident sunlight was not perpendicular to the surface of the soil samples.
A gradual decrease in reflectance with increasing salt contents was evident, as shown in Figure 5, especially for the wavelength above 1400 nm since NaHCO 3 and Na 2 CO 3 dominants the salt minerals in Soda saline-alkali soils. Howari et al. [62] and Farifteh et al. [63,64] showed that due to the unique resonance caused by stretching and bending, the sharp decline in reflectance of NaHCO 3 and Na 2 CO 3 as the wavelength increases can be considered diagnostic since no spectra of other salt minerals were investigated. Similar finding was reported by Wang et al. [65] that the Na 2 CO 3 -affected soil has the lowest reflectance shape comparing with saline soils affected by Na 2 SO 4 , NaCl and the soil without any salt, their results also indicate that a clear decreasing pattern of reflectance with increasing content of Na 2 CO 3 can be distinguished especially for the wavelength above 1400 nm, which is quite in agreement with the results in this study. Although the residual water contents are very limited (ranging from 0.18% to 2.34% with variance only 0.03%, indicated in Table 5), hygroscopic water can still generate spectral response difference among the block soil samples to a certain extent. This is because from the results studied by Huang et al. [66], salt content can increase the water-holding capacity of soil and thus further reduce the spectral reflectance of soil samples with higher salinities. In addition, according to Tedeschi and Dell' Aquila [67] and Huang et al., the soil aggregate ability decreases with increasing salinity as a result of the deflocculation of clay colloids induced by the salt, which maybe another reason causing the salinity inversely proportional to the spectral reflectance of the soil samples. Figure 5 also shows that the soil samples under different surface conditions have the same absorption bands and similar corresponding absorption characteristics. Due to the influence of water content in the atmosphere, the reflectance curves exhibit obvious absorption feature at 1900 nm. The reflectance curves also exhibit a weak absorption characteristics of NaCl at 1190 nm, 1150 nm and three weak absorption valleys of NaHCO 3 at 1470 nm, 1990 nm and 2170 nm, these weak absorption features are caused by the resonance by the stretching and bending between the ions of the soil salt minerals. The water evaporation during the drying process caused the salt ions to precipitate on the surface of soil samples, making the spectral response to the main soil properties of the block soil samples quite high. The regression results given in Tables 2 and 3 indicate that when considering the CON texture feature of crack images, the modeling accuracies of main soil properties can be greatly improved compared with those models developed based on the original reflectance. This is due to the fact that the desiccation cracking process of the soil samples can make the salt separate out on the surface of soil samples, which makes the spectral response of the block soil samples to soil properties stronger than that of the comparison soil samples. Since the reflectance of the crack region is close to zero and the cracking degree is proportional to the salinity of soil samples, a more spectral difference of the block soil samples can be found, and the estimation accuracy of soil properties is, therefore, improved.
In order to quantitatively evaluate the effectiveness of the regression models considering the CON texture feature, the measured soil properties and the estimated soil properties of all the validation soil samples were linearly fitted to a function of y = x shown in Figure 7. In addition, the root-mean-square error (RMSE), the relative root-mean-square error (RMSE%), the mean absolute error (MAE), the relative mean absolute error (MAE%) and the predicted ratio of deviation (PRD) were calculated as the evaluation indexes, and the results are shown in Table 6. According to the criteria given by Farifteh et al. [68], both univariate and multivariate regression models developed from the block soil samples provide good accuracy for EC, salinity and Na + (0.82 < R 2 < 0.9 and RPD>2), the calculated RMSE, RMSE%, MAE, and MAE% also show that the regression models of EC, salinity and Na + have quite high reliability and stability. In addition, the univariate and multivariate prediction models developed from the comparison soil samples also show an acceptable prediction accuracy for EC, salinity and Na + (0.5 < R 2 < 0.65 and RPD > 1.5), but high RMSE, RMSE%, MAE The MAE% indicate that the reliability and stability are poor. However, although CON texture features considered, both univariate and multivariate regression models are inaccurate for pH, whether based on the block soil samples or the comparison soil samples. prediction models developed from the comparison soil samples also show an acceptable prediction accuracy for EC, salinity and Na + (0.5<R 2 <0.65 and RPD>1.5), but high RMSE, RMSE%, MAE The MAE% indicate that the reliability and stability are poor. However, although CON texture features considered, both univariate and multivariate regression models are inaccurate for pH, whether based on the block soil samples or the comparison soil samples.  For the soda saline-alkali soil with surface cracks under natural conditions, the regression models developed in this study can be well implemented to estimate the main soil properties such as EC, salinity and the cation of Na + with relatively high accuracy and less work intensity. In detail, only soil samples need to be photographed and then collected firstly in the field. Secondly, the spectral measurements can be carried out based on the obtained soil samples (the soil samples collected may be soil blocks or soil powders, corresponding to the block soil samples and the comparison soil samples mentioned in this study) under laboratory-controlled conditions, which usually takes less time than the spectral measurement under natural world. Thirdly, the CON texture features of soil samples with surface cracks can be calculated individually from the crack images obtained in the field. After that, the developed regression models in this study can, therefore, be used for soil properties prediction using the mixed spectral reflectance of the soil samples computed from Equation (6) in Section 3.5. However, the method and the regression linear models developed in this study are still limited for only the salt-affected soil with relatively high clay content, which generates surface cracks commonly under natural conditions. Moreover, other anticipated problems, such as the height of the photograph, the CON scale, the humidity effects, and the atmosphere attenuation, also need to be considered in the field measurement.

Conclusions
In this study, a cracking experiment was carried out on 57 samples of soda saline-alkali soil in the Songnen Plain of China. A procedure was studied to improve soil property estimation of soda saline-alkali soil based on spectral measurement considering the CON texture feature of the soil surface with cracks. The spectral response of the block soil samples separated by crack regions and the 2 mm comparison soil samples with different soil properties was analyzed and compared. After that, univariate and multivariate linear regression models were developed for four main soil properties of interest including Na + , EC, pH and total salinity based on the original reflectance and the mixed reflectance considering the CON texture feature for the block soil samples and the comparison soil samples respectively. The main conclusions in this investigation are as follows: 1.
The soil salinity can dominate the cracking process of soda saline-alkali soil to a certain extent for the soil with similar clay content and clay mineral composition.

2.
The spectral characteristics of the block soil samples are more correlated with soil properties than the comparison soil samples, the mixed reflectance of soil samples considering the CON texture feature can largely improve the prediction accuracy of the main soil properties.

3.
Linear regression models based on mixed reflectance considering the CON texture feature of block soil samples show good potential for predicting the main soil properties.