Estimating soil profile salinity under vegetation cover based on UAV multi-source remote sensing

Estimating soil profile salinity under vegetation cover based on UAV multi-source remote sensing | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Article Estimating soil profile salinity under vegetation cover based on UAV multi-source remote sensing Zhenhai Luo, Meihua Deng, Min Tang, Rui Liu, Shaoyuan Feng, Chao Zhang, and 1 more This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-4971758/v1 This work is licensed under a CC BY 4.0 License Status: Published Journal Publication published 21 Jan, 2025 Read the published version in Scientific Reports → Version 1 posted 11 You are reading this latest preprint version Abstract Soil salinization is the most common land degradation problem in arid, semi-arid and coastal areas of China, which seriously affects local crop yield, economic development, and environmental sustainability. There are few studies on estimating soil salinity at different depths under vegetation cover. In this study, field soil control experiments were employed to collect multi-source remote sensing data under barley growth, and soil salt content (SSC) with various depths. Three types of feature variables were built based on images and were filtered by the boosting decision tree (BDT) method. Besides, four machine learning algorithms coupling with seven variable combination groups were used to comprehensively establish soil salinity estimation model. Finally, the performances of estimation model for different crop over ratios were evaluated. The results showed that the gaussian process regression (GPR) model based on the full variable group at the depths of 0 ~ 10 cm and 30 ~ 40 cm is more accurate than other models. The validation R 2 is 0.774 and 0.705, and the RMSE is 0.185% and 0.31%;The random forest (RF) models based on spectral index and texture data at 10 ~ 20 cm and 20 ~ 30 cm depths are more accurate, with R 2 of 0.666 and 0.714. SSC may be quantitatively inverted at various depths using the machine learning model based on multi-source remote sensing, which also serves as a guide for monitoring soil salinization. Earth and environmental sciences/Biogeochemistry Earth and environmental sciences/Environmental sciences Physical sciences/Engineering soil salt content UAV vegetation cover soil depth multi-source remote sensing data Figures Figure 1 Figure 2 Figure 3 Figure 4 Figure 5 Figure 6 Figure 7 1. Introduction In China, the soil salinization area extends 9.21 million hectares, accounting for 6.62% of the total cultivated land. The process of saline soil formation is intricate, posing significant challenging for the detection and dynamic monitoring of soil salinization. Conventionally, soil salinization was assessed through field sampling and chemical analysis, these methods that were both time consuming, and labor-intensive 1 . In contrast, remote sensing technology, a rapid and expansive means to gather data on ground objects across different temporal and spatial, is an ideal tool for soil salinity monitoring 2 . The spectral response of soils varied with salt content, with high-salinity soils exhibiting stronger responses in the visible and near-infrared bands than those with lower salinity 3 . Leveraging remote sensing for the dynamic monitoring of soil salinization is crucial for the efficient utilization of soil and water resources, providing essential insights into the timing, patterns, and locations of potential changes in soil salinity, thus facilitating improved resource management and planning. In recent years, Unmanned Aerial Vehicle (UAV) and other aerial remote sensing platforms have experienced rapid advancements, increasingly integrating into civilian applications and becoming prominent in agricultural research. UAVs boast advantages like portability, high flexibility, and customizable flight durations. Zhao et al. 4 used multispectral remote sensing data from three research locations to establish soil salinity inversion models based on support vector machine (SVM), RF, backpropagation neural network (BPNN), and extreme learning machine (ELM). The results showed that the four spectral index-based models could get a high level of inversion accuracy. Similarly, Wei et al. 5 used a UAV equipped with Micro-MCA multispectral sensors to capture pictures for evaluating the soil salinity in a small area of the Hetao Irrigation District. Chen et al. 6 developed estimation models of SSC in sunflower fields at different soil depths during the budding and blooming stages using UAV multispectral data. IVUSHKIN et al. 7 found that UAVs equipped with multiple sensors of hyperspectral, multispectral, thermal infrared, and LiDAR cameras have great potential for monitoring soil salinization. Feature indices like the soil salinity index and vegetation index, derived from spectral band reflectance transformations, serve as important variables for estimating soil salinity content 8 . Qi et al. 9 collected spectrum reflectance and spectral indices based on a UAV ground cooperative system, and then used machine learning algorithms like BPNN to build a salinity inversion model. The findings indicated that the constructed model accurately captured the level of salinization in the research area. The relationship between soil salinity and remote sensing feature variables was frequently nonlinear due to the combination of complex factors such as soil, vegetation, and atmospheric signatures 10 . A common method for estimating SSC was to use the mathematical statistical model, especially the linear regression model including partial least squares regression (PLSR) 11 , 12 . But the relationship between spectral covariates and soil properties in nature was rarely linear 13 . Machine learning (ML) algorithms, renowned for their ability to address nonlinear problems and process high dimensional data, often outperform statistical regression models in predicting soil salinity. Algorithms like RF、SVM and BPNN could capture these intricate nonlinear relationships through multi-layered results, making them particularly effective for soil salinity estimation 14 . Hu 15 compared PLSR and RF methods using hyperspectral first order differentiation, broad band spectral indices, and narrow band spectral indices as independent variables. It was reported that the RF model could better predict soil salinity and the model for the bare soil area had the highest prediction accuracy compared to the other sample areas. Zhang et al. 16 found that the SVM method based on texture features had an improved effect on the monitoring accuracy of soil salinity information. Wei et al. 5 established various machine learning models to evaluate and optimize the best salt estimate model. However, the prediction accuracy of a single machine learning algorithm could fluctuate under varying conditions. Therefore, evaluating the performance of multiple machine learning regression algorithms is essential for constructing reliable soil salinity prediction models that could adapt to diverse environmental factors. There is a strong correlation between vegetation growth and soil salinity, as evidenced spectrally in two primary ways: differences in leaf spectral reflectance and significant variations in the texture features of spectral images influenced by changes in soil salinity or leaf characteristics 17 , 18 Studies have shown that texture features are extensively utilized to reveal vegetation characteristics variations, combining them with spectral information could effectively improve the accuracy of predictive models 19 , 20 , 21 , 22 , 23 . Huang et al. 24 used Sentinel-2 image combined with texture features to significantly improve the classification effect of moderate saline soil in the Yellow River Delta region. Nevertheless, the prevalent use of remote sensing data for soil salinity detection primarily relied on spectral indices 5 , with limited studies investigating the use of texture features in soil salinity estimation 24 . In the case of bare soil, the spectrum could directly determine the salt content of the soil surface. Conversely under vegetation coverage, soil salinity could be indirectly assessed by vegetation canopy spectral signatures 25 . Few investigations on estimating soil salinity under vegetation cover have been conducted in the past. Therefore, the purpose of this study is to construct soil salinity estimation models at various depths beneath barley growth. The main research contents were to:(1) evaluate the potential and feasibility of UAV remote sensing for SSC estimation under barley coverage; (2) optimize the feature variables of spectral band, spectral index and texture data based on BDT method; (3) Validate the accuracy of soil salt estimation model for different vegetation coverage conditions. This study presented a novel method for the dynamic monitoring of soil salinity in agricultural production helped to achieving precise irrigation and fertilization practices, and improved agricultural productivity efficiency and sustainable utilization of soil resources. 2. Materials and methods 2.1 Study area and experimental design The study area is located at the ecological experimental station of Yangzhou University, situated in the Jianghuai Plain of Jiangsu Province in eastern China (119°24′E, 32°21′N), at an altitude of 5 meters. This region is characterized by a subtropical monsoon climate, with an average annual frostfree period, precipitation, evaporation, and air temperature of 223 days 937 mm, 1063 mm, and 14.8°C, respectively. Soil samples for the experiment were collected from Tiaozini reclamation area in Dongtai City, Jiangsu Province. This area was subject to marine intrusion and groundwater topdressing, with high salinity in the soil tillage layer. The crop studied was barley (Hordeum vulgare L.). Four different soil salinity treatments were set: control (no salt), low salinity (3‰), medium salinity (5‰), and high salinity (10‰). The salinity experiment was carried out using a box planting setup with dimensions of 100 cm in length, 40 cm in width, and 40 cm in height, each containing 120 kg of base soil. Each treatment was replicated thrice, resulting in a total of 12 experimental units. The cultivation adhered to local management practices, encompassing weeding, pest, and disease prevention and control. Barley was planted annually at the early of November, with each box accommodating two rows and was harvested at the late of May of the following year. To maintain the soil salinity status of each treatment, the experimental box devices were designed as sealed containers to prevent salt leaching by precipitation and irrigation. 2.2 Data collection and acquisition 2.2.1 Soil salt measurement During the barley growth stages of reviving-jointing, jointing-filling, and grain filling maturity, soil electrical conductivity data were measured every 7–10 days, in conjunction with the acquisition of multispectral imagery. A total of 84 dataset were collected during barley growth. The soil electrical conductivity was measured using the EC450 conductivity meter (Spectrum Technologies Co., Ltd., Chicago, IL, USA). First, the electrode was calibrated using a calibration solution (conductivity: 1413 µS/cm). After calibration, the electrode was inserted into the soil profile to measure conductivity at soil depths of 0ཞ10 cm, 10ཞ20 cm, 20ཞ30 cm, and 30ཞ40 cm. The soil conductivity values were directly recorded by the handheld meter. Subsequently, the SSC was derived using the empirical formula (SSC = 0.2882EC + 0.0183, %). 2.2.2 Multispectral data acquisition and processing This study leveraged remote sensing data from multiple UAV based sensors to enhance monitoring of soil properties through data fusion. The DJI Inspire 2 UAV platform (DJI Inc., Shenzhen, China) used in this study, which was equipped with a Altum multi-spectral and infrared camera (MicaSense, Inc., Seattle, WA, USA). This camera was capable of capturing images across six spectral bands (blue, green, red, red edge, near-infrared, and thermal infrared) simultaneously. Variables The UAV flight operations were conducted under optimal conditions of clear skies and calm winds, between 11 a.m. and 2 p.m. local time. To ensure high-quality data, the sensor was prewarmed for 5 minutes and the reference plate was used for radiometric calibration before each flight. The flight altitude and the cross-track overlap were maintained to 25 meters and 75%, respectively. The camera was oriented vertically downward to achieve a ground resolution of 1.1 cm per pixel. After preprocessing the collected multispectral images by radiation correction, geometric correction, and image mosaic, the reflectance data of each pixel were generated during crop growth period. To extract canopy reflectance, a region of interest (ROI) was preset near the center of each plot. 2.2.3 Vegetation coverage calculation After collecting multispectral imagery, the canopy light interception of barley was measured using the AccuPAR LP-80 ceptometer (Decagon Devices Inc, Pullman, WA, USA). The LP-80 is equipped with 80 independent sensors, each spaced uniformly at 1 cm intervals, capable of measuring solar radiation within the 400–700 nm band in different modes. During the measurement, the sensor was strategically positioned in the center of the plot, aligned with the row direction. With the acquired photosynthetically active radiation at the top and bottom of the canopy, the instrument automatically calculated the crop LAI through a built-in algorithm with the other variables. In accordance with the reports and employing a classification method that considers similar geographical features and vegetation types 26 , vegetation cover was classified according to the LAI value as follows: low coverage (0, 0.45), medium coverage (0.45, 0.75) and high coverage (0.75, 1). 2.3 Construction of feature variables Spectral indices were classical variables that integrated the spectral characteristics of each band of ground objects and enhanced specific information through mathematical transformations and combinations of reflectance values from different bands. Salinity indices, which were frequently employed for the rapid assessment of soil salinization, exhibited a strong correlation with bare soil salinity 27 . Similarly, vegetation indices were frequently used for the quantitative assessment of vegetation growth 28 . In this study, a selection of widely recognized spectral indicators were employed for the monitoring of soil salinization, including 15 salinity indices and 15 vegetation indices. The selected salinity indices include: Normalized Difference Salinity Index (NDSI), R-edge Normalized Difference Salinity Index (NDSI-reg), Brightness Index (BI), Salinity Index 1 (SI1), R-edge Salinity Index 1 (SI1-reg), Salinity Index 2 (SI2), R-edge Salinity Index 2 (SI2-reg), Salinity Index 3 (SI3), R-edge Salinity Index 3 (SI3-reg), Salinity Index (SI-T), Salinity Index S1, Salinity Index S2, Salinity Index S3, Salinity Index S5, Salinity Index SI; The selected vegetation indices include: Normalized Difference Vegetation Index (NDVI), R-edge Normalized Difference Vegetation Index (NDVI-reg), Difference Vegetation Index (DVI), R-edge Difference Vegetation Index (DVI-reg), Enhanced Vegetation Index (EVI), R-edge Difference Vegetation Index (EVI-reg), Triangular Vegetation Index (TVI), Normalized Greenness Index (NDGI), Simple Ratio Index (SR), Modified Soil Adjusted Vegetation Index (MSAVI), Optimized Soil Adjusted Vegetation Index (OSAVI), Soil Adjusted Vegetation Index (SAVI), Visible Light Band Difference Vegetation Index (VDVI), Visible Atmospherically Resistant Index (VARI), Green Normalized Difference Vegetation Index (GNDVI), and a normalized relative canopy temperature (NRCT). The calculation formulas for these indices were shown in Table 1 . Table 1 Spectral index and calculation formula Feature Formulation Reference NDSI NDSI= \(\:\frac{R-NIR}{R+NIR}\) 29 NDSI-reg NDSI-reg= \(\:\frac{RedEdge-NIR}{RedEdge+NIR}\) 30 BI BI= \(\:\sqrt{{R}^{2}+{NIR}^{2}}\) 31 SI1 SI1= \(\:\sqrt{G\times\:R}\) 28,32 SI1-reg SI1-reg= \(\:\sqrt{G\times\:RedEdge}\) 30 SI2 SI2= \(\:\sqrt{{G}^{2}+{R}^{2}+{NIR}^{2}}\) 28,32 SI2-reg SI2-reg= \(\:\sqrt{{G}^{2}+{RedEdge}^{2}+{NIR}^{2}}\) 30 SI3 SI3= \(\:\sqrt{{G}^{2}+{R}^{2}}\) 28,32 SI3-reg SI3-reg= \(\:\sqrt{{G}^{2}+{RedEdge}^{2}}\) 30 SI-T SI-T= \(\:100\times\:R/NIR\) 33 S1 S1= \(\:B/R\) 34 S2 S2= \(\:\frac{B-R}{B+R}\) 34 S3 S3= \(\:G\times\:R/B\) 34 S5 S5= \(\:B\times\:R/G\) 34 SI SI= \(\:\sqrt{B\times\:R}\) 35 NDVI NDVI= \(\:\frac{NIR-R}{NIR+R}\) 36 NDVI-reg NDVI-reg= \(\:\frac{NIR-RedEdge}{NIR+RedEdge}\) 30 DVI DVI= \(\:NIR-R\) 28 DVI-reg DVI-reg= \(\:NIR-RedEdge\) 30 EVI EVI= \(\:2.5\times\:\frac{NIR-R}{NIR+6R-7.5B+1}\) 37 EVI-reg EVI-reg= \(\:2.5\times\:\frac{NIR-RedEdge}{NIR+6RedEdge-7.5B+1}\) 30 TVI TVI= \(\:0.5\times\:\frac{120\times\:\left(NIR-G\right)}{200\times\:\:\left(R-G\right)}\) 9 NDGI NDGI= \(\:\frac{G-R}{G+R}\) 38 SR SR= \(\:NIR/R\) 39 MSAVI MSAVI= \(\:\frac{\left(2NIR-1\right)-\sqrt{{\left(2NIR+1\right)}^{2}-8\left(NIR-R\right)}}{2}\) 30 OSAVI OSAVI= \(\:1.16\times\:\frac{NIR-R}{NIR+R+0.16}\) 30 SAVI SAVI= \(\:1.5\times\:\frac{NIR-R}{NIR+R+0.5}\) 40 VDVI VDVI= \(\:\frac{2G-\left(R+B\right)}{2G+\left(R+B\right)}\) 41 VARI VARI= \(\:\frac{G-R}{G+R+B}\) 42 GNDVI GNDVI= \(\:\frac{NIR-G}{NIR+G}\) 9 NRCT NRCT= \(\:\frac{Ti-Tmin}{Tmax-Tmin}\) 43 Note: B, R, G, NIR, and Rededge represent the reflectance of blue, red, green, near-infrared, and red edge bands, respectively. The image texture was indicative of variations in soil surface’s color and gray level, which was closely related to the soil's salt and water content. In this study, the statistical Gray Level Cooccurrence Matrix (GLCM) method was used to extract the textural features of the images. GLCM is a prevalent method widely used for image feature extraction, texture analysis, and quality evaluation 44 , describing the correlation between pixel gray levels within images. Eight characteristic variables including mean (MEA), variance (VAR), uniformity (HOM), contrast (CON), difference (DIS), entropy (ENT), second-order moment (SEC) and correlation (COR) were calculated for each band of UAV multispectral images utilizing the second-order statistical filtering tool. Variables In total, this study incorporated 5 band reflectance values, canopy temperature, 31 spectral indices, and 40 texture feature data, amounting to 77 feature variables serving as independent variables. 2.4 Feature variable selection and set construction 2.4.1 Optimization of characteristic variables To optimize the complexity of the model input variables, this study employed the BDT method to refine the selection of the 77 characteristic variables. BDT is one of the embedded type methods which used to compute estimates of Predictor Importance for the tree by summing changes in the mean squared error (MSE) due to splits on every predictor and dividing the sum by the number of branch nodes 45 . The selection steps involved in this optimization process include: (1) the data were divided into twenty groups; (2) importance score of feature variables in each set of data was calculated based on BDT; (3) normalizing the important scores for various data categories. (4) A statistical analysis and sorting were performed on the normalized data of 20 groups. The method used in this study was as follows: The obtained 84 datasets were randomly divided, with 70% allocated to the training subset and the remaining 30% reserved for the test subset. This process was repeated 20 times. The results from each iteration were normalized across different groups (spectral band reflectance group, spectral index group, and texture data group). Based on the contribution degree of characteristic variables within each group, the threshold was set to 0.1 to select the pivotal variables for model training and establishing. 2.4.2 Construction of model training The construction of a predictive model necessitates the careful extraction of input features and the judicious selection of appropriate machine learning algorithms. To explore the influence of different features on model performance, the feature combinations were categorized as follows: Single variable groups: (1) Band reflectance. Composed of the original reflectance values of each spectral band. These values were directly reflected the physical and chemical properties of the soil, serving as the primary features for predicting SSC. (2) Spectral Index. Consists of various spectral indices (such as NDVI, RVI, etc.) calculated from different spectral bands. These indices were designed to enhance specific soil or vegetation characteristics. (3) Texture data. Extracting from spectral data, such as the gray level cooccurrence matrix. These features captured subtle changes and spatial structure of the soil surface. Pairwise variable groups: (1) Spectral band reflectance and indices combination: This combination retained the original spectral information while incorporating enhanced information, potentially providing more comprehensive responses of soil salinity changes. (2) Spectral band reflectance and texture feature combination: By integrating the complementary advantages of both types of data, this combination revealed spectral characteristics and captured spatial heterogeneity of the soil surface. (3) Spectral indices and texture feature combination: This combination enriched the enhanced information with spatial structure data, potentially improving the model's generalization. Full Variable Group: This group contained all types data of spectral bands reflectance, spectral indices, and texture features. The comprehensive inclusion of all potentially relevant information was intended to enhance the model's generalization capabilities and predictive performance. 2.5 Machine learning models Based on the SSC and the selected feature variables, four machine learning methods of RF, SVM, GPR, and BPNN methods was used to construct soil salt estimation models. RF was an ensemble learning method that used to establish prediction models for classification and regression problems by using randomly unrelated decision trees 46 . In recent years, RF has been widely used in vegetation growth parameters estimate and the estimation of soil physical and chemical parameters. For instance, Huang et al. 47 established several soil salinity estimation models based on Landsat-8 OLI images in a study of oasis soil salinity in arid areas, and found that the estimation accuracy of the RF modeling method was higher than that of classical statistical models. Similarly, Sui et al. 48 developed a soil salinity estimation model based on original observations and satellite data in a study of coastal soil salinity, utilizing hydrological connectivity measurements and the RF algorithm. SVM was a method that implemented the concept of structural risk minimization, effectively addressing small sample sizes, nonlinearity, and high-dimensional data. SVM offered strong expression capability, generalization ability, and learning efficiency. It easily integrated with multi-source information, thereby achieving higher estimation accuracy 49 . For example, Cai et al. (2010) combined multispectral and texture features and used the SVM classifier to identify soil affected by salt, confirming that the SVM classifier effectively extracted soil salinization distribution information in the Yinchuan Plain. Similarly, Guan et al. 50 introduced SVM theory into the dynamic prediction of soil EC value, constructing a dynamic prediction model for soil salinity aimed at managing irrigation water in salinized areas. GPR was a nonparametric Bayesian regression method that predicted by assuming the data distributed as a multivariate Gaussian distribution, providing an estimate of prediction uncertainty. It was flexible and efficient for small datasets but can be computationally expensive and difficult to scale for larger ones. Additionally, selecting and adjusting the appropriate kernel function required considerable experience and testing. BPNN was a feedforward network composed of multiple neurons capable of learning and identifying nonlinear relationships in complex systems. It exhibited strong self-learning ability, adaptability, and resistance to interference, making it highly promising for the estimation of soil physical and chemical parameters. For instance, Wang et al. 51 successfully established a prediction model for soil moisture and salinity using BPNN with Landsat-8 satellite data. 2.6 Technical workflow This study aims to identify the most effective variables set for accurately predicting SSC, thereby improving the overall performance and applicability of the model in diverse soil depths. SSC was designated as the dependent variable, while different variable groups serving as the independent variables. The sample data were randomly divided into two groups, with 70% allocated for model training and 30% for validation. Four different machine learning methods (RF, SVM, GPR, and BPNN) were employed to estimate SSC. Each method constructed corresponding prediction models, optimizing the performance by adjusting the model hyperparameters. To further evaluate model accuracy, 10-fold cross-validation was used to construct and test the SSC estimation models (Fig. 1 ). 3. Results 3.1 Statistical of soil salinity distribution The salt content of sampling points at different soil depths was categorized as follows: non-saline soil ( 1.0%). The non-salt treatment had an average salt content of 0.214%, the low-salt treatment had an average salt content of 0.257%, the medium-salt treatment had an average salt content of 0.418%, and the high-salt treatment had an average value of 1.353%. The statistical of obtained salt content data were as following: in the 0 to 10 cm soil depth, the measured SSC were relatively lower than those of other depths, which varied from 0.031% to1.04%, with an average of 0.428%. At a depth of 10 to 20 cm, the SSC ranged from 0.087–1.406%, with an average of 0.561%. The SSC in the 20 to 30 cm depth were the highest, with the average SSC of 0.633%, varying between 0.055 ~ 1.806%. The SSC distribution in 30 ~ 40 cm was similar to that in the 10 ~ 20 cm soil layer, with an average of 0.619%. According to these criteria, the measured salt grade distribution in the study area was shown in Fig. 2 . 3.2 Contribution and selection of feature variables In this study, the important factor of 77 feature variables of three types data for different soil depths were calculated (Fig. 3 ). After normalizing the feature variables important factor of different groups and setting a threshold of 0.1, the feature variables were screened. Common variables across different depths were selected for the ML model training. Finally, the selected variables for different groups were as follow: Band reflectance: B, R-edge, and NIR; spectral Index: SI2, SI2-reg, BI, SI-T, DVI, EVI, SAVI, OSAVI, DVI-reg, EVI-reg, and MSAVI; Texture Data B-MEA, B-VAR, B-CON, B-ENT, B-SEC, G-VAR, RE-MEA, NIR-MEA, NIR-HOM, NIR-ENT, and NIR-SEC. Altogether, 25 feature variables were filtered for modeling. 3.3 Comparision of model performance for various ML methods and data groups The optimal variables of the three variables types of band reflectance, spectral index, and texture data and their combinations were used as the independent variables and the SSC as the target variable imputing into the machine learning (RF, SVM, GPR, and BPNN) model to establish the soil salt prediction model. The prediction accuracies of these models based on different variable combinations were shown in Figs. 4 and 5 . Figures 4 and 5 demonstrated that for each algorithm, the band reflectance and spectral index variables performed well in estimating SSC at a depth of 0 to 10 cm within the single variable groups. Among these, the BPNN model, based on the spectral index, achieved the best performance with an R² of 0.74 and RMSE of 0.26%. The whole variable combination-based GPR model was the best in the multivariate combination group. It exhibited a stableR² value around 0.77, an RMSE of 0.24%, and great overall stability. The RF model identified the band reflectance as the most single variable for SSC at a depth of 10 to 20 cm, whereas the SVM model produced the worst estimation. Across the four approaches tested, texture data yielded suboptimal results, with an average RMSE exceeding 0.34%. However, The RF model that incorporated both texture and spectral index data outperformed others in the multivariate combination, achieving a mean RMSE of 0.31% and an R² of approximately 0.61, making it the most effective model overall. Additionally, the GPR model's results outperformed those of other approaches. In contrast, the SVM and BPNN models consistently delivered less accurate predictions across all groups. At a depth of 20 to 30 cm, the RF model built using band reflectance and spectral index demonstrated superior performance, with average R² values of 0.62 and 0.67 and mean RMSE of 0.31% and 0.29%, respectively. This model outperformed than the models built by the other three algorithm. Within the multivariate combination, the spectral index and texture data performed better outcomes in the RF model compared to other combinations, yielding an R² of 0.67 and an RMSE of 0.29%. Similarly, the GPR algorithm, when applied to the texture data showed higher accuracy than the other three algorithms. The SVM model achieved its highest accuracy when all variables were combined, while the accuracy of models based on band reflectance and spectral index data in the GPR and BPNN algorithms was marginally higher than that of other combinations, with R² values of 0.64 and 0.58, and RMSE values of 0.30% and 0.35%, respectively. The RF model built using textural data as a single variable had the lowest accuracy at a depth of 30 to 40 cm, with an average RMSE of 0.39% and R² of 0.37. In contrast, the RF model constructed using the spectral index showed higher accuracy than the other three models. Within the multi-variable combination group, the SVM model retained consistent accuracy, with the R² of approximately 0.55. The GPR model, when applied to the full variable set, demonstrated superior accuracy compared to other combinations. Notably, the BPNN model exhibited better accuracy when based on band reflectance and spectral index data than other models. In conclusion, soil salinity estimation at different depths were improved by applying a variety of machine learning techniques, with the effectiveness of these methods differing across modeling groups. The most effective method for estimating soil salinity was the whole variable group based GPR model, particularly for both surface (0ཞ10 cm) and deep (30ཞ40 cm) measurements. However, for middle depth soils (10–20 cm and 20–30 cm), the RF model based on spectral index and texture data displayed the best outcomes. The RF model effectively reduced the variance of a single tree by integrating multiple decision trees, thereby improving the model's stability and estimation accuracy. While the SVM and BPNN models performed marginally worse overall than GPR and RF, they were still fairly good at predicting within particular categories. 3.4 Outer-sample validation and analysis The performances of optimal estimation models have been evaluated 20 times using the outer-samples (not training) data (Fig. 6 ). At a soil depth of 0ཞ10 cm, the R² values of the validation set fluctuated between 0.6 and 0.9, with an average R²of 0.77. The RMSE values typically ranged from 0.1 to 0.3%, with an average RMSE of 0.185%. This indicated that the GPR model had high prediction accuracy and stability in predicting shallow soil salinity. For the 10ཞ20 cm depth, the RF model showed favorable results, with an average R² of 0.67. The RMSE values ranged between 0.2 and 0.4%, with an average RMSE of 0.325%. At the 20 ~ 30 cm depth, the RF model achieved a maximum R 2 of 0.87 and minimum RMSE of 0.212%. At a deepest soil layer, the model's R 2 ranged from a minimum of t 0.52 to a maximum of 0.92, indicating high overall accuracy and stability. The RMSE values fluctuated between 0.169% and 0.485%, with an average RMSE of 0.31%. The modeling and validation accuracies of four soil salinity estimation models at different depths were compared. Models were selected at each depth to achieve better estimation results in both modeling and validation. The average R 2 values in the modeling set at each depth exceeded 0.6, and the average RMSE were below 0.4%, indicating that the model had high accuracy and good stability. 3.5 Performance of salinity estimations under different vegetation coverage In this study, the LAI was used as a proxy for vegetation cover, which is divided into three different levels: low, medium and high. Correspondingly, soil salinity data are stratified based on these cover classes. The results of SSC estimation under various covers were presented in Fig. 7 . The accuracy of the salt estimation models at various depths was assessed based on the salt data collected under different vegetation coverage levels, and corresponding scatter plots were created. As shown in Fig. 5 , all models achieved R² values greater than 0.4, with the highest value reaching 0.83, indicating that the model's accuracy was highest for soil surface with intermediate vegetation covering. The models demonstrated good accuracy for soil depths of 0 to 10 cm under each vegetation coverage level, with R² values exceeding 0.7 and low RMSE values. At a depth of 10 to 20 cm, the model accuracy was highest under low vegetation covering, while the accuracy under medium vegetation coverage was lowest. However, at a depth of 20 to 30 cm, the maximum validation accuracy was achieved with medium coverage, with R² values exceeding 0.7 across all three coverages levels. With considerable vegetation coverage, the maximum accuracy at 30 to 40 cm was achieved, with an RMSE of 0.241% and an R² of 0.78. Generally, when using the data selected under each vegetation coverage level as the validation set, the favored models at different depths showed good estimation performance. 4. Discussion Soil salinity estimation by UAV multispectral remote sensing data plays a positive role in salinity monitoring and management. In this study, the spectral data and salinity data acquired during the vegetation cover period were used to construct and validate the estimation model for soil salinity at different depths. In addition, the model performances were validated under different vegetation cover ratios and the impact of cover on salt estimation in different soil layers was also explored. Feature selection is an important part of constructing the soil salinity estimation model. Through the BDT method, 25 key features were finally identified, including bands reflectance (e.g., blue, red-edge, and near-infrared), spectral indices (e.g., SI2-reg, EVI, and DVI, etc.), and texture data (e.g., the mean, variance, and second-order moments of the GLCM features, etc.). These features efficiently captured the spectral properties of the soil and the spatial structural information, providing the model with sufficient input variables. Numerous studies have demonstrated that the red edge band is more sensitive to soil salinity and is strongly correlated the spectral properties of the vegetation canopy, thereby improving the accuracy of soil salinity estimation 30 , 52 . To facilitate modeling and prediction, Hu et al. 53 employed UAV to gather remote sensing images of areas with both vegetation cover and barren ground. The regions with vegetation cover showed the most accurate prediction results. As a result, the spectral index produced by incorporating the red edge band significantly increases the accuracy of the soil salinity estimation in saline agricultural lands with vegetation cover. The red edge and near-infrared bands showed high importance at all depths in this study, as they sensitively revealed the changes in soil moisture and salt content, serving as important indicators for soil salinity prediction. This result aligned with the reported of Taghadosi et al. 54 . Spectral indices such as EVI, EVI-reg and BI indirectly indicated the salinity status of the soil by enhancing the vegetation characteristics. Lobell et al. 55 showed that EVI is more reliable than NDVI for salinity monitoring. Spectral index groups, such as SI2-reg, EVI-reg, and DVI-reg, involved the calculation of the red-edge band, further proving the importance of the red-edge band in soil salinity monitoring, consistent with the results of Ma et al. 56 Textural features, such as the mean, variance, and contrast in the GLCM, were instrumental in capturing subtle soil surface changes and spatial structure information. Tai 57 showed that incorporating texture features could effectively improve the accuracy of soil salinity estimation models under vegetation cover conditions. Different modelling groups also exerted an influence on estimation model performance. In the single-variable modelling group, spectral indices had a greater contribution to soil salinity estimation, with the RF model demonstrating high accuracy across all four depths (Fig. 4 ). Bian et al. 58 found that there is a certain correlation between spectral index and soil salinity, which has a great contribution to the application of estimation. Although the accuracy of the model built by band reflectance was slightly lower than that of spectral indices, the overall difference was not significant, which was in agreement with Wang et al. 59 , who evaluated the prediction accuracy of different variable groups for soil salinity in oasis. In contrast, texture data performed poorly among the single variables. However, when combined with other variables, the texture information from multispectral images could improve the accuracy of soil salinity estimation at different depths of vegetation cover. Zheng et al. 18 showed that the combination of texture data and spectral information significantly improved the accuracy of rice biomass estimation. The texture features of the UAV multispectral images provided rich information 60 , making them applicable for estimating soil salinity at different depths. Incorporating texture data significantly enhances model performance, particularly in the RF model. When texture data is combined with other data, the accuracy of the model is notably improved. This study also demonstrates that using sensitive bands and spectral indices as the input variables leads to better modeling and validation outcomes in soil salinity estimation models. Nevertheless, in multi-variate groups, variables may affect and constrain each other, and combining variables of high importance not always achieve optimal results. For instance, the accuracy of all groups in the BPNN model was not optimal compared to other groups. Introducing an excessive number of independent variables could lead to information redundancy, overfitting, and a decrease in model accuracy 61 . Previous studies have highlighted the superiority of machine learning methods in estimating soil salinity content due to the complexity and indirect relationships between variables 62 . In this study, RF, SVM, GPR and BPNN were used to model soil salinity at different depths and to select the most effective model. The evaluation of model criteria revealed that RF and GPR have advantages in estimating soil salinity. Although SVM has strong generalization capability for handling nonlinear problems, it is sensitive to noisy data 63 . The models exhibited varying levels of accuracy in soil salinity estimation at different depths. The GPR model demonstrates better prediction accuracy than other methods for surface and deep soils, while the RF model performed well in intermediate soils. Due to the complexity of soil salinization mechanism and the intricate nonlinear relationship between soil spectral characteristics and soil salinity, machine learning is well-suited for revealing these relationships. Their strong nonlinear fitting and generalization capabilities, make them ideal for simulating the complex interactions among variables. Ma et al. 56 found that the RF model had better estimation results when predicting the soil salinity in Ebinur Lake wetland using multispectral and Digital Elevation Model (DEM) data. Therefore, machine learning algorithms have shown superior performance in both modeling -and validation of soil salinity estimation model. For instance, Wei et al. 5 developed SSC estimation models using RF, SVM, and BPNN algorithms, with the RF model achieving the highest prediction accuracy (R 2 = 0.84). Additionally, Zhu et al. 64 and Yu et al. 65 concluded that the RF based SSC estimation model had a high accuracy. The essence of estimating soil salinity during the vegetation cover period is to indirectly obtain information on soil salinity through the crop spectral response to soil salinity. The vegetation cover ratio of the given crop in the same period may correlate with the degree of soil salinity. The optimal estimation model for various soil depth varied with vegetation cover conditions, with the optimal model on depths ranging from 0ཞ10 cm under both low and medium vegetation cover conditions, and from 20ཞ30 cm under high vegetation cover. During the vegetation cover period, the absorption of soil water by crops was mainly done by the lateral root. Consequently, the soil salinity in the soil layer where the lateral roots were located will significantly affect crop water absorption 66 , causing crop growth stress. Higher soil salinity levels result in more severe stress and poorer vegetation growth, which would be indirectly expressed in the vegetation canopy. 5. Conclusion In this study, SSC at different depths was estimated and analyzed based on UAV multi-spectral data by composing seven type groups from band reflectance, spectral index and texture feature. Moreover, feature variables selection based on BDT method and four machine learning algorithms (RF, SVM, GPR and BPNN) were incorporated in estimation model establishment. (1) The estimation accuracy of the RF and GPR models surpasses that of the SVM and BPNN models across all four soil depths in the SSC estimation models. (2) The GPR model, based on the entire variable group, is the most accurate for estimating soil salt content at depths of 0 ~ 10 cm and 30 ~ 40 cm. The R 2 values for the validation set were 0.774 and 0.705, and the modeling set, 0.77 and 0.62, with corresponding RMSE of 0.185% and 0.31%, respectively. The RF model, based on the spectral index and texture data set, is the most accurate estimation model for depths of 10 to 20 cm and 20 to 30 cm. (3) The model was further validated using salt data under varying vegetation covers. The highest validation accuracy was obtained with medium vegetation coverage at depths of 0 ~ 10 cm and 20 ~ 30 cm, as well as high vegetation coverage at depths of 10 ~ 20 cm and 30 ~ 40 cm. (4) The model imputes were screened using the BDT approach to select the independent variables at each depth. The selected results demonstrated that the model constructed using only texturing data group had low accuracy, but the accuracy increased dramatically when additional variables were included. There are still some problems and limitations in this study. This experiment was conducted under box planting conditions, and the experimental results have some limitations, which may not fully reflect the soil salinity conditions under actual field conditions. In addition, the study has not yet considered the effects of crop type and soil water content on the estimation model, which may have a significant impact on model accuracy in practical applications. Therefore, future studies need to consider more influencing factors in order to further optimize the estimation model and improve the accuracy and stability of soil salinity prediction. Declarations Competing Interests Statement The authors declare no competing interests. Funding This research was funded by the National Natural Science Foundation of China (52379049, 52209071), the Postgraduate Research & Practice Innovation Program of Jiangsu Province (Yangzhou University) (No. SJCX23_1945), the China Postdoctoral Science Foundation (2023T160552, 2020M671623), “Chunhui Plan” Cooperative Scientific Research Project of Ministry of Education of China (HZKY20220115), the “Blue Project” of Yangzhou University and the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD). Author Contribution Conceptualization, C.Z.; methodology, Z.L.; software, Z.L.; validation, M.D. and M.T.; formal analysis, Z.L.; resources, M.D.; data curation, M.D. and R.L.; Writing—original preparation, Z.L.; Writing—review and editing, C.Z. and M.T.; visualization, Z.Z.; supervision, C.Z. and S.F.; project administration, C.Z.; funding acquisition, C.Z. and S.F. All authors have read and agreed to the published version of the manuscript. Acknowledgement The authors would like to acknowledge the Agricultural Water and Hydrological Ecology Experiment Station, Yangzhou University, for providing the experimental facilities. Data Availability The datasets used and analyzed during the current study are available from the corresponding author upon reasonable request. References Seifi, M., Ahmadi, A., Neyshabouri, M. R., Taghizadeh-Mehrjardi, R. & Bahrami, H. A. Remote and Vis-NIR spectra sensing potential for soil salinization estimation in the eastern coast of Urmia hyper saline lake, Iran 20100398 (Society and Environment, 2020). Wang, L. W. & Wei, Y. X. Estimating the total nitrogen and total phosphorus content of wetland soils using hyperspectral models. Acta Ecol. Sin . 36 (16), 5116–5125 (2016). Rao, B. R. M. et al. Spectral behaviour of salt-affected soils. Int. J. Remote Sens. 16 (12), 2125–2136 (1995). Zhao, W., Zhou, C., Zhou, C., Ma, H. & Wang, Z. Soil salinity inversion model of oasis in arid area based on UAV multispectral remote sensing. Remote Sens. 14 (8), 1804 (2022). Wei, G. et al. Estimation of soil salt content by combining UAV-borne multispectral sensor and machine learning algorithms. PeerJ , 8 , e9087. (2020). Chen, J. et al. UAV remote sensing inversion of soil salinity in field of sunflower. Trans. Chin. Soc. Agricultural Mach. 51 (7), 178–191 (2020). Ivushkin, K. et al. UAV based soil salinity assessment of cropland. Geoderma . 338 , 502–512 (2019). Nicolas, H. & Walter, C. Detecting salinity hazards within a semiarid context by means of combining soil and remote-sensing data. Geoderma . 134 (1–2), 217–230 (2006). Qi, G., Zhao, G. & Xi, X. Soil salinity inversion of winter wheat areas based on satellite-unmanned aerial vehicle-ground collaborative system in coastal of the Yellow River Delta. Sensors . 20 (22), 6521 (2020). Wang, N. et al. Integrating remote sensing and landscape characteristics to estimate soil salinity using machine learning methods: A case study from Southern Xinjiang, China. Remote Sens. 12 (24), 4118 (2020). Farifteh, J., Van der Meer, F., Atzberger, C. & Carranza, E. J. M. Quantitative analysis of salt-affected soil reflectance spectra: A comparison of two adaptive methods (PLSR and ANN). Remote Sens. Environ. 110 (1), 59–78 (2007). Xu, L. & Wang, Q. Retrieval of soil water content in saline soils from emitted thermal infrared spectra using partial linear squares regression. Remote Sens. 7 (11), 14646–14662 (2015). Ge, X. et al. Combining UAV-based hyperspectral imagery and machine learning algorithms for soil moisture content monitoring. PeerJ , 7 , e6926. (2019). Poblete, T., Ortega-Farías, S., Moreno, M. A. & Bardeen, M. Artificial neural network to predict vine water status spatial variability using multispectral information obtained from an unmanned aerial vehicle (UAV). Sensors . 17 (11), 2488 (2017). Hu, J. Estimation of Soil Salinity in Arid Area Based on Multi-Source Remote Sensing (Doctoral dissertation, Zhejiang University). (2019). Zhang Fei, Z. F. et al. I., Extracting of soil salinization by SVM and accuracy evaluation based on texture characteristic. (2009). Zhang, Z. et al. Inversion of soil salt content by UAV remote sensing under different vegetation coverage. Trans. Chin. Soc. Agric. Mach. 53 , 220–230 (2022). Zheng, H. et al. Enhancing the nitrogen signals of rice canopies across critical growth stages through the integration of textural and spectral information from unmanned aerial vehicle (UAV) multispectral imagery. Remote Sens. 12 (6), 957 (2020). Yang, H., Wang, Z., Cao, J., Wu, Q. & Zhang, B. Estimating soil salinity using Gaofen-2 imagery: A novel application of combined spectral and textural features. Environ. Res. 217 , 114870 (2023). Ren, J., Li, X., Zhao, K., Fu, B. & Jiang, T. Study of an on-line measurement method for the salt parameters of soda-saline soils based on the texture features of cracks. Geoderma . 263 , 60–69 (2016). Andrade, R. et al. Proximal sensor data fusion and auxiliary information for tropical soil property prediction: Soil texture. Geoderma . 422 , 115936 (2022). Drielsma, M. J., Love, J., Taylor, S., Thapa, R. & Williams, K. J. General Landscape Connectivity Model (GLCM): a new way to map whole of landscape biodiversity functional connectivity for operational planning and reporting. Ecol. Model. 465 , 109858 (2022). Cai, S., Zhang, R., Liu, L. & Zhou, D. A method of salt-affected soil information extraction based on a support vector machine with texture features. Math. Comput. Model. 51 (11–12), 1319–1325 (2010). Huang, C. et al. Land cover mapping in cloud-prone tropical areas using Sentinel-2 data: Integrating spectral features with Ndvi temporal dynamics. Remote Sens. 12 (7), 1163 (2020). Jia, K. L. & Zhang, J. H. Impacts of different alkaline soil on canopy spectral characteristics of overlying vegetation. Spectrosc. Spectr. Anal. 34 (3), 782–786 (2014). Zhang, J. et al. Estimating soil salinity with different fractional vegetation cover using remote sensing. Land. Degrad. Dev. 32 (2), 597–612 (2021). Wang, F. et al. Sensitivity analysis of soil salinity and vegetation indices to detect soil salinity variation by using Landsat series images: applications in different oases in Xinjiang, China. Acta Ecol. Sin. 37 (15), 5007–5022 (2017). Chi, Y., Sun, J., Liu, W., Wang, J. & Zhao, M. Mapping coastal wetland soil salinity in different seasons using an improved comprehensive land surface factor system. Ecol. Ind. 107 , 105517 (2019). Khan, N. M., Rastoskuev, V. V., Shalina, E. V. & Sato, Y. Mapping salt-affected soils using remote sensing indicators-a simple approach with the use of GIS IDRISI. In 22nd Asian conference on remote sensing (Vol. 5, No. 9). (2001), November. Zhang, S. et al. Integrated satellite, unmanned aerial vehicle (UAV) and ground inversion of the SPAD of winter wheat in the reviving stage. Sensors . 19 (7), 1485 (2019). Song, C. Y., Ren, H. X. & Huang, C. Estimating soil salinity in the Yellow River Delta, Eastern China—an integrated approach using spectral and terrain indices with the generalized additive model. Pedosphere . 26 (5), 626–635 (2016). Abbas, A., Khan, S., Hussain, N., Hanjra, M. A. & Akbar, S. Characterizing soil salinity in irrigated agriculture using a remote sensing approach 5543–52 (Physics and chemistry of the Earth, 2013). Tripathi, N. K., Rai, B. K. & Dwivedi, P. Spatial modeling of soil alkalinity in GIS environment using IRS data. In Proceedings of the 18th Asian Conference on Remote Sensing, Kuala Lumpur, Malaysia (Vol. 20, No. 25, pp. 81–86). (1997), October. Allbed, A., Kumar, L. & Aldakheel, Y. Y. Assessing soil salinity using soil salinity and vegetation indices derived from IKONOS high-spatial resolution imageries: Applications in a date palm dominated region. Geoderma . 230 , 1–8 (2014). Zhou, L. et al. Research on quick dynamic monitoring of soil salinization in Qaidam Basin[J]. Science of Surveying and Mapping, 46(7): 99 ~ 106, 114. (2021). Khan, N. M., Rastoskuev, V. V., Sato, Y. & Shiozawa, S. Assessment of hydrosaline land degradation by using a simple approach of remote sensing indicators. Agric. Water Manage. 77 (1–3), 96–109 (2005). Liu, H. Q. & Huete, A. A feedback based modification of the NDVI to minimize canopy background and atmospheric noise. IEEE Trans. Geosci. Remote Sens. 33 (2), 457–465 (1995). Lyon, J. G., Yuan, D., Lunetta, R. S. & Elvidge, C. D. A change detection experiment using vegetation indices. Photogram. Eng. Remote Sens. 64 (2), 143–150 (1998). Birth, G. S. & McVey, G. R. Measuring the color of growing turf with a reflectance spectrophotometer 1. Agron. J. 60 (6), 640–643 (1968). Arachchi, M. H., Field, D. J. & McBratney, A. B. Quantification of soil carbon from bulk soil samples to predict the aggregate-carbon fractions within using near-and mid-infrared spectroscopic techniques. Geoderma . 267 , 207–214 (2016). Ma, S., He, B., Ge, X. & Luo, X. Spatial prediction of soil salinity based on the Google Earth Engine platform with multitemporal synthetic remote sensing images. Ecol. Inf. 75 , 102111 (2023). Gitelson, A. A., Kaufman, Y. J., Stark, R. & Rundquist, D. Novel algorithms for remote estimation of vegetation fraction. Remote Sens. Environ. 80 (1), 76–87 (2002). Elsherbiny, O., Zhou, L., Feng, L. & Qiu, Z. Integration of visible and thermal imagery with an artificial neural network approach for robust forecasting of canopy water content in rice. Remote Sens. 13 (9), 1785 (2021). Xiao, Y., Dong, Y., Huang, W., Liu, L. & Ma, H. Wheat fusarium head blight detection using UAV-based spectral and texture features in optimal window size. Remote Sens. 13 (13), 2437 (2021). Drucker, H. & Cortes, C. Boosting decision trees. Advances in neural information processing systems, 8. (1995). Mutanga, O., Adam, E. & Cho, M. A. High density biomass estimation for wetland vegetation using WorldView-2 imagery and random forest regression algorithm. Int. J. Appl. Earth Obs. Geoinf. 18 , 399–406 (2012). Huang, X. Y., Wang, X. M. & KAWUQIATI, B. Inversion of soil salinity of an oasis in an arid area based on Landsat8 OLI images. Remote Sens. Nat. Resour. , 35 (1). (2023). Sui, H. et al. Soil salinity estimation over coastal wetlands based on random Forest algorithm and hydrological connectivity metric. Front. Mar. Sci. 9 , 895172 (2022). Tan, K., Ye, Y. Y., Du, P. J. & Zhang, Q. Q. Estimation of heavy metal concentrations in reclaimed mining soils using reflectance spectroscopy. Spectrosc. Spectr. Anal. 34 (12), 3317–3322 (2014). Guan, X., Wang, S., Gao, Z. & Lv, Y. Dynamic prediction of soil salinization in an irrigation district based on the support vector machine. Math. Comput. Model. 58 (3–4), 719–724 (2013). Wang, J., Wang, W., Hu, Y., Tian, S. & Liu, D. Soil moisture and salinity inversion based on new remote sensing index and neural network at a salina-alkaline wetland. Water . 13 (19), 2762 (2021). Kaplan, G. & Avdan, U. Evaluating the utilization of the red edge and radar bands from sentinel sensors for wetland classification. Catena . 178 , 109–119 (2019). Hu, J. et al. Quantitative estimation of soil salinity using UAV-borne hyperspectral and satellite multispectral images. Remote Sens. 11 (7), 736 (2019). Taghadosi, M. M., Hasanlou, M. & Eftekhari, K. Retrieval of soil salinity from Sentinel-2 multispectral imagery. Eur. J. Remote Sens. 52 (1), 138–154 (2019). Lobell, D. B. et al. Regional-scale assessment of soil salinity in the Red River Valley using multi‐year MODIS EVI and NDVI. J. Environ. Qual. 39 (1), 35–41 (2010). Ma, G., Ding, J., Han, L. & Zhang, Z. Digital mapping of soil salinization in arid area wetland based on variable optimized selection and machine learning. Trans. Chin. Soc. Agric. Eng. 36 , 124–131 (2020). Tai, X. Soil Salinity Monitoring Model on UAV Multi Spectral Remote Sensing Vegetation Cover[D]. Northwest Agricultural and Forestry University, DOI: (2022). 10.27409/d Bian, L., Wang, J., Guo, B., Cheng, K. & Wei, H. Remote sensing extraction of soil salinity in Yellow River Delta Kenli County based on feature space. Remote Sens. Technol. Application . 35 (1), 211–218 (2020). Wang, F. et al. Environmental sensitive variable optimization and machine learning algorithm using in soil salt prediction at oasis. Trans. Chin. Soc. Agric. Eng. 34 , 102–110 (2018). Guo, Y. et al. Integrating spectral and textural information for identifying the tasseling date of summer maize using UAV based RGB images. Int. J. Appl. Earth Obs. Geoinf. 102 , 102435 (2021). Morellos, A. et al. Machine learning based prediction of soil total nitrogen, organic carbon and moisture content by using VIS-NIR spectroscopy. Biosyst. Eng. 152 , 104–116 (2016). Cui, X. et al. Estimating soil salinity under sunflower cover in the Hetao Irrigation District based on unmanned aerial vehicle remote sensing. Land. Degrad. Dev. 34 (1), 84–97 (2023). Shieh, H. L. & Kuo, C. C. A reduced data set method for support vector regression. Expert Syst. Appl. 37 (12), 7781–7787 (2010). Zhu, C., Ding, J., Zhang, Z. & Wang, Z. Exploring the potential of UAV hyperspectral image for estimating soil salinity: Effects of optimal band combination algorithm and random forest. Spectrochim. Acta Part A Mol. Biomol. Spectrosc. 279 , 121416 (2022). Yu, X., Chang, C., Song, J., Zhuge, Y. & Wang, A. Precise monitoring of soil salinity in China’s Yellow River Delta using UAV-borne multispectral imagery and a soil salinity retrieval index. Sensors . 22 (2), 546 (2022). Zhang, Z. T. et al. UAV Multispectral Remote Sensing Soil Salinity Inversion Based on Different Fractional Vegetation Coverages[J]. Trans. Chin. Soc. Agricultural Mach. , 53 (08): 220–230. (2022). Additional Declarations No competing interests reported. Cite Share Download PDF Status: Published Journal Publication published 21 Jan, 2025 Read the published version in Scientific Reports → Version 1 posted Editorial decision: Revision requested 30 Oct, 2024 Reviewers agreed at journal 19 Oct, 2024 Reviews received at journal 19 Oct, 2024 Reviews received at journal 09 Oct, 2024 Reviewers agreed at journal 02 Oct, 2024 Reviewers agreed at journal 28 Sep, 2024 Reviewers invited by journal 26 Sep, 2024 Editor assigned by journal 26 Sep, 2024 Editor invited by journal 05 Sep, 2024 Submission checks completed at journal 03 Sep, 2024 First submitted to journal 25 Aug, 2024 You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-4971758","acceptedTermsAndConditions":true,"allowDirectSubmit":false,"archivedVersions":[],"articleType":"Article","associatedPublications":[],"authors":[{"id":361580570,"identity":"cfc92871-e3d9-415c-a33d-95b1ec8fa328","order_by":0,"name":"Zhenhai Luo","email":"","orcid":"","institution":"Yangzhou University","correspondingAuthor":false,"prefix":"","firstName":"Zhenhai","middleName":"","lastName":"Luo","suffix":""},{"id":361580571,"identity":"8fae637d-bf2d-45c0-a0ea-f7f5e240aaab","order_by":1,"name":"Meihua Deng","email":"","orcid":"","institution":"Yangzhou University","correspondingAuthor":false,"prefix":"","firstName":"Meihua","middleName":"","lastName":"Deng","suffix":""},{"id":361580572,"identity":"693defc4-1d39-4a43-b861-d37c44d4fd3a","order_by":2,"name":"Min Tang","email":"","orcid":"","institution":"Yangzhou University","correspondingAuthor":false,"prefix":"","firstName":"Min","middleName":"","lastName":"Tang","suffix":""},{"id":361580573,"identity":"899529a8-1338-4cdf-8f29-e4d62f6c8f42","order_by":3,"name":"Rui Liu","email":"","orcid":"","institution":"Yangzhou University","correspondingAuthor":false,"prefix":"","firstName":"Rui","middleName":"","lastName":"Liu","suffix":""},{"id":361580574,"identity":"c6e2412c-fc7b-4f65-af57-32993963334d","order_by":4,"name":"Shaoyuan Feng","email":"","orcid":"","institution":"Yangzhou University","correspondingAuthor":false,"prefix":"","firstName":"Shaoyuan","middleName":"","lastName":"Feng","suffix":""},{"id":361580575,"identity":"ababf508-dfa5-41df-995d-093bfaea73c4","order_by":5,"name":"Chao Zhang","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAA9ElEQVRIiWNgGAWjYDACCQhZ38bMfAAicoA4LRaM/extCSRpqWCc2XPGgDgt8rObnz38UiPBbHAj55vkzzYGOb4bCYyfC/BoYZxzzNxY5pgEm8GN3G0Skm0MxpI3EpilZ+DRwiyRYCYt2SDBA9Zi2MaQuOFGAhszDx4tbBLp30BaJIAOeyaR2MZQT1ALj0SOmeTHBgkDyZ4zbBIH2xgSDAhpkZDIKZNmOCaRwM/eZmzZcE7CcOaZh83S+LTIz0jfJvmjpg5oMvPDmz/KbOT5jicf/IxPCwjAnMEiAYkmxgYCGoBKfkC1fiCodBSMglEwCkYkAAD09UfNz0w9CwAAAABJRU5ErkJggg==","orcid":"","institution":"Yangzhou University","correspondingAuthor":true,"prefix":"","firstName":"Chao","middleName":"","lastName":"Zhang","suffix":""},{"id":361580576,"identity":"aac3e805-eead-4179-90b0-22f2b1ba900a","order_by":6,"name":"Zhen Zheng","email":"","orcid":"","institution":"Jiangsu University","correspondingAuthor":false,"prefix":"","firstName":"Zhen","middleName":"","lastName":"Zheng","suffix":""}],"badges":[],"createdAt":"2024-08-25 08:04:10","currentVersionCode":1,"declarations":"","doi":"10.21203/rs.3.rs-4971758/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-4971758/v1","draftVersion":[],"editorialEvents":[{"content":"https://doi.org/10.1038/s41598-024-82868-9","type":"published","date":"2025-01-21T15:57:29+00:00"}],"editorialNote":"","failedWorkflow":false,"files":[{"id":66953778,"identity":"d1843fa6-c3ca-41e5-ad88-e6f9139e8641","added_by":"auto","created_at":"2024-10-18 10:48:54","extension":"png","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":175647,"visible":true,"origin":"","legend":"\u003cp\u003eTechnical workflow for soil salinity estimation model under vegetation cover based on multi-source remote sensing\u003c/p\u003e","description":"","filename":"1.png","url":"https://assets-eu.researchsquare.com/files/rs-4971758/v1/bb13f255c53d6b4e52fee06e.png"},{"id":66953347,"identity":"4285d8cc-0782-4d7b-9765-742077b4be3c","added_by":"auto","created_at":"2024-10-18 10:40:54","extension":"png","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":42689,"visible":true,"origin":"","legend":"\u003cp\u003eStatistics of soil salinity in different soil depths (a) experimental treatments (b).\u003c/p\u003e","description":"","filename":"2.png","url":"https://assets-eu.researchsquare.com/files/rs-4971758/v1/11d273094c3cb89e96df4589.png"},{"id":66954526,"identity":"333e26e5-e1c5-4a52-b674-0831e7ed3183","added_by":"auto","created_at":"2024-10-18 10:56:54","extension":"png","order_by":3,"title":"Figure 3","display":"","copyAsset":false,"role":"figure","size":111841,"visible":true,"origin":"","legend":"\u003cp\u003eContributions of feature variables of three data types for SSC variations under different depths.\u003c/p\u003e","description":"","filename":"3.png","url":"https://assets-eu.researchsquare.com/files/rs-4971758/v1/2197333dfa4e3b17de20b8c6.png"},{"id":66953344,"identity":"5d7788b1-8a88-4743-974a-df3189d61165","added_by":"auto","created_at":"2024-10-18 10:40:54","extension":"png","order_by":4,"title":"Figure 4","display":"","copyAsset":false,"role":"figure","size":92770,"visible":true,"origin":"","legend":"\u003cp\u003ePerformance of different ML methods and data groups for SSC estimation. (a ~ d) Summary of 0~10cm SSC R² estimation results; (e ~ h) Summary of 10~20cm SSC R² estimation results; (i ~ l) Summary of 20~30cm SSC R² estimation results; (m ~ p) Summary of 30~40cm SSC R² estimation results.\u003c/p\u003e","description":"","filename":"4.png","url":"https://assets-eu.researchsquare.com/files/rs-4971758/v1/010a36150e00147d2947092d.png"},{"id":66953350,"identity":"685d4a89-74da-446e-93ae-fbbd2645dcc8","added_by":"auto","created_at":"2024-10-18 10:40:54","extension":"png","order_by":5,"title":"Figure 5","display":"","copyAsset":false,"role":"figure","size":97063,"visible":true,"origin":"","legend":"\u003cp\u003ePerformance of different ML methods and data groups for SSC estimation. (a ~ d) Summary of 0~10cm SSC RMSE estimation results; (e ~ h) Summary of 10~20cm SSC RMSE estimation results; (i ~ l) Summary of 20~30cm SSC RMSE estimation results; (m ~ p) Summary of 30~40cm SSC RMSE estimation results.\u003c/p\u003e","description":"","filename":"5.png","url":"https://assets-eu.researchsquare.com/files/rs-4971758/v1/00e428b0931b332decb89133.png"},{"id":66953349,"identity":"b0e36c3f-8507-48a2-bfcc-d8297d001017","added_by":"auto","created_at":"2024-10-18 10:40:54","extension":"png","order_by":6,"title":"Figure 6","display":"","copyAsset":false,"role":"figure","size":54692,"visible":true,"origin":"","legend":"\u003cp\u003eAccuracy of optimal SSC estimation model for different soil layers\u003c/p\u003e","description":"","filename":"6.png","url":"https://assets-eu.researchsquare.com/files/rs-4971758/v1/83f4a39dd869b45c00b36522.png"},{"id":66953348,"identity":"d2fa0683-ef18-4649-9d96-05a21bc483fa","added_by":"auto","created_at":"2024-10-18 10:40:54","extension":"png","order_by":7,"title":"Figure 7","display":"","copyAsset":false,"role":"figure","size":122679,"visible":true,"origin":"","legend":"\u003cp\u003e(a ~c) Scatter plots of measured and predicted data at different coverage levels of 0~10 cm; (d ~f) Scatter plots of measured and predicted data at different coverage levels of 10~20 cm; (g ~i) Scatter plots of measured and predicted data at different coverage levels of 20~30 cm; (j ~l) Scatter plots of measured and predicted data at different coverage levels of 20~30 cm.\u003c/p\u003e","description":"","filename":"7.png","url":"https://assets-eu.researchsquare.com/files/rs-4971758/v1/b1332b1d6b7f7400828df5a9.png"},{"id":74858381,"identity":"ecc41da5-c8d4-42a3-8047-293e0826ca07","added_by":"auto","created_at":"2025-01-27 16:08:48","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":1565566,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-4971758/v1/02d99017-af93-46bf-bf1c-05dfe0a240ff.pdf"}],"financialInterests":"No competing interests reported.","formattedTitle":"Estimating soil profile salinity under vegetation cover based on UAV multi-source remote sensing","fulltext":[{"header":"1. Introduction","content":"\u003cp\u003eIn China, the soil salinization area extends 9.21\u0026nbsp;million hectares, accounting for 6.62% of the total cultivated land. The process of saline soil formation is intricate, posing significant challenging for the detection and dynamic monitoring of soil salinization. Conventionally, soil salinization was assessed through field sampling and chemical analysis, these methods that were both time consuming, and labor-intensive\u003csup\u003e\u003cspan citationid=\"CR1\" class=\"CitationRef\"\u003e1\u003c/span\u003e\u003c/sup\u003e. In contrast, remote sensing technology, a rapid and expansive means to gather data on ground objects across different temporal and spatial, is an ideal tool for soil salinity monitoring\u003csup\u003e\u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2\u003c/span\u003e\u003c/sup\u003e. The spectral response of soils varied with salt content, with high-salinity soils exhibiting stronger responses in the visible and near-infrared bands than those with lower salinity\u003csup\u003e\u003cspan citationid=\"CR3\" class=\"CitationRef\"\u003e3\u003c/span\u003e\u003c/sup\u003e. Leveraging remote sensing for the dynamic monitoring of soil salinization is crucial for the efficient utilization of soil and water resources, providing essential insights into the timing, patterns, and locations of potential changes in soil salinity, thus facilitating improved resource management and planning.\u003c/p\u003e \u003cp\u003eIn recent years, Unmanned Aerial Vehicle (UAV) and other aerial remote sensing platforms have experienced rapid advancements, increasingly integrating into civilian applications and becoming prominent in agricultural research. UAVs boast advantages like portability, high flexibility, and customizable flight durations. Zhao et al.\u003csup\u003e\u003cspan citationid=\"CR4\" class=\"CitationRef\"\u003e4\u003c/span\u003e\u003c/sup\u003e used multispectral remote sensing data from three research locations to establish soil salinity inversion models based on support vector machine (SVM), RF, backpropagation neural network (BPNN), and extreme learning machine (ELM). The results showed that the four spectral index-based models could get a high level of inversion accuracy. Similarly, Wei et al.\u003csup\u003e\u003cspan citationid=\"CR5\" class=\"CitationRef\"\u003e5\u003c/span\u003e\u003c/sup\u003e used a UAV equipped with Micro-MCA multispectral sensors to capture pictures for evaluating the soil salinity in a small area of the Hetao Irrigation District. Chen et al.\u003csup\u003e\u003cspan citationid=\"CR6\" class=\"CitationRef\"\u003e6\u003c/span\u003e\u003c/sup\u003e developed estimation models of SSC in sunflower fields at different soil depths during the budding and blooming stages using UAV multispectral data. IVUSHKIN et al.\u003csup\u003e\u003cspan citationid=\"CR7\" class=\"CitationRef\"\u003e7\u003c/span\u003e\u003c/sup\u003e found that UAVs equipped with multiple sensors of hyperspectral, multispectral, thermal infrared, and LiDAR cameras have great potential for monitoring soil salinization. Feature indices like the soil salinity index and vegetation index, derived from spectral band reflectance transformations, serve as important variables for estimating soil salinity content\u003csup\u003e\u003cspan citationid=\"CR8\" class=\"CitationRef\"\u003e8\u003c/span\u003e\u003c/sup\u003e. Qi et al.\u003csup\u003e\u003cspan citationid=\"CR9\" class=\"CitationRef\"\u003e9\u003c/span\u003e\u003c/sup\u003e collected spectrum reflectance and spectral indices based on a UAV ground cooperative system, and then used machine learning algorithms like BPNN to build a salinity inversion model. The findings indicated that the constructed model accurately captured the level of salinization in the research area.\u003c/p\u003e \u003cp\u003eThe relationship between soil salinity and remote sensing feature variables was frequently nonlinear due to the combination of complex factors such as soil, vegetation, and atmospheric signatures\u003csup\u003e\u003cspan citationid=\"CR10\" class=\"CitationRef\"\u003e10\u003c/span\u003e\u003c/sup\u003e. A common method for estimating SSC was to use the mathematical statistical model, especially the linear regression model including partial least squares regression (PLSR)\u003csup\u003e\u003cspan citationid=\"CR11\" class=\"CitationRef\"\u003e11\u003c/span\u003e,\u003cspan citationid=\"CR12\" class=\"CitationRef\"\u003e12\u003c/span\u003e\u003c/sup\u003e. But the relationship between spectral covariates and soil properties in nature was rarely linear\u003csup\u003e\u003cspan citationid=\"CR13\" class=\"CitationRef\"\u003e13\u003c/span\u003e\u003c/sup\u003e. Machine learning (ML) algorithms, renowned for their ability to address nonlinear problems and process high dimensional data, often outperform statistical regression models in predicting soil salinity. Algorithms like RF、SVM and BPNN could capture these intricate nonlinear relationships through multi-layered results, making them particularly effective for soil salinity estimation\u003csup\u003e\u003cspan citationid=\"CR14\" class=\"CitationRef\"\u003e14\u003c/span\u003e\u003c/sup\u003e. Hu\u003csup\u003e\u003cspan citationid=\"CR15\" class=\"CitationRef\"\u003e15\u003c/span\u003e\u003c/sup\u003e compared PLSR and RF methods using hyperspectral first order differentiation, broad band spectral indices, and narrow band spectral indices as independent variables. It was reported that the RF model could better predict soil salinity and the model for the bare soil area had the highest prediction accuracy compared to the other sample areas. Zhang et al.\u003csup\u003e\u003cspan citationid=\"CR16\" class=\"CitationRef\"\u003e16\u003c/span\u003e\u003c/sup\u003e found that the SVM method based on texture features had an improved effect on the monitoring accuracy of soil salinity information. Wei et al.\u003csup\u003e\u003cspan citationid=\"CR5\" class=\"CitationRef\"\u003e5\u003c/span\u003e\u003c/sup\u003e established various machine learning models to evaluate and optimize the best salt estimate model. However, the prediction accuracy of a single machine learning algorithm could fluctuate under varying conditions. Therefore, evaluating the performance of multiple machine learning regression algorithms is essential for constructing reliable soil salinity prediction models that could adapt to diverse environmental factors.\u003c/p\u003e \u003cp\u003eThere is a strong correlation between vegetation growth and soil salinity, as evidenced spectrally in two primary ways: differences in leaf spectral reflectance and significant variations in the texture features of spectral images influenced by changes in soil salinity or leaf characteristics\u003csup\u003e\u003cspan citationid=\"CR17\" class=\"CitationRef\"\u003e17\u003c/span\u003e,\u003cspan citationid=\"CR18\" class=\"CitationRef\"\u003e18\u003c/span\u003e\u003c/sup\u003e Studies have shown that texture features are extensively utilized to reveal vegetation characteristics variations, combining them with spectral information could effectively improve the accuracy of predictive models\u003csup\u003e\u003cspan citationid=\"CR19\" class=\"CitationRef\"\u003e19\u003c/span\u003e,\u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e20\u003c/span\u003e,\u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e21\u003c/span\u003e,\u003cspan citationid=\"CR22\" class=\"CitationRef\"\u003e22\u003c/span\u003e,\u003cspan citationid=\"CR23\" class=\"CitationRef\"\u003e23\u003c/span\u003e\u003c/sup\u003e. Huang et al.\u003csup\u003e\u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e24\u003c/span\u003e\u003c/sup\u003e used Sentinel-2 image combined with texture features to significantly improve the classification effect of moderate saline soil in the Yellow River Delta region. Nevertheless, the prevalent use of remote sensing data for soil salinity detection primarily relied on spectral indices\u003csup\u003e\u003cspan citationid=\"CR5\" class=\"CitationRef\"\u003e5\u003c/span\u003e\u003c/sup\u003e, with limited studies investigating the use of texture features in soil salinity estimation\u003csup\u003e\u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e24\u003c/span\u003e\u003c/sup\u003e. In the case of bare soil, the spectrum could directly determine the salt content of the soil surface. Conversely under vegetation coverage, soil salinity could be indirectly assessed by vegetation canopy spectral signatures\u003csup\u003e\u003cspan citationid=\"CR25\" class=\"CitationRef\"\u003e25\u003c/span\u003e\u003c/sup\u003e.\u003c/p\u003e \u003cp\u003eFew investigations on estimating soil salinity under vegetation cover have been conducted in the past. Therefore, the purpose of this study is to construct soil salinity estimation models at various depths beneath barley growth. The main research contents were to:(1) evaluate the potential and feasibility of UAV remote sensing for SSC estimation under barley coverage; (2) optimize the feature variables of spectral band, spectral index and texture data based on BDT method; (3) Validate the accuracy of soil salt estimation model for different vegetation coverage conditions. This study presented a novel method for the dynamic monitoring of soil salinity in agricultural production helped to achieving precise irrigation and fertilization practices, and improved agricultural productivity efficiency and sustainable utilization of soil resources.\u003c/p\u003e"},{"header":"2. Materials and methods","content":"\u003cdiv id=\"Sec3\" class=\"Section2\"\u003e \u003ch2\u003e2.1 Study area and experimental design\u003c/h2\u003e \u003cp\u003eThe study area is located at the ecological experimental station of Yangzhou University, situated in the Jianghuai Plain of Jiangsu Province in eastern China (119\u0026deg;24\u0026prime;E, 32\u0026deg;21\u0026prime;N), at an altitude of 5 meters. This region is characterized by a subtropical monsoon climate, with an average annual frostfree period, precipitation, evaporation, and air temperature of 223 days 937 mm, 1063 mm, and 14.8\u0026deg;C, respectively. Soil samples for the experiment were collected from Tiaozini reclamation area in Dongtai City, Jiangsu Province. This area was subject to marine intrusion and groundwater topdressing, with high salinity in the soil tillage layer.\u003c/p\u003e \u003cp\u003eThe crop studied was barley (Hordeum vulgare L.). Four different soil salinity treatments were set: control (no salt), low salinity (3\u0026permil;), medium salinity (5\u0026permil;), and high salinity (10\u0026permil;). The salinity experiment was carried out using a box planting setup with dimensions of 100 cm in length, 40 cm in width, and 40 cm in height, each containing 120 kg of base soil. Each treatment was replicated thrice, resulting in a total of 12 experimental units. The cultivation adhered to local management practices, encompassing weeding, pest, and disease prevention and control. Barley was planted annually at the early of November, with each box accommodating two rows and was harvested at the late of May of the following year. To maintain the soil salinity status of each treatment, the experimental box devices were designed as sealed containers to prevent salt leaching by precipitation and irrigation.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec4\" class=\"Section2\"\u003e \u003ch2\u003e2.2 Data collection and acquisition\u003c/h2\u003e \u003cdiv id=\"Sec5\" class=\"Section3\"\u003e \u003ch2\u003e2.2.1 Soil salt measurement\u003c/h2\u003e \u003cp\u003eDuring the barley growth stages of reviving-jointing, jointing-filling, and grain filling maturity, soil electrical conductivity data were measured every 7\u0026ndash;10 days, in conjunction with the acquisition of multispectral imagery. A total of 84 dataset were collected during barley growth. The soil electrical conductivity was measured using the EC450 conductivity meter (Spectrum Technologies Co., Ltd., Chicago, IL, USA). First, the electrode was calibrated using a calibration solution (conductivity: 1413 \u0026micro;S/cm). After calibration, the electrode was inserted into the soil profile to measure conductivity at soil depths of 0ཞ10 cm, 10ཞ20 cm, 20ཞ30 cm, and 30ཞ40 cm. The soil conductivity values were directly recorded by the handheld meter. Subsequently, the SSC was derived using the empirical formula (SSC\u0026thinsp;=\u0026thinsp;0.2882EC\u0026thinsp;+\u0026thinsp;0.0183, %).\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec6\" class=\"Section3\"\u003e \u003ch2\u003e2.2.2 Multispectral data acquisition and processing\u003c/h2\u003e \u003cp\u003eThis study leveraged remote sensing data from multiple UAV based sensors to enhance monitoring of soil properties through data fusion. The DJI Inspire 2 UAV platform (DJI Inc., Shenzhen, China) used in this study, which was equipped with a Altum multi-spectral and infrared camera (MicaSense, Inc., Seattle, WA, USA). This camera was capable of capturing images across six spectral bands (blue, green, red, red edge, near-infrared, and thermal infrared) simultaneously. Variables The UAV flight operations were conducted under optimal conditions of clear skies and calm winds, between 11 a.m. and 2 p.m. local time. To ensure high-quality data, the sensor was prewarmed for 5 minutes and the reference plate was used for radiometric calibration before each flight. The flight altitude and the cross-track overlap were maintained to 25 meters and 75%, respectively. The camera was oriented vertically downward to achieve a ground resolution of 1.1 cm per pixel. After preprocessing the collected multispectral images by radiation correction, geometric correction, and image mosaic, the reflectance data of each pixel were generated during crop growth period. To extract canopy reflectance, a region of interest (ROI) was preset near the center of each plot.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec7\" class=\"Section3\"\u003e \u003ch2\u003e2.2.3 Vegetation coverage calculation\u003c/h2\u003e \u003cp\u003eAfter collecting multispectral imagery, the canopy light interception of barley was measured using the AccuPAR LP-80 ceptometer (Decagon Devices Inc, Pullman, WA, USA). The LP-80 is equipped with 80 independent sensors, each spaced uniformly at 1 cm intervals, capable of measuring solar radiation within the 400\u0026ndash;700 nm band in different modes. During the measurement, the sensor was strategically positioned in the center of the plot, aligned with the row direction. With the acquired photosynthetically active radiation at the top and bottom of the canopy, the instrument automatically calculated the crop LAI through a built-in algorithm with the other variables. In accordance with the reports and employing a classification method that considers similar geographical features and vegetation types\u003csup\u003e\u003cspan citationid=\"CR26\" class=\"CitationRef\"\u003e26\u003c/span\u003e\u003c/sup\u003e, vegetation cover was classified according to the LAI value as follows: low coverage (0, 0.45), medium coverage (0.45, 0.75) and high coverage (0.75, 1).\u003c/p\u003e \u003c/div\u003e \u003c/div\u003e \u003cdiv id=\"Sec8\" class=\"Section2\"\u003e \u003ch2\u003e2.3 Construction of feature variables\u003c/h2\u003e \u003cp\u003eSpectral indices were classical variables that integrated the spectral characteristics of each band of ground objects and enhanced specific information through mathematical transformations and combinations of reflectance values from different bands. Salinity indices, which were frequently employed for the rapid assessment of soil salinization, exhibited a strong correlation with bare soil salinity\u003csup\u003e\u003cspan citationid=\"CR27\" class=\"CitationRef\"\u003e27\u003c/span\u003e\u003c/sup\u003e. Similarly, vegetation indices were frequently used for the quantitative assessment of vegetation growth\u003csup\u003e\u003cspan citationid=\"CR28\" class=\"CitationRef\"\u003e28\u003c/span\u003e\u003c/sup\u003e. In this study, a selection of widely recognized spectral indicators were employed for the monitoring of soil salinization, including 15 salinity indices and 15 vegetation indices.\u003c/p\u003e \u003cp\u003eThe selected salinity indices include: Normalized Difference Salinity Index (NDSI), R-edge Normalized Difference Salinity Index (NDSI-reg), Brightness Index (BI), Salinity Index 1 (SI1), R-edge Salinity Index 1 (SI1-reg), Salinity Index 2 (SI2), R-edge Salinity Index 2 (SI2-reg), Salinity Index 3 (SI3), R-edge Salinity Index 3 (SI3-reg), Salinity Index (SI-T), Salinity Index S1, Salinity Index S2, Salinity Index S3, Salinity Index S5, Salinity Index SI; The selected vegetation indices include: Normalized Difference Vegetation Index (NDVI), R-edge Normalized Difference Vegetation Index (NDVI-reg), Difference Vegetation Index (DVI), R-edge Difference Vegetation Index (DVI-reg), Enhanced Vegetation Index (EVI), R-edge Difference Vegetation Index (EVI-reg), Triangular Vegetation Index (TVI), Normalized Greenness Index (NDGI), Simple Ratio Index (SR), Modified Soil Adjusted Vegetation Index (MSAVI), Optimized Soil Adjusted Vegetation Index (OSAVI), Soil Adjusted Vegetation Index (SAVI), Visible Light Band Difference Vegetation Index (VDVI), Visible Atmospherically Resistant Index (VARI), Green Normalized Difference Vegetation Index (GNDVI), and a normalized relative canopy temperature (NRCT). The calculation formulas for these indices were shown in Table\u0026nbsp;\u003cspan refid=\"Tab1\" class=\"InternalRef\"\u003e1\u003c/span\u003e.\u003c/p\u003e \u003cp\u003e \u003cdiv class=\"gridtable\"\u003e\u003ctable float=\"Yes\" id=\"Tab1\" border=\"1\"\u003e \u003ccaption language=\"En\"\u003e \u003cdiv class=\"CaptionNumber\"\u003eTable 1\u003c/div\u003e \u003cdiv class=\"CaptionContent\"\u003e \u003cp\u003eSpectral index and calculation formula\u003c/p\u003e \u003c/div\u003e \u003c/caption\u003e \u003ccolgroup cols=\"3\"\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c1\" colnum=\"1\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c2\" colnum=\"2\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c3\" colnum=\"3\"\u003e\u003c/div\u003e \u003cthead\u003e \u003ctr\u003e \u003cth align=\"left\" colname=\"c1\"\u003e \u003cp\u003eFeature\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c2\"\u003e \u003cp\u003eFormulation\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c3\"\u003e \u003cp\u003eReference\u003c/p\u003e \u003c/th\u003e \u003c/tr\u003e \u003c/thead\u003e \u003ctbody\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eNDSI\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eNDSI= \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\frac{R-NIR}{R+NIR}\\)\u003c/span\u003e\u003c/span\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e29\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eNDSI-reg\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eNDSI-reg=\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\frac{RedEdge-NIR}{RedEdge+NIR}\\)\u003c/span\u003e\u003c/span\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e30\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eBI\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eBI=\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\sqrt{{R}^{2}+{NIR}^{2}}\\)\u003c/span\u003e\u003c/span\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e31\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eSI1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eSI1=\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\sqrt{G\\times\\:R}\\)\u003c/span\u003e\u003c/span\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e28,32\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eSI1-reg\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eSI1-reg=\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\sqrt{G\\times\\:RedEdge}\\)\u003c/span\u003e\u003c/span\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e30\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eSI2\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eSI2=\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\sqrt{{G}^{2}+{R}^{2}+{NIR}^{2}}\\)\u003c/span\u003e\u003c/span\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e28,32\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eSI2-reg\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eSI2-reg=\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\sqrt{{G}^{2}+{RedEdge}^{2}+{NIR}^{2}}\\)\u003c/span\u003e\u003c/span\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e30\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eSI3\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eSI3=\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\sqrt{{G}^{2}+{R}^{2}}\\)\u003c/span\u003e\u003c/span\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e28,32\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eSI3-reg\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eSI3-reg=\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\sqrt{{G}^{2}+{RedEdge}^{2}}\\)\u003c/span\u003e\u003c/span\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e30\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eSI-T\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eSI-T=\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:100\\times\\:R/NIR\\)\u003c/span\u003e\u003c/span\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e33\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eS1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eS1=\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:B/R\\)\u003c/span\u003e\u003c/span\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e34\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eS2\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eS2=\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\frac{B-R}{B+R}\\)\u003c/span\u003e\u003c/span\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e34\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eS3\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eS3=\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:G\\times\\:R/B\\)\u003c/span\u003e\u003c/span\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e34\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eS5\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eS5=\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:B\\times\\:R/G\\)\u003c/span\u003e\u003c/span\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e34\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eSI\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eSI=\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\sqrt{B\\times\\:R}\\)\u003c/span\u003e\u003c/span\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e35\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eNDVI\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eNDVI=\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\frac{NIR-R}{NIR+R}\\)\u003c/span\u003e\u003c/span\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e36\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eNDVI-reg\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eNDVI-reg=\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\frac{NIR-RedEdge}{NIR+RedEdge}\\)\u003c/span\u003e\u003c/span\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e30\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eDVI\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eDVI= \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:NIR-R\\)\u003c/span\u003e\u003c/span\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e28\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eDVI-reg\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eDVI-reg= \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:NIR-RedEdge\\)\u003c/span\u003e\u003c/span\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e30\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eEVI\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eEVI=\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:2.5\\times\\:\\frac{NIR-R}{NIR+6R-7.5B+1}\\)\u003c/span\u003e\u003c/span\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e37\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eEVI-reg\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eEVI-reg=\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:2.5\\times\\:\\frac{NIR-RedEdge}{NIR+6RedEdge-7.5B+1}\\)\u003c/span\u003e\u003c/span\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e30\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eTVI\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eTVI=\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:0.5\\times\\:\\frac{120\\times\\:\\left(NIR-G\\right)}{200\\times\\:\\:\\left(R-G\\right)}\\)\u003c/span\u003e\u003c/span\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e9\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eNDGI\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eNDGI=\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\frac{G-R}{G+R}\\)\u003c/span\u003e\u003c/span\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e38\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eSR\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eSR=\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:NIR/R\\)\u003c/span\u003e\u003c/span\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e39\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eMSAVI\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eMSAVI=\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\frac{\\left(2NIR-1\\right)-\\sqrt{{\\left(2NIR+1\\right)}^{2}-8\\left(NIR-R\\right)}}{2}\\)\u003c/span\u003e\u003c/span\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e30\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eOSAVI\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eOSAVI=\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:1.16\\times\\:\\frac{NIR-R}{NIR+R+0.16}\\)\u003c/span\u003e\u003c/span\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e30\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eSAVI\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eSAVI=\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:1.5\\times\\:\\frac{NIR-R}{NIR+R+0.5}\\)\u003c/span\u003e\u003c/span\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e40\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eVDVI\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eVDVI=\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\frac{2G-\\left(R+B\\right)}{2G+\\left(R+B\\right)}\\)\u003c/span\u003e\u003c/span\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e41\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eVARI\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eVARI=\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\frac{G-R}{G+R+B}\\)\u003c/span\u003e\u003c/span\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e42\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eGNDVI\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eGNDVI=\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\frac{NIR-G}{NIR+G}\\)\u003c/span\u003e\u003c/span\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e9\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eNRCT\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eNRCT=\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\frac{Ti-Tmin}{Tmax-Tmin}\\)\u003c/span\u003e\u003c/span\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e43\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003c/tbody\u003e \u003c/colgroup\u003e \u003ctfoot\u003e \u003ctr\u003e\u003ctd colspan=\"3\"\u003eNote: B, R, G, NIR, and Rededge represent the reflectance of blue, red, green, near-infrared, and red edge bands, respectively.\u003c/td\u003e\u003c/tr\u003e \u003c/tfoot\u003e \u003c/table\u003e\u003c/div\u003e \u003c/p\u003e \u003cp\u003eThe image texture was indicative of variations in soil surface\u0026rsquo;s color and gray level, which was closely related to the soil's salt and water content. In this study, the statistical Gray Level Cooccurrence Matrix (GLCM) method was used to extract the textural features of the images. GLCM is a prevalent method widely used for image feature extraction, texture analysis, and quality evaluation\u003csup\u003e\u003cspan citationid=\"CR44\" class=\"CitationRef\"\u003e44\u003c/span\u003e\u003c/sup\u003e, describing the correlation between pixel gray levels within images. Eight characteristic variables including mean (MEA), variance (VAR), uniformity (HOM), contrast (CON), difference (DIS), entropy (ENT), second-order moment (SEC) and correlation (COR) were calculated for each band of UAV multispectral images utilizing the second-order statistical filtering tool. Variables In total, this study incorporated 5 band reflectance values, canopy temperature, 31 spectral indices, and 40 texture feature data, amounting to 77 feature variables serving as independent variables.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec9\" class=\"Section2\"\u003e \u003ch2\u003e2.4 Feature variable selection and set construction\u003c/h2\u003e \u003cdiv id=\"Sec10\" class=\"Section3\"\u003e \u003ch2\u003e2.4.1 Optimization of characteristic variables\u003c/h2\u003e \u003cp\u003eTo optimize the complexity of the model input variables, this study employed the BDT method to refine the selection of the 77 characteristic variables. BDT is one of the embedded type methods which used to compute estimates of Predictor Importance for the tree by summing changes in the mean squared error (MSE) due to splits on every predictor and dividing the sum by the number of branch nodes\u003csup\u003e\u003cspan citationid=\"CR45\" class=\"CitationRef\"\u003e45\u003c/span\u003e\u003c/sup\u003e. The selection steps involved in this optimization process include: (1) the data were divided into twenty groups; (2) importance score of feature variables in each set of data was calculated based on BDT; (3) normalizing the important scores for various data categories. (4) A statistical analysis and sorting were performed on the normalized data of 20 groups. The method used in this study was as follows: The obtained 84 datasets were randomly divided, with 70% allocated to the training subset and the remaining 30% reserved for the test subset. This process was repeated 20 times. The results from each iteration were normalized across different groups (spectral band reflectance group, spectral index group, and texture data group). Based on the contribution degree of characteristic variables within each group, the threshold was set to 0.1 to select the pivotal variables for model training and establishing.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec11\" class=\"Section3\"\u003e \u003ch2\u003e2.4.2 Construction of model training\u003c/h2\u003e \u003cp\u003eThe construction of a predictive model necessitates the careful extraction of input features and the judicious selection of appropriate machine learning algorithms. To explore the influence of different features on model performance, the feature combinations were categorized as follows:\u003c/p\u003e \u003cp\u003eSingle variable groups: (1) Band reflectance. Composed of the original reflectance values of each spectral band. These values were directly reflected the physical and chemical properties of the soil, serving as the primary features for predicting SSC. (2) Spectral Index. Consists of various spectral indices (such as NDVI, RVI, etc.) calculated from different spectral bands. These indices were designed to enhance specific soil or vegetation characteristics. (3) Texture data. Extracting from spectral data, such as the gray level cooccurrence matrix. These features captured subtle changes and spatial structure of the soil surface.\u003c/p\u003e \u003cp\u003ePairwise variable groups: (1) Spectral band reflectance and indices combination: This combination retained the original spectral information while incorporating enhanced information, potentially providing more comprehensive responses of soil salinity changes. (2) Spectral band reflectance and texture feature combination: By integrating the complementary advantages of both types of data, this combination revealed spectral characteristics and captured spatial heterogeneity of the soil surface. (3) Spectral indices and texture feature combination: This combination enriched the enhanced information with spatial structure data, potentially improving the model's generalization.\u003c/p\u003e \u003cp\u003eFull Variable Group: This group contained all types data of spectral bands reflectance, spectral indices, and texture features. The comprehensive inclusion of all potentially relevant information was intended to enhance the model's generalization capabilities and predictive performance.\u003c/p\u003e \u003c/div\u003e \u003c/div\u003e \u003cdiv id=\"Sec12\" class=\"Section2\"\u003e \u003ch2\u003e2.5 Machine learning models\u003c/h2\u003e \u003cp\u003eBased on the SSC and the selected feature variables, four machine learning methods of RF, SVM, GPR, and BPNN methods was used to construct soil salt estimation models.\u003c/p\u003e \u003cp\u003eRF was an ensemble learning method that used to establish prediction models for classification and regression problems by using randomly unrelated decision trees\u003csup\u003e\u003cspan citationid=\"CR46\" class=\"CitationRef\"\u003e46\u003c/span\u003e\u003c/sup\u003e. In recent years, RF has been widely used in vegetation growth parameters estimate and the estimation of soil physical and chemical parameters. For instance, Huang et al.\u003csup\u003e\u003cspan citationid=\"CR47\" class=\"CitationRef\"\u003e47\u003c/span\u003e\u003c/sup\u003e established several soil salinity estimation models based on Landsat-8 OLI images in a study of oasis soil salinity in arid areas, and found that the estimation accuracy of the RF modeling method was higher than that of classical statistical models. Similarly, Sui et al.\u003csup\u003e\u003cspan citationid=\"CR48\" class=\"CitationRef\"\u003e48\u003c/span\u003e\u003c/sup\u003e developed a soil salinity estimation model based on original observations and satellite data in a study of coastal soil salinity, utilizing hydrological connectivity measurements and the RF algorithm.\u003c/p\u003e \u003cp\u003eSVM was a method that implemented the concept of structural risk minimization, effectively addressing small sample sizes, nonlinearity, and high-dimensional data. SVM offered strong expression capability, generalization ability, and learning efficiency. It easily integrated with multi-source information, thereby achieving higher estimation accuracy\u003csup\u003e\u003cspan citationid=\"CR49\" class=\"CitationRef\"\u003e49\u003c/span\u003e\u003c/sup\u003e. For example, Cai et al. (2010) combined multispectral and texture features and used the SVM classifier to identify soil affected by salt, confirming that the SVM classifier effectively extracted soil salinization distribution information in the Yinchuan Plain. Similarly, Guan et al.\u003csup\u003e\u003cspan citationid=\"CR50\" class=\"CitationRef\"\u003e50\u003c/span\u003e\u003c/sup\u003e introduced SVM theory into the dynamic prediction of soil EC value, constructing a dynamic prediction model for soil salinity aimed at managing irrigation water in salinized areas.\u003c/p\u003e \u003cp\u003eGPR was a nonparametric Bayesian regression method that predicted by assuming the data distributed as a multivariate Gaussian distribution, providing an estimate of prediction uncertainty. It was flexible and efficient for small datasets but can be computationally expensive and difficult to scale for larger ones. Additionally, selecting and adjusting the appropriate kernel function required considerable experience and testing.\u003c/p\u003e \u003cp\u003eBPNN was a feedforward network composed of multiple neurons capable of learning and identifying nonlinear relationships in complex systems. It exhibited strong self-learning ability, adaptability, and resistance to interference, making it highly promising for the estimation of soil physical and chemical parameters. For instance, Wang et al.\u003csup\u003e\u003cspan citationid=\"CR51\" class=\"CitationRef\"\u003e51\u003c/span\u003e\u003c/sup\u003e successfully established a prediction model for soil moisture and salinity using BPNN with Landsat-8 satellite data.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec13\" class=\"Section2\"\u003e \u003ch2\u003e2.6 Technical workflow\u003c/h2\u003e \u003cp\u003eThis study aims to identify the most effective variables set for accurately predicting SSC, thereby improving the overall performance and applicability of the model in diverse soil depths. SSC was designated as the dependent variable, while different variable groups serving as the independent variables. The sample data were randomly divided into two groups, with 70% allocated for model training and 30% for validation. Four different machine learning methods (RF, SVM, GPR, and BPNN) were employed to estimate SSC. Each method constructed corresponding prediction models, optimizing the performance by adjusting the model hyperparameters. To further evaluate model accuracy, 10-fold cross-validation was used to construct and test the SSC estimation models (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003e).\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e"},{"header":"3. Results","content":"\u003cdiv id=\"Sec15\" class=\"Section2\"\u003e \u003ch2\u003e3.1 Statistical of soil salinity distribution\u003c/h2\u003e \u003cp\u003eThe salt content of sampling points at different soil depths was categorized as follows: non-saline soil (\u0026lt;\u0026thinsp;0.2%), mild salinization (salt content 0.2\u0026thinsp;~\u0026thinsp;0.5%), severe salinization (0.5\u0026thinsp;~\u0026thinsp;1.0%), and saline soil (\u0026gt;\u0026thinsp;1.0%). The non-salt treatment had an average salt content of 0.214%, the low-salt treatment had an average salt content of 0.257%, the medium-salt treatment had an average salt content of 0.418%, and the high-salt treatment had an average value of 1.353%. The statistical of obtained salt content data were as following: in the 0 to 10 cm soil depth, the measured SSC were relatively lower than those of other depths, which varied from 0.031% to1.04%, with an average of 0.428%. At a depth of 10 to 20 cm, the SSC ranged from 0.087\u0026ndash;1.406%, with an average of 0.561%. The SSC in the 20 to 30 cm depth were the highest, with the average SSC of 0.633%, varying between 0.055\u0026thinsp;~\u0026thinsp;1.806%. The SSC distribution in 30\u0026thinsp;~\u0026thinsp;40 cm was similar to that in the 10\u0026thinsp;~\u0026thinsp;20 cm soil layer, with an average of 0.619%. According to these criteria, the measured salt grade distribution in the study area was shown in Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003e.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec16\" class=\"Section2\"\u003e \u003ch2\u003e3.2 Contribution and selection of feature variables\u003c/h2\u003e \u003cp\u003eIn this study, the important factor of 77 feature variables of three types data for different soil depths were calculated (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003e). After normalizing the feature variables important factor of different groups and setting a threshold of 0.1, the feature variables were screened. Common variables across different depths were selected for the ML model training. Finally, the selected variables for different groups were as follow: Band reflectance: B, R-edge, and NIR; spectral Index: SI2, SI2-reg, BI, SI-T, DVI, EVI, SAVI, OSAVI, DVI-reg, EVI-reg, and MSAVI; Texture Data B-MEA, B-VAR, B-CON, B-ENT, B-SEC, G-VAR, RE-MEA, NIR-MEA, NIR-HOM, NIR-ENT, and NIR-SEC. Altogether, 25 feature variables were filtered for modeling.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec17\" class=\"Section2\"\u003e \u003ch2\u003e3.3 Comparision of model performance for various ML methods and data groups\u003c/h2\u003e \u003cp\u003eThe optimal variables of the three variables types of band reflectance, spectral index, and texture data and their combinations were used as the independent variables and the SSC as the target variable imputing into the machine learning (RF, SVM, GPR, and BPNN) model to establish the soil salt prediction model. The prediction accuracies of these models based on different variable combinations were shown in Figs.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003e and \u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003e.\u003c/p\u003e \u003cp\u003eFigures\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003e and \u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003e demonstrated that for each algorithm, the band reflectance and spectral index variables performed well in estimating SSC at a depth of 0 to 10 cm within the single variable groups. Among these, the BPNN model, based on the spectral index, achieved the best performance with an R\u0026sup2; of 0.74 and RMSE of 0.26%. The whole variable combination-based GPR model was the best in the multivariate combination group. It exhibited a stableR\u0026sup2; value around 0.77, an RMSE of 0.24%, and great overall stability. The RF model identified the band reflectance as the most single variable for SSC at a depth of 10 to 20 cm, whereas the SVM model produced the worst estimation. Across the four approaches tested, texture data yielded suboptimal results, with an average RMSE exceeding 0.34%. However, The RF model that incorporated both texture and spectral index data outperformed others in the multivariate combination, achieving a mean RMSE of 0.31% and an R\u0026sup2; of approximately 0.61, making it the most effective model overall. Additionally, the GPR model's results outperformed those of other approaches. In contrast, the SVM and BPNN models consistently delivered less accurate predictions across all groups.\u003c/p\u003e \u003cp\u003eAt a depth of 20 to 30 cm, the RF model built using band reflectance and spectral index demonstrated superior performance, with average R\u0026sup2; values of 0.62 and 0.67 and mean RMSE of 0.31% and 0.29%, respectively. This model outperformed than the models built by the other three algorithm. Within the multivariate combination, the spectral index and texture data performed better outcomes in the RF model compared to other combinations, yielding an R\u0026sup2; of 0.67 and an RMSE of 0.29%. Similarly, the GPR algorithm, when applied to the texture data showed higher accuracy than the other three algorithms. The SVM model achieved its highest accuracy when all variables were combined, while the accuracy of models based on band reflectance and spectral index data in the GPR and BPNN algorithms was marginally higher than that of other combinations, with R\u0026sup2; values of 0.64 and 0.58, and RMSE values of 0.30% and 0.35%, respectively.\u003c/p\u003e \u003cp\u003eThe RF model built using textural data as a single variable had the lowest accuracy at a depth of 30 to 40 cm, with an average RMSE of 0.39% and R\u0026sup2; of 0.37. In contrast, the RF model constructed using the spectral index showed higher accuracy than the other three models. Within the multi-variable combination group, the SVM model retained consistent accuracy, with the R\u0026sup2; of approximately 0.55. The GPR model, when applied to the full variable set, demonstrated superior accuracy compared to other combinations. Notably, the BPNN model exhibited better accuracy when based on band reflectance and spectral index data than other models.\u003c/p\u003e \u003cp\u003eIn conclusion, soil salinity estimation at different depths were improved by applying a variety of machine learning techniques, with the effectiveness of these methods differing across modeling groups. The most effective method for estimating soil salinity was the whole variable group based GPR model, particularly for both surface (0ཞ10 cm) and deep (30ཞ40 cm) measurements. However, for middle depth soils (10\u0026ndash;20 cm and 20\u0026ndash;30 cm), the RF model based on spectral index and texture data displayed the best outcomes. The RF model effectively reduced the variance of a single tree by integrating multiple decision trees, thereby improving the model's stability and estimation accuracy. While the SVM and BPNN models performed marginally worse overall than GPR and RF, they were still fairly good at predicting within particular categories.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec18\" class=\"Section2\"\u003e \u003ch2\u003e3.4 Outer-sample validation and analysis\u003c/h2\u003e \u003cp\u003eThe performances of optimal estimation models have been evaluated 20 times using the outer-samples (not training) data (Fig.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003e). At a soil depth of 0ཞ10 cm, the R\u0026sup2; values of the validation set fluctuated between 0.6 and 0.9, with an average R\u0026sup2;of 0.77. The RMSE values typically ranged from 0.1 to 0.3%, with an average RMSE of 0.185%. This indicated that the GPR model had high prediction accuracy and stability in predicting shallow soil salinity. For the 10ཞ20 cm depth, the RF model showed favorable results, with an average R\u0026sup2; of 0.67. The RMSE values ranged between 0.2 and 0.4%, with an average RMSE of 0.325%. At the 20\u0026thinsp;~\u0026thinsp;30 cm depth, the RF model achieved a maximum R\u003csup\u003e\u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2\u003c/span\u003e\u003c/sup\u003e of 0.87 and minimum RMSE of 0.212%. At a deepest soil layer, the model's R\u003csup\u003e\u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2\u003c/span\u003e\u003c/sup\u003e ranged from a minimum of t 0.52 to a maximum of 0.92, indicating high overall accuracy and stability. The RMSE values fluctuated between 0.169% and 0.485%, with an average RMSE of 0.31%.\u003c/p\u003e \u003cp\u003eThe modeling and validation accuracies of four soil salinity estimation models at different depths were compared. Models were selected at each depth to achieve better estimation results in both modeling and validation. The average R\u003csup\u003e\u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2\u003c/span\u003e\u003c/sup\u003e values in the modeling set at each depth exceeded 0.6, and the average RMSE were below 0.4%, indicating that the model had high accuracy and good stability.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec19\" class=\"Section2\"\u003e \u003ch2\u003e3.5 Performance of salinity estimations under different vegetation coverage\u003c/h2\u003e \u003cp\u003eIn this study, the LAI was used as a proxy for vegetation cover, which is divided into three different levels: low, medium and high. Correspondingly, soil salinity data are stratified based on these cover classes. The results of SSC estimation under various covers were presented in Fig.\u0026nbsp;\u003cspan refid=\"Fig7\" class=\"InternalRef\"\u003e7\u003c/span\u003e.\u003c/p\u003e \u003cp\u003eThe accuracy of the salt estimation models at various depths was assessed based on the salt data collected under different vegetation coverage levels, and corresponding scatter plots were created. As shown in Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003e, all models achieved R\u0026sup2; values greater than 0.4, with the highest value reaching 0.83, indicating that the model's accuracy was highest for soil surface with intermediate vegetation covering. The models demonstrated good accuracy for soil depths of 0 to 10 cm under each vegetation coverage level, with R\u0026sup2; values exceeding 0.7 and low RMSE values. At a depth of 10 to 20 cm, the model accuracy was highest under low vegetation covering, while the accuracy under medium vegetation coverage was lowest. However, at a depth of 20 to 30 cm, the maximum validation accuracy was achieved with medium coverage, with R\u0026sup2; values exceeding 0.7 across all three coverages levels. With considerable vegetation coverage, the maximum accuracy at 30 to 40 cm was achieved, with an RMSE of 0.241% and an R\u0026sup2; of 0.78. Generally, when using the data selected under each vegetation coverage level as the validation set, the favored models at different depths showed good estimation performance.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e"},{"header":"4. Discussion","content":"\u003cp\u003eSoil salinity estimation by UAV multispectral remote sensing data plays a positive role in salinity monitoring and management. In this study, the spectral data and salinity data acquired during the vegetation cover period were used to construct and validate the estimation model for soil salinity at different depths. In addition, the model performances were validated under different vegetation cover ratios and the impact of cover on salt estimation in different soil layers was also explored.\u003c/p\u003e \u003cp\u003eFeature selection is an important part of constructing the soil salinity estimation model. Through the BDT method, 25 key features were finally identified, including bands reflectance (e.g., blue, red-edge, and near-infrared), spectral indices (e.g., SI2-reg, EVI, and DVI, etc.), and texture data (e.g., the mean, variance, and second-order moments of the GLCM features, etc.). These features efficiently captured the spectral properties of the soil and the spatial structural information, providing the model with sufficient input variables. Numerous studies have demonstrated that the red edge band is more sensitive to soil salinity and is strongly correlated the spectral properties of the vegetation canopy, thereby improving the accuracy of soil salinity estimation\u003csup\u003e\u003cspan citationid=\"CR30\" class=\"CitationRef\"\u003e30\u003c/span\u003e,\u003cspan citationid=\"CR52\" class=\"CitationRef\"\u003e52\u003c/span\u003e\u003c/sup\u003e. To facilitate modeling and prediction, Hu et al.\u003csup\u003e\u003cspan citationid=\"CR53\" class=\"CitationRef\"\u003e53\u003c/span\u003e\u003c/sup\u003e employed UAV to gather remote sensing images of areas with both vegetation cover and barren ground. The regions with vegetation cover showed the most accurate prediction results. As a result, the spectral index produced by incorporating the red edge band significantly increases the accuracy of the soil salinity estimation in saline agricultural lands with vegetation cover. The red edge and near-infrared bands showed high importance at all depths in this study, as they sensitively revealed the changes in soil moisture and salt content, serving as important indicators for soil salinity prediction. This result aligned with the reported of Taghadosi et al.\u003csup\u003e\u003cspan citationid=\"CR54\" class=\"CitationRef\"\u003e54\u003c/span\u003e\u003c/sup\u003e. Spectral indices such as EVI, EVI-reg and BI indirectly indicated the salinity status of the soil by enhancing the vegetation characteristics. Lobell et al.\u003csup\u003e\u003cspan citationid=\"CR55\" class=\"CitationRef\"\u003e55\u003c/span\u003e\u003c/sup\u003e showed that EVI is more reliable than NDVI for salinity monitoring. Spectral index groups, such as SI2-reg, EVI-reg, and DVI-reg, involved the calculation of the red-edge band, further proving the importance of the red-edge band in soil salinity monitoring, consistent with the results of Ma et al.\u003csup\u003e\u003cspan citationid=\"CR56\" class=\"CitationRef\"\u003e56\u003c/span\u003e\u003c/sup\u003e Textural features, such as the mean, variance, and contrast in the GLCM, were instrumental in capturing subtle soil surface changes and spatial structure information. Tai\u003csup\u003e\u003cspan citationid=\"CR57\" class=\"CitationRef\"\u003e57\u003c/span\u003e\u003c/sup\u003e showed that incorporating texture features could effectively improve the accuracy of soil salinity estimation models under vegetation cover conditions.\u003c/p\u003e \u003cp\u003eDifferent modelling groups also exerted an influence on estimation model performance. In the single-variable modelling group, spectral indices had a greater contribution to soil salinity estimation, with the RF model demonstrating high accuracy across all four depths (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003e). Bian et al.\u003csup\u003e\u003cspan citationid=\"CR58\" class=\"CitationRef\"\u003e58\u003c/span\u003e\u003c/sup\u003e found that there is a certain correlation between spectral index and soil salinity, which has a great contribution to the application of estimation. Although the accuracy of the model built by band reflectance was slightly lower than that of spectral indices, the overall difference was not significant, which was in agreement with Wang et al.\u003csup\u003e\u003cspan citationid=\"CR59\" class=\"CitationRef\"\u003e59\u003c/span\u003e\u003c/sup\u003e, who evaluated the prediction accuracy of different variable groups for soil salinity in oasis. In contrast, texture data performed poorly among the single variables. However, when combined with other variables, the texture information from multispectral images could improve the accuracy of soil salinity estimation at different depths of vegetation cover. Zheng et al.\u003csup\u003e\u003cspan citationid=\"CR18\" class=\"CitationRef\"\u003e18\u003c/span\u003e\u003c/sup\u003e showed that the combination of texture data and spectral information significantly improved the accuracy of rice biomass estimation. The texture features of the UAV multispectral images provided rich information\u003csup\u003e\u003cspan citationid=\"CR60\" class=\"CitationRef\"\u003e60\u003c/span\u003e\u003c/sup\u003e, making them applicable for estimating soil salinity at different depths. Incorporating texture data significantly enhances model performance, particularly in the RF model. When texture data is combined with other data, the accuracy of the model is notably improved. This study also demonstrates that using sensitive bands and spectral indices as the input variables leads to better modeling and validation outcomes in soil salinity estimation models. Nevertheless, in multi-variate groups, variables may affect and constrain each other, and combining variables of high importance not always achieve optimal results. For instance, the accuracy of all groups in the BPNN model was not optimal compared to other groups. Introducing an excessive number of independent variables could lead to information redundancy, overfitting, and a decrease in model accuracy\u003csup\u003e\u003cspan citationid=\"CR61\" class=\"CitationRef\"\u003e61\u003c/span\u003e\u003c/sup\u003e.\u003c/p\u003e \u003cp\u003ePrevious studies have highlighted the superiority of machine learning methods in estimating soil salinity content due to the complexity and indirect relationships between variables\u003csup\u003e\u003cspan citationid=\"CR62\" class=\"CitationRef\"\u003e62\u003c/span\u003e\u003c/sup\u003e. In this study, RF, SVM, GPR and BPNN were used to model soil salinity at different depths and to select the most effective model. The evaluation of model criteria revealed that RF and GPR have advantages in estimating soil salinity. Although SVM has strong generalization capability for handling nonlinear problems, it is sensitive to noisy data\u003csup\u003e\u003cspan citationid=\"CR63\" class=\"CitationRef\"\u003e63\u003c/span\u003e\u003c/sup\u003e. The models exhibited varying levels of accuracy in soil salinity estimation at different depths. The GPR model demonstrates better prediction accuracy than other methods for surface and deep soils, while the RF model performed well in intermediate soils. Due to the complexity of soil salinization mechanism and the intricate nonlinear relationship between soil spectral characteristics and soil salinity, machine learning is well-suited for revealing these relationships. Their strong nonlinear fitting and generalization capabilities, make them ideal for simulating the complex interactions among variables. Ma et al.\u003csup\u003e\u003cspan citationid=\"CR56\" class=\"CitationRef\"\u003e56\u003c/span\u003e\u003c/sup\u003e found that the RF model had better estimation results when predicting the soil salinity in Ebinur Lake wetland using multispectral and Digital Elevation Model (DEM) data. Therefore, machine learning algorithms have shown superior performance in both modeling -and validation of soil salinity estimation model. For instance, Wei et al.\u003csup\u003e\u003cspan citationid=\"CR5\" class=\"CitationRef\"\u003e5\u003c/span\u003e\u003c/sup\u003e developed SSC estimation models using RF, SVM, and BPNN algorithms, with the RF model achieving the highest prediction accuracy (R\u003csup\u003e2\u003c/sup\u003e\u0026thinsp;=\u0026thinsp;0.84). Additionally, Zhu et al.\u003csup\u003e\u003cspan citationid=\"CR64\" class=\"CitationRef\"\u003e64\u003c/span\u003e\u003c/sup\u003e and Yu et al.\u003csup\u003e\u003cspan citationid=\"CR65\" class=\"CitationRef\"\u003e65\u003c/span\u003e\u003c/sup\u003e concluded that the RF based SSC estimation model had a high accuracy.\u003c/p\u003e \u003cp\u003eThe essence of estimating soil salinity during the vegetation cover period is to indirectly obtain information on soil salinity through the crop spectral response to soil salinity. The vegetation cover ratio of the given crop in the same period may correlate with the degree of soil salinity. The optimal estimation model for various soil depth varied with vegetation cover conditions, with the optimal model on depths ranging from 0ཞ10 cm under both low and medium vegetation cover conditions, and from 20ཞ30 cm under high vegetation cover. During the vegetation cover period, the absorption of soil water by crops was mainly done by the lateral root. Consequently, the soil salinity in the soil layer where the lateral roots were located will significantly affect crop water absorption\u003csup\u003e\u003cspan citationid=\"CR66\" class=\"CitationRef\"\u003e66\u003c/span\u003e\u003c/sup\u003e, causing crop growth stress. Higher soil salinity levels result in more severe stress and poorer vegetation growth, which would be indirectly expressed in the vegetation canopy.\u003c/p\u003e"},{"header":"5. Conclusion","content":"\u003cp\u003eIn this study, SSC at different depths was estimated and analyzed based on UAV multi-spectral data by composing seven type groups from band reflectance, spectral index and texture feature. Moreover, feature variables selection based on BDT method and four machine learning algorithms (RF, SVM, GPR and BPNN) were incorporated in estimation model establishment. (1) The estimation accuracy of the RF and GPR models surpasses that of the SVM and BPNN models across all four soil depths in the SSC estimation models. (2) The GPR model, based on the entire variable group, is the most accurate for estimating soil salt content at depths of 0\u0026thinsp;~\u0026thinsp;10 cm and 30\u0026thinsp;~\u0026thinsp;40 cm. The R\u003csup\u003e\u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2\u003c/span\u003e\u003c/sup\u003e values for the validation set were 0.774 and 0.705, and the modeling set, 0.77 and 0.62, with corresponding RMSE of 0.185% and 0.31%, respectively. The RF model, based on the spectral index and texture data set, is the most accurate estimation model for depths of 10 to 20 cm and 20 to 30 cm. (3) The model was further validated using salt data under varying vegetation covers. The highest validation accuracy was obtained with medium vegetation coverage at depths of 0\u0026thinsp;~\u0026thinsp;10 cm and 20\u0026thinsp;~\u0026thinsp;30 cm, as well as high vegetation coverage at depths of 10\u0026thinsp;~\u0026thinsp;20 cm and 30\u0026thinsp;~\u0026thinsp;40 cm. (4) The model imputes were screened using the BDT approach to select the independent variables at each depth. The selected results demonstrated that the model constructed using only texturing data group had low accuracy, but the accuracy increased dramatically when additional variables were included. There are still some problems and limitations in this study. This experiment was conducted under box planting conditions, and the experimental results have some limitations, which may not fully reflect the soil salinity conditions under actual field conditions. In addition, the study has not yet considered the effects of crop type and soil water content on the estimation model, which may have a significant impact on model accuracy in practical applications. Therefore, future studies need to consider more influencing factors in order to further optimize the estimation model and improve the accuracy and stability of soil salinity prediction.\u003c/p\u003e"},{"header":"Declarations","content":"\u003cp\u003e\u003cstrong\u003eCompeting Interests Statement\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe authors declare no competing interests.\u003c/p\u003e\u003ch2\u003eFunding\u003c/h2\u003e \u003cp\u003eThis research was funded by the National Natural Science Foundation of China (52379049, 52209071), the Postgraduate Research \u0026amp; Practice Innovation Program of Jiangsu Province (Yangzhou University) (No. SJCX23_1945), the China Postdoctoral Science Foundation (2023T160552, 2020M671623), \u0026ldquo;Chunhui Plan\u0026rdquo; Cooperative Scientific Research Project of Ministry of Education of China (HZKY20220115), the \u0026ldquo;Blue Project\u0026rdquo; of Yangzhou University and the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD).\u003c/p\u003e\u003ch2\u003eAuthor Contribution\u003c/h2\u003e\u003cp\u003eConceptualization, C.Z.; methodology, Z.L.; software, Z.L.; validation, M.D. and M.T.; formal analysis, Z.L.; resources, M.D.; data curation, M.D. and R.L.; Writing\u0026mdash;original preparation, Z.L.; Writing\u0026mdash;review and editing, C.Z. and M.T.; visualization, Z.Z.; supervision, C.Z. and S.F.; project administration, C.Z.; funding acquisition, C.Z. and S.F. All authors have read and agreed to the published version of the manuscript.\u003c/p\u003e\u003ch2\u003eAcknowledgement\u003c/h2\u003e\u003cp\u003eThe authors would like to acknowledge the Agricultural Water and Hydrological Ecology Experiment Station, Yangzhou University, for providing the experimental facilities.\u003c/p\u003e\u003ch2\u003eData Availability\u003c/h2\u003e\u003cp\u003eThe datasets used and analyzed during the current study are available from the corresponding author upon reasonable request.\u003c/p\u003e"},{"header":"References","content":"\u003col\u003e\u003cli\u003e\u003cspan\u003eSeifi, M., Ahmadi, A., Neyshabouri, M. R., Taghizadeh-Mehrjardi, R. \u0026amp; Bahrami, H. A. \u003cem\u003eRemote and Vis-NIR spectra sensing potential for soil salinization estimation in the eastern coast of Urmia hyper saline lake, Iran\u003c/em\u003e20100398 (Society and Environment, 2020).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWang, L. W. \u0026amp; Wei, Y. X. Estimating the total nitrogen and total phosphorus content of wetland soils using hyperspectral models. \u003cem\u003eActa Ecol. Sin\u003c/em\u003e. \u003cb\u003e36\u003c/b\u003e (16), 5116\u0026ndash;5125 (2016).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eRao, B. R. M. et al. Spectral behaviour of salt-affected soils. \u003cem\u003eInt. J. Remote Sens.\u003c/em\u003e \u003cb\u003e16\u003c/b\u003e (12), 2125\u0026ndash;2136 (1995).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eZhao, W., Zhou, C., Zhou, C., Ma, H. \u0026amp; Wang, Z. Soil salinity inversion model of oasis in arid area based on UAV multispectral remote sensing. \u003cem\u003eRemote Sens.\u003c/em\u003e \u003cb\u003e14\u003c/b\u003e (8), 1804 (2022).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWei, G. et al. Estimation of soil salt content by combining UAV-borne multispectral sensor and machine learning algorithms. \u003cem\u003ePeerJ\u003c/em\u003e, \u003cb\u003e8\u003c/b\u003e, e9087. (2020).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eChen, J. et al. UAV remote sensing inversion of soil salinity in field of sunflower. \u003cem\u003eTrans. Chin. Soc. Agricultural Mach.\u003c/em\u003e \u003cb\u003e51\u003c/b\u003e (7), 178\u0026ndash;191 (2020).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eIvushkin, K. et al. UAV based soil salinity assessment of cropland. \u003cem\u003eGeoderma\u003c/em\u003e. \u003cb\u003e338\u003c/b\u003e, 502\u0026ndash;512 (2019).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eNicolas, H. \u0026amp; Walter, C. Detecting salinity hazards within a semiarid context by means of combining soil and remote-sensing data. \u003cem\u003eGeoderma\u003c/em\u003e. \u003cb\u003e134\u003c/b\u003e (1\u0026ndash;2), 217\u0026ndash;230 (2006).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eQi, G., Zhao, G. \u0026amp; Xi, X. Soil salinity inversion of winter wheat areas based on satellite-unmanned aerial vehicle-ground collaborative system in coastal of the Yellow River Delta. \u003cem\u003eSensors\u003c/em\u003e. \u003cb\u003e20\u003c/b\u003e (22), 6521 (2020).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWang, N. et al. Integrating remote sensing and landscape characteristics to estimate soil salinity using machine learning methods: A case study from Southern Xinjiang, China. \u003cem\u003eRemote Sens.\u003c/em\u003e \u003cb\u003e12\u003c/b\u003e (24), 4118 (2020).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eFarifteh, J., Van der Meer, F., Atzberger, C. \u0026amp; Carranza, E. J. M. Quantitative analysis of salt-affected soil reflectance spectra: A comparison of two adaptive methods (PLSR and ANN). \u003cem\u003eRemote Sens. Environ.\u003c/em\u003e \u003cb\u003e110\u003c/b\u003e (1), 59\u0026ndash;78 (2007).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eXu, L. \u0026amp; Wang, Q. Retrieval of soil water content in saline soils from emitted thermal infrared spectra using partial linear squares regression. \u003cem\u003eRemote Sens.\u003c/em\u003e \u003cb\u003e7\u003c/b\u003e (11), 14646\u0026ndash;14662 (2015).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eGe, X. et al. Combining UAV-based hyperspectral imagery and machine learning algorithms for soil moisture content monitoring. \u003cem\u003ePeerJ\u003c/em\u003e, \u003cb\u003e7\u003c/b\u003e, e6926. (2019).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003ePoblete, T., Ortega-Far\u0026iacute;as, S., Moreno, M. A. \u0026amp; Bardeen, M. Artificial neural network to predict vine water status spatial variability using multispectral information obtained from an unmanned aerial vehicle (UAV). \u003cem\u003eSensors\u003c/em\u003e. \u003cb\u003e17\u003c/b\u003e (11), 2488 (2017).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHu, J. Estimation of Soil Salinity in Arid Area Based on Multi-Source Remote Sensing (Doctoral dissertation, Zhejiang University). (2019).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eZhang Fei, Z. F. et al. I., Extracting of soil salinization by SVM and accuracy evaluation based on texture characteristic. (2009).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eZhang, Z. et al. Inversion of soil salt content by UAV remote sensing under different vegetation coverage. \u003cem\u003eTrans. Chin. Soc. Agric. Mach.\u003c/em\u003e \u003cb\u003e53\u003c/b\u003e, 220\u0026ndash;230 (2022).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eZheng, H. et al. Enhancing the nitrogen signals of rice canopies across critical growth stages through the integration of textural and spectral information from unmanned aerial vehicle (UAV) multispectral imagery. \u003cem\u003eRemote Sens.\u003c/em\u003e \u003cb\u003e12\u003c/b\u003e (6), 957 (2020).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eYang, H., Wang, Z., Cao, J., Wu, Q. \u0026amp; Zhang, B. Estimating soil salinity using Gaofen-2 imagery: A novel application of combined spectral and textural features. \u003cem\u003eEnviron. Res.\u003c/em\u003e \u003cb\u003e217\u003c/b\u003e, 114870 (2023).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eRen, J., Li, X., Zhao, K., Fu, B. \u0026amp; Jiang, T. Study of an on-line measurement method for the salt parameters of soda-saline soils based on the texture features of cracks. \u003cem\u003eGeoderma\u003c/em\u003e. \u003cb\u003e263\u003c/b\u003e, 60\u0026ndash;69 (2016).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eAndrade, R. et al. Proximal sensor data fusion and auxiliary information for tropical soil property prediction: Soil texture. \u003cem\u003eGeoderma\u003c/em\u003e. \u003cb\u003e422\u003c/b\u003e, 115936 (2022).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eDrielsma, M. J., Love, J., Taylor, S., Thapa, R. \u0026amp; Williams, K. J. General Landscape Connectivity Model (GLCM): a new way to map whole of landscape biodiversity functional connectivity for operational planning and reporting. \u003cem\u003eEcol. Model.\u003c/em\u003e \u003cb\u003e465\u003c/b\u003e, 109858 (2022).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eCai, S., Zhang, R., Liu, L. \u0026amp; Zhou, D. A method of salt-affected soil information extraction based on a support vector machine with texture features. \u003cem\u003eMath. Comput. Model.\u003c/em\u003e \u003cb\u003e51\u003c/b\u003e (11\u0026ndash;12), 1319\u0026ndash;1325 (2010).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHuang, C. et al. Land cover mapping in cloud-prone tropical areas using Sentinel-2 data: Integrating spectral features with Ndvi temporal dynamics. \u003cem\u003eRemote Sens.\u003c/em\u003e \u003cb\u003e12\u003c/b\u003e (7), 1163 (2020).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eJia, K. L. \u0026amp; Zhang, J. H. Impacts of different alkaline soil on canopy spectral characteristics of overlying vegetation. \u003cem\u003eSpectrosc. Spectr. Anal.\u003c/em\u003e \u003cb\u003e34\u003c/b\u003e (3), 782\u0026ndash;786 (2014).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eZhang, J. et al. Estimating soil salinity with different fractional vegetation cover using remote sensing. \u003cem\u003eLand. Degrad. Dev.\u003c/em\u003e \u003cb\u003e32\u003c/b\u003e (2), 597\u0026ndash;612 (2021).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWang, F. et al. Sensitivity analysis of soil salinity and vegetation indices to detect soil salinity variation by using Landsat series images: applications in different oases in Xinjiang, China. \u003cem\u003eActa Ecol. Sin.\u003c/em\u003e \u003cb\u003e37\u003c/b\u003e (15), 5007\u0026ndash;5022 (2017).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eChi, Y., Sun, J., Liu, W., Wang, J. \u0026amp; Zhao, M. Mapping coastal wetland soil salinity in different seasons using an improved comprehensive land surface factor system. \u003cem\u003eEcol. Ind.\u003c/em\u003e \u003cb\u003e107\u003c/b\u003e, 105517 (2019).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eKhan, N. M., Rastoskuev, V. V., Shalina, E. V. \u0026amp; Sato, Y. Mapping salt-affected soils using remote sensing indicators-a simple approach with the use of GIS IDRISI. In 22nd Asian conference on remote sensing (Vol. 5, No. 9). (2001), November.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eZhang, S. et al. Integrated satellite, unmanned aerial vehicle (UAV) and ground inversion of the SPAD of winter wheat in the reviving stage. \u003cem\u003eSensors\u003c/em\u003e. \u003cb\u003e19\u003c/b\u003e (7), 1485 (2019).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSong, C. Y., Ren, H. X. \u0026amp; Huang, C. Estimating soil salinity in the Yellow River Delta, Eastern China\u0026mdash;an integrated approach using spectral and terrain indices with the generalized additive model. \u003cem\u003ePedosphere\u003c/em\u003e. \u003cb\u003e26\u003c/b\u003e (5), 626\u0026ndash;635 (2016).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eAbbas, A., Khan, S., Hussain, N., Hanjra, M. A. \u0026amp; Akbar, S. \u003cem\u003eCharacterizing soil salinity in irrigated agriculture using a remote sensing approach\u003c/em\u003e5543\u0026ndash;52 (Physics and chemistry of the Earth, 2013).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eTripathi, N. K., Rai, B. K. \u0026amp; Dwivedi, P. Spatial modeling of soil alkalinity in GIS environment using IRS data. In Proceedings of the 18th Asian Conference on Remote Sensing, Kuala Lumpur, Malaysia (Vol. 20, No. 25, pp. 81\u0026ndash;86). (1997), October.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eAllbed, A., Kumar, L. \u0026amp; Aldakheel, Y. Y. Assessing soil salinity using soil salinity and vegetation indices derived from IKONOS high-spatial resolution imageries: Applications in a date palm dominated region. \u003cem\u003eGeoderma\u003c/em\u003e. \u003cb\u003e230\u003c/b\u003e, 1\u0026ndash;8 (2014).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eZhou, L. et al. Research on quick dynamic monitoring of soil salinization in Qaidam Basin[J]. Science of Surveying and Mapping, 46(7): 99\u0026thinsp;~\u0026thinsp;106, 114. (2021).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eKhan, N. M., Rastoskuev, V. V., Sato, Y. \u0026amp; Shiozawa, S. Assessment of hydrosaline land degradation by using a simple approach of remote sensing indicators. \u003cem\u003eAgric. Water Manage.\u003c/em\u003e \u003cb\u003e77\u003c/b\u003e (1\u0026ndash;3), 96\u0026ndash;109 (2005).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLiu, H. Q. \u0026amp; Huete, A. A feedback based modification of the NDVI to minimize canopy background and atmospheric noise. \u003cem\u003eIEEE Trans. Geosci. Remote Sens.\u003c/em\u003e \u003cb\u003e33\u003c/b\u003e (2), 457\u0026ndash;465 (1995).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLyon, J. G., Yuan, D., Lunetta, R. S. \u0026amp; Elvidge, C. D. A change detection experiment using vegetation indices. \u003cem\u003ePhotogram. Eng. Remote Sens.\u003c/em\u003e \u003cb\u003e64\u003c/b\u003e (2), 143\u0026ndash;150 (1998).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eBirth, G. S. \u0026amp; McVey, G. R. Measuring the color of growing turf with a reflectance spectrophotometer 1. \u003cem\u003eAgron. J.\u003c/em\u003e \u003cb\u003e60\u003c/b\u003e (6), 640\u0026ndash;643 (1968).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eArachchi, M. H., Field, D. J. \u0026amp; McBratney, A. B. Quantification of soil carbon from bulk soil samples to predict the aggregate-carbon fractions within using near-and mid-infrared spectroscopic techniques. \u003cem\u003eGeoderma\u003c/em\u003e. \u003cb\u003e267\u003c/b\u003e, 207\u0026ndash;214 (2016).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMa, S., He, B., Ge, X. \u0026amp; Luo, X. Spatial prediction of soil salinity based on the Google Earth Engine platform with multitemporal synthetic remote sensing images. \u003cem\u003eEcol. Inf.\u003c/em\u003e \u003cb\u003e75\u003c/b\u003e, 102111 (2023).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eGitelson, A. A., Kaufman, Y. J., Stark, R. \u0026amp; Rundquist, D. Novel algorithms for remote estimation of vegetation fraction. \u003cem\u003eRemote Sens. Environ.\u003c/em\u003e \u003cb\u003e80\u003c/b\u003e (1), 76\u0026ndash;87 (2002).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eElsherbiny, O., Zhou, L., Feng, L. \u0026amp; Qiu, Z. Integration of visible and thermal imagery with an artificial neural network approach for robust forecasting of canopy water content in rice. \u003cem\u003eRemote Sens.\u003c/em\u003e \u003cb\u003e13\u003c/b\u003e (9), 1785 (2021).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eXiao, Y., Dong, Y., Huang, W., Liu, L. \u0026amp; Ma, H. Wheat fusarium head blight detection using UAV-based spectral and texture features in optimal window size. \u003cem\u003eRemote Sens.\u003c/em\u003e \u003cb\u003e13\u003c/b\u003e (13), 2437 (2021).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eDrucker, H. \u0026amp; Cortes, C. Boosting decision trees. Advances in neural information processing systems, 8. (1995).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMutanga, O., Adam, E. \u0026amp; Cho, M. A. High density biomass estimation for wetland vegetation using WorldView-2 imagery and random forest regression algorithm. \u003cem\u003eInt. J. Appl. Earth Obs. Geoinf.\u003c/em\u003e \u003cb\u003e18\u003c/b\u003e, 399\u0026ndash;406 (2012).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHuang, X. Y., Wang, X. M. \u0026amp; KAWUQIATI, B. Inversion of soil salinity of an oasis in an arid area based on Landsat8 OLI images. \u003cem\u003eRemote Sens. Nat. Resour.\u003c/em\u003e, \u003cb\u003e35\u003c/b\u003e(1). (2023).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSui, H. et al. Soil salinity estimation over coastal wetlands based on random Forest algorithm and hydrological connectivity metric. \u003cem\u003eFront. Mar. Sci.\u003c/em\u003e \u003cb\u003e9\u003c/b\u003e, 895172 (2022).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eTan, K., Ye, Y. Y., Du, P. J. \u0026amp; Zhang, Q. Q. Estimation of heavy metal concentrations in reclaimed mining soils using reflectance spectroscopy. \u003cem\u003eSpectrosc. Spectr. Anal.\u003c/em\u003e \u003cb\u003e34\u003c/b\u003e (12), 3317\u0026ndash;3322 (2014).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eGuan, X., Wang, S., Gao, Z. \u0026amp; Lv, Y. Dynamic prediction of soil salinization in an irrigation district based on the support vector machine. \u003cem\u003eMath. Comput. Model.\u003c/em\u003e \u003cb\u003e58\u003c/b\u003e (3\u0026ndash;4), 719\u0026ndash;724 (2013).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWang, J., Wang, W., Hu, Y., Tian, S. \u0026amp; Liu, D. Soil moisture and salinity inversion based on new remote sensing index and neural network at a salina-alkaline wetland. \u003cem\u003eWater\u003c/em\u003e. \u003cb\u003e13\u003c/b\u003e (19), 2762 (2021).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eKaplan, G. \u0026amp; Avdan, U. Evaluating the utilization of the red edge and radar bands from sentinel sensors for wetland classification. \u003cem\u003eCatena\u003c/em\u003e. \u003cb\u003e178\u003c/b\u003e, 109\u0026ndash;119 (2019).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHu, J. et al. Quantitative estimation of soil salinity using UAV-borne hyperspectral and satellite multispectral images. \u003cem\u003eRemote Sens.\u003c/em\u003e \u003cb\u003e11\u003c/b\u003e (7), 736 (2019).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eTaghadosi, M. M., Hasanlou, M. \u0026amp; Eftekhari, K. Retrieval of soil salinity from Sentinel-2 multispectral imagery. \u003cem\u003eEur. J. Remote Sens.\u003c/em\u003e \u003cb\u003e52\u003c/b\u003e (1), 138\u0026ndash;154 (2019).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLobell, D. B. et al. Regional-scale assessment of soil salinity in the Red River Valley using multi‐year MODIS EVI and NDVI. \u003cem\u003eJ. Environ. Qual.\u003c/em\u003e \u003cb\u003e39\u003c/b\u003e (1), 35\u0026ndash;41 (2010).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMa, G., Ding, J., Han, L. \u0026amp; Zhang, Z. Digital mapping of soil salinization in arid area wetland based on variable optimized selection and machine learning. \u003cem\u003eTrans. Chin. Soc. Agric. Eng.\u003c/em\u003e \u003cb\u003e36\u003c/b\u003e, 124\u0026ndash;131 (2020).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eTai, X. Soil Salinity Monitoring Model on UAV Multi Spectral Remote Sensing Vegetation Cover[D]. Northwest Agricultural and Forestry University, DOI: (2022). \u003cdiv class=\"ExternalRefDOI\"\u003e10.27409/d\u003c/div\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eBian, L., Wang, J., Guo, B., Cheng, K. \u0026amp; Wei, H. Remote sensing extraction of soil salinity in Yellow River Delta Kenli County based on feature space. \u003cem\u003eRemote Sens. Technol. Application\u003c/em\u003e. \u003cb\u003e35\u003c/b\u003e (1), 211\u0026ndash;218 (2020).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWang, F. et al. Environmental sensitive variable optimization and machine learning algorithm using in soil salt prediction at oasis. \u003cem\u003eTrans. Chin. Soc. Agric. Eng.\u003c/em\u003e \u003cb\u003e34\u003c/b\u003e, 102\u0026ndash;110 (2018).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eGuo, Y. et al. Integrating spectral and textural information for identifying the tasseling date of summer maize using UAV based RGB images. \u003cem\u003eInt. J. Appl. Earth Obs. Geoinf.\u003c/em\u003e \u003cb\u003e102\u003c/b\u003e, 102435 (2021).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMorellos, A. et al. Machine learning based prediction of soil total nitrogen, organic carbon and moisture content by using VIS-NIR spectroscopy. \u003cem\u003eBiosyst. Eng.\u003c/em\u003e \u003cb\u003e152\u003c/b\u003e, 104\u0026ndash;116 (2016).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eCui, X. et al. Estimating soil salinity under sunflower cover in the Hetao Irrigation District based on unmanned aerial vehicle remote sensing. \u003cem\u003eLand. Degrad. Dev.\u003c/em\u003e \u003cb\u003e34\u003c/b\u003e (1), 84\u0026ndash;97 (2023).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eShieh, H. L. \u0026amp; Kuo, C. C. A reduced data set method for support vector regression. \u003cem\u003eExpert Syst. Appl.\u003c/em\u003e \u003cb\u003e37\u003c/b\u003e (12), 7781\u0026ndash;7787 (2010).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eZhu, C., Ding, J., Zhang, Z. \u0026amp; Wang, Z. Exploring the potential of UAV hyperspectral image for estimating soil salinity: Effects of optimal band combination algorithm and random forest. \u003cem\u003eSpectrochim. Acta Part A Mol. Biomol. Spectrosc.\u003c/em\u003e \u003cb\u003e279\u003c/b\u003e, 121416 (2022).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eYu, X., Chang, C., Song, J., Zhuge, Y. \u0026amp; Wang, A. Precise monitoring of soil salinity in China\u0026rsquo;s Yellow River Delta using UAV-borne multispectral imagery and a soil salinity retrieval index. \u003cem\u003eSensors\u003c/em\u003e. \u003cb\u003e22\u003c/b\u003e (2), 546 (2022).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eZhang, Z. T. et al. UAV Multispectral Remote Sensing Soil Salinity Inversion Based on Different Fractional Vegetation Coverages[J]. \u003cem\u003eTrans. Chin. Soc. Agricultural Mach.\u003c/em\u003e, \u003cb\u003e53\u003c/b\u003e(08): 220\u0026ndash;230. (2022).\u003c/span\u003e\u003c/li\u003e\u003c/ol\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":true,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":false,"hideJournal":false,"highlight":"","institution":"","isAcceptedByJournal":true,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":false,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"scientific-reports","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"scirep","sideBox":"Learn more about [Scientific Reports](http://www.nature.com/srep/)","snPcode":"","submissionUrl":"","title":"Scientific Reports","twitterHandle":"","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"stoa","reportingPortfolio":"Scientific Reports","inReviewEnabled":true,"inReviewRevisionsEnabled":true},"keywords":"soil salt content, UAV, vegetation cover, soil depth, multi-source remote sensing data","lastPublishedDoi":"10.21203/rs.3.rs-4971758/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-4971758/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003cp\u003eSoil salinization is the most common land degradation problem in arid, semi-arid and coastal areas of China, which seriously affects local crop yield, economic development, and environmental sustainability. There are few studies on estimating soil salinity at different depths under vegetation cover. In this study, field soil control experiments were employed to collect multi-source remote sensing data under barley growth, and soil salt content (SSC) with various depths. Three types of feature variables were built based on images and were filtered by the boosting decision tree (BDT) method. Besides, four machine learning algorithms coupling with seven variable combination groups were used to comprehensively establish soil salinity estimation model. Finally, the performances of estimation model for different crop over ratios were evaluated. The results showed that the gaussian process regression (GPR) model based on the full variable group at the depths of 0\u0026thinsp;~\u0026thinsp;10 cm and 30\u0026thinsp;~\u0026thinsp;40 cm is more accurate than other models. The validation R\u003csup\u003e\u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2\u003c/span\u003e\u003c/sup\u003e is 0.774 and 0.705, and the RMSE is 0.185% and 0.31%;The random forest (RF) models based on spectral index and texture data at 10\u0026thinsp;~\u0026thinsp;20 cm and 20\u0026thinsp;~\u0026thinsp;30 cm depths are more accurate, with R\u003csup\u003e\u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2\u003c/span\u003e\u003c/sup\u003e of 0.666 and 0.714. SSC may be quantitatively inverted at various depths using the machine learning model based on multi-source remote sensing, which also serves as a guide for monitoring soil salinization.\u003c/p\u003e","manuscriptTitle":"Estimating soil profile salinity under vegetation cover based on UAV multi-source remote sensing","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2024-10-18 10:40:49","doi":"10.21203/rs.3.rs-4971758/v1","editorialEvents":[{"type":"communityComments","content":0},{"type":"decision","content":"Revision requested","date":"2024-10-30T05:12:16+00:00","index":"","fulltext":""},{"type":"reviewerAgreed","content":"96636350359983593172556678222506695281","date":"2024-10-19T14:54:23+00:00","index":"hide","fulltext":""},{"type":"editorInvitedReview","content":"","date":"2024-10-19T06:09:37+00:00","index":"hide","fulltext":""},{"type":"editorInvitedReview","content":"","date":"2024-10-09T06:59:55+00:00","index":"hide","fulltext":""},{"type":"reviewerAgreed","content":"158749014807589711203283651998850751268","date":"2024-10-02T15:59:23+00:00","index":"hide","fulltext":""},{"type":"reviewerAgreed","content":"189564883993777159164584403853638307588","date":"2024-09-29T00:35:34+00:00","index":"hide","fulltext":""},{"type":"reviewersInvited","content":"","date":"2024-09-26T18:24:58+00:00","index":"","fulltext":""},{"type":"editorAssigned","content":"","date":"2024-09-26T18:23:44+00:00","index":"","fulltext":""},{"type":"editorInvited","content":"","date":"2024-09-05T17:18:47+00:00","index":"","fulltext":""},{"type":"checksComplete","content":"","date":"2024-09-03T09:50:30+00:00","index":"","fulltext":""},{"type":"submitted","content":"Scientific Reports","date":"2024-08-25T08:02:47+00:00","index":"","fulltext":""}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"scientific-reports","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"scirep","sideBox":"Learn more about [Scientific Reports](http://www.nature.com/srep/)","snPcode":"","submissionUrl":"","title":"Scientific Reports","twitterHandle":"","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"stoa","reportingPortfolio":"Scientific Reports","inReviewEnabled":true,"inReviewRevisionsEnabled":true}}],"origin":"","ownerIdentity":"e2c052a0-d84c-4aa5-86d6-cac94a000814","owner":[],"postedDate":"October 18th, 2024","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"published-in-journal","subjectAreas":[{"id":38468838,"name":"Earth and environmental sciences/Biogeochemistry"},{"id":38468839,"name":"Earth and environmental sciences/Environmental sciences"},{"id":38468840,"name":"Physical sciences/Engineering"}],"tags":[],"updatedAt":"2025-01-27T16:01:18+00:00","versionOfRecord":{"articleIdentity":"rs-4971758","link":"https://doi.org/10.1038/s41598-024-82868-9","journal":{"identity":"scientific-reports","isVorOnly":false,"title":"Scientific Reports"},"publishedOn":"2025-01-21 15:57:29","publishedOnDateReadable":"January 21st, 2025"},"versionCreatedAt":"2024-10-18 10:40:49","video":"","vorDoi":"10.1038/s41598-024-82868-9","vorDoiUrl":"https://doi.org/10.1038/s41598-024-82868-9","workflowStages":[]},"version":"v1","identity":"rs-4971758","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-4971758","identity":"rs-4971758","version":["v1"]},"buildId":"qtupq5eGEP_6zYnWcrvyt","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}