Many studies have attempted to predict soil organic matter (SOM), whereas mapping high-precision and high-resolution SOM maps remains a challenge due to the difficulty of selecting appropriate satellite data sources and prediction algorithms. This study aimed to investigate the influence of different remotely sensed images and machine learning algorithms on SOM prediction. We constructed two comparative experiments, i.e., full-band and common-band variable datasets of Sentinel-2A and MODIS images using Google Earth Engine (GEE). The predictive performances of random forest (RF), artificial neural network (ANN), and support vector regression (SVR) algorithms were evaluated, and the SOM map was generated for the Songnen Plain. Results showed that the model based on the full-band Sentinel-2A dataset achieved the best performance. The application of Sentinel-2A data resulted in mean relative improvements (RIs) of 7.67% and 5.87%, respectively. The RF achieved a lower root mean squared error (RMSE = 0.68%) and a higher coefficient of determination (R² = 0.67) in all of the predicted scenarios than ANN and SVR. The resultant SOM map accurately characterized the SOM spatial distribution. Therefore, the Sentinel-2A data have obvious advantages over MODIS due to their higher spectral and spatial resolutions, and the combination of the RF algorithm and GEE is an effective approach to SOM mapping.

Soil organic matter (SOM) plays a critical role in agroecosystems and the terrestrial carbon cycle. Thus, accurately mapping SOM promotes sustainable agriculture and estimations of soil carbon pools. However, few studies have analyzed the changing trends in multi-period SOM prediction accuracies for single cropland soil types and mapped their spatial SOM patterns. Using time series 7 MOD09A1 images during the bare soil period, we combined the pixel dates of training samples and precipitation data to explore the variation in SOM accuracy for two typical cropland soil types. The advantage of using single soil type data versus the total dataset was evaluated, and SOM maps were drawn for the northern Songnen Plain. When almost no precipitation occurred on or near the optimal pixel date, the accuracies increased, and vice versa. SOM models of the two soil types achieved a lower root mean squared error (RMSE = 0.55%, 0.79%) and mean absolute error (MAE = 0.39%, 0.58%) and a higher coefficient of determination (R² = 0.65, 0.75) than the model using the total dataset and resulted in a mean relative improvement (RI) of 30.21%. The SOM decreased from northeast to southwest. The results provide reference data for the accurate management of cultivated soil and determining carbon sequestration.