Land surface models are important tools for understanding and predicting the land components of the carbon cycle. Atmospheric CO₂ concentration data are a valuable source of information that can be used to improve the accuracy of these models. In this study, we present a statistical ensemble-variational data assimilation method named EnVarDA to calibrate parameters of a land surface model using these data. We show that this method is easy to implement and more efficient and accurate than traditional methods.