Search

The leaf economics spectrum 1,2 and the global spectrum of plant forms and functions 3 revealed fundamental axes of variation in plant traits, which represent different ecological strategies that are shaped by the evolutionary development of plant species 2 . Ecosystem functions depend on environmental conditions and the traits of species that comprise the ecological communities 4 . However, the axes of variation of ecosystem functions are largely unknown, which limits our understanding of how ecosystems respond as a whole to anthropogenic drivers, climate and environmental variability 4,5 . Here we derive a set of ecosystem functions 6 from a dataset of surface gas exchange measurements across major terrestrial biomes. We find that most of the variability within ecosystem functions (71.8%) is captured by three key axes. The first axis reflects maximum ecosystem productivity and is mostly explained by vegetation structure. The second axis reflects ecosystem water-use strategies and is jointly explained by variation in vegetation height and climate. The third axis, which represents ecosystem carbon-use efficiency, features a gradient related to aridity, and is explained primarily by variation in vegetation structure. We show that two state-of-the-art land surface models reproduce the first and most important axis of ecosystem functions. However, the models tend to simulate more strongly correlated functions than those observed, which limits their ability to accurately predict the full range of responses to environmental changes in carbon, water and energy cycling in terrestrial ecosystems 7,8 ., (© 2021. The Author(s).)

Leaf dark respiration (Rdark ) is an important yet poorly quantified component of the global carbon cycle. Given this, we analyzed a new global database of Rdark and associated leaf traits. Data for 899 species were compiled from 100 sites (from the Arctic to the tropics). Several woody and nonwoody plant functional types (PFTs) were represented. Mixed-effects models were used to disentangle sources of variation in Rdark . Area-based Rdark at the prevailing average daily growth temperature (T) of each site increased only twofold from the Arctic to the tropics, despite a 20°C increase in growing T (8-28°C). By contrast, Rdark at a standard T (25°C, Rdark (25) ) was threefold higher in the Arctic than in the tropics, and twofold higher at arid than at mesic sites. Species and PFTs at cold sites exhibited higher Rdark (25) at a given photosynthetic capacity (Vcmax (25) ) or leaf nitrogen concentration ([N]) than species at warmer sites. Rdark (25) values at any given Vcmax (25) or [N] were higher in herbs than in woody plants. The results highlight variation in Rdark among species and across global gradients in T and aridity. In addition to their ecological significance, the results provide a framework for improving representation of Rdark in terrestrial biosphere models (TBMs) and associated land-surface components of Earth system models (ESMs)., (© 2015 The Authors. New Phytologist © 2015 New Phytologist Trust.)

We describe the Max Planck Institute Carbon Cycle Data Assimilation System (MPI-CCDAS) built around the tangent-linear version of the land surface scheme of the MPI-Earth System Model v1 (JSBACH). The simulated terrestrial biosphere processes (phenology and carbon balance) were constrained by observations of the fraction of photosynthetically active radiation (TIP-FAPAR product) and by observations of atmospheric CO2 at a global set of monitoring stations for the years 2005-2009. The system successfully, and computationally efficiently, improved average foliar area and northern extra-tropical seasonality of foliar area when constrained by TIP-FAPAR. Global net and gross carbon fluxes were improved when constrained by atmospheric CO2, although the system tended to underestimate tropical productivity. Assimilating both data streams jointly allowed the MPI-CCDAS to match both observations (TIP-FAPAR and atmospheric CO2) equally well as the single data stream assimilation cases, therefore overall increasing the appropriateness of the resultant parameter values and biosphere dynamics. Our study thus highlights the role of the TIP-FAPAR product in stabilising the underdetermined atmospheric inversion problem and demonstrates the value of multiple-data stream assimilation for the simulation of terrestrial biosphere dynamics. The constraint on regional gross and net CO2 flux patterns is limited through the parametrisation of the biosphere model. We expect improvement on that aspect through a refined initialisation strategy and inclusion of further biosphere observations as constraints. [ABSTRACT FROM AUTHOR]

Since 70% of global forests are managed and forests impact the global carbon cycle and the energy exchange with the overlying atmosphere, forest management has the potential to mitigate climate change. Yet, none of the land-surface models used in Earth system models, and therefore none of today's predictions of future climate, accounts for the interactions between climate and forest management. We addressed this gap in modelling capability by developing and parametrising a version of the ORCHIDEE land-surface model to simulate the biogeochemical and biophysical effects of forest management. The most significant changes between the new branch called ORCHIDEE-CAN (SVN r2290) and the trunk version of ORCHIDEE (SVN r2243) are the allometric-based allocation of carbon to leaf, root, wood, fruit and reserve pools; the transmittance, absorbance and reflectance of radiation within the canopy; and the vertical discretisation of the energy budget calculations. In addition, conceptual changes were introduced towards a better process representation for the interaction of radiation with snow, the hydraulic architecture of plants, the representation of forest management and a numerical solution for the photosynthesis formalism of Farquhar, von Caemmerer and Berry. For consistency reasons, these changes were extensively linked throughout the code. Parametrisation was revisited after introducing 12 new parameter sets that represent specific tree species or genera rather than a group of often distantly related or even unrelated species, as is the case in widely used plant functional types. Performance of the new model was compared against the trunk and validated against independent spatially explicit data for basal area, tree height, canopy structure, gross primary production (GPP), albedo and evapotranspiration over Europe. For all tested variables, ORCHIDEE-CAN outperformed the trunk regarding its ability to reproduce large-scale spatial patterns as well as their inter-annual variability over Europe. Depending on the data stream, ORCHIDEE-CAN had a 67 to 92% chance to reproduce the spatial and temporal variability of the validation data. [ABSTRACT FROM AUTHOR]

Terrestrial productivity in semi-arid woodlands is strongly susceptible to changes in precipitation, and semi- arid woodlands constitute an important element of the global water and carbon cycles. Here, we use the Carbon Cycle Data Assimilation System (CCDAS) to investigate the key parameters controlling ecological and hydrological activities for a semi-arid savanna woodland site in Maun, Botswana. Twenty-four eco-hydrological process parameters of a terres-trial ecosystem model are optimized against two data streams separately and simultaneously: daily averaged latent heat flux (LHF) derived from eddy covariance measurements, and decadal fraction of absorbed photosynthetically active radiation (FAPAR) derived from the Sea-viewing Wide Field-of- view Sensor (SeaWiFS). Assimilation of both data streams LHF and FAPAR for the years 2000 and 2001 leads to improved agreement between measured and simulated quantities not only for LHF and FAPAR, but also for photosynthetic CO2 uptake. The mean uncertainty reduction (relative to the prior) over all parameters is 14.9% for the simultaneous assimilation of LHF and FAPAR, 8.5% for assimilating LHF only, and 6.1% for assimilating FAPAR only. The set of parameters with the highest uncertainty reduction is similar between assimilating only FAPAR or only LHF. The highest uncertainty reduction for all three cases is found for a parameter quantifying max- imum plant-available soil moisture. This indicates that not only LHF but also satellite-derived FAPAR data can be used to constrain and indirectly observe hydrological quantities. [ABSTRACT FROM AUTHOR]

The vertical distribution of soil organic matter (SOM) in the profile may constitute an important factor for soil carbon cycling. However, the formation of the SOM profile is currently poorly understood due to equifinality, caused by the entanglement of several processes: input from roots, mixing due to bioturbation, and organic matter leaching. In this study we quantified the contribution of these three processes using Bayesian parameter estimation for the mechanistic SOM profile model SOMPROF. Based on organic carbon measurements, 13 parameters related to decomposition and transport of organic matter were estimated for two temperate forest soils: an Arenosol with a mor humus form (Loobos, the Netherlands), and a Cambisol with mull-type humus (Hainich, Germany). Furthermore, the use of the radioisotope 210Pbex as tracer for vertical SOM transport was studied. For Loobos, the calibration results demonstrate the importance of organic matter transport with the liquid phase for shaping the vertical SOM profile, while the effects of bioturbation are generally negligible. These results are in good agreement with expectations given in situ conditions. For Hainich, the calibration offered three distinct explanations for the observations (three modes in the posterior distribution). With the addition of 210Pbex data and prior knowledge, as well as additional information about in situ conditions, we were able to identify the most likely explanation, which indicated that root litter input is a dominant process for the SOM profile. For both sites the organic matter appears to comprise mainly adsorbed but potentially leachable material, pointing to the importance of organo-mineral interactions. Furthermore, organic matter in the mineral soil appears to be mainly derived from root litter, supporting previous studies that highlighted the importance of root input for soil carbon sequestration. The 210Pbex measurements added only slight additional constraint on the estimated parameters. However, with sufficient replicate measurements and possibly in combination with other tracers, this isotope may still hold value as tracer for SOM transport. [ABSTRACT FROM AUTHOR]

Terrestrial carbon (C) cycle models applied for climate projections simulate a strong increase in net primary productivity (NPP) due to elevated atmospheric CO2 concentration during the 21st century. These models usually neglect the limited availability of nitrogen (N) and phosphorus (P), nutrients that commonly limit plant growth and soil carbon turnover. To investigate how the projected C sequestration is altered when stoichiometric constraints on C cycling are considered, we incorporated a P cycle into the land surface model JSBACH (Jena Scheme for Biosphere-Atmosphere Coupling in Hamburg), which already includes representations of coupled C and N cycles. The model reveals a distinct geographic pattern of P and N limitation. Under the SRES (Special Report on Emissions Scenarios) A1B scenario, the accumulated land C uptake between 1860 and 2100 is 13% (particularly at high latitudes) and 16% (particularly at low latitudes) lower in simulations with N and P cycling, respectively, than in simulations without nutrient cycles. The combined effect of both nutrients reduces land C uptake by 25% compared to simulations without N or P cycling. Nutrient limitation in general may be biased by the model simplicity, but the ranking of limitations is robust against the parameterization and the inflexibility of stoichiometry. After 2100, increased temperature and high CO2 concentration cause a shift from N to P limitation at high latitudes, while nutrient limitation in the tropics declines. The increase in P limitation at high-latitudes is induced by a strong increase in NPP and the low P sorption capacity of soils, while a decline in tropical NPP due to high autotrophic respiration rates alleviates N and P limitations. The quantification of P limitation remains challenging. The poorly constrained processes of soil P sorption and biochemical mineralization are identified as the main uncertainties in the strength of P limitation. Even so, our findings indicate that global land C uptake in the 21st century is likely overestimated in models that neglect P and N limitations. In the long term, insufficient P availability might become an important constraint on C cycling at high latitudes. Accordingly, we argue that the P cycle must be included in global models used for C cycle projections. [ABSTRACT FROM AUTHOR]

The vertical distribution of soil organic matter (SOM) in the profile may constitute a significant factor for soil carbon cycling. However, the formation of the SOM profile is currently poorly understood due to equifinality, caused by the entanglement of several processes: input from roots, mixing due to bioturbation, and organic matter leaching. In this study we quantified the contribution of these three processes using Bayesian parameter estimation for the mechanistic SOM profile model SOMPROF. Based on organic carbon measurements, 13 parameters related to decomposition and transport of organic matter were estimated for two temperature forest soils: an Arenosol with a mor humus form (Loobos, The Netherlands), and a Cambisol with mull type humus (Hainich, Germany). Furthermore, the use of the radioisotope 210Pbex as tracer for vertical SOM transport was studied. For Loobos the calibration results demonstrate the importance of liquid phase transport for shaping the vertical SOM profile, while the effects of bioturbation are generally negligible. These results are in good agreement with expectations given in situ conditions. For Hainich the calibration offered three distinct explanations for the observations (three modes in the posterior distribution). With the addition of 210Pbex data and prior knowledge, as well as additional information about in situ conditions, we were able to identify the most likely explanation, which identified root litter input as the dominant process for the SOM profile. For both sites the organic matter appears to comprise mainly adsorbed but potentially leachable material, pointing to the importance of organo-mineral interactions. Furthermore, organic matter in the mineral soil appears to be mainly derived from root litter, supporting previous studies that highlighted the importance of root input for soil carbon sequestration. The 210Pbex measurements added only slight additional constraint on the estimated parameters. However, with sufficient replicate measurements and possibly in combination with other tracers, this isotope may still hold value as tracer for a SOM transport. [ABSTRACT FROM AUTHOR]

Terrestrial productivity in semi-arid woodlands is strongly susceptible to changes in precipitation, and semi-arid woodlands constitute an important element of the global water and carbon cycles. Here, we use the Carbon Cycle Data Assimilation System (CC-DAS) to investigate the mechanisms controlling ecological and hydrogical activities for a semi-arid savanna woodland site in Maun, Botswana. Twenty-four eco-hydrological process parameters of a terrestrial ecosystem model are optimized against two data streams either separately or simultaneously: daily averaged latent heat flux (LHF) derived from eddy covariance measurement, and decadal fraction of absorbed photosynthetically active radiation (FAPAR) derived from Sea-viewingWide Field-of-view Sensor (SeaWiFS). Assimilation of both LHF and FAPAR for the years 2000 and 2001 leads to improved agreement between measured and simulated quantities not only for LHF and FAPAR, but also for photosynthetic CO2 uptake. The closest agreement is found for each observed data stream when only the same data stream is assimilated. The mean uncertainty reduction (relative to the prior) over all parameters is 16.1% for the simultaneous assimilation of LHF and FAPAR, 9.2% for assimilating LHF only, and 7.8% for assimilating FAPAR only. Furthermore, the set of parameters with the highest uncertainty reduction is similar between assimilating only FAPAR or only LHF. The highest uncertainty reduction is found for a parameter describing maximum plant-available soil moisture for all three cases. This indicates that not only LHF but also satellite-derived FAPAR data can be used to constrain and indirectly observe hydrological quantities. [ABSTRACT FROM AUTHOR]