Die neuen Ertragstafeln der NW-FVA (Albert et al. 2021) bieten für die fünf Hauptbaumarten Eiche, Buche, Fichte, Douglasie und Kiefer ein Planungsinstrument, welches das aktuelle Zuwachsniveau der Wälder Nordwestdeutschlands widerspiegelt und das empfohlene waldbauliche Behandlungskonzept einer gestaffelten Hochdurchforstung unterstellt.
Neben den Tafelwerten bietet et.nwfva
sowohl die Möglichkeit zu
bonitieren als auch Bestandeshöhen und Tafelwerte für beliebige
Bonitäten und Alter zu bestimmen. Für die Inter- und Extrapolation
kommen standardmäßig funktionale Bonitätsfächer zum Einsatz, aber die
klassische Methode mittels Dreisatz steht ebenfalls zur Verfügung. Die
funktionale Interpolation kann den Bereich von der -3. bis 7.
Ertragsklasse und den gesamten Altersbereich bedienen. Die klassische
Interpolation ist auf den Bonitätsbereich -2. bis 4. Ertragsklasse und
den Altersbereich 5 bis max. zulässiges Alter (Ei 220, Bu 180 und Fi,
Dgl, Ki 160) beschränkt.
Die aktuelle Version des Paketes kann für die gängigen Betriebssysteme (Windows, Linux, macOS) von CRAN
install.packages("et.nwfva")
und aus dem R-Universe der NW-FVA installiert werden.
# NW-FVA Universe in R bekannt machen
options(repos = c(
nwfva = "https://nw-fva.r-universe.dev",
CRAN = "https://cloud.r-project.org"))
# und Paket wie gewohnt installieren
install.packages("et.nwfva")
Die Entwicklerversion kann mit dem Paket remotes
direkt aus
Github kompiliert und installiert
werden.
remotes::install_github("rnuske/et.nwfva")
Einen Überblick über die vorhandenen Tafeln liefert et_liste()
und
Detailinformationen zu einer Tafel et_info()
. Ganze Tafeln einer
Baumart oder Auszüge aus einer Tafel, z.B. für bestimmte Bonitäten oder
Alter, bietet et_tafel()
.
library(et.nwfva)
# Ertragstafelwerte für Buche 1. Ertragsklasse ausgeben (hier: Zeilen 1 bis 6)
head(et_tafel("Buche", bon=1))
#> Ekl Alter N Hg H100 G Dg Dw V N_aus G_aus Dg_aus V_aus iV GWL dGZ
#> 1 1 35 3028 9.7 12.9 13.2 7.5 11.2 54 1518 5.5 6.8 24 9.5 111 3.2
#> 2 1 40 2083 12.3 15.3 15.2 9.7 14.4 82 945 4.9 8.1 28 11.1 167 4.2
#> 3 1 45 1487 14.7 17.5 16.9 12.0 17.7 112 596 4.4 9.7 31 12.2 228 5.1
#> 4 1 50 1104 16.9 19.5 18.3 14.5 21.1 143 382 4.1 11.7 34 13.0 293 5.9
#> 5 1 55 879 18.9 21.4 19.6 16.8 24.5 174 225 3.9 14.9 36 13.5 361 6.6
#> 6 1 60 721 20.7 23.1 20.8 19.1 27.8 206 159 3.7 17.2 38 13.9 430 7.2
# Eine Ertragstafelzeile interpoliert für Fichte 2.4. Ekl. und Alter 82
et_tafel("Fi", bon=2.4, alter=82)
#> Ekl Alter N Hg H100 G Dg Dw V N_aus G_aus Dg_aus V_aus iV GWL dGZ
#> 1 2.4 82 672 25.3 26.3 39.1 27.4 38.9 436.8 42 3 30.1 32.8 12.9 768.6 9.4
Die Bonitierung von Beständen erfolgt mit et_bonitaet()
anhand von
Bestandesalter und -höhe standardmäßig über funktionale Bonitätsfächer.
Auf Wunsch ist auch die klassische Bonitierung mittels Dreisatz
verfügbar. Im Fall der Interpolation zwischen den Tafelwerten
unterscheiden sich die Verfahren kaum, aber im Extrapolationsbereich
können die Abweichungen durchaus beachtlich sein (siehe Abbildung oben).
# Ertragsklasse eines 80 jährigen Eichenbestandes mit der Mittelhöhe 16m
et_bonitaet("Eiche", alter=80, hoehe=16)
#> [1] 3.8
et_bonitaet("Eiche", alter=80, hoehe=16, methode="klassisch")
#> [1] 3.9
Sollte im Hinblick auf das Alter und/oder die Bestandeshöhe eine
Extrapolation der Tafelwerte notwendig sein, gibt die klassische Methode
einen Hinweis aus. Falls dabei der zulässige Extrapolationsbereich
verlassen würde, geben beide Methoden eine Warnung und je nach Wert des
Parameters kapp_na
den jeweiligen Grenzwert oder NA
zurück.
# Ertragsklasse eines 120 jährigen Kiefernbestandes mit Oberhöhe 32m
et_bonitaet("Kiefer", alter=120, hoehe=45, hoehe_typ="ober", methode="funktional")
#> [1] -2.4
et_bonitaet("Kiefer", alter=120, hoehe=45, hoehe_typ="ober", methode="klassisch")
#> Für die Ertragsklasse(n) '-1' ist das höchste Alter in der Tafel '90'. Es wurde auf das Alter '120' extrapoliert.
#> [1] -1.8
# Anfrage außerhalb des zulässigen Bonitätsbereichs
et_bonitaet("Kiefer", alter=70, hoehe=40, hoehe_typ="ober")
#> Warning: Die Bestandeshöhe 40 im Alter 70 ergibt eine Bonität außerhalb des Intervalls
#> [-3,7]. Da kapp_na=TRUE, wurde die Bonität auf NA gesetzt.
#> [1] NA
Die Bestandeshöhenwerte für die Abbildungen oben wurden mit der Funktion
et_hoehe()
berechnet. Umrechnungen von relativen Ertragsklassen in
absolute Oberhöhenbonitäten (Oberhöhe H100 im Alter 100 in Metern) und
umgekehrt ermöglicht die Funktion et_bontrans()
.
# Bestandesmittelhöhe für Buche im Alter 75 und 1.25. Ekl.
et_hoehe("bu", alter=75, bon=1.25, hoehe_typ="mittel")
#> [1] 24.4
et_hoehe("bu", alter=75, bon=1.25, hoehe_typ="mittel", methode="klass")
#> [1] 24.4
# Umrechnung von der Oberhöhenbonität 27 zur relativen Ertragsklasse für Kiefer
et_bontrans("Kiefer", bon=27, richtung="abs_zu_rel")
#> [1] 1.5
Albert M., Nagel J., Schmidt M., Nagel R.-V., Spellmann H. (2021): Eine neue Generation von Ertragstafeln für Eiche, Buche, Fichte, Douglasie und Kiefer [Datensatz]. Version 1.0. Zenodo. https://doi.org/10.5281/zenodo.6343906