Inhaltsverzeichnis

Gänseblümchen

Um diesen Abschnitt verstehen zu können, sollte das Kapitel Einführung in die Numerik zuerst behandelt werden! Dieser Abschnitt eignet sich auch als fächerübergreifende Einheit in Mathematik, Biologie (theoretische Ökologie) oder Physik.

Die schwarze Variante ist relativ selten, da sie in der Natur nicht vorkommt. Dennoch findet man mitunter ihr Bild im Internet. :-)
Ein Gänseblümchen in der weißen Variante. Ein Gänseblümchen in der schwarzen Variante.

Was ist die Gaia-Theorie?

Daisyworld (Gänseblümchen-Welt) ist ein theoretisches Modell, das von Lovelock und Watson entwickelt wurde, um die Gaia-Theorie zu untermauern. Die Gaia-Theorie wurde von der Mikrobiologin Lynn Margulis und dem Physiker James Lovelock entwickelt.

Hören wir uns an, was MEPHISTOPHELES dazu sagt:

Wer will was Lebendiges erkennen und beschreiben,
Sucht erst den Geist heraus zu treiben,
Dann hat er die Teile in seiner Hand,
Fehlt, leider! nur das geistige Band.
Goethe - Faust I


Genau diese Vorgehensweise, die ein zu untersuchendes System in seine Einzelteile zerlegt, um eben diese dann auch einzeln zu untersuchen, kritisiert Lovelock. Die Erforschung emergenter Eigenschaften - das sind Eigenschaften, die nur im komplexen Zusammenspiel aller Einzelteile in Erscheinung treten - ist mit diesem Ansatz nicht möglich. Lovelock und Margulis Augenmerk galt deswegen dem Gesamtsystem. Unter Gesamtsystem wird dabei gleich die Erde inklusiv allem Leben auf der Erde verstanden. Nach der Gaia-Theorie ist die Erde selbst ein lebender Organismus. Die Erde bildet zusammen mit allen Lebewesen ein komplex interagierendes Netzwerk. Die moderne Physik hat erst kürzlich damit begonnen, Eigenschaften interagierender Systeme systematisch zu untersuchen. Forschung im Bereich der Chaostheorie, selbst organisierte Kritizität, Fraktale, neuronale Netze, nicht lineare Systeme etc… oder einfach die Physik der komplexen Systeme sind hoch aktuell. Lovelock ist mit seiner Theorie teilweise auf erhebliche Kritik seiner Kollegen gestoßen. Um genauer zu erklären, was Lovelock eigentlich meinte, entwickelte er die Daisyworld.

Daisyworld

Die Daisyworld ist ein relativ einfaches komplexes System. Es besteht lediglich aus vier Komponenten:

  1. Eine Sonne, die dank ihrer Wärme (abgestrahlte Energie) Leben spenden kann.
  2. Ein Planet, der die Sonne umkreist.
  3. Eine Population von weißen Gänseblumen.
  4. Eine Population von schwarzen Gänseblumen.


Der Modellcharakter kommt schon darin zum Ausdruck, dass es schwarze Gänseblumen gar nicht gibt. Dennoch gelang es Lovelock und Watson mit diesem Modell große Aufmerksamkeit auf sich zu ziehen.

Damit Leben auf dem Planeten existieren kann, ist es notwendig, dass die Temperatur des Planeten in einem Bereich liegt, der es möglich macht, dass sich Leben entwickeln kann (Abbildung 1). Ein Planet ohne Leben wird umso heißer, je größer die Sonneneinstrahlung wird.
In der Entwicklungsgeschichte der Sonne, hat die Sonnenleuchtkraft schon um ca: 25% zugenommen, ohne dass sich die Temperatur der Erde dabei wesentlich erhöht hat. Wie lässt sich das erklären?
Lovelock u. Watson zeigten in ihrem Modell, dass es möglich ist, die Temperatur in Abhängigkeit von der Sonneneinstrahlung in einem großen Bereich konstant zu halten. Notwendig dafür ist die Existenz zweier Gänseblumenarten (Abbildung 2), die Aufgrund ihrer unterschiedlichen Eigenschaften in der Lage sind, die Gesamttemperatur des Planeten in der Nähe ihrer eigenen optimalen Wachstumstemperatur zu stabilisieren. Mit lediglich einer Blumenart war es nicht möglich, die Temperatur des Planeten in Abhängigkeit der sich veränderten Sonneneinstrahlung so kontrolliert konstant zu halten. Als Lovelock dann zusätzliche Gänseblumenarten und Lebewesen (Füchse, Kaninchen…) hinzufügte, machte er die Entdeckung, dass sich das System sehr stabil verhält. Dieses wertete er als starkes Plädoyer für die Artenvielfalt.
Bevor Lovelock zeigen konnte, dass die Regulierung der Planetentemperatur, eine wie er sagte System emergente Eigenschaft ist und sich auf dem interagierenden System ergibt, wurde der Gaia-Theorie vorgeworfen, sie sei teleologisch.
Als teleologisch könnte hier z.b. gelten, dass wie durch ein Wunder die Temperatur der Erde gerade derjenigen entspricht, die notwendig ist, damit sich Leben entwickeln kann.

Abbildung 1: Die Grafik zeigt den Einfluss der Sonneneinstrahlung auf die Temperatur eines Planeten im Daisyworld Modell. Die hellgraue Kurve ist der Temperaturverlauf ohne Leben. Hier gilt. umso höher die Sonneneinstrahlung ist, desto höher ist auch die Temperatur des Planeten. Damit findet in diesem Szenario keine Selbstregulation statt. Gänzlich anders verhält es sich (grüne Kurve), wenn wir Leben zulassen. Leben heißt hier, dass zwei Gänseblumenarten, die um den zur Verfügung stehenden Boden in Konkurrenz stehen, existieren. Mit Leben entwickelt sich nun das Gesamtsystem (Planet + Blumenpopulationen) in einem Zustand, der es dem Planet ermöglicht, fast unabhängig von der Sonneneinstrahlung die Temperatur nahe der optimalen Wachstumstemperatur der beiden Blumenarten zu halten. Bei einer Sonneneinstrahlung von ≈1.6 ist es für die Blumen zu heiß und sie sterben aus. Ab diesem Wert entspricht der Temperaturverlauf daher dem ohne Leben.
Abbildung 2: Dargestellt ist die Änderung der beiden Gänseblumenpopulationen in Abhängigkeit der Sonneneinstrahlung. Wenn es kalt ist, ist der Anteil der schwarzen Gänseblumen wesentlich höher als der Anteil der weißen Blumen. Mit zunehmender Sonneneinstrahlung schwindet der Anteil der schwarzen Gänseblumen zugunsten der weißen Gänseblumen. Ab einer Sonneneinstrahlung von 1.6 ist es für beide Populationen nicht mehr möglich zu überleben. Damit ist die Selbstregulation des Gesamtsystems (Planet + Blumenarten) zusammengebrochen.

Der nächste Abschnitt erklärt die physikalischen Grundlagen, die notwendig sind, um das Modell zu verstehen.

Physikalische Grundlagen

Dieser Abschnitt erklärt die Albedo und das Stefan-Bolzmann-Gesetz. Beides spielt in der Daisyworld eine wichtige Rolle.

Die Albedo

Als Albedo (Rückstrahlvermögen) bezeichnet man die Fähigkeit eines Körpers, auf ihn einfallende Strahlung wieder zu reflektieren. Die Albedo ist damit ein einfacher Zahlenwert zwischen 0 und 1.



Frisch gefallener Schnee hat eine Albedo von 0.8 ~ 0.90, d.h. bis zu 90% der einfallenden Strahlung werden wieder reflektiert. Dadurch heizt sich der Körper (in diesem Fall der Schnee) nicht so schnell auf.

Dunkler Asphalt hat eine Albedo von ~ 0.15, d.h lediglich 15% der einfallenden Strahlung wird reflektiert. Deswegen ist der Asphalt im Sommer natürlich auch recht warm.
Allgemein gilt:

Dunkle Körper haben eine niedrige Albedo und sind wärmer, als helle Körper.
$\longleftarrow$ Albedo fällt ab (Körper werden wärmer)
$\longrightarrow$ Albedo steigt an ( Körper werden kälter)

Im Daisyworld-Modell müssen die Albedos der weißen Gänseblumen, der noch unbewachsenen Fläche der Erde und der schwarzen Blumen berücksichtigt werden.

Bezeichnen wir die Albedo der weißen Gänseblumen mit $A_w$ , die der unbewachsenen Fläche mit $A_g$ und die der schwarzen Blumen mit $A_b$ dann folgt unabhängig von den tatsächlichen Werten für die Albedo aus dem bereits gesagten: \begin{equation} A_b < A_g < A_w \end{equation}

Die Indizes stehen für (b=black, g=ground und w=white).

Das Stefan-Bolzmann-Gesetz

  • E = Energie, hat die Einheit J (Joule)
  • P = Leistung (engl. Power) ist Energie/Zeit, hat die Einheit Watt: 1 W = 1 J/s.
  • I = Intensität, die auf eine Fläche bezogene Leistung (Energiefluss) , hat die Einheit $W/m^2$ bzw. $ J/(s\cdot m^2)$.

Jeder Körper strahlt Energie ab. Dabei gilt:
Umso heißer der Körper ist, desto mehr Energie wird abgestrahlt.
Das Stefan-Bolzmann-Gesetz beschreibt nun den genauen Zusammenhang zwischen der Energieabstrahlung eines Körpers und die dafür notwendige Temperatur. Interessanterweise hängt die abgestrahlte Energie sehr empfindlich von der Temperatur ab. Verdoppelt man z.B. die Temperatur eines Körpers (streng genommen werden ideale Körper betrachtet, worauf wir aber hier nicht eingehen wollen) dann wird die abgestrahlte Energie 16-mal so groß.

Das Stefan-Bolzmann-Gesetz lautet: \begin{equation} I = \sigma \cdot T^4 \end{equation}

Hier ist I die abgestrahlte Energie pro Zeit und Fläche (Einheit: $W/m^2$). T ist die Temperatur in Kelvin, sie geht mit der vierten Potenz ein. Und σ eine Konstante mit dem Zahlenwert σ=$5.67032\cdot 10^{-8}\; W/(K^4m^2)$. Diese Konstante heißt auch Stefan-Bolzmann-Konstante.

Berechnung der Oberflächen Temperatur der Sonne.

Auf einen Quadratmeter Erde erreicht uns in jeder Sekunde eine Sonnenenergie von 1367 J. Dies bezeichnet man als Solarkonstante S= 1,367 $kW/m^2$. Hierbei handelt es sich natürlich um einen Mittelwert. Der Abstand Erde-Sonne beträgt $150\cdot10^6$ km, der Sonnenradius ist R = $6,963\cdot 10^8$m,
Wie groß ist die Temperatur der Sonne?
Lösung

Das Modell

In Daisyworld konkurrieren zwei Blumenpopulationen um die noch frei zur Verfügung stehenden Platz, dabei beziehen sie ihre Energie von der Sonne.

Die Sonne liefert die Energie die zum Überleben der zwei Gänseblumenarten notwendig ist. Aber die Sonneneinstrahlung ist nicht konstant. Im laufe ihre Entwicklung nimmt die Sonneneinstrahlung immer mehr zu. Ein Planet, auf dem Leben herrscht, muss mit dieser Herausforderung fertig werden. Die Temperatur des Planeten hängt von der Albedo A des Planeten ab. Diese Albedo setzt sich aus den drei Albedos der weißen Blumen, der schwarzen Blumen und der grauen Fläche zusammen. Die graue Fläche ist die für die Pflanzen noch frei zur Verfügung stehenden Fläche. Die Anteile dieser Flächen können sich verändern, je nachdem wie sich die Blumenpopulationen entwickeln.

In den folgenden Abschnitten werden alle Gleichungen im Detail beschrieben, die notwendig sind, um das System zu beschreiben. Im darauf folgenden Abschnitt wird dann ein Algorithmus angegeben, um die Daisyworld zu simulieren.

Populationsgleichungen

Daisyworld beherbergt schwarze und weiße Gänseblumen, die ihre optimale Wachstumstemperatur bei 22.5 °C haben. Die Blumen unterscheiden sich in ihrer Albedo. Die Albedo der unbewachsenen Fläche des Planeten ist größer als die der schwarzen, aber kleiner als die der weißen Blumen. Die Gleichungen, die die Änderungen der Gänseblumenpopulationen beschreiben lauten:

\begin{equation} \frac{d\alpha_w}{dt} = \alpha_w (\alpha_g \beta (T_w) - \gamma) \end{equation} \begin{equation} \frac{d\alpha_b}{dt} = \alpha_b (\alpha_g \beta (T_b) - \gamma) \end{equation}

Die Bezeichnungen sind der Tabelle 1 zu entnehmen. In Worten lautet Gleichung (3), dass die zeitliche Änderung der weißen Gänseblumen ( $\frac{d\alpha_w}{dt}$ ) proportional zur jeweils vorhandenen Anzahl ( $\alpha_w$) der Blumen ist. $\beta(T_w)$ ist die Wachstumsrate; sie ist temperaturabhängig. Die Abhängigkeit von der noch unbewachsenen Fläche $\alpha_g$ ist ebenfalls berücksichtigt.$\gamma$ ist eine Konstante, die die Sterberate modelliert. Es gilt dann, dass der Faktor $-\gamma\cdot \alpha_w$ die Abnahme der Population beschreibt. Der Faktor $\alpha_w \alpha_g \beta(T_w)$ bestimmt, wie die Population in jedem Zeitschritt ($dt$) zunimmt. (Entsprechendes gilt für Gleichung (4)).

Wir fassen zusammen: Die Änderung der Blumenpopulation werden durch ein Differentialgleichungssystem beschrieben.
Die Blumen sind abhängig von der noch frei zur Verfügung stehenden Fläche $\alpha_g$ und einer temperaturabhängigen Wachstumsrate $\beta$. Die gesamte zur Verfügung stehende Fläche ist auf 1 normiert, d.h. \begin{equation} \alpha_g + \alpha_b + \alpha_w = 1 \end{equation}

Die Wachstumsrate

Die Wachstumsraten der Gänseblumen hängt nur von der Temperatur ab. Wenn es zu kalt ist (ab 5 ° C abwärts), kann kein Wachstum mehr stattfinden. Wenn es auf dem Planeten 40° C und heißer wird, ist es für die Gänseblumen zu lebensfeindlich und es findet abermals kein Wachstum mehr statt. Die optimale Wachstumstemperatur liegt für die Blumen bei 22.5°C. Die Wachstumsrate $\beta(T)$ ist eine Zahl zwischen 0 und 1. 0 bedeutet kein Wachstum und 1 bedeutet optimale Wachstumsbedingungen.

Der Graph der Funktion β(T). Die Wachstumsrate nimmt Werte zwischen 0 (bei 5 °C und 40 °C) und 1 (bei 22.5 °C) an.

Aus dem oben erwähnten Bedingungen folgt für die Wachstumsrate β(T) : \begin{equation} \beta(T) = 1 - \frac{1}{17.5^2} \; (T-22.5)^2\end{equation}

In dieser Form muss T in °Celsius eingegeben werden.

In der Simulation verwenden wir alle Gleichungen in einer Form, die als Eingabe die Temperatur in Kelvin erwarten, als Umrechnungsfaktor verwenden wir 0 °C = 273 Kelvin.

Mit dieser Konvention lautet die Formel dann: \begin{equation} \beta(T) = 1 - \frac{1}{17.5^2} \; (T-295.5)^2\end{equation}

mit $\frac{1}{17.5^2}$ ≈ 0.003265

Leiten Sie ausgehend von den oben erwähnten Bedingungen für die Wachstumsrate die Gleichung für β(T) her.

aber erst einmal selber probieren ;-)

Planetentemperatur und die Rolle der Albedo

In diesem Abschnitt leiten wir die Formel (Gleichung 8) her. Diese Formel ermöglicht die Berechnung der Temperatur des Planeten in Abhängigkeit der Albedo A.

\begin{equation} T^4 = L \frac{S_0}{4\sigma} (1-A) \end{equation}

T ist die Temperatur des Planeten. $S_0$ und $\sigma$ sind Konstanten. L ist eine Zahl, (0.6 ≤ L ≤ 1.65) die es uns ermöglicht, die Änderung der Sonneneinstrahlung zu modellieren. Auch unsere Sonne wird ihre Sonneneinstrahlung in den nächsten Millionen Jahren weiter erhöhen. Für L = 1 beschreiben wir die Sonneneinstrahlung zum jetzigen Zeitpunkt, dementsprechend beschreiben wir mit L > 1 die zukünftige Sonneneinstrahlung. A ist die Albedo des Planeten, auf die weiter unten genauer eingegangen wird.

Die Temperatur des Planeten bleibt dann konstant, wenn die pro Zeit eingestrahlte Energie ($P_{EIN}$) genauso groß ist wie die pro Zeit abgestrahlte Energie ($P_{AUS}$) (Abbildung 3).


\begin{equation}P_{EIN} = P_{AUS}\end{equation}

Die von dem Planeten abgestrahlte Intensität können wir mithilfe des Stefan-Bolzman-Gesetzes angeben, dabei gilt $I_{AUS}= \frac{P_{AUS}}{m^2}$: \begin{equation} I_{AUS} = \sigma\cdot T^4\end{equation}

Hierbei ist T die Temperatur des Planeten.

Im dynamischen Gleichgewicht muss die pro Zeiteinheit (das ist die Leistung) eingestrahlte Energie genau so groß sein, wie die pro Zeiteinheit abgestrahlte Energie, denn wäre z.B. die eingestrahlte Leistung größer, würde sich der Planet immer mehr erwärmen.
Abbildung 3 Wenn die Erde die gesamte von der Sonne abgestrahlte Energie aufnehmen soll, dann müssen alle Strahlen senkrecht auf die Erde fallen. Das ist aber nicht der Fall und somit kann die Erde letztlich nur mit seiner Querschnittfläche ($\pi \cdot R^2$) die von der Sonne abgestrahlte Energie aufnehmen. Ein Teil dieser aufgenommen Energie wird unmittelbar reflektiert (rote Pfeile). Der andere Teil (blaue Pfeile) wird durch Wärmestrahlung abgegeben. Die Wärmestrahlung erfolgt von der gesamten Oberfläche der Erde.

Um einen Ausdruck für die eingestrahlte Leistung ($P_{EIN}$) zu erhalten, wird die Solarkonstante $S_0$ benötigt. Die Solarkonstante taucht jetzt als Produkt mit dem zusätzlichen Parameter L auf. Der Parameter L ermöglicht es, die veränderte Sonneneinstrahlung zu modellieren.

\begin{equation}I_{EIN} =(1-A) L S_0\end{equation}

Hierbei wurde die Albedo A des Planeten eingeführt. Da die Albedo A den Anteil beschreibt der reflektiert wird, beschreibt der Faktor (1-A) denjenigen Anteil, der nicht reflektiert wird und somit auf dem Planeten trifft. Die eingestrahlte Intensität ist also derjenige Teil der Sonneneinstrahlung, der nicht reflektiert wird.

Um die Albedo A des Planeten berechnen zu können, ist es notwendig zu wissen, wie groß der jeweilige Anteil von weißen Blumen, schwarzen Blumen bzw. der noch unbewachsen Fläche ist. Die Albedo des Planeten ist dann das gewichtete Mittel, das sich aus der jeweiligen Flächenanteile multipliziert mit der zugehörigen Albedo berechnet.

\begin{equation} A = \alpha_g \; A_g + \alpha_b \; A_b + \alpha_w \; A_w\end{equation}

Die Albedo A bestimmt damit also indirekt die Planetentemperatur. Da die Albedo A von $\alpha_g$,$\alpha_b$ und $\alpha_w$ abhängig ist, können auf diese Weise die Blumenpopulationen auf die Planetentemperatur einwirken. Wir erinnern an den Zusammenhang zwischen Leistung P und Intensität I: P = I · Fläche.
Aus Abbildung (3) folgt nun:

\begin{equation} P_{EIN} = I_{EIN} \cdot \pi \cdot R^2\end{equation} \begin{equation} P_{AUS} = I_{AUS}\cdot 4 \cdot \pi \cdot R^2 \end{equation}

Durch das Gleichsetzen der obigen Ausdrücke und das Einsetzen der Intensitäten, erhalten wir schließlich Gleichung (15).
(bitte unbedingt nachrechnen):
$$ I_{EIN} = 4\; I_{AUS} \Rightarrow (1-A) L S_0 = 4 \sigma \cdot T^4 $$ \begin{equation} \Rightarrow T^4 = L \frac{S_0}{4\sigma} (1-A) \end{equation}

Die lokalen Temperaturen und der Transportparameter R

Die lokalen Temperaturen $T_b$ , $T_w$ und $T_g$ sind die jeweiligen Temperaturen der weißen und schwarzen Blumen, sowie der unbewachsenen Fläche. Schwarze Blumen haben aufgrund ihrer Albedo eine höhere Temperatur als die weißen Blumen. Da die Wachstumsrate β(T) von der Temperatur abhängt, wirken sich die lokalen Temperaturen auf die Entwickung der Populationen aus. Der Anteil der Blumenpopulationen wird für die Berechnung der Albedo A des Planeten benötigt. Die Albedo A bestimmt wiederum die Temperatur des Planeten. Die verschiedenen (lokalen) Temperaturen können sich mit der Zeit angleichen, da sie durch Wärmetransport miteinander verbunden sind. Wie gut (schnell) dieser Temperaturausgleich ist, wird dabei durch den Parameter R bestimmt. Der Parameter R kann Werte zwischen 0 und 1 annehmen.
Betrachten wir die folgende Gleichung um die Wirkung des Parameters R auf die lokale Temperatur $T_w$ zu untersuchen:

\begin{equation} T_w^4 = R( 1-A_w) \cdot \frac{L\;S_0}{4\;\sigma} + (1-R)(1-A)\cdot\frac{L\; S_0}{4\;\sigma}\end{equation} Für $T_b$ und $T_g$ verläuft die Argumentation völlig analog, deswegen gehen wir in diesem Abschnitt immer von Gl. (16) aus. Die Gleichungen für $T_b$ und $T_g$ lauten:

\begin{equation} T_b^4 = R( 1-A_b) \cdot \frac{L\;S_0}{4\;\sigma} + (1-R)(1-A)\cdot\frac{L\; S_0}{4\;\sigma}\end{equation}

\begin{equation} T_g^4 = R( 1-A_g) \cdot \frac{L\;S_0}{4\;\sigma} + (1-R)(1-A)\cdot\frac{L\; S_0}{4\;\sigma}\end{equation} Um sich Gleichungen die durch einen Parameter gesteuert werden zu veranschaulichen, betrachtet man häufig wie sich die Gleichung für die Randwerte des Parameters verhalten. Für Gleichung 16 bedeutet das, dass der Fall $R=0$ und $R=1$ von besonderem Interesse ist:
Wir untersuchen zuerst den Fall $R = 0$, aus der Gleichung folgt dann:

\begin{equation} T_w^4 = ( 1-A) \cdot \frac{L\;S_0}{4\;\sigma} \end{equation}

Diese Gleichung ist identisch mit Gleichung 8. Es gibt jetzt keinen Unterschied zwischen den Temperaturen mehr: $T_w =T_b = T$. Wenn es keinen Unterschied zwischen den Temperaturen gibt, bedeutet das, dass Temperaturunterschiede sofort ausgeglichen werden. Damit beschreibt $R=0$ den Fall des vollständigen Wärmetransport.
Betrachten wir jetzt den Fall R=1 aus Gleichung 16 folgt nun: \begin{equation} T_w^4 = ( 1-A_w) \cdot \frac{L\;S_0}{4\;\sigma} \end{equation}

Dies ist die entgegengesetzte Situation: Die Temperaturen können sich jetzt nicht mehr ausgleichen. Vielmehr ist es so, als existiere jede Population nur für sich. Dieser Fall beschreibt den Fall ohne Wärmetransport. Damit liegt vollständige Isolation vor.


Wenn der Parameter R nun zwischen 0 ≤ R ≤ 1 variiert wird, haben wir eine Möglichkeit, den Wärmetransport des Planeten vorzugeben. In der Simulation setzen wir $R=0.12$, d.h. die verschiedenen Temperaturen werden relativ schnell einander angeglichen.

Leiten sie die Gleichung 21a aus der Gleichung 16 her. Setzen Sie hierfür die Gl.(8) in die Gl.(16) ein um formen sie um.

...aber erst einmal selber probieren ;-)

Eine Implementierung in Java

Wir fassen alle benötigten Parameter und Gleichungen noch einmal zusammen, bevor wir einen Algorithmus zur Modellierung von Daisyworld angeben.

Die verwendete Parameter und Gleichungen

ParameterBedeutungParameterbezeichnung in JavaTypische Anfangswerte
$A$Albedo des PlanetenA wird berechnet
$A_g$Albedo der unbewachsenen FlächeAg0.5
$A_b$Albedo der schwarzen GänseblumenAb0.25
$A_w$Albedo der weißen GänseblumenAw0.75
$\alpha_g$Anteil der ungewachsenen Fläche (des Untergrundes)ag0.8
$\alpha_b$Anteil (Population) der schwarzen Gänseblumenab0.1
$\alpha_w$Anteil (Population) der weißen Gänseblumenaw0.1
TTemperatur des PlanetenTwird berechnet
$T_w$Temperatur der weißen Gänseblumen (w=white)Tw wird berechnet
$T_b$Temperatur der schwarzen Gänseblumen (b=black)Tb wird berechnet
$T_g$Temperatur der unbewachsenen Fläche (g=ground)Tg wird berechnet
$\beta$(…)Wachstumsfaktor der weißen/schwarzen Blumen beta(…)wird jeweils berechnet
$\gamma$Sterberate der Blumengamma0.3
LKontrollfaktor um die Veränderung der Sonneneinstrahlung ($S_0$) zu modellieren: 0.6 ≤ L ≤ 1.65 L = 1.0 entspricht dem heutigen Wert für die Sonneneinstrahlung.L0.6
RTransportparamer (0 ≤ R ≤ 1), zur Modellierung des Wärmetransportes.R=1 bedeutet vollständige Isolation und R=0 bedeutet vollständiger Wärmetransport.R0.12
$S_0$Solarkonstante (für die Erde gilt $S_0 = 1.367\cdot kW/m^2$ wir verwenden hier einen höheren Wert. Für L=1 gilt dann in der Simulation S=3668 )S3668
$\sigma$Stefan Bolzmann Konstante $\sigma=5.67032\cdot 10^{-8} \frac{W}{K^4\; m^2}$ sigma $5.67032\cdot10^{-8}$
Tabelle 1:

Übersicht über die verwendeten Gleichungen (Temperaturen in Kelvin)

BedeutungGleichungNr.Beschchreibung
Planetarische Albedo $A = \alpha_g\;A_g + \alpha_b\;A_b +\alpha_w\;A_w$(12)Die planetarische Albedo ist in etwa ein Mittelwert aus der Albedo der schwarzen und weißen Gänseblumen, sowie der Albedo $A_g$ des unbewachsenen Untergrundes.
Populationsgrößen
Änderungsrate der weißen Gänseblümchen $\frac{d\alpha_w}{dt} = \alpha_w (\alpha_g\; \beta(T_w) -\gamma)$(3)Diese Gleichung beschreibt die Entwicklung der weißen Gänseblümchen.
Änderungsrate der schwarzen Gänseblümchen $\frac{d\alpha_b}{dt} = \alpha_b (\alpha_g\; \beta(T_w) -\gamma)$(4) Diese Gleichung beschreibt die Entwicklung der schwarzen Gänseblümchen.
Anteil der unbewachsenen Fläche (des Untergrundes)$\alpha_g = 1- \alpha_b -\alpha_w$ (5)
Wachstumsraten
Wachstumsrate (weiße Blumen)$\beta(T_w) =1 - 0.003265\cdot(T_w - 295.5)^2$(7a) Für die weißen Blumen wird $T_w$ verwendet: $\beta(T_w)$
Wachstumsrate (schwarze Blumen)$\beta(T_b) =1 - 0.003265\cdot(T_b - 295.5)^2$(7b) Für die schwarzen Blumen wird $T_b$ verwendet: $\beta(T_b)$
Temperaturen
Temperatur des Planeten$T^4 = L\frac{S_0}{4\sigma}\cdot(1-A)$ (8)
Temperatur der weißen Gänseblumen$T^4_w = R\;L\frac{S_0}{4\sigma}\cdot(A-A_w)+ T^4$ (21a)
Temperatur der schwarzen Gänseblumen$T^4_b = R\;L\frac{S_0}{4\sigma}\cdot(A-A_b)+ T^4$ (21b)
Termperatur der unbewachsenen Fläche$T^4_g = R\;L\frac{S_0}{4\sigma}\cdot(A-A_g)+ T^4$ (21c)
Tabelle 2:

Der Algorithmus

SchritteBeschreibung
Step 1.Startwerte festlegen (Starte mit L=0.6 - dieser Parameter regelt die Sonneneinstrahlung).
Step 2.Berechne die Albedo des Planeten A
Step 3.

Berechne sämtliche Temperaturen:

  • Berechne die Temperatur des Planeten $T$.
  • Berechne die Temperaturen der weißen $T_w$ Blumen.
  • Berechne die Temperaturen der schwarzen $T_w$ Blumen.
  • Berechne die Temperatur der unbewachsenen Fläche $T_g$.
Step 4.

Berechne die Wachstumsrate β der weißen Gänseblumen.
Berechne die Wachstumsrate β der schwarzen Gänseblumen.

Step 5.

Berechne $\alpha_w$ den neuen Anteil der weißen Gänseblumen.
Berechne $\alpha_b$ den neuen Anteil der schwarzen Gänseblumen.
Berechne $\alpha_g$ den neuen Anteil der unbewachsenen Fläche.
Hinweis: Es gilt $\alpha_w + \alpha_b+\alpha_g = 1$.
Falls eine Gänseblumenart < 0.01 ist ($\alpha_b$ < 0.01 oder $\alpha_w$ < 0.01), wird sie zurück auf 0.01 gesetzt.

Step 6.Wiederhole von Step 2., solange bis die Populationsgrößen (Flächenanteile) $\alpha_w$ und $\alpha_b$ konvergieren. Oder eine bestimmte Schrittzahl erreicht wurde.
Step 7.

Zeichne die neuen Anteile $\alpha_w$ und $\alpha_b$ in Abhängigkeit der Sonneneinstrahlung in ein Koordinatensystem.
Alternativ kann auch die Temperatur des Planeten in Abhängigkeit der Sonneneinstrahlung in ein Koordinatensystem eingetragen werden. Zum Vergleich sollte dann auch die Temperatur des Planeten ohne Leben eingetragen werden.

Step 8.Erhöhe L und wiederhole von Step 2. um die nächsten Wertepaare zu erhalten.
Tabelle 3:

Die Änderungen von $\alpha_w$ und $\alpha_b$ werden durch Differentialgleichungen beschrieben, dies muss in der Numerik gesondert berücksichtigt werden. Wir verwenden den Euler-Algorithmus zur numerischen Lösung der Differentialgleichungen mit einer Schrittweite (timeStep) Δt = 0.01 und 1000 Iterationen (Step 6). Step 6 bedeutet ferner, dass das System jeweils genügend Zeit hat, sich in einem stabilen Systemzustand zu entwickeln, bevor L inkrementiert wird. Zu beachten ist auch, dass die Werte $\alpha_w$ und $\alpha_b$ nur einmal am Anfang initialisiert werden und dann über die gesamte Simulation nicht mehr verändert werden. Lovelock und Watson setzten $\alpha_w$ und $\alpha_b$ gleich 0.01, wenn diese Größen einen Wert unter 0.01 annehmen bevor es mit Step 2 weitergeht, das machen wir natürlich auch so. Der Algorithmus ist robust, d.h. initialisiert man z.B. das System neu bevor man jeweils mit Step 2 beginnt, dann ändert sich die Dynamik nicht grundsätzlich, allerdings hängt sie dann wesentlich stärker von den Populationsgrößen ab. Die Modellvariante, bei dem die Populationen nicht immer neu initialisiert werden, ist die realistischere Variante.
Eine Java Bibliothek Daisyworld (inkl. Sourcecode) steht (falls gewünscht) zur Verfügung. Die Abbildungen 1 und 2 wurden mit dieser Java-Klasse und dem freien Tool xmgrace erstellt.

Lesetipps

Das Original von James Lovelock. Diesen Buch bietet einen unterhaltenden und informativen Einstieg in das Thema. …für diejenigen die es ganz genau wissen wollen hier der englische ArtikelExterner Link, indem die Daisyworld genau beschrieben wird.

Dieses interessante Modell hat noch nicht ausgedient. Hier eine Liste mit Folgeartikeln, die aufbauend auf dem einfachen 'Daisywold' Modell viele weitere Aspekte der theoretischen Ökologie behandeln:

weitere Literatur
  • Ackland (2004)
  • Ackland et al (2003)
  • Akagi 2006
  • Adams et al (2003)
  • Andrew (1996)
  • Andrew (1997)
  • Baldocchi et al (2005)
  • Cohen & Rich (2000)
  • Dagg (2002)
  • Dewar (2003)
  • Dewar (2003b)
  • Dewar (2005)
  • Dewar (2005b)
  • Downing (2002)
  • Downing & Zvirinsky (2000)
  • Dyke & Harvey (2005)
  • Dyke & Harvey (2005b)
  • Franck et al (2000)
  • Harding & Lovelock (1996)
  • Harvey (2004)
  • Jascourt & Raymond (1991)
  • Kleidon (2004)
  • Koelsag et al (1999)
  • Lansing et al (1998)
  • Lapenis (2002)
  • Lenton (1998)
  • Lenton (2002)
  • Lenton & Lovelock (2000)
  • Lenton & Lovelock (2001)
  • Lenton & van Oijen (2002)
  • Lenton & Wilkinson (2003)
  • Lovelock (1992)
  • McCormack(2003)
  • Nature (2005)
  • Nevinson et al (1999)
  • Pujol (2002)
  • Pujol et al (2005)
  • Robertson & Robinson (1998)
  • Saunders (1994)
  • Saunders (1998)
  • Saunders et al (2000)
  • Seto & Akagi (2005)
  • Staley (2002)
  • Stocker (1995)
  • Sugimoto (2002)
  • Toniazzo et al (2004)
  • Von Bloh (1997)
  • Von Bloh (1999)
  • Weber (2001)
  • Wilkinson (2003)
  • Williams & Nobel (2005)
  • Zeng (1990)