Browsing by Author "Reichstein, Markus"

The energy balance at most surface-atmosphere flux research sites remains unclosed. The mechanisms underlying the discrepancy between measured energy inputs and outputs across the global FLUXNET tower network are still under debate. Recent reviews have identified exchange processes and turbulent motions at large spatial and temporal scales in heterogeneous landscapes as the primary cause of the lack of energy balance closure at some intensively-researched sites, while unmeasured storage terms cannot be ruled out as a dominant contributor to the lack of energy balance closure at many other sites. We analyzed energy balance closure across 173 ecosystems in the FLUXNET database and explored the relationship between energy balance closure and landscape heterogeneity using MODIS products and GLOBEstat elevation data. Energy balance closure per research site (C-EBS)averaged 0.84 +/- 0.20, with best average closures in evergreen broadleaf forests and savannas (0.91-0.94) and worst average closures in crops, deciduous broadleaf forests, mixed forests and wetlands (0.70-0.78). Half-hourly or hourly energy balance closure on a percent basis increased with friction velocity (u.) and was highest on average under near-neutral atmospheric conditions. C-EBS was significantly related to mean precipitation, gross primary productivity and landscape-level enhanced vegetation index (EVI) from MODIS, and the variability in elevation, MODIS plant functional type, and MODIS EVI. A linear model including landscape-level variability in both EVI and elevation, mean precipitation, and an interaction term between EVI variability and precipitation had the lowest Akaike's information criterion value. C-EBS in landscapes with uniform plant functional type approached 0.9 and C-EBS in landscapes with uniform EVI approached 1. These results suggest that landscape-level heterogeneity in vegetation and topography cannot be ignored as a contributor to incomplete energy balance closure at the flux network level, although net radiation measurements, biological energy assimilation, unmeasured storage terms, and the importance of good practice including site selection when making flux measurements should not be discounted. Our results suggest that future research should focus on the quantitative mechanistic relationships between energy balance closure and landscape-scale heterogeneity, and the consequences of mesoscale circulations for surface-atmosphere exchange measurements. (C) 2012 Elsevier B.V. All rights reserved.

We apply and compare three widely applicable methods for estimating ecosystem transpiration (T) from eddy covariance (EC) data across 251 FLUXNET sites globally. All three methods are based on the coupled water and carbon relationship, but they differ in assumptions and parameterizations. Intercomparison of the three daily T estimates shows high correlation among methods (R between .89 and .94), but a spread in magnitudes of T/ET (evapotranspiration) from 45% to 77%. When compared at six sites with concurrent EC and sap flow measurements, all three EC-based T estimates show higher correlation to sap flow-based T than EC-based ET. The partitioning methods show expected tendencies of T/ET increasing with dryness (vapor pressure deficit and days since rain) and with leaf area index (LAI). Analysis of 140 sites with high-quality estimates for at least two continuous years shows that T/ET variability was 1.6 times higher across sites than across years. Spatial variability of T/ET was primarily driven by vegetation and soil characteristics (e.g., crop or grass designation, minimum annual LAI, soil coarse fragment volume) rather than climatic variables such as mean/standard deviation of temperature or precipitation. Overall, T and T/ET patterns are plausible and qualitatively consistent among the different water flux partitioning methods implying a significant advance made for estimating and understanding T globally, while the magnitudes remain uncertain. Our results represent the first extensive EC data-based estimates of ecosystem T permitting a data-driven perspective on the role of plants’ water use for global water and carbon cycling in a changing climate.

Extreme climatic events, such as droughts and heat stress, induce anomalies in ecosystem–atmosphere CO2 fluxes, such as gross primary production (GPP) and ecosystem respiration (Reco), and, hence, can change the net ecosystem carbon balance. However, despite our increasing understanding of the underlying mechanisms, the magnitudes of the impacts of different types of extremes on GPP and Reco within and between ecosystems remain poorly predicted. Here we aim to identify the major factors controlling the amplitude of extreme-event impacts on GPP, Reco, and the resulting net ecosystem production (NEP). We focus on the impacts of heat and drought and their combination. We identified hydrometeorological extreme events in consistently downscaled water availability and temperature measurements over a 30-year time period. We then used FLUXNET eddy covariance flux measurements to estimate the CO2 flux anomalies during these extreme events across dominant vegetation types and climate zones. Overall, our results indicate that short-term heat extremes increased respiration more strongly than they downregulated GPP, resulting in a moderate reduction in the ecosystem's carbon sink potential. In the absence of heat stress, droughts tended to have smaller and similarly dampening effects on both GPP and Reco and, hence, often resulted in neutral NEP responses. The combination of drought and heat typically led to a strong decrease in GPP, whereas heat and drought impacts on respiration partially offset each other. Taken together, compound heat and drought events led to the strongest C sink reduction compared to any single-factor extreme. A key insight of this paper, however, is that duration matters most: for heat stress during droughts, the magnitude of impacts systematically increased with duration, whereas under heat stress without drought, the response of Reco over time turned from an initial increase to a downregulation after about 2 weeks. This confirms earlier theories that not only the magnitude but also the duration of an extreme event determines its impact. Our study corroborates the results of several local site-level case studies but as a novelty generalizes these findings on the global scale. Specifically, we find that the different response functions of the two antipodal land–atmosphere fluxes GPP and Reco can also result in increasing NEP during certain extreme conditions. Apparently counterintuitive findings of this kind bear great potential for scrutinizing the mechanisms implemented in state-of-the-art terrestrial biosphere models and provide a benchmark for future model development and testing.

Across a landscape, aquatic-terrestrial interfaces within and between ecosystems are hotspots of organic matter (OM) mineralization. These interfaces are characterized by sharp spatio-temporal changes in environmental conditions, which affect OM properties and thus control OM mineralization and other transformation processes. Consequently, the extent of OM movement at and across aquatic-terrestrial interfaces is crucial in determining OM turnover and carbon (C) cycling at the landscape scale. Here, we propose expanding current concepts in aquatic and terrestrial ecosystem sciences to comprehensively evaluate OM turnover at the landscape scale. We focus on three main concepts toward explaining OM turnover at the landscape scale: the landscape spatiotemporal context, OM turnover described by priming and ecological stoichiometry, and anthropogenic effects as a disruptor of natural OM transfer magnitudes and pathways. A conceptual framework is introduced that allows for discussing the disparities in spatial and temporal scales of OM transfer, changes in environmental conditions, ecosystem connectivity, and microbial–substrate interactions. The potential relevance of priming effects in both terrestrial and aquatic systems is addressed. For terrestrial systems, we hypothesize that the interplay between the influx of OM, its corresponding elemental composition, and the elemental demand of the microbial communities may alleviate spatial and metabolic thresholds. In comparison, substrate level OM dynamics may be substantially different in aquatic systems due to matrix effects that accentuate the role of abiotic conditions, substrate quality, and microbial community dynamics. We highlight the disproportionate impact anthropogenic activities can have on OM cycling across the landscape. This includes reversing natural OM flows through the landscape, disrupting ecosystem connectivity, and nutrient additions that cascade across the landscape. This knowledge is crucial for a better understanding of OM cycling in a landscape context, in particular since terrestrial and aquatic compartments may respond differently to the ongoing changes in climate, land use, and other anthropogenic interferences.

Fundamental axes of variation in plant traits result from trade-offs between costs and benefits of resource-use strategies at the leaf scale. However, it is unclear whether similar trade-offs propagate to the ecosystem level. Here, we test whether trait correlation patterns predicted by three well-known leaf- and plant-level coordination theories – the leaf economics spectrum, the global spectrum of plant form and function, and the least-cost hypothesis – are also observed between community mean traits and ecosystem processes. We combined ecosystem functional properties from FLUXNET sites, vegetation properties, and community mean plant traits into three corresponding principal component analyses. We find that the leaf economics spectrum (90 sites), the global spectrum of plant form and function (89 sites), and the least-cost hypothesis (82 sites) all propagate at the ecosystem level. However, we also find evidence of additional scale-emergent properties. Evaluating the coordination of ecosystem functional properties may aid the development of more realistic global dynamic vegetation models with critical empirical data, reducing the uncertainty of climate change projections.

Forests play a crucial role in the global carbon (C) cycle by storing and sequestering a substantial amount of C in the terrestrial biosphere. Due to temporal dynamics in climate and vegetation activity, there are significant regional variations in carbon dioxide (CO2) fluxes between the biosphere and atmosphere in forests that are affecting the global C cycle. Current forest CO2 flux dynamics are controlled by instantaneous climate, soil, and vegetation conditions, which carry legacy effects from disturbances and extreme climate events. Our level of understanding from the legacies of these processes on net CO2 fluxes is still limited due to their complexities and their long-term effects. Here, we combined remote sensing, climate, and eddy-covariance flux data to study net ecosystem CO2 exchange (NEE) at 185 forest sites globally. Instead of commonly used non-dynamic statistical methods, we employed a type of recurrent neural network (RNN), called Long Short-Term Memory network (LSTM) that captures information from the vegetation and climate's temporal dynamics. The resulting data-driven model integrates interannual and seasonal variations of climate and vegetation by using Landsat and climate data at each site. The presented LSTM algorithm was able to effectively describe the overall seasonal variability (Nash-Sutcliffe efficiency, NSE = 0.66) and across-site (NSE = 0.42) variations in NEE, while it had less success in predicting specific seasonal and interannual anomalies (NSE = 0.07). This analysis demonstrated that an LSTM approach with embedded climate and vegetation memory effects outperformed a non-dynamic statistical model (i.e. Random Forest) for estimating NEE. Additionally, it is shown that the vegetation mean seasonal cycle embeds most of the information content to realistically explain the spatial and seasonal variations in NEE. These findings show the relevance of capturing memory effects from both climate and vegetation in quantifying spatio-temporal variations in forest NEE.

In the version of this Article originally published, the wrong Supplementary Information pdf was uploaded, in which the figures did not correspond with those mentioned in the main text and the R code was not presented properly. This has now been replaced.

Forests dominate carbon (C) exchanges between the terrestrial biosphere and the atmosphere on land. In the long term, the net carbon flux between forests and the atmosphere has been significantly impacted by changes in forest cover area and structure due to ecological disturbances and management activities. Current empirical approaches for estimating net ecosystem productivity (NEP) rarely consider forest age as a predictor, which represents variation in physiological processes that can respond differently to environmental drivers, and regrowth following disturbance. Here, we conduct an observational synthesis to empirically determine to what extent climate, soil properties, nitrogen deposition, forest age and management influence the spatial and interannual variability of forest NEP across 126 forest eddy-covariance flux sites worldwide. The empirical models explained up to 62% and 71% of spatio-temporal and across-site variability of annual NEP, respectively. An investigation of model structures revealed that forest age was a dominant factor of NEP spatio-temporal variability in both space and time at the global scale as compared to abiotic factors such as nutrient availability, soil characteristics, and climate. These findings emphasize the importance of forest age in quantifying spatio-temporal variation in NEP using empirical approaches.

The total uptake of carbon dioxide by ecosystems via photosynthesis (gross primary productivity, GPP) is the largest flux in the global carbon cycle. A key ecosystem functional property determining GPP is the photosynthetic capacity at light saturation (GPPsat), and its interannual variability (IAV) is propagated to the net land–atmosphere exchange of CO2. Given the importance of understanding the IAV in CO2 fluxes for improving the predictability of the global carbon cycle, we have tested a range of alternative hypotheses to identify potential drivers of the magnitude of IAV in GPPsat in forest ecosystems. Our results show that while the IAV in GPPsat within sites is closely related to air temperature and soil water availability fluctuations, the magnitude of IAV in GPPsat is related to stand age and biodiversity (R2 = 0.55, P < 0.0001). We find that the IAV of GPPsat is greatly reduced in older and more diverse forests, and is higher in younger forests with few dominant species. Older and more diverse forests seem to dampen the effect of climate variability on the carbon cycle irrespective of forest type. Preserving old forests and their diversity would therefore be beneficial in reducing the effect of climate variability on Earth's forest ecosystems.

Characterization of state-dependent model biases in land surface models can highlight model deficiencies, and provide new insights into model development. In this study, artificial neural networks (ANNs) are used to estimate the state-dependent biases of a land surface model (ORCHIDEE: ORganising Carbon and Hydrology in Dynamic EcosystEms). To characterize state-dependent biases in ORCHIDEE, we use multi-year flux measurements made at 125 eddy covariance sites that cover 7 different plant functional types (PFTs) and 5 climate groups. We determine whether the state-dependent model biases in five flux variables (H: sensible heat, LE: latent heat, NEE: net ecosystem exchange, GPP: gross primary productivity and R-eco: ecosystem respiration) are transferable within and between three different timescales (diurnal, seasonal-annual and interannual), and between sites (categorized by PFTs and climate groups). For each flux variable at each site, the spectral decomposition method (singular system analysis) was used to reconstruct time series on the three different timescales. At the site level, we found that the share of state-dependent model biases (hereafter called "error transferability") is larger for seasonal-annual and interannual timescales than for the diurnal timescale, but little error transferability was found between timescales in all flux variables. Thus, performing model evaluations at multiple timescales is essential for diagnostics and future development. For all PFTs, climate groups and timescale components, the state-dependent model biases are found to be transferable between sites within the same PFT and climate group, suggesting that specific model developments and improvements based on specific eddy covariance sites can be used to enhance the model performance at other sites within the same PFT-climate group. This also supports the legitimacy of upscaling from the ecosystem scale of eddy covariance sites to the regional scale based on the similarity of PFT and climate group. However, the transferability of state-dependent model biases between PFTs or climate groups is not always found on the seasonal-annual and interannual timescales, which is contrary to transferability found on the diurnal timescale and the original time series. (C) 2012 Elsevier B.V. All rights reserved.

Abstract 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 .