Representing soil organic carbon (SOC) dynamics in Earth system models (ESMs) is a key source of uncertainty in predicting carbon–climate feedbacks. Using machine learning, we develop and compare predictive relationships in observations (Obs) and ESMs. We find different relationships between environmental factors and SOC stocks in Obs and ESMs. SOC prediction in ESMs may be improved by representing the functional relationships of environmental controllers in a way consistent with observations.