As researchers who have published over recent years on the issue of comparing the climate effects of different greenhouse gases, we would like to highlight a simple innovation that would enhance the transparency of stocktakes of progress towards achieving any multi-decade-timescale global temperature goal. In addition to specifying targets for total CO2-equivalent emissions of all greenhouse gases, governments and corporations could also indicate the separate contribution to these totals from greenhouse gases with lifetimes around 100 years or longer, notably CO2 and nitrous oxide, and the contribution from Short-Lived Climate Forcers (SLCFs), notably methane and some hydrofluorocarbons. This separate indication would support an objective assessment of the implications of aggregated emission targets for global temperature, in alignment with the UNFCCC Parties’ Decision (4/CMA.1)1 to provide “information necessary for clarity, transparency and understanding” in nationally determined contributions (NDCs) and long-term low-emission development strategies (LT-LEDSs).

While differences remain between us regarding how best to set fair yet ambitious targets for individual emitters2,3,4,5, including how any additional information might be used, and the interpretation of the Paris Agreement, it is important to emphasise the high level of agreement on the underlying science of how different greenhouse gases affect global temperature. The 2018 IPCC Special Report on 1.5 °C (SR1.5)6 stated “Reaching and sustaining net-zero global anthropogenic CO2 emissions and declining net non-CO2 radiative forcing (Planetary energy imbalance resulting directly from human-induced changes.) would halt anthropogenic global warming on multi-decadal timescales (high confidence). The maximum temperature reached is then determined by cumulative net global anthropogenic CO2 emissions up to the time of net zero CO2 emissions (high confidence) and the level of non-CO2 radiative forcing in the decades prior to the time that maximum temperatures are reached (medium confidence)”. The IPCC 6th Assessment Report (AR6)7 confirmed “limiting human-induced global warming to a specific level requires limiting cumulative CO2 emissions, reaching at least net zero CO2 emissions, along with strong reductions in other greenhouse gas emissions”.

Parties to the Paris Agreement agreed in Katowice in 2018 (Decision 18/CMA.1)1 to report past emissions of individual gases separately and use 100-year Global Warming Potentials (GWP100) when aggregating them to CO2-equivalent (we refer to these here as CO2-e100 emissions). The separate specification of individual gases minimises ambiguity in determining the climate impact of past emissions. NDCs and other future targets are, however, almost always expressed in terms of aggregate CO2-e100 emissions only, for which the implications for global temperature are ambiguous8,9. Separate specification of the contribution from CO2 helps, but ambiguity in global temperature outcomes remains if targets for non-CO2 gases comprise a mixture of long-lived climate forcers (LLCFs), such as nitrous oxide, with atmospheric lifetimes around 100 years or longer, and SLCFs, such as methane, most of which have lifetimes shorter than 20 years10.

Specifying the contributions of all gases individually in future targets as well as the reporting of past emissions would resolve the ambiguity in global temperature outcomes, and would also help quantify non-climate benefits of emission reductions, especially for methane11. Governments and particularly corporations may, however, wish to retain some level of aggregation across gases to allow flexibility in how they achieve their targets. Fortunately, a much less restrictive approach delivers almost all the transparency benefits from a climate perspective. The climate system responds similarly over a broad range of timescales to equal emissions expressed in tonnes of CO2-e100 of all LLCFs, including CO212. Likewise, the net radiative forcing due to SLCFs on multi-decadal timescales is similar to the aggregated rate of SLCF emissions expressed in tonnes of CO2-e100 per year multiplied by the 100-year Absolute Global Warming Potential (AGWP100) of CO213. With this additional information, it is straightforward to express the SR1.5 statement quoted above in terms of CO2-e100 emissions: human-induced warming over any multi-decade time-interval is approximately the sum of (i) aggregate CO2-e100 emissions of LLCFs, including CO2, multiplied by a constant parameter, the Transient Climate Response to cumulative CO2 Emissions, or TCRE14 (the TCRE can alternatively be thought of as the Absolute Global Temperature-Change Potential for a sustained emission of CO2 divided by the time-horizon, AGTPs/H13); (ii) any change in decadal-average radiative forcing due to SLCFs multiplied by another constant parameter, the Transient Climate Response to Forcing, or TCRF, another name for the “fast” component(s) of the climate response15; and (iii) a gradual adjustment to average SLCF forcing16, all evaluated over the same time-interval.

Hence a separate indication of the contributions of LLCFs and SLCFs in emission targets, or equivalently the LLCF contribution to total CO2-e100 emissions, is required to allow for the global temperature outcome to be calculated relatively unambiguously. It is important to note, however, that the evaluation of emission targets at the national or corporate level cannot be undertaken from a physical science perspective alone, but also depends on economic, social, equity and political considerations2,3,4,5,17, including responsibility for past warming, capacity for and costs of abatement, and non-climate impacts. Separate specification would also facilitate the use of alternate or flexible emission metrics, which may be useful for achieving a cost-effective emission trajectory over time18 or addressing specific policy goals such as limiting near-term rates of warming19. Indicative contributions from LLCF and SLCF abatement would not preclude trade-offs between them, but would clarify the need to monitor the temperature impacts of any such trade-offs over a range of timescales20.

It has long been accepted21 that stringent mitigation of both LLCFs and SLCFs is needed to meet any ambitious temperature goal, but making progress on two fronts necessitates monitoring progress on two fronts. Some countries (but very few companies) already specify the contribution of LLCFs and/or SLCFs to total CO2-e100 emissions in NDCs, LT-LEDSs and science-based targets (https://sciencebasedtargets.org/) communicated under the Greenhouse Gas Protocol. Quantifying the aggregated implications of these targets for future global temperature simply requires a much wider uptake of this practice, representing a simple and achievable innovation that would enhance the transparency of any stocktake of progress towards any global temperature outcome. Separate indication of LLCF and/or SLCF contributions could be communicated by countries as additional information consistent with Decision 4/CMA.1. This does not have to affect any existing or planned NDCs or long-term net zero strategies22 communicated using aggregate CO2-e100.

Why separate specification is so useful

To quantify the SR1.5 and AR6 statements quoted above, human-induced global temperature change over a multi-decade time-interval ∆t, relative to the level of human-induced warming at the beginning of that interval (e.g. the present day or pre-industrial), can be decomposed using the framework articulated above as follows:

$${\Delta} T = \kappa _E\overline {E_C} {\Delta} t + \kappa _F\left( {{\Delta} F_N + \rho \overline {F_N} {\Delta} t} \right),$$
(1)

where \(\overline {E_C}\) and \(\overline {F_N}\) are globally aggregated average CO2 emission-rates and non-CO2 radiative forcing, respectively (so \(\overline {E_C} {\Delta} t\) is cumulative CO2 emissions), and ∆FN is the change in decadal-average non-CO2 forcing, all evaluated over that interval (the geophysical “Zero Emissions Commitment” is expected to be relatively small over a multi-decade time-interval23, but this may not be the case on longer timescales). The coefficients κE (the TCRE) and κF (the TCRF, or “fast” component of the climate response to any forcing change, denoted c1 in ref. 12, or sum of fast components24: see supplementary material), are both scenario-independent in the absence of strongly non-linear carbon cycle feedbacks or climate response. The only scenario-dependent coefficient is ρ, the fractional Rate of Adjustment to Constant Forcing (RACF), or the relatively small fractional rate at which forcing needs to decline to maintain stable temperatures. It depends on how fast and how recently FN has increased (this term represents the delayed adjustment to past forcing increases, so is larger for more recent and rapid increases). If FN varies only on multi-decadal timescales, ρ = c2/(κFs2), where c2 is the “slow” (multi-century) component of the climate sensitivity, and s2 the deep ocean thermal adjustment timescale. For representative12 coefficient values, ρ 0.3% per year, making this third term usually small.

Aggregate CO2-e100 emissions cannot be used to calculate FN if these comprise a mixture of LLCFs and SLCFs. Aggregate CO2-e100 emissions of LLCFs, EL, can, however, be combined unambiguously and have the same impact on global temperature on decade to century timescales as the corresponding quantity of CO2. Likewise, aggregate CO2-e100 emissions of SLCFs, ES, multiplied by the AGWP100 of CO2, A100, give SLCF radiative forcing, FS (A100 normally includes a first-order estimate of the impact of carbon cycle feedbacks25 so, for consistency, this should also be included in the GWP100 values used to compute ES).

For emissions reported as CO2-e100 the above expression can therefore be re-written (now grouping all LLCFs with CO2):

$${\Delta} T = \kappa _E\overline {E_L} {\Delta} t + \kappa _F\left( {{\Delta} F_S + \rho \overline {F_S} {\Delta} t} \right),$$
(2)

or equivalently, using FS=A100ES on multi-decadal timescales,

$${\Delta} T = \kappa _E\overline {E_L} {\Delta} t + \kappa _FA_{100}\left( {{\Delta} E_S + \rho \overline {E_S} {\Delta} t} \right).$$
(3)

Hence ∆T can be estimated directly using well-known (albeit uncertain) climate system properties if, and only if, total CO2-e100 emissions of long-lived climate forcers, EL, are specified in emission targets together with total CO2-e100 emissions, EL+ES; or, equivalently, EL and ES are specified separately. ∆T cannot be calculated from the sum of EL+ ES alone.

This is illustrated by Fig. 1, which shows the impact of LLCF and SLCF emissions, expressed as CO2-e100, on global temperature change over a multi-decade period, relative to the level of warming at the beginning of that period, calculated with a simple climate model12. Stylised cases of constant (darker shades) and step-change (+10%, lighter shades, and −50%, dotted lines) emissions are shown in panels a and c. Warming due to LLCF emissions (the term \(\kappa _E\overline {E_L} {\Delta} t\) in Eq. (3)) increases linearly with cumulative emissions in all three cases (panel b). Warming due to an ongoing constant emission of an SLCF that started decades before the beginning of this period (the \(\kappa _FA_{100}\rho \overline {E_S} {\Delta} t\) term) also increases linearly (panel d, darker blue) but at a slower rate per tCO2-e100 emitted (by a factor of about 4, because κE ≈ 4 × κFA100ρ): global temperatures have already partially equilibrated with this constant emission (by how much depends on how long ago these SLCF emissions began, which is why ρ is the only scenario-dependent coefficient in these expressions). Finally, warming due to an increase in SLCF emissions (the κFA100ES term, panel d, lighter blue) is 4–5 times greater than would be expected from the same increase in tCO2-e100 emissions of an LLCF (panel b, lighter red) over the 20 years following the increase (κFA100 ≈ 4.5 × κE × 20 years). Hence the AR6 statement “expressing methane emissions as CO2 equivalent emissions using GWP100 overstates the effect of constant methane emissions on global surface temperature by a factor of 3–4 … while understating the effect of any new methane emission source by a factor of 4–5 over the 20 years following the introduction of the new source”26 applies to the impact of global emissions of any SLCF. Any decrease in SLCF emissions also has a much greater impact on temperatures over a multi-decade period per tCO2-e100 avoided than a corresponding decrease in LLCF emissions (red and blue dotted lines) (Fig. 1).

Fig. 1: Stylised LLCF and SLCF emissions and resulting global temperature change ∆T over a multi-decade period.
figure 1

Darker bands in panels a and c show, respectively, constant LLCF and SLCF emissions of 1 tCO2-e100 per year starting some decades before the interval shown. Pale bands show a 10% increase one-quarter of the way through the interval shown, while dotted lines show a 50% decrease. Resulting temperature changes relative to the start of this interval shown in panels b and d, calculated using a simple climate model: vertical axes in b and d are scaled identically to illustrate smaller rate of warming due to constant SLCF emissions and much larger warming impact of any change in SLCF emissions relative to the warming due to identical CO2-e100 LLCF emissions. Vertical arrows in the right show predicted contributions to ∆T from the individual terms in Eq. (3): three arrows in panel b show cumulative LLCF emissions over this interval multiplied by the TCRE for the three scenarios shown; the lower and upper arrows in panel d show, respectively, the predicted warming due to ongoing constant SLCF emissions and additional warming due to the 10% increase. The figure illustrates that Eq. (3) allows reliable, if approximate, prediction of multi-decade warming ∆T if, and only if, LLCF and SLCF emissions are specified separately.

Temperature changes in the figure are calculated using a particular model, LLCF, SLCF and scenario. The figure would, however, appear similar if another model, combination of gases or scenario of prior emissions were used, provided emissions do not change rapidly immediately before the beginning or end of the period shown, because the relationship between emissions and warming expressed in Eq. (3) is generic. Individual terms in Eq. (3), assuming constant coefficients, are shown by the arrows on the right of panels b and d. These match the warming calculated by the explicit simple climate model within modelling uncertainties. The figure shows temperature change relative to the start of the period rather than absolute warming because the latter is not determined by Eq. (3) but depends on the prior LLCF and SLCF emissions history (the specific scenario used to generate this figure is shown in full in the Supplementary Information).

Temperature change ∆T over a multi-decade period depends, to first order, only on cumulative emissions of LLCFs \(\overline {E_L} {\Delta} t\), cumulative emissions of SLCFs \(\overline {E_S} {\Delta} t\), and net change in total SLCF emission rates ∆ES, over that period alone. As the SR1.5 and AR6 emphasised, future warming depends on future emissions. Making use of this information, however, requires both EL and ES to be specified: only specifying the sum EL + ES introduces an ambiguity in temperature outcome.

Separate specification also facilitates assessing the implications of different metrics. For example, aggregate CO2-equivalent emissions using the 20-year Global Warming Potential (GWP20) can be approximated by EL+ 3ES if both EL and ES are reported as CO2-e100, with a slightly higher multiplicative factor (up to 4) if ES is dominated by forcers with lifetimes of order one year (Table 8.A.1 of ref. 12 shows that GWP20 values are similar to GWP100 values for LLCFs and 3 or 4 times GWP100 values for gases with lifetimes of order a decade or a year, respectively). Finally, we re-emphasise that these expressions capture our physical understanding of how global emissions of LLCFs and SLCFs collectively determine global temperature change, and illustrate the utility of separate specification of EL and ES. How this understanding is used to inform the assessment of the adequacy of individual emission targets depends on other considerations listed above and cannot be argued from a physical science perspective alone. There will be several other advantages to the additional communication such as being able to estimate air quality co-benefits of mitigation.