The FLUXNET2015 dataset provides ecosystem-scale data on CO₂, water, and energy exchange between the biosphere and the atmosphere, and other meteorological and biological measurements, from 212 sites around the globe (over 1500 site-years, up to and including year 2014). These sites, independently managed and operated, voluntarily contributed their data to create global datasets. Data were quality controlled and processed using uniform methods, to improve consistency and intercomparability across sites. The dataset is already being used in a number of applications, including ecophysiology studies, remote sensing studies, and development of ecosystem and Earth system models. FLUXNET2015 includes derived-data products, such as gap-filled time series, ecosystem respiration and photosynthetic uptake estimates, estimation of uncertainties, and metadata about the measurements, presented for the first time in this paper. In addition, 206 of these sites are for the first time distributed under a Creative Commons (CC-BY 4.0) license. This paper details this enhanced dataset and the processing methods, now made available as open-source codes, making the dataset more accessible, transparent, and reproducible.

Managing environmental resources under conditions of climate change and extreme climate events remains among the most challenging research tasks in the field of sustainable development. A particular challenge in many regions such as East Africa is often the lack of sufficiently long-term and spatially representative observed climate data. To overcome this data challenge we used a combination of accessible data sources based on station data, earth observations by remote sensing, and regional climate models. The accuracy of the Africa Rainfall Climatology version 2.0 (ARC2), Climate Hazards Group InfraRed Precipitation (CHIRP), CHIRP with Station data (CHIRPS), Observational-Reanalysis Hybrid (ORH), and regional climate models (RCMs) are evaluated against station data obtained from the respective national weather services and international databases. We did so by performing a comparison in three ways: point to pixel, point to area grid cell average, and stations' average to area grid cell average over 21 regions of East Africa: 17 in Ethiopia, 2 in Kenya, and 2 in Tanzania. We found that the latter method provides better correlation and significantly reduces biases and errors. The correlations were analysed at daily, dekadal (10 days), and monthly resolution for rainfall and maximum and minimum temperature (T_max and T_min) covering the period of 1983–2005. At a daily timescale, CHIRPS, followed by ARC2 and CHIRP, is the best performing rainfall product compared to ORH, individual RCMs (I-RCM), and RCMs' mean (RCMs). CHIRPS captures the daily rainfall characteristics well, such as average daily rainfall, amount of wet periods, and total rainfall. Compared to CHIRPS, ARC2 showed higher underestimation of the total (−30 %) and daily (−14 %) rainfall. CHIRP, on the other hand, showed higher underestimation of the average daily rainfall (−53 %) and duration of dry periods (−29 %). Overall, the evaluation revealed that in terms of multiple statistical measures used on daily, dekadal, and monthly timescales, CHIRPS, CHIRP, and ARC2 are the best performing rainfall products, while ORH, I-RCM, and RCMs are the worst performing products. For T_max and T_min, ORH was identified as the most suitable product compared to I-RCM and RCMs. Our results indicate that CHIRPS (rainfall) and ORH (T_max and T_min), with higher spatial resolution, should be the preferential data sources to be used for climate change and hydrological studies in areas of East Africa where station data are not accessible.

A Correction to this paper has been published: https://doi.org/10.1038/s41597-021-00851-9.

Peatlands store substantial amounts of carbon and are vulnerable to climate change. We present a modified version of the Organising Carbon and Hydrology In Dynamic Ecosystems (ORCHIDEE) land surface model for simulating the hydrology, surface energy, and CO₂ fluxes of peatlands on daily to annual timescales. The model includes a separate soil tile in each 0.5° grid cell, defined from a global peatland map and identified with peat-specific soil hydraulic properties. Runoff from non-peat vegetation within a grid cell containing a fraction of peat is routed to this peat soil tile, which maintains shallow water tables. The water table position separates oxic from anoxic decomposition. The model was evaluated against eddy-covariance (EC) observations from 30 northern peatland sites, with the maximum rate of carboxylation (V_cmax) being optimized at each site. Regarding short-term day-to-day variations, the model performance was good for gross primary production (GPP) (r² = 0.76; Nash–Sutcliffe modeling efficiency, MEF = 0.76) and ecosystem respiration (ER, r² = 0.78, MEF = 0.75), with lesser accuracy for latent heat fluxes (LE, r² = 0.42, MEF = 0.14) and and net ecosystem CO₂ exchange (NEE, r² = 0.38, MEF = 0.26). Seasonal variations in GPP, ER, NEE, and energy fluxes on monthly scales showed moderate to high r² values (0.57–0.86). For spatial across-site gradients of annual mean GPP, ER, NEE, and LE, r² values of 0.93, 0.89, 0.27, and 0.71 were achieved, respectively. Water table (WT) variation was not well predicted (r² < 0.1), likely due to the uncertain water input to the peat from surrounding areas. However, the poor performance of WT simulation did not greatly affect predictions of ER and NEE. We found a significant relationship between optimized V_cmax and latitude (temperature), which better reflects the spatial gradients of annual NEE than using an average V_cmax value.

The response of crops to changing climate and atmospheric CO₂ concentration ([CO₂]) could have large effects on food production, and impact carbon, water, and energy fluxes, causing feedbacks to the climate. To simulate the response of temperate crops to changing climate and [CO₂], which accounts for the specific phenology of crops mediated by management practice, we describe here the development of a process-oriented terrestrial biogeochemical model named ORCHIDEE-CROP (v0), which integrates a generic crop phenology and harvest module, and a very simple parameterization of nitrogen fertilization, into the land surface model (LSM) ORCHIDEEv196, in order to simulate biophysical and biochemical interactions in croplands, as well as plant productivity and harvested yield. The model is applicable for a range of temperate crops, but is tested here using maize and winter wheat, with the phenological parameterizations of two European varieties originating from the STICS agronomical model. We evaluate the ORCHIDEE-CROP (v0) model against eddy covariance and biometric measurements at seven winter wheat and maize sites in Europe. The specific ecosystem variables used in the evaluation are CO₂ fluxes (net ecosystem exchange, NEE), latent heat, and sensible heat fluxes. Additional measurements of leaf area index (LAI) and aboveground biomass and yield are used as well. Evaluation results revealed that ORCHIDEE-CROP (v0) reproduced the observed timing of crop development stages and the amplitude of the LAI changes. This is in contrast to ORCHIDEEv196 where, by default, crops have the same phenology as grass. A halving of the root mean square error for LAI from 2.38 ± 0.77 to 1.08 ± 0.34 m² m⁻² was obtained when ORCHIDEEv196 and ORCHIDEE-CROP (v0) were compared across the seven study sites. Improved crop phenology and carbon allocation led to a good match between modeled and observed aboveground biomass (with a normalized root mean squared error (NRMSE) of 11.0–54.2 %), crop yield, daily carbon and energy fluxes (with a NRMSE of ∼ 9.0–20.1 and ∼ 9.4–22.3 % for NEE), and sensible and latent heat fluxes. The simulated yields for winter wheat and maize from ORCHIDEE-CROP (v0) showed a good match with the simulated results from STICS for three sites with available crop yield observations, where the average NRMSE was ∼ 8.8 %. The model data misfit for energy fluxes were within the uncertainties of the measurements, which themselves showed an incomplete energy balance closure within the range 80.6–86.3 %. The remaining discrepancies between the modeled and observed LAI and other variables at specific sites were partly attributable to unrealistic representations of management events by the model. ORCHIDEE-CROP (v0) has the ability to capture the spatial gradients of carbon and energy-related variables, such as gross primary productivity, NEE, and sensible and latent heat fluxes across the sites in Europe, which is an important requirement for future spatially explicit simulations. Further improvement of the model, with an explicit parameterization of nutritional dynamics and management, is expected to improve its predictive ability to simulate croplands in an Earth system model.