Dynamic changes in soil organic matter content affects the sustainable supply of soil water and fertilizer and impacts the stability of soil ecological function. Understanding the spatial distribution characteristics of soil organic matter will help deepen our understanding of the differences in soil organic matter content, soil formation law; such understanding would be useful for rational land use planning. Taking terrain data, meteorological data, and remote sensing data as auxiliary variables and the ordinary Kriging (OK) method as a control, this study compares the spatial prediction accuracies and mapping effects of various models (MLR, RK, GWR, GWRK, MGWR, and MGWRK) on soil organic matter. Our results show that the spatial distribution trend of soil organic matter predicted by each model is similar, but the prediction of composite models can reflect more mapping details than that of unitary models. The OK method can provide better support for spatial prediction when the sampling points are dense; however, the local models are superior in dealing with spatial non-stationarity. Notably, the MGWR model is superior to the GWR model, but the MGWRK model is inferior to the GWRK model. As a new method, the prediction accuracy of MGWRK reached 47.72% for the OK and RK methods and 40.08% for the GWRK method. The GWRK method achieved a better prediction accuracy. The influence mechanism of soil organic matter is complex, but the MGWR model more clearly reveals the complex nonlinear relationship between soil organic matter content and factors influencing it. This research can provide reference methods and mapping technical support to improve the spatial prediction accuracy of soil organic matter.