- Open Access
Terrestrial vegetation redistribution and carbon balance under climate change
Carbon Balance and Management volume 1, Article number: 6 (2006)
Dynamic Global Vegetation Models (DGVMs) compute the terrestrial carbon balance as well as the transient spatial distribution of vegetation. We study two scenarios of moderate and strong climate change (2.9 K and 5.3 K temperature increase over present) to investigate the spatial redistribution of major vegetation types and their carbon balance in the year 2100.
The world's land vegetation will be more deciduous than at present, and contain about 125 billion tons of additional carbon. While a recession of the boreal forest is simulated in some areas, along with a general expansion to the north, we do not observe a reported collapse of the central Amazonian rain forest. Rather, a decrease of biomass and a change of vegetation type occurs in its northeastern part. The ability of the terrestrial biosphere to sequester carbon from the atmosphere declines strongly in the second half of the 21st century.
Climate change will cause widespread shifts in the distribution of major vegetation functional types on all continents by the year 2100.
The distribution of the world's vegetation has changed with past changes in climate and will continue to do so in the future. Due to rapidly increasing greenhouse gas concentrations, climate changes now more quickly than it has been doing for a long time  – but the pattern is irregular due to the complex changes in weather patterns, warming and rainfall change. So how much will vegetation change and where will it change most dramatically? Being able to answer these questions, even roughly, is important for two reasons. First, much of human well-being depends on ecosystems, due to the many services they provide . Second, land ecosystems contain large amounts of carbon which could be released as a consequence of major changes – they therefore may accelerate or slow down climate change substantially [3–5] (accelerate, e.g. due to increasing carbon emissions from organic soils, wildfires or forest die-back, or slow down, e.g. through increased vegetation growth and storage in dry or cold soils).
Terrestrial vegetation responds to climate change on several levels. Changes in temperature, precipitation, light and nutrient availability, and in atmospheric CO2 concentration influence plant biochemistry and physiology as well as the allocation of carbon to long- or short-lived plant parts such as leaves, stems and roots. Additionally, plants have evolved different functional strategies to cope with adverse conditions such as drought, cold or inundation (for example, the evergreen and the deciduous strategies of trees), therefore changes in these conditions eventually lead to changes in the species composition of an ecosystem – even if several decades may be needed for the process.
Mapping the outcome of the complex interplay of these processes has become possible due to the development of Dynamic Global Vegetation Models (DGVMs) , which simulate the terrestrial balances of carbon and water as well as the temporal development of vegetation in response to changing climate. The geographical pattern of vegetation emerges as a result of different responses of plant functional types to climate, with respect to productivity, bioclimatic constraints, access to resources and space, and sensitivity to natural disturbances such as fire.
Presently, the land biosphere is a net sink of carbon . Most simulations of the land biosphere's response to future climate change (as simulated by climate models) show a decline in this sink beginning around the middle of the 21st century [4, 8–10], with some scenarios even showing a net carbon loss by the end of the century [3, 9, 11]. The magnitude of this terrestrial feedback on climate is projected to be an additional increase in atmospheric CO2 concentration of between 20 and 200 ppm, implying an additional increase in temperature of between 0.1 K and 1.5 K .
Looking beyond these global numbers, all these simulations contain dramatic regional changes in vegetation structure and composition, in some cases of catastrophic extent. In some model simulations, for example, a collapse of parts of the Amazonian rain forest occurs by the year 2100 due to strongly decreasing rainfall . A decline in boreal forest area due to increasing heat stress on boreal trees has been reported . Another example is a transition from temperate savannah to subtropical woodland for a highland location in Africa . Changes such as these would imply a significant change in the composition and structure of the respective ecosystems – however, they differ depending on the greenhouse gas emission scenario, climate scenario and biosphere simulation model used.
As a step towards the identification of a more robust assessment, we use the state-of-the-art LPJ-DGVM [14, 15] and present results for change from the present in natural vegetation for two scenarios of climate change, selected to represent a wide range of potential futures from moderate (though not weak) to strong change by 2100. The consequences of moderate climate change (temperature increase over land: 2.9 K) were computed for ECHAM5 climate model projections under the SRES-B1 emission scenario (rising to 550 ppm CO2 in 2100). The consequences of strong climate change (temperature increase over land: 5.3 K) were computed for HadCM3 climate model projections under the SRES-A2 emission scenario (rising to 856 ppm CO2 in 2100).
Changes of vegetation patterns
Figure 1 shows simulated changes in the cover fraction of major vegetation types. Patterns of change are generally similar for moderate and for strong climate change but are more pronounced and extensive for stronger warming and atmospheric CO2 concentration.
Arctic tundra disappears in northern Eurasia, where warming is disproportionately strong, as deciduous woodlands expand. In northern Canada, new tundra is formed on the Arctic fringe that at present is only sparsely or not vegetated.
Boreal forests are affected by a number of large-scale changes. Evergreen vegetation increases its presence along the northern edge of the boreal zone. Further to the south, although in central western Eurasia and in Canada forests remain largely mixed, a widespread shift toward more deciduousness occurs. At the southern edge of the boreal forest, where it borders the steppe in central Asia and Canada, a recession of forest cover is observed in several areas due to increasing drought. In the stronger climate change scenario, the boreal forest collapses or changes into an open woodland in southern eastern Siberia, central western Siberia, and to the southwest of Hudson Bay in Canada due to increased mortality caused by heat stress. One mechanism by which this occurs is peak-summer heat stress on boreal species in areas where continentally cold winters prevent an invasion of temperate species.
In the temperate zone, evergreen plant functional types increase their fractions at the expense of deciduous vegetation in several regions, for example in the southeastern US, in Europe and in parts of eastern China. The savannahs and woodlands of South America and in southern Africa, on the other hand, show increased deciduousness.
The vegetation composition of the tropical evergreen zone is largely unaffected by climate change on the level of functional strategies. The simulations do not show a collapse of the Amazonian rain forest even for strong climate change. A decline in forest area is observed, however, along its eastern fringe.
In several regions with savannah, for example in southern Africa, grasslands experience woody encroachment. In many semi-arid regions, increased water use efficiency due to higher atmospheric CO2 concentrations leads to an increase in the herbaceous or grassy cover. Grass cover also increases where forests recede. Regional changes in precipitation, however, may lead to a net decline of non-woody vegetation cover in the affected regions. Increases and decreases of grass cover show a strongly differentiated spatial pattern around the globe.
Changes in terrestrial carbon storage and exchange
In these simulations, the biosphere contains 228 and 205 GtC (1 GtC = 1 billion tons of carbon) more carbon in the year 2100 than it did in the year 2000, for the moderate and the strong scenario, respectively. Vegetation biomass increases by 131 and 125 GtC, and the sum of carbon contained in litter and soils by 97 and 80 GtC. Net primary production increases by 18 and 28 GtC/yr, soil respiration by 17 and 24 GtC. Emissions from fires increase by 3.3 and 5.3 GtC/yr, respectively. As a consequence, the terrestrial biosphere's carbon sink is 2.4 and 1.6 GtC/yr smaller in 2100 than it is today (current sinks, 1971–2000: 3.0 and 1.9 GtC/yr; future sinks 2071–2100: 0.6 and 0.3 GtC/yr).
Figure 2 shows that these numbers are the cumulative result of a differentiated spatial pattern that is more pronounced for stronger climate change. The high northern latitudes act as a considerable sink of carbon as warming and carbon fertilisation stimulate vegetation growth. The boreal and temperate mid-latitudes lose carbon as biomass declines regionally. For strong climate warming, heterotrophic soil respiration is additionally stimulated. The African tropical and subtropical zone acts as a sink as carbon fertilisation and increased water use efficiency support more vegetation. In the strong climate change scenario, biomass strongly declines in the eastern reaches of the tropical rainforests of South America, in a region that is considerably larger than the region in which forest type changes substantially.
Figure 3 shows the temporal evolution of net carbon exchange between the terrestrial land surface and the atmosphere between 1900 and 2100 for the two scenarios. The terrestrial carbon sink currently achieved will increase until around 2025, then decline. At the end of the century, periods characterised by net terrestrial carbon release become more frequent and the average sink strength drops to close to zero. Major excursions toward the end of the century are related to periods with larger-than-average extremes in simulated climate. Exceptionally dry years with increased wildfire activity (carbon source) are first followed by re-growth during a period with slightly more precipitation (sink), then by renewed fire activity (source).
A major unknown in projecting future vegetation change is the disputed magnitude of carbon fertilisation effects on plant growth. In keeping with short-term experimental data [16, 17], DGVMs generally assume a diminishing but persistent stimulation of carbon assimilation with rising atmospheric carbon dioxide concentration and largely invariant allocation of carbon gains to plant compartments. Both of these assumptions remain under dispute in terms of their magnitude, persistence, and general applicability under nutrient and water constraints .
A source of uncertainty is also the spatial distribution of future changes in precipitation. Climate models differ considerably in the spatial pattern of their projections, but these changes are a first-order driver of biome structure. Uncertainty also remains due to gaps in generalised knowledge about the extent and mechanisms by which plants and ecosystems are able to adapt to change, and the thresholds beyond which compensatory mechanisms fail. The existence of such thresholds is well-known from the occasionally devastating regional effects on vegetation of years with anomalously extreme climate.
Future changes in nutrient limitations of vegetation growth, particularly with respect to nitrogen and phosphorous, may play an important role in some areas. The LPJ-DGVM assumes leaf nitrogen contents that optimise photosynthesis, but does not simulate allocation shifts away from leaves in response to nitrogen limitations. Simulated biomass increases are, however, not of a magnitude that would likely lead to general nutrient limitation on the global scale, though regional effects are to be expected. Increased water demand of vegetation at higher air temperatures is partly compensated in the LPJ-DGVM simulations by decreases in transpiration under higher ambient CO2 concentration. This effect is closely intertwined with concurrent shifts in vegetation abundance and composition . Projections of future fire frequency and severity in LPJ-DGVM simulations depend on the evolution of fuel load, fuel moisture and climatic fire probability, though fire occurrence in ecosystems under environmental stress may be influenced by additional factors.
In our simulations we do not observe a full-scale die-back of the Amazonian rain forest even under strong warming. The pattern of biomass reduction differs from that produced by the vegetation component of the HadCM3 simulation , but is similar to a pattern observed with a different DGVM . We suggest two explanations. First, HadCM3's simulated precipitation for the present is lower in some areas than observed. Our normalisation of climate model data to observed averages increases precipitation to observed levels, and the rainforest survives the subsequently simulated reduction, though biomass is lost in the eastern part of the region. Second, we do not compute feedbacks between climate and vegetation, which have been shown to enhance desiccation by 20% .
One of the mechanisms by which vegetation biomass is reduced in our simulations of Amazonia is an increase in fire frequency due to reduced precipitation and enhanced litter production. This causes a younger forest age with forest biomass reduced to nearly half of its present value. Where the forest manages to re-establish after fires, the biome type remains unchanged. Forest is replaced by non-woody vegetation in those regions where the frequency and severity of fires prevents its re-establishment. Carbon emissions from fire increase most strongly in northwestern South America, large parts of southern Africa, the Sahel, Australia and India.
The simulated northward expansion of the boreal forest is generally accepted, and there is evidence for a recession of the southern boreal drought-induced tree line . Little substantial discussion exists of the potential causes of boreal forest die-back due to heat stress, possibly due to high peak tissue temperatures. Mortality of stressed boreal trees is increased due to insect infestation and subsequent fire. However, overall vegetation dynamics are complicated by changes in permafrost thawing depth and snow cover, and by potential effects of carbon fertilisation.
The spatial pattern of the world's vegetation will change under climate change [22–24]. But by how much and where? The answer given by current research has to remain tentative, but all indications are that the changes will be wide-spread if atmospheric CO2 increase and warming are not limited to a small magnitude. It would be unwise to treat these changes lightly. Shifts in the spatial distribution of vegetation types signal severe change in the underlying ecosystems, with effects on a large number of species. The modelling, observational and political challenges that follow for biosphere and ecosystem research are as large, if not larger, than those encountered in climate system research. The biosphere will continue to remove anthropogenic carbon from the atmosphere but with diminishing strength in the second half of the 21st century. Large regions of the globe can be identified that are projected to be either long-term sinks or sources of carbon.
Vegetation models of the next generation will be challenged to provide an improved basis for determining the impact of climate change on vegetation distribution, structure and properties. The processes to be modelled, their physiological basis and interlinkages on the level of ecosystems will require review, discussion in terms of relevance and priority, and integration into new model versions. Close links between biosphere modellers, field ecologists, plant physiologists, biogeographers and experimentalists in biogeochemistry will be required for achieving this advance .
Climate change is an important but by far not the only pressure on the terrestrial biosphere. Still expanding land use and deforestation as well as chemical pollution, interregional exchange of species and, perhaps, more widespread genetic engineering in the future will in all likelihood rival or even overshadow the effects of climate change. Today's tropical forests, for example, are under pressure mainly from socioeconomic drivers that are unrelated to but will interact with climate change.
The implications of these findings for policy are threefold. First, biome change is an important dimension in which the degree of acceptable climate change has to be evaluated. Second, the saturation and subsequent decline of terrestrial carbon storage is a factor to be taken into account in discussions of mitigation and adaptation policies. Third, a substantial focus on climatic and non-climatic causes of ecosystem change is required in the next several years to advance the reliability of available projections of land biosphere change.
The LPJ-DGVM [14, 15] is a process-based model of key ecosystem processes that govern terrestrial biogeochemistry and biogeography. A Farquhar-Collatz photosynthesis scheme is coupled to a two-layered soil hydrological scheme for daily simulations of gross primary production and plant respiration, including effects of drought stress on assimilation and evapotranspiration, depending on tissue-specific C:N ratios, biomass and phenology, and using a modified Arrhenius temperature-dependent formulation. Assimilated carbon is allocated annually to four pools (leaves, sapwood, heartwood and fine roots) to satisfy a set of allometric and functional relations. Leaf and root turnover, as well as plant mortality, feed a litter pool, a slow and a fast soil carbon pool that decay depending on soil temperature through a modified Arrhenius formulation, and soil moisture. Vegetation functional differences are represented by seven woody and two herbaceous plant functional types (PFTs) differentiated by their physiological, physiognomic and phenological attributes. These may co-exist at any location, depending on plant competition for resources and space, and a set of environmental constraints. Fire disturbance is simulated as a function of a threshold litter load and surface soil moisture, permafrost was not modelled. A selection of model validation is contained in [14, 19, 26–28].
Climate model data were interpolated to 0.5 degrees resolution and normalised in each pixel to the observed 1961–90 CRU2002 climatology  mean value. For Echam5 , the global mean temperature over land remained nearly unchanged as a consequence, while for HadCM3  it was increased by about 1 K. Individual pixels will, however, have experienced different amounts of normalisation. Precipitation was normalised using a relative factor. A spinup of LPJ was performed lasting 900 years using CRU data to place carbon pools into a realistic initial state. This was followed by a second shorter spin-up with normalised GCM data to allow adjustments of vegetation distribution to the characteristics of GCM climate. Transient runs from 1860 to 2100 followed, of which the averages of the periods 1971–2000 and 2071–2100 were evaluated. The two combinations of scenario and model used were selected to represent a moderate (but not weak) and a strong example of potential futures from among the range of results produced by a larger number of combinations of scenarios and models.
Karl TR, Trenbreth KE: Modern global climate change. Science 2003, 302: 1719–1723. 10.1126/science.1090228
Alcamo J, van Vuuren D, Ringler DC, Cramer W, Masui T, Alder J, Schulze K: Changes in nature's balance sheet: model-based estimates of future worldwide ecosystem services. Ecology and Society 2005, 10: 19.
Cox PM, Betts RA, Jones CD, Spall SA, Totterdell IJ: Acceleration of global warming due to carbon-cycle feedbacks in a coupled climate model. Nature 2000, 408: 184–187. 10.1038/35041539
Cramer W, Bondeau A, Woodward FI, Prentice IC, Betts RA, Brovkin V, Cox PM, Fisher V, Foley JA, Friend AD, Kucharik C, Lomas MR, Ramankutty N, Sitch S, Smith B, White A, Young-Molling C: Global response of terrestrial ecosystem structure and function to CO2 and climate change: Results from six dynamic global vegetation models. Global Change Biology 2001, 357–373. 10.1046/j.1365-2486.2001.00383.x
Friedlingstein P, Dufresne J-L, Cox PM, Rayner P: How positive is the feedback between climate change and the carbon cycle? Tellus 2003, 55B: 692–700.
Prentice IC, Bondeau A, Cramer W, Harrison SP, Hickler T, Lucht W, Sitch S, Smith B, Sykes MT: Dynamic global vegetation modelling: quantifying terrestrial ecosystem responses to large-scale environmental change. In Terrestrial Ecosystems in a Changing World. Edited by: Canadell JD, Pataki E, Pitelka LF. Berlin, The IGBP Series, Springer-Verlag; 2006.
House JI, Prentice IC, Ramankutty N, Houghton RA, Heimann M: Reconciling apparent inconsistencies in estimates of terrestrial CO2 sources and sinks. Tellus 2003, 55B: 345–363.
Berthelot M, Friedlingstein P, Ciais P, Dufresne J-L, Monfray P: How uncertainties in future climate change predictions translate into future terrestrial carbon fluxes. Glob Change Biol 2005, 1: 959–970. 10.1111/j.1365-2486.2005.00957.x
Schaphoff S, Lucht W, Gerten D, Sitch S, Cramer W, Prentice IC: Terrestrial biosphere carbon storage under alternative climate projections. Climatic Change 2006. doi: 10.1007/s10584–005–9002–5. doi: 10.1007/s10584-005-9002-5.
Woodward FI, Lomas MR: Vegetation dynamics – simulating responses to climatic change. Biol Rev 2004, 79: 643–670. 10.1017/S1464793103006419
Joos F, Prentice IC, Sitch S, Meyer R, Hooss G, Plattner G-K, Gerber S, Hasselmann K: Global warming feedbacks on terrestrial carbon uptake under the Intergovernmental Panel on Climate Change (IPCC) emission scenarios. Global Biogeochemical Cycles 2001, 15: 891–907. 10.1029/2000GB001375
Friedlingstein P, Cox P, Betts R, Bopp L, von Bloh W, Brovkin V, Cadule P, Doney S, Eby M, Fung I, Bala G, John J, Jones C, Joos F, Kato T, Kawamiya M, Knorr W, Lindsay K, Matthews HD, Raddatz T, Rayner P, Reick C, Roeckner E, Snitzler K-G, Schnur R, Strassmann K, Weaver AJ, Yoshikawa C, Zeng N: Climate-carbon cycle feedback analysis: Results from the C4MIP model intercomparison. J Clim 2006, 19: 3337–3353. 10.1175/JCLI3800.1
Cox PM, Betts RA, Collins M, Harris PP, Huntingford C, Jones CD: Amazonian forest dieback under climate-carbon cycle projections for the 21st century. Theoretical and Applied Climatology 2004, 78: 137–156. 10.1007/s00704-004-0049-4
Sitch S, Smith B, Prentice IC, Arneth A, Bondeau A, Cramer W, Kaplan JO, Levis S, Lucht W, Sykes MT, Thonicke K, Venevsky S: Evaluation of ecosystem dynamics, plant geography and terrestrial carbon cycling in the LPJ dynamic global vegetation model. Global Change Biol 2003, 9: 161–185. 10.1046/j.1365-2486.2003.00569.x
Gerten D, Schaphoff S, Haberlandt U, Lucht W, Sitch S: Terrestrial vegetation and water balance: Hydrological evaluation of a dynamic global vegetation model. Journal of Hydrology 2004, 286: 249–270. 10.1016/j.jhydrol.2003.09.029
Norby RJ, Luo YQ: Evaluating ecosystem responses to rising atmospheric CO2 and global warming in a multi-factor world. New Phytologist 2004, 162: 281–293. 10.1111/j.1469-8137.2004.01047.x
Ellsworth DS, Reich PB, Naumburg ES, Koch GW, Kubiske ME, Smith SD: Photosynthesis, carboxylation and leaf nitrogen responses of 16 species to elevated pCO2 across four free-air CO2 enrichment experiments in forest, grassland and desert. Global Change Biology 2004, 10: 2121–2138. 10.1111/j.1365-2486.2004.00867.x
Körner C, Asshoff R, Bignucolo O, Hättenschwiler S, Keel SG, Peléz-Riedl S, Pepin S, Siegwolf RTW, Zotz G: Carbon flux and growth in mature deciduous forest trees exposed to elevated CO2. Science 2005, 309: 1360–1362. 10.1126/science.1113977
Gerten D, Lucht W, Schaphoff S, Cramer W, Hickler T, Wagner W: Hydrologic resilience of the terrestrial biosphere. Geophys Res Lett 2005., 32: doi: 10.1029/2005GL024247. doi: 10.1029/2005GL024247.
Betts RA, Cox PM, Harris PP, Huntingford C, Jones CD: The role of ecosystem-atmosphere interactions in simulated Amazon forest dieback under global climate warming. Theoretical and Applied Climatology 2004, 78: 157–175. 10.1007/s00704-004-0050-y
Hogg EH, Schwarz AG: Regeneration of planted conifers across climatic moisture gradients on the Canadian prairies: implications for distribution and climate change. Journal of Biogeography 1997, 24: 527–534. 10.1111/j.1365-2699.1997.00138.x
Neilson RP, Drapek RJ: Potentially complex biosphere responses to transient global warming. Global Change Biology 1998, 4: 505–521. 10.1046/j.1365-2486.1998.t01-1-00202.x
Leemans R, Eichhout B: Another reason for concern: regional and global impacts on ecosystems for different levels of climate change. Global Environmental Change 2004, 14: 219–228. 10.1016/j.gloenvcha.2004.04.009
Gerber S, Joos F, Prentice IC: Sensitivity of a dynamic global vegetation model to climate and atmospheric CO 2 . Global Change Biology 2004, 10: 1223–1239. 10.1111/j.1529-8817.2003.00807.x
Moorcroft PR: How close are we to a predictive science of the biosphere? Trends in Ecology and Evolution 2006, 21: 400–407. 10.1016/j.tree.2006.04.009
Peylin P, Bousquet P, Le Quéré C, Sitch S, Friedlingstein P, McKinley G, Gruber N, Rayner P, Ciais P: Multiple constraints on regional CO2 flux variations over land and oceans. Global Biogeochemical Cycles 2005., 19: doi:10.1029/2003GB002214. doi:10.1029/2003GB002214.
Lucht W, Prentice IC, Myneni RB, Sitch S, Friedlingstein P, Cramer W, Bousquet P, Buermann W, Smith B: Climatic control of the high-latitude vegetationgreening trend and Pinatubo effect. Science 2002, 296: 1687–1689. 10.1126/science.1071828
Beer C, Lucht W, Schmullius C, Shvidenko A: Small net carbon dioxide uptake by Russian forests during 1981–1999. Geophys Res Lett, in press.
New MG, Hulme M, Jones PD: Representing twentieth-century space-time climate variability, Part II, Development of 1901–1996 monthly grids of terrestrial surface climate. Journal of Climate 2000, 13: 2217–2238. Publisher Full Text 10.1175/1520-0442(2000)013<2217:RTCSTC>2.0.CO;2
Roeckner E, Bauml G, Bonaventura L, Brokopf R, et al.: The general circulation model ECHAM5. Part I: Model description. Hamburg, Max-Planck-Institute for Meteorology; 2003. Report 349.
Gordon C, Cooper C, Senior CA, Banks HT, Gregory JM, Johns TC, Mitchell JFB, Wood RA: The simulation of SST, sea ice extents and ocean heat transports in a version of the Hadley Centre coupled model without flux adjustments. Climate Dynamics 2000, 16: 147–168. 10.1007/s003820050010
We acknowledge project funding from the Pakt für Forschung of the Leibniz Association (SS, UH) and by the State of Brandenburg (TE). We thank Werner von Bloh, Christoph Müller, Stefanie Jachner and Alberte Bondeau for valuable contributions, and Andreas Fischlin, Guy Midgley, Ron Neilson, Colin Prentice and Stephen Sitch for discussions. We acknowledge the international modeling groups for providing IPCC climate model scenario simulations for analysis, the Program for Climate Model Diagnosis and Intercomparison (PCMDI) for collecting and archiving the model data, the JSC/CLIVAR Working Group on Coupled Modelling (WGCM) and their Coupled Model Intercomparison Project (CMIP) and Climate Simulation Panel for organizing the model data analysis activity, and the IPCC WG1 TSU for technical support. The IPCC Data Archive at Lawrence Livermore National Laboratory is supported by the Office of Science, U.S. Department of Energy. We thank Will Steffen for having been the Handling Editor and three anonymous reviewers for their excellent comments.
The author(s) declare that they have no competing interests.
WL designed and supervised the analysis and wrote the manuscript. SS carried out the analysis and co-devised the design. TE and UH contributed data handling and analysis, WC analysis, writing and supervision.
Wolfgang Lucht, Sibyll Schaphoff contributed equally to this work.