iForest - Biogeosciences and Forestry

iForest - Biogeosciences and Forestry

Incorporating management history into forest growth modelling

iForest - Biogeosciences and Forestry, Volume 4, Issue 5, Pages 212-217 (2011)
doi: https://doi.org/10.3832/ifor0597-004
Published: Nov 03, 2011 - Copyright © 2011 SISEF

Research Articles

Collection/Special Issue: COST Action FP0903 (2010) - Rome (Italy)
Research, monitoring and modelling in the study of climate change and air pollution impacts on forest ecosystems
Guest Editors: E Paoletti, J-P Tuovinen, N Clarke, G Matteucci, R Matyssek, G Wieser, R Fischer, P Cudlin, N Potocic

Mechanistic modelling is an important tool for understanding the impacts of climate change and pollutants on forest growth. One of the common practical limitations of these models is a lack of specific information regarding management activities such as thinning or harvesting, which can have a very strong influence on the accuracy of results. The use of inventory data for model parameterization and calibration is also problematic, as inventories are designed to have large volumes of data amalgamated to give accurate mean results across large areas. The precision of single point estimates is often quite low.This study uses BIOME-BGC to model forest growth on 1133 sites of the Austrian National Forest Inventory, and develops a method to estimate timber removal patterns prior to the commencement of record keeping on the sites. Recognizing the poor precision of individual point estimates in the data, we do not seek to precisely calibrate the model to the data on each point. Rather, we assume that the point-wise inventory estimates will be normally distributed around the true values. We then model each site assuming no management interventions, and compare this with inventory results. Plotting the “error” between model results and NFI data shows a strong right-skew, reflecting the modelled lack of timber removals. A Box-Cox transformation of the error plot, centred on zero, would represent an unbiased model estimate of the data, thus we can determine the historic timber removals as the difference between the original error curve and its Box-Cox transformation. Calibrating the model with this information allow us to represent forest volume with greater accuracy than would otherwise be possible.

BIOME-BGC, Inventory, Uncertainty, Thinning, Model initialisation


Mechanistic modelling is an increasingly popular method of tracking carbon fluxes in forest ecosystems, important for carbon balance studies, climate and pollution studies and others. Properly calibrated and validated, such models can also be used as diagnostic tools, to separate the influence of forest growth impactors such as climate change and Nitrogen deposition ([8]). Such models have however been sometimes criticised in the past for a poor representation of forest mass growth ([29], [26]). This paper presents part of a wider study that will compare aboveground biomass estimates obtained with the BIOME-BGC model ([33], [23]) with field data from the Austrian National Forest Inventory ([9]). We detail here the methodology developed to estimate forest management activities that may have occurred prior to the starting date of the available inventory data.

BIOME-BGC simulates fluxes and pools of carbon, water and nitrogen, beginning with a spin-up run of many forest growth cycles, to estimate a natural forest stand for the site without human influence ([24]). We follow this with two cycles of clear cutting to represent early industrial forest use in Central Europe, with the final start date of the present forest established in accordance with the mean forest age for each site from the NFI data. In managed forests it is likely that some biomass removals have occurred in the current forest rotation, but if these were prior to the commencement of the inventory programme then the timing and extent of those thinnings is unknown. As forest management is one of the major influences on forest carbon fluxes and pools ([30], [11], [17]), this represents a potentially serious source of uncertainty or bias in model outputs ([38]). In a recent Europe-wide comparison of forest models, Tupek et al. ([34]) did not include management in their BGC simulation, yet recognise management as a major factor.

Direct site-by-site comparisons of model outputs with inventory measurements is difficult, due to the fundamentally different approaches to scale. Models assume homogeneity at the scale of their input data ([14]), where single sites are assumed to be representative of a broader area. Inventories on the other hand are a statistical measure, designed so that a large number of results are aggregated to give a mean figure for a region or strata. Although at a broad scale the results should be comparable, individual site-by-site comparisons are likely to have little meaning. In terms of using inventory data for model initialisation and validation, this means that some liberties must be taken with the strictly mechanistic operation of biogeochemical process models, as a large component of random error in the individual plot observations must be incorporated.

National forest inventories (NFIs) contain a vast body of information, however the use of this data for research purposes has to date been limited ([18]), and few biogeochemical analyses have been applied to inventory data ([27]). This is likely due to two complications, the first to do with the data needs of complex biogeochemical models and the second with the statistical nature of forest inventories. Comprehensive inventories such as the Austrian NFI ([9]), species-specific model parameterisation ([23]) and modern spatial interpolation methods ([21]) may fulfil most of the data needs, but the uncertainty of prior management activities remains.

Comparison of biogeochemical modelling and inventory results in managed forests requires solving two simultaneous problems: the statistical uncertainty of the inventory data at the plot scale and the lack of knowledge of historical management. Our approach in this paper is to assume that the errors in a sophisticated NFI and a mature, well-tested model will both be normally distributed around zero, and thus, if assumptions of plot management history are correct, the site-by-site disagreements between the two approaches will also be normally distributed around zero.

  Data and methods 

The Austrian National Forest Inventory ([9]) is organized into tracts each of 4 points on a 200m square. 5600 such tracts are arranged in a square grid pattern across the country, including over areas that are not currently forested. An angle-count with a basal area factor of 4 is performed on each point, recording all trees of over 10.5cm dbh. A fixed-area plot of radius 2.60 metres is installed to record all trees of dbh between 5.0 and 10.5cm. Twenty percent of plots were measured each year, and re-measured after 5 years. Timber volumes are calculated here from measured diameters and heights using the allometric equations of Pollanschütz ([25]). Using one point from each forested tract, a total of 2224 points is obtained. We use here a subset of the database for the measurement periods 1981-85 and 1986-91, and limit the analysis to those plots with forest present in both periods, a species homogeneity of over 80 percent by basal area and less than 30 percent timber removal from the plot between the measurement periods, as our current version of BIOME-BGC has not been tested for mixed-species forests and a major removal within the increment period could confuse the results. This limits our analysis to 1133 plots.

Separate species parameterisations are available for Norway spruce (Picea abies) depending on site elevation above or below 1000m asl ([23]), but as no parameters are available for Silver fir (Abies alba) these were modelled as low elevation spruce giving a total of 386 plots (354 spruce dominated, 32 fir) in this group. Some 439 plots were high elevation spruce, 80 Pinus (67 Scots pine Pinus sylvestris and 13 Pinus nigra), 47 Larch (Larix decidua), 28 combining Oak (Quercus robur/petraea) and Ash (Fraxinus excelsior), 146 beech (Fagus sylvatica) and other broadleaves and seven stone pine (Pinus cembra).

BIOME-BGC is a mechanistic model developed to simulate the ecosystem processes of a forest stand on a daily timestep. Storages and fluxes of water, carbon and nitrogen are tracked throughout various pools in the vegetation, litter, and soil. The model has been widely applied across a range of forested ecosystems around the world, has been extensively compared with flux measurements ([32]), and was specifically parameterized and validated for major central European species by Pietsch et al. ([23]). Norway spruce in particular has been widely studied with BIOME-BGC ([4], [15], [22], [12]), giving confidence that the model can accurately represent spruce growth across a range of different circumstances. Forest management operations were included in BIOME-BGC by Cienciala & Tatarinov ([4]) and Petritsch et al. ([22]), based on prescribed management regimes.

Net primary production of the 1133 plots in our dataset is simulated using the Austrian incarnation of the DAYMET climate interpolation ([19], [7]). Species and age information is obtained from the NFI, and other inputs comprise of an interpolation of Austrian soil data ([21]), a 1 km² nationwide raster of nitrogen deposition in 1995 ([28]) and species-specific parameters determined by Pietsch et al. ([23]). BIOME-BGC operates by simulating the development of an unmanaged forest over many centuries, until a point where carbon pools are meta-stable over the course of successive mortality cycles ([24]). This is followed by two cycles of clearcutting and replanting to simulate early-industrial forest exploitation in Central Europe, and the planting date for the current stand is determined from the NFI data.

The primary output of BIOME-BGC is net primary production (NPP). Internal routines allocate this in fixed proportions to the various tree compartments, of which we are here interested in aboveground wood carbon. Without consideration of current rotation management, the model output would apply to a potential, unmanaged forest condition. In most cases this gives a biomass estimate substantially higher than the recorded data from the NFI. At the plot level the NFI data is extremely imprecise, so a four step procedure was developed to determine the appropriate removal assumptions to include in the modelling. These assumptions are then used in a second simulation of forest growth to determine final results (Fig. 1).

Fig. 1 - Schematic diagram of modelling procedure. The first run of the BIOME-BGC model is conducted without any management assumptions. From the comparison of those outputs with the measured standing volumes from the National Forest Inventory, it is possible to estimate the management history on each site. This information is then applied in the second run of the model.

  Enlarge/Shrink   Download   Full Width  Open in Viewer

(1) The true standing biomass in Period 1 was estimated as the mean of the NFI results for Periods 1 and 2, less half of the plot increment calculated from the NFI data using the “Starting Value” method ([10], [31]). This “temporal averaging” reduces some of the extreme variability in single-point angle count measurements ([6]).

(2) The “unmanaged” model outputs may be compared with the estimated Period 1 biomass, and the mismatch (“error”) calculated. A density plot of this mismatch shows a strong right-skew (most BGC estimates were much too high, fewer were too low). Plots were then stratified by dominant species group and data/model mismatches calculated.

(3) We assume that the mismatch between the NFI biomass estimations and the model outputs should be normally distributed. A Box-Cox power transformation ([2]) was applied to the mismatch data for each species group, and the thinning requirement for each plot calculated as the difference between the mismatch curve and the transformed curve, centred on zero (Fig. 2), less a minor adjustment to compensate for the extra soil nitrogen use in a thinned stand.

Fig. 2 - Relationship of model output to NFI data, by plot. The original “error” curve (closed symbols, dark blue) shows a clear positive skew, reflecting a model overestimation of volume with no assumed thinnings. Applying a Box-Cox power transformation to this curve increases the normalcy of the error distribution (open symbols, pink). The required thinning for any plot (in m³/ha) is the lateral distance between that plot’s corresponding points on the two curves. Note that the transformed curve has been scaled vertically down for this figure, to increase the clarity of the concept.

  Enlarge/Shrink   Download   Full Width  Open in Viewer

(4) Removals calculated in step 3 are modelled as thinnings taken place some time prior to Period 1 (although they also may have been partial harvests or significant disturbance events). To avoid these having a biophysical impact on modelling post-Period 1 due to increased nitrogen consumption by soil bacteria, we apply a series of steps to determine reasonable assumptions for the timing of the removals.

  1. If the stand is over 50 years old, first removal is assumed to be at age 40 ([20]).
  2. If insufficient biomass can be removed in a single step (required removals are greater than modelled standing biomass), subsequent removals occur at ages 70, 100 and 130 and 5 years prior to the Period 1 inventory date.
  3. If a subsequent removal is required, the prior removal is limited to 40 percent of biomass.
  4. A maximum of 95% of biomass may be removed in the final step.

To convert biomass figures to more readily visualised timber volume per hectare values, we compiled a regression of all 33 779 tree measurements in the available dataset, calculating conversion factors between biomass estimates made with allometry and biomass expansion factors developed for Central European species (Lexer, pers. comm. - 10 February 2010) and merchantable timber volume estimates consistent with Pollanschütz ([25]) - Tab. 1. All species are included here, as some of our 1133 plots contained a minor component of other species. It should be noted however that for species other than those in our major groups often few data points were available for the biomass/timber volume regression, and thus the multipliers may be unreliable for those species. The Lexer equations were used for consistency with concurrent studies ([35]), as the only available comprehensive estimates for Austrian conditions.

Tab. 1 - Multipliers for conversion of above-ground woody biomass (kg/m2) to volume (m3/ha) consistent with Pollanschütz ([25]).

Species Multiplier Species Multiplier
Spruce (Picea abies), low elevation 18.91 Chestnut (Castanea sativa) 12.31
Fir (Abies alba) 18.28 Black Locust (Robinia pseudoacacia) 11.11
Larch (Larix decidua) 21.46 Sorbus/Prunus spp. 9.94
Scots Pine (Pinus sylvestris) 22.12 Birch (Betula spp.) 11.01
Black Pine (Pinus nigra) 18.29 Black Alder (Alnus glutinosa) 16.93
Stone Pine (Pinus cembra) 7.45 White Alder (Alnus incana) 10.66
Douglas Fir (Pseudotsuga menziesii) 6.25 Linden (Tilia cordata) 15.01
Other conifers 6.81 Aspen, Poplars (Populus spp.) 13.74
Beech (Fraxinus excelsior) 13.17 Black Poplar (Populus nigra) 12.58
Oak (Quercus robur/petraea) 12.53 Hybrid Poplars (Populus x spp.) 11.83
Hornbeam (Carpinus betulus) 11.61 Willow (Salix spp.) 8.95
Ash (Fraxinus excelsior) 13.61 Other broadleafs 9.06
Maple (Acer spp.) 11.41 Spruce, high elevation 15.9
Elm (Ulmus spp.) 11.53 - -

  Enlarge/Reduce  Open in Viewer

Given the very poor precision of angle-count estimates of volume on single points, a regression of observed versus predicted standing volume is inappropriate. We therefore apply the methods of Bland & Altman ([1]), which for each data point plot the difference between values determined through each method against their mean for each point.

  Results and discussion  

The initial simulation of forest growth on the NFI plots was performed without any management interventions. A comparison of model outputs with observed data displayed a remarkably even normal curve for those plots where the model appears to underestimate the NFI (the blue line below zero in Fig. 2). This gives some confidence that if management were appropriately considered, the variation between the simulations and the NFI data would also be normally distributed for apparent overestimations.

Overall, the procedure enabled a close simulation of standing volume (Fig. 3). Some problems exist with the Larch simulation, probably due to poor species data resulting from the low number of samples used in the initial species parameterisation ([23], Pietsch, pers. comm. - 25 June 2010). Results for the Oak/Ash and Stone Pine groups are rather uncertain due to the low numbers of plots available in this study for those species. The Bland-Altman plots (Fig. 4) show that the errors are consistent over the range of standing volumes, with little evidence of heteroscedasticity.

Fig. 3 - Standing volume in Period 1. Error bars show one standard deviation above and below the mean.

  Enlarge/Shrink   Download   Full Width  Open in Viewer

Fig. 4 - Bland-Altman plot of standing volume determined with inventory and with the BIOME-BGC model. In each plot the “x” axis is the mean of the volume assessed with each method, while the “y” axis is the difference between the methods (BGC subtract inventory). Dashed lines indicate the mean difference.

  Enlarge/Shrink   Download   Full Width  Open in Viewer

The Lexer biomass equations are based only on dbh. Several authors have presented allometry for major European species in different regions ([13] and [37] for Fagus sylvatica, [36] for Picea abies and [3], [5] for Pinus sylvestris and Quercus spp. respectively). All of these authors concluded that allometric models can be considerably improved by including height and other variables, rather than basing results solely on dbh. The Austrian NFI however contains many records where height is calculated from dbh rather than directly measured, which would appear to negate any possible advantage of using more complicated allometrics.

Besides natural variability, several other factors may account for the differences between model output and NFI records. The version of BIOME-BGC we use in this study is not designed to simulate mixed-species forest, so the up to 20 percent of non-homogenous species on any plot may have some minor effect. Uncertainty in biomass expansion factors ([16]) is also possibly biasing on some plots, where the by-plot ratio of volume to biomass is not the same as the mean by-tree relationship for the dominant species (Tab. 1).

Although it may seem that the results presented here are an artificial “forcing” of model assumptions to match the NFI data, in most cases the simulated thinnings are timed to occur several decades prior to the inventory measurement. The model is thus required to accurately simulate growth over long periods before the results in Fig. 3 are extracted. Alternative approaches such as thinning to precisely match field data ([28]) would not be possible for apparent model underestimations and thus would give an underestimation in the mean. Thinning all plots equally so that the simulation and the field data have equal means would simply shift the blue curve in Fig. 1 to the left, giving a non-normal error distribution. An “expert opinion” approach ([20]) would ignore the standing volume information available in the NFI dataset. Although single plot data in the NFI is extremely imprecise, in the aggregate it can be assumed to provide accurate values.

By thinning to target a normal error curve, apparent model underestimations are little affected. In cases where a plot simulated without thinning greatly overestimates standing volume a heavy or a series of thinnings is indicated. There are any number of possible normal error curves that could have been selected as the “target” error curve, but the procedure using the Box-Cox transformation was chosen to maximise the use of the available data. The timing of the thinning events at 40, 70, 100 and 130 years of age appears to be reasonably consistent with records from the NFI (Fig. 5). A breakdown of thinning volume removals (Tab. 2) suggests that the assumption procedure tends to underestimate the number of thinnings, but overestimate the volume removed in each thinning event. The results shown in Fig. 3 however suggest that this has a minimal effect on the final simulated standing volumes.

Fig. 5 - Frequency of removal events in each age class between the first two periods of the modern Austrian National Forest Inventory. Tick marks on the “x” axis correspond to the assumed timing of thinnings for BIOME-BGC model initialization.

  Enlarge/Shrink   Download   Full Width  Open in Viewer

Tab. 2 - Comparison of determined thinning assumptions with NFI records between the first and second modern Austrian National Forest Inventory periods. Inventory values for the number of thinning events and the sum of volume removed are the annual mean of the 5-year inter-period, while the assumptions values are the annual mean of the 20 years from 1961-1980.

Group Inventory Thinning assumptions
n yr-1 mean
(m3 ha-1)
(m3 ha-1yr-1)
n yr-1 mean
(m3 ha-1)
(m3 ha-1yr-1)
All 46.00 55.27 2542.33 22.95 69.24 1589.10
Spruce < 1000m and fir 21.20 55.22 1170.69 12.60 84.01 1058.55
Spruce > 1000m 15.00 64.25 963.75 6.55 53.10 347.80
Larch 0.20 15.17 3.03 0.00 0.00 0.00
Scots pine 4.00 42.01 168.02 0.15 23.67 3.55
Beech 5.00 42.04 210.20 3.65 49.10 179.20
Oak/ash 0.40 51.47 20.59 0.00 0.00 0.00
Stone pine 0.20 30.24 6.05 0.00 0.00 0.00

  Enlarge/Reduce  Open in Viewer


The incorporation of unknown forest management histories into ecosystem modelling is a challenging issue, yet the ability of models to accurately mimic forest volume growth is dependant on management being appropriately considered. As we will rarely know the history of a large number of sites in sufficient detail to use real figures for historic thinnings, supportable procedures must be developed for appropriate assumptions. While National forest inventories provide large datasets of potential calibration data, the statistical nature of inventories means that we cannot assume a single inventory point to be representative of a wider area and directly comparable with model outputs for that point. The statistical estimation of management history in this work does not attempt to describe the real, unknown thinning history of each stand. Rather, it establishes a set of assumptions that can mimic the true management history in such a way that we can produce unbiased estimates of standing volume which, at a sufficiently large scale, are compatible with National Forest Inventories.

The procedure developed here has provided us with model-driven timber volume estimates for 1133 sites across the whole of Austria that are overall a close match to NFI records at the national scale, are of comparable precision, and maximise the use of available data. Future work will examine how well model simulations compare with NFI records for subsequent inventory periods, where statistical indications are available regarding management interventions on each plot.


This work was supported by the two projects:Auswirkungen des Klimawandels auf Österreichs Wälder - Entwicklung und vergleichende Evaluierung unterschiedlicher Prognosemodelle (WAMOD) and Comparing satellite versus inventory driven carbon estimates for Austrian Forests (MOTI), both funded by the Energy Fund of the Federal State of Austria and managed by Kommunalkredit Public Consulting GmbH. We thank the two anonymous reviewers and the editor of the journal for helpful suggestions to improve this manuscript.


Bland JM, Altman DG (1986). Statistical methods for assessing agreement between two methods of clinical measurement. Lancet 327 (8476): 307-310.
CrossRef | Gscholar
Box GEP, Cox DR (1964). An analysis of transformations. Journal of the Royal Statistical Society Series B 26 (2): 211-252.
Cienciala E, Cerný M, Tatarinov F, Apltauer J, Exnerová Z (2006). Biomass functions applicable to Scots pine. Trees 20: 483-495.
CrossRef | Gscholar
Cienciala E, Tatarinov F (2006). Application of BIOME-BGC model to managed forests 2. Comparison with long-term observations of stand production for major tree species. Forest Ecology and Management 237: 252-256.
CrossRef | Gscholar
Cienciala E, Apltauer J, Exnerová Z, Tatarinov F (2008). Biomass functions applicable to oak trees grown in Central-European forestry. Journal of Forest Science 54: 109-120.
Eastaugh CS, Hasenauer H (2010). The usefulness of time series angle-count forest inventory data in assessing forest growth model accuracy. Forestry Ideas 16 (2) 171-180.
Eastaugh CS, Petritsch R, Hasenauer H (2010). Climate characteristics across the Austrian forest estate from 1960 to 2008. Austrian Journal of Forest Science 127 (3): 133-146.
Eastaugh CS, Pötzelsberger E, Hasenauer H (2011). Assessing the impacts of climate change and nitrogen deposition on Norway spruce (Picea abies L. Karst) growth in Austria with BIOME-BGC. Tree Physiology 31 (3): 262-274.
CrossRef | Gscholar
Gabler K, Schadauer K (2006). Methoden der Österreichischen Waldinventur 2000/02. BFW-Berichte 135, Federal Research and Training Centre for Forests, Natural Hazards and Landscape, Vienna, Austria, pp. 135.
Grosenbaugh LR (1958). Point-sampling and line-sampling: probability theory, geometric implications, synthesis. USDA Forest Service, South Forest Experimental Station, Occasional Paper 160, pp. 34.
Jandl R, Lindner M, Vesterdal L, Bauwens B, Baritz R, Hagedorn F, Johnson DW, Minnkkinen K, Byrne K (2007). How strongly can forest management influence soil carbon sequestration? Geoderma 17: 253-268.
CrossRef | Gscholar
Jochheim H, Puhlmann M, Beese F, Berthold D, Einert P, Kallweit R, Konopatzky A, Meesenburg H, Meiwes KJ, Raspe S, Schulte-Bisping H, Schulz C (2009). Modelling the carbon budget of intensive forest monitoring sites in Germany using the simulation model BIOME-BGC. iForest 2: 7-10.
CrossRef | Gscholar
Joosten R, Schumacher J, Wirth C, Schulte A (2004). Evaluating tree carbon predictions for beech (Fagus sylvatica L.) in western Germany. Forest Ecology and Management 189: 87-96.
CrossRef | Gscholar
Landsberg JJ, Kaufmann MR, Binkley D, Isebrands J, Jarvis PG (1991). Evaluating progress towards closed forest models based on fluxes of carbon, water and nutrients. Tree Physiology 9: 1-15.
Lange H (2007). Modelling carbon dynamics in forest ecosystems using Biome-BGC. In: “Greenhouse gas budget of soils under changing climate and land use” (Jandl R, Olsson M eds). BFV, Wien, Austria, pp. 63-70.
Lehtonen A, Cienciala E, Tatarinov F, Mäkäpää R (2007). Uncertainty estimation of biomass expansion factors for Norway spruce in the Czech Republic. Annals of Forest Science 64: 133-140.
CrossRef | Gscholar
Lindner M, Green T, Woodall CW, Perry CH, Nabuurs GJ, Sanz MJ (2008). Impacts of forest ecosystem management on greenhouse gas budgets. Forest Ecology and Management 256: 191-193.
CrossRef | Gscholar
Nabuurs GJ, Hengeveld GM, van der Werf DC, Heidema AH (2010). European forest carbon balance assessed with inventory based methods - an introduction to a special section. Forest Ecology and Management 260: 239-240.
CrossRef | Gscholar
Petritsch R (2002). Anvendung und Validierung des Klimainterpolationsmodells DAYMET in Österreich. MSc thesis, University of Natural Resources and Applied Life Sciences, Institute of Silviculture, Vienna, Austria, pp. 95.
Petritsch R (2008). Large scale mechanistic ecosystem modeling in Austria. Ph.D. thesis, University of Natural Resources and Applied Life Sciences, Institute of Silviculture, Vienna, Austria, pp. 135.
Petritsch R, Hasenauer H (2007). Interpolating input parameters for large scale ecosystem models. Austrian Journal of Forest Science 124 (2): 135-151.
Petritsch R, Hasenauer H, Pietsch SA (2007). Incorporating forest growth response to thinning within biome-BGC. Forest Ecology and Management 242 (1-2): 324-336.
CrossRef | Gscholar
Pietsch SA, Hasenauer H, Thornton PE (2005). BGC-model parameters for tree species growing in central European forests. Forest Ecology and Management 211: 264-295.
CrossRef | Gscholar
Pietsch SA, Hasenauer H (2006). Evaluating the self-initialization procedure for large-scale ecosystem models. Global Change Biology 12: 1658-1669.
CrossRef | Gscholar
Pollanschütz J (1974). Formzahlfunktionen der Hauptbaumarten Österreichs. Allgemeine Forstzeitung 85: 341-343.
Pretzsch H, Grote R, Reinking B, Rötzer TH, Seifert ST (2008). Models for forest ecosystem management: A European perspective. Annals of Botany 101: 1065-1087.
CrossRef | Gscholar
Purves D, Pacala S (2008). Predictive models of forest dynamics. Science 320: 1452-1453.
CrossRef | Gscholar
Pötzelsberger E (2008). Assessing the productivity and water regime in the Schmittenbach catchment area. Diploma thesis, University of Natural Resources and Applied Life Sciences, Institute of Silviculture,Vienna, Austria, pp. 75.
Rennolls K, Tomé M, McRoberts RE, Vanclay J, LeMay V, Guan BT, Gertner G (2007). Potential contributions of statistics and modelling to sustainable forest management: review and synthesis. In: “Sustainable forestry: from monitoring and modelling to knowledge management and policy science” (Reynolds KM, Thomson AJ, Kohl M, Shannon MA, Ray D, Rennolls K eds). CABI, pp. 525.
Rypdal K, Baritz R (2002). Estimating and managing uncertainties in order to detect terrestrial greenhouse gas removals. CICERO Working Paper 2002:07, pp. 10.
Online | Gscholar
Schieler K (1997). Methode der Zuwachsberechnung der Österreichischen Waldinventur. Ph.D. Thesis, University of Natural Resources and Applied Life Sciences Vienna, Institute for Forest Growth, Vienna, Austria, pp. 99.
Schimel D, Melillo J, Tian H, McGuire AD, Kicklighter D, Kittel T, Rosenbloom N, Running S, Thornton P, Ojima D, Parton W, Kelly R, Sykes M, Neilson R, Rizzo B (2000). Contribution of increasing CO2 and climate to carbon storage by ecosystems in the United States. Science 287 (5460): 2004-2006.
CrossRef | Gscholar
Thornton PE (1998). Description of a numerical simulation model for predicting the dynamics of energy, water carbon and nitrogen in a terrestrial ecosystem. PhD thesis, University of Montana, Missoula, USA, pp. 280.
Tupek B, Zanchi G, Verkerk PJ, Churkina G, Viovy N, Hughes JK, Lindner M (2010). A comparison of alternative modelling approaches to evaluate the European forest carbon fluxes. Forest Ecology and Management 260: 241-251.
CrossRef | Gscholar
WAMOD (2010). Auswirkungen des Klimawandels auf Österreichs Wälder - Entwicklung und vergleichende Evaluierung unterschiedlicher Prognosemodelle (WAMOD). Projektbericht 2010 (Antragsnummer A760631). BOKU, Institut für Waldbau; Institut für Waldwachstumsforschung, BFW, Institut für Waldinventur, Wien, Österreich (in German).
Wirth C, Schumacher J, Schultze ED (2004). Generic biomass functions for Norway spruce in Central Europe - a meta-analysis approach towards prediction and uncertainty estimation. Tree Physiology 24: 121-139.
Online | Gscholar
Wutzler T, Wirth C, Schumacher J (2008). Generic biomass functions for Common Beech (Fagus sylvatica L.) in Central Europe - predictions and components of uncertainty. Canadian Journal of Forest Research 38: 1661-1675.
CrossRef | Gscholar
Zaehle, S, Sitch, S, Prentice IC, Liski J, Cramer W, Erhard M, Hickler T, Smith B (2006). The importance of age-related decline in forest NPP for modelling regional carbon balances. Ecological Applications 16: 1555-1574.
CrossRef | Gscholar

Authors’ Affiliation

CS Eastaugh
H Hasenauer
Institute of Silviculture, University of Natural Resources and Life Sciences (BOKU), Vienna (Austria).

Corresponding author



Eastaugh CS, Hasenauer H (2011). Incorporating management history into forest growth modelling. iForest 4: 212-217. - doi: 10.3832/ifor0597-004

Paper history

Received: Nov 30, 2010
Accepted: Aug 22, 2011

First online: Nov 03, 2011
Publication Date: Nov 03, 2011
Publication Time: 2.43 months

© SISEF - The Italian Society of Silviculture and Forest Ecology 2011

  Open Access

This article is distributed under the terms of the Creative Commons Attribution-Non Commercial 4.0 International (https://creativecommons.org/licenses/by-nc/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

Creative Commons Licence

Breakdown by View Type

(Waiting for server response...)

Article Usage

Total Article Views: 49573
(from publication date up to now)

Breakdown by View Type
HTML Page Views: 42644
Abstract Page Views: 2272
PDF Downloads: 3553
Citation/Reference Downloads: 28
XML Downloads: 1076

Web Metrics
Days since publication: 4608
Overall contacts: 49573
Avg. contacts per week: 75.31

Article citations are based on data periodically collected from the Clarivate Web of Science web site
(last update: Feb 2023)

Total number of cites (since 2011): 8
Average cites per year: 0.62


Publication Metrics

by Dimensions ©

List of the papers citing this article based on CrossRef Cited-by.


iForest Similar Articles

iForest Database Search

Google Scholar Search

PubMed Search


This website uses cookies to ensure you get the best experience on our website. More info