How Do You Model Soil Carbon Change Over Time?

A soil carbon model takes a description of a system (climate, soil, management) and returns a trajectory: how the carbon stock is expected to evolve under that description. That trajectory is how organisations justify capital allocation to regenerative programmes, forecast Scope 3 emissions and removals under the GHG Protocol Land Sector and Removals Guidance, and, where the project structure warrants it, file claims on voluntary carbon markets.

The question “which model should we use?” has no universal answer. The right question is narrower: which model is fit for the decision, on the system at hand, given the data the project can realistically get? Answering it is the single most important judgment call in a soil carbon modelling engagement.

What modelling actually delivers, and what it doesn't

A soil carbon model is a structured hypothesis about how carbon moves through the system. Its outputs are conditional: given these inputs, here is the expected trajectory and its uncertainty. A model does not measure the soil; it predicts it.

This matters because much of the disagreement between project developers and verifiers comes from confusing the two1. A modelled trajectory substituted for field data in annual reporting is a compliance risk. A modelled trajectory used alongside field data, with uncertainty ranges explicit, is standard practice under the GHG Protocol Land Sector and Removals Guidance and under most voluntary market methodologies.

The candidate models

Four frameworks dominate practical work. Each was designed for a specific class of question; using one outside its domain is possible but expensive in credibility.

All four have been benchmarked against long-term experimental data across dozens of sites. RothC and Century have the longest international track record; AMG is the strongest option for French and other European cropping systems; DayCent is the specialist tool when N₂O emissions are in scope2.

Choosing the right model for the question

Four practical questions determine the choice.

Is nitrogen in scope?

If the decision depends on N₂O emissions, fertilizer-use efficiency, or the interaction between nitrogen management and carbon dynamics, RothC is structurally the wrong choice; it has no N pool. Century or DayCent are the honest answers. If the decision is purely about soil carbon under stable nitrogen management, RothC is defensible and simpler.

Is the system water-limited?

RothC's default water modifier underperforms in dryland and Mediterranean systems, where biological activity is limited by moisture for much of the year3. A modified RothC calibration (Farina et al.) exists and is the practical choice in those contexts. Century and DayCent handle water limitation natively through their submodels.

What time resolution is needed?

Monthly is sufficient for a 20-year decadal trajectory. Daily is required when the outcome of interest is a rapid flux: trace-gas emissions after a fertilizer application, for instance. RothC, AMG, and Century operate monthly; DayCent operates daily.

What data can the project realistically produce?

A model that needs inputs the project cannot supply is a model that will run on assumptions. If daily weather is not available, DayCent will operate on a downscaling, which is defensible if the downscaling is disclosed. If residue inputs and management history are unknown, every model is in the same position: the uncertainty from missing inputs dominates the uncertainty from the model.

Initialization: the silent source of error

Every soil carbon model partitions the stock into pools with different turnover rates: fast, slow, and near-inert. The initial split across pools is rarely measurable and is usually estimated. For RothC, the Falloon equation estimates the inert pool from total SOC; it is a statistical approximation with a material effect on 30-year projections4.

The equation above is a statistical fit, not a physical measurement. Its parameters carry uncertainty that compounds over decadal projections. The more defensible alternative is spin-up: run the model to equilibrium under a reconstructed historical management regime, then use the resulting pool distribution as the starting point. The reconstruction is the hard part, and is also the step most commonly glossed over in project documentation. A 50-year historical back-cast under a realistic rotation reduces the effect of initialization on forecast trajectories by an order of magnitude compared with default-pool initialization5.

Reading outputs as envelopes, not lines

A soil carbon trajectory with a tight 95 percent confidence interval either reflects a well-calibrated, data-rich system, or underreporting. Monte Carlo propagation of input uncertainty (climate variability, clay content, residue inputs, initialization) produces a trajectory envelope that is usually substantially wider than the mean line. That envelope is the honest deliverable. Strip it out and the remaining number is, at best, a plausible central estimate.

For scenario comparison work (baseline vs. intervention, for instance), the useful output is the difference between scenarios, with its own propagated uncertainty. Much of the input noise is shared between scenarios and cancels in the difference, which is often tighter than either individual trajectory. This is the right framing for nearly all carbon project feasibility questions.

Key takeaways