Representing soil organic carbon (SOC) dynamics in Earth system models (ESMs) is a key source of uncertainty in predicting carbon climate feedbacks. We used machine learning to develop and compare predictive relationships in observations and ESMs. We found different relationships between environmental factors and SOC stocks in observations and ESMs. SOC predictions in ESMs may be improved by representing the functional relationships of environmental controllers consistent with observations.