Radiative forcing of carbon dioxide, methane, and nitrous oxide: A significant revision of the methane radiative forcing

New calculations of the radiative forcing (RF) are presented for the three main well‐mixed greenhouse gases, methane, nitrous oxide, and carbon dioxide. Methane's RF is particularly impacted because of the inclusion of the shortwave forcing; the 1750–2011 RF is about 25% higher (increasing from 0.48 W m−2 to 0.61 W m−2) compared to the value in the Intergovernmental Panel on Climate Change (IPCC) 2013 assessment; the 100 year global warming potential is 14% higher than the IPCC value. We present new simplified expressions to calculate RF. Unlike previous expressions used by IPCC, the new ones include the overlap between CO2 and N2O; for N2O forcing, the CO2 overlap can be as important as the CH4 overlap. The 1750–2011 CO2 RF is within 1% of IPCC's value but is about 10% higher when CO2 amounts reach 2000 ppm, a value projected to be possible under the extended RCP8.5 scenario.


Introduction
The radiative forcing (RF) due to changes in concentrations of the relatively well mixed greenhouse gases (WMGHGs) is the largest component of total RF due to human activity over the past century [Myhre et al., 2013a].
The headline RF values for CO 2 , CH 4 , and N 2 O presented in recent Intergovernmental Panel on Climate Change (IPCC) assessments [e.g., Myhre et al., 2013a] are calculated using simplified expressions presented in the IPCC Third Assessment Report [Ramaswamy et al., 2001, section 6.3.5]. These were largely based on the work of Myhre et al. [1998, henceforth MHSS98]. MHSS98 used updated RF calculations to modify expressions adopted by IPCC in its first assessment [Shine et al., 1990] which originated from earlier work of Hansen et al. [1988] and Wigley [1987]. The MHSS98 expressions have been widely used elsewhere, for example, for calculating WMGHG RF in simple climate models and in NOAA's Annual Greenhouse Gas Index (http://www. esrl.noaa.gov/gmd/aggi/aggi.html (accessed 11 November 2016)). Other similar simplified expressions have been proposed [e.g., Shi, 1992;Byrne and Goldblatt, 2014], and Hansen et al. [2000] presented updates to the Hansen et al. [1988] fits.
The purpose of this letter is to update these expressions in a number of important ways. MHSS98 used broadband and narrowband radiation transfer calculations, which had been evaluated against clear-sky line-by-line (LBL) calculations, as a basis for the simplified expressions. Here all-sky LBL calculations are used directly. Since MHSS98, there have been many updates to spectroscopic data (both spectral lines and continua) which need to be incorporated. In addition, MHSS98 incorporated the shortwave component of the CO 2 RF; more recently, the potential importance of the shortwave bands of CH 4 has started to be recognized [Collins et al., 2006;Li et al., 2010] and we demonstrate the significance of these bands for RF here. The fits are also extended to cover a wider range of concentrations, including those relevant to the glacial-interglacial cycles over the past 800,000years. Finally, we now incorporate the effect of the CO 2 -N 2 O spectral overlap [see, e.g., Byrne and Goldblatt, 2014] which was not included in MHSS98. Byrne and Goldblatt [2014] have also presented simplified expressions based on high-spectral-resolution calculations, but there are restrictions in using their expressions. First, they reported instantaneous forcings (see Myhre et al. [2013a] for detailed definitions) and so neglected the role of stratospheric temperature adjustment on RF (see below). Second, the form of their expressions (which are referenced to preindustrial concentrations) means that they are less easy to apply to arbitrary changes in WMGHG concentrations over their range of applicability. They do, however, extend to a wider range of concentrations, making them potentially applicable to gas concentrations several million years ago.
The RF calculations here include adjustment to stratospheric temperatures, which is known to improve the utility of RF as a predictor of surface temperature change in climate models, compared to the instantaneous forcing [e.g., Myhre et al., 2013a]. Myhre et al. [2013a] presented the concept of effective radiative forcing (ERF), which includes tropospheric adjustments (of, for example, temperatures and clouds) with surface temperatures held fixed. ERF is a conceptually more applicable predictor of surface temperature response across a wide range of forcing mechanisms. However, it would be a formidable job to calculate ERFs for the range of cases presented here; it would require the use of climate models, with low-spectral-resolution radiative transfer codes, and the tropospheric adjustment would likely be significantly model dependent. Myhre et al. [2013a] concluded that the central estimates of ERFs for the WMGHGs considered here were the same as the stratosphere-adjusted RFs, but with a larger uncertainty. Unlike MHSS98, we do not include RF due to halogenated gases, as these have been covered in detail by Hodnebrog et al. [2013].

Methods
The RF calculations use the Oslo line-by-line (OLBL) code which was described in Myhre et al. [2006]. OLBL has contributed, as a benchmark model, to several radiation code intercomparisons including Forster et al. [2005Forster et al. [ , 2011, Myhre et al. [2009], and Randles et al. [2013]. We use the Myhre et al. [2006] two-atmosphere (tropical/ extratropical) approach, which includes cloud data from the International Satellite Cloud Climatology Project, to generate global forcings. OLBL code is based on the GENLN2 LBL code, which was used to generate the clear-sky forcings presented in MHSS98, but it has been enhanced in a number of important ways. It is coupled to a 16-stream Discrete Ordinate code [Stamnes et al., 1988] to compute irradiances and now includes clouds and stratospheric temperature adjustment using the standard fixed-dynamical heating methodology. The shortwave RF part of OLBL is now updated to be representative for global simulations. The solar irradiance spectrum is taken from Lean et al. [2005]. Solar radiative transfer simulations are performed for five solar zenith angles, and day length fraction is calculated for the tropical and extratropical profiles. Present-day natural and anthropogenic aerosols are included using the OsloCTM2 simulations for AeroCom Phase II [Myhre et al., 2013b] with broadband surface albedos taken from the same simulations. Absorption data from HITRAN 2008 edition [Rothman et al., 2009] are adopted both for the longwave and shortwave RF calculations. During this work the water vapor continuum was also updated (see section 3.2 for details).
We use 48 combinations (see supporting information Table S1) of CO 2 , N 2 O, and CH 4 concentrations covering the range from the lowest typical values in the 800,000year ice core record (180ppm of CO 2 , 200ppb of N 2 O, and 340ppb of CH 4 ) [Masson-Delmotte et al., 2013] to the highest anticipated amount in the year 2300 (2000 ppm of CO 2 , 525ppb of N 2 O, and 3500ppb of CH 4 ) [Meinshausen et al., 2011].

Results
The RF results for all 48 cases are presented in Table S1 and include the shortwave RF, the instantaneous clear-sky and all-sky longwave RF, the adjusted longwave RF (which includes the influence of the shortwave forcing on stratospheric temperature), and the net forcing. The row numbers in Table S1 will be used to refer to particular cases; forcing values for the cases where only one gas changes, and the impact of overlap, are shown in Figure 2, which will be discussed in section 4. The change, relative to the MHSS98 expressions, is also shown in Table S1. Updated expressions are discussed in section 4. The new calculations are, for the most part, within about 5% of the old fits, but with a few notable exceptions. In particular, for cases where only methane changes (Table S1, rows 4, 7, and 10), the new calculations are 17-27% higher than the old fitsthe reasons for this will be discussed in section 3.2. For high CO 2 cases (2000ppm-Table S1, rows 37-48) the new calculations are typically 10% higher than the old fits, indicating that the CO 2 forcing increases more rapidly than expected from a purely logarithmic dependence (as noted by others [e.g., Hansen et al., 1988;Zhong and Haigh, 2013], and see section 4). Such high CO 2 amounts were outside the range considered by MHSS98. The ice age CO 2 concentrations (180ppm-Table S1, rows 13-24) were also outside the range considered in MHSS98, but nevertheless, the old fits are within about 5% of the new calculations, which indicates that the pure logarithmic dependence was more appropriate for lower CO 2 levels. A few cases in Table S1 Geophysical Research Letters 10.1002/2016GL071930 (rows 6, 8, and 12) have a high error (exceeding 50%) where the forcing is a small residual due to an increase in CH 4 and a decrease in N 2 O, or vice versa.

Explanation of Methane Forcing Changes
The most striking feature of the results is the enhancement of the methane RF by about 17-27% compared to the old expressions; this is beyond the nominal uncertainty estimates (about 10%) given in successive IPCC reports for WMGHG forcings. Two mechanisms are responsible-it is primarily due to the role of methane's shortwave bands, which were not included in MHSS98, with a secondary effect of an update to the water vapor continuum strength relative to that used in MHSS98.
Considering, as an example, the CH 4 change, from 750 to 1800ppb (Table S1, row 4), the instantaneous shortwave forcing due to methane (0.03Wm −2 ) is about 6% of the total methane forcing (0.58Wm −2 ). This is further enhanced because the shortwave forcing affects stratospheric temperature adjustment; this effect is manifested by the change between the instantaneous and adjusted longwave forcings. Without the shortwave effect, the adjustment is found to be small and negative; the instantaneous RF of 0.516Wm −2 decreases by 2% to 0.504Wm −2 after adjustment for the 750-1800ppb CH 4 change. Rind and Lacis [1993] report a 0.5% decrease after stratosphere temperature adjustment for a 280-560ppb change in CH 4 , using a radiativeconvective model. When methane's shortwave absorption is included, the longwave adjustment is positive, increasing the longwave forcing from its instantaneous RF value of 0.516 to 0.548Wm −2 after adjustment, a 6% increase (or 9% relative to adjusted forcing without the shortwave effect). Hence, in total, the shortwave effect increases the methane forcing by 15% from its adjusted longwave-only value of 0.504 to 0.582Wm −2 .
The fact that the shortwave forcing for CH 4 is positive is notable, in the context of earlier work. For CO 2 , the shortwave bands cause a forcing that opposes and so acts to decrease the magnitude of the CO 2 longwaveonly forcing by about 5% (e.g., MHSS98 and Table S1, rows 13, 25, and 37); this is due to the combined effect of the direct absorption of the solar radiation in the stratosphere and its subsequent effect on temperature adjustment. For CH 4 , the instantaneous forcing calculations of Collins et al. [2006] and Forster et al. [2011] also imply a negative, rather than a positive shortwave forcing; in Collins et al. [2006] the shortwave forcing is 15% of, but the opposite sign to, the longwave irradiance change at 200hPa for a doubling of CH 4 from preindustrial concentrations. However, these calculations were for an idealized clear-sky case with a fixed solar zenith angle of 53°and a surface albedo of 0.1. The shortwave bands of N 2 O contribute less than 1% of the forcing; this forcing is included here but discussed no further.
The contrasting signs of the shortwave forcings of CH 4 and CO 2 can be explained by reference to Figure 1. This shows the spectral variation of global-mean shortwave net forcing for the two gases (180ppm to 389 ppm for CO 2 and 750 to 1800ppb for CH 4 ). Supporting information Figure S1 shows the upward and downward forcing components. For both gases (Figures 1a and 1b), the sign of the forcing varies with wavelength, the net impact being the residual of these. In the case of CO 2 the negative forcing due to its 2.7µm band dominates. For CH 4 , the positive forcing due to its 1.7 and 2.3µm bands dominates. This contrasting behavior is driven by two processes. One is the stratospheric opacity of the gases, and the other is the degree of overlap of absorption bands with the near-IR bands of water vapor. Figures 1c and 1d show the sum of line strengths for each 1nm interval in the OLBL for water vapor (in both plots) and CO 2 and CH 4 .
For all bands, the downward shortwave flux at the tropopause is always decreased by the increased concentrations, due to increased absorption in the stratosphere (Figures S1a and S1b). The sign of the forcing depends on whether this negative contribution dominates over the increased absorption by these gases in the troposphere, which contributes a positive forcing. For CO 2 , the extremely strong band at 4.3µm makes little contribution to forcing at band center; absorption is almost complete in the stratosphere at unperturbed concentrations, and it only starts to contribute to forcing at the band edges. The CO 2 2.7µm band lies toward the center of a strong water vapor band (Figure 1d). The change in the downward forcing is strongly negative ( Figure S1b), but there is not a compensating increase in tropospheric absorption because the heavy spectral overlap with water vapor strongly mutes the impact of CO 2 increases. By contrast, the weaker bands of CO 2 at 1.6 and 2.0µm lie in, or toward the edges of, windows between the main water vapor bands ( Figure 1d); hence, they are more able to increase tropospheric absorption causing a positive forcing that dominates over the negative stratospheric component.

10.1002/2016GL071930
Methane has a strong band at 3.3µm (Figure 1a), which lies within a region of relatively strong water vapor absorption (Figure 1c), leading to a negative net forcing. The weaker bands at 1.6 and 2.3µm lie toward the center of the windows in the water vapor spectrum. They cause a positive forcing which more than compensates for the 3.3µm negative forcing. In summary, the magnitude and sign of the solar RF depend on band strength, gas concentration, and overlap with water vapor.
The contrast with the implied negative forcing in earlier work originates from clouds. For the Collins et al. [2006] clear-sky case, there is little tropospheric upward scatter of radiation (no clouds and low surface albedo), so that photons largely take a single pass through the troposphere. When clouds are added, the added CH 4 absorbs not only downwelling radiation but also upward scattered radiation, greatly increasing the tropospheric absorption (see Table S2). For both the extratropical and tropical cases, the inclusion of clouds causes the shortwave RF to change from a negative to a positive forcing of almost the same size; the negative clear-sky forcing of about −0.04Wm −2 is consistent with the Collins et al. [2006] value of −0.13Wm −2 (for approximately the same methane change) after accounting for the higher incoming solar irradiance for their (fixed-Sun) case. The effect of clouds on the shortwave RF from CH 4 is similar to the effect of clouds on the shortwave RF of black carbon aerosols [Haywood and Shine, 1997].
In trying to reconcile the present (longwave) results with those in MHSS98, we identified a further influence on the forcing due to CH 4 (and, to a lesser extent, the N 2 O). MHSS98 used the Clough-Kneizys-Davies (CKD) water vapor continuum version 0 in the OLBL calculations. In the region of the CH 4 and N 2 O bands that are most responsible for the longwave forcing (around 1300cm −1 ), the foreign continuum was weakened by about a factor of 3.5 in subsequent versions of CKD [e.g., Mlawer et al., 1998] and successor versions (the Mlawer-Tobin-CKD) [Mlawer et al., 2012]. This reduces the effect of water vapor overlap with CH 4 , increasing its RF. Changes in the self-continuum in this spectral region were much smaller during these updates. To isolate the effect of the continuum changes, the instantaneous cloudy-sky longwave forcings were examined when CKD version 0 is updated to version 2.4.1. The update caused the CH 4 forcing for a 1800 to 3500ppb change to increase by 4.2%. For comparison, the N 2 O forcing increased by about 1.7% for a 323 to 525ppb change; for CO 2 , RF changed by less than 0.2%.
In total, comparing the new RFs to those calculated with the MHSS98 expressions, the CH 4 forcing increases by 17-27% for cases 4, 7, and 10 shown in Table S1. This is mostly due to the shortwave forcing, with a smaller contribution from the decreased strength of the foreign continuum in the region of methane's longwave bands. We note that contemporary radiation codes used in climate models often neglect the shortwave band of methane (e.g., all those models participating in the Collins et al. [2006] and CCMval [Fomichev and Forster, 2010] intercomparisons); by contrast, it is likely that more recent versions of the water vapor continuum in the mid-infrared are already in use in these codes, and so they will be less affected by this component of our update to the MHSS98 calculations.
Based on our findings, the neglect of the CH 4 shortwave bands would underestimate the total historical anthropogenic RF by about 0.1Wm −2 and hence underestimate temperature change. However, many climate models have too high longwave RF of methane, which could (inadvertently) compensate, to some extent, for the lack of solar absorption by methane [Collins et al., 2006]. The same is true for models in the CCMval comparison [Fomichev and Forster, 2010], where the majority of model configurations gave too high a methane longwave forcing, some by as much as 30% (E. Rozanov, personal communication, 2016). Solar absorption by methane has not been included in observationally based estimates of climate sensitivity [e.g., Skeie et al., 2014;Johansson et al., 2015], and thus, these have used a too low total historical RF; however, the change in methane RF is relatively small compared to uncertainties in many of the non-greenhouse gas RFs and in particular those related to aerosol RF. Finally, we note that the strong spectral variation of the sign of the CH 4 forcing (Figure 1), and its strong dependence on the overlap with water vapor, may pose a challenge to its accurate representation in climate model radiation codes, with relatively low spectral resolution.

CO 2 -N 2 O Overlaps
Nitrous oxide possesses a moderately strong band centered at 17µm, which overlaps with the high-wavelength side of the 15µm complex of CO 2 absorption bands; hence, N 2 O forcing can be affected by CO 2 changes and vice versa. This overlap was neglected in most previous simplified expressions of CO 2 [e.g., Shi, 1992;Hansen et al.,1988;MHSS98;Ramaswamy et al., 2001]. It was included (but not highlighted) in Byrne and Goldblatt's [2014] calculations of cloudy-sky instantaneous RF. Table S1, the N 2 O forcing for a change from 200 to 525ppb is diminished by about 14% when using a CO 2 value of 2000ppm (0.83Wm −2 , differencing rows 38 and 39) rather than 180ppm (0.97Wm −2 , differencing rows 14 and 15), assuming CH 4 at 1800ppb. This is almost the same as the effect on N 2 O RF but using a CH 4 value of 3500ppb (0.89Wm −2 , differencing rows 11 and 12) rather than 340ppb (1.02W m −2 , differencing rows 8 and 9), assuming CO 2 at 389ppm. Hence, when calculating the N 2 O forcing, coincident changes in CO 2 can be as important as coincident changes in CH 4 . The CO 2 forcing for a 180 to 2000ppm change is diminished by about 1% if N 2 O values of 525ppb (differencing rows 15 and 39) are used, rather than 200ppb (differencing rows 14 and 38), assuming CH 4 at 1800ppb. By contrast, the CO 2 forcing changes by less than 0.1% across the range of CH 4 concentrations investigated here (for example, differencing rows 46 and 19, versus rows 43 and 22).

Uncertainties
Successive IPCC assessments [see Myhre et al., 2013a, section 8.3.1] have drawn on a range of evidence to estimate an uncertainty range for the WMGHG RFs of ±10%, when calculated using line-by-line codes. Hodnebrog et al. [2013] (see especially their section 3.6) presented an extended discussion of sources of uncertainty in forcing, based on the prior literature, including intercomparison studies. Although this was in the context of halocarbon radiative forcing, much of the discussion carries over to the longwave component of the WMGHG RFs discussed here.
There are distinct sources of uncertainty, including the underlying spectroscopic data, the calculation methodology, and the specification of atmospheric conditions (and the averaging of these conditions). Assessment of the impact of the then recent updates of the HITRAN database on longwave forcing [e.g., Kratz, 2008], and a more detailed assessment specifically for CO 2 [Mlynczak et al., 2016], indicates uncertainties of less than 1%. Hodnebrog et al. [2013] estimate that uncertainties from the radiation code and the method of calculating stratospheric adjustment are 5 and 4%, respectively, while specification of atmospheric temperature, clouds, and tropopause position each contributes 3%, 5%, and 4%, respectively, and that due to temporal and spatial averaging is about 1%. Their overall conclusion was that the total (5-95%) uncertainty was ±13% for halogenated gases; this is slightly higher than the Myhre et al. [2013a] Geophysical Research Letters 10.1002/2016GL071930 estimate of 10% for all WMGHGs, but the increase is specific to the halocarbons where uncertainties in spectroscopic data and the vertical profile are higher than the gases considered here.
We are not aware of such a detailed assessment for WMGHG forcing in the shortwave. Collins et al. [2006] show good agreement among four line-by-line codes for the shortwave forcing of greenhouse gases, for a specified clear-sky atmospheric state, and a single zenith angle; for example, at 200hPa, the standard deviation was about 2% of the mean for CO 2 and a water vapor perturbation and no more than 4% (given the precision that their results were reported) for CH 4 . However, now that our calculations demonstrate an enhanced role for methane in this spectral region, a more detailed assessment becomes important. Rothman et al. [2013] note incomplete understanding of CH 4 line shapes particularly in the near-infrared; the impact of future HITRAN versions on calculations of the CH 4 shortwave forcing will need to be assessed. Toon et al.
[2016] note a number of instances in this wavelength region where balloon-borne spectral measurements are in better agreement with versions of HITRAN prior to the 2012 edition. In addition, there is ongoing uncertainty in the strength of the water vapor continuum in near-IR windows, particularly at atmospheric temperatures, [e.g., Shine et al., 2016], where the CH 4 shortwave forcing is important (Figure 1).
Other error sources for computing shortwave forcing, beyond those considered for the longwave, include those resulting from the specification of surface albedo and aerosol concentrations and the integration over day length to produce a day-averaged forcing; these will also be affected by the temporal and spatial averaging used. A number of sensitivity tests of the solar forcing were performed here for CH 4 changes from 750 to 1800ppb, with all other gases held fixed. The complete removal of aerosol reduces the shortwave forcing by 10%, as the lower albedo means less reflected solar radiation to be absorbed by CH 4 . Similarly, increasing the surface albedo by 1 percentage point (from a global-mean value of 12.3 to 13.3%) increased the shortwave forcing by 6.7%. These sensitivity studies indicate that there is some interdependence between the CH 4 shortwave forcing and other forcings, which could be examined in more detail in future work. Increasing the number of solar zenith angles at which the forcing is calculated from 5 to 15 changed the forcing by only 0.2%.
We retain the Myhre et al. [2013a] ±10% uncertainty range for the longwave forcing and hence use this for the total uncertainty for CO 2 and N 2 O, given the small contribution from the shortwave forcing for these gases. For CH 4 we tentatively adopt a shortwave forcing uncertainty of about double this value (±25%), but given that the shortwave forcing is about 15% of the total forcing (including its contribution to the longwave adjustment discussed in section 3.2), this yields a total uncertainty in CH 4 forcing of about ±14%.

Simplified Expressions
Based on the results presented in section 3, the MHSS98 simplified expressions are updated. While the basic forms of the equations are retained, important modifications are introduced, in order to achieve improved agreement with the detailed calculations, over the range of concentration changes. The new expressions, derived using polynomial fitting, are shown in Table 1 and are valid for the CO 2 , N 2 O, and CH 4 concentrations given in section 2.
For CO 2 , the logarithmic form used in MHSS98 (RF=α ln(C/C o ), where C o and C are the initial and final CO 2 concentrations) is retained but the nature of α term is changed; instead of being a constant, it is now a function of the CO 2 and N 2 O concentrations. Note that the absolute value of CO 2 concentration change |C-C o | is adopted in the α term, to ensure that the forcing is symmetric for increases or decreases in CO 2 . That is, RF(C, C o )=−RF(C o ,C). The C-C o terms only become important for large changes in CO 2 ; as shown below, for historical forcings, ignoring these terms and using the simple α ln(C/C o ) form, and a midrange N 2 O concentration, yields a value of α within 1% of the value of 5.35Wm −2 given in MHSS98.
For N 2 O and CH 4 , the basic square root dependence of RF used in MHSS98 (e.g., for N 2 O, RF=β(√N−√N o ) plus an overlap term) is retained, but with two changes: the β term is no longer a constant, and it also now incorporates the overlap terms. MHSS98 adopted the Hansen et al. [1988] overlap term; this had the useful feature that it was the same for both N 2 O and CH 4 , but it is a complicated expression including the logarithm of products of N 2 O and CH 4 raised to noninteger powers. We find that a different approach, with separate N 2 O and CH 4 overlap terms, enables the overlap to be included by modifying the β term. In addition, as noted in section 3, it is necessary to include a N 2 O-CO 2 overlap term.

10.1002/2016GL071930
In all cases, the overlap term is treated differently to previous practice. In most previous work including MHSS98, the overlap term has been treated by considering the final concentration of the gas for which the RF is being calculated (e.g., M, in the case of CH 4 ) and the initial concentration of overlapping gas (e.g., N o when calculating the CH 4 RF and M o when calculating the N 2 O RF). By contrast Myhre et al. [2013a] used the final concentrations of the overlapping gases (e.g., N when calculating CH 4 RF). One undesirable aspect of both approaches is that the calculated RF is not symmetric about increases and decreases in the gas concentration (e.g., RF(M,M o )≠−RF(M o ,M)) if the overlapping gas changes in concentration. To alleviate this problem, the mean concentration (e.g., Þ ) is adopted here for all terms (other than the (C-C o ) terms) within the square brackets in Table 1 equations-the quality of the fits is, in the majority of cases, superior to using the final concentrations.
The final columns of Table S1 show the RF calculated using the new expressions, for each gas separately, and their sum, and the error in the sum relative to the OLBL calculations. Except for two cases where the CH 4 and N 2 O forcings strongly oppose each other (Table S1, rows 6 and 8, and see section 3.1), the new expressions reproduce the OLBL forcings to within better than 5% across the range of concentration changes. In a few cases, the old fits are better, but the most severe cases (Table S1, rows 16, 18, and 24) are for unlikely combinations with very low paleo-CO 2 and present or future concentrations of N 2 O. Even in the worst of these cases, the new fits reproduce the OLBL calculations to within better than 4%; however, taken as a whole, the new fits are considerably better.
The results are summarized in Figure 2. Figures 2a-2c show the individual forcings for CO 2 , CH 4 , and N 2 O, comparing the new expressions with the OLBL results. Figure 2d illustrates the strength of the overlap terms, and their representation by the new expressions, by showing the difference in forcing for a given gas, between the highest and lowest concentrations of the overlapping gases; since this range covers the entire ice age to maximum projected 2300 values, it represents an extreme test. For CO 2 the difference in forcing between the maximum and minimum N 2 O concentrations is at the 1% level, and for CH 4 the effect of N 2 O overlap is at less than the 10% level. For N 2 O, both CO 2 and CH 4 overlaps contribute about equally, for the specified concentration changes, with the difference exceeding 10% for each gas.

Impact on Historical and Future Forcings and Emission Metrics
The new expressions increase the IPCC AR5 [Myhre et al., 2013a, Table 8.2] RFs for CO 2 , CH 4 , and N 2 O for the period 1750-2011 from 1.82, 0.48, and 0.17Wm −2 to 1.83, 0.61, and 0.17Wm −2 or by 0.5%, 25%, and 2%, respectively; the difference in the sum of the three forcings is 0.14Wm −2 . Provisional values for 2015 global-annual-mean concentrations (from http://www.esrl.noaa.gov/gmd/-accessed 3 October 2016) for CO 2 (399ppm), CH 4 (1834ppb), and N 2 O (328ppb) (using the same preindustrial values of 278ppm, 722ppb,  and 270ppbv, respectively) yield RFs of 1.95, 0.62, and 0.18Wm −2 ; this gives an increased RF for these three gases of 0.14Wm −2 relative to 2011 concentrations, using the new expressions. Figure 3a compares the RF derived from old and updated expressions for 1750-2011; for the new expressions, the shaded area shows the radiative uncertainty discussed in section 3.4. It shows that for CO 2 the RF using the old and new expressions differs by only about 1%. As anticipated from the results in section 3, this difference is larger for the other two gases; in 2011 it is about 2% for N 2 O forcing and 25% for CH 4 forcing. Figure 3b shows RF using the extended RCP8.5 out to 2300. For CO 2 , the percentage difference between RF values using old and new expressions increases to about 9%, illustrating the importance of the nonconstant α term. For N 2 O and CH 4 RF, the difference between the new and old expressions in 2300 reaches −21% and 12%, respectively. Thus, in future scenarios the new expressions suggest that forcing will be even higher than current estimates due to the increase in CO 2 and CH 4 RF, which is slightly offset by the decrease in N 2 O forcing, which is due to CO 2 overlap at high CO 2 concentrations.
The calculations of metrics (global warming potential (GWP) and global temperature change potential (GTP)) presented in Myhre et al. [2013a] are also impacted. These utilize radiative efficiencies (REs) for small perturbations around present-day concentrations. For CO 2 and N 2 O these are affected by about 1% or less compared to the values presented in Table 8.A1 of Myhre et al. [2013a]. By contrast, the CH 4 radiative efficiency (RE), for small perturbations about present-day concentrations, increases from 3.63×10 −4 to 4.48×10 −4 Wm −2 ppb −1 , an increase of 23%. This percentage difference is slightly lower than for the industrial era CH 4 change because differences in the updated expression relative to MHSS98 are largest for low CH 4 concentrations (see Table S1). Since the GWP and GTP values for CH 4 in Myhre et al. [2013a] include indirect effects due to ozone and stratospheric water vapor change, and the absolute contribution of these is unchanged by the increase in the RE, the metrics themselves The forcings are relative to 389ppm for CO 2 , 1800ppb for CH 4 , and 323ppb for N 2 O; the effect of overlapping gases is included using these same mixing ratios. (d) Illustration of the effect of the overlap on forcings for OLBL (symbols) and simple expressions (lines); the difference between the forcing for the highest and lowest concentrations of the overlapping gas is plotted: N 2 O overlap (525 minus 200ppb), CH 4 overlap (3500 minus 340ppb), and CO 2 overlap (2000 minus 180ppm). The effect of CH 4 and CO 2 overlaps of N 2 O are almost the same for these choices of mixing ratio, so that these lines are almost coincident on the plot.

10.1002/2016GL071930
increase by about 14%. The GWP for the 100year time horizon, the most commonly used metric, increases from 28 to 32.

Conclusions
New line-by-line calculations of RF due to CO 2 , CH 4 , and N 2 O are used as a basis to update the simple expressions in MHSS98 that have been used widely, including in IPCC assessments, to calculate WMGHG forcings. While the changes to CO 2 and N 2 O RF between the old and new fits are within the 10% uncertainty estimate given in Myhre et al. [2013a], the changes in the CH 4 RF are significantly larger. We show that inclusion of the shortwave (near-infrared) bands of methane, which were neglected in many previous analyses, is most responsible for this change; there is a smaller contribution due to an update of the water vapor foreign continuum relative to that used in MHSS98. The new expressions reproduce the line-by-line calculations to better than 5%, but the quality of the fits must be distinguished from the uncertainty in the forcing calculations themselves-based on previous work and calculations performed here, we assess the uncertainty to be ±10% for CO 2 and N 2 O and ±14% for CH 4 . The higher uncertainty for CH 4 arises from the greater importance of the shortwave forcing. The results emphasize the need for continuing assessment and refinement of RF calculations used in deriving such simple expressions. Further assessment of the methane shortwave effect is particularly warranted, as is continued examination of the fundamental spectroscopic basis of greenhouse gas radiative forcing calculations [e.g., Mlynczak et al., 2016].