Climate Warming in Response to Emission Reductions Consistent with the Paris Agreement

Department of Climate and Environment Change, Key Laboratory of Land Surface Pattern and Simulation, Institute of Geographic Sciences and Natural Resources Research, Chinese Academy of Sciences, 11A, Datun Road, Chaoyang District, Beijing 100101, China School of Earth and Ocean Sciences, University of Victoria, 3800 Finnerty Road, Victoria, BC, Canada, V8W 3V6 Department of Climate and Environment Change, Institute of Geographic Sciences and Natural Resources Research, Chinese Academy of Sciences, Beijing, China


Introduction
To limit global mean warming to well below 2 °C, in accordance with the Paris Agreement [1], 190 countries submitted their Intended Nationally Determined Contributions (INDCs), which outline the intended post-2020 emission plans of each country [2].INDCs became the first target of greenhouse gas (GHG) mitigation reached through a bottom-up approach by nationally intended efforts, so it is easier to monitor a level of commitment than before through a top-down system.However, the impacts of these emission-reduction efforts on climate warming are poorly understood and their adequacy to meet the long-term goal of stabilizing the global mean temperature to 1.5 °C or 2 °C above the preindustrial level is still unknown.
To simulate climate response under INDC scenarios, running a full suite of comprehensive Earth system models (ESMs), such as the CMIP5 models (Coupled Model Intercomparison Project), is unrealistic due to the high computational cost, while running only one certain model is not representative of climate response of the Earth system due to potential model biases.Recent studies have shown a near linear relationship between cumulative carbon emissions and temperature change [3][4][5][6][7][8][9][10]; thereby providing a way to evaluate climate response under INDC scenarios without the need of running additional simulations by comprehensive Earth system models.
In the Fifth Assessment Report (AR5) of the Intergovernmental Panel on Climate Change (IPCC) [11], future climate change was projected under a set of Representative Concentration Pathway (RCP) scenarios, using a model ensemble of comprehensive Earth system models (CMIP5) [12] and a reduced-complexity carbon-cycle and climate model (Model for the Assessment of Greenhouse Gas Induced Climate Change (MAGICC)) [13,14].e ratio of temperature increase to cumulative carbon emissions, also referred to as the transient climate response to cumulative carbon emissions (TCRE), is relatively constant over time and is independent of the CO 2 emissions pathway [4,15,16].Expert judgements [4,5,7,10] based on multiple lines of evidence estimate TCRE to be likely between 0.8 °C and 2.5 °C per 1000 GtC (5 to 95%), of cumulative CO 2 emissions.Most of the above studies presented results only for CO 2 emissions, without considering the effects of non-CO 2 forcing.When non-CO 2 GHG forcing is considered, the additional net warming from non-CO 2 forcings contributes to lower levels of emissions allowed to reach the given temperature target [16][17][18].erefore, estimates of the climate response to cumulative carbon emissions provide a useful benchmark for assessing the climate response under different emission scenarios.In this study, we make use of the available data of climate responses to cumulative emissions based on CMIP5 and MAGICC simulations, to estimate the global mean warming in response to the INDC scenarios.is study is structured as follows: Section 2 describes the methods and data sources, Section 3.1 presents the conditional and unconditional INDC committed emissions, Section 3.2 examines climate responses to cumulative emissions in CMIP5 and MAGICC models and presents an estimate of warming under INDC emission scenarios, while Section 4 provides the conclusions and further discussion.

INDC Data.
e INDC dataset is continuously updated, and includes 192 countries (165 INDCs) that submitted their pledges up to 2017 (July).Most countries have reported their composite targets, such as emission targets, energy targets, forest targets, and adaptation targets [2].e emission targets reported by countries vary from absolute emission target (e.g., reducing emissions by a given amount of GtC/year) to emission target relative to the base year level (e.g., reducing emissions back to 2010 or 2000 levels), or emission-reduction target relative to the baseline emission scenario (e.g., reducing emissions compared with the business-as-usual scenario (BAU, 2030 levels)).e base years of emission data of each country were obtained from the UNFCCC national inventories [19].e baseline scenario data of each country were calculated according to the predicted emissions from the Stockholm Environment Institute [20].e reported emission targets of each country were extracted and, subsequently, total cumulative emissions were calculated (see Section 2.2). e 165 INDCs cover 192 countries, which together account for more than 99% of global GHG emissions.In addition, three countries (Libya, Nicaragua, and Syria) still had not submitted INDC reports at the time of this analysis.e INDC emission targets consist of both the unconditional targets (voluntary emission reductions, irrespective of the actions of other countries) and conditional targets (more aggressive mitigation actions if certain conditions are met regarding the provision of finance or technological assistance from developed countries).Supplementary Table S1 shows the INDCs of the 165 countries (up to July 2017) analyzed in this study.

Calculating Cumulative Global INDC Emissions.
Global cumulative INDC emissions were estimated by summing national INDC emissions for each year.National INDC commitments provide the emission level of the pre-2030 period (Figure 1).For the post-2030 period, annual emission pathways were obtained by linear interpolation based on the expected level of development.
e global emissions are assumed to peak in 2030 and then decline steadily (Figure 2).We assume that the continued action of national emission reduction was adopted in the scenario over the 21st century, and the relatively constant decarbonisation rates were followed for the period after 2030.
Cumulative carbon emissions and temperature responses from RCP experiments were used to construct a function of climate responses to cumulative emissions referred to as TCRE all (Figure 3): where TCRE all represents the climate response for all GHGs, both for CO 2 effect and for non-CO 2 GHG effect.ΔI represents the cumulative anthropogenic GHG emissions above the current level in 2012, including CO 2 and other non-CO 2 GHG emissions.All non-CO 2 emissions were converted into a unified unit of CO 2 equivalent emissions, according to the global warming potential of each gas [21].ΔT is the corresponding change of global temperature, subject to decadal smoothing.ese data were obtained from the RCP simulation experiments for CMIP5 and MAGICC models.e warming above the current level (ΔT INDC ) under the INDC scenarios was estimated by the following equation: where ΔI INDC represents cumulative emissions (from year 2012) under the INDC scenarios, which was calculated by summing national INDC emissions for each year.en, the warming level above the preindustrial level was estimated based on the sum of the INDC warming above the current level (ΔT INDC ) and the current warming above the preindustrial level.e current warming in 2012 was estimated to be about 0.85 ± 0.14 °C [22].

2
Advances in Meteorology e spatial pattern of global warming under the INDC scenario was estimated based on the time-slice approach [23][24][25][26], where the spatial state at a specific warming point related to ΔT INDC is taken from the decadal time slices with the respective mean warming for each model separately.
is study uses the spatial output from 12 comprehensive Earth system models from the CMIP5 project [27].ese models include BCC_CSM 1.1 (China), CanESM2 (Canada), CESM1 (BGC) (USA), GFDL-ESM2G (USA), GFDL-ESM2M (USA), INM-CM4 (Russia), IPSL-CM5A-LR (France), IPSL-CM5A-MR (France), IPSL-CM5B-LR (France), MIROC-ESM (Japan), and MPI-ESM-LR and MPI-ESM-MR (Germany).We make use of RCP 4.5, 6.0, and 8.5 scenarios and identify the respective warming patterns corresponding to INDC warming for each model, followed by computing a multimodel average state of the spatial warming pattern based on all model simulations.e simulations are regridded into a common 144 × 72 grid (2.5 °× 2.5 °). e CMIP5 models considered in this study are comprehensive Earth system models (ESMs) with coupled carbon-climate system responses, where terrestrial and ocean carbon-cycle processes are coupled with atmosphere-ocean general circulation models [27,28].In addition to CMIP5 ESMs, we also make use of the MAGICC scenario database.e MAGICC model consists of reduced-complexity carbon-cycle and climate models and emulates the global and annual mean behavior of significantly more complex CMIP5 models [13,14,29].

INDC Emissions Rate and Cumulative INDC Emissions.
Figure 1 shows the emission level of each country in 2030 for unconditional INDC pledges, where emission levels are expressed as a percentage of total global emissions.e countries including China, India, United States, and European Union-28 account for the largest proportions of annual global emissions in 2030 and most underdeveloped countries from Africa, Latin America, and southwest Asia account for a small proportion.
On a global scale, for the unconditional and conditional pledges (Figure 2, red and pink lines), the total INDC  S1).Emission level for each country is expressed as a percentage of total global emissions in 2030.e member states of European Union (EU) are shown as a whole as their emission target is submitted for the whole region.emissions target was about 14.2∼14.9GtC/year for 2030, 14.1∼14.6GtC/year for 2025, and 14.0∼14.3GtC/year for 2020 (Figure 2, solid lines).e rate of annual emission from 2012 to 2030 increased on average by 0.7% per year.e global INDC emissions target under conditional pledge is about 0.7 GtC/year less than that of the unconditional pledges in 2030.In terms of the RCP emissions, the INDC emissions level was intermediate between the emission levels of RCP 4.5 and RCP 8.5.To extend the INDC emissions to 2100, the continued action of emission reduction was adopted in the scenario over the 21st century.We assumed global emissions peaked in 2030, as this condition is essential for the control of warming to meet long-term targets [30].In 2050, the estimated GHG emissions were about 11.3∼14.0GtC/year, while in 2100 they were about 6.4∼9.0GtC/year.From 2030 to 2100, the INDC emission level was intermediate between the emission levels of RCP 4.5 and RCP 6.0.
Cumulative emissions for the unconditional and conditional pledges are estimated to be 263∼270 GtC (for the 2012 to 2030 period) and 940∼1120 GtC (for the 2012 to 2100 period).e INDC cumulative emissions during 2012-2100 were higher than those of RCP 2.6 by about 400∼580 GtC but lower than those of RCP 6.0 and RCP 8.5 by about 380∼560 GtC and 1380∼1560 GtC, respectively.

Global Mean Temperature Estimates in Response to INDC Emission Reductions.
Global mean temperature is proportional to cumulative carbon emissions for a range of emission scenarios considered here (Figure 3(a)).e black line in Figure 3(a) shows the historical values and the colored lines are the results for different RCP scenarios.It was found that the relationship between cumulative carbon emissions and temperature increase does not differ much for different RCP scenarios for low warming targets (such as 1.5 °C or 2.0 °C) and was nearly constant for each RCP pathway, with only a small and stable change when cumulative emissions approached 2000 GtC. e results from MAGICC were well aligned with the CMIP5 results for the RCP pathways.erefore, for lower temperature targets such as 1.5 °C and 2.0 °C, this relationship could be used as an approximation of the projection of the climate response to INDC scenarios.We estimated the ratio of a median TCRE all (1) of 2.12 °C per 1000 GtC using CMIP5 results (blue cross in Figure (3b)).
e value was about 2.06 °C per 1000 GtC using MAGICC results (pink cross in Figure 3(b)).e uncertainty range is likely 1.63∼2.59°C per 1000 GtC (5 to 95%).Note that these values apply to cumulative emissions of up to 600 GtC (about 600 GtC have been emitted at the present time).If only CO 2 -induced temperature response is considered to estimate climate warming (gray line in Figure 3(a), based on the 1PCT simulations, where atmospheric CO 2 concentration increases at a rate of 1%•yr −1 using CMIP5 models), TCRE is lower, due to lack of the net warming from non-CO2 forcings that are present in the RCP scenarios.
Estimates of global mean warming under INDC emission scenarios are based on the above definition of the climate response to cumulative emissions (TCRE all ).For the unconditional INDC pledges (Figure 4(a)), the global mean   Advances in Meteorology temperature change in 2030 is projected to be 0.57 °C (median) above the 2012 baseline for the CMIP5 simulations and 0.56 °C (median) for the MAGICC simulations.e likely range is 0.44∼0.70°C (5 to 95%) above the 2012 baseline.Relative to the preindustrial levels, the increase is projected to be 1.42 °C and 1.41 °C for the CMIP5 and MAGICC models, respectively (likely range, 1.29∼1.55°C) (Figure 4(c)).By the end of this century, the global temperature is projected to be 2.36 °C (CMIP5) and 2.30 Correspondingly, for the conditional INDC pledges (Figure 4(b)), the global temperature increase in 2030 is projected to be 0.56 °C above the 2012 level for the CMIP5 simulations and 0.54 °C for the MAGICC simulations (range, 0.43∼0.68°C).In 2100, the global temperature is projected to be 1.98 °C (CMIP5) and 1.93 °C (MAGICC) above the 2012 level (range, 1.53∼2.42°C) and 2.83 °C (CMIP5) and 2.78 °C (MAGICC) above the preindustrial level (Figure 4(c)).
Figures 5(a) and 5(b) show multimodel mean regional patterns of surface temperature changes for unconditional INDC scenario in 2030 and 2100, respectively.e Arctic

Advances in Meteorology
warming is projected to be most profound, exceeding the global average by a factor of three, with about 3∼5 °C Arctic warming in 2030 and 8∼10 °C in 2100 relative to the preindustrial level.e warming in midlatitude is nearly a factor of two greater than the global average in both 2030 and 2100.e south oceans and parts of North Atlantic exhibit lowest warming (Figures 5(a

Discussion and Conclusions
Our study estimated the global mean temperature increase under the INDC commitments in 2030 to range from 1.29 to 1.55 °C (median of 1.41 °C) above the preindustrial level, reaching 2.67∼3.74°C (median of 3.17 °C) in 2100.Our best estimates were within the reported ranges from other studies (e.g., UNEP [31], CI [32], and CAT [33]).Our best estimate for global mean warming is a little higher than that of Rogelj et al. [30] (the median of 2.6∼3.1 °C).We also compared the results of this study with the temperature increase resulting from the RCP scenarios, reported in IPCC AR5 [16,34].Table 1 gives the warming estimates for each scenario in 2030 and 2100.In 2030, the INDC level of warming will be higher than that estimated from RCP 2.6 and 6.0, but lower than that estimated in RCP 8.5.However, the temperature differences between the INDC scenarios and other scenarios are very small (in the order of 0.01 °C).In 2100, the INDC warming will be higher than that estimated form RCP 4.5 and 2.6 scenarios and lower than that estimated from RCP 6.0 and 8.5, but closer to the warming estimated in RCP 4.5 and RCP 6.0 (Table 1).Sanderson et al. [35] proposed a set of idealized emission pathways consistent with reaching the 2 °C temperature target, which showed that if the INDCs for 2030 remain the same as committed, only net zero GHG emissions by 2085 and negative emissions implemented later   Advances in Meteorology in the century are required to limit the global mean warming below the 2 °C warming.However, such rapid emission reductions would be quite difficult to implement [30].
Compared with the long-term temperature-limiting targets, the global mean warming resulting from INDC emission reductions exceeds the 2 °C or 1.5 °C stabilization goals set by UNFCCC [36].erefore, even the conditional INDCs and resulting emission reductions are insufficient to meet the globally agreed long-term goal of limiting global mean warming to well below 2 °C or 1.5 °C.To achieve the desired temperature levels in the future, it is therefore necessary to adjust and strengthen the INDC emissionreduction commitments.
Regional climate change prediction is more relevant for assessing impact-related temperature-controlled goals.e results indicate that by 2030 the high-latitude regions show higher warming level under INDC emission pledges, while by 2100 the significant warming trends will reach across the globe.In particular, the intense warming of the Arctic and other high-latitude land regions may have serious implications for other components of the Earth system (e.g., Arctic summer sea ice retreat, Alpine glaciers loss, and coral reefs bleaching).Future regional extreme climate change should be given more attention.
e results presented in this study are sensitive to future non-CO 2 emission trajectories, which are not fully included in the INDC-committed emission reductions.For example, strong aerosol emissions could have negative effects on warming, even by several tenths of a degree in the short term [37,38].Further research is needed regarding the probability distribution of warming in response to INDC scenarios and its sensitivity to different non-CO 2 emission trajectories.
e global mean temperature under different INDC emission scenarios was estimated using estimates of TCRE based on CMIP5 and MAGICC responses.e advantage of this approach is the avoidance of high costs of running all comprehensive ESMs to simulate the climate response under INDC scenarios.is approach also incorporates the overall structural uncertainty from ranges of model responses from model ensembles, rather than the results based only on one model.Note that the MAGICC TCRE all values are slightly lower than that of CMIP5 mean value (as shown in Figures 3  (a) and 3(b)). is is why the estimated future climate from IPCC WG1 is slightly warmer than that from the IPCC WG3.
e results presented here are also subject to uncertainties due to different representation of carbon cycle processes in climate models.Also, since permafrost carbon cycle feedbacks and ice sheets are not currently represented in CMIP5 models considered here, they could lead to even higher warming levels and associated feedbacks.
INDC emission pledges are nonbinding and will be evaluated every five years, with the pre-evaluation by UNFCCC in 2018, and the first formal evaluation in 2023.e outcome will be used as the input for new INDCs.erefore, it is necessary to evaluate the climate impacts of the new INDCs in a timely manner.Our study shows that climate warming under current INDC scenarios is projected to greatly exceed the long-term Paris Agreement goal of stabilizing the global mean temperature at 2 °C or 1.5 °C above the preindustrial level, suggesting that rapid emission reductions in conjunction with negative emissions [30] may be necessary to ensure temperature stabilization consistent with the Paris Agreement.

Figure 1 :
Figure 1: National emissions under Intended Nationally Determined Contribution (INDC) scenarios for unconditional pledges in 2030.e emission target data of each country are calculated based on national INDCs in this study (see Section 2.1 and Supplementary TableS1).Emission level for each country is expressed as a percentage of total global emissions in 2030.e member states of European Union (EU) are shown as a whole as their emission target is submitted for the whole region.

5 Figure 2 :
Figure 2: Global INDC emissions compared with representative concentration pathway (RCP) scenarios.e black line shows the historical observed emissions.e blue lines show future RCP emission scenarios.Colored lines show Intended Nationally Determined Contribution (INDC) emissions under the unconditional pledge (red line) and conditioned pledge (pink line).Solid (red and pink) lines represent emissions before 2030 and dashed lines represent emissions after 2030.

Figure 3 :
Figure 3: Global temperature change as a function of cumulative carbon emissions from various lines of evidence.(a) Simulated GHG-attributable warming as a function of cumulative emissions based on representative concentration pathway (RCP) and 1%•yr −1 CO 2 increase (1PCT) simulations from the fifth phase of the Coupled Model Intercomparison Project (CMIP5) and Model for the Assessment of Greenhouse Gas Induced Climate Change (MAGICC) experiments.e MAGICC was used in order to compare the results of CMIP5 models used in IPCC Working Group 1 to the MAGICC model results in Working Group 3. (b) e ratio of GHG-attributable warming to cumulative carbon emissions (TCRE).e ranges of numbers in panel (b) indicate the amounts of cumulative emissions from the horizontal axis of (a) for each simulation.ese values indicate different ratios of TCRE at different amounts of cumulative emissions from various RCP simulations.

Figure 4 :
Figure 4: Global temperature increases under Intended Nationally Determined Contribution (INDC) scenarios.(a) unconditional INDC, (b) conditional INDC, and (c) temperature change relative to the preindustrial level (CMIP5: the fifth phase of the Coupled Model Intercomparison Project; MAGICC: Model for the Assessment of Greenhouse Gas Induced Climate Change).e method of calculating temperature response to INDC emissions is presented in Section 2.3.e INDC-induced warming is calculated based on (2) by multiplying ΔT INDC (cumulative carbon emission under the INDC scenarios) by TCRE all (the ratio of GHG-attributable warming to cumulative emissions).

Figure 5 :
Figure 5: Simulated model mean temperature changes in response to Intended Nationally Determined Contribution (INDC) emissions for unconditional pledges in (a) 2030 and (b) 2100.e temperature anomalies are relative to preindustrial levels.e simulation data from 12 CMIP5 models were used to produce an average state of warming pattern for INDC (see Section 2.3).e 2030 spatial pattern for INDC was estimated based on RCP 4.5 scenario experiment due to the similar emission levels of RCP 4.5 and INDC for 2030.e 2100 spatial pattern for INDC was estimated based on RCP 6.0 and 8.5 scenarios experiments due to the available simulations of the two scenarios for larger temperature increases.