Multimodel emission metrics for regional emissions of short lived climate forcers

Introduction Conclusions References Tables Figures


Introduction
Climate is impacted by various emitted gases and particles with a range of radiative efficiencies, lifetimes, and climate efficacies (e.g., Myhre et al., 2013).Emissions of CO 2 , N 2 O, and some of the gases included in the Kyoto Protocol are defined as long-lived greenhouse gases (LLGHGs).In addition, emissions of BC, OC, SO 2 , NH 3 , NO x , CO, and VOC cause changes in atmospheric levels of short lived climate forcers (SLCFs), such as ozone and aerosols (BC, OC, sulphate and nitrate).CH 4 is a well-mixed gas as the LLGHGs, but is often categorized together with the SLCFs since its lifetime is shorter than a realistic time scale for stabilizing anthropogenic influence on climate.There has recently been increased interest by policy makers to mitigate these SLCFs, for instance as advocated by the Climate and Clean Air Coalition (CCAC) motivated by co-benefits to climate and air quality (Schmale et al., 2014).The climate benefits for the next few decades of reducing SLCFs today are comparable to a climate policy on LL-GHGs (Smith and Mizrahi, 2013).However, Myhre et al. (2011) point out that reducing emissions of SLCFs today might potentially result in a delay in CO 2 mitigation, which may give unwanted long-term consequences (Pierrehumbert, 2014).Studies show that climate change in the long term is mainly governed by CO 2 emissions; however, mitigation of SLCFs may temporarily decrease the rate of warming (Shoemaker et al., 2013;Bowerman et al., 2013).Rogelj et al. (2014) argue that quantifying the climate impact of actual mitigating policies targeted on SLCFs is difficult, as the sources are common for a range of SLCFs and LLGHGs; thus, these linkages should be considered.
The impact of emissions of different SLCFs may be measured with the use of emission metrics which quantify an idealized climate impact per unit mass of emissions of a given species.Various applications exist (Fuglestvedt et al., 2003;Tanaka et al., 2010;Aamaas et al., 2013), the main ones are to (1) provide an "exchange rate" be-Figures

Back Close
Full tween different emitted species used in mitigation policies, (2) compare different activities and technologies that emit species over time such as in Life Cycle Assessment (LCA), and (3) compare in a simplified manner the climate responses of various emissions to gain and communicate scientific understanding.The most common emission metrics are time integrated radiative forcing (Absolute Global Warming Potential, AGWP) (IPCC, 1990) and temperature perturbation (Absolute Global Temperature change Potential, AGTP) (Shine et al., 2005(Shine et al., , 2007)), which normalized to CO 2 become GWP and GTP, respectively.Physically based metrics evaluate the idealized climate impact (integrated global mean RF for GWP or global mean temperature change for the GTP) over a certain time period (for the GWP) or at a given time after the emissions (for the GTP).This time period is called the time horizon and this choice is inevitably influenced by value judgments.Here we present metric values for a range of time horizons.Among the value choices are for instance looking at either temperature or forcing and what time horizon to pick (Fuglestvedt et al., 2003;Tol et al., 2012;Myhre et al., 2013).The Kyoto Protocol used GWP with a time horizon of 100 years.
Emissions metrics have normally been calculated for global emissions.However, due to short lifetimes compared to atmospheric mixing times, and because the chemistry and radiative effects on climate depends on the regional physical conditions, even the global mean radiative forcing depends on the region of emissions (Fuglestvedt et al., 1999;Wild et al., 2001;e.g., Berntsen et al., 2005;Naik et al., 2005).Then, the emission metric values will vary for different emission locations (Fuglestvedt et al., 2010).In addition, a distinct pattern in the temperature response appear from all forcings (Boer and Yu, 2003;Shindell et al., 2010).A growing literature investigates how the weights of the emission metrics change as emissions from different regions of the world are considered.Collins et al. (2013) assessed variations in emission metrics for four different regions (East Asia, Europe, North America, and South Asia) for aerosols and ozone precursors, based on radiative forcings from consistent multimodel experiments from the Hemispheric Transport of Air Pollution (HTAP) experiments given by Yu et al. (2013); Fry et al. (2012).Collins et al. (2010)

Conclusions References
Tables Figures

Back Close
Full ric values differ between regions, including vegetation responses.Bond et al. (2011) quantified differences in RFs for BC and OC emissions from different locations and types of emissions.For SLCFs, the impact of emissions depends also on the seasons.As the chemistry and radiative effects vary between summer and winter, the RF per unit emissions will differ between the seasons.An additional factor is that the magnitude of emissions fluctuates between the seasons, which can also be the case for LLGHGs.E.g., emissions of certain species from wood burning for domestic heating will be much larger in winter than summer (Streets et al., 2003).Bellouin et al. (2015) et al. (2015) to calculate emission metrics for the different regions and seasons.We produce emission metrics for standard pulse emissions, but also for an emission profile consisting of a ramp up period of 15 years followed by a sustained case, which can illustrate a gradual implementation of technology standards.As the study is based on several models running the same experiments, this data allows us to investigate the robustness in our findings.We analyze the robustness for individual species, as well as for hypothetical policy mitigation packages.Finally, we discuss how the emission metrics presented here can be used in mitigation policies.Introduction

Conclusions References
Tables Figures

Back Close
Full 2 Material and methods

Radiative forcing
An overview of the 4 different coupled-chemistry climate models or chemical-transport models presented by Bellouin et al. (2015), their resolution and species investigated (SO 2 , BC, OC, NH 3 , NO x , CO, VOC, and CH 4 ) is given in Table 1.Not all models have calculated RF for all species.While all four models give RFs for BC, OC, and SO 2 , only the OsloCTM2 calculated RF for NH 3 .Three models (OsloCTM2, HadGEM3, NorESM) have calculated RFs for the ozone precursors and CH 4 .
The calculations are based on different processes that affect RF, see Bellouin et al. (2015).For aerosols and aerosol precursors, all four models calculate the aerosol direct and 1st indirect (cloud-albedo) effect, except ECHAM6 which only diagnosed direct RF.For BC, OsloCTM2 estimated in addition the RF from BC deposition on the snow and semi-direct effect.Only a few previous studies, such as Bond et al. (2013), have included the semi-direct effect in emission metrics.For the ozone precursors and CH 4 , the total RF consists of the aerosol direct and 1st indirect effects, short-lived ozone effect, methane effect, and methane-induced ozone effect.Only OsloCTM2 includes nitrate aerosols, but nitrate aerosol RF has been used to complement the estimates by other models.
The best estimate of a species' RF is given as the sum of all the processes, in which the average across the models is used for each process.Not all processes, nor species, have been modeled by all models, and hence, the average for a process can be based on anything from only one model to four models.ECHAM6 is not included in the best estimate for OC and SO 2 , since this model does not diagnose the 1st indirect effect.As this 1st indirect effect is significant compared to the direct effect in OsloCTM2 and NorESM for OC and SO 2 , but to a much smaller degree for BC, we exclude ECHAM6 for the OC and SO 2 averages, but not for the BC average.
For the high and low estimate, we sum the highest and lowest value, respectively, for each individual process.Introduction

Conclusions References
Tables Figures

Back Close
Full

Emission metrics
In this study, we use the emission metrics GWP and GTP with varying time horizons.In all perturbations, RF is annually and globally averaged, thus, the responses are also annually averaged.AGWP for species i at time horizon H is defined as where RF is the time varying radiative forcing following a unit mass pulse emission at time zero.For aerosols, the radiative forcing values (RF ss ) (W m −2 (kg yr −1 ) −1 ) calculated by Bellouin et al. (2015) are based on assuming that the emissions are sustained for a year and hence the concentrations are close to equilibrium values because of their short lifetimes.These RF ss values have been converted into RF values (W m −2 kg −1 ) for an instantaneous emission for BC, OC, SO 2 , and NH 3 by the formula (Aamaas et al., 2013):

Conclusions References
Tables Figures

Back Close
Full where τ is the lifetime (yr) of the species.This conversion is only applicable when the adjustment time of the species is significantly less than one year.The AGTP is given as where IRF T (H − t) is the temperature response at time H to a unit radiative forcing at time t.These emissions metrics (AGWP, AGTP) are given in absolute forms.They can be normalized to the corresponding effect of CO 2 , where M is GWP or GTP, given as To calculate the time-varying RF for a pulse emission of CO 2 an impulse response function (IRF C ) for CO 2 is needed.Here we use the IRF C based on the Bern Carbon Cycle Model (Joos et al., 2013) as reported in Myhre et al. (2013).The IRF T is based on simulations with the Hadley Centre CM3 climate model (Boucher and Reddy, 2008).These parameterizations have uncertainties, and the spread due to IRF T is larger for SLCFs than for species with longer lifetimes (Olivié and Peters, 2013) The emission metric parameterizations for CH 4 comes from Myhre et al. (2013).
Emission metrics for pulse emissions are in principle the most useful metrics, even though emissions follow a given temporal profile.A pulse can be seen as an instantaneous emission, or constant emission during a short period ( H), followed by no emissions.In real life, implementing mitigation can be a gradual process where emissions are gradually reduced over some period, followed by a sustained level of emission reduction.This reflects regulations or technical improvements that are phased in over a given period and then sustained indefinitely.Such an emission profile, or mitigation profile, can be called a "ramp up".These different types of emission profiles are shown Introduction

Conclusions References
Tables Figures

Back Close
Full in Fig. 1.For ramp up or any other emissions scenarios, the emission metric can be calculated by a convolution.A temperature response is calculated as (5) E is the emission scenario and AGTP gives the temporal temperature perturbation for a unit of emissions.The absolute metrics for compound i for the ramp up scenarios where AM p i H is the corresponding absolute pulse metric (e.g., AGWP or AGTP) for time horizon H, and E(t e ) is the emission at time t e .The integral in Eq. ( 5) is the general notation, while we apply this in our calculations with the sum in Eq. ( 6).Note that the sustained case is a special case where E (t) = 1 for all t.For a ramp up period of mitigation of TH years, change in emissions in year t in the first TH years are t×E (t) TH and after that E (t).We show results only for a ramp up period of TH = 15 years, but we have also investigated other implementation rates.The total response for a scenario is found by multiplying Eq. ( 6) with the total emission change.Note that since emission metric values for SLCFs increase with decreasing time horizon (because they are short lived), their "ramp up" emission metrics values are significantly higher than the standard pulse based values.Emission metrics normalized to the corresponding absolute emission metric for ramp up emissions of CO 2 (M R i (H)) are calculated by Note that since the pulse metrics are given by region and season, so are the ramp up metrics (M R i (H)).For policymakers to apply this concept to compare different (n) sets of mitigation options (all following the same ramp up profiles over time, but with different mix of species, regions, and seasons) the net impacts (I n (H)) (i.e., AGWP or AGTP) for all options must be calculated according to Here ∆E i (j ) denotes the total mitigation (at the end of the ramp up period) of component i emitted in region j .

Best estimates
First, we present the best estimate of emission metric values for pulse emissions, see Table 2 for GTP(20) values.Additional values for GWP and for other selected time horizons are given in Table S1 in the Supplement.Due to space constraints, we can only present values for a few time horizons.What emission metric and time horizon to use depends on the application, and a range of different justified choices are possible (Aamaas et al., 2013).If the focus is on temperature change in the next few decades, GTP(20) can be applied.In Fig. 2, GTP(20) values are given for the different species, decomposed by a range of processes.Figure 3 presents results for GWP(100) for the ozone precursors.We first focus on a few selected time horizons, but Sect.3.1.5shows how emission metrics evolve for a range of time horizons.Introduction

Conclusions References
Tables Figures

Back Close
Full The uncertainties in Figs. 2 and 3 are given as the range across all contributing models.The uncertainty is in general larger than the variation between different regions and seasons.Thus, when including the uncertainty, it is less clear which region and season give the largest and smallest emission metric values.However, we will show in Sect.3.1.3and 3.1.4that the best estimate is more robust than the uncertainty bars indicate.
The emission metric values for the shipping sector are based on only two models (OsloCTM2 and NorESM).We do not provide uncertainty ranges for shipping due to the low numbers of models.Further, the robustness of these values presented is lower than for the other regions for the same reason.
We find distinct differences between regions and seasons for all species.For the aerosols BC, OC, and SO 2 , the magnitude of the total GTP(20) values are higher for emissions during summer than winter and larger for Europe than for East Asia.Collins et al. (2013) also estimated higher values for Europe than East Asia, while Fuglestvedt et al. (2010) based on earlier calculations in the literature gave partly conflicting results.
As a significant share of the emissions from the shipping sector, as well for global emissions, are occurring in the Northern Hemisphere, the seasonal variation is similar for these two categories except for BC for shipping.Seasonal variations are mainly driven by aerosol RF, which is mainly located in the shortwave spectrum.Greater sunlight duration in local summer yields stronger RFs (Bellouin et al., 2015).Seasonal differences in atmospheric lifetimes, caused by seasonality in precipitation, will also contribute.
For BC, the elevated aerosol direct effect in summer is partially cancelled out by a cooling effect by the semi-direct effect (see Fig. 2).The semi-direct effect is due to the absorption of solar radiation of particles, which affects the atmospheric static stability, and impacts on clouds.The impact of BC deposition on snow is largest for emissions during winter and larger for Europe than East Asia.The BC surface albedo effect is governed by the extent of snow and ice covered surface areas, but depends also on the availability of solar radiation where the BC is deposited.For Europe, the snow effect is 65 % of the direct effect in winter and 3.1 % in summer, while the corresponding Introduction

Conclusions References
Tables Figures

Back Close
Full percentages are 26 and 1.5 % for East Asia.The shares are similar for the shipping and global, with lowest shares for global emissions.As explained by Bellouin et al. (2015), this is due to atmospheric transport: according to the models, European emissions of BC are preferentially transported to the Arctic, where they modify the albedo of snow.Seasonality is driven by snow cover, which is larger in winter and early spring.
In Europe, the semi-direct effect is −32 % of the direct effect in summer and −3.3 % in winter, while it is −51 and 14 %, respectively, for East Asia.As the other regions are a mix of summer and winter because both hemispheres are included, the semidirect effect is smeared out on the two seasons, but largest in absolute value for NH summer.For NH 3 , the GTP(20) value is larger for Europe than East Asia, in summer but not for winter, as explained by Bellouin et al. (2015).Ammonium nitrate aerosol formation is strongly dependent on relative humidity and temperature, and competes for ammonium with sulphate aerosols, which has larger concentrations in local summer (Bellouin et al., 2011).Those complex interactions may explain different seasonalities in different regions, and will contribute to model diversity.
For the ozone precursors, the variability between regions and seasons is smallest for CO.For CO, GTP(20) values are higher for winter than summer.Due to the longer lifetime of CO during winter, a large fraction of the CO emitted during winter will undergo long-range transport and will be oxidized in relatively clean low-NO x environments.There CO-oxidation will reduce OH and thus increase the methane lifetime.As can be seen in Fig. 2, it is the indirect methane effect that leads to higher metric values for wintertime emissions.Furthermore, GTP(20) values of CO are slightly larger for East Asia than Europe.For VOC, the seasonal variability is opposite with highest GTP(20) values for summer.Further, GTP(20) values are higher for Europe than East Asia.The overall picture is a bit more complex for NO x .The seasonal difference is very small for GTP(20) values in East Asia.However, for Europe, the GTP(20) value is more negative for summer and less negative for winter.Shipping has the largest GTP(20) values in magnitude for all ozone precursors, driven by a large methane effect.Collins et al. (2013) observed the same annual pattern for Europe and East Asia as we do.Introduction

Conclusions References
Tables Figures

Back Close
Full For aerosols emissions and the major aerosols precursors, the relative ratios between the different regions and seasons are constant while varying the emission metric and time horizon applied.On the other hand, the relative ratios between different emission metric values for the ozone precursors differ with varying emission metrics and time horizons.The ratios for the aerosols are fixed since the aerosols have little effect on perturbations of atmospheric composition and components with long adjustment times.By contrast, the ozone precursors affect processes with longer time constants.By causing a change in OH-levels, methane with an adjustment time of about 10 years is perturbed.Hence, we also show GWP(100) values for the ozone precursors (Fig. 3), while similar figures for the other species are provided in the Supplement (Fig. S1).For the ozone precursors, the aerosol direct and indirect effect and the short-lived ozone effect are given relatively more weight for GWP(100) than GTP(20) than the methane effect and methane-induced ozone effect, since GWP integrates the RF up to the time horizon, while GTP is an end-point indicator.As the time horizon increases, the relative contribution from methane and methane-induced ozone increases and the contribution from aerosols and short-lived ozone decreases.The overall picture presented here for GTP(20) and GWP(100) is mostly similar.But for NO x , no significant seasonal difference was observed in GTP(20) values for East Asia, the value in winter is almost twice as negative as the summer values for GWP(100).

Comparison with literature
As already observed, the variations with respect to regional emissions for emission metric values are in line with Collins et al. (2013).Fuglestvedt et al. (2010) also presented emission metrics with respect to regional emissions based on earlier calculations in the literature, but with some conflicting results between available studies.Our findings are therefore somewhat in line with Fuglestvedt et al. (2010).In general, the specific emission metric values are also comparable with Collins et al. (2013).However, a complete comparison is not possible as we have included the effect of aerosols for the ozone precursors and the semi-direct and deposition on snow effect for BC.Introduction

Conclusions References
Tables Figures

Back Close
Full The findings are also generally similar to previous estimates for emission metrics of global emissions (e.g., Fuglestvedt et al., 2010), with some discrepancies we will discuss here.A comparison of modeled GWP and GTP values with a selection from the literature for some selected time horizons is given in Table S1.
For BC, Bond et al. (2011Bond et al. ( , 2013) ) presented about 50-60 % higher emission metric values, while other studies (Fuglestvedt et al., 2010;Collins et al., 2013) are in line with this study and Hodnebrog et al. (2014) give significant lower values.As discussed in Hodnebrog et al. (2014), the atmospheric lifetime of BC may be shorter and the BC emissions may be larger than previously thought (e.g., Fuglestvedt et al., 2010) leading to emission metric values almost halved compared to previous estimates (−44 % for the example given in Hodnebrog et al., 2014).The OC values in our study are almost 200 % higher in magnitude than the literature, driven by the high values in one of the models (NorESM).The OC values from NorESM are driven by a strong indirect effect.When this indirect effect is excluded, the NorESM value is similar to the others as well as the literature.For SO 2 , the emission metric values for the winter season are similar to the literature, while more than doubled for summer.As for OC, the more negative emission metric values for SO 2 are driven by the inclusion of the indirect effect.The one study (Shindell et al., 2009) we found on NH 3 gave emission metric values that are about the double of our annual average.The literature shows a wide range in the emission metric values for NO x depending on the source and region.Our estimates are within this range but, on the more negative side within the range, about 80 % stronger than the values used for land-based emissions in Myhre et al. (2013).The emission metric values for CO are roughly 0-30 % higher than in the literature, partly driven by the additional positive impact of including the aerosol effect.For VOC, the emission metric values are roughly the double or more than those found in the literature, even with a negative contribution from the aerosol effect (Bellouin et al., 2015).

Robustness for individual species
The differences in the emission metric values between the emission regions and seasons of emissions, seen for the best estimate holds generally in each model, which strengthens our confidence in the modeled variations between regions and seasons.
For emissions of aerosols and their precursors, the magnitude of GTP(20) values is higher in summer than winter in 88 % of the model cases and another 8 % are marginally the other way.The consistency between the individual models and our best estimate based on the models is 100 % for SO 2 .The metric values for European emissions are larger in magnitude for most cases than East Asia.In summer, this is true for 92 % of the cases and 58 % in winter in addition to 17 % that are marginally the opposite.Yu et al. (2013) also observed that the regional dependency in RF was robust for a number of models with the same regional pattern as in our study.
For the ozone precursors, the variation in GTP(20) values observed for the best estimate also holds for most of the models.For both regional and seasonal variability, 83 % of the model cases agree with the best estimate.For CO, all cases agree that the GTP(20) values are larger for East Asian emissions than European emissions and for winter than summer, even though the relative differences in GTP(20) values between Europe and East Asia in summer and winter are relatively small.The findings for NO x and VOC are also relatively robust, where the model cases agree 83 % for NO x and 67 % for VOC.The same tendencies in the regional pattern were also found by Collins et al. (2013).metrics, because there might be significant correlations between species.By robustness here, we mean how uncertain is the total climate impact of a given set of emission changes (changes of multiple species, seasons and regions) and related to this how robust would a ranking (in terms of net climate impact) of possible mitigation measures be, given the individual uncertainties shown in Figs. 2 and 3.

Robustness in total climate impact
Models with more efficient vertical transport and/or slow removal of aerosols by wet scavenging will tend to give longer lifetimes for the aerosols and thus stronger RF per unit emission for all aerosol species, and thus emission metric values for the individual species and seasons would be correlated.This means that the ranking of measures and the net impact of measures that lead to reduction in emissions of coemitted species that cause a cooling effect could be more robust.Similar effects can be expected across ozone precursors due to non-linear chemistry effects and removal efficiencies; for instance, such correlations across models were observed for the climate effect of NO x emissions from aviation by Holmes et al. (2011).To investigate this we first focus on the correlation.To put all species on a common scale we calculate the normalized variability (across species, regions and seasons) for the best estimate (NV BE ) and for the individual model estimates (NV m ) NV BE (r, s, i ) = M BE (r, s, i ) − M BE,min (i ) and M BE (r, s, i ) denotes the best estimate for the emission metric value for species i , region r and season s, while M m (r, s, i ) denotes the emission metric value from a single model m for species i , region r and season s.
The values of NV BE are numbers between 0 and 1. Figure 4  of the symbol indicate component.Since the processes that could lead to correlations are somewhat different for aerosols and ozone precursors (e.g.non-linear chemistry effects for the latter) the species are split into two separate panels.Figure 4 clearly shows the correlation between the species for the individual model emission metrics.For the aerosols, HadGEM and particularly NorESM tend to give higher (in absolute terms, i.e. more negative for cooling agents) emission metric values compared to the best estimate, while ECHAM gives much lower values.For the ozone precursors, the picture is the opposite, with NorESM being lower than the BE while the OsloCTM is higher.This indicates that for both aerosols and ozone precursors there are generic features in the models related to representation of key processes (e.g.vertical mixing, wet scavenging, ozone production efficiency etc.) that systematically affects the emission metric values.
These correlations between the estimates for the individual species have to be taken into account when the uncertainty in the net effect of a multi-component mitigation policy is estimated.Since different SLCFs are often co-emitted, most mitigation options will affect emissions of several species at the same time.The uncertainty in the estimate of the net effect depends on the composition of the mitigation, i.e. mix of species, regions, and sectors.To be useful for policymaking, the emission metrics should be robust enough so that there is trust in the sign of the net effect of a mitigation measure and that the uncertainty in the emission metrics does not hinder a ranking of different measures when cost-efficiency is considered.Figure 5 shows the estimates of the net effect (here in terms of temperature change after 20 years, i.e. using AGTP(20) for pulse emissions) when using the sets of emission metrics from the individual models.First, we consider a global mitigation of a 10 % reduction in emissions of all SLCFs for which the best estimate is positive for the AGTP(20) (BC, OC, and VOC -labelled B on Fig. 5), and then a 10 % global reduction of all SLCFs (an extreme case of also reducing the co-emitted cooling species OC, SO 2 , and NO x -case A on Fig. 5).The shipping sector is not included in this sensitivity test as the best estimate is only based on two models.ECHAM6 did not calculate RFs for the ozone precursors, therefore, values for Introduction

Conclusions References
Tables Figures

Back Close
Full the best estimate is given for those species.NH 3 is not included, as only OsloCTM2 provided RF estimates of that.These scenario estimates are based on emission inventories for 2008 (Klimont et al., 2015).For a 10 % reduction in emissions of the warming SLCFs (BC, CO, and VOC), the best estimate gives a global reduction in temperature of 0.55 mK 20 years after a pulse, with a spread of −0.39 to −0.83 mK.When the cooling components are included, the best estimate gives a global warming of 0.48 mK, with models ranging from 0.05 to 0.46 mK.Hence, all models agree that a reduction of those six SLCFs will cause warming, but for two of the models only a marginal warming.The black bars in Fig. 5 give the uncertainty in the net global temperature effect assuming all the metric values are independent.This gives a similar or narrower uncertainty interval than the spread of the estimates using the individual model metrics, again showing that there is considerable correlation between in the model estimates.However, if the difference between the models were 100 % systematic (i.e. one model always giving the lowest estimates by magnitude and another model giving the highest), then the model based interval would be given by the blue bar in Fig. 5. From this analysis, we conclude that the uncertainty for an estimate of the net temperature effect of multi-component emission change is enhanced due to the correlations; however, for mitigation measures that mainly change emissions of species with positive GTPs, the sign of the global temperature signal is robust.
Since not all processes are included in all the models, the average of all models in Fig. 5 will differ from the best estimate.This deviation is observed in both scenarios, but clearest for a mitigation scenario including both warming and cooling SLCFs, as the net climate impact is a sum of large positive and negative numbers.The processes not included are dominated by cooling.Three out of four models do not include the cooling from the semi-direct effect of BC, as well as the mainly cooling from nitrate for the ozone precursors and SO 2 .As a consequence, the individual models tend towards more cooling or less warming than the best estimate for a mitigation scenario of SLCFs.
Our findings show that the robustness is largest for individual species, i. a subgroup of species that correlates, such as aerosols.Lowest robustness is given for mitigation for all SLCFs.However, we observe that all models agree whether two hypothetical mitigation scenarios give warming or cooling.

Variations with time horizon
We have until now presented emission metric values at certain fixed time horizons; however, these values vary greatly with time horizon.SLCFs impact the atmosphere for a short time, as aerosols and aerosol precursors have atmospheric lifetimes of about a week.Methane, however, has an atmospheric perturbation lifetime of about 12.4 years (IPCC, 2013).Due to the inertia in the climate system, the climate is impacted for at least 10-20 years from a radiative forcing lasting only a week (Peters et al., 2011;Solomon et al., 2010;Fuglestvedt et al., 2010).The denominator in the emission metrics is CO 2 , which impact the atmosphere for centuries (IPCC, 2013).However, aerosols are very strong at perturbing the radiative balance of the Earth while they are situated in the atmosphere; for instance, the radiative efficiency (W m −2 kg −1 ) of black carbon is about a million times larger than the radiative efficiency of CO 2 .Thus, the magnitude of the normalized emission metric values is very high for short time horizons, but decreases rapidly with increasing time horizon.The aerosols have the highest emission metric values in magnitude for the shortest time horizons, see Figs. 6 and 7 for GTP and GWP values in the first 50 years after a pulse emission.Additional figures are provided in the Supplement.NO x often has a positive emission metric value for the first 5-10 years, followed by negative numbers, as the sum of the short-lived effects are positive and the longer lived effect negative.However (see Fig. 7), we find cooling already from year one for emissions in Europe during all seasons and East Asia during winter as the cooling from the aerosol effect is as large as or larger than the short-lived ozone effect.This aerosol effect is cooling for all regions except shipping.
The time dimension is especially important for NO the first 10 years, but becoming the second most negative after 20 years.For a specific region and season, the weighting between the aerosols and ozone precursors is also changing with variable time horizon.

Global temperature response
We have applied the emission metrics on an emission dataset for year 2008(Klimont et al., 2015)).The variability discussed in the previous section is also found in the global temperature response for regional and seasonal emissions (Fig. 8).A seasonal profile is included in the emissions, with typically largest emissions in the winter season.As for the emission metric GTP, the temperature response drops rapidly off due to the short lifetimes of the SLCFs.In general, the total temperature response is governed by the SO 2 emissions.Hence, the total climate impact is a cooling for all regions and seasons, but largest for emissions in summer.The emission mix is different between the regions.For instance, SO 2 and NO x generally dominate for shipping.Europe and East Asia have a broader mix of SLCFs that impact the climate.The temperature perturbation, dominated by cooling, is in agreement with Aamaas et al. (2013), who also showed that the warming of CO 2 is larger than the net cooling from the SLCFs after only 15 years.We have presented the global temperature response, while regional variations will occur beyond this global mean response (e.g., Lund et al., 2011).

Gradual implementation of mitigation
We have calculated emission metrics for pulse emissions, which is the typical way of presentation.However, changes in emissions are often gradual in real life.In this section, we present how the emission metric values differ based on a gradual implementation of mitigation policy (see Fig. 9), which is calculated by convolution as given in Eq. ( 6).We show results only for a ramp up period of 15 years, but we have also looked at other implementation rates.The emission metric values presented here are Introduction

Conclusions References
Tables Figures

Back Close
Full for Europe in the summer season.We complement with emissions of CH 4 , based on parameterizations in Myhre et al. (2013).
For species that have a shorter influence on the climate system than CO 2 , the normalized emission metric values will almost always be larger in magnitude for sustained emissions than pulse emissions.The only exception is for species with competing processes on different timescales, such as for NO x in Fig. 9.The growth scenarios give slightly higher normalized emission metrics than the sustained case, but those emission metric values approach each other in the long term.The longer the ramp up period lasts, the larger the emission metric value becomes, but the value converges to the sustained emission case for time horizons beyond the ramp up period.The normalized emission metric values are higher in the growth scenarios than the sustained case since the impact of the shorter lived effects are given more weight than CO 2 .Hence, a mitigation scenario that will have a gradually increasing effect over several years will for most species have a higher metric value than for mitigation that instantly takes effect.What this means is that one obtains the benefit of mitigating SLCFs (i.e., higher CO 2 equivalent emission reduction and thus higher value on an emission trading market or in a cost-effectiveness analysis) now.The reason is the planned emission reductions of the shorter lived species close to the time horizon has a large impact.Hence, these emission metrics for ramp up scenarios should be used with care.If there is a chance that the emission reductions are reversable and will not be kept in place (or replaced by even stronger reductions) until the time horizon, the ramp up metrics will overestimate the effects.
We also present the temporal evolution for all the regions and seasons for BC and NO x in Fig. 10, while Fig. 9 only showed these emission metric values for Europe in summer.While the regions and seasons are ranked the same for all time horizons for the aerosols, the ranking may differ for the ozone precursors for different emission metrics and different time horizons due to competing processes on different timescales, especially for NO x .Introduction

Conclusions References
Tables Figures

Back Close
Full The other significant difference between emission metrics based on pulse and ramp up emissions is the sign switch for NO x (see Fig. 9).In the pulse case, the GTP values are negative or turn negative within the first 6 years, with NO x from shipping during summer taking 10 years.The sign switch is much slower for the ramp up scenario emission metrics.Even after 10 years, 5 out of the 8 cases give positive GTP values (see Fig. 10).In the long run (> 22 years), all the GTP values with the exception of shipping in summer, are negative.Thus, if a time horizon of 10 years is picked, the mitigation policies of NO x will depend highly on the assumed emission scenario.

Policy implications
Emission metrics can be applied as an "exchange rate" between different emissions in climate polices, such as for different LLGHGs in the Kyoto Protocol.While the calculations of how emissions impact the climate build on scientific knowledge, how to use the emission metrics is given by political choices.There is no particular reason why there should be one and only one goal for our climate policy (Fuglestvedt et al., 2000;Rypdal et al., 2005;Daniel et al., 2012;Sarofim, 2012;Jackson, 2009;Victor and Kennel, 2014).In particular there may be harmful impacts of exceeding a longterm temperature constraint (e.g., 2 • C), while at the same time there is more immediate concern about short term effects over the next decade or so.The rationale behind a climate policy focusing on SLCFs must be that there are potential harmful effects of climate change over the next few decades.However, CO 2 and other LLGHGs should also be included in evaluation of possible mitigation measures under a short-term goal.Historically, emission metrics within international climate policy have been applied to emissions of LLGHGs.However, as the uncertainty for the emission metrics of SLCFs is reduced and the values become more robust, this opens up for regimes that also include non-methane SLCFs beyond CH 4 , e.g., CCAC.Recently, Mexico included BC in their Intended Nationally Determined Contribution (INDC) (Mexico, 2015).
A general difference between LLGHGs and SLCFs is that the location of the LLGHG emissions does not matter, while we have shown that different locations, as well as tim-Introduction

Conclusions References
Tables Figures

Back Close
Full ing of emissions, will cause different impacts of SLCFs (Fuglestvedt et al., 1999;Naik et al., 2005;Berntsen et al., 2005Berntsen et al., , 2006;;Shindell and Faluvegi, 2009).In addition to differences in the total global response, the spatial distribution of the impact depends on the location and timing of the SLCFs emissions.Further, we have shown that individual models may give significantly different emission metric values than other models.

Conclusions
We have presented emission metrics for regional emissions of several SLCFs (BC, OC, SO 2 , NH 3 , NO x , CO, and VOC) based on four different models.We have focused on the emission regions Europe and East Asia, but also given numbers for the global shipping sector and total emissions from all countries.Values have been estimated for emissions in both the summer and winter seasons.For the aerosols, the magnitude of the emission metric values is larger for Europe than East Asia and for summer than winter.The variability between the models is generally larger than the variations between regions and seasons.However, most models agree that specific regions and seasons have larger emission metric values than others.Hence, the robustness of this ranking is better than can be interpreted from the variability between models.The co-variability between models is also seen for the ozone precursors.For CO, the emission metric values are larger for East Asia than Europe and for winter than summer.The pattern is the opposite for VOC with larger emission metric values in Europe and in summer.NO x is more complex with more negative values in summer than winter for Europe.
In East Asia, we model no significant difference between the seasons for GTP(20) for NO x , while the GWP(100) for winter emissions is more negative.
We have also calculated emission metrics for transient scenarios where we consider a ramp up of the emission over time.This emission metric will better represent the effect of imposing a mitigation measure (i.e. a new technology standard) that is known to give a long term change of emissions.For species that have a shorter influence on the atmosphere than CO 2 , the magnitude of the emission metric value is larger Figures for a mitigation scenario with a gradually increasing effect over several years than for the standard pulse based emission metric.The only exception is species that have competing short and longer lived effects that are positive and negative, notably for NO x .We observe variability in the emission metrics between different regions and seasons, however, with varying robustness between the models.As the certainties in the numbers increases, the regional and seasonal differences may be accounted for in mitigation policies, agreements and potential trading schemes involving SLCFs.One robust finding in our study is that, per unit mass of emissions, emissions of aerosols and their precursors in Europe should likely be given more weight than emissions in East Asia, as well as emissions in summer likely more weight than in winter.When emission metrics are applied, the selection of the specific emission metric and time horizon is of significance.The emission metric values for SLCFs drop quickly with time horizon.For the ozone precursors, the ranking between different regions and seasons can vary with different time horizon.Thus, emission metrics must be used based on careful consideration of these factors.
The Supplement related to this article is available online at doi:10.5194/acpd-15-26089-2015-supplement.Introduction

Conclusions References
Tables Figures

Back Close
Full Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | The emission metric values for CH 4 are mostly lower than those in Myhre et al. (2013) (29 and 19 % lower for GTP(20) and GWP(100), respectively), mainly due to a shorter methane atmospheric lifetime, as well as a smaller contribution from the indirect effect on ozone.Discussion Paper | Discussion Paper | Discussion Paper | Emission metrics are used to quantify the climate impacts of different sets of emission changes following either mitigation policies or changes caused by some other mechanisms (e.g.technological development).However, the uncertainties given by the model ranges for individual regions, seasons and species shown in Figs. 2 and 3 do not provide a good indication for the robustness of the total impacts estimated by the emission Discussion Paper | Discussion Paper | Discussion Paper | is a scatter plot between NV BE and all the individual NV m values, where the colors indicate model and shapes Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | e., what region and season of emissions to mitigate for an individual species.Next follows Discussion Paper | Discussion Paper | Discussion Paper | x and the other ozone precursors, as different regions and seasons are given different weights with different time horizons.For instance, shipping in summer has most positive GTP values for NO x of all cases in Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper |

Figure 10 .
Figure 10.The emission metric values for ramp up scenario emissions.GTP (top) and GWP (bottom) are given for BC (left) and NO x (right).
These global-mean RFs of various species were calculated for emissions in different regions.The three regions, following tier 1 HTAP regions, are Europe (Western and Eastern Europe up to 66• N including Turkey), East Asia (China, Korea, and Japan),