Zuwachsmodell


  1. Überblick
    1. Name
    2. Modellentwickler
    3. Zielsetzung/Zweck
    4. Modelltyp
    5. Zustandsvariablen
    6. Skalenebenen
    7. Prozessüberblick und Ablaufplan
    8. Stochastizität
  2. Details
    1. Input
  3. Literatur


Überblick

Name des Modells:

TreeGrOSS, Version 7.5.7

Modellentwickler:

Prof. Dr. J. Nagel, Dr. H. Duda, J. Hansen, Dr. S. Hentschel, Dr. M. Schmidt, Dr. M. Albert

Zielsetzung/Zweck:

Das Waldwachstumsmodell prognostiziert das Einzelbaumwachstum mit den Modelleinstellungen für Nordwestdeutschland (Nagel, 2009, Döbbeler et al., 2002).

Modelltyp:

statistisches, distanzunabhängiges Einzelbaumwachtumsmodell

Zustandsvariablen:
Name Einheit Beschreibung
ig cm2 Grundflächenzuwachs des Einzelbaumes
ih m Höhenzuwachs des Einzelbaumes
Skaleneben:

Einzelbaum, 5-jähige Prognoseperioden

Prozessüberblick und Ablaufplan:

Das Waldwachstumsmodell BWINPro wurde für 5-jährige Zeitschritte parametrisiert. Dies hängt im Wesentlichen damit zusammen, dass die meisten Versuchsflächen in einem 5-jährigen Turnus aufgenommen wurden. Theoretisch lassen sich aber mit der Software auch kleinere Zeitintervalle simulieren. Ein Zeitschritt wird mit der Methode grow() der Klasse Stand ausgelöst und es werden von der Software die folgende Schritte durchlaufen:

  1. Überprüfung der altersbedingten Mortalität
  2. Überprüfung der konkurrenzbedingten Mortalität
  3. Schätzung von Höhen- und Grundflächenzuwachs
  4. Erhöhung des Alters um 5 Jahre
  5. Bestimmung des BHD und der Höhe nach 5 Jahren
  6. Anpassung von Kronenansatz und -breite an die neuen Dimensionen
  7. Update der ertragskundlichen Werte durch die Methode descspecies()
Stochastizität:

Zufallseffekte sind in der Durchmesserzuwachs- und Höhenzuwachsschätzung enthalten.

Details

Die Simulationssoftware TreeGrOSS wurde in erster Linie für das Waldwachstumsmodell BWINPro erstellt (Nagel et al., 2002, Nagel, 2005, Nagel et al., 2006). Bei dem Wachstumsmodell handelt es sich um ein statistisches Modell, bei dem jeder einzelne Baum eines Bestandes in seiner Entwicklung beschrieben wird. Dieser sogenannte Einzelbaumansatz ermöglicht es, nahezu je Bestandesstruktur und –zusammensetzung zu simulieren. In BWINPro wird das Baumwachstum stark abstrahiert und ist auf den BHD- und Höhenzuwachs reduziert. Die Funktionen zur Zuwachsschätzung wurden mit Hilfe des Versuchsflächenmatrials baumartenweise parametrisiert. Als Zuwachs beeinflussende Größen werden das Alter, die Kronenmantelfläche, der Kronenkonkurrenzindex und dessen Veränderung bei Durchforstungen verwendet. Über den spezifischen Standraumbedarf der Baumarten, der sich hauptsächlich aus der Kronengröße ergibt, werden die Mischbestandseffekte im Modell realisiert. Die Kronengröße wird im Modell aus der Höhe, dem Kronenansatz und der Kronenbreite für einen unterstellten Paraboliden berechnet, wobei die für die Kronenmantelfläche und den Kronenkonkurrenzindex notwendigen Kroneninformationen über statische Funktionen aus dem BHD, der Höhe und der Bestandesoberhöhe geschätzt werden.

Der Simulator lässt sich theoretisch nur einsetzen, wenn ein kompletter Datensatz für einen Bestand zur Verfügung steht. Im praktischen Forstbetrieb, der Forsteinrichtung und auch bei Inventuren ist es aus Zeit- und Kostengründen meist nicht möglich, eine derart aufwendige Datenerhebung durchzuführen. Aus diesem Grund ist die Simulationssoftware mit einigen Datenergänzungsroutinen und einer Routine zur Erzeugung von Durchmesserverteilungen ausgestattet. Mit dem Aufruf der Methode missingData() der Klasse Stand werden fehlende Informationen zu den Einzelbäumen ergänzt.

Höhenzuwachsschätzung

In den meisten Einzelbaummodellen ist die Schätzung des Höhenzuwachses mit der größten Modellunsicherheit behaftet. Selten haben Höhenzuwachsmodelle Bestimmtheitsmasse, die größer als 0,4 betragen. Dies hängt zu einem großen Teil damit zusammen, dass die Baumhöhen schwierig zu messen und daher mit einem nicht unerheblichem Fehler belastet sind. Gerade in Laubholzbeständen ist der Höhenzuwachs älterer Bäume selbst bei fünfjährigen Messintervallen zum Teil geringer als der Messfehler. So haben zum Beispiel Buchen der II. Ertragsklasse nach Schober in einem Alter von 120 Jahren einen durchschnittlichen Höhenzuwachs von 0,6m bei einer Mittelhöhe von 30,5m. Der Fehler in der Höhenmessung wird bei Kramer und Akca (1982) mit 1% angegeben. Nach den Erfahrungen der Nordwestdeutschen Forstlichen Versuchsanstalt dürfte er bei Laubholz aber eher bis zu 3% betragen. Das bedeutet, dass der Höhenmessfehler bereits der einmaligen Messung mit 0,92m erheblich über dem fünfjährigen Zuwachs liegt. Um den Höhenwachstumsmodellen eine gewisse Robustheit zu verleihen, werden daher häufig Höhenwachstumsentwicklungen verwendet, die aus Bestandeswerten hergeleitet werden.

Der Höhenzuwachs wird in vielen statistischen Waldwachstumssimulatoren über zwei Teilmodelle berechnet. Das erste Teilmodell schätzt den potentiellen Höhenzuwachs. Der potentielle Höhenzuwachs kann entweder aus Beobachtungswerten, um die eine Umhüllende gelegt wird hergeleitet werden (Pretzsch 2001; Guericke 2001) oder er kann aus der Oberhöhenentwicklung von Bonitätsrahmen von Ertragstafeln oder Bestandesmodellen berechnet werden.

Der individuelle Höhenzuwachs eines Baumes wird in den meisten statistischen Modellen geschätzt, in dem der potentielle Höhenzuwachs mit einem Faktor modifiziert wird (Hasenauer 2006).

Durchmesserzuwachsschätzung

Im Modell BWINPro wird der Durchmesserzuwachs aus der Schätzung des Grundflächenzuwachses errechnet. Der Grundflächenzuwachses wird über die Variablen Alter, Kronenmantelfläche, Kronenkonkurrenzindex und die Veränderung des Kronenkonkurrenzindexes bestimmt. Die Kronenmantelfläche wird aus dem Kronenansatz und der Kronenbreite nach der Formel eines quadratischen Paraboliden hergeleitet.

Die Kronenveränderung erfolgt in BWINPro über statischen Funktionen für den Kronenansatz und die Kronenbreite. Für den Kronenansatz gilt die Nebenbedingung, dass der Kronenansatz nicht kleiner werden kann.

Bestandesbehandlung

Die Wuchsleistung der Bäume und auch des gesamten Bestandes kann durch forstliche Eingriffe erhebliche beeinflusst werden. Aus diesem Grund muss ein Waldwachstumssimulator auch Eingriffmöglichkeiten verfügen. Die Eingriffe erfolgen hier automatisiert. Baumartenspezifisch wird ab einer erreichten Oberhöhe und bei entsprechender Bestandesgrundfläche eine mäßige Durchforstung durchgeführt. Die Intensität des Eingriffs orientiert sich an der Absenkung des natürlichen Bestockungsgrades.

Input:

Baumart Eiche (Quercus spec.)

Variable Funktion
Durchmessergenerierung [cm]

t.d = (-1.937+1.082*sp.dg)*(((6.9/(-1.937+1.082*sp.dg))^ (4.669+0.366*sp.dg-0.234*dmax))-ln(1.0-random))^ (1.0/(4.669+0.366*sp.dg-0.234*dmax))
(NAGEL u. BIGING 1995)

Einheitshöhenkurve [m]

t.h = 1.3+(sp.hg-1.3)*exp(0.14657227*(1.0-(sp.dg/t.d)))*exp(3.78686023* ((1.0/sp.dg)-(1.0/t.d)))
(NAGEL 1999)

Höhenvariabilität [m]

t.hv = -0.1944676+0.3535610*ln(sp.dg)
(ALBERT 2000)

Volumenfunktion [m³]

t.v = 3.141592*t.h*(t.d/200)^2*(0.4786-(1.011176/t.d)+(2.10428/t.h)- (203.1997/(t.d*t.h*t.h)))
(BERGEL 1974)

Kronenbreite [m]

t.cw = (2.6618+0.1152*t.d)*(1.0-exp(-exp(ln(t.d/8.3381)*1.4083)))
(DÖBBELER ET. AL. 2001)

Kronenansatz [m]

t.cb = t.h*(1.0-exp(-abs((-0.5268+0.2287*t.h/t.d-0.00453*t.d+0.4712*ln(sp.h100)))))
(DÖBBELER ET. AL. 2001)

Höhe entsprechend der Mittelhöhenbonität des Standort-Leistungs-Modells hg100 [m]

Höhe = 1.2164*t.hg100*(1.0-exp(-0.0194*t.age))^1.1344
(DÖBBELER ET. AL. 2001)

potenzieller Höhenzuwachs [m]

ihpot = ((1.2164*t.hg100*(1.0-exp(-0.0194*(t.age+5.0)))^1.1344)
-(1.2164*t.hg100*(1.0-exp(-0.0194*t.age))^1.1344))
(DÖBBELER ET. AL. 2001)

Höhenzuwachs [m]

t.hinc = t.h*((((1.2164*t.hg100*(1.0-exp(-0.0194*(t.age+5.0)))^1.1344)
-(1.2164*t.hg100*(1.0-exp(-0.0194*t.age))^1.1344))/sp.h100)+ (0.01676*(t.hinc^1.3349)))
(DÖBBELER ET. AL. 2001)

Höhenzuwachsstreuung [m]

herror = 0.082

Grundflächenzuwachs [cm²]

t.dinc = exp(-6.5350+1.3260*ln((3.141593*(t.cw/2.0)/(6.0*(t.h-t.cb)^2.0))* (((4.0*(t.h-t.cb)^2.0+(t.cw/2.0)^2.0)^1.5)
-(t.cw/2.0)^3.0))-0.8437*ln(t.age)-0.9373*t.c66xy+0.1239* t.c66cxy-0.1263*ln(5.0))
(DÖBBELER ET. AL. 2001)

Durchmesserzuwachsstreuung [cm]

derror = 0.617
(DÖBBELER ET. AL. 2001);

Maximale Dichte der Grundfläche [m²/ha]

MaxDichte = 0.0001*3.141592/(16*0.000002807*0.5814*(t.h^(0.9082-1.1830)))
(DÖBBELER 2004)

Maximales Alter [Jahre] MaxAlter = 600
Zielstärkendurchmesser [cm] 80.0
Höhe der 1. Durchforstung [m] 14.0
Mäßige Durchforstung [m; NB°] 14.0;0.75;18.0;18.0;0.80;24.0;24.0;0.85;100.0

Baumart Buche (Fagus silvatica)

Variable Funktion
Durchmessergenerierung [cm]

t.d = (-4.282+1.132*sp.dg)*(((6.9/(-4.282+1.132*sp.dg))^ (4.518+0.317*sp.dg-0.200*dmax))-ln(1.0-random))^ (1.0/(4.518+0.317*sp.dg-0.200*dmax))
(NAGEL u. BIGING 1995)

Einheitshöhenkurve [m]

t.h = 1.3+(sp.hg-1.3)*exp(0.20213328*(1.0-(sp.dg/t.d)))*exp(5.64023296* ((1.0/sp.dg)-(1.0/t.d)))
(NAGEL 1999)

Höhenvariabilität [m]

t.hv = 1.1217150+0.2203473*ln(sp.BHD_STD)
(ALBERT 2000)

Volumenfunktion [m³]

t.v = 3.141592*t.h*(t.d/200)^2*(0.4039+0.0017335*t.h+1.1267/t.h-118.188/ (t.d*t.d*t.d)+0.0000042*t.d*t.d)
(BERGEL 1973)

Kronenbreite [m]

t.cw = (2.0837+0.15*t.d)*(1.0-exp(-exp(ln(t.d/5.7292)*1.3341)))
(DÖBBELER ET. AL. 2001)

Kronenansatz [m]

t.cb = t.h*(1.0-exp(-abs((0.25704+0.11819*t.h/t.d-0.002065*t.d+0.13831* ln(sp.h100)))))
(DÖBBELER ET. AL. 2001)

Höhe entsprechend der Mittelhöhenbonität des Standort-Leistungs-Modells hg100 [m]

Höhe = -75.65900+23.19200*ln(t.age)-1.46800*(ln(t.age)^2)+0.0*t.hg100+0.21520* t.hg100*ln(t.age)
(NAGEL 1999

potenzieller Höhenzuwachs [m]

ihpot = ((-75.65900+23.19200*ln(t.age+5)-1.46800*((ln(t.age+5))^2.0)+0.21520* t.hg100*(ln(t.age+5)))-(-75.65900+23.19200*ln(t.age)-1.46800* ((ln(t.age))^2.0)+0.21520* t.hg100*(ln(t.age))))
(NAGEL 1999

Höhenzuwachs [m]

t.hinc = t.h*((((-75.65900+23.19200*ln(t.age+5)-1.46800*((ln(t.age+5))^2.0)+ 0.21520*t.hg100*(ln(t.age+5)))-(-75.65900+23.19200*ln(t.age)-1.46800* ((ln(t.age))^2.0)+0.21520*t.hg100*(ln(t.age))))/sp.h100)+ (0.00159*(t.hinc^1.9086))
(DÖBBELER ET. AL. 2001)

Höhenzuwachsstreuung [m]

herror = 0.082

Grundflächenzuwachs [cm²]

t.dinc = exp(-7.393+1.375*ln(3.14159265359*(t.cw/2.0)/(6.0*(t.h-t.cb)^2)* ((4.0*(t.h-t.cb)^2+(t.cw/2)^2)^1.5-(t.cw/2)^3))-0.791*ln(t.age) -0.793*t.c66xy+0.809*t.c66cxy-0.0*ln(5.0)))
(DÖBBELER ET. AL. 2001)

Durchmesserzuwachsstreuung [cm]

derror = 0.762

Maximale Dichte der Grundfläche [m²/ha]

MaxDichte = 0.0001*3.141592/(16*0.00000010829*8.3652*(t.h^(1.5374-1.7365)))
(DÖBBELER 2004) 

Maximales Alter [Jahre]

MaxAlter = 300

Zielstärkendurchmesser [cm]

60.0

Höhe der 1. Durchforstung [m]

12.0

Mäßige Durchforstung [m;NB°]

12.0;0.7;22.0;22.0;0.65;28.0;28.0;0.75;100.0


Baumart Fichte (Picea abies)

Variable Funktion
Durchmessergenerierung [cm]

t.d = (-2.492+1.104*sp.dg)*(((6.9/(-2.492+1.104*sp.dg))^(3.418+0.353*sp.dg- 0.192*dmax))-ln(1.0-random))^(1.0/(3.418+0.353*sp.dg-0.192*dmax))
(NAGEL u. BIGING 1995)

Einheitshöhenkurve [m]

t.h = 1.3+(sp.hg-1.3)*exp(0.18290951*(1.0-(sp.dg/t.d)))* exp(5.68789430*((1.0/sp.dg)-(1.0/t.d)))
(NAGEL 1999)

Höhenvariabilität [m]

t.hv = 0.1441427+0.5552640*ln(sp.BHD_STD)
(ALBERT 2000)

Volumenfunktion [m³]

t.v = 3.141592*t.h*(t.d/200)^2*(0.04016-27.56211/(t.d*t.d)+1.36195/ln(t.d) +0.057654*t.h/t.d)
(BERGEL 1987)

Kronenbreite [m]

t.cw = (1.2644+0.1072*t.d)*(1.0-exp(-exp(ln(t.d/0.000001)*1.0)))
(NAGEL 1999)

Kronenansatz [m]

t.cb = t.h*(1.0-exp(-abs((2.0417-0.3335*t.h/t.d+0.00906*t.d-0.9004* ln(sp.h100)))))
(DÖBBELER ET AL 2002)

Höhe entsprechend der Mittelhöhenbonität des Standort-Leistungs-Modells hg100 [m]

Höhe = -49.87200+7.33090*ln(t.age)+0.77338*(ln(t.age)^2)+0.52684*t.hg100+ 0.10542*t.hg100*ln(t.age)
(NAGEL 1999)

potenzieller Höhenzuwachs [m]

ihpot = ((1.2164*t.hg100*(1.0-exp(-0.0194*(t.age+5.0)))^1.1344) -(1.2164*t.hg100*(1.0-exp(-0.0194*t.age))^1.1344))
(NAGEL 1999)

Höhenzuwachs [m]

t.hinc = t.h*((((-49.87200+7.33090*ln(t.age+5)+0.77338*((ln(t.age+5))^2.0)+ 0.52684*t.hg100+0.10542*t.hg100*(ln(t.age+5)))-(-49.87200+7.33090* ln(t.age)+0.77338*((ln(t.age))^2.0)+0.52684*t.hg100+0.10542* t.hg100*(ln(t.age))))/sp.h100)+(0.00271*(t.hinc^2.1725)))
(DÖBBELER ET AL 2002)

Höhenzuwachsstreuung [m]

herror = 0.082

Grundflächenzuwachs [cm²]

t.dinc = exp(-6.2018+1.2984*ln(3.14159265359*(t.cw/2.0)/(6.0*(t.h-t.cb)^2)* ((4.0*(t.h-t.cb)^2+(t.cw/2)^2)^1.5-(t.cw/2)^3))-0.9366*ln(t.age)- 1.2835*t.c66xy+0.2962*t.c66cxy+0.2926*1.6094)
(DÖBBELER ET AL 2002)

Durchmesserzuwachsstreuung [cm]

derror = 0.638

Maximale Dichte der Grundfläche [m²/ha]

MaxDichte = 0.0001*3.141592/(16*0.0000012874*1.2842*(t.h^(0.7148-1.1914)))
(DÖBBELER 2004)

Maximales Alter [Jahre]

MaxAlter = 240

Zielstärkendurchmesser [cm]

45.0

Höhe der 1. Durchforstung [m]

12.0

Mäßige Durchforstung [m;NB°]

12.0;0.7;20.0;20.0;0.75;26.0;26.0;0.8;100.0


Baumart Douglasie (Pseudotsuga menziesii)

Variable Funktion
Durchmessergenerierung [cm]

t.d = (-0.621+1.060*sp.dg)*(((6.9/(-0.621+1.060*sp.dg))^(4.380+0.236*sp.dg- 0.141*dmax))-ln(1.0-random))^(1.0/(4.380+0.236*sp.dg-0.141*dmax))
(NAGEL u. BIGING 1995)

Einheitshöhenkurve [m]

t.h = 1.3+(sp.hg-1.3)*exp(0.19965100*(1.0-(sp.dg/t.d)))*exp(4.63277655* ((1.0/sp.dg)-(1.0/t.d)))
(NAGEL 1999)

Höhenvariabilität [m]

t.hv = 0.2071047+0.5843520*ln(sp.BHD_STD)
(ALBERT 2000)

Volumenfunktion [m³]

t.v = 3.141592*t.h*(t.d/200)^2*(-200.31914/(t.h*t.d*t.d)+0.8734/t.d-0.0052* ln(t.d)*ln(t.d)+7.3594/(t.h*t.d)+0.46155)
(BERGEL 1987)

Kronenbreite [m]

t.cw = (2.919+0.0939*t.d)*(1.0-exp(-exp(ln(t.d/10.0161)*1.362)))
(DÖBBELER ET AL 2002)

Kronenansatz [m]

t.cb = t.h*(1.0-exp(-abs((-1.8796+0.34056*t.h/t.d-0.00610*t.d+0.8262* ln(sp.h100)))))
(DÖBBELER ET AL 2002)

Höhe entsprechend der Mittelhöhenbonität des Standort-Leistungs-Modells hg100 [m]

Höhe = -47.09070+11.4322*ln(t.age)-0.0*(ln(t.age)^2)-0.0*t.hg100+0.20063* t.hg100*ln(t.age)
(NAGEL 1999)

potenzieller Höhenzuwachs [m]

ihpot = ((-47.09070+11.4322*ln(t.age+5)-0.0*((ln(t.age+5))^2.0)-0.0*t.hg100+ 0.20063*t.hg100*(ln(t.age+5)))-(-47.09070+11.4322*ln(t.age)-0.0* ((ln(t.age))^2.0)-0.0*t.hg100+0.20063* t.hg100*(ln(t.age))))
(NAGEL 1999)

Höhenzuwachs [m]

t.hinc = t.h*((((-47.09070+11.4322*ln(t.age+5)-0.0*((ln(t.age+5))^2.0) -0.0*t.hg100+0.20063*t.hg100*(ln(t.age+5)))-(-47.09070+11.4322*ln(t.age) -0.0*((ln(t.age))^2.0)-0.0*t.hg100+0.20063*t.hg100*(ln(t.age))))/sp.h100) +(0.00159*(t.hinc^2.5255)))
(DÖBBELER ET AL 2002)

Höhenzuwachsstreuung [m]

herror = 0.082

Grundflächenzuwachs [cm²]

t.dinc = exp(-7.9766+1.5135*ln(3.14159265359*(t.cw/2.0)/(6.0*(t.h-t.cb)^2)* ((4.0*(t.h-t.cb)^2 +(t.cw/2)^2)^1.5 (t.cw/2)^3))-1.0009*ln(t.age)-0.4481*t.c66xy+ 0.5099*t.c66cxy+0.3038*1.6094)
(DÖBBELER ET AL 2002)

Durchmesserzuwachsstreuung [cm]

derror = 0.725

Maximale Dichte der Grundfläche [m²/ha]

MaxDichte = 0.0001*3.141592/(16.0*0.00000089306*2.4088*(t.h^(0.7726-1.3555)))
(DÖBBELER 2004)

Maximales Alter [Jahre]

MaxAlter = 280

Zielstärkendurchmesser [cm]

65.0

Höhe der 1. Durchforstung [m]

14.0

Mäßige Durchforstung [m;NB°]

14.0;0.65;20.0;20.0;0.70;26.0;26.0;0.7;100.0


Baumart Kiefer (Pinus sylvestris)

Variable Funktion
Durchmessergenerierung [cm]

t.d = (-0.047+1.047*sp.dg)*(((6.9/(-0.047+1.047*sp.dg))^ (3.640+0.332*sp.dg-0.180*dmax))-ln(1.0-random))^ (1.0/(3.640+0.332*sp.dg-0.180*dmax))
(NAGEL u. BIGING 1995)

Einheitshöhenkurve [m]

t.h = 1.3+(sp.hg-1.3)*exp(0.25963741*(1.0-(sp.dg/t.d)))*exp(1.30645374* ((1.0/sp.dg)-(1.0/t.d)))
(NAGEL 1999)

Höhenvariabilität [m]

t.hv = -1.8315300+0.9701583*ln(sp.dg)
(ALBERT 2000)

Volumenfunktion [m³]

t.v = 3.141592*t.h*(t.d/200)^2*(0.40804-318.3342/(t.h*t.d*t.d)+36.90522/ (t.h*t.d)-4.05292/(t.d*t.d))
(BERGEL 1987)

Kronenbreite [m]

t.cw = (1.2783+0.11388*t.d)*(1.0-exp(-exp(ln(t.d/8.705220)*1.33944)))
(DÖBBELER ET AL 2002)

Kronenansatz [m]

t.cb = t.h*(1.0-exp(-abs((1.2085-0.2392*t.h/t.d+0.00742*t.d-0.7897* ln(sp.h100)))))
(DÖBBELER ET AL 2002)

Höhe entsprechend der Mittelhöhenbonität des Standort-Leistungs-Modells hg100 [m]

Höhe = -31.67480+11.64500*ln(t.age)-1.04989*(ln(t.age)^2)-0.43221*t.hg100+ 0.31253*t.hg100*ln(t.age)
(NAGEL 1999)

potenzieller Höhenzuwachs [m]

ihpot = ((-31.67480+11.64500*ln(t.age+5)-1.04989*((ln(t.age+5))^2.0)- 0.43221*t.hg100 +0.31253*t.hg100*(ln(t.age+5)))-(-31.67480+11.64500* ln(t.age)-1.04989*((ln(t.age))^2.0)-0.43221*t.hg100+0.31253*t.hg100* (ln(t.age))))
(NAGEL 1999)

Höhenzuwachs [m]

t.hinc = t.h*((((-31.67480+11.64500*ln(t.age+5)-1.04989*((ln(t.age+5))^2.0)- 0.43221*t.hg100+0.31253*t.hg100*(ln(t.age+5)))-(-31.67480+11.64500* ln(t.age)-1.04989*((ln(t.age))^2.0)-0.43221*t.hg100+0.31253*t.hg100* (ln(t.age))))/sp.h100)+(0.0000*(t.hinc^1.0)))
(DÖBBELER ET AL 2002)

Höhenzuwachsstreuung [m]

herror = 0.082

Grundflächenzuwachs [cm²]

t.dinc = exp(-5.0479+0.9508*ln(3.14159265359*(t.cw/2.0)/(6.0*(t.h-t.cb)^2)* ((4.0*(t.h-t.cb)^2+(t.cw/2)^2)^1.5-(t.cw/2)^3))-0.7835*ln(t.age) -0.7639*t.c66xy+0.7113*t.c66cxy-0.1891*1.6094)
(DÖBBELER ET AL 2002)

Durchmesserzuwachsstreuung [cm]

derror = 0.649

Maximale Dichte der Grundfläche [m²/ha]

MaxDichte = 0.0001*3.141592/(16*0.0000025729*0.2838*(t.h^(0.6277-0.7621)))
(DÖBBELER 2004)

Maximales Alter [Jahre]

MaxAlter = 280

Zielstärkendurchmesser [cm]

45.0

Höhe der 1. Durchforstung [m]

12.0

Mäßige Durchforstung [m;NB°]

12.0;0.7;18.0;18.0;0.75;24.0;24.0;0.8;100.0

Literatur

Albert, M., (2000): Ein funktionalisierter Höhenergänzungsalgorithmus für Einzelbaumwachstumsmodelle. Jahrestagung des DVFFA -Sektion Ertragskunde- in Kaiserslautern.

Bergel, D. (1973): Formzahluntersuchungen an Buche, Fichte, europäischer Lärche und japanischer Lärche zur Aufstellung neuer Massentafeln. Allg. Forst- u. J. Ztg. 144 (5/6): 117-124

Bergel, D. (1974): Massentafeln II Eiche Roteiche Kiefer. Niedersächsische Forstliche Versuchsanstalt, Göttingen

Bergel, D. (1985): Douglasienertragstafel für Nordwestdeutschland. Niedersächsische Forstliche Versuchsanstalt.

Bergel, D. (1987): Derbholz-Massentafeln III (Nordwestdeutschland). Niedersächsische Forstliche Versuchsanstalt.

Döbbeler, H. (2004): Simulation und Bewertung von Nutzungsstrategien unter heutigen und veränderten Klimabedingungen mit dem Wuchsmodell SILVA 2.2. Dissertation Univ. Göttingen. S. 231

Döbbeler, H.; Albert, M.; Schmidt, M.; Nagel, J. (2002): BWINPro – Programm zur Bestandesanalyse und Prognose. Handbuch zur Version 6.1. Niedersächsische Forstliche Versuchsanstalt Abteilung Waldwachstum, Göttingen, S. 121

Guericke, M. (2001): Versuchsflächenanalyse, Modellparametrisierung und waldbauliche Konsequenzen für die Behandlung von Buchen-Lärchen-Mischbeständen im Südniedersächsischen Bergland. Dissertation Universität Göttingen (http://webdoc.sub.gwdg.de/diss/2001/guericke/index.html)

Hausenauer (2006): Sustainable Forest Management, Growth Models for Europe. Springer, Berlin, Heidelberg

Kramer, H.; Akça; A. (1982): Leitfaden für Dendrometrie und Bestandesinventur, J.D. Sauerländer's Verlag, Frankfurt a.M.

Nagel, J.; Biging, G.S. (1995): Schätzung der Parameter der Weibullfunktion zur Generierung von Durchmesserverteilungen. Allg. Forst- u. J. Ztg. 166 (9/10): 185-189

Nagel, J. (1999): Konzeptionelle Überlegungen zum schrittweisen Aufbau eines waldwachstumskundlichen Simulationssystems für Nordwestdeutschland. Schriften aus der Forstlichen Fakultät der Universität Göttingen und der Nieders. Forstl. Versuchsanstalt, Band 128, J.D. Sauerländer's Verlag, Frankfurt a.M., S.122

Nagel, J.; Albert, M.; Schmidt, M. (2002): Das waldbauliche Prognose- und Entscheidungsmodell BWINPro 6.1. Forst u. Holz 57, (15/16) 486-493

Nagel, J. (2005): TreeGrOSS eine Java basierte Softwarekomponente zur Waldwachstumsmodellierung für Forschung. Lehre und Praxis. Deutscher Verband Forstlicher Forschungsanstalten - Sektion Forstliche Biometrie und Informatik, 15. Tagung Freiburg 9.-10. Oktober 2003, 33-37.

Nagel, J.; Duda, H.; Hansen, J. (2006): Forest Simulator BWINPro7. Forst und Holz 61, Heft 10, S.427-429

Pretzsch, H. 2001: Modellierung des Waldwachstums. Parey Buchverlag Berlin. 341 S.