Quantification of soil organic carbon (SOC) and pH, and their spatial variations at regional scales, is a foundation to adequately assess agriculture, pollution control, or environmental health and ecosystem functioning, so as to establish better practices for land use and land management. In this study, we used the random forest (RF) model to map the distribution of SOC and pH in the topsoil (0–20 cm) and estimate SOC and pH changes from 1982 to 2012 in Liaoning Province, Northeast China. A total of 10 covariates (elevation, slope gradient, topographic wetness index (TWI), mean annual temperature (MAT), mean annual precipitation (MAP), visible-red band 3 (B3), near-infrared band 4 (B4), short-wave infrared band 5 (B5), normalized difference vegetation index (NDVI), and land-use data) and a set of 806 (in 1982) and 973 (in 2012) soil samples were selected. Cross-validation technology was used to test the performance and uncertainty of the RF model. We found that the prediction R2 of SOC and pH was 0.69 and 0.54 for 1982, and 0.63 and 0.48 for 2012, respectively. Elevation, NDVI, and land use are the main environmental variables affecting the spatial variability of SOC in both periods. Correspondingly, the topographic wetness index and mean annual precipitation were the two most critical environmental variables affecting the spatial variation of pH. The mean SOC and pH decreased from 18.6 to 16.9 kg−1 and 6.9 to 6.6, respectively, over a 30-year period. SOC distribution generated using the RF model showed a decreasing SOC trend from east to west across the city in the two periods. In contrast, the spatial distribution of pH showed an opposite trend in both periods. This study provided important information of spatial variations in SOC and pH to agencies and communities in this region, to evaluate soil quality and make decisions on remediation and prevention of soil acidification and salinization.

Accurately mapping the spatial distribution information of soil organic carbon (SOC) stocks is a key premise for soil resource management and environment protection. Rapid development of satellite remote sensing provides a great opportunity for monitoring SOC stocks at a large scale. In this study, based on 12 environmental variables of multispectral remote sensing, topography and climate and 236 soil sampling data, three different boosted regression tree (BRT) models were compared to obtain the most accurate map of SOC stocks covering the forest area of Lvshun District in the Northeast China. Four validation indexes, including mean absolute error (MAE), root mean square error (RMSE), coefficient of determination (R2), and Lin’s concordance correlation coefficient (LCCC) were calculated to evaluate the performance of the three models. The results showed that the full variable model performed the best, except the model using multispectral remote sensing variables. In the full variable model, the regional SOC stocks are primarily determined by multispectral remote sensing variables, followed by topographic and climatic variables, with the relative importance of variables in the model being 63%, 28%, and 9%, respectively. The average prediction results of full variables model and only multispectral remote sensing variables model were 8.99 and 9.32 kg m−2, respectively. Our results indicated that there is a strong dependence of SOC stocks on multispectral remote sensing data when forest ecosystems have dense natural vegetation. Our study suggests that the multispectral remote sensing variables should be used to map SOC stocks of forest ecosystems in our study region.