Wie modelliert man die Veränderung des Bodenkohlenstoffs über die Zeit?

Ein Bodenkohlenstoffmodell nimmt eine Beschreibung eines Systems (Klima, Boden, Bewirtschaftung) entgegen und liefert eine Trajektorie: wie sich der Kohlenstoffvorrat unter dieser Beschreibung voraussichtlich entwickelt. Mit dieser Trajektorie begründen Organisationen die Kapitalallokation für regenerative Programme, prognostizieren Scope-3-Emissionen und -Entnahmen nach der GHG Protocol Land Sector and Removals Guidance und machen, wo die Projektstruktur es rechtfertigt, Ansprüche auf freiwilligen Kohlenstoffmärkten geltend.

Die Frage „Welches Modell sollen wir verwenden?“ hat keine allgemeingültige Antwort. Die richtige Frage ist enger gefasst: Welches Modell ist für die Entscheidung geeignet, für das vorliegende System, angesichts der Daten, die das Projekt realistischerweise beschaffen kann? Sie zu beantworten ist die wichtigste Ermessensentscheidung in einem Mandat zur Bodenkohlenstoffmodellierung.

Hinter den Akronymen teilen die gängigen Arbeitsmodelle einen Mechanismus: Der Bodenkohlenstoff wird in Pools aufgeteilt, die jeweils mit ihrer eigenen Rate zerfallen, gespeist von einem Strom frischer Einträge aus Pflanzenmaterial und Wirtschaftsdünger. Wer diesen Mechanismus versteht, für den beantworten sich die meisten praktischen Fragen von selbst, einschließlich der beiden, die eine Prognose über Jahrzehnte stillschweigend bestimmen: wo man den Ausgangspunkt der Pools annimmt und wie breit die ehrliche Antwort ist. Die drei Werkzeuge in diesem Beitrag lassen eine kompakte Version dieses Mechanismus laufen, sodass man sein Verhalten selbst bewegen kann, statt nur darüber zu lesen.

Was Modellierung tatsächlich leistet, und was nicht

Ein Bodenkohlenstoffmodell ist eine strukturierte Hypothese darüber, wie sich Kohlenstoff durch das System bewegt. Seine Ergebnisse sind bedingt: Bei diesen Eingangsdaten ist dies die erwartete Trajektorie und ihre Unsicherheit. Ein Modell misst den Boden nicht; es sagt ihn vorher.

Das ist wichtig, weil ein großer Teil der Meinungsverschiedenheiten zwischen Projektentwicklern und Verifizierern daher rührt, dass beides verwechselt wird1. Eine modellierte Trajektorie, die in der jährlichen Berichterstattung an die Stelle von Felddaten tritt, ist ein Compliance-Risiko. Eine modellierte Trajektorie, die neben Felddaten und mit ausgewiesenen Unsicherheitsspannen verwendet wird, ist gängige Praxis nach der GHG Protocol Land Sector and Removals Guidance und nach den meisten Methodologien des freiwilligen Marktes.

Wie diese Modelle tatsächlich funktionieren

Bodenkohlenstoff ist keine einzelne Substanz, die mit einer einzigen Rate zerfällt. Eine frische Wurzel und ein jahrhundertealter Huminstoffkomplex sind beide „organischer Bodenkohlenstoff“, doch der eine setzt sich in Monaten um, der andere überdauert die Person, die ihn misst. Jedes gebräuchliche Modell geht damit auf dieselbe Weise um: Es teilt den Kohlenstoff in eine Handvoll Pools mit jeweils eigener Umsatzzeit auf und verfolgt sie getrennt2.

Innerhalb jedes Pools folgt der Abbau einer Kinetik erster Ordnung: Die in einem bestimmten Monat verlorene Menge ist proportional zur vorhandenen Menge, sodass ein sich selbst überlassener Pool entlang einer Exponentialkurve zerfällt, die durch eine Ratenkonstante k bestimmt wird.

Die Ratenkonstante k ist eine Eigenschaft des Pools: hoch für frische Rückstände, niedrig für humifizierte Substanz, praktisch null für die inerte Fraktion. Die Funktion f skaliert k nach oben oder unten, je nach den Bedingungen, in denen sich der Boden tatsächlich befindet. Wärme beschleunigt den Abbau, Feuchtestress verlangsamt ihn, und eine lebende Pflanzendecke bremst ihn im Vergleich zu offenem Boden. Deshalb baut derselbe Rückstandseintrag unter einer Zwischenfrucht in einer kühlen gemäßigten Saison Kohlenstoff auf und verliert ihn in einer heißen, trockenen.

Nicht der gesamte abgebaute Kohlenstoff verlässt den Boden. Ein Teil wird als CO₂ in die Atmosphäre veratmet; der Rest wird in mikrobieller Biomasse und humifizierten Verbindungen stabilisiert, die wieder in die langsamen Pools einfließen. Der Tongehalt steuert diese Aufteilung, weil feiner texturierte Böden einen größeren Teil des umgesetzten Kohlenstoffs schützen. Deshalb ist die Tonfraktion in jedem der folgenden Modelle eine erforderliche Eingangsgröße3.

Setzt man diese Teile zusammen, bewegt sich ein Boden unter gleichbleibender Bewirtschaftung auf einen Gleichgewichtszustand zu, den Vorrat, bei dem sich die jährlichen Einträge und die jährlichen Verluste genau die Waage halten. Erhöht man die Einträge, steigt der Gleichgewichtszustand; erwärmt sich das Klima oder sinkt der Tongehalt, fällt er. Die Annäherung ist langsam, getaktet von der Umsatzzeit des langsamsten aktiven Pools, weshalb glaubwürdige Bodenkohlenstoffgewinne über Jahrzehnte statt über Jahre angegeben werden. Das Werkzeug unten lässt diesen Mechanismus vorwärts laufen. Verändern Sie Eintrag, Tongehalt und Temperatur, und beobachten Sie, wo sich der Vorrat einpendelt und wie lange es dauert, bis er dort ankommt.

Zwei Eigenschaften dieser Kurve sind für ein Projekt wichtig. Das Ziel wird durch die Einträge und den Standort bestimmt, nicht durch den Ausgangszustand des Bodens, sodass zwei Flächen unter derselben Bewirtschaftung aus unterschiedlichen Richtungen demselben Gleichgewichtszustand zustreben. Und der Weg ist lang: Selbst eine große Veränderung braucht Jahrzehnte, um sich auszuwirken, und die ersten Jahre bewegen sich langsam. Eine modellierte Trajektorie, die in fünf Jahren ein neues Gleichgewicht erreicht, beschreibt ein Artefakt ihrer Konfiguration, keinen Boden.

Die vier infrage kommenden Modelle

Vier Ansätze dominieren die praktische Arbeit. Jeder wurde für eine bestimmte Klasse von Fragen entwickelt, und jeder trägt den Fingerabdruck des Systems, an dem er aufgebaut wurde. RothC ist aus den Rothamsted-Langzeitversuchen an gemäßigten Ackeroberböden hervorgegangen2. AMG begann als Drei-Pool-Modell, das seine INRA-Autoren zuerst an Anbauversuchen in der argentinischen Pampa erprobten, und wurde in seiner heutigen Form auf Grundlage von sechzig französischen Langzeitfeldversuchen an zwanzig Standorten neu aufgebaut45. Century wurde entwickelt, um die Gehalte an organischer Substanz in den Grasländern der Great Plains zu erklären, mit an den Kohlenstoffkreislauf gekoppeltem Stickstoff6 (Phosphor und Schwefel folgten in späteren Versionen), und DayCent ist sein Gegenstück mit täglichem Zeitschritt, gebaut, um die Spurengasflüsse aufzulösen, die ein Monatsmodell nicht erfassen kann7. Eines davon außerhalb seines angestammten Bereichs zu verwenden, ist möglich, kostet aber Glaubwürdigkeit. Zur Prüfung, ob ein neuer Standort überhaupt innerhalb dieses Bereichs liegt, siehe Repräsentativ wofür?

Alle vier wurden an experimentellen Langzeitdaten von Dutzenden Standorten getestet. Der klassische Vergleich stellte neun Modelle der organischen Bodensubstanz sieben Langzeitdatensätzen gegenüber und fand keinen eindeutigen Sieger: Die Leistung hing vom Standort und von der Fragestellung ab8. Neuere Ensemble-Arbeiten kommen zum selben Urteil und fügen eine Warnung hinzu: Die Streuung zwischen den Modellen ist selbst eine Unsicherheitsquelle, die eine Projektion mit nur einem Modell verbirgt11. RothC und Century haben die längste internationale Erfolgsbilanz; AMG ist die stärkste Option für französische und andere europäische Ackerbausysteme; DayCent ist das Spezialwerkzeug, wenn N₂O-Emissionen im Untersuchungsrahmen liegen.

Das richtige Modell für die Fragestellung wählen

Vier praktische Fragen bestimmen die Wahl.

Liegt Stickstoff im Untersuchungsrahmen?

Hängt die Entscheidung von N₂O-Emissionen, der Düngernutzungseffizienz oder der Wechselwirkung zwischen Stickstoffmanagement und Kohlenstoffdynamik ab, ist RothC strukturell die falsche Wahl; es hat keinen N-Pool. Century oder DayCent sind die ehrlichen Antworten. Geht es bei der Entscheidung ausschließlich um Bodenkohlenstoff unter stabilem Stickstoffmanagement, ist RothC vertretbar und einfacher.

Ist das System wasserlimitiert?

Der Standard-Wassermodifikator von RothC schneidet in Trockengebieten und mediterranen Systemen schwächer ab, in denen die biologische Aktivität einen großen Teil des Jahres durch Feuchtigkeit begrenzt ist9. Eine modifizierte RothC-Kalibrierung (Farina et al.) existiert und ist in diesen Kontexten die praktische Wahl. Century und DayCent bilden die Wasserlimitierung über ihre Teilmodelle von Haus aus ab.

Welche zeitliche Auflösung wird benötigt?

Ein monatlicher Schritt genügt für eine Trajektorie über 20 Jahre. Ein täglicher Schritt ist erforderlich, wenn die interessierende Zielgröße ein schneller Fluss ist: zum Beispiel Spurengasemissionen nach einer Düngung. RothC, AMG und Century arbeiten monatlich; DayCent arbeitet täglich.

Welche Daten kann das Projekt realistischerweise erheben?

Ein Modell, das Eingangsdaten benötigt, die das Projekt nicht liefern kann, ist ein Modell, das auf Annahmen laufen wird. Sind tägliche Wetterdaten nicht verfügbar, arbeitet DayCent mit einem Downscaling, was vertretbar ist, sofern das Downscaling offengelegt wird. Sind Rückstandseinträge und Bewirtschaftungshistorie unbekannt, sind alle Modelle in derselben Lage: Die Unsicherheit aus fehlenden Eingangsdaten überwiegt die Unsicherheit aus dem Modell.

Kalibrierung, Validierung und wann ein Modell passt

Ein Modell, das nie an gemessenem Kohlenstoff in einem vergleichbaren System geprüft wurde, ist eine Hypothese, kein Beleg. Die Kalibrierung stimmt die Parameter eines Modells so ab, dass es einen bekannten Datensatz reproduziert; die Validierung stellt die schwierigere Frage, ob es anschließend Daten vorhersagt, die es nie gesehen hat. Beides wird regelmäßig verwechselt, und ein Modell, das seine Kalibrierungsdaten hervorragend abbildet, kann an einem unabhängigen Standort dennoch scheitern.

Die Instrumente dafür sind die Langzeitversuche, von denen einige seit mehr als einem Jahrhundert laufen und in denen der Bodenkohlenstoff unter festen Behandlungen verfolgt wurde. Die Modellgüte wird an ihnen mit vertrauten Statistiken angegeben: der Wurzel des mittleren quadratischen Fehlers für die Größe der Abweichung und der Modelleffizienz dafür, ob das Modell besser abschneidet, als einfach den Mittelwert vorherzusagen. Der Vergleich von neun Modellen, der für diese Arbeit den Maßstab setzte, fand an einigen Standorten Übereinstimmung und an anderen deutliche Abweichungen, ohne dass ein Ansatz überall dominierte8. Ein Ensemble aus 26 Modellen auf Schwarzbrache-Böden kam zu einem noch schärferen Schluss: Es gibt keine vertretbare Möglichkeit, das beste Modell im Voraus zu bestimmen, und es braucht ein Ensemble von etwa zehn Modellen, bevor sich die Schätzung stabilisiert11.

Daraus folgen zwei praktische Regeln. Validieren Sie an Zeitreihendaten aus der Region und dem System, das Sie modellieren, nicht am globalen Durchschnitt, denn ein an gemäßigten Ackerböden kalibriertes Modell bietet keine Garantie für tropische Agroforstsysteme. Und behandeln Sie den Anwendbarkeitsbereich als reale Grenze: Ein Modell, das über die Bedingungen hinaus extrapolieren soll, unter denen es getestet wurde, liefert eine Zahl mit falscher Zuversicht, worum es in Repräsentativ wofür? geht. Braucht ein Projekt eine modellbasierte Bilanzierung auf dem Niveau, das nationale Inventare verwenden, spricht das IPCC von Tier 3 und verlangt genau diesen Beleg: ein Modell, das an für das System repräsentativen Messungen validiert wurde12. Die Übersichtsarbeiten, die die Entwicklung des Fachgebiets verfolgen, kommen von der Forschungsseite her zum selben Schluss10.

Dieser Abschnitt ist eine Einführung in ein Thema, das mehr umfasst, als ein einzelner Abschnitt fassen kann: Das Wort Kalibrierung trägt vier verschiedene Bedeutungen, die Güte innerhalb der Stichprobe schmeichelt, und das gewählte Kreuzvalidierungsschema ist selbst eine Aussage über die Vorhersageaufgabe. Ein Feldleitfaden zur Kalibrierung und Validierung eines Modells behandelt sie einzeln.

Initialisierung: die stille Fehlerquelle

Jedes Bodenkohlenstoffmodell teilt den Vorrat in Pools mit unterschiedlichen Umsatzraten auf: schnell, langsam und nahezu inert. Die anfängliche Aufteilung auf die Pools ist selten messbar und wird meist geschätzt. Für RothC schätzt die Falloon-Gleichung den inerten Pool aus dem gesamten SOC13; das ist eine statistische Näherung, und die angenommene inerte Fraktion wirkt sich erheblich auf Projektionen über 30 Jahre aus, wie das Werkzeug unten zeigt.

Die obige Gleichung ist eine statistische Anpassung, keine physikalische Messung. Ihre Parameter tragen eine Unsicherheit, die sich über Projektionen von Jahrzehnten aufschaukelt. Die besser vertretbare Alternative ist der Spin-up: das Modell unter einem rekonstruierten historischen Bewirtschaftungsregime bis zum Gleichgewicht laufen lassen und die resultierende Poolverteilung dann als Ausgangspunkt verwenden. Die Rekonstruktion ist der schwierige Teil und zugleich der Schritt, der in der Projektdokumentation am häufigsten übergangen wird. Die Initialisierung über eine historische Rückrechnung unter einer realistischen Fruchtfolge statt über Standardpools ist der etablierte Weg, die anfängliche Aufteilung vertretbar zu machen14.

Das Werkzeug unten zeigt, warum dies keine Formalität ist. Legen Sie einen gemessenen Vorrat fest, legen Sie die künftige Bewirtschaftung fest, und ändern Sie nur die Annahme über den inerten Pool. Die inerte Fraktion gewinnt oder verliert nie Kohlenstoff; legt man mehr vom Vorrat dorthin, bleibt weniger übrig, das reagieren kann. Dieselbe Messung, auf zwei Arten aufgeteilt, kann unter der einen Annahme einen projizierten Gewinn und unter der anderen einen projizierten Verlust ergeben. Das ist die gesamte Entscheidung, entschieden durch eine Zahl, die niemand gemessen hat.

Ergebnisse als Bänder lesen, nicht als Linien

Eine Bodenkohlenstofftrajektorie mit einem engen 95-Prozent-Konfidenzintervall spiegelt entweder ein gut kalibriertes, datenreiches System wider oder eine zu geringe Ausweisung der Unsicherheit. Die Monte-Carlo-Fortpflanzung der Unsicherheit der Eingangsdaten (Klimavariabilität, Tongehalt, Rückstandseinträge, Initialisierung) ergibt ein Trajektorienband, das meist deutlich breiter ist als die Mittelwertlinie. Dieses Band ist das ehrliche Ergebnis. Entfernt man es, ist die verbleibende Zahl bestenfalls eine plausible zentrale Schätzung.

Für Szenarienvergleiche (zum Beispiel Baseline gegenüber Maßnahme) ist das nützliche Ergebnis die Differenz zwischen den Szenarien, mit ihrer eigenen fortgepflanzten Unsicherheit. Ein großer Teil des Rauschens in den Eingangsdaten ist den Szenarien gemeinsam und hebt sich in der Differenz auf, die oft enger ist als jede der einzelnen Trajektorien. Das ist die richtige Betrachtungsweise für nahezu alle Machbarkeitsfragen von Kohlenstoffprojekten.

Das Werkzeug unten macht das greifbar. Wechseln Sie zwischen der getrennten Darstellung der beiden Szenarien, jedes mit einem Band, das breit genug ist, um das Signal zu verschlucken, und der Darstellung ihrer Differenz, in der sich das gemeinsame Klima- und Bodenrauschen aufhebt und eine vertretbare Zahl hervortritt. Je mehr Baseline und Maßnahme dieselben Bedingungen teilen, desto enger wird diese Differenz. Das ist der statistische Grund, warum seriöse MRV-Rahmenwerke gepaarte Designs und Kontrollflächen gegenüber absoluten Vorher-Nachher-Aussagen bevorzugen15.

Die Forschungsfront: was die Poolmodelle auslassen

Die vier gängigen Arbeitsmodelle gehen alle auf dieselbe Idee zurück: konzeptionelle Pools, die durch ihre Umsatzrate definiert sind. Das ist eine wirkungsvolle Abstraktion, doch die Pools sind keine Dinge, die man messen gehen kann. Man kann einen Boden nicht in seine abbaubaren und humifizierten Fraktionen aussieben; es sind buchhalterische Kompartimente, die aus dem Verhalten erschlossen werden. Seit zwanzig Jahren treibt diese Lücke eine zweite Generation von Modellen an, die auf messbaren Größen aufbauen.

Die Verschiebung in der zugrunde liegenden Wissenschaft führt von der Chemie zur Biologie und Mineralogie. Persistenz hängt nach heutiger Auffassung weniger davon ab, dass ein Molekül an sich schwer abbaubar ist, als davon, ob Mikroben es erreichen können und ob Mineraloberflächen es festhalten18. Mikrobiell explizite Modelle wie MIMICS und Millennial stellen mikrobielle Biomasse und mineralassoziierten Kohlenstoff in den Mittelpunkt und ersetzen erschlossene Pools durch Fraktionen, die sich grundsätzlich im Labor messen lassen1617. Sie reagieren unter Erwärmung anders als Modelle erster Ordnung, und der Unterschied ist nicht klein.

Nichts davon schickt RothC oder Century für die Projektarbeit in den Ruhestand. Die neueren Modelle sind schwerer zu parametrisieren, weniger an Langzeitfelddaten getestet und noch nicht in MRV-Standards verankert. Doch sie zeigen, wo die Unsicherheit einer Projektion über Jahrzehnte tatsächlich liegt, und wer die Poolstruktur als physikalische Wahrheit statt als nützliche Fiktion behandelt, wird vom nächsten Jahrzehnt der Literatur überrascht werden.

Von der Trajektorie zur belastbaren Aussage

Eine Trajektorie wird erst dann zu einem Vermögenswert, wenn ein Standard sie akzeptiert, und die hier maßgeblichen Rahmenwerke laufen auf eine Haltung hinaus: modellieren, um zu planen und zwischen Messungen zu interpolieren, aber messen, um Ansprüche geltend zu machen. Die GHG Protocol Land Sector and Removals Guidance betrachtet modellierte Entnahmen nur dann als berichtsfähig, wenn sie auf Primärdaten beruhen und mit einem laufenden Monitoringprogramm einhergehen, und erwartet, dass die Unsicherheit bewertet und berichtet wird1. Nationale Inventare, die eine modellbasierte Tier-3-Bilanzierung verwenden, unterliegen derselben Anforderung der Validierung an repräsentativen Daten12. Die Referenzarbeit zum MRV von Bodenkohlenstoff stellt ausdrücklich fest, dass Modelle die Beprobungskosten senken, die Notwendigkeit zu messen aber nicht aufheben15.

Das ehrliche Ergebnis eines Modellierungsmandats ist daher selten eine einzelne Zahl. Es ist eine Trajektorie mit ihrem Band, eine klare Angabe, welches Modell verwendet wurde und warum, eine dokumentierte Initialisierung und ein Plan für die Messungen, die sie bestätigen oder korrigieren werden. Studien im kontinentalen Maßstab modellieren Millionen Hektar nach genau dieser Logik19, und sie gilt auch im Maßstab der einzelnen Fläche: Das Modell trägt die Projektion zwischen den Punkten, an denen der Boden tatsächlich beprobt wird.

Kernaussagen