4.5.1 Grundlagen
4.5.2 Wasser- und Nährstoffbilanzierung
4.5.3 Assimilation und Allokation
4.5.4 Phänologie und Mortalität
4.5.5 Regeneration, Management und Störungen
4.5.6 Parameter, Initialisierung und Triebkräfte
4.5.7 Validierung und Anwendung

Das Modell ist eine Entwicklung des Potsdam-Instituts für Klimafolgenforschung. Eine ausführliche Modellbeschreibung findet sich auf www.pik-potsdam.de/research/climate-impacts-and-vulnerabilities/models/4c.

Abb. 4‑5: Modellschema 4C

Das Waldwachstumsmodell 4C (FORESEE – FORESt Ecosystems in a changing Environment) beschreibt die Sukzessionsdynamik von Waldbeständen, die entweder auf Basis von simulierter natürlicher Regeneration oder ausgehend von einem durch eine Forstinventur definierten Anfangszustand aufwachsen (Schaber et al., 1999, Suckow et al., 2001).

Die Baumindividuen einer Art, die gleiches Alter und identische Baum­dimensionen besitzen, sind in Kohorten zusammengefasst. Produktion und Wachstum werden für jede Kohorte berechnet. Die Konkurrenz der Kohorten um Licht, Wasser und Nährstoffe beeinflusst ihr Wachstum, ihre Mortalität und die Ver­jüngung im Bestand. Die Positionen der Individuen der Kohorten im Bestand sind nicht bekannt; es wird angenommen, dass die Bäume gleichmäßig auf der Bestandesfläche verteilt sind. Es werden keine Unterschiede in den Wachstums­bedingungen einzelner Individuen derselben Kohorte berücksichtigt. Der Kronenraum des Bestandes ist in 0.5 m hohe Schichten eingeteilt. Das Blattwerk ist auf die Kronen­schichten zwischen Kronenansatzhöhe und Baumhöhe verteilt. Abb. 4‑5 beschreibt das prinzipielle Zusammenwirken der Teilmodelle zur Vege­tations- und Bodendynamik in 4C.

Die Berechnung der Flüsse und die Änderung der Zustandsvariablen erfolgt prozess­abhängig mit unterschiedlichen zeitlichen Schrittweiten:

• Wasser- und Nährstoffflüsse sowie Phänologie: im Tagestakt,
• Photosynthetische Produktion: Schrittweiten zwischen wahlweise einem Tag bis zu einer Woche,
• Allokation, Wachstum, Mortalität und Regeneration: im Jahrestakt.

Das Wasser- und Nährstoffangebot (Stickstoff, Kohlenstoff) wird in Abhängigkeit von Boden, Bestand und Wetter bilanziert (Einzelmodell 4C: Grote et al., 1999; 4C im Modellverbund ArcEGMO-PSCN: Kap. 5 und 7). 4C berechnet zunächst für jede Kohorte einen aktuellen Transpirationsentzug in Abhängigkeit von der durch die Interzeptionsverdunstung reduzierten potenziellen Evapotranspiration und der aktuellen Bodenfeuchte der durchwurzelten Schichten, der als Wasserentzug pro Schicht an das Bodenwassermodell übergeben wird.

Durch den täglichen Entzug von Wasser und Nährstoffen einerseits sowie die jährliche Bilanzierung des Streufalls und die Zufuhr zum Bodenkompartiment andererseits wird der Nährstoffkreislauf im System Pflanze – Boden geschlossen. Zusätzlichen Eintrag erhält das System durch Deposition, Verluste treten durch Auswaschung auf.

Die photosynthetische Nettoassimilationsleistung wird nach einem Ansatz von Haxeltine & Prentice (1996) als Funktion von absorbierter photosynthetisch aktiver Strahlung, Lufttemperatur sowie Bodenwasser- und Nährstoffverfügbarkeit berechnet. Die Produktion der einzelnen Kohorten ist artenspezifisch und abhängig vom Anteil der Kohorte an der vom Bestand absorbierten Strahlung. Aus der jähr­lichen Bruttoassimilationsleistung wird nach Abzug der Respiration das Wachstum der Kompartimente Feinwurzeln, Stamm, Blattwerk, sowie Äste und Grobwurzeln berechnet. Die Allokationskoeffizienten für die einzelnen Kompartimente werden dabei so bestimmt, dass sie einem vorgeschriebenen Verhältnis des Querschnitts der leitenden Gewebe im Holz zu den zu versorgenden Blattmassen und einer aus­gewogenen Leistung des Aufnahmevermögens der Feinwurzeln und der Produk­tionskapazität der Blätter entsprechen (functional balance und pipe model theory, siehe z.B. Mäkelä, 1986). Zusätzlich gehen in die Bestimmung der Allokations­koeffizienten eine Höhenwachstums­funktion ein, die vom relativen Lichtgenuss und dem Blattmassenzuwachs abhängig ist, sowie die zu ersetzenden Streuproduktions­flüsse. Die Kronenansatzhöhe wächst, wenn die Nettoproduktion der untersten Kronenschicht negativ wird. Damit sind die Baumdimensionen Höhe, Brusthöhen­durchmesser und Kronenansatzhöhe im Modell berechenbar.

Das Modell enthält ein Modul zur Simulation des Blattaustriebstages der Buche, das mit dem phänologischen Datensatz des Deutschen Wetterdienstes parametrisiert wurde (Schaber & Badeck, 2003). Bei der Berechnung des Blattaustriebs wird davon ausgegangen, dass die Regulation der ontogenetischen Entwicklung durch die Interaktion von wachstumsfördernden (Promotoren) und wachstumshemmenden (Inhibitoren) Verbindungen (Phytohormone) erfolgt und die Balance zwischen Promotoren und Inhibitoren den physiologischen Entwicklungszustand der Pflanze und ihre Reaktion auf externe Bedingungen bestimmt (Schaber, 2002). Wesentliche Einflussfaktoren sind dabei die Lufttemperatur und Tageslänge.

Der Blattwurf wird anhand eines festen durchschnittlichen Datums bestimmt. Für die Laubbaumarten wurden die erforderlichen Parameter an Daten des DWD angepasst. Für die immergrünen Baumarten werden gegenwärtig im Modell noch keine phänologischen Stadien unterschieden.

Die Mortalität der Individuen innerhalb der Kohorten wird aus der Kohlenstoffbilanz bestimmt (stressbedingte Mortalität) oder nach einer vorgegebenen altersbedingten Mortalitätsrate. Beide Ansätze können auch kombiniert werden. Stressbedingte Mortalität tritt auf, wenn die Kohlenstoffbilanz über einen bestimmten Zeitraum negativ wird und daher z.B. die Blattmasse der Individuen einer Kohorte eines Jahres geringer ist als die des Vorjahres. Sie führt damit zur Verringerung der Individuenzahl bzw. zum vollständigen Absterben der Kohorte. Die intrinsische artenspezifische Mortalität wird über ein maximales Alter der Baumart definiert (Botkin & Nisbet, 1992).

Auf Bestandesebene können Verjüngungs- und Bewirtschaftungsmaßnahmen simuliert werden. Die Verjüngung beschreibt das Aufwachsen von gepflanzten Setzlingen oder von Sämlingen, die jährlich in Abhängigkeit von Umweltbedingungen im Bestand aufkeimen können.

Das Modell erlaubt die Simulation einer Vielzahl von Bewirtschaftungsvarianten, die durch die Kombination verschiedener implementierter Methoden möglich sind. Ein Bewirtschaftungsplan steuert die verschiedenen Eingriffe. In Abhängigkeit von der Bestandeshöhe werden Jungwuchspflege, Läuterung und Jungbestandespflege simuliert. Für die anschließenden Bestandeseingriffe kann Hoch- oder Nieder­durchforstung gewählt werden. Die Durchforstung wird durch Durchforstungsstärke und Durchforstungsintervall beschrieben. Zur Ernte des Bestandes, gesteuert über die Umtriebszeit, ermöglicht das Modell den Kahlschlag mit anschließender Neu­anpflanzung oder natürlicher Verjüngung oder einen Schirmschlag, ebenfalls kombiniert mit einer Unterpflanzung. Für die Pflanzung kann eine Artenspezifikation und Pflanzdichte vorgegeben werden. Die Generierung von Mischbeständen ist möglich. Der Bewirtschaftungsplan erlaubt damit auch langfristige Simulationen über mehrere Umtriebszeiten.

Störungen wie Windwurf oder Borkenkäferbefall im Bestand können im Modell durch Verarbeitung einer Störungszeitreihe simuliert werden. Diese Störungen wirken insbesondere durch eine zusätzliche Mortalität in den Baumkohorten auf die Bestandesdynamik.

Die Störungsdynamik durch Borkenkäferbefall wird durch eine Zeitreihe der Störungsereignisse mit Information über befallene Fläche (Koordinaten), Störungstyp, Störungsjahr, Anteil der befallenen Fläche an der Gesamtfläche des Bestandes und Typ des Bestandes (bewirtschaftet oder unbewirtschaftet) beschrieben. Im Fall von Borkenkäferbefall wird die Mortalität der Fichten, die älter als 50 sind, angenommen. Bei den bewirtschafteten Beständen werden alle befallenen und damit gestorbenen Fichten vollständig geerntet, Wurzeln und Benadlung verbleiben im Ökosystem. Bei den unbewirtschafteten Beständen werden diese Bäume komplett in die Zersetzungspools des Bodens überführt, für die Stämme, wird eine verzögerte Über­führung in den Zersetzungspool für Stammholz angenommen, da die toten Bäume noch einige Jahre im Bestand stehen. Sowohl im bewirtschafteten als auch im unbewirtschafteten Fall, wird bei einem Anteil der Störungsfläche, akkumuliert über mehrere aneinanderfolgende Störungsjahre, von über 30% der Bestandesfläche, eine natürliche Verjüngung der vorkommenden Baumarten simuliert.

Die Artenzusammensetzung und Bestandesstruktur können auf der Basis von Inventurdaten initialisiert werden. Die Kohorten werden dann entsprechend der gemessenen Verteilungen von Brusthöhendurchmesser sowie, falls vorhanden, Baumhöhe und Kronenansatzhöhe definiert.

Als treibende klimatische Größen werden Tagesmittel bzw. -summen der Tempe­ratur, des Niederschlags, der Luftfeuchte und der Strahlung benötigt, die wahlweise aus Messreihen eingelesen oder aus Monatswerten mit Hilfe eines Wettergenerators hergeleitet werden können. Außerdem wird eine Zeitreihe der atmosphärischen CO₂-Konzentration und Stickstoff-Deposition vorgegeben.

Das Modell liegt in einer ersten, hier umrissenen Ausbaustufe vor und wird an langen Zeitreihen von Dauerbeobachtungsflächen hinsichtlich der Bestandesdynamik getestet (Schaber 2002, Schaber et al. 1999; Mäkelä et al., 2000). Es ist gegenwärtig für die Baum­arten Buche (Fagus sylvatica L.), Fichte (Picea abies L.), Kiefer (Pinus sylvestris L.), Eiche (Quercus robur L.) und Birke (Betula pendula Roth) parametrisiert. Einzelne Komponenten des Waldwachstumsmodells 4C wurden für mehrere Baumarten an verschiedenen Standorten (v.a. an Level II Standorten in Brandenburg, Sachsen und Thüringen) validiert.

Im Rahmen des BMBF-Projektes „Wälder und Forstwirtschaft im Globalen Wandel: Strategien für eine integrierte Wirkungsanalyse und –bewertung“ (LK9528-9533) und des EU-Projektes SilviStrat“: Silvicultural Response Strategies to Climatic Change in Management of European Forests” (EVK2-2000-00723) wurde das Modell für Analysen eingesetzt (Lasch et al., 2002; Lasch et al., 2005).