Due to the importance of soil organic carbon (SOC) in supporting ecosystem services, accurate SOC assessment is vital for scientific research and decision making. However, most previous studies focused on single soil depth, leading to a poor understanding of SOC in multiple depths. To better understand the spatial distribution pattern of SOC in Northeast and North China Plain, we compared three machine learning algorithms (i.e., Cubist, Extreme Gradient Boosting (XGBoost) and Random Forest (RF)) within the digital soil mapping framework. A total of 386 sampling sites (1584 samples) following specific criteria covering all dryland districts and counties and soil types in four depths (i.e., 0–10, 10–20, 20–30 and 30–40 cm) were collected in 2017. After feature selection from 249 environmental covariates by the Genetic Algorithm, 29 variables were used to fit models. The results showed SOC increased from southern to northern regions in the spatial scale and decreased with soil depths. From the result of independent verification (validation dataset: 80 sampling sites), RF (R²: 0.58, 0.71, 0.73, 0.74 and RMSE: 3.49, 3.49, 2.95, 2.80 g kg⁻¹ in four depths) performed better than Cubist (R²: 0.46, 0.63, 0.67, 0.71 and RMSE: 3.83, 3.60, 3.03, 2.72 g kg⁻¹) and XGBoost (R²: 0.53, 0.67, 0.70, 0.71 and RMSE: 3.60, 3.60, 3.00, 2.83 g kg⁻¹) in terms of prediction accuracy and robustness. Soil, parent material and organism were the most important covariates in SOC prediction. This study provides the up-to-date spatial distribution of dryland SOC in Northeast and North China Plain, which is of great value for evaluating dynamics of soil quality after long-term cultivation.

Timely monitoring and early warning of soil salinity are crucial for saline soil management. Environmental variables are commonly used to build soil salinity prediction model. However, few researches have been done to summarize the environmental sensitive variables for soil electrical conductivity (EC) estimation systematically. Additionally, the performance of Multiple Linear Regression (MLR), Geographically Weighted Regression (GWR), and Random Forest regression (RFR) model, the representative of current main methods for soil EC prediction, has not been explored. Taking the north of Yinchuan plain irrigation oasis as the study area, the feasibility and potential of 64 environmental variables, extracted from the Landsat 8 remote sensed images in dry season and wet season, the digital elevation model, and other data, were assessed through the correlation analysis and the performance of MLR, GWR, and RFR model on soil salinity estimation was compared. The results showed that: 1) 10 of 15 imagery texture and spectral band reflectivity environmental variables extracted from Landsat 8 image in dry season were significantly correlated with soil EC, while only 3 of these indices extracted from Landsat 8 image in wet season have significant correlation with soil EC. Channel network base level, one of the terrain attributes, had the largest absolute correlation coefficient of 0.47 and all spatial location factors had significant correlation with soil EC. 2) Prediction accuracy of RFR model was slightly higher than that of the GWR model, while MLR model produced the largest error. 3) In general, the soil salinization level in the study area gradually increased from south to north. In conclusion, the remote sensed imagery scanned in dry season was more suitable for soil EC estimation, and topographic factors and spatial location also play a key role. This study can contribute to the research on model construction and variables selection for soil salinity estimation in arid and semiarid regions.

The effects of environmental factors on topsoil nutrient distribution have been extensively discussed, but it remains unclear how they affect spatial characteristics of soil carbon (C), nitrogen (N), and phosphorus (P) stoichiometry at different depths. We collected 184 soil samples in the typical black soil region of northeast China. Ordinary kriging was performed to describe the spatial distribution of soil C, N, and P eco-stoichiometry. Redundancy analysis was used to explore relationships between C:N:P ratios and physicochemical characteristics. The soil classification was studied by hierarchical cluster analysis. The mean C, N, and P contents ranged from 15.67 to 20.08 g·kg⁻¹, 1.15 to 1.51 g·kg⁻¹, and 0.80 to 0.90 g·kg⁻¹ within measured depths. C, N, and P concentrations and stoichiometry increased from southwest to northeast, and the Songhua River was identified as an important transition zone. At 0–20 cm, soil water content explained most of the C, N, and P content levels and ratios in cluster 1, while latitude had the highest explanatory ability in cluster 2. For 20–40 cm, soil bulk density was the main influencing factor in both clusters. Our findings contribute to an improved knowledge of the balance and ecological interactions of C, N, and P in northeast China for its sustainability.