Die nachhaltige Bewirtschaftung des Waldes stellt eine der wichtigsten Anforderungen an das forstwirtschaftliche Handeln dar. In diesem Kontext befasst sich die vorliegende Arbeit mit der modellhaften Abbildung des Einzelbaumwachstums der Baumart Pinus caribaea var. hondurensis und bildet damit die Grundlage zur Entwicklung nachhaltiger Bewirtschaftungskonzepte in der Dominikanischen Republik.
Die Datenbasis für die Berechnung von Ertragskennwerten und für die sich anschließende Modellierung beruht auf der Anlage von 30 Probeflächen, die in den Jahren 2008, 2009 und 2010 aufgenommen wurden. 16 Probeflächen entfallen dabei auf die Region Jarabacoa und 14 Probeflächen auf die Region Restauración. Nach der Markierung der Probeflächen und der Einzelbäume wurden u.a. der Durchmesser d1,3, die Baumhöhen h, die Kronenbreite kb, der Kronenansatz ka sowie die Baumpositionen erfasst. Aus diesen Messwerten konnten zunächst Ertragskenngrößen berechnet werden. Das Bestandesalter t deckt hierbei einen Bereich von 3 bis 38 Jahren ab. Der Mitteldurchmesser dg variiert zwischen einem Minimalwert von 8,1 cm und einem Maximum von 38,4 cm. Des Weiteren konnten für alle Aufnahmen Mittelhöhen hg mit Werten von 6,7 m bis 33,7 m ermittelt werden. Die Grundfläche G nimmt Werte zwischen 11,3 m² und 63,0 m² an.
Aufbauend auf den Ergebnissen der Basisauswertung und dem Einzelbaumwachstumssimulator BWINPro (NAGEL 1999) wurde anschließend ein regionales Wachstumsmodell angepasst. Zu diesem Zweck waren insbesondere Änderungen in den Bereichen der Datenergänzung, der Zuwachsprognose und zur Nachbildung der Mortalität notwendig, wobei sich an der vorhandenen Modellstruktur orientiert wurde. Im Bereich des Datenergänzungsmoduls war die Anpassung von Funktionen zur Erzeugung einer Durchmesserliste und zur Generierung von Einzelbaumhöhen notwendig. Die Erzeugung einzelner Durchmesser erfolgte auf der Basis der Weibull-Verteilung und linearer Gleichungen, die den Maßstabsparameter b und den Formparameter c dieser Verteilung beschreiben.
In diesem Zusammenhang kamen die Prädiktoren Mitteldurchmesser dg und Maximaldurchmesser dmax zum Einsatz. Zur Generierung von Einzelbaumhöhen wurden zwei Einheitshöhenkurvenmodelle angepasst und auf ihre Eignung hin verglichen. Das Modell von SLOBODA et al. (1993) erwies sich gegenüber dem Modell, das von ARCANGELI et al. (2013) verwendet wurde, als besser geeignet. Zur Erzeugung einer gewissen Höhenvariation bei den generierten Einzelbaumhöhen wurde zudem ein Modell in Anlehnung an ALBERT (2000) angepasst. Zur Berechnung der Standardabweichung der Höhenresiduen um die Bestandeshöhenkurve nach MICHAILOFF (1943) erscheint in diesem Zusammenhang das Modell mit der unabhängigen Variable dg als geeignet.
Kronenbreiten- und Kronenansatzfunktionen werden sowohl im Bereich der Datenergänzung als auch in der Zuwachsprognose benötigt. Bei der Ermittlung eines geeigneten Kronenansatzmodells wurde auf die Funktionen von SCHMIDT (2001), VAN DEUSEN und BIGING (1985) sowie auf die Funktion von PRETZSCH (2001) zurückgegriffen. Die Funktion von SCHMIDT (2001), die auch in der niedersächsischen Version von BWINPro verwendet wird, weist die beste Anpassung an das Datenmaterial auf. In die Findung einer passenden Kronenbreitenfunktion wurden die Modelle aus den Wachstumssimulatoren Silva (PRETZSCH 2001) und BWINPro (DÖBBELER et al. 2011) einbezogen. Hierbei erwies sich die Funktion aus Silva als besser geeignet.
Für die Anpassung der Zuwachsprognose der Baumart P. caribaea var. hondurensis unter den Wuchsbedingungen in der Dominikanischen Republik mussten Funktionen zur Berechnung des Durchmesser- und Höhenzuwachses aktualisiert werden. Das Durchmesserzuwachsmodell beruht hierbei weitestgehend auf dem Ansatz von NAGEL (1999). Als unabhängige Variablen wurden die Kronenmantelfläche km, das Bestandesalter t und der positionsunabhängige Konkurrenzindex C66 nach WENSEL et al. (1987) verwendet. Die Berechnung erfolgte ausschließlich unter Einbeziehung von Bäumen, von denen der Durchmesser d1,3, die Baumhöhe h sowie die Kronenbreite kb und der Kronenansatz ka als Messwerte zur Verfügung standen. Aufgrund dieser Einschränkung umfasste der Stichprobenumfang lediglich 475 Bäume. Mit einem adjustierten Bestimmtheitsmaß von 0,6409 wurde im Vergleich zu anderen Untersuchungen dennoch ein gutes Anpassungsergebnis erzielt.
Eine direkte Schätzung von Einzelbaumhöhenzuwächsen analog zur Durchmesserzuwachsschätzung konnte aufgrund fehlender Höhenzuwächse nicht realisiert werden. Aus diesem Grund wird hierfür ein Verfahren zur Berechnung von Höhenzuwächsen auf der Basis von Einheitshöhenkurven empfohlen. Die Umsetzung des Ansatzes kann in Anlehnung an den Ansatz von KAHN und ĎURSKÝ (1999) erfolgen. Die Voraussetzung zur Verwendung dieses Ansatzes stellt ein Einheitshöhenkurvenmodell und ein Modell zur Beschreibung der Bestandeshöhenentwicklung dar. Um Konsistenz zwischen beiden notwendigen Modellen herzustellen, erfolgte die Parametrisierung auf der Basis der Oberhöhe h100 und des Oberdurchmessers d100. Durch die Verwendung dieser Werte sind die Funktionen zudem stabiler gegenüber Durchforstungen.
Für die Anpassung einer Alters-Höhen-Kurve zur Beschreibung der Oberhöhenentwicklung wurden verschiedene Wachstumsfunktionen getestet. Hierzu zählen die CHAPMAN-RICHARDS-Funktion (RICHARDS 1959), die Funktionen von GOMPERTZ (1825), KORF (1939), SCHUMACHER (1939), SLOBODA (1971) sowie eine logistische Funktion. Nach abschließender Prüfung ergab die SLOBODA-Funktion (1971) den besten Kompromiss zwischen statistischer Anpassung und biologischer Plausibilität. Für das notwendige Einheitshöhenkurvenmodell wurde wie im Bereich der Datenergänzung der Ansatz von SLOBODA et al. (1993) verwendet, wobei als Prädiktoren die Oberhöhe h100 und der Oberdurchmesser d100 in die Berechnung eingingen.
Die Berechnung des Höhenzuwachses eines Einzelbaumes erfolgte durch Differenzbildung zwischen den Einheitshöhen am Anfang und Ende eines Simulationsschrittes.
Zur Gewährleistung der Modellplausibilität und um sicherzustellen, dass die Bestandesdichte nicht ungehindert ansteigt, war die Parametrisierung eines Mortalitätsmodells notwendig. Da innerhalb des Untersuchungszeitraums keine Einzelbaummortalität beobachtet werden konnte und somit die Umsetzung komplexer Modelle nicht möglich war, wurde auf die Bestandesdichte-Regel nach REINEKE (1933) zurückgegriffen. Mithilfe des Mitteldurchmessers dg und der ermittelten Grenzbeziehung können die maximale Stammzahl Nmax und somit auch die maximale Grundfläche Gmax für einen Bestand berechnet werden. Erreicht die Grundfläche eines Bestandes den berechneten Maximalwert, so sterben Bäume mit dem höchsten Konkurrenzindex ab. / The sustainable management of forest ecosystems is one of the most important requirements in forestry decisions. In this context, this thesis deals with the model-based description of individual tree growth of the tree species Pinus caribaea var. hondurensis and provides the basis for the development of sustainable management concepts in the Dominican Republic.
The data for the calculation of stand yield values and for the following modeling process was based on 30 sample plots. The measurements in the pine stands were made in the years 2008, 2009 and 2010. Sixteen sample plots were placed in the region of Jarabacoa and 14 plots were located close to the city of Restauración. After the selection of specific plot size in the stand, all trees were marked and the diameter d1,3, tree height h, crown width kb, crown base ka and tree positions were measured. Using these measurements, yield parameters could be calculated. The stand age t covered a range from 3 to 38 years. The mean quadratic diameter dg varied between a minimum value of 8,1 cm and a maximum of 38,4 cm. Furthermore, mean heights hg with values of 6,7 m to 33,7 m could be determined. The stand basal area G were measured with values between 11,3 m² and 63,0 m².
Based on the results of the calculation of yield parameters and the growth simulator BWINPro (NAGEL 1999), a regional growth model was adjusted. For this objective, adaptations in model routines for data generation, growth prediction and mortality were necessary. The Adaption of routines for the data generation required changes in the algorithms of tree diameter estimation and individual tree height prediction. For the generation of diameter values, the Weibull distribution was adjusted to each diameter distribution of the different stands. In a second step, linear equations which describe the scale parameter b and the shape parameter c of this distribution were fitted by using the independent variables of mean quadratic diameter dg and the maximum diameter dmax. In case where no height measurements taken, a generalized diameter-height equation is parameterized. In this context the two models from SLOBODA et al. (1993) and ARCANGELI et al. (2013) were used and compared. The equation of SLOBODA et al. (1993) showed a better suitability. In order to produce heights with a certain variation, the model after ALBERT (2000) was adjusted. With this model it is now possible, to calculate the standard deviation of the height residuals around the height curve of MICHAILOFF (1943) with the predictor dg.
Crown width and crown base functions are required in the model routines for data generation and growth prediction. For the determination of a suitable crown base model, the functions of SCHMIDT (2001), of VAN DEUSEN und BIGING (1985) and from PRETZSCH (2001) were included in the study. The function of SCHMIDT (2001), which is also used in the version of BWINPro in Lower Saxony, showed the best adaption to the empirical data. For the estimation of crown width values the functions of the growth simulator Silva (PRETZSCH 2001) and BWINPro (DÖBBELER et al. 2011) were tested. Here, the function from Silva proved to be better suited.
To adapt the growth prediction routines for the tree species P. caribaea var. hondurensis under the conditions in the Dominican Republic, the actualization of diameter and height growth models were necessary. The diameter growth model is based on the approach of NAGEL (1999). As independent variables the crown surface area km, stand age t and the position independent competition index C66 by WENSEL et al. (1987) were used. The calculation was performed only with trees, of which the diameter d1,3, tree height h, crown width kb and crown base ka existed as empirical measurements. Because of this limitation, the sample size consists of only 475 trees. A good result of the statistical fit was achieved with an adjusted coefficient of determination of 0,6409. The comparison to other studies supports this assertion.
A direct estimation of individual tree height growth similar to the diameter increment model could not be implemented due to lack of height measurements. For this reason a method for calculating the height growth of single trees on the basis of a generalized diameter-height equation is recommended. The implementation of the described approach can be carried out in accordance to KAHN and ĎURSKÝ (1999). The requirements for using this approach are a generalized diameter-height equation and a growth function for the stand height. To ensure consistency between the two necessary models, the parameterization is based on the dominant height h100 and the dominant diameter d100. By using these values, the functions are also more stable with regard to thinning. For the adaptation of an age-height curve to describe the top height development, different growth functions were tested. These include the CHAPMAN-RICHARDS function (RICHARDS 1959), the functions of GOMPERTZ (1825), KORF (1939), SCHUMACHER (1939), SLOBODA (1971) and a logistic function. After a comprehensive evaluation, the SLOBODA function (1971) represents the best compromise between statistical adaptation and biological plausibility. For the necessary generalized diameter-height equation, again the approach of SLOBODA et al. (1993) was used, with the predictors dominant height h100 and the dominant diameter d100. The calculation of the height increment of a single tree is set by the difference between the heights from the generalized diameter-height equation at the beginning and the end of one simulation step. The age-height curve is used to extract the required information for the generalized diameter-height equation.
To ensure the plausibility of the whole growth model and to ensure that the stand density does not rise unlimited, the adjustment of a mortality model was necessary. Within the research period no individual tree mortality was observed and the implementation of complex models was not possible. For this reason the stand density rule by REINEKE (1933) was implemented. By using the mean quadratic diameter dg and the estimated equation, it is possible to calculate the maximum tree number Nmax. Furthermore the maximum basal area Gmax can be calculated by multiplying Nmax with the basal area of the dg. If the basal area G of the simulated forest stand exceeds the maximum basal area Gmax, the trees with the highest competition index die. The elimination of trees stops if G is lower than Gmax. / El uso sostenible del bosque representa uno de los requisitos más importantes de la gestión forestal. En este contexto, el presente trabajo expone un modelo de crecimiento individual de la especie de árbol Pinus caribaea var. hondurensis, el mismo que pretende contribuir a sentar las bases para el desarrollo de conceptos de aprovechamiento sostenible en la República Dominicana.
La base para el cálculo de valores de rendimiento y para el modelaje la constituyeron los datos tomados en 30 parcelas de muestreo en los años 2008, 2009 y 2010. La ubicación de las parcelas fue georeferenciada; 16 de ellas se encuentran en la región de Jarabacoa y 14 parcelas están en la región Restauración. Cada árbol fue marcado, se registró su ubicación dentro de la parcela y se tomaron las respectivas medidas de diámetro (d1,3), altura (h), diámetro de copa (kb) y altura de fuste limpio (ka). La edad de los árboles en las parcelas de muestreo cubrió un rango desde los tres hasta los 38 años. El diámetro a 1,3 m correspondiente al árbol de área basal promedio (dg) varió entre un mínimo de 8,1 cm hasta un máximo de 38,4cm. También se determinó que la altura de los árboles con área basal promedio (hg) fue de entre 6,7 m hasta 33,7 m. Mientras que el área basal (G) presentó valores entre 11,3 m² y 63,0 m².
A partir de estos resultados de análisis básico y a través de la utilización del programa de simulación del crecimiento individual de árboles BWINPro (NAGEL 1999), se elaboró un modelo de crecimiento específico para la región. Con este propósito, fue necesario realizar modificaciones para calcular valores faltantes, los pronósticos de incremento y para replicar los valores de mortalidad. En cuanto al cálculo de los valores faltantes, se adaptaron las funciones para generar una lista de diámetros y de alturas de cada uno de los árboles. El cálculo del diámetro de cada árbol se basó en la distribución de Weibull y en ecuaciones lineales que describen el parámetro de escala b y el parámetro de forma c de dicha distribución. En este sentido, se utilizó el diámetro a 1,3 m correspondiente al árbol de área basal promedio (dg) y el diámetro máximo dmax como predictores. Con el fin de generar las alturas de cada árbol se adaptaron dos modelos generalizados de diámetro-altura y se avaluó su aptitud. El modelo de SLOBODA et al. (1993) demostró ser más apto para simular la altura de los árboles comparado con el de ARCANGELI et al. (2013). Para producir una cierta variación en la distribución de alturas generadas se ajustó adicionalmente un modelo según ALBERT (2000). El modelo con la variable independiente dg es el que resultó ser más adecuado para calcular las desviaciones estándar de los residuos de las alturas de los árboles con respecto a la ecuacion altura-diámetro de MICHAILOFF (1943).
Tanto para el cálculo de los valores faltantes, como para realizar los pronósticos de incremento, se requieren funciones para determinar la altura de fuste limpio y el diámetro de copa. Con el fin de identificar un modelo adecuado para calcular la altura de fuste limpio se utilizaron las funciones de SCHMIDT (2001), VAN DEUSEN y BIGING (1985), así como también la de PRETZSCH (2001). La función de SCHMIDT (2001), que también es utilizada en la versión de BWINPro de Baja Sajonia, es la que presentó el mejor ajuste a los datos disponibles. En la búsqueda de una función adecuada para calcular el diámetro de copa, se consideraron los modelos de simulación de crecimiento de Silva (PRETZSCH 2001) y BWINPro (DÖBBELER et al. 2011). Entre ellos, el de Silva resultó ser el más adecuado.
Las funciones para calcular el incremento del diámetro y la altura fueron actualizadas para ajustar el pronóstico de incremento de la especie P. caribaea var. hondurensis bajo las condiciones ambientales de la República Dominicana. El modelo para calcular el incremento del diámetro se basa principalmente en el enfoque de NAGEL (1999). Como variables independientes se consideraron el área superficial de copa (km), la edad de la plantación (t) y el índice de competencia (C66) según WENSEL et al. (1987). El cálculo consideró exclusivamente aquellos árboles cuyos valores de diámetro (d1,3), altura (h), diámetro de copa (kb) y altura de fuste limpio (ka) fueron medidos. A causa de esta limitación el tamaño de la muestra fue de únicamente 475 árboles. A pesar de esto, y en comparación con otras investigaciones se logró un buen resultado con un coeficiente de determinación ajustado de 0,6409.
Una estimación directa del incremento de la altura individual análoga a la estimación del incremento del diámetro, no fue posible debido a que faltaban datos del incremento de altura. Por esta razón, el presente trabajo recomienda utilizar un procedimiento para calcular el incremento de la altura basado en una función diámetro-altura. La aplicación de tal procedimiento puede darse de conformidad al enfoque de KAHN und ĎURSKÝ (1999), que utilizaron un modelo generalizado de diámetro-altura y un modelo para describir el desarrollo de la altura del rodal. La parameterización de se basó en la altura dominante (h100) y el diámetro dominante (d100), logrando de tal manera consistencia entre ambos modelos. Además, el uso de estos valores permite que las funciones sean más estables en caso de que se realice un aclareo del bosque. Con el propósito de ajustar de una curva de altura-edad para describir el desarrollo de la altura dominante, se evaluaron varias funciones de crecimiento. Entre ellas se consideraron las de CHAPMAN-RICHARDS (RICHARDS 1959), GOMPERTZ (1825), KORF (1939), SCHUMACHER (1939), SLOBODA (1971), así como también una función logística. Tras evaluar las mencionadas funciones, la de SLOBODA (1971) fue la que presentó el mejor equilibrio entre el ajuste estadístico y la plausibilidad biológica. Para el modelo generalizado diámetro-altura se utilizó, al igual que para el cálculo de valores faltantes, el enfoque de SLOBODA et al. (1993), en el cuyos cálculos se incluyen como predictores la altura dominante (h100) y el diámetro dominante (d100). El incremento de altura de un individuo resultó de la diferencia entre la altura de cada árbol al inicio y al final de cada iteración de la simulación.
Con el objeto de garantizar la plausibilidad del modelo y, para evitar que la densidad del rodal crezca libremente fue necesario parameterizar un modelo de mortalidad. Debido a que durante el tiempo de estudio no se observó la muerte de ningún árbol, y por lo cual no fue posible utilizar un modelo más complejo, se empleó la ley de densidad del rodal de Reineke, según REINEKE (1933). Con ayuda del diámetro a 1,3 m correspondiente al árbol de área basal promedio (dg) y a una relación lineal, se puede calcular el número máximo de árboles por hectárea y así también el área basal máxima de un rodal. Cuando el área basal del rodal alcanza el valor máximo calculado, mueren los árboles que poseen los índices de competencia más altos.
| Identifer | oai:union.ndltd.org:DRESDEN/oai:qucosa:de:qucosa:29288 |
| Date | 03 August 2015 |
| Creators | Körner, Michael |
| Contributors | Röhle, Heinz, Guericke, Martin, Kapp, Gerald, Technische Universität Dresden |
| Source Sets | Hochschulschriftenserver (HSSS) der SLUB Dresden |
| Language | German |
| Detected Language | German |
| Type | doc-type:doctoralThesis, info:eu-repo/semantics/doctoralThesis, doc-type:Text |
| Rights | info:eu-repo/semantics/openAccess |
Page generated in 0.0045 seconds