[latexpage]

Es gibt Forschungsansätze zum CO$_2$ Senkeneffekt, die von der unbestrittenen kurzen Residenzzeit des CO$_2$ ausgehen, die in der Größenordnung von 3 Jahren liegt. Das prominenteste dieser Modelle ist das von Prof. Harde, beschrieben in einem Artikel aus 2019. Der Artikel geht insbesondere auf Gegenargumente ein, die gegen eine frühere Publikation von ihm vorgebracht wurde.

Die Hauptgleichung des Modells von Prof. Harde ist (Gleichung 23) – es werden die im Artikel verwendeten und erklärten Original-Bezeichnungen verwendet:

$ \frac{dC(t)}{dt}=e_{N0}+\beta_e\cdot\Delta T(t) + e_A(t) – \frac{C(t)}{\tau_{R0}+\beta_{\tau}\cdot\Delta T(t)} $

Prof. Harde schreibt, dass die Bestimmung aller Parameter mit den vorhandenen Daten unterbestimmt ist, und dass nur eine Kombination zwischen $\beta_e$ und $\beta_{\tau}$ bestimmt werden kann, wobei im Prinzip auch einer der beiden Koeffizienten gegen 0 gehen kann. Daher nehmen wir vorläufig an, dass $\beta_{\tau}=0$ ist, was bedeutet, dass sich die konzentrationsabhängige Absorptionskonstante nicht mit der Zeit ändert.

Damit vereinfacht sich die Differentialgleichung:
$ \frac{dC(t)}{dt}= e_A(t) + e_{N0}+\beta_e\cdot\Delta T(t) – \frac{1}{\tau_{R0}}\cdot C(t) $
Auf den ersten Blick ist das eine direkte „Übersetzung“ der diskreten Gleichung des temperaturabhängigen linearen Modells (Gleichung 2) – auch hier werden die dort verwendeten Symbole angewandt:
$C_{i+1} – C_i = E_i – a\cdot C_i – b\cdot T_i – c $
Auf den ersten Blick haben alle Terme klare Entsprechungen:

• Konzentrationswachstum $C_{i+1} – C_i \Longleftrightarrow \frac{dC(t)}{dt}$
• Anthropogene Emissionen $E_i \Longleftrightarrow e_A(t)$
• Konzentrationsabhängige Senkenwirkung $a\cdot C_i\Longleftrightarrow \frac{1}{\tau_{R0}}\cdot C(t) $
• Temperaturabhängige natürliche Emissionen $-b\cdot T_i\Longleftrightarrow \beta_e\cdot\Delta T(t)$
• Als konstant angenommene natürliche Emissionen $-c\Longleftrightarrow e_{N0}$

Aber die Konstanten $c$ und $e_{N0}$ passen nicht zusammen. $c$ wird als Schätzwert der jährlichen natürlichen Netto-Emissionen mit 13.6 ppm/Jahr ermittelt, während $e_{N0}$ mit 93.3 ppm/Jahr als Vorwissen aller natürlichen Emissionen vorgegeben wird. Demzufolge beschreibt $\tau_{R0}$ auch, wie beschrieben, die Residenzzeit von etwa 3 Jahren, während $\frac{1}{a} \approx 23 a $ in der Größenordnung der Zeitkonstante der Absorptionszeit der Bombentest-Daten (etwa 16 Jahre) liegt. Neuere Berechnungen mit monatsgenauen Messdaten ergeben mit dem temperaturabhängigen linearen Modell $a=0.06$, also $\frac{1}{a} \approx 17 a $.

Um die Modelle inhaltlich aufeinander abzustimmen müssen die „Brutto“ natürlichen Emissionen $e_{N0}$ aufgeteilt werden in die Jahres-„Netto“-Emissionen $-c$ und den restlichen, innerhalb des Jahres re-absorbierten natürlichen Emissionen $e_{Rest}$ (dass dabei ein Austausch mit anderem CO$_2$ stattfindet, steht außer Frage, es geht hier nur um die Menge)
$e_{N0} = e_{Rest} – c $
sowie deren zugehörige Temperatur-Koeffizienten
$\beta_e = \beta_{Rest} – b$
Ebenso müssen die Kehrwerte der Zeitkonstanten entsprechend angepasst werden.
$\frac{1}{\tau_{R0}} = \frac{1}{\tau_{Rest}} + a$
Um die beiden Ansätze kompatibel zu bekommen, muß diese Bedingung erfüllt sein:
$e_{Rest} +\beta_{Rest}\cdot\Delta T(t)-\frac{1}{\tau_{Rest}}\cdot C(t) = 0$
und die Temperaturkoeffizienten sollten im gleichen Verhältnis stehen wie die Emissionen:
$\frac{beta_{Rest}}{\beta_e} = \frac{e_{Rest}}{e_{N0}}$

Was hilft uns das?

Es erscheint mir wichtig, zu einem grundsätzlich gemeinsamen Verständnis zu kommen, um nicht unnötige Fronten in der Wissenschaft aufzurichten, denn leider gibt es sehr viele Mißverständnisse um die Modelle des Kohlenstoffzyklus.

Beide Konzepte, Residenzzeit und Ausgleichszeit, haben in der Diskussion eine große Rolle gespielt. Sie werden häufig gegeneinander ausgespielt, und daher kann es zum Dialog wesentlich beitragen, die Zusammenhänge beider Konzepte zu verstehen. Dann können die Mosaiksteine des Kohlenstoffzyklus zusammengefügt werden.

Die wichtigste Konsequenz ist zu zeigen, dass die Gleichungen zur Behandlung der der Residenzzeit kompatibel ist zur Gleichung zur Bestimmung der Ausgleichszeit, und dass beide Fragestellungen nicht gegeneinander ausgespielt werden dürfen:

• Es ist legitim zu fragen, wie sich anthropogene Emissionen auf die Konzentration auswirken, und es ist wenig hilfreich, diese Frage damit herunterzuspielen, dass es vielfach größere natürliche Emissionen gibt
• Es ist auch legitim zu sagen, dass trotz des jährlichen Netto-Senkeneffekts eine Temperaturabhängigkeit der natürlichen Emissionen gibt
• Es ist wert zu betonen, dass sich die lineare Temperaturabhängigkeit der Emissionen und die lineare Konzentrationsabhängigkeit der Absorptionen teilweise kompensieren, dass die Gesamtwirkung wie eine Temperaturunabhängigkeit des Senkeneffekt bedeuten, solange die Temperatur linear mit der CO$_2$ Konzentration ansteigt.

[latexpage]

In der Klimadiskussion wird das Argument der Rückkopplung durch Wasserdampf dazu herangezogen, um die Klimawirkung der Treibhausgase — die Sensitivität bei Verdoppelung von deren Konzentration in der Atmosphäre — , die nach Strahlungstransportgleichung und generellem Konsens maximal 0,8° beträgt, um einen angeblichen Faktor 2-6 zu verstärken. Allerdings wird das gewöhnlich nicht genauer quantifiziert, es werden in der Regel nur Formeln mit dem „finalen Feedback“ angegeben.

Vor kurzem haben David Coe, Walter Fabinski und Gerhard Wiegleb in der Publikation „The Impact of CO₂, H₂O and Other ‚Greenhouse Gases‘ on Equilibrium Earth Temperatures“ unter anderem genau diese Rückkopplung beschrieben und analysiert. In Anlehnung an ihre Publikation wird dieser Effekt mit dem teils gleichen, teils leicht unterschiedlichen Ansatz im folgenden hergeleitet. Die Ergebnisse sind fast identisch.

Dabei wird hier von allen anderen Effekten, die bei der Bildung von Wasserdampf auftreten, wie z.B. Wolkenbildung, abgesehen.

Der grundsätzliche Mechanismus der Wasserdampfrückkopplung

Ausgangspunkt ist eine Temperaturerhöhung der Atmosphäre um ∆T₀, ungeachtet deren Ursache. Typischerweise wird der Treibhauseffekt als primäre Ursache angenommen. Die Argumentation ist nun, dass die erwärmte Atmosphäre mehr Wasserdampf aufnehmen kann, d.h. der Sättigungsdampfdruck (SVP = „Saturation Water Pressure“) erhöht sich und es wird angenommen, dass sich konsequenterweise auch die Wasserdampfkonzentration ∆H₂O erhöht, und zwar als lineare Funktion der Temperaturänderung. (Die Temperaturänderung ist so klein, dass eine Linearisierung auf jeden Fall legitim ist):
$\Delta H_2O = j\cdot \Delta T_0 $
Dabei ist $j$ die Proportionalitätskonstante für die Wasserdampfkonzentration.
Eine erhöhte Wasserdampfkonzentration bewirkt wiederum aufgrund der Treibhauswirkung von Wasserdampf eine Temperaturerhöhung, die linear von der Wasserdampfkonzentration abhängt:
$\Delta T_1 = k\cdot \Delta H_2O $
Zusammengefaßt bewirkt also die auslösende Temperaturerhöhung ∆T₀ eine Folgeerhöhung der Temperatur ∆T₁:
$\Delta T_1 = j\cdot k\cdot \Delta T_0 $
Da die Voraussetzung des Verfahrens ist, dass die Ursache der auslösenden Temperaturerhöhung unerheblich ist, bewirkt die Erhöhung um ∆T₁ natürlich ebenfalls wieder einen Rückkopplungszyklus:
$\Delta T_2 = j\cdot k\cdot \Delta T_1 = (j\cdot k)^2\cdot \Delta T_0$
Dies wiederholt sich rekursiv. Die finale Temperaturänderung ist demnach eine geometrische Reihe:
$\Delta T = \Delta T_0\sum_{n=0}^\infty(j\cdot k)^n = \Delta T_0\cdot \frac{1}{1-j\cdot k} $
Wäre $j\cdot k\ge 1$, würde die Reihe divergieren und die Temperatur über alle Grenzen wachsen. Daher ist es wichtig, sich über die Größe dieser beiden Rückkopplungsfaktoren Klarheit zu verschaffen.

Abhängigkeit der möglichen Wasserdampfkonzentration von der Temperatur

Die maximal mögliche Wasserdampfkonzentratio in Abhängigkeit von der Temperatur T (in °C) ist durch den Sättigungsdampfdruck SVP (englisch „saturation vapour pressure“, SVP) begrenzt. Dieser wird durch die Arden Buck Gleichung, (eine moderne, überarbeitete Version der Magnus-Formel) sehr genau beschrieben:
$ SVP = 0.61121\cdot \exp{((18.678-\frac{T}{234.5})(\frac{T}{257.14+T}))} $
Es wird hier die Standard-Atmosphäre mit 15°C Boden- bzw. Wasseroberflächentemperatur und adiabatischem Temperaturgradient von -6.5°C/km betrachtet.

Die absolute Differenz $\frac{\Delta (SVP(T))}{\Delta T}$ ist naturgemäß bei höheren Temperaturen, also in Bodennähe, am größten:

Die relative Differenz $\frac{\frac{\Delta (SVP(T))}{\Delta T}}{SVP(T)}$ wird mit zunehmender Höhe größer, bewegt sich zwischen 4% und 8%.

Die mögliche Zunahme der relativen Luftfeuchtigkeit – das Verhältnis des tatsächlichen Dampfdrucks im Vergleich zum Sättigungsdampfdrucks – als Folge der globalen Temperaturerhöhung $T_0$ ist durch diese relative Änderung des Sättigungsdampfdrucks begrenzt.

Da die mittlere, dominante Infrarot-Abstrahlung der Erde etwa in der Höhe 5000m stattfindet, und sich oberhalb davon kaum mehr Wasserdampf befindet, ist es sinnvoll, 6% als oberes Limit der Änderung der relativen Luftfeuchtigkeit infolge einer Temperaturerhöhung um 1°C anzunehmen. Demzufolge ergibt sich die Konstante $j$ als $j=0.06$. Dieser Wert ist etwas kleiner als die üblicherweise genannten (aber gewöhnlich nicht belegten) 7%. Nach dem obigen Diagramm wären 7% Erhöhung der Luftfeuchtigkeit erst oberhalb von 8000 m ü.d.M. möglich.

Abhängigkeit des Treibhauseffekts von der Änderung der relativen Luftfeuchtigkeit

Der Infrarot Strahlungstransport in der Atmosphäre ist von der relativen Luftfeuchtigkeit abhängig. Dies wird in dem bekannten und bewährten Simulationsprogram MODTRAN berücksichtigt. Mit zunehmender Luftfeuchtigkeit sinkt infolge des Treibhauseffektes des Wasserdampfes die ausgehende Infrarotstrahlung.

Zwischen der Luftfeuchtigkeit 60% und 100% ist die Strahlungsabnahme linear. Daher wird zur Ermittlung der Abnahme der Strahlungsleistung und der zur Kompensation notwendigen Temperaturerhöhung die Zunahme der relativen Luftfeuchtigkeit von 80% auf 86% betrachtet.

Dazu stellen wir die Parameter der MODTRAN Simulation auf

• die aktuelle CO2-Konzentration von 420 ppm,
• eine relative Luftfeuchtigkeit von 80%,
• und eine Wolkenkonstellation, die der mittleren IR Abstrahlungsleistung von 240 $\frac{W}{m^2}$ nahe kommt.

Der Temperatur-Offset wird nun so lange vergrößert, bis die reduzierte iR-Abstrahlung von 0.6 \frac{W}{m^2} durch Temperaturerhöhung wieder ausgeglichen ist. Dies ist bei einer Erhöhung der Bodentemperatur um 0.185 °C der Fall.

Eine 6% höhere relative Luftfeuchtigkeit bewirkt also einen Treibhauseffekt, der durch eine Temperaturerhöhung von 0.185°C ausgeglichen wird, Auf eine Änderung um (theoretische) 100% Luftfeuchtigkeit hochgerechnet sind das $k=3.08$°C/100% .

Der finale Rückkopplungsfaktor und der gesamte Treibhauseffekt

Damit bewirkt eine um 1 Grad höhere Temperatur in einem Rückkopplungszyklus eine zusätzliche Temperaturerhöhung um $k\cdot j = 0.06*3.08= 0.185$.

Die geometrische Reihe führt zu einem Verstärkungsfaktor $f$ des reinen CO$_2$ Treibhauseffekts um
$f=\frac{1}{1-0.185} = 1.23 $

Damit ist die um die Wasserdampfrückkopplung verstärkte Sensitivität bei Verdopplung der CO$_2$ Konzentration $\Delta T$ nicht mehr $\Delta T_0=0.8$°C, sondern
$\Delta T = 1.23\cdot 0.8$ °C = 0.98°C $\approx$ 1°C

Würde man die vom „Mainstream“ postulierten 7% maximaler Erhöhung der Luftfeuchtigkeit zugrunde legen, wäre die zusätzliche Temperaturerhöhung 0.215 °C und demzufolge die Verstärkung des Treibhauseffekts
$f=\frac{1}{1-0.215} = 1.27 $. Die Sensitivität für CO$_2$ Verdoppelung wäre dann
$\Delta T = 1.27\cdot 0.8$ °C = 1.02°C $\approx$ 1°C

Dieses Ergebnis berücksichtigt nicht die um die durch höhere Wasserdampfkonzentration stärkere Wolkenbildung und deren Abschirmung des einfallenden Sonnenlichts, die eine negative Rückkopplung bewirkt.