Comment modéliser l’évolution du carbone du sol dans le temps ?

Un modèle de carbone du sol prend une description d’un système (climat, sol, gestion) et renvoie une trajectoire : la façon dont le stock de carbone devrait évoluer selon cette description. Cette trajectoire est ce qui permet aux organisations de justifier l’allocation de capital à des programmes régénératifs, de prévoir les émissions et absorptions du Scope 3 au titre de la GHG Protocol Land Sector and Removals Guidance et, lorsque la structure du projet le justifie, de faire valoir des crédits sur les marchés volontaires du carbone.

La question « quel modèle devrions-nous utiliser ? » n’a pas de réponse universelle. La bonne question est plus circonscrite : quel modèle est adapté à la décision, sur le système considéré, compte tenu des données que le projet peut raisonnablement obtenir ? Y répondre constitue le jugement le plus déterminant de toute mission de modélisation du carbone du sol.

Sous les acronymes, les modèles usuels partagent un même mécanisme : le carbone du sol réparti en compartiments qui se décomposent chacun à leur propre rythme, alimentés par un flux d’apports frais de matière végétale et de fumier. Comprenez ce mécanisme et la plupart des questions pratiques se résolvent d’elles-mêmes, y compris les deux qui décident discrètement d’une prévision décennale, à savoir l’endroit où l’on suppose que démarrent les compartiments et la largeur de la réponse honnête. Les trois outils de cet article font tourner une version compacte de ce mécanisme : le comportement devient ainsi quelque chose que l’on peut manipuler, et pas seulement lire.

Ce que la modélisation apporte réellement, et ce qu’elle n’apporte pas

Un modèle de carbone du sol est une hypothèse structurée sur la façon dont le carbone se déplace dans le système. Ses sorties sont conditionnelles : étant donné ces données d’entrée, voici la trajectoire attendue et son incertitude. Un modèle ne mesure pas le sol ; il le prédit.

Cela compte, car une grande part du désaccord entre développeurs de projets et vérificateurs vient de la confusion entre les deux1. Une trajectoire modélisée substituée aux données de terrain dans le reporting annuel est un risque de conformité. Une trajectoire modélisée utilisée aux côtés de données de terrain, avec des plages d’incertitude explicites, est une pratique courante au titre de la GHG Protocol Land Sector and Removals Guidance et de la plupart des méthodologies des marchés volontaires.

Comment ces modèles fonctionnent réellement

Le carbone du sol n’est pas une substance unique qui se décompose à une vitesse unique. Une racine fraîche et un complexe humique vieux de plusieurs siècles sont tous deux du « carbone organique du sol », et pourtant l’un se renouvelle en quelques mois quand l’autre survit à celui qui le mesure. Tous les modèles d’usage courant traitent cela de la même façon : ils répartissent le carbone en une poignée de compartiments, chacun doté de son propre temps de renouvellement, et les suivent séparément2.

Au sein de chaque compartiment, la décomposition suit une cinétique du premier ordre : la quantité perdue au cours d’un mois donné est proportionnelle à la quantité présente, de sorte qu’un compartiment laissé à lui-même décroît selon une courbe exponentielle fixée par une constante de vitesse k.

La constante de vitesse k est une propriété du compartiment : élevée pour les résidus frais, faible pour la matière humifiée, quasi nulle pour la fraction inerte. La fonction f module k à la hausse ou à la baisse selon les conditions réelles du sol. La chaleur accélère la décomposition, le stress hydrique la ralentit, et un couvert vivant la freine par rapport à un sol nu. C’est pourquoi un même apport de résidus accumule du carbone sous une culture de couverture lors d’une saison tempérée fraîche, et en perd lors d’une saison chaude et sèche.

Tout le carbone décomposé ne quitte pas le sol. Une fraction est respirée vers l’atmosphère sous forme de CO₂ ; le reste se stabilise en biomasse microbienne et en composés humifiés qui réalimentent les compartiments lents. La teneur en argile gouverne ce partage, car les sols à texture plus fine protègent une plus grande part du carbone transformé. C’est pourquoi la fraction argileuse est une donnée d’entrée requise dans chacun des modèles ci-dessous3.

Assemblez ces éléments et un sol sous gestion stable dérive vers un état stationnaire, le stock auquel les apports annuels équilibrent exactement les pertes annuelles. Augmentez les apports et l’état stationnaire s’élève ; réchauffez le climat ou abaissez la teneur en argile et il baisse. L’approche est lente, cadencée par le temps de renouvellement du compartiment actif le plus lent, ce qui explique pourquoi les gains crédibles de carbone du sol se comptent en décennies plutôt qu’en années. L’outil ci-dessous fait avancer ce mécanisme. Faites varier l’apport, l’argile et la température, et observez où le stock se stabilise et le temps qu’il met pour y parvenir.

Deux traits de cette courbe importent pour un projet. La destination est fixée par les apports et le site, non par le point de départ du sol, de sorte que deux parcelles conduites sous la même gestion convergent vers le même état stationnaire par des chemins différents. Et le trajet est long : même un changement important met des décennies à se manifester, et les premières années avancent lentement. Une trajectoire modélisée qui atteint un nouvel équilibre en cinq ans décrit un artefact de son paramétrage, pas un sol.

Les quatre modèles candidats

Quatre cadres dominent la pratique. Chacun a été conçu pour une classe de questions précise, et chacun porte l’empreinte du système sur lequel il a été bâti. RothC est né des essais de longue durée de Rothamsted sur les horizons de surface arables tempérés2. AMG a d’abord été calibré sur des essais de grandes cultures français, puis reconstruit, dans sa forme actuelle, sur plus de quarante sites européens de longue durée45. Century a été conçu pour expliquer les niveaux de matière organique des prairies des Grandes Plaines, avec l’azote, le phosphore et le soufre couplés au cycle du carbone6, et DayCent en est l’homologue au pas de temps journalier, bâti pour résoudre les flux de gaz traces qu’un modèle mensuel ne peut saisir7. En utiliser un hors de son domaine d’origine est possible mais coûteux en crédibilité. Sur la manière de vérifier qu’un nouveau site relève bien de ce domaine, voir Représentatif de quoi ?

Les quatre ont été confrontés à des données expérimentales de longue durée sur des dizaines de sites. La comparaison classique a opposé neuf modèles de matière organique du sol à sept jeux de données de longue durée sans dégager de vainqueur unique : la performance dépendait du site et de la question8. Des travaux de modélisation d’ensemble plus récents parviennent au même verdict et ajoutent un avertissement : l’écart entre modèles est lui-même une source d’incertitude qu’une projection à modèle unique dissimule11. RothC et Century ont le plus long historique international ; AMG est la meilleure option pour les systèmes de grandes cultures français et autres européens ; DayCent est l’outil spécialisé lorsque les émissions de N₂O sont dans le périmètre.

Choisir le bon modèle pour la question

Quatre questions pratiques déterminent le choix.

L’azote est-il dans le périmètre ?

Si la décision dépend des émissions de N₂O, de l’efficience d’utilisation des engrais ou de l’interaction entre la gestion de l’azote et la dynamique du carbone, RothC est structurellement le mauvais choix ; il n’a pas de compartiment azote. Century ou DayCent sont les réponses honnêtes. Si la décision porte purement sur le carbone du sol sous une gestion stable de l’azote, RothC est défendable et plus simple.

Le système est-il limité par l’eau ?

Le modificateur hydrique par défaut de RothC donne de moins bons résultats dans les systèmes arides et méditerranéens, où l’activité biologique est limitée par l’humidité pendant une grande partie de l’année9. Une calibration modifiée de RothC (Farina et al.) existe et constitue le choix pratique dans ces contextes. Century et DayCent gèrent nativement la limitation hydrique grâce à leurs sous-modèles.

Quelle résolution temporelle est nécessaire ?

Le pas mensuel suffit pour une trajectoire décennale sur 20 ans. Le pas journalier est requis lorsque le résultat recherché est un flux rapide : des émissions de gaz traces après une application d’engrais, par exemple. RothC, AMG et Century fonctionnent au pas mensuel ; DayCent fonctionne au pas journalier.

Quelles données le projet peut-il réellement produire ?

Un modèle qui a besoin de données que le projet ne peut fournir est un modèle qui tournera sur des hypothèses. Si les données météorologiques journalières ne sont pas disponibles, DayCent fonctionnera sur une désagrégation, ce qui est défendable si la désagrégation est explicitée. Si les apports de résidus et l’historique de gestion sont inconnus, tous les modèles sont dans la même situation : l’incertitude issue des données manquantes domine celle issue du modèle.

Calibration, validation, et savoir quand un modèle convient

Un modèle qui n’a jamais été confronté à du carbone mesuré sur un système comparable est une hypothèse, pas une preuve. La calibration ajuste les paramètres d’un modèle pour reproduire un jeu de données connu ; la validation pose la question plus difficile de savoir s’il prédit ensuite des données qu’il n’a jamais vues. On confond couramment les deux, et un modèle qui épouse parfaitement ses données de calibration peut malgré tout échouer sur un site indépendant.

Les instruments de cette vérification sont les essais de longue durée, dont certains fonctionnent depuis plus d’un siècle, où le carbone du sol a été suivi sous traitements fixes. On rend compte de la qualité d’un modèle à l’aune de statistiques familières : la racine de l’erreur quadratique moyenne pour l’ampleur de l’écart, et l’efficience du modèle pour savoir s’il fait mieux que de simplement prédire la moyenne. La comparaison de neuf modèles qui a fait référence en la matière a trouvé un accord sur certains sites et une nette divergence sur d’autres, aucun cadre ne dominant partout8. Un ensemble de 26 modèles sur des sols en jachère nue a abouti à une conclusion plus tranchée encore : sur de longues projections, le choix du modèle peut peser autant que le choix de la gestion11.

Deux règles pratiques en découlent. Validez sur des séries temporelles issues de la région et du système que vous modélisez, et non sur la moyenne mondiale, car un modèle calibré sur des sols arables tempérés n’offre aucune garantie en agroforesterie tropicale. Et traitez le domaine d’applicabilité comme une frontière réelle : un modèle sommé d’extrapoler au-delà des conditions sur lesquelles il a été testé renvoie un chiffre à la confiance trompeuse, ce qui fait l’objet de Représentatif de quoi ? Lorsqu’un projet a besoin d’une comptabilité fondée sur un modèle au niveau où l’emploient les inventaires nationaux, le GIEC parle alors de niveau 3 (Tier 3) et exige précisément cette preuve : un modèle validé face à des mesures représentatives du système12. Les articles de synthèse qui suivent l’évolution du domaine font le même constat du côté de la recherche10.

L’initialisation : la source d’erreur silencieuse

Tout modèle de carbone du sol répartit le stock en compartiments aux vitesses de renouvellement différentes : rapide, lent et quasi inerte. La répartition initiale entre compartiments est rarement mesurable, et le plus souvent estimée. Pour RothC, l’équation de Falloon estime le compartiment inerte à partir du COS total ; c’est une approximation statistique dont l’effet sur les projections à 30 ans est loin d’être négligeable13.

L’équation ci-dessus est un ajustement statistique, pas une mesure physique. Ses paramètres portent une incertitude qui s’amplifie au fil des projections décennales. L’alternative plus défendable est le spin-up : faire tourner le modèle jusqu’à l’équilibre sous un régime de gestion historique reconstruit, puis utiliser la distribution de compartiments obtenue comme point de départ. La reconstruction est la partie difficile, et c’est aussi l’étape la plus souvent escamotée dans la documentation de projet. Une rétro-simulation historique sur 50 ans sous une rotation réaliste réduit l’effet de l’initialisation sur les trajectoires prévues d’un ordre de grandeur par rapport à une initialisation fondée sur les compartiments par défaut14.

L’outil ci-dessous montre pourquoi il ne s’agit pas d’un détail technique. Fixez un stock mesuré, fixez la gestion future, et ne faites varier que l’hypothèse retenue sur le compartiment inerte. La fraction inerte ne gagne ni ne perd jamais de carbone : en y plaçant davantage de stock, on en laisse moins qui puisse réagir. La même mesure, répartie de deux façons, peut produire un gain projeté sous une hypothèse et une perte projetée sous l’autre. C’est là toute la décision, tranchée par un nombre que personne n’a mesuré.

Lire les sorties comme des enveloppes, pas des lignes

Une trajectoire de carbone du sol assortie d’un étroit intervalle de confiance à 95 % reflète soit un système bien calibré et riche en données, soit une sous-déclaration. La propagation Monte-Carlo de l’incertitude des données d’entrée (variabilité climatique, teneur en argile, apports de résidus, initialisation) produit une enveloppe de trajectoire généralement bien plus large que la ligne moyenne. Cette enveloppe est le livrable honnête. Retirez-la et le chiffre qui reste n’est, au mieux, qu’une estimation centrale plausible.

Pour un travail de comparaison de scénarios (référence contre intervention, par exemple), la sortie utile est la différence entre scénarios, avec sa propre incertitude propagée. Une grande part du bruit d’entrée est partagée entre scénarios et s’annule dans la différence, souvent plus étroite que l’une ou l’autre trajectoire individuelle. C’est le bon cadrage pour la quasi-totalité des questions de faisabilité de projets carbone.

L’outil ci-dessous rend cela concret. Basculez entre l’affichage des deux scénarios séparément, chacun assorti d’une enveloppe assez large pour engloutir le signal, et l’affichage de leur différence, où le bruit climatique et pédologique partagé s’annule et où émerge un nombre défendable. Plus la référence et l’intervention partagent les mêmes conditions, plus cette différence se resserre. C’est la raison statistique pour laquelle les cadres de MRV sérieux récompensent les dispositifs appariés et les sites témoins plutôt que les revendications absolues de type avant-après15.

La frontière : ce que les modèles à compartiments laissent de côté

Les quatre modèles usuels descendent tous de la même idée : des compartiments conceptuels définis par leur vitesse de renouvellement. C’est une abstraction puissante, mais ces compartiments ne sont pas des objets que l’on peut aller mesurer. On ne peut pas tamiser un sol pour en séparer les fractions décomposables et humifiées ; ce sont des compartiments comptables déduits du comportement. Depuis vingt ans, cet écart nourrit une seconde génération de modèles bâtis sur des grandeurs mesurables.

Le déplacement de la science sous-jacente va de la chimie vers la biologie et la minéralogie. La persistance, selon la vision actuelle, tient moins à ce qu’une molécule soit intrinsèquement difficile à dégrader qu’à la capacité des microbes à l’atteindre et à celle des surfaces minérales à la retenir18. Les modèles à microbiologie explicite comme MIMICS et Millennial placent la biomasse microbienne et le carbone associé aux minéraux au centre, et remplacent les compartiments déduits par des fractions que l’on peut, en principe, mesurer en laboratoire1617. Ils réagissent différemment des modèles du premier ordre sous l’effet du réchauffement, et l’écart est loin d’être négligeable.

Rien de tout cela ne met RothC ou Century à la retraite pour le travail de projet. Les modèles plus récents sont plus difficiles à paramétrer, moins éprouvés face aux données de terrain de longue durée, et pas encore inscrits dans les référentiels de MRV. Mais ils indiquent où réside réellement l’incertitude d’une projection décennale, et le praticien qui prend la structure en compartiments pour une vérité physique plutôt que pour une fiction utile se fera surprendre par la prochaine décennie de littérature.

De la trajectoire à la revendication défendable

Une trajectoire ne devient un actif que lorsqu’un référentiel l’accepte, et les cadres qui comptent ici convergent vers une même posture : modéliser pour planifier et pour interpoler entre les mesures, mais mesurer pour revendiquer. La GHG Protocol Land Sector and Removals Guidance ne considère les absorptions modélisées comme déclarables que si elles s’accompagnent d’un programme de mesure et de suivi et d’une incertitude explicite1. Les inventaires nationaux qui recourent à une comptabilité de niveau 3 (Tier 3) fondée sur un modèle sont soumis à la même exigence de validation face à des données représentatives12. L’article de référence sur le MRV du carbone du sol est explicite : les modèles réduisent le coût d’échantillonnage mais ne suppriment pas le besoin de mesurer15.

Le livrable honnête d’une mission de modélisation est donc rarement un nombre unique. C’est une trajectoire assortie de son enveloppe, un énoncé clair du modèle employé et des raisons de ce choix, une initialisation documentée, et un plan des mesures qui la confirmeront ou la corrigeront. Des études à l’échelle continentale modélisent des millions d’hectares selon cette logique même19, et cela vaut aussi à l’échelle de la parcelle : le modèle porte la projection entre les points où le sol est effectivement échantillonné.

Points clés