Max Planck Institute for Biogeochemistry

A large body of information on emissions from the various types of biomass burning has been accumulated over the past decade, to a large extent as a result of International Geosphere‐Biosphere Programme/International Global Atmospheric Chemistry research activities. Yet this information has not been readily accessible to the atmospheric chemistry community because it was scattered over a large number of publications and reported in numerous different units and reference systems. We have critically evaluated the presently available data and integrated these into a consistent format. On the basis of this analysis we present a set of emission factors for a large variety of species emitted from biomass fires. Where data were not available, we have proposed estimates based on appropriate extrapolation techniques. We have derived global estimates of pyrogenic emissions for important species emitted by the various types of biomass burning and compared our estimates with results from inverse modeling studies.

There is growing recognition that classifying terrestrial plant species on the basis of their function (into 'functional types') rather than their higher taxonomic identity, is a promising way forward for tackling important ecological questions at the scale of ecosystems, landscapes or biomes. These questions include those on vegetation responses to and vegetation effects on, environmental changes (e.g. changes in climate, atmospheric chemistry, land use or other disturbances). There is also growing consensus about a shortlist of plant traits that should underlie such functional plant classifications, because they have strong predictive power of important ecosystem responses to environmental change and/or they themselves have strong impacts on ecosystem processes. The most favoured traits are those that are also relatively easy and inexpensive to measure for large numbers of plant species. Large international research efforts, promoted by the IGBP–GCTE Programme, are underway to screen predominant plant species in various ecosystems and biomes worldwide for such traits. This paper provides an international methodological protocol aimed at standardising this research effort, based on consensus among a broad group of scientists in this field. It features a practical handbook with step-by-step recipes, with relatively brief information about the ecological context, for 28 functional traits recognised as critical for tackling large-scale ecological questions.

Abstract This paper discusses the advantages and disadvantages of the different methods that separate net ecosystem exchange (NEE) into its major components, gross ecosystem carbon uptake (GEP) and ecosystem respiration ( R eco ). In particular, we analyse the effect of the extrapolation of night‐time values of ecosystem respiration into the daytime; this is usually done with a temperature response function that is derived from long‐term data sets. For this analysis, we used 16 one‐year‐long data sets of carbon dioxide exchange measurements from European and US‐American eddy covariance networks. These sites span from the boreal to Mediterranean climates, and include deciduous and evergreen forest, scrubland and crop ecosystems. We show that the temperature sensitivity of R eco , derived from long‐term (annual) data sets, does not reflect the short‐term temperature sensitivity that is effective when extrapolating from night‐ to daytime. Specifically, in summer active ecosystems the long‐term temperature sensitivity exceeds the short‐term sensitivity. Thus, in those ecosystems, the application of a long‐term temperature sensitivity to the extrapolation of respiration from night to day leads to a systematic overestimation of ecosystem respiration from half‐hourly to annual time‐scales, which can reach >25% for an annual budget and which consequently affects estimates of GEP. Conversely, in summer passive (Mediterranean) ecosystems, the long‐term temperature sensitivity is lower than the short‐term temperature sensitivity resulting in underestimation of annual sums of respiration. We introduce a new generic algorithm that derives a short‐term temperature sensitivity of R eco from eddy covariance data that applies this to the extrapolation from night‐ to daytime, and that further performs a filling of data gaps that exploits both, the covariance between fluxes and meteorological drivers and the temporal structure of the fluxes. While this algorithm should give less biased estimates of GEP and R eco , we discuss the remaining biases and recommend that eddy covariance measurements are still backed by ancillary flux measurements that can reduce the uncertainties inherent in the eddy covariance data.

Abstract The Lund–Potsdam–Jena Dynamic Global Vegetation Model (LPJ) combines process‐based, large‐scale representations of terrestrial vegetation dynamics and land‐atmosphere carbon and water exchanges in a modular framework. Features include feedback through canopy conductance between photosynthesis and transpiration and interactive coupling between these ‘fast’ processes and other ecosystem processes including resource competition, tissue turnover, population dynamics, soil organic matter and litter dynamics and fire disturbance. Ten plants functional types (PFTs) are differentiated by physiological, morphological, phenological, bioclimatic and fire‐response attributes. Resource competition and differential responses to fire between PFTs influence their relative fractional cover from year to year. Photosynthesis, evapotranspiration and soil water dynamics are modelled on a daily time step, while vegetation structure and PFT population densities are updated annually. Simulations have been made over the industrial period both for specific sites where field measurements were available for model evaluation, and globally on a 0.5°° × 0.5°° grid. Modelled vegetation patterns are consistent with observations, including remotely sensed vegetation structure and phenology. Seasonal cycles of net ecosystem exchange and soil moisture compare well with local measurements. Global carbon exchange fields used as input to an atmospheric tracer transport model (TM2) provided a good fit to observed seasonal cycles of CO 2 concentration at all latitudes. Simulated inter‐annual variability of the global terrestrial carbon balance is in phase with and comparable in amplitude to observed variability in the growth rate of atmospheric CO 2 . Global terrestrial carbon and water cycle parameters (pool sizes and fluxes) lie within their accepted ranges. The model is being used to study past, present and future terrestrial ecosystem dynamics, biochemical and biophysical interactions between ecosystems and the atmosphere, and as a component of coupled Earth system models.

Abstract Eleven coupled climate–carbon cycle models used a common protocol to study the coupling between climate change and the carbon cycle. The models were forced by historical emissions and the Intergovernmental Panel on Climate Change (IPCC) Special Report on Emissions Scenarios (SRES) A2 anthropogenic emissions of CO2 for the 1850–2100 time period. For each model, two simulations were performed in order to isolate the impact of climate change on the land and ocean carbon cycle, and therefore the climate feedback on the atmospheric CO2 concentration growth rate. There was unanimous agreement among the models that future climate change will reduce the efficiency of the earth system to absorb the anthropogenic carbon perturbation. A larger fraction of anthropogenic CO2 will stay airborne if climate change is accounted for. By the end of the twenty-first century, this additional CO2 varied between 20 and 200 ppm for the two extreme models, the majority of the models lying between 50 and 100 ppm. The higher CO2 levels led to an additional climate warming ranging between 0.1° and 1.5°C. All models simulated a negative sensitivity for both the land and the ocean carbon cycle to future climate. However, there was still a large uncertainty on the magnitude of these sensitivities. Eight models attributed most of the changes to the land, while three attributed it to the ocean. Also, a majority of the models located the reduction of land carbon uptake in the Tropics. However, the attribution of the land sensitivity to changes in net primary productivity versus changes in respiration is still subject to debate; no consensus emerged among the models.

Carbon Cycle and Climate Change As climate change accelerates, it is important to know the likely impact of climate change on the carbon cycle (see the Perspective by Reich ). Gross primary production (GPP) is a measure of the amount of CO 2 removed from the atmosphere every year to fuel photosynthesis. Beer et al. (p. 834 , published online 5 July) used a combination of observation and calculation to estimate that the total GPP by terrestrial plants is around 122 billion tons per year; in comparison, burning fossil fuels emits about 7 billion tons annually. Thirty-two percent of this uptake occurs in tropical forests, and precipitation controls carbon uptake in more than 40% of vegetated land. The temperature sensitivity (Q10) of ecosystem respiratory processes is a key determinant of the interaction between climate and the carbon cycle. Mahecha et al. (p. 838 , published online 5 July) now show that the Q10 of ecosystem respiration is invariant with respect to mean annual temperature, independent of the analyzed ecosystem type, with a global mean value for Q10 of 1.6. This level of temperature sensitivity suggests a less-pronounced climate sensitivity of the carbon cycle than assumed by recent climate models.

A high-resolution deuterium profile is now available along the entire European Project for Ice Coring in Antarctica Dome C ice core, extending this climate record back to marine isotope stage 20.2, approximately 800,000 years ago. Experiments performed with an atmospheric general circulation model including water isotopes support its temperature interpretation. We assessed the general correspondence between Dansgaard-Oeschger events and their smoothed Antarctic counterparts for this Dome C record, which reveals the presence of such features with similar amplitudes during previous glacial periods. We suggest that the interplay between obliquity and precession accounts for the variable intensity of interglacial periods in ice core records.

Abstract Plant traits – the morphological, anatomical, physiological, biochemical and phenological characteristics of plants and their organs – determine how primary producers respond to environmental factors, affect other trophic levels, influence ecosystem processes and services and provide a link from species richness to ecosystem functional diversity. Trait data thus represent the raw material for a wide range of research from evolutionary biology, community and functional ecology to biogeography. Here we present the global database initiative named TRY, which has united a wide range of the plant trait research community worldwide and gained an unprecedented buy‐in of trait data: so far 93 trait databases have been contributed. The data repository currently contains almost three million trait entries for 69 000 out of the world's 300 000 plant species, with a focus on 52 groups of traits characterizing the vegetative and regeneration stages of the plant life cycle, including growth, dispersal, establishment and persistence. A first data analysis shows that most plant traits are approximately log‐normally distributed, with widely differing ranges of variation across traits. Most trait variation is between species (interspecific), but significant intraspecific variation is also documented, up to 40% of the overall variation. Plant functional types (PFTs), as commonly used in vegetation models, capture a substantial fraction of the observed variation – but for several traits most variation occurs within PFTs, up to 75% of the overall variation. In the context of vegetation models these traits would better be represented by state variables rather than fixed parameter values. The improved availability of plant trait data in the unified global database is expected to support a paradigm shift from species to trait‐based ecology, offer new opportunities for synthetic plant trait research and enable a more realistic and empirically grounded representation of terrestrial vegetation in Earth system models.

Aerosols are of central importance for atmospheric chemistry and physics, the biosphere, climate, and public health. The airborne solid and liquid particles in the nanometer to micrometer size range influence the energy balance of the Earth, the hydrological cycle, atmospheric circulation, and the abundance of greenhouse and reactive trace gases. Moreover, they play important roles in the reproduction of biological organisms and can cause or enhance diseases. The primary parameters that determine the environmental and health effects of aerosol particles are their concentration, size, structure, and chemical composition. These parameters, however, are spatially and temporally highly variable. The quantification and identification of biological particles and carbonaceous components of fine particulate matter in the air (organic compounds and black or elemental carbon, respectively) represent demanding analytical challenges. This Review outlines the current state of knowledge, major open questions, and research perspectives on the properties and interactions of atmospheric aerosols and their effects on climate and human health.

Abstract. Accurate assessment of anthropogenic carbon dioxide (CO2) emissions andtheir redistribution among the atmosphere, ocean, and terrestrial biospherein a changing climate – the “global carbon budget” – is important tobetter understand the global carbon cycle, support the development ofclimate policies, and project future climate change. Here we describe andsynthesize data sets and methodology to quantify the five major componentsof the global carbon budget and their uncertainties. Fossil CO2emissions (EFOS) are based on energy statistics and cement productiondata, while emissions from land-use change (ELUC), mainlydeforestation, are based on land use and land-use change data andbookkeeping models. Atmospheric CO2 concentration is measured directlyand its growth rate (GATM) is computed from the annual changes inconcentration. The ocean CO2 sink (SOCEAN) and terrestrialCO2 sink (SLAND) are estimated with global process modelsconstrained by observations. The resulting carbon budget imbalance(BIM), the difference between the estimated total emissions and theestimated changes in the atmosphere, ocean, and terrestrial biosphere, is ameasure of imperfect data and understanding of the contemporary carboncycle. All uncertainties are reported as ±1σ. For the lastdecade available (2010–2019), EFOS was 9.6 ± 0.5 GtC yr−1 excluding the cement carbonation sink (9.4 ± 0.5 GtC yr−1 when the cement carbonation sink is included), andELUC was 1.6 ± 0.7 GtC yr−1. For the same decade, GATM was 5.1 ± 0.02 GtC yr−1 (2.4 ± 0.01 ppm yr−1), SOCEAN 2.5 ± 0.6 GtC yr−1, and SLAND 3.4 ± 0.9 GtC yr−1, with a budgetimbalance BIM of −0.1 GtC yr−1 indicating a near balance betweenestimated sources and sinks over the last decade. For the year 2019 alone, thegrowth in EFOS was only about 0.1 % with fossil emissions increasingto 9.9 ± 0.5 GtC yr−1 excluding the cement carbonation sink (9.7 ± 0.5 GtC yr−1 when cement carbonation sink is included), and ELUC was 1.8 ± 0.7 GtC yr−1, for total anthropogenic CO2 emissions of 11.5 ± 0.9 GtC yr−1 (42.2 ± 3.3 GtCO2). Also for 2019, GATM was5.4 ± 0.2 GtC yr−1 (2.5 ± 0.1 ppm yr−1), SOCEANwas 2.6 ± 0.6 GtC yr−1, and SLAND was 3.1 ± 1.2 GtC yr−1, with a BIM of 0.3 GtC. The global atmospheric CO2concentration reached 409.85 ± 0.1 ppm averaged over 2019. Preliminarydata for 2020, accounting for the COVID-19-induced changes in emissions,suggest a decrease in EFOS relative to 2019 of about −7 % (medianestimate) based on individual estimates from four studies of −6 %, −7 %,−7 % (−3 % to −11 %), and −13 %. Overall, the mean and trend in thecomponents of the global carbon budget are consistently estimated over theperiod 1959–2019, but discrepancies of up to 1 GtC yr−1 persist for therepresentation of semi-decadal variability in CO2 fluxes. Comparison ofestimates from diverse approaches and observations shows (1) no consensusin the mean and trend in land-use change emissions over the last decade, (2)a persistent low agreement between the different methods on the magnitude ofthe land CO2 flux in the northern extra-tropics, and (3) an apparentdiscrepancy between the different methods for the ocean sink outside thetropics, particularly in the Southern Ocean. This living data updatedocuments changes in the methods and data sets used in this new globalcarbon budget and the progress in understanding of the global carbon cyclecompared with previous publications of this data set (Friedlingstein et al.,2019; Le Quéré et al., 2018b, a, 2016, 2015b, a, 2014,2013). The data presented in this work are available at https://doi.org/10.18160/gcp-2020 (Friedlingstein et al., 2020).

Chapter 3 of Managing the Risks of Extreme Events and Disasters to Advance Climate Change Adaptation. A Special Report of Working Groups I and II of the Intergovernmental Panel on Climate Change (IPCC). This chapter addresses changes in weather and climate events relevant to extreme impacts and disasters.

The global distribution of dust aerosol is simulated with the Georgia Tech/Goddard Global Ozone Chemistry Aerosol Radiation and Transport (GOCART) model. In this model all topographic lows with bare ground surface are assumed to have accumulated sediments which are potential dust sources. The uplifting of dust particles is expressed as a function of surface wind speed and wetness. The GOCART model is driven by the assimilated meteorological fields from the Goddard Earth Observing System Data Assimilation System (GEOS DAS) which facilitates direct comparison with observations. The model includes seven size classes of mineral dust ranging from 0.1–6 μm radius. The total annual emission is estimated to be between 1604 and 1960 Tg yr −1 in a 5‐year simulation. The model has been evaluated by comparing simulation results with ground‐based measurements and satellite data. The evaluation has been performed by comparing surface concentrations, vertical distributions, deposition rates, optical thickness, and size distributions. The comparisons show that the model results generally agree with the observations without the necessity of invoking any contribution from anthropogenic disturbances to soils. However, the model overpredicts the transport of dust from the Asian sources to the North Pacific. This discrepancy is attributed to an overestimate of small particle emission from the Asian sources.

Abstract. Although the definition and measurement techniques for atmospheric "black carbon" ("BC") or "elemental carbon'' ("EC") have long been subjects of scientific controversy, the recent discovery of light-absorbing carbon that is not black ("brown carbon, Cbrown") makes it imperative to reassess and redefine the components that make up light-absorbing carbonaceous matter (LAC) in the atmosphere. Evidence for the atmospheric presence of Cbrown comes from (1) spectral aerosol light absorption measurements near specific combustion sources, (2) observations of spectral properties of water extracts of continental aerosol, (3) laboratory studies indicating the formation of light-absorbing organic matter in the atmosphere, and (4) indirectly from the chemical analogy of aerosol species to colored natural humic substances. We show that brown carbon may severely bias measurements of "BC" and "EC" over vast parts of the troposphere, especially those strongly polluted by biomass burning, where the mass concentration of Cbrown is high relative to that of soot carbon. Chemical measurements to determine "EC" are biased by the refractory nature of Cbrown as well as by complex matrix interferences. Optical measurements of "BC" suffer from a number of problems: (1) many of the presently used instruments introduce a substantial bias into the determination of aerosol light absorption, (2) there is no unique conversion factor between light absorption and "EC" or "BC" concentration in ambient aerosols, and (3) the difference in spectral properties between the different types of LAC, as well as the chemical complexity of Cbrown, lead to several conceptual as well as practical complications. We also suggest that due to the sharply increasing absorption of Cbrown towards the UV, single-wavelength light absorption measurements may not be adequate for the assessment of absorption of solar radiation in the troposphere. We discuss the possible consequences of these effects for our understanding of tropospheric processes, including their influence on UV-irradiance, atmospheric photochemistry and radiative transfer in clouds.

At eight European field sites, the impact of loss of plant diversity on primary productivity was simulated by synthesizing grassland communities with different numbers of plant species. Results differed in detail at each location, but there was an overall log-linear reduction of average aboveground biomass with loss of species. For a given number of species, communities with fewer functional groups were less productive. These diversity effects occurred along with differences associated with species composition and geographic location. Niche complementarity and positive species interactions appear to play a role in generating diversity-productivity relationships within sites in addition to sampling from the species pool.

Plant traits-the morphological, anatomical, physiological, biochemical and phenological characteristics of plants-determine how plants respond to environmental factors, affect other trophic levels, and influence ecosystem properties and their benefits and detriments to people. Plant trait data thus represent the basis for a vast area of research spanning from evolutionary biology, community and functional ecology, to biodiversity conservation, ecosystem and landscape management, restoration, biogeography and earth system modelling. Since its foundation in 2007, the TRY database of plant traits has grown continuously. It now provides unprecedented data coverage under an open access data policy and is the main plant trait database used by the research community worldwide. Increasingly, the TRY database also supports new frontiers of trait-based plant research, including the identification of data gaps and the subsequent mobilization or measurement of new data. To support this development, in this article we evaluate the extent of the trait data compiled in TRY and analyse emerging patterns of data coverage and representativeness. Best species coverage is achieved for categorical traits-almost complete coverage for 'plant growth form'. However, most traits relevant for ecology and vegetation modelling are characterized by continuous intraspecific variation and trait-environmental relationships. These traits have to be measured on individual plants in their respective environment. Despite unprecedented data coverage, we observe a humbling lack of completeness and representativeness of these continuous traits in many aspects. We, therefore, conclude that reducing data gaps and biases in the TRY database remains a key challenge and requires a coordinated approach to data mobilization and trait measurements. This can only be achieved in collaboration with other initiatives.

Summary The possible responses of ecosystem processes to rising atmospheric CO 2 concentration and climate change are illustrated using six dynamic global vegetation models that explicitly represent the interactions of ecosystem carbon and water exchanges with vegetation dynamics. The models are driven by the IPCC IS92a scenario of rising CO 2 ( Wigley et al . 1991 ), and by climate changes resulting from effective CO 2 concentrations corresponding to IS92a, simulated by the coupled ocean atmosphere model HadCM2‐SUL. Simulations with changing CO 2 alone show a widely distributed terrestrial carbon sink of 1.4–3.8 Pg C y −1 during the 1990s, rising to 3.7–8.6 Pg C y −1 a century later. Simulations including climate change show a reduced sink both today (0.6–3.0 Pg C y −1 ) and a century later (0.3–6.6 Pg C y −1 ) as a result of the impacts of climate change on NEP of tropical and southern hemisphere ecosystems. In all models, the rate of increase of NEP begins to level off around 2030 as a consequence of the ‘diminishing return’ of physiological CO 2 effects at high CO 2 concentrations. Four out of the six models show a further, climate‐induced decline in NEP resulting from increased heterotrophic respiration and declining tropical NPP after 2050. Changes in vegetation structure influence the magnitude and spatial pattern of the carbon sink and, in combination with changing climate, also freshwater availability (runoff). It is shown that these changes, once set in motion, would continue to evolve for at least a century even if atmospheric CO 2 concentration and climate could be instantaneously stabilized. The results should be considered illustrative in the sense that the choice of CO 2 concentration scenario was arbitrary and only one climate model scenario was used. However, the results serve to indicate a range of possible biospheric responses to CO 2 and climate change. They reveal major uncertainties about the response of NEP to climate change resulting, primarily, from differences in the way that modelled global NPP responds to a changing climate. The simulations illustrate, however, that the magnitude of possible biospheric influences on the carbon balance requires that this factor is taken into account for future scenarios of atmospheric CO 2 and climate change.

We present new reference values for the NIST SRM 610–617 glasses following ISO guidelines and the International Association of Geoanalysts’ protocol. Uncertainties at the 95% confidence level (CL) have been determined for bulk‐ and micro‐analytical purposes. In contrast to former compilation procedures, this approach delivers data that consider present‐day requirements of data quality. New analytical data and the nearly complete data set of the GeoReM database were used for this study. Data quality was checked by the application of the Horwitz function and by a careful investigation of analytical procedures. We have determined quantitatively possible element inhomogeneities using different test portion masses of 1, 0.1 and 0.02 μg. Although avoiding the rim region of the glass wafers, we found moderate inhomogeneities of several chalcophile/siderophile elements and gross inhomogeneities of Ni, Se, Pd and Pt at small test portion masses. The extent of inhomogeneity was included in the determination of uncertainties. While the new reference values agree with the NIST certified values with the one exception of Mn in SRM 610, they typically differ by as much as 10% from the Pearce et al. (1997) values in current use. In a few cases (P, S, Cl, Ta, Re) the discrepancies are even higher.

Abstract. Accurate assessment of anthropogenic carbon dioxide (CO2) emissions andtheir redistribution among the atmosphere, ocean, and terrestrial biospherein a changing climate is critical to better understand the global carboncycle, support the development of climate policies, and project futureclimate change. Here we describe and synthesize data sets and methodologies toquantify the five major components of the global carbon budget and theiruncertainties. Fossil CO2 emissions (EFOS) are based on energystatistics and cement production data, while emissions from land-use change(ELUC), mainly deforestation, are based on land use and land-use changedata and bookkeeping models. Atmospheric CO2 concentration is measureddirectly, and its growth rate (GATM) is computed from the annualchanges in concentration. The ocean CO2 sink (SOCEAN) is estimatedwith global ocean biogeochemistry models and observation-baseddata products. The terrestrial CO2 sink (SLAND) is estimated withdynamic global vegetation models. The resulting carbon budget imbalance(BIM), the difference between the estimated total emissions and theestimated changes in the atmosphere, ocean, and terrestrial biosphere, is ameasure of imperfect data and understanding of the contemporary carboncycle. All uncertainties are reported as ±1σ. For the year 2021, EFOS increased by 5.1 % relative to 2020, withfossil emissions at 10.1 ± 0.5 GtC yr−1 (9.9 ± 0.5 GtC yr−1 when the cement carbonation sink is included), and ELUC was 1.1 ± 0.7 GtC yr−1, for a total anthropogenic CO2 emission(including the cement carbonation sink) of 10.9 ± 0.8 GtC yr−1(40.0 ± 2.9 GtCO2). Also, for 2021, GATM was 5.2 ± 0.2 GtC yr−1 (2.5 ± 0.1 ppm yr−1), SOCEAN was 2.9 ± 0.4 GtC yr−1, and SLAND was 3.5 ± 0.9 GtC yr−1, with aBIM of −0.6 GtC yr−1 (i.e. the total estimated sources were too low orsinks were too high). The global atmospheric CO2 concentration averaged over2021 reached 414.71 ± 0.1 ppm. Preliminary data for 2022 suggest anincrease in EFOS relative to 2021 of +1.0 % (0.1 % to 1.9 %)globally and atmospheric CO2 concentration reaching 417.2 ppm, morethan 50 % above pre-industrial levels (around 278 ppm). Overall, the meanand trend in the components of the global carbon budget are consistentlyestimated over the period 1959–2021, but discrepancies of up to 1 GtC yr−1 persist for the representation of annual to semi-decadalvariability in CO2 fluxes. Comparison of estimates from multipleapproaches and observations shows (1) a persistent large uncertainty in theestimate of land-use change emissions, (2) a low agreement between thedifferent methods on the magnitude of the land CO2 flux in the northernextratropics, and (3) a discrepancy between the different methods on thestrength of the ocean sink over the last decade. This living data updatedocuments changes in the methods and data sets used in this new globalcarbon budget and the progress in understanding of the global carbon cyclecompared with previous publications of this data set. The data presented inthis work are available at https://doi.org/10.18160/GCP-2022 (Friedlingstein et al., 2022b).

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

Abstract. Accurate assessment of anthropogenic carbon dioxide (CO2) emissions andtheir redistribution among the atmosphere, ocean, and terrestrial biosphere– the “global carbon budget” – is important to better understand theglobal carbon cycle, support the development of climate policies, andproject future climate change. Here we describe data sets and methodology toquantify the five major components of the global carbon budget and theiruncertainties. Fossil CO2 emissions (EFF) are based on energystatistics and cement production data, while emissions from land use change(ELUC), mainly deforestation, are based on land use and land use changedata and bookkeeping models. Atmospheric CO2 concentration is measureddirectly and its growth rate (GATM) is computed from the annual changesin concentration. The ocean CO2 sink (SOCEAN) and terrestrialCO2 sink (SLAND) are estimated with global process modelsconstrained by observations. The resulting carbon budget imbalance(BIM), the difference between the estimated total emissions and theestimated changes in the atmosphere, ocean, and terrestrial biosphere, is ameasure of imperfect data and understanding of the contemporary carboncycle. All uncertainties are reported as ±1σ. For the lastdecade available (2009–2018), EFF was 9.5±0.5 GtC yr−1,ELUC 1.5±0.7 GtC yr−1, GATM 4.9±0.02 GtC yr−1 (2.3±0.01 ppm yr−1), SOCEAN 2.5±0.6 GtC yr−1, and SLAND 3.2±0.6 GtC yr−1, with a budgetimbalance BIM of 0.4 GtC yr−1 indicating overestimated emissionsand/or underestimated sinks. For the year 2018 alone, the growth in EFF wasabout 2.1 % and fossil emissions increased to 10.0±0.5 GtC yr−1, reaching 10 GtC yr−1 for the first time in history,ELUC was 1.5±0.7 GtC yr−1, for total anthropogenicCO2 emissions of 11.5±0.9 GtC yr−1 (42.5±3.3 GtCO2). Also for 2018, GATM was 5.1±0.2 GtC yr−1 (2.4±0.1 ppm yr−1), SOCEAN was 2.6±0.6 GtC yr−1, and SLAND was 3.5±0.7 GtC yr−1, with a BIM of 0.3 GtC. The global atmospheric CO2 concentration reached 407.38±0.1 ppm averaged over 2018. For 2019, preliminary data for the first 6–10 months indicate a reduced growth in EFF of +0.6 % (range of−0.2 % to 1.5 %) based on national emissions projections for China, theUSA, the EU, and India and projections of gross domestic product correctedfor recent changes in the carbon intensity of the economy for the rest ofthe world. Overall, the mean and trend in the five components of the globalcarbon budget are consistently estimated over the period 1959–2018, butdiscrepancies of up to 1 GtC yr−1 persist for the representation ofsemi-decadal variability in CO2 fluxes. A detailed comparison amongindividual estimates and the introduction of a broad range of observationsshows (1) no consensus in the mean and trend in land use change emissionsover the last decade, (2) a persistent low agreement between the differentmethods on the magnitude of the land CO2 flux in the northernextra-tropics, and (3) an apparent underestimation of the CO2variability by ocean models outside the tropics. This living data updatedocuments changes in the methods and data sets used in this new globalcarbon budget and the progress in understanding of the global carbon cyclecompared with previous publications of this data set (Le Quéré etal., 2018a, b, 2016, 2015a, b, 2014, 2013). The data generated bythis work are available at https://doi.org/10.18160/gcp-2019 (Friedlingsteinet al., 2019).