Global environmental change is rapidly altering the dynamics of terrestrial vegetation, with consequences for the functioning of the Earth system and provision of ecosystem services 1,2 . Yet how global vegetation is responding to the changing environment is not well established. Here we use three long-term satellite leaf area index (LAI) records and ten global ecosystem models to investigate four key drivers of LAI trends during 1982-2009. We show a persistent and widespread increase of growing season integrated LAI (greening) over 25% to 50% of the global vegetated area, whereas less than 4% of the globe shows decreasing LAI (browning). Factorial simulations with multiple global ecosystem models suggest that CO 2 fertilization e ects explain 70% of the observed greening trend, followed by nitrogen deposition (9%), climate change (8%) and land cover change (LCC) (4%). CO 2 fertilization e ects explain most of the greening trends in the tropics, whereas climate change resulted in greening of the high latitudes and the Tibetan Plateau. LCC contributed most to the regional greening observed in southeast China and the eastern United States. The regional e ects of unexplained factors suggest that the next generation of ecosystem models will need to explore the impacts of forest demography, di erences in regional management intensities for cropland and pastures, and other emerging productivity constraints such as phosphorus availability.Changes in vegetation greenness have been reported at regional and continental scales on the basis of forest inventory and satellite measurements 3-8 . Long-term changes in vegetation greenness are driven by multiple interacting biogeochemical drivers and land-use effects 9 . Biogeochemical drivers include the fertilization effects of elevated atmospheric CO 2 concentration (eCO 2 ), regional climate change (temperature, precipitation and radiation), and varying rates of nitrogen deposition. Land-use-related drivers involve changes in land cover and in land management intensity, including fertilization, irrigation, forestry and grazing 10 . None of these driving factors can be considered in isolation, given their strong interactions with one another. Previously, a few studies had investigated the drivers of global greenness trends 6,7,11 , with a limited number of models and satellite observations, which prevented an appropriate quantification of uncertainties 12 .Here, we investigate trends of leaf area index (LAI) and their drivers for the period 1982 to 2009 using three remotely sensed data sets (GIMMS3g, GLASS and GLOMAP) and outputs from ten ecosystem models run at global extent (see Supplementary Information). We use the growing season integrated leaf area index (hereafter, LAI; Methods) as the variable of our study. We first analyse global and regional LAI trends for the study period and differences between the three data sets. Using modelling results, we then quantify the contributions of CO 2 fertilization, climatic factors, nitrogen deposition and LCC to the observed trends...

Abstract. Accurate assessment of anthropogenic carbon dioxide (CO2) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere in a changing climate – the “global carbon budget” – is important to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe and synthesize data sets and methodology to quantify the five major components of the global carbon budget and their uncertainties. Fossil CO2 emissions (EFOS) are based on energy statistics and cement production data, while emissions from land-use change (ELUC), mainly deforestation, are based on land use and land-use change data and bookkeeping models. Atmospheric CO2 concentration is measured directly and its growth rate (GATM) is computed from the annual changes in concentration. The ocean CO2 sink (SOCEAN) and terrestrial CO2 sink (SLAND) are estimated with global process models constrained by observations. The resulting carbon budget imbalance (BIM), the difference between the estimated total emissions and the estimated changes in the atmosphere, ocean, and terrestrial biosphere, is a measure of imperfect data and understanding of the contemporary carbon cycle. All uncertainties are reported as ±1σ. For the last decade 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), and ELUC 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 budget imbalance BIM of −0.1 GtC yr−1 indicating a near balance between estimated sources and sinks over the last decade. For the year 2019 alone, the growth in EFOS was only about 0.1 % with fossil emissions increasing to 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 was 5.4 ± 0.2 GtC yr−1 (2.5 ± 0.1 ppm yr−1), SOCEAN was 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 CO2 concentration reached 409.85 ± 0.1 ppm averaged over 2019. Preliminary data for 2020, accounting for the COVID-19-induced changes in emissions, suggest a decrease in EFOS relative to 2019 of about −7 % (median estimate) based on individual estimates from four studies of −6 %, −7 %, −7 % (−3 % to −11 %), and −13 %. Overall, the mean and trend in the components of the global carbon budget are consistently estimated over the period 1959–2019, but discrepancies of up to 1 GtC yr−1 persist for the representation of semi-decadal variability in CO2 fluxes. Comparison of estimates from diverse approaches and observations shows (1) no consensus in 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 of the land CO2 flux in the northern extra-tropics, and (3) an apparent discrepancy between the different methods for the ocean sink outside the tropics, particularly in the Southern Ocean. This living data update documents changes in the methods and data sets used in this new global carbon budget and the progress in understanding of the global carbon cycle compared 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).

NATURE CLIMATE CHANGE | ADVANCE ONLINE PUBLICATION | www.nature.com/natureclimatechange 1 D espite two decades of effort to curb emissions of CO 2 and other greenhouse gases (GHGs), emissions grew faster during the 2000s than in the 1990s 1 , and by 2010 had reached ~50 Gt CO 2 equivalent (CO 2 eq) yr −1 (refs 2,3). The continuing rise in emissions is a growing challenge for meeting the international goal of limiting warming to less than 2 °C relative to the pre-industrial era, particularly without stringent climate policies to decrease emissions in the near future 2-4 . As negative emissions technologies (NETs) seem ever more necessary 3,[5][6][7][8][9][10] To have a >50% chance of limiting warming below 2 °C, most recent scenarios from integrated assessment models (IAMs) require large-scale deployment of negative emissions technologies (NETs). These are technologies that result in the net removal of greenhouse gases from the atmosphere. We quantify potential global impacts of the different NETs on various factors (such as land, greenhouse gas emissions, water, albedo, nutrients and energy) to determine the biophysical limits to, and economic costs of, their widespread application. Resource implications vary between technologies and need to be satisfactorily addressed if NETs are to have a significant role in achieving climate goals.options, to be able to decide which pathways are most desirable for dealing with climate change.There are distinct classes of NETs, such as: (1) bioenergy with carbon capture and storage (BECCS) 11,12 ; (2) direct air capture of CO 2 from ambient air by engineered chemical reactions (DAC) 13,14 ; (3) enhanced weathering of minerals (EW) 15 , where natural weathering to remove CO 2 from the atmosphere is accelerated and the products stored in soils, or buried in land or deep ocean [16][17][18][19] ; (4) afforestation and reforestation (AR) to fix atmospheric carbon in biomass and soils [20][21][22] ; (5) manipulation of carbon uptake by the ocean, either

The growth rate of atmospheric carbon dioxide (CO2) concentrations since industrialization is characterized by large interannual variability, mostly resulting from variability in CO2 uptake by terrestrial ecosystems (typically termed carbon sink). However, the contributions of regional ecosystems to that variability are not well known. Using an ensemble of ecosystem and land-surface models and an empirical observation-based product of global gross primary production, we show that the mean sink, trend, and interannual variability in CO2 uptake by terrestrial ecosystems are dominated by distinct biogeographic regions. Whereas the mean sink is dominated by highly productive lands (mainly tropical forests), the trend and interannual variability of the sink are dominated by semi-arid ecosystems whose carbon balance is strongly associated with circulation-driven variations in both precipitation and temperature.

Accurate assessment of anthropogenic carbon dioxide (CO 2 ) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere -the "global carbon budget" -is important to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe data sets and methodology to quantify the five major components of the global carbon budget and their uncertainties. Fossil CO 2 emissions (E FF ) are based on energy statistics and cement production data, while emissions from land use and land-use change (E LUC ), mainly deforestation, are based on land use and land-use change data and bookkeeping models. Atmospheric CO 2 concentration is measured directly and its growth rate (G ATM ) is computed from the annual changes in concentration. The ocean CO 2 sink (S OCEAN ) and terrestrial CO 2 sink (S LAND ) are estimated with global process models constrained by observations. The resulting carbon budget imbalance (B IM ), the difference between the estimated total emissions and the estimated changes in the atmosphere, ocean, and terrestrial biosphere, is a measure of imperfect data and understanding of the contemporary carbon cycle. All uncertainties are reported as ±1σ . For the last decade available (2008-2017), E FF was 9.4 ± 0.5 GtC yr −1 , E LUC 1.5 ± 0.7 GtC yr −1 , G ATM 4.7 ± 0.02 GtC yr −1 , S OCEAN 2.4 ± 0.5 GtC yr −1 , and S LAND 3.2 ± 0.8 GtC yr −1 , with a budget imbalance B IM of 0.5 GtC yr −1 indicating overestimated emissions and/or underestimated sinks. For the year 2017 alone, the growth in E FF was about 1.6 % and emissions increased to 9.9 ± 0.5 GtC yr −1 . Also for 2017, E LUC was 1.4 ± 0.7 GtC yr −1 , G ATM was 4.6 ± 0.2 GtC yr −1 , S OCEAN was 2.5 ± 0.5 GtC yr −1 , and S LAND was 3.8 ± 0.8 GtC yr −1 , with a B IM of 0.3 GtC. The global atmospheric CO 2 concentration reached 405.0 ± 0.1 ppm averaged over 2017. For 2018, preliminary data for the first 6-9 months indicate a renewed growth in E FF of +2.7 % (range of 1.8 % to 3.7 %) based on national emission projections for China, the US, the EU, and India and projections of gross domestic product corrected for recent changes in the carbon intensity of the economy for the rest of the world. The analysis presented here shows that the mean and trend in the five components of the global carbon budget are consistently estimated over the period of 1959-2017, but discrepancies of up to 1 GtC yr −1 persist for the representation of semi-decadal variability in CO 2 fluxes. A detailed comparison among individual estimates and the introduction of a broad range of observations show (1) no consensus in the mean and trend in land-use change emissions, (2) a persistent low agreement among the different methods on the magnitude of the land CO 2 flux in the northern extra-tropics, and (3) an apparent underestimation of the CO 2 variability by ocean models, originating outside the tropics. This living data update documents changes in the methods and data sets used in this new global c...

Abstract. Accurate assessment of anthropogenic carbon dioxide (CO2) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere – the “global carbon budget” – is important to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe data sets and methodology to quantify the five major components of the global carbon budget and their uncertainties. Fossil CO2 emissions (EFF) are based on energy statistics and cement production data, while emissions from land use change (ELUC), mainly deforestation, are based on land use and land use change data and bookkeeping models. Atmospheric CO2 concentration is measured directly and its growth rate (GATM) is computed from the annual changes in concentration. The ocean CO2 sink (SOCEAN) and terrestrial CO2 sink (SLAND) are estimated with global process models constrained by observations. The resulting carbon budget imbalance (BIM), the difference between the estimated total emissions and the estimated changes in the atmosphere, ocean, and terrestrial biosphere, is a measure of imperfect data and understanding of the contemporary carbon cycle. All uncertainties are reported as ±1σ. For the last decade 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 budget imbalance BIM of 0.4 GtC yr−1 indicating overestimated emissions and/or underestimated sinks. For the year 2018 alone, the growth in EFF was about 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 anthropogenic CO2 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, the USA, the EU, and India and projections of gross domestic product corrected for recent changes in the carbon intensity of the economy for the rest of the world. Overall, the mean and trend in the five components of the global carbon budget are consistently estimated over the period 1959–2018, but discrepancies of up to 1 GtC yr−1 persist for the representation of semi-decadal variability in CO2 fluxes. A detailed comparison among individual estimates and the introduction of a broad range of observations shows (1) no consensus in 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 of the land CO2 flux in the northern extra-tropics, and (3) an apparent underestimation of the CO2 variability by ocean models outside the tropics. This living data update documents changes in the methods and data sets used in this new global carbon budget and the progress in understanding of the global carbon cycle compared with previous publications of this data set (Le Quéré et al., 2018a, b, 2016, 2015a, b, 2014, 2013). The data generated by this work are available at https://doi.org/10.18160/gcp-2019 (Friedlingstein et al., 2019).

Atmospheric vapor pressure deficit (VPD) is a critical variable in determining plant photosynthesis. Synthesis of four global climate datasets reveals a sharp increase of VPD after the late 1990s. In response, the vegetation greening trend indicated by a satellite-derived vegetation index (GIMMS3g), which was evident before the late 1990s, was subsequently stalled or reversed. Terrestrial gross primary production derived from two satellite-based models (revised EC-LUE and MODIS) exhibits persistent and widespread decreases after the late 1990s due to increased VPD, which offset the positive CO2 fertilization effect. Six Earth system models have consistently projected continuous increases of VPD throughout the current century. Our results highlight that the impacts of VPD on vegetation growth should be adequately considered to assess ecosystem responses to future climate conditions.

Abstract. Accurate assessment of anthropogenic carbon dioxide (CO 2 ) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere -the "global carbon budget" -is important to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe data sets and methodology to quantify the five major components of the global carbon budget and their uncertainties. CO 2 emissions from fossil fuels and industry (E FF ) are based on energy statistics and cement production data, respectively, while emissions from land-use change (E LUC ), mainly deforestation, are based on land-cover change data and bookkeeping models. The global atmospheric CO 2 concentration is measured directly and its rate of growth (G ATM ) is computed from the annual changes in concentration. The ocean CO 2 sink (S OCEAN ) and terrestrial CO 2 sink (S LAND ) are estimated with global process models constrained by observations. The resulting carbon budget imbalance (B IM ), the difference between the estimated total emissions and the estimated changes in the atmosphere, ocean, and terrestrial biosphere, is a measure of imperfect data and understanding of the contemporary carbon cycle. All uncertainties are reported as ±1σ . For the last decade available (2007)(2008)(2009)(2010)(2011)(2012)(2013)(2014)(2015)(2016), E FF was 9.4 ± 0.5 GtC yr −1 , E LUC 1.3 ± 0.7 GtC yr −1 , G ATM 4.7 ± 0.1 GtC yr −1 , S OCEAN 2.4 ± 0.5 GtC yr −1 , and S LAND 3.0 ± 0.8 GtC yr −1 , with a budget imbalance B IM of 0.6 GtC yr −1 indicating overestimated emissions and/or underestimated sinks. For year 2016 alone, the growth in E FF was approximately zero and emissions remained at 9.9 ± 0.5 GtC yr −1 . Also for 2016, E LUC was 1.3 ± 0.7 GtC yr −1 , G ATM was 6.1 ± 0.2 GtC yr −1 , S OCEAN was 2.6 ± 0.5 GtC yr −1 , and S LAND was 2.7 ± 1.0 GtC yr −1 , with a small B IM of −0.3 GtC. G ATM continued to be higher in 2016 compared to the past decade (2007)(2008)(2009)(2010)(2011)(2012)(2013)(2014)(2015)(2016), reflecting in part the high fossil emissions and the small S LAND consistent with El Niño conditions. The global atmospheric CO 2 concentration reached 402.8 ± 0.1 ppm averaged over 2016. For 2017, preliminary data for the first 6-9 months indicate a renewed growth in E FF of +2.0 % (range of 0.8 to 3.0 %) based on national emissions projections for China, USA, and India, and projections of gross domestic product (GDP) corrected for recent changes in the carbon intensity of the economy for the rest of the world. This living data update documents changes in the methods and data sets used in this new global carbon budget compared with previous publications of this data set (Le Quéré et al