We describe a new method for evaluating the radiative forcing, the climate feedback parameter (W m −2 K −1 ) and hence the effective climate sensitivity from any GCM experiment in which the climate is responding to a constant forcing. The method is simply to regress the top of atmosphere radiative flux against the global average surface air temperature change. This method does not require special integrations or off‐line estimates, such as for stratospheric adjustment, to obtain the forcing, and eliminates the need for double radiation calculations and tropopause radiative fluxes. We show that for CO 2 and solar forcing in a slab model and an AOGCM the method gives results consistent with those obtained by conventional methods. For a single integration it is less precise but since it does not require a steady state to be reached its precision could be improved by running an ensemble of short integrations.

Abstract The Atlantic thermohaline circulation (THC) is an important part of the earth's climate system. Previous research has shown large uncertainties in simulating future changes in this critical system. The simulated THC response to idealized freshwater perturbations and the associated climate changes have been intercompared as an activity of World Climate Research Program (WCRP) Coupled Model Intercomparison Project/Paleo-Modeling Intercomparison Project (CMIP/PMIP) committees. This intercomparison among models ranging from the earth system models of intermediate complexity (EMICs) to the fully coupled atmosphere–ocean general circulation models (AOGCMs) seeks to document and improve understanding of the causes of the wide variations in the modeled THC response. The robustness of particular simulation features has been evaluated across the model results. In response to 0.1-Sv (1 Sv ≡ 106 m3 s−1) freshwater input in the northern North Atlantic, the multimodel ensemble mean THC weakens by 30% after 100 yr. All models simulate some weakening of the THC, but no model simulates a complete shutdown of the THC. The multimodel ensemble indicates that the surface air temperature could present a complex anomaly pattern with cooling south of Greenland and warming over the Barents and Nordic Seas. The Atlantic ITCZ tends to shift southward. In response to 1.0-Sv freshwater input, the THC switches off rapidly in all model simulations. A large cooling occurs over the North Atlantic. The annual mean Atlantic ITCZ moves into the Southern Hemisphere. Models disagree in terms of the reversibility of the THC after its shutdown. In general, the EMICs and AOGCMs obtain similar THC responses and climate changes with more pronounced and sharper patterns in the AOGCMs.

We quantify forcing and feedbacks across available CMIP5 coupled atmosphere‐ocean general circulation models (AOGCMs) by analysing simulations forced by an abrupt quadrupling of atmospheric carbon dioxide concentration. This is the first application of the linear forcing‐feedback regression analysis of Gregory et al. (2004) to an ensemble of AOGCMs. The range of equilibrium climate sensitivity is 2.1–4.7 K. Differences in cloud feedbacks continue to be important contributors to this range. Some models show small deviations from a linear dependence of top‐of‐atmosphere radiative fluxes on global surface temperature change. We show that this phenomenon largely arises from shortwave cloud radiative effects over the ocean and is consistent with independent estimates of forcing using fixed sea‐surface temperature methods. We suggest that future research should focus more on understanding transient climate change, including any time‐scale dependence of the forcing and/or feedback, rather than on the equilibrium response to large instantaneous forcing.

Geophysical Research LettersVolume 38, Issue 18 OceansOpen Access Revisiting the Earth's sea-level and energy budgets from 1961 to 2008 Correction(s) for this article Correction to “Revisiting the Earth's sea-level and energy budgets from 1961 to 2008” John A. Church, Neil J. White, Leonard F. Konikow, Catia M. Domingues, J. Graham Cogley, Eric Rignot, Jonathan M. Gregory, Michiel R. Broeke, Andrew J. Monaghan, Isabella Velicogna, Volume 40Issue 15Geophysical Research Letters pages: 4066-4066 First Published online: August 8, 2013 John A. Church, John A. Church [email protected] Centre for Australian Weather and Climate Research and Wealth from Oceans Flagship, CSIRO Marine and Atmospheric Research, Hobart, Tasmania, AustraliaSearch for more papers by this authorNeil J. White, Neil J. White Centre for Australian Weather and Climate Research and Wealth from Oceans Flagship, CSIRO Marine and Atmospheric Research, Hobart, Tasmania, AustraliaSearch for more papers by this authorLeonard F. Konikow, Leonard F. Konikow U.S. Geological Survey, Reston, Virginia, USASearch for more papers by this authorCatia M. Domingues, Catia M. Domingues Antarctic Climate and Ecosystems Cooperative Research Centre, CSIRO Marine and Atmospheric Research, Aspendale, Victoria, AustraliaSearch for more papers by this authorJ. Graham Cogley, J. Graham Cogley Department of Geography, Trent University, Peterborough, Ontario, CanadaSearch for more papers by this authorEric Rignot, Eric Rignot Department of Earth System Science, University of California, Irvine, California, USA Jet Propulsion Laboratory, California Institute of Technology, Pasadena, California, USASearch for more papers by this authorJonathan M. Gregory, Jonathan M. Gregory NCAS-Climate, Department of Meteorology, University of Reading, Reading, UK Met Office Hadley Centre, Exeter, UKSearch for more papers by this authorMichiel R. van den Broeke, Michiel R. van den Broeke Institute for Marine and Atmospheric Research, Utrecht University, Utrecht, NetherlandsSearch for more papers by this authorAndrew J. Monaghan, Andrew J. Monaghan National Centre for Atmospheric Research, Boulder, Colorado, USASearch for more papers by this authorIsabella Velicogna, Isabella Velicogna Department of Earth System Science, University of California, Irvine, California, USA Jet Propulsion Laboratory, California Institute of Technology, Pasadena, California, USASearch for more papers by this author John A. Church, John A. Church [email protected] Centre for Australian Weather and Climate Research and Wealth from Oceans Flagship, CSIRO Marine and Atmospheric Research, Hobart, Tasmania, AustraliaSearch for more papers by this authorNeil J. White, Neil J. White Centre for Australian Weather and Climate Research and Wealth from Oceans Flagship, CSIRO Marine and Atmospheric Research, Hobart, Tasmania, AustraliaSearch for more papers by this authorLeonard F. Konikow, Leonard F. Konikow U.S. Geological Survey, Reston, Virginia, USASearch for more papers by this authorCatia M. Domingues, Catia M. Domingues Antarctic Climate and Ecosystems Cooperative Research Centre, CSIRO Marine and Atmospheric Research, Aspendale, Victoria, AustraliaSearch for more papers by this authorJ. Graham Cogley, J. Graham Cogley Department of Geography, Trent University, Peterborough, Ontario, CanadaSearch for more papers by this authorEric Rignot, Eric Rignot Department of Earth System Science, University of California, Irvine, California, USA Jet Propulsion Laboratory, California Institute of Technology, Pasadena, California, USASearch for more papers by this authorJonathan M. Gregory, Jonathan M. Gregory NCAS-Climate, Department of Meteorology, University of Reading, Reading, UK Met Office Hadley Centre, Exeter, UKSearch for more papers by this authorMichiel R. van den Broeke, Michiel R. van den Broeke Institute for Marine and Atmospheric Research, Utrecht University, Utrecht, NetherlandsSearch for more papers by this authorAndrew J. Monaghan, Andrew J. Monaghan National Centre for Atmospheric Research, Boulder, Colorado, USASearch for more papers by this authorIsabella Velicogna, Isabella Velicogna Department of Earth System Science, University of California, Irvine, California, USA Jet Propulsion Laboratory, California Institute of Technology, Pasadena, California, USASearch for more papers by this author First published: 16 September 2011 https://doi.org/10.1029/2011GL048794Citations: 378AboutSectionsPDF ToolsRequest permissionExport citationAdd to favoritesTrack citation ShareShare Give accessShare full text accessShare full-text accessPlease review our Terms and Conditions of Use and check box below to share full-text version of article.I have read and accept the Wiley Online Library Terms and Conditions of UseShareable LinkUse the link below to share a full-text version of this article with your friends and colleagues. Learn more.Copy URL Share a linkShare onFacebookTwitterLinkedInRedditWechat Abstract [1] We review the sea-level and energy budgets together from 1961, using recent and updated estimates of all terms. From 1972 to 2008, the observed sea-level rise (1.8 ± 0.2 mm yr−1 from tide gauges alone and 2.1 ± 0.2 mm yr−1 from a combination of tide gauges and altimeter observations) agrees well with the sum of contributions (1.8 ± 0.4 mm yr−1) in magnitude and with both having similar increases in the rate of rise during the period. The largest contributions come from ocean thermal expansion (0.8 mm yr−1) and the melting of glaciers and ice caps (0.7 mm yr−1), with Greenland and Antarctica contributing about 0.4 mm yr−1. The cryospheric contributions increase through the period (particularly in the 1990s) but the thermosteric contribution increases less rapidly. We include an improved estimate of aquifer depletion (0.3 mm yr−1), partially offsetting the retention of water in dams and giving a total terrestrial storage contribution of −0.1 mm yr−1. Ocean warming (90% of the total of the Earth's energy increase) continues through to the end of the record, in agreement with continued greenhouse gas forcing. The aerosol forcing, inferred as a residual in the atmospheric energy balance, is estimated as −0.8 ± 0.4 W m−2 for the 1980s and early 1990s. It increases in the late 1990s, as is required for consistency with little surface warming over the last decade. This increase is likely at least partially related to substantial increases in aerosol emissions from developing nations and moderate volcanic activity. Key Points The sea level budget is closed from 1972 Oceans are continuing to warm and sea level toi rise Aerosol loading has increased since 2000 1. Introduction [2] Sea-level rise is a high-profile and important aspect of climate change. However, in the last two Intergovernmental Panel on Climate Change (IPCC) assessments, the sum of observed contributions to sea-level rise has consistently been less than the observed rise over multi-decadal periods, thus reducing confidence in the sea-level projections. Domingues et al. [2008], and more recently Moore et al. [2011], made significant progress in closing the sea-level budget from 1961 to 2003. However, to do so, they had to estimate multi-decadal trends for the ice sheets, rather than use observed or calculated time series, and as a result they failed to reproduce the observed decadal variability in sea levels. [3] Closing the sea-level budget requires accurate estimates of ocean thermal expansion and the related ocean heat-content change, by far the largest storage of heat in the climate system [Bindoff et al., 2007]. Despite the Earth's energy budget being one of the most fundamental aspects of the climate system there have been failures to adequately resolve the energy budget for recent years [Trenberth, 2010; Trenberth and Fasullo, 2010]. [4] Clearly, the Earth's sea-level and energy budgets are closely related and must be solved in a consistent manner. We therefore consider them together, using new and updated estimates of all terms for the past five decades. By updating the time series of sea-level rise and its contributions, including a new estimate of groundwater depletion, we find an improved closure of the sea-level budget from 1972 to the present, both in the mean trends from 1972 to 2008, and also the variability. We also use these estimates to look at the implications for the energy budget calculation of Murphy et al. [2009]. 2. The Components of Sea-Level Change and Earth's Energy Storage [5] Time series from 1961 to 2008 of all of the contributions to the sea-level and energy budgets are given in Figures 1– 3. Because of the uncertainty affecting estimates of the contributions during the 1960s, linear trends in the sea-level budget (Table 1) are given for 1972 to 2008 and 1993 to 2008. Similarly, integrals of components of the energy budget from 1972 to 2008 and 1993 to 2008 are shown in Table 2. Figure 1Open in figure viewerPowerPoint Cumulative net groundwater depletion, 1961–2008. Figure 2Open in figure viewerPowerPoint The global sea-level budget from 1961 to 2008. (a) The observed sea level using coastal and island tide gauges (solid black line with grey shading indicating the estimated uncertainty) and using TOPEX/Poseidon/Jason-1&2 satellite altimeter data (dashed black line). The two estimates have been matched at the start of the altimeter record in 1993. Also shown are the various components as described in the text. (b) The observed sea level and the sum of components. The estimated uncertainties (consistent with Table 1) are indicated by the shading. The two time series are plotted such that they have the same average over 1972 to 2008. Figure 3Open in figure viewerPowerPoint The Earth's energy budget from 1961 through 2008. (a) The top line indicates the integrated radiative forcing from well-mixed greenhouse gases, ozone and solar changes. This radiative forcing is offset by the cooling due to stratospheric aerosols from volcanic eruptions (cyan), the warming of the upper (0 to 700 m) ocean (blue), the deep ocean (red) and energy absorbed in the melting of ice and warming the atmosphere and the land (green), and the change in outgoing radiation inferred from temperature change (brown). The integrated aerosol (and other) forcing is estimated as a residual (grey). (b) The inferred global averaged aerosol radiative cooling (red) is calculated for three different values of λ [1.25 W m−2 K−1 (solid) and 0.75 W m−2 K−1 and 1.75 W m−2 K−1 (dashed)]. The aerosol coolings for the central value of λ are also given for the Ishii and Kimoto [2009] and Levitus et al. [2009] heat content estimates (green and orange respectively). The equivalent results for the Gregory and Forster [2008] volcanic forcing and the present heat content estimates are shown by the blue line and shading. Results from previous aerosol forcing estimates of Gregory and Forster [2008] and GISS forcing are shown in black solid and dashed respectively. Table 1. The Sea-Level Budgeta Component 1972 → 2008 1993 → 2008 Total s.l. (t.g. only) 1.83 ± 0.18b 2.61 ± 0.55 Total s.l. (t.g. + sat) 2.10 ± 0.16 3.22 ± 0.41 Shallow thermal (0–700m) 0.63 ± 0.09 0.71 ± 0.31 Deep thermal (700–3000m) 0.07 ± 0.10 0.07 ± 0.10 Abyssal thermal (3000m-bottom) 0.10 ± 0.06 0.10 ± 0.06 Total thermal (full depth) 0.80 ± 0.15 0.88 ± 0.33 Glaciers & Ice Caps 0.67 ± 0.03 0.99 ± 0.04 Greenland Ice Sheet 0.12 ± 0.17 0.31 ± 0.17 Antarctic Ice Sheet 0.30 ± 0.20 0.43 ± 0.20 Land ice (G&IC, GIS, AIS) 1.09 ± 0.26 1.73 ± 0.27 Thermal (full depth) + Land ice 1.89 ± 0.30 2.61 ± 0.42 Dam retention −0.44 ± 0.15 −0.30 ± 0.15 Groundwater depletion 0.26 ± 0.07 0.35 ± 0.07 Natural terrestrial storage 0.07 ± 0.10 −0.14 ± 0.10 Total terrestrial storage −0.11 ± 0.19 −0.08 ± 0.19 Total mass contributions 0.98 ± 0.33 1.66 ± 0.33 Total thermal + Mass 1.78 ± 0.36 2.54 ± 0.46 Residual (t.g. only) 0.05 ± 0.40 0.08 ± 0.72 Residual (t.g. + sat) 0.32 ± 0.39 0.69 ± 0.62 a Linear trends in mm yr−1 (with one standard deviation error estimates) are shown for the sea-level time series from the reconstructed tide-gauge data (t.g.) and from joining the altimeter data to the reconstructed data in 1993 (t.g. + sat) and for each term in the sea-level budget for two separate time intervals. Uncertainty estimates for the sea level, shallow thermal expansion and the ground water depletion are from our analysis and all other uncertainty estimates are from the relevant publications as cited in the text. b Bold numbers indicate sum of other rows, as indicated in first column. Table 2. The Earth's Heat Budgeta Component 1972–2008 1993–2008 Shallow ocean (0–700m) 112.6 45.9 Deep ocean (700–3000m) 49.7 20.7 Abyssal ocean (3000m-bottom) 30.7 12.8 Total ocean storage 193.0b 79.4 Glaciers (Latent only) 3.0 1.7 Antarctica (Latent only) 1.4 0.8 Greenland (Latent only) 0.7 0.6 Sea ice 2.5 1.0 Continents 4.7 2.0 Atmosphere 2.0 1.2 Total other storage 14.2 7.3 Total storage 207.2 86.7 Solar + Ozone + well-mixed GHGs 1461.7 709.6 Energy consumption 13.0 6.5 Volcanic (GISS) −207.8 −14.6 Outgoing radiation −343.4 −199.3 Total forcing 923.6 502.2 Total forcing - Total storage 716.4 415.5 a The integrated changes in the heat storage and the radiative forcing are in units of 10+21 J. The total forcing minus the total storage is the amount of energy that must be balanced by the aerosol cooling (or other climate forcing). b Bold numbers indicate sum of other rows, as indicated in first column. Changes in the Terrestrial Storage of Water [6] Previous studies of sea-level change lacked robust estimates of groundwater depletion. We use results from a recent study of Konikow [2011]. Data on the long-term net depletion of groundwater (the net reduction in the volume of water stored in the subsurface saturated zone), computed directly on a volumetric basis, are most widely available in the United States, but generally sparse or unavailable elsewhere. Depletion volumes and rates in the U.S.A. [Konikow, 2011] (Figure 1) were calculated using comprehensive calibrated groundwater simulation models, analytical approaches [e.g., Konikow and Neuzil, 2007], and/or groundwater budget analyses applied to 41 separate aquifer systems or geographic subareas of aquifer systems. [7] Global depletion is calculated where direct estimates are available (Figure 1 and Table S1 in the auxiliary material), and estimated for unmeasured areas outside of the U.S.A. by correlations with global pumpage estimates. In the U.S.A., total groundwater withdrawals in 2000 were 115 km3 [Hutson et al., 2004], and the total depletion that year was 17.7 km3, or 15.4% of withdrawals. The total global withdrawals that year [Margat and Custodio, 2004; Shah et al., 2007] were 800 km3. Assuming that the global depletion represents about the same fraction of global withdrawals as in the U.S.A., then global depletion would be approximately 123 km3. Data from the U.S.A. and several aquifer systems outside the U.S.A. indicate a net depletion in 2000 of 86 km3, so the unaccounted depletion in the rest of the world would be 37 km3 —30% of the total and slightly more than twice that in the U.S.A. The temporal rate of depletion from 1960–2008 was estimated in 46 separate areas by computations from calibrated simulation models or by correlation with pumpage records or water-level records. The average time rate was then applied (assuming a constant value) to the uncounted depletion in the rest of the world, yielding the rates shown in Figure 1. The uncertainty in each approach was estimated, and the weighted average indicated that the uncertainty in the global estimate is ±27%. [8] Konikow's [2011] calculations indicate that during 1961–2008, groundwater was depleted by 727 km3 in the U.S.A. (averaging 15.2 km3 yr−1), equivalent to a global sea-level rise of 2.0 mm in the 48-year period. Global groundwater depletion during 1961–2008 was 3,960 km3, and the rate of depletion averaged 82 ± 22 km3 yr−1 (0.23 ± 0.06 mm yr−1) during the 48-year period. The depletion rate was increasing and was 122 km3 yr−1 (0.34 mm yr−1) during 2000 and 140 km3 yr−1 (0.39 mm yr−1) during 2008, which can account for 12% of the reported recent rate of sea-level rise of 3.2 mm yr−1, estimated from satellite altimeter data [Church and White, 2011]. [9] The method used by Konikow [2011] contrasts with the approach followed by Wada et al. [2010], who estimated depletion indirectly from the difference between global-model based estimates of natural recharge and indirect estimates of withdrawal fluxes. Their approach assumes that depletion equals the difference between natural recharge and (anthropogenic) withdrawals, and may inadequately account for surface water-groundwater interactions, artificial recharge, and the fact that withdrawals are also compensated by reductions in natural discharge and increases in recharge, as described by Theis [1940]. Consequently, the Wada et al. [2010] estimates are almost a factor of two larger than the Konikow estimates used here. [10] The retention of water in man-made reservoirs, with associated seepage into the surrounding soil [Chao et al., 2008], reduced sea level by around 30 mm from 1900 to 2008, mostly after 1950 (slightly over 0.5 mm yr−1; no error bars are given and we assign an uncertainty of 30% based on the difference between the nominal and actual reservoir capacity and that the reservoirs are assumed to be 85% full). Natural exchanges with terrestrial reservoirs average out to near-zero over multi-decadal time periods [Ngo-Duc et al., 2005; Milly and Shmakin, 2002] but can have larger rates over shorter periods. The total terrestrial storage term partially offsets other contributions to sea-level rise (Figure 2a). Ocean Changes [11] For 1961 to 2008, upper (0–700 m) ocean thermosteric sea level and heat content were updated from Domingues et al. [2008]. Ocean temperature profiles (1960–2004) are from the ENACT/ENSEMBLES version 3b (EN3b) data base [Ingleby and Huddleston, 2007] for bottles, Conductivity-Temperature-Depth (CTDs) and eXpendable BathyThermographs (XBTs). To minimise biases in XBT data [Gouretski and Koltermann, 2007], we use a time-dependent (but depth constant) fall-rate correction [Wijffels et al., 2008]. To complement the above EN3b data selection, we use the most recent version of quality-controlled Argo data (2000–2008), corrected for pressure biases [Barker et al., 2011]. [12] Thermal expansion and heat content estimates of the upper 700 m of the global ocean fall by about 7 mm and 40 × 1021 J from 1961 to 1969 and then rise steadily, oscillating (less than 5 mm and 30 × 1021 J) around the long-term trends, contributing almost 18 mm and 124 × 1021 J by 2008, with respect to 1961. Because of a lack of data, we use a linear trend for the deep-ocean contribution. For the depth range 700–3000 m for 1955 to 2003, the trends are 0.07 ± 0.1 mm yr−1 and 1 × 1021 J yr−1 [Levitus et al., 2005; Antonov et al., 2005], but these estimates may be biased low because of inadequate deep-ocean sampling, particularly in the early part of the record. Below 3000 m and since the 1980s, the trends are 0.09 ± 0.06 mm yr−1 and 1 × 1021 J yr−1 [Purkey and Johnson, 2010]. There are few long-term observations to indicate if this is a recent or a longer-term trend. Either way, it would contribute less than 5 mm to sea level since 1961 and we use this trend for the full period. The full-depth thermal expansion is about 25 mm over 1961–2008 (Figure 2a; 0.8 ± 0.1 mm yr−1 since 1972). Glaciers and Ice Caps [13] The contributions from Glaciers and Ice Caps (G&IC) were updated from Cogley [2009] by interpolation of direct (in-situ) and geodetic (based on repeated mapping and leading to larger estimates than previous results) measurements of mass balance for the entire glacierized area outside the Greenland Ice Sheet (GIS) and Antarctica [Cogley, 2009]. The contributions from glaciers around the Antarctic Ice Sheet (AIS), where measurements are few, were obtained by upscaling [Cogley, 2009]. The G&IC contributions, with pentadal resolution from 1951–2008, are interpolated to yearly values. G&IC melting contributed over 30 mm for the full period and increased in the mid 1970s and the late 1990s. Comparatively little energy is required to melt this quantity of glacial ice (4 × 1021 J). Greenland and Antarctica [14] We use an improved estimate of the GIS contribution defined by the ice discharge minus the surface mass balance [Rignot et al., 2008, 2011; van den Broeke et al., 2009]. It is positive from 1961 to the early 1970s, then essentially zero until the 1990s when it becomes positive once more and contributes to sea-level at a rate similar to ocean thermal expansion. [15] The mass balance of the AIS is less well constrained and remains controversial (for example, see Velicogna [2009] versus Zwally and Giovinetto [2011]). We use the estimates of Rignot et al. [2011] which are based on atmospheric climate models constrained by re-analysis data and ice discharge estimates from 1980. The latter are derived from velocities determined from satellite synthetic aperture radar data and ice thickness, mostly from radar echo sounding. The Antarctic contribution increases through the 1990s and is consistent with GRACE gravity data since 2003 [Velicogna, 2009]. Prior to 1980, we assume a value of 0.1 mm yr−1 in 1950 (based on the work by Lemke et al. [2007, p. 365]) and then interpolate to yearly values. The 1950 value is not well determined and could be in the range from zero to perhaps as much as 0.4 mm yr−1; the value chosen here is also consistent with the IPCC Fourth Assessment Report (AR4) estimates. Using either zero or 0.2 mm yr−1 in 1950 would alter the integrated sea-level change from 1950 to the end of the record by less than 2 mm. The AIS contributes slightly more than Greenland over the full period. Melting of the Greenland and Antarctic ice over 1972 to 2008 requires about 2 × 1021 J and melting of sea ice [Levitus et al., 2005] a similar amount. [16] Note that from 2004 to the end of the record, the ice sheets contribute to sea level at a combined rate of 1.3 ± 0.4 mm yr−1, approaching the estimated combined rate from upper ocean thermal expansion and the melting of G&ICs. Other Contributions [17] Changes in atmospheric water content [Trenberth et al., 2005] and the melting of sea ice and ice shelves [Shepherd et al., 2010] make insignificant contributions to sea level change. The additional heat storage in the atmosphere is computed from the mass of the atmosphere (5.3 × 1018 kg) [Gill, 1982] and the heat capacity Cp of 1000 J kg−1 K−1. This amounts to 2.0 × 1021 J since 1972. Huang [2006] estimates the additional heat stored in the ground from 1950 to 2000 as 6.7 × 1021 J and we assume this increase to be linear. The geothermal heat flux of 0.087 W m−2 [1.4 × 1021 J yr−1; Pollack et al., 1993] is not included in the budget as this remains relatively constant over time and is not a component of the perturbations to the energy budget since preindustrial times that we focus on here. For comparison, we do include the anthropogenic consumption of energy using a 1980–2006 estimate from the US Energy Information Administration (http://www.eia.gov/pub/international/iealf/tablee1.xls), linearly extrapolated to 1972 and 2008. Although increasing rapidly, this term (0.4 × 1021 J yr−1, less than 0.01 W m−2) is more than 100 times smaller than the forcing from greenhouse gases and more than an order of magnitude less than the energy required for the observed heat storage increases. 3. Observed Sea-Level Rise and the Sea-Level Budget [18] We use an estimate of global averaged sea level from Church and White [2011]. From 1993 to 2009, they estimated global-averaged sea level with satellite altimeter data using standard techniques. For the period from 1961, they ‘reconstructed’ sea level using coastal and island tide gauges and empirical orthogonal functions of sea level calculated with satellite altimeter data. Fiedler and Conrad [2010] indicate that there can be an important impact on tide-gauge data from terrestrial loading and gravitational changes resulting from dam storage. However, for the large number of tide gauges used here (over 200), we estimate the impact of this correction (which we include) on global mean sea level is only about 0.05 mm yr−1. Similar corrections for changes in the cryosphere are smaller, at least until recently. The satellite and the in situ sea-level time series have also been corrected for glacial isostatic adjustment [Davis and Mitrovica, 1996; M. Tamisiea, personal communication, 2004] so that the sea-level estimates are a measure of changes in ocean volume. [19] The reconstructed sea level falls more than 10 mm from 1961 to 1966 and then rises to the end of the record with small oscillations of about 5 mm partly related to volcanic eruptions [Domingues et al., 2008]. The greater rate of rise after 1993 in the tide-gauge data is confirmed by satellite altimeter observations. The reconstructed sea level mostly agrees with the more accurate satellite data within uncertainty limits and over the satellite record the total rise is almost identical but the linear trend for the tide gauges over the 1993–2008 period is smaller (by about 20%). [20] The sum of the contributions agrees well with the observed rate of sea-level rise, and certainly within error bars (Figure 2b and Table 1) from 1972. Error estimates (Table 1) come from weighted least-squares estimates of the trends for observed sea level and upper-ocean thermal expansion, and for the remaining time series from uncertainty estimates provided by the originators of the data sets (where given). Thermal expansion contributes about 40% of the observed rise since 1972, the proportion decreasing as the ice contributions increase. The sum of the glacier and thermosteric contributions explains approximately 75% of the observed rise since 1972. It is also clear that both the Greenland and Antarctic ice sheets are an important component in closing the sea-level budget, particularly since 1993, as found for example by Cazenave and Llovel [2010]. The larger differences between the sum of contributions and the observed sea-level change occur in the 1960s when sea level, thermal expansion and terrestrial storage estimates fall but the poorly constrained cryospheric contributions are positive. This poorer agreement may be a result of the poorer ocean (the number of ocean observations increased significantly in the early 1970s) and cryospheric data coverage in the 1960s. As a result, it is likely we cannot adequately resolve the response to and the recovery from the 1963 volcanic eruption of Mount Agung. The response of the climate system to volcanic eruptions in the 1980s and 1990s is likely to be better observed but perhaps still incompletely resolved. [21] From 1993, the sum of contributions is not significantly different to the observed rise but the linear trend in the altimeter observations is larger. A larger terrestrial storage component, as inferred by Wada et al. [2010], would reduce these differences. Note that similar sea-level budget closures for 1955 to 2005 [e.g., Moore et al., 2011] and 1993 to 2007 [e.g., Cazenave and Llovel, 2010] have recently been obtained using slightly different contributions from those estimated here. A significant issue remains the identification of systematic errors in the observational data base. For example, the IPCC AR4 sea-level budget for 1993–2003 [Bindoff et al., 2007] has an upper-ocean thermosteric sea-level rise of 1.6 mm yr−1. This is about twice the current estimates (Table 1), inconsistent with ocean warming estimated using XBT data corrected for fall-rate errors [Wijffels et al., 2008]. 4. Implications for the Earth's Energy Budget [22] We follow the approach of Murphy et al. [2009] to examine the Earth's energy balance. We use updates of radiative forcing (F) from well-mixed greenhouse gases and fluctuations in solar radiation [Gregory and Forster, 2008] (with the stratospheric ozone scaled as by Murphy et al. [2009]) and volcanic aerosols from the Goddard Institute for Space Studies (GISS) [Sato et al., 1993], updated at http://data.giss.nasa.gov/modelforce/RadF.txt, and Gregory and Forster [2008]. Integrated from 1961 through to the end of 2008, well-mixed greenhouse gases, ozone and solar forcing would result in the storage of 1,735 × 1021 J of additional heat in the Earth's system if it were all retained (Figure 3). However, as discussed in section 2, only 250 × 1021 J (234 from 1961–2008, 207 from 1972–2008) of the forcing is stored in the Earth system (Figure 3a and Table 2), over 90% in the ocean with minor storage in the atmosphere and the solid earth and consumption in the melting of ice. Negative radiative forcing from volcanic aerosols offset about 320 × 1021 J (317 for GISS, 398 for Gregory and Forster [2008]) of the greenhouse gas forcing. The remaining portion (1200 × 1021 J; 1190 for the whole period, 1060 for 1972–2008) must be balanced by the sum of reflection of solar radiation by tropospheric aerosols and increased radiation from a warming Earth. [23] The rate of storage of heat by the climate system equals the net heat flux into the climate system N = F − λΔT, where F is the radiative forcing, ΔT is the global averaged temperature change and the climate feedback parameter λ

As part of the Coupled Model Intercomparison Project, integrations with a common design have been undertaken with eleven different climate models to compare the response of the Atlantic thermohaline circulation (THC) to time‐dependent climate change caused by increasing atmospheric CO 2 concentration. Over 140 years, during which the CO 2 concentration quadruples, the circulation strength declines gradually in all models, by between 10 and 50%. No model shows a rapid or complete collapse, despite the fairly rapid increase and high final concentration of CO 2 . The models having the strongest overturning in the control climate tend to show the largest THC reductions. In all models, the THC weakening is caused more by changes in surface heat flux than by changes in surface water flux. No model shows a cooling anywhere, because the greenhouse warming is dominant.

A warming signal has penetrated into the world's oceans over the past 40 years. The signal is complex, with a vertical structure that varies widely by ocean; it cannot be explained by natural internal climate variability or solar and volcanic forcing, but is well simulated by two anthropogenically forced climate models. We conclude that it is of human origin, a conclusion robust to observational sampling and model differences. Changes in advection combine with surface forcing to give the overall warming pattern. The implications of this study suggest that society needs to seriously consider model predictions of future climate change.