The FLUXCOM ensemble of global land-atmosphere energy fluxes.

Although a key driver of Earth's climate system, global land-atmosphere energy fluxes are poorly constrained. Here we use machine learning to merge energy flux measurements from FLUXNET eddy covariance towers with remote sensing and meteorological data to estimate global gridded net radiation, latent and sensible heat and their uncertainties. The resulting FLUXCOM database comprises 147 products in two setups: (1) 0.0833° resolution using MODIS remote sensing data (RS) and (2) 0.5° resolution using remote sensing and meteorological data (RS + METEO). Within each setup we use a full factorial design across machine learning methods, forcing datasets and energy balance closure corrections. For RS and RS + METEO setups respectively, we estimate 2001-2013 global (±1 s.d.) net radiation as 75.49 ± 1.39 W m-2 and 77.52 ± 2.43 W m-2, sensible heat as 32.39 ± 4.17 W m-2 and 35.58 ± 4.75 W m-2, and latent heat flux as 39.14 ± 6.60 W m-2 and 39.49 ± 4.51 W m-2 (as evapotranspiration, 75.6 ± 9.8 × 103 km3 yr-1 and 76 ± 6.8 × 103 km3 yr-1). FLUXCOM products are suitable to quantify global land-atmosphere interactions and benchmark land surface model simulations.

Memory effects of climate and vegetation affecting net ecosystem CO2 fluxes in global forests.

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.

Compensatory water effects link yearly global land CO2 sink changes to temperature.

Large interannual variations in the measured growth rate of atmospheric carbon dioxide (CO2) originate primarily from fluctuations in carbon uptake by land ecosystems. It remains uncertain, however, to what extent temperature and water availability control the carbon balance of land ecosystems across spatial and temporal scales. Here we use empirical models based on eddy covariance data and process-based models to investigate the effect of changes in temperature and water availability on gross primary productivity (GPP), terrestrial ecosystem respiration (TER) and net ecosystem exchange (NEE) at local and global scales. We find that water availability is the dominant driver of the local interannual variability in GPP and TER. To a lesser extent this is true also for NEE at the local scale, but when integrated globally, temporal NEE variability is mostly driven by temperature fluctuations. We suggest that this apparent paradox can be explained by two compensatory water effects. Temporal water-driven GPP and TER variations compensate locally, dampening water-driven NEE variability. Spatial water availability anomalies also compensate, leaving a dominant temperature signal in the year-to-year fluctuations of the land carbon sink. These findings help to reconcile seemingly contradictory reports regarding the importance of temperature and water in controlling the interannual variability of the terrestrial carbon balance. Our study indicates that spatial climate covariation drives the global carbon cycle response.

Diagnosing the Dynamics of Observed and Simulated Ecosystem Gross Primary Productivity with Time Causal Information Theory Quantifiers.

Data analysis and model-data comparisons in the environmental sciences require diagnostic measures that quantify time series dynamics and structure, and are robust to noise in observational data. This paper investigates the temporal dynamics of environmental time series using measures quantifying their information content and complexity. The measures are used to classify natural processes on one hand, and to compare models with observations on the other. The present analysis focuses on the global carbon cycle as an area of research in which model-data integration and comparisons are key to improving our understanding of natural phenomena. We investigate the dynamics of observed and simulated time series of Gross Primary Productivity (GPP), a key variable in terrestrial ecosystems that quantifies ecosystem carbon uptake. However, the dynamics, patterns and magnitudes of GPP time series, both observed and simulated, vary substantially on different temporal and spatial scales. We demonstrate here that information content and complexity, or Information Theory Quantifiers (ITQ) for short, serve as robust and efficient data-analytical and model benchmarking tools for evaluating the temporal structure and dynamical properties of simulated or observed time series at various spatial scales. At continental scale, we compare GPP time series simulated with two models and an observations-based product. This analysis reveals qualitative differences between model evaluation based on ITQ compared to traditional model performance metrics, indicating that good model performance in terms of absolute or relative error does not imply that the dynamics of the observations is captured well. Furthermore, we show, using an ensemble of site-scale measurements obtained from the FLUXNET archive in the Mediterranean, that model-data or model-model mismatches as indicated by ITQ can be attributed to and interpreted as differences in the temporal structure of the respective ecological time series. At global scale, our understanding of C fluxes relies on the use of consistently applied land models. Here, we use ITQ to evaluate model structure: The measures are largely insensitive to climatic scenarios, land use and atmospheric gas concentrations used to drive them, but clearly separate the structure of 13 different land models taken from the CMIP5 archive and an observations-based product. In conclusion, diagnostic measures of this kind provide data-analytical tools that distinguish different types of natural processes based solely on their dynamics, and are thus highly suitable for environmental science applications such as model structural diagnostics.

Physical Limits of Solar Energy Conversion in the Earth System.

Solar energy provides by far the greatest potential for energy generation among all forms of renewable energy. Yet, just as for any form of energy conversion, it is subject to physical limits. Here we review the physical limits that determine how much energy can potentially be generated out of sunlight using a combination of thermodynamics and observed climatic variables. We first explain how the first and second law of thermodynamics constrain energy conversions and thereby the generation of renewable energy, and how this applies to the conversions of solar radiation within the Earth system. These limits are applied to the conversion of direct and diffuse solar radiation - which relates to concentrated solar power (CSP) and photovoltaic (PV) technologies as well as biomass production or any other photochemical conversion - as well as solar radiative heating, which generates atmospheric motion and thus relates to wind power technologies. When these conversion limits are applied to observed data sets of solar radiation at the land surface, it is estimated that direct concentrated solar power has a potential on land of up to 11.6 PW (1 PW=10(15) W), whereas photovoltaic power has a potential of up to 16.3 PW. Both biomass and wind power operate at much lower efficiencies, so their potentials of about 0.3 and 0.1 PW are much lower. These estimates are considerably lower than the incoming flux of solar radiation of 175 PW. When compared to a 2012 primary energy demand of 17 TW, the most direct uses of solar radiation, e.g., by CSP or PV, have thus by far the greatest potential to yield renewable energy requiring the least space to satisfy the human energy demand. Further conversions into solar-based fuels would be reduced by further losses which would lower these potentials. The substantially greater potential of solar-based renewable energy compared to other forms of renewable energy simply reflects much fewer and lower unavoidable conversion losses when solar radiation is directly converted into renewable energy.

Two methods for estimating limits to large-scale wind power generation.

Wind turbines remove kinetic energy from the atmospheric flow, which reduces wind speeds and limits generation rates of large wind farms. These interactions can be approximated using a vertical kinetic energy (VKE) flux method, which predicts that the maximum power generation potential is 26% of the instantaneous downward transport of kinetic energy using the preturbine climatology. We compare the energy flux method to the Weather Research and Forecasting (WRF) regional atmospheric model equipped with a wind turbine parameterization over a 10(5) km2 region in the central United States. The WRF simulations yield a maximum generation of 1.1 We⋅m(-2), whereas the VKE method predicts the time series while underestimating the maximum generation rate by about 50%. Because VKE derives the generation limit from the preturbine climatology, potential changes in the vertical kinetic energy flux from the free atmosphere are not considered. Such changes are important at night when WRF estimates are about twice the VKE value because wind turbines interact with the decoupled nocturnal low-level jet in this region. Daytime estimates agree better to 20% because the wind turbines induce comparatively small changes to the downward kinetic energy flux. This combination of downward transport limits and wind speed reductions explains why large-scale wind power generation in windy regions is limited to about 1 We⋅m(-2), with VKE capturing this combination in a comparatively simple way.

Cross-modulated amplitudes and frequencies characterize interacting components in complex systems.

The dynamics of complex systems is characterized by oscillatory components on many time scales. To study the interactions between these components we analyze the cross modulation of their instantaneous amplitudes and frequencies, separating synchronous and antisynchronous modulation. We apply our novel technique to brain-wave oscillations in the human electroencephalogram and show that interactions between the alpha wave and the delta or beta wave oscillators as well as spatial interactions can be quantified and related with physiological conditions (e.g., sleep stages). Our approach overcomes the limitation to oscillations with similar frequencies and enables us to quantify directly nonlinear effects such as positive or negative frequency modulation.