Saturation effects limit the application of vegetation indices (VIs) in dense vegetation areas. The possibility to mitigate them by adopting a negative soil adjustment factor X is addressed. Two leaf area index (LAI) data sets are analyzed using the Google Earth Engine (GEE) for validation. The first one is derived from observations of MODerate resolution Imaging Spectroradiometer (MODIS) from 16 April 2013, to 21 October 2020, in the Apiacás area. Its corresponding VIs are calculated from a combination of Sentinel-2 and Landsat-8 surface reflectance products. The second one is a global LAI dataset with VIs calculated from Landsat-5 surface reflectance products. A linear regression model is applied to both datasets to evaluate four VIs that are commonly used to estimate LAI: normalized difference vegetation index (NDVI), soil adjusted vegetation index (SAVI), transformed SAVI (TSAVI), and enhanced vegetation index (EVI). The optimal soil adjustment factor of SAVI for LAI estimation is determined using an exhaustive search. The Dickey-Fuller test indicates that the time series of LAI data are stable with a confidence level of 99%. The linear regression results stress significant saturation effects in all VIs. Finally, the exhaustive searching results show that a negative soil adjustment factor of SAVI can mitigate the SAVIs’ saturation in the Apiacás area (i.e., X = −0.148 for mean LAI = 5.35), and more generally in areas with large LAI values (e.g., X = −0.183 for mean LAI = 6.72). Our study further confirms that the lower boundary of the soil adjustment factor can be negative and that using a negative soil adjustment factor improves the computation of time series of LAI.

Yan’an new district (YND) is one of the largest civil engineering projects for land creation in Loess Plateau, of which the amount of earthwork exceeds 600 million m³, to create 78.5 km² of flat land. Such mega-scale engineering activities and complex geological characteristics have induced wide land deformation in the region. Small baseline subset synthetic aperture radar interferometry (SBAS-InSAR) method and 55 Sentinel-1A (S-1A) images were utilized in the present work to investigate the urban surface deformation in the Yan’an urban area and Yan’an new airport (YNA) from 2015 to 2019. The results were validated by the ground leveling measurements in the YNA. It is found that significant uneven surface deformation existed in both YND and YNA areas with maximum accumulative subsidence of 300 and 217 mm, respectively. Moreover, the average subsidence rate of the YND and YNA areas ranged from −70 to 30 mm/year and −50 to 25 mm/year, respectively. The present work shows that the land deformation suffered two periods (from 2015 to 2017 and from 2017 to 2019) and expanded from urban center to surrounding resettlement area, which are highly relevant with urban earthwork process. It is found that more than 60% of land subsidence occurs at filled areas, while more than 65% of surface uplifting occurs at excavation areas. The present work shows that the subsidence originates from the earth filling and the load of urban buildings, while the release of stress is the major factor for the land uplift. Moreover, it is found that the collapsibility of loess and concentrated precipitation deteriorates the degree of local land subsidence. The deformation discovered by this paper shows that the city may suffer a long period of subsidence, and huge challenges may exist in the period of urban maintaining buildings and infrastructure facilities.

The suitability evaluation of agricultural land at the regional scale is of great significance for protecting land and water resources and building sustainable agricultural systems. Based on climate, soil, topographical, and surface water resources, land suitability index (LSI) data for maize, rice, and soybeans are established using an analytical hierarchy process and matter element analysis (AHP–MEA) model in Jilin Province, China. The results show that there is a significant positive linear correlation between the LSI and the measured yield, which indicates that the model has an ideal effect and certain reference and extension significance. The main limiting factors for maize and soybean planting are pH, total nitrogen (TN), available phosphorus (AP), and soil texture, while water shortage limits rice planting. Different spatial structure optimization schemes for planting are established using the LSI and measured yield, along with economic indices. This study shows that the scheme that integrates policy and cost can make full use of land and water resources and promote the economic growth of agriculture. After optimization, the planting areas of maize, rice, and soybeans were 7.22, 2.44, and 0.71 million ha, respectively, representing an increase of 15.71 billion yuan over the agricultural GDP for the existing planting structure. It is expected that this study will provide a basis for follow-up studies on crop cultivation suitability.

With the development of industrialization and urbanization, heavy metal contamination in agricultural soils tends to accumulate rapidly and harm human health. Visible and near-infrared (Vis-NIR) spectroscopy provides the feasibility of fast monitoring of the variation of heavy metals. This study explored the potential of fractional-order derivative (FOD), the optimal band combination algorithm and different mathematical models in estimating soil heavy metals with Vis-NIR spectroscopy. A total of 80 soil samples were collected from an agriculture area in Suzi river basin, Liaoning Province, China. The spectra for mercury (Hg), chromium (Cr), and copper (Cu) of the samples were obtained in the laboratory. For spectral preprocessing, FODs were allowed to vary from 0 to 2 with an increment of 0.2 at each step, and the optimal band combination algorithm was applied to the spectra after FOD. Then, four mathematical models, namely, partial least squares regression (PLSR), adaptive neural fuzzy inference system (ANFIS), random forest (RF) and generalized regression neural network (GRNN), were used to estimate the concentration of Hg, Cr and Cu. Results showed that high-order FOD had an excellent effect in highlighting hidden information and separating minor absorbing peaks, and the optimal band combination algorithm could remove the influence of spectral noise caused by high-order FOD. The incorporation of the optimal band combination algorithm and FOD is able to further mine spectral information. Furthermore, GRNN made an obvious improvement to the estimation accuracy of all studied heavy metals compared to ANFIS, PLSR, and RF. In summary, our results provided more feasibility for the rapid estimation of Hg, Cr, Cu and other heavy metal pollution areas in agricultural soils.

Black soil in northeast China is gradually degraded and soil organic matter (SOM) content decreases at a rate of 0.5% per year because of the long-term cultivation. SOM content can be obtained rapidly by visible and near-infrared (Vis–NIR) spectroscopy. It is critical to select appropriate preprocessing techniques for SOM content estimation through Vis–NIR spectroscopy. This study explored three categories of preprocessing techniques to improve the accuracy of SOM content estimation in black soil area, and a total of 496 ground samples were collected from the typical black soil area at 0–15 cm in Hai Lun City, Heilongjiang Province, northeast of China. Three categories of preprocessing include denoising, data transformation and dimensionality reduction. For denoising, Svitzky-Golay filter (SGF), wavelet packet transform (WPT), multiplicative scatter correction (MSC), and none (N) were applied to spectrum of ground samples. For data transformation, fractional derivatives were allowed to vary from 0 to 2 with an increment of 0.2 at each step. For dimensionality reduction, multidimensional scaling (MDS) and locally linear embedding (LLE) were introduced and compared with principal component analysis (PCA), which was commonly used for dimensionality reduction of soil spectrum. After spectral pretreatments, a total of 132 partial least squares regression (PLSR) models were constructed for SOM content estimation. Results showed that SGF performed better than the other three denoising methods. Low-order derivatives can accentuate spectral features of soil for SOM content estimation; as the order increases from 0.8, the spectrum were more susceptible to spectral noise interferences. In most cases, 0.2–0.8 order derivatives exhibited the best estimation performance. Furthermore, PCA yielded the optimal predictability, the mean residual predictive deviation (RPD) and maximum RPD of the models using PCA were 1.79 and 2.60, respectively. The application of appropriate preprocessing techniques could improve the efficiency and accuracy of SOM content estimation, which is important for the protection of ecological and agricultural environment in black soil area.