It is shown that by calibrating the simulation model BIOME-BGC with mandatory and optional Level II data, within the ICP Forest programme, a well-founded calculation of the carbon budget of forest stands is achievable and, based on succeeded calibration, the modified BIOME-BGC model is a useful tool to assess the effect of climate change on forest ecosystems.
The main focus of the intensive monitoring of forests (Level II) within the ICP Forest programme of the UN-ECE is the investigation of effects of atmospheric deposition on element budgets, soil chemistry, crown condition, and tree growth of forest ecosystems. Investigations on forest carbon budget are limited to the sampling of stem increments within a 5 years interval. In a limited number of stands some supplementary measurements are carried out concerning the above and below ground carbon turnover.
The aim of this study was the assessment of additional measurements (litter fall, permanent girth measurement, tree ring analyses, soil respiration, long term development of soil carbon, soil physics, soil moisture) to improve the initialization, calibration and validation of dynamic simulation models with respect to the carbon budget of forest ecosystems. The validated simulation models can be applied for the assessment of the impact of climate change on the carbon budget of forests.
Methods
Sixteen German Level II plots with Fagus sylvatica (3), Pinus sylvestris (8) and Picea abies (5) were selected, where the measurements of the mandatory and optional programmes (meteorology, stand precipitation, soil moisture, deposition, soil solution, soil chemistry, stem biomass (5 years cycle), stem girth measurements, litter fall, crown condition, foliar analysis) and in some cases data on soil physics, tree ring analyses, and sap flow were available. At these plots soil respiration was measured, additionally (Badeck et al. 2007). For differentiation between heterotrophic and root respiration the trenching technique was used (Hanson et al. 2000).
In our investigation the data were used for the initialization and calibration of BIOME-BGC version 4.2 (Thornton et al. 2002 - http://www.ntsg.umt.edu/models/bgc/), a simulation model for calculating the energy, water, carbon, and nitrogen budget of biomes. For the purpose of the project the programme code of BIOME-BGC was extended by a forest management module, a multi-layered soil water, C, and temperature module, and a root growth module. The programme code concerning phenology, decomposition, canopy evaporation and transpiration modules was modified. At the present status the model shows a more species and stand specific behavior. The number of ecophysiological model parameters (including management) was enhanced from 37 to 80.
The carbon and water budgets of these 16 Level-II sites were calculated for the simulation period 1996-2005 using meteorological data of the open field measurements of the plots.
For analyzing the effects of climate change, the climate scenarios of the statistical model STAR (Orlowsky 2007) for the DWD stations and the years 2046-2055, based on the SRES emission scenario A2 of the IPCC and the global climate model ECHAM5-MPI-OM were used. For the base and climate change scenarios identical initial model values were used. The climate change scenarios consider the predicted atmospheric CO2 concentrations, too, but did not use changing nitrogen deposition.
Results
Comparisons of simulation results with measurements on stand precipitation, soil moisture, and sap flow demonstrate, that the model BIOME-BGC (including modifications) allows a precise simulation of water budget (not depicted). The same applies to soil temperature.
The long term averages of leaf litter fall and stem growth can be simulated with satisfying precision, if the plot number 906 is disregarded (Fig. 1, Fig. 2). In some stands basal area increment from tree ring analyses are weakly correlated to simulated stem increment, mainly in pine stands on sandy soils of Brandenburg (Fig. 3). The annual variation of measured values of leaf litter fall and stem growth exceeds the annual variation of simulation results.
Repeated measurements of soil carbon are only available for two stands (Solling beech and spruce). In these stands, the increasing soil carbon stocks cannot be found in the current simulation results. Instead, a slow decrease is simulated (not depicted).
The seasonality of soil respiration (heterotrophic + root resp.) can be reproduced by the model predictions, with the exception of some peaks in summer (Fig. 4). Major deviations occur if measurements at soil sites with trenched roots are interpreted as heterotrophic respiration and compared to simulation results (not depicted).
Simulations using climate change scenarios showed an increasing productivity (GPP, NPP, stem and leaf growth - Fig. 5). On the other hand, the plant and heterotrophic respiration is rising, too. The carbon sequestration of the stands (NEP) is increased under the climate change scenarios (Fig. 6). Due to reduced precipitation, a weak reduction of evaporation and a considerable reduction of soil water outflow (43 % of base scenario) are predicted (not depicted). Transpiration reacts different depending on tree species and site conditions.
Discussion
Correct simulations of hydrology and the soil temperature of the forest ecosystem, as realized in this investigation, are essential preconditions for estimating the carbon budget of forest ecosystems with dynamic simulation models.
Even if simulated leaf litter fall and stem growth more or less agree with measurements, the higher annual variations of observations compared to those of simulations indicates, that the model does not yet consider all essential factors determining stem and leaf growth (e.g., fructification, air pollutants). In case of the plot number 906 successful model calibration is not possible. This is may caused by data errors on driving forces, initial values, or on data for calibration.
The divergent simulation results on development of soil carbon pools compared to measurements may have different causes. First, the spatial heterogeneity of soil carbon leads to a wide range of data (Lloyd & McKee 1983, Klinck et al. 2008) that complicate model initialization and calibration. Additionally, measurements on coarse woody debris are lacking, and data on root turnover, one of the main input-path of organic matter into the soil, are uncertain. Finally, seed production as a further pathway of organic matter into the litter layer is neglected by the model.
With respect to these uncertainties the additional soil respiration measurements represent valuable data for model initialization and calibration. Measured soil respiration with living roots could be simulated satisfactorily. In contrast, the root trenching technique for separation of root and heterotrophic respiration may produce some artifacts, at least under dry summer conditions. Thus, the latter measurements are an uncertain database for model calibration aiming on separated consideration of root and heterotrophic respiration.
The climate change scenarios applied in this investigation are characterized by rising temperature and atmospheric CO2 concentrations and decreasing precipitation. The application of these climate change scenarios are connected to some artifacts caused by systematic climate differences between the meteorological stations of the national weather service (used for climate change scenarios) and open land stations belonged to Level II plots (used for base scenarios).
The simulation results on climate change effects on water budget with decreasing evaporation and water outflow are in line with other investigations (Gerstengarbe et al. 2003). The results of higher wood production under climate change conditions are in agreement with analyses on growth trends of the past (Spiecker 1999, Boisvenue & Running 2006) and simulation studies for future development (Lasch et al. 2002, Eggers et al. 2008). Obviously, the simulated rising water stress at some Level II plots is counteracted by the fertilizing effect of rising CO2 concentration. By contrast, a simulation study on pine stands in Brandenburg (Germany) reports reduced productivity (Lasch & Suckow 2007). In Europe different gradients (north-south, altitude, maritime-continental) of forest growth reaction on climate change are expected. In the boreal zone the limiting low temperature effect will be reduced, whereas in the Mediterranean zone the limiting drought effect is expected to be enhanced (Kellomäki et al. 2005, De Vries et al. 2007).
Conclusions
The monitoring data from the current mandatory and optional Level II programme provide valuable information for calibration of dynamic simulation models for calculating the carbon budget of forest ecosystems. Additionally, the data offer some hints for further model development. By calibrating the simulation model BIOME-BGC with mandatory and optional Level II data, a well-founded calculation of the carbon budget of forest stands is achievable. Furthermore, based on succeeded calibration, the modified BIOME-BGC model is a useful tool to assess climate change effects on forest ecosystems.
Acknowledgements
These results are part of a project that was financially supported by the European Commission, Contract No. Forest Focus - DE 2003/2004 BB 5, DE 2003/2004 BY 4, DE 2003/2004 NI 6.
ReferencesBadeck FW, Beese F, Berthold D, Einert P, Jochheim H, Kallweit R, Konopatzky A, Lasch P, Meesenburg H, Meiwes KJ, Puhlmann M, Raspe S, Schulte-Bisping H, Schulz C, Suckow FParametrisierung, Kalibrierung und Validierung von Modellen des Kohlenstoffumsatzes in Waldökosystemen und deren Böden. Abschlussbericht Forest Focus C2-Projekte DE 2003/2004 BB 5, DE 2003/2004 BY 4, DE 2003/2004 NI 6, pp. 110.2007Boisvenue C, Running SWImpacts of climate change on natural forest productivity - evidence since the middle of the 20th century. Global Change Biology 12: 862-8822006De Vries W, Wamelink GWW, Reinds GJ, Wieggers HJJ, Mol-Dijkstra JP, Kros H, Nabuurs GJ, Pus-sinen A, Solberg B, Dobbertin M, Laubhann D, Sterba H, Van Oijen MAssessment of the relative importance of nitrogen deposition, climate change and forest management on the sequestration of carbon by forests in Europe. Alterra Rapport 1538: pp. 302.2007Eggers J, Lindner M, Zudin S, Zaehle S, Liski JImpact of changing wood demand, climate and land use on European forest resources and carbon stocks during the 21st century. Global Change Biology 14: 2288-2303.2008Gerstengarbe FW, Badeck F, Hattermann F, Krysanova V, Lahmer W, Lasch P, Stock M, Suckow F, Wechsung F, Werner PCStudie zur klimatischen Entwicklung im Land Brandenburg bis 2055 und deren Auswirkungen auf den Wasserhaushalt, die Forst- und Landwirtschaft sowie die Ableitung erster Perspektiven. PIK Report 83: pp. 77.2003Hanson PJ, Edwards NT, Garten CT, Andrews JASeparating root and soil microbial contributions to soil respiration: A review of methods and observations. Biogeochemistry 48: 115-146.2000Kellomäki S, Peltola H, Bauwens B, Dekker M, Mohren F, Badeck F-W, Gracia C, Sánchez A, Pla E, Sabaté SEuropean mitigation and adaptation potentials: conclusions and recommendation. In: “Management of European forests under changing climatic conditions” (Kellomäki S, Leinonen S eds). Final Report of the Project: “Silvicultural response strategies to climatic change in management of european forests”, pp. 402-427.2005Klinck U, Fröhlich D, Meiwes KJZur Problematik der Quantifizierung von C- und N-Gehalten, C/N-Verhältnissen und C- und N-Vorräten der Humusauflage und des mineral-ischen Oberbodens. Allgemeine Forst und Jagdzeitung 179: 1-7.2008Lasch P, Suckow FReaktion von Kiefernbeständen unter Klimaänderungen - eine Analyse mit dem Waldwachstumsmodell 4C. In: “Kiefer im nordostdeutschen Tiefland - Ökologie und Bewirtschaftung” (Kätzel R, Möller K, Löffler S, Engel J, Liero K eds). Eberswalder Forstliche Schriftenreihe 32: 230-237.2007Lasch P, Badeck F-W, Lindner M, Suckow FSensitivity of simulated forest growth to changes in climate and atmospheric CO2. Forstwissenschaftliches Centralblatt 121 (Suppl. 1): 155-171.2002Lloyd FT, McKee WHReplications and subsamples needed to show treatment responses on forest soils of the coastal plain. Soil Science Society of America Journal 47: 587-590.1983Orlowsky BSetzkasten Vergangenheit - ein kombinatorischer Ansatz für regionale Klimasimulationen. PhD thesis, Department of Geoscience, University Hamburg, Germany, pp. 187.2007Spiecker HOverview of recent growth trends in European forests. Water Air and Soil Pollution 116: 33-46.1999Thornton PE, Law BE, Gholz HL, Clark KL, Falge E, Ellsworth DS, Golstein AH, Monson RK, Hollinger D, Falk M, Chen J, Sparks JPModeling and measuring the effects of disturbance history and climate on carbon and water budgets in evergreen needleleaf forests. Agricultural and Forest Meteorology 113: 185-222.2002
Simulation results compared to measurements: litter fall.
Simulation results compared to measurements: wood increment.
Simulation results compared to measurements: deviation of the measured basal area increment from hyperbolic trend compared to the deviation of the simulated current annual stem increment from the linear trend.