Max Planck Institute for Meteorology

Abstract. By coordinating the design and distribution of global climate model simulations of the past, current, and future climate, the Coupled Model Intercomparison Project (CMIP) has become one of the foundational elements of climate science. However, the need to address an ever-expanding range of scientific questions arising from more and more research communities has made it necessary to revise the organization of CMIP. After a long and wide community consultation, a new and more federated structure has been put in place. It consists of three major elements: (1) a handful of common experiments, the DECK (Diagnostic, Evaluation and Characterization of Klima) and CMIP historical simulations (1850–near present) that will maintain continuity and help document basic characteristics of models across different phases of CMIP; (2) common standards, coordination, infrastructure, and documentation that will facilitate the distribution of model outputs and the characterization of the model ensemble; and (3) an ensemble of CMIP-Endorsed Model Intercomparison Projects (MIPs) that will be specific to a particular phase of CMIP (now CMIP6) and that will build on the DECK and CMIP historical simulations to address a large range of specific questions and fill the scientific gaps of the previous CMIP phases. The DECK and CMIP historical simulations, together with the use of CMIP data standards, will be the entry cards for models participating in CMIP. Participation in CMIP6-Endorsed MIPs by individual modelling groups will be at their own discretion and will depend on their scientific interests and priorities. With the Grand Science Challenges of the World Climate Research Programme (WCRP) as its scientific backdrop, CMIP6 will address three broad questions: – How does the Earth system respond to forcing? – What are the origins and consequences of systematic model biases? – How can we assess future climate changes given internal climate variability, predictability, and uncertainties in scenarios? This CMIP6 overview paper presents the background and rationale for the new structure of CMIP, provides a detailed description of the DECK and CMIP6 historical simulations, and includes a brief introduction to the 21 CMIP6-Endorsed MIPs.

Abstract ERA‐40 is a re‐analysis of meteorological observations from September 1957 to August 2002 produced by the European Centre for Medium‐Range Weather Forecasts (ECMWF) in collaboration with many institutions. The observing system changed considerably over this re‐analysis period, with assimilable data provided by a succession of satellite‐borne instruments from the 1970s onwards, supplemented by increasing numbers of observations from aircraft, ocean‐buoys and other surface platforms, but with a declining number of radiosonde ascents since the late 1980s. The observations used in ERA‐40 were accumulated from many sources. The first part of this paper describes the data acquisition and the principal changes in data type and coverage over the period. It also describes the data assimilation system used for ERA‐40. This benefited from many of the changes introduced into operational forecasting since the mid‐1990s, when the systems used for the 15‐year ECMWF re‐analysis (ERA‐15) and the National Centers for Environmental Prediction/National Center for Atmospheric Research (NCEP/NCAR) re‐analysis were implemented. Several of the improvements are discussed. General aspects of the production of the analyses are also summarized. A number of results indicative of the overall performance of the data assimilation system, and implicitly of the observing system, are presented and discussed. The comparison of background (short‐range) forecasts and analyses with observations, the consistency of the global mass budget, the magnitude of differences between analysis and background fields and the accuracy of medium‐range forecasts run from the ERA‐40 analyses are illustrated. Several results demonstrate the marked improvement that was made to the observing system for the southern hemisphere in the 1970s, particularly towards the end of the decade. In contrast, the synoptic quality of the analysis for the northern hemisphere is sufficient to provide forecasts that remain skilful well into the medium range for all years. Two particular problems are also examined: excessive precipitation over tropical oceans and a too strong Brewer‐Dobson circulation, both of which are pronounced in later years. Several other aspects of the quality of the re‐analyses revealed by monitoring and validation studies are summarized. Expectations that the ‘second‐generation’ ERA‐40 re‐analysis would provide products that are better than those from the firstgeneration ERA‐15 and NCEP/NCAR re‐analyses are found to have been met in most cases. © Royal Meteorological Society, 2005. The contributions of N. A. Rayner and R. W. Saunders are Crown copyright.

The Synthesis Report (SYR) distils and integrates the findings of the three Working Group contributions to the Fifth Assessment Report (AR5) of the Intergovernmental Panel on Climate Change (IPCC), the most comprehensive assessment of climate change undertaken thus far by the IPCC: Climate Change 2013: The Physical Science Basis; Climate Change 2014: Impacts, Adaptation, and Vulnerability; and Climate Change 2014: Mitigation of Climate Change. The SYR also incorporates the findings of two Special Reports on Renewable Energy Sources and Climate Change Mitigation (2011) and on Managing the Risks of Extreme Events and Disasters to Advance Climate Change Adaptation (2011).

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.

We present a first description and evaluation of GEOS‐CHEM, a global threedimensional (3‐D) model of tropospheric chemistry driven by assimilated meteorological observations from the Goddard Earth Observing System (GEOS) of the NASA Data Assimilation Office (DAO). The model is applied to a 1‐year simulation of tropospheric ozone‐NO x ‐hydrocarbon chemistry for 1994, and is evaluated with observations both for 1994 and for other years. It reproduces usually to within 10 ppb the concentrations of ozone observed from the worldwide ozonesonde data network. It simulates correctly the seasonal phases and amplitudes of ozone concentrations for different regions and altitudes, but tends to underestimate the seasonal amplitude at northern midlatitudes. Observed concentrations of NO and peroxyacetylnitrate (PAN) observed in aircraft campaigns are generally reproduced to within a factor of 2 and often much better. Concentrations of HNO 3 in the remote troposphere are overestimated typically by a factor of 2–3, a common problem in global models that may reflect a combination of insufficient precipitation scavenging and gas‐aerosol partitioning not resolved by the model. The model yields an atmospheric lifetime of methylchloroform (proxy for global OH) of 5.1 years, as compared to a best estimate from observations of 5.5 +/− 0.8 years, and simulates H 2 O 2 concentrations observed from aircraft with significant regional disagreements but no global bias. The OH concentrations are ∼20% higher than in our previous global 3‐D model which included an UV‐absorbing aerosol. Concentrations of CO tend to be underestimated by the model, often by 10–30 ppb, which could reflect a combination of excessive OH (a 20% decrease in model OH could be accommodated by the methylchloroform constraint) and an underestimate of CO sources (particularly biogenic). The model underestimates observed acetone concentrations over the South Pacific in fall by a factor of 3; a missing source from the ocean may be implicated.

A new high-resolution regional climate change ensemble has been established for Europe within the World Climate Research Program Coordinated Regional Downscaling Experiment (EURO-CORDEX) initiative. The first set of simulations with a horizontal resolution of 12.5 km was completed for the new emission scenarios RCP4.5 and RCP8.5 with more simulations expected to follow. The aim of this paper is to present this data set to the different communities active in regional climate modelling, impact assessment and adaptation. The EURO-CORDEX ensemble results have been compared to the SRES A1B simulation results achieved within the ENSEMBLES project. The large-scale patterns of changes in mean temperature and precipitation are similar in all three scenarios, but they differ in regional details, which can partly be related to the higher resolution in EURO-CORDEX. The results strengthen those obtained in ENSEMBLES, but need further investigations. The analysis of impact indices shows that for RCP8.5, there is a substantially larger change projected for temperature-based indices than for RCP4.5. The difference is less pronounced for precipitation-based indices. Two effects of the increased resolution can be regarded as an added value of regional climate simulations. Regional climate model simulations provide higher daily precipitation intensities, which are completely missing in the global climate model simulations, and they provide a significantly different climate change of daily precipitation intensities resulting in a smoother shift from weak to moderate and high intensities.

Abstract. Aerosols affect the climate system by changing cloud characteristics in many ways. They act as cloud condensation and ice nuclei, they may inhibit freezing and they could have an influence on the hydrological cycle. While the cloud albedo enhancement (Twomey effect) of warm clouds received most attention so far and traditionally is the only indirect aerosol forcing considered in transient climate simulations, here we discuss the multitude of effects. Different approaches how the climatic implications of these aerosol effects can be estimated globally as well as improvements that are needed in global climate models in order to better represent indirect aerosol effects are discussed in this paper.

Understanding and quantifying the global methane (CH4) budget is important for assessing realistic pathways to mitigate climate change. Atmospheric emissions and concentrations of CH4 continue to increase, making CH4 the second most important human-influenced greenhouse gas in terms of climate forcing, after carbon dioxide (CO2). The relative importance of CH4 compared to CO2 depends on its shorter atmospheric\nlifetime, stronger warming potential, and variations in atmospheric growth rate over the past decade, the causes of which are still debated. Two major challenges in reducing uncertainties in the atmospheric growth rate arise from the variety of geographically overlapping CH4 sources and from the destruction of CH4 by short-lived hydroxyl radicals (OH). To address these challenges, we have established a consortium of multidisciplinary scientists under the umbrella of the Global Carbon Project to synthesize and stimulate new research aimed at improving and regularly updating the global methane budget. Following Saunois et al. (2016), we present here the second version of the living review paper dedicated to the decadal methane budget, integrating results of top-down studies (atmospheric observations within an atmospheric inverse-modelling framework) and bottom-up estimates (including process-based models for estimating land surface emissions and atmospheric chemistry, inventories of anthropogenic emissions, and data-driven extrapolations).\nFor the 2008–2017 decade, global methane emissions are estimated by atmospheric inversions (a top-down approach) to be 576 TgCH4 yr-1 (range 550–594, corresponding to the minimum and maximum estimates of the model ensemble). Of this total, 359 TgCH4 yr-1 or 60% is attributed to anthropogenic sources, that is emissions caused by direct human activity (i.e. anthropogenic emissions; range 336–376 TgCH4 yr-1 or 50 %–65 %). The mean annual total emission for the new decade (2008–2017) is 29 TgCH4 yr-1 larger than our estimate for the previous decade (2000–2009), and 24 TgCH4 yr-1 larger than the one reported in the previous budget for 2003–2012 (Saunois et al., 2016). Since 2012, global CH4 emissions have been tracking the warmest scenarios assessed by the Intergovernmental Panel on Climate Change. Bottom-up methods suggest almost 30% larger global emissions (737 TgCH4 yr-1, range 594–881) than top-down inversion methods. Indeed, bottom-up estimates for natural sources such as natural wetlands, other inland water systems, and geological sources are higher than top-down estimates. The atmospheric constraints on the top-down budget suggest that at least some of these bottom-up emissions are overestimated. The latitudinal distribution of atmospheric observation-based emissions indicates a predominance of tropical emissions (∼65% of the global budget, <30◦N) compared to mid-latitudes (∼30 %, 30–60◦ N) and high northern latitudes (∼4 %, 60–90◦N). The most important source of uncertainty in the methane budget is attributable to natural emissions, especially those from wetlands and other inland waters.\nSome of our global source estimates are smaller than those in previously published budgets (Saunois et al., 2016; Kirschke et al., 2013). In particular wetland emissions are about 35 TgCH4 yr-1 lower due to improved partition wetlands and other inland waters. Emissions from geological sources and wild animals are also found to be smaller by 7 TgCH4 yr-1 by 8 TgCH4 yr-1, respectively. However, the overall discrepancy between bottom-up and top-down estimates has been reduced by only 5% compared to Saunois et al. (2016), due to a higher estimate of emissions from inland waters, highlighting the need for more detailed research on emissions factors. Priorities for improving the methane budget include (i) a global, high-resolution map of water-saturated soils and inundated areas emitting methane based on a robust classification of different types of emitting habitats; (ii) further development of process-based models for inland-water emissions; (iii) intensification of methane observations at local scales (e.g., FLUXNET-CH4 measurements) and urban-scale monitoring to constrain bottom-up land surface models, and at regional scales (surface networks and satellites) to constrain atmospheric inversions; (iv) improvements of transport models and the representation of photochemical sinks in top-down inversions; and (v) development of a 3D variational inversion system using isotopic and/or co-emitted species such as ethane to improve source partitioning.\nThe data presented here can be downloaded from https://doi.org/10.18160/GCP-CH4-2019 (Saunois et al.,\n2020) and from the Global Carbon Project

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

Abstract. The Model for Ozone and Related chemical Tracers, version 4 (MOZART-4) is an offline global chemical transport model particularly suited for studies of the troposphere. The updates of the model from its previous version MOZART-2 are described, including an expansion of the chemical mechanism to include more detailed hydrocarbon chemistry and bulk aerosols. Online calculations of a number of processes, such as dry deposition, emissions of isoprene and monoterpenes and photolysis frequencies, are now included. Results from an eight-year simulation (2000–2007) are presented and evaluated. The MOZART-4 source code and standard input files are available for download from the NCAR Community Data Portal (http://cdp.ucar.edu).

Addresses both fundamental and applied aspects of ocean waves including the use of wave observations made from satellites. More specifically it describes the WAM model, its scientific basis, its actual implementation, and its many applications. This model has been developed by an international group (the Wave Modelling group), and is based on a detailed physical description of air/sea interactions. It is widely used for wave forecasting for meteorological and oceanographic purposes. The three sections of the volume describe the basic statistical theory and the relevant physical processes; the numerical model and its global and regional applications; and satellite observations, their interpretation and use in data assimilation. Written by leading experts, it is a comprehensive guide and reference for researchers and advanced students in physical oceanography, meteorology, fluid dynamics, coastal engineering and physics.

Water scarcity severely impairs food security and economic prosperity in many countries today. Expected future population changes will, in many countries as well as globally, increase the pressure on available water resources. On the supply side, renewable water resources will be affected by projected changes in precipitation patterns, temperature, and other climate variables. Here we use a large ensemble of global hydrological models (GHMs) forced by five global climate models and the latest greenhouse-gas concentration scenarios (Representative Concentration Pathways) to synthesize the current knowledge about climate change impacts on water resources. We show that climate change is likely to exacerbate regional and global water scarcity considerably. In particular, the ensemble average projects that a global warming of 2 °C above present (approximately 2.7 °C above preindustrial) will confront an additional approximate 15% of the global population with a severe decrease in water resources and will increase the number of people living under absolute water scarcity (<500 m(3) per capita per year) by another 40% (according to some models, more than 100%) compared with the effect of population growth alone. For some indicators of moderate impacts, the steepest increase is seen between the present day and 2 °C, whereas indicators of very severe impacts increase unabated beyond 2 °C. At the same time, the study highlights large uncertainties associated with these estimates, with both global climate models and GHMs contributing to the spread. GHM uncertainty is particularly dominant in many regions affected by declining water resources, suggesting a high potential for improved water resource projections through hydrological model development.

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

The new Max‐Planck‐Institute Earth System Model (MPI‐ESM) is used in the Coupled Model Intercomparison Project phase 5 (CMIP5) in a series of climate change experiments for either idealized CO 2 ‐only forcing or forcings based on observations and the Representative Concentration Pathway (RCP) scenarios. The paper gives an overview of the model configurations, experiments related forcings, and initialization procedures and presents results for the simulated changes in climate and carbon cycle. It is found that the climate feedback depends on the global warming and possibly the forcing history. The global warming from climatological 1850 conditions to 2080–2100 ranges from 1.5°C under the RCP2.6 scenario to 4.4°C under the RCP8.5 scenario. Over this range, the patterns of temperature and precipitation change are nearly independent of the global warming. The model shows a tendency to reduce the ocean heat uptake efficiency toward a warmer climate, and hence acceleration in warming in the later years. The precipitation sensitivity can be as high as 2.5% K −1 if the CO 2 concentration is constant, or as small as 1.6% K −1 , if the CO 2 concentration is increasing. The oceanic uptake of anthropogenic carbon increases over time in all scenarios, being smallest in the experiment forced by RCP2.6 and largest in that for RCP8.5. The land also serves as a net carbon sink in all scenarios, predominantly in boreal regions. The strong tropical carbon sources found in the RCP2.6 and RCP8.5 experiments are almost absent in the RCP4.5 experiment, which can be explained by reforestation in the RCP4.5 scenario.

ABSTRACT Thawing permafrost and the resulting microbial decomposition of previously frozen organic carbon (C) is one of the most significant potential feedbacks from terrestrial ecosystems to the atmosphere in a changing climate. In this article we present an overview of the global permafrost C pool and of the processes that might transfer this C into the atmosphere, as well as the associated ecosystem changes that occur with thawing. We show that accounting for C stored deep in the permafrost more than doubles previous high-latitude inventory estimates, with this new estimate equivalent to twice the atmospheric C pool. The thawing of permafrost with warming occurs both gradually and catastrophically, exposing organic C to microbial decomposition. Other aspects of ecosystem dynamics can be altered by climate change along with thawing permafrost, such as growing season length, plant growth rates and species composition, and ecosystem energy exchange. However, these processes do not appear to be able to compensate for C release from thawing permafrost, making it likely that the net effect of widespread permafrost thawing will be a positive feedback to a warming climate.

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

Abstract. Ocean ecosystems are increasingly stressed by human-induced changes of their physical, chemical and biological environment. Among these changes, warming, acidification, deoxygenation and changes in primary productivity by marine phytoplankton can be considered as four of the major stressors of open ocean ecosystems. Due to rising atmospheric CO2 in the coming decades, these changes will be amplified. Here, we use the most recent simulations performed in the framework of the Coupled Model Intercomparison Project 5 to assess how these stressors may evolve over the course of the 21st century. The 10 Earth system models used here project similar trends in ocean warming, acidification, deoxygenation and reduced primary productivity for each of the IPCC's representative concentration pathways (RCPs) over the 21st century. For the "business-as-usual" scenario RCP8.5, the model-mean changes in the 2090s (compared to the 1990s) for sea surface temperature, sea surface pH, global O2 content and integrated primary productivity amount to &amp;amp;plus;2.73 (±0.72) °C, −0.33 (±0.003) pH unit, −3.45 (±0.44)% and −8.6 (±7.9)%, respectively. For the high mitigation scenario RCP2.6, corresponding changes are +0.71 (±0.45) °C, −0.07 (±0.001) pH unit, −1.81 (±0.31)% and −2.0 (±4.1)%, respectively, illustrating the effectiveness of extreme mitigation strategies. Although these stressors operate globally, they display distinct regional patterns and thus do not change coincidentally. Large decreases in O2 and in pH are simulated in global ocean intermediate and mode waters, whereas large reductions in primary production are simulated in the tropics and in the North Atlantic. Although temperature and pH projections are robust across models, the same does not hold for projections of subsurface O2 concentrations in the tropics and global and regional changes in net primary productivity. These high uncertainties in projections of primary productivity and subsurface oxygen prompt us to continue inter-model comparisons to understand these model differences, while calling for caution when using the CMIP5 models to force regional impact models.

Abstract. Accurate assessment of anthropogenic carbon dioxide(CO2) emissions and their redistribution among the atmosphere,ocean, and terrestrial biosphere – the “global carbon budget” – isimportant to better understand the global carbon cycle, support thedevelopment of climate policies, and project future climate change. Here wedescribe data sets and methodology to quantify the five major components ofthe global carbon budget and their uncertainties. Fossil CO2emissions (EFF) are based on energy statistics and cementproduction data, while emissions from land use and land-use change (ELUC),mainly deforestation, are based on land use and land-use change data andbookkeeping models. Atmospheric CO2 concentration is measureddirectly and its growth rate (GATM) is computed from the annualchanges in concentration. The ocean CO2 sink (SOCEAN)and terrestrial CO2 sink (SLAND) are estimated withglobal process models constrained by observations. The resulting carbonbudget imbalance (BIM), the difference between the estimatedtotal emissions and the estimated changes in the atmosphere, ocean, andterrestrial biosphere, is a measure of imperfect data and understanding ofthe contemporary carbon cycle. All uncertainties are reported as ±1σ. For the last decade available (2008–2017), EFF was9.4±0.5 GtC yr−1, ELUC 1.5±0.7 GtC yr−1, GATM 4.7±0.02 GtC yr−1,SOCEAN 2.4±0.5 GtC yr−1, and SLAND 3.2±0.8 GtC yr−1, with a budget imbalance BIM of0.5 GtC yr−1 indicating overestimated emissions and/or underestimatedsinks. For the year 2017 alone, the growth in EFF was about 1.6 %and emissions increased to 9.9±0.5 GtC yr−1. Also for 2017,ELUC was 1.4±0.7 GtC yr−1, GATM was 4.6±0.2 GtC yr−1, SOCEAN was 2.5±0.5 GtC yr−1, and SLAND was 3.8±0.8 GtC yr−1,with a BIM of 0.3 GtC. The global atmosphericCO2 concentration reached 405.0±0.1 ppm averaged over 2017.For 2018, preliminary data for the first 6–9 months indicate a renewedgrowth in EFF of +2.7 % (range of 1.8 % to 3.7 %) basedon national emission projections for China, the US, the EU, and India andprojections of gross domestic product corrected for recent changes in thecarbon intensity of the economy for the rest of the world. The analysispresented here shows that the mean and trend in the five components of theglobal carbon budget are consistently estimated over the period of 1959–2017,but discrepancies of up to 1 GtC yr−1 persist for the representationof semi-decadal variability in CO2 fluxes. A detailed comparisonamong individual estimates and the introduction of a broad range ofobservations show (1) no consensus in the mean and trend in land-use changeemissions, (2) a persistent low agreement among the different methods onthe magnitude of the land CO2 flux in the northern extra-tropics,and (3) an apparent underestimation of the CO2 variability by oceanmodels, originating outside the tropics. This living data update documentschanges in the methods and data sets used in this new global carbon budgetand the progress in understanding the global carbon cycle compared withprevious publications of this data set (Le Quéré et al., 2018, 2016,2015a, b, 2014, 2013). All results presented here can be downloaded fromhttps://doi.org/10.18160/GCP-2018.

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 (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 datasets and methodology 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 firsttime, an approach is shown to reconcile the difference in our ELUCestimate with the one from national greenhouse gas inventories, supportingthe assessment of collective countries' climate progress. For the year 2020, EFOS declined by 5.4 % relative to 2019, withfossil emissions at 9.5 ± 0.5 GtC yr−1 (9.3 ± 0.5 GtC yr−1 when the cement carbonation sink is included), and ELUC was 0.9 ± 0.7 GtC yr−1, for a total anthropogenic CO2 emission of10.2 ± 0.8 GtC yr−1 (37.4 ± 2.9 GtCO2). Also, for2020, GATM was 5.0 ± 0.2 GtC yr−1 (2.4 ± 0.1 ppm yr−1), SOCEAN was 3.0 ± 0.4 GtC yr−1, and SLANDwas 2.9 ± 1 GtC yr−1, with a BIM of −0.8 GtC yr−1. Theglobal atmospheric CO2 concentration averaged over 2020 reached 412.45 ± 0.1 ppm. Preliminary data for 2021 suggest a rebound in EFOSrelative to 2020 of +4.8 % (4.2 % to 5.4 %) globally. Overall, the mean and trend in the components of the global carbon budgetare consistently estimated over the period 1959–2020, but discrepancies ofup to 1 GtC yr−1 persist for the representation of annual tosemi-decadal variability in CO2 fluxes. Comparison of estimates frommultiple approaches and observations shows (1) a persistent largeuncertainty in the estimate of land-use changes emissions, (2) a lowagreement between the different methods on the magnitude of the landCO2 flux in the northern extra-tropics, and (3) a discrepancy betweenthe different methods on the strength of the ocean sink over the lastdecade. This living data update documents changes in the methods and datasets used in this new global carbon budget and the progress in understandingof the global carbon cycle compared with previous publications of this dataset (Friedlingstein et al., 2020, 2019; LeQuéré et al., 2018b, a, 2016, 2015b, a, 2014, 2013). Thedata presented in this work are available at https://doi.org/10.18160/gcp-2021 (Friedlingstein et al., 2021).