
Fig. 1
GHG forcing as presented by the Goddard Institute for Space Studies (GISS; http://data.giss.nasa.gov/modelforce/RadF.txt), National Oceanic and Atmospheric Administration (NOAA; http://www.esrl.noaa.gov/gmd/aggi/AGGI_Table.csv) and the Representative Concentration Pathways group (RCP; http://www.pik-potsdam.de/~mmalte/rcps/data/20THCENTURY_MIDYEAR_RADFORCING.xls). All forcings are set equal at 1982 to permit comparison.

Fig. 2
Correlation of global temperature and GHG forcing. Temperature anomaly data are HadCrut3 (Brohan et al., 2006, as extended at http://www.cru.uea.ac.uk/cru/data/temperature/). Forcing (blend of RCP and NOAA as discussed in the text) is relative to preindustrial. Ratio of scales of two vertical axes was set by slope of graph of ΔT s vs. forcing.

Fig. 3
Forcing by LLGHGs and non-LLGHG forcing over the time period 1970–2010 as given by the GISS and blended RCP-NOAA data sets.

Fig. 4
Graph of temperature anomaly vs. forcing by LLGHGs for the years 1970–2010 (indicated by colour). Forcing is the average of GISS and blended RCP-NOAA, relative to 1970, Fig. 3. Slope S tr=0.39±0.03 K/(W m−2), where the 1−σ uncertainty is based only on the uncertainty in the fit; forcing is relative to 1970; temperature anomaly HadCrut3 is relative to base period 1961–1990. Correlation coefficient r 2=0.80.
Table 1. Calculation of lower-bound transient and equilibrium sensitivities
[i] LLGHG forcing over period 1970–2010 ΔF is based on the mean of GISS and blended RCP-NOAA forcing data sets; uncertainty in forcing is taken as ±16% as discussed in text. Column 3 presents values for forcing by LLGHGs only and thus yields a best estimate for lower-bound transient and equilibrium sensitivity. Uncertainty in S tr reflects uncertainties in ΔF and dΔT s/dΔF. Heating rate N and associated uncertainty are from Table 2. Time range for ΔT s is for middle of time period examined relative to assumed steady state at beginning of twentieth century. Last column shows lower bounds of the 5–95% uncertainty range, evaluated as the best-estimate value of the lower bound minus 1.64 times the 1−σ uncertainty for the probability distribution function for the quantity taken as normally distributed. Values of S tr and S eq expressed in the unit K (3.7 W m−2)−1 are shown to permit comparison with commonly reported CO2 doubling temperature ΔT 2×.

Fig. 5
Heat content of the world ocean to depth of 2000 m. Slope (0.48±0.02×1022 J yr−1) of linear fit (blue) to data for years 1970–2008, indicated by arrows, corresponds to heating rate relative to the area of the planet N=0.30±0.01 W m−2. Data from Levitus et al. (2012).
Table 2. Contributions to planetary heating rate
