Global Mean Surface Temperature Response Research Articles

National contributions to climate change due to historical emissions of carbon dioxide, methane, and nitrous oxide since 1850

Anthropogenic emissions of carbon dioxide (CO2), methane (CH4) and nitrous oxide (N2O) have made significant contributions to global warming since the pre-industrial period and are therefore targeted in international climate policy. There is substantial interest in tracking and apportioning national contributions to climate change and informing equitable commitments to decarbonisation. Here, we introduce a new dataset of national contributions to global warming caused by historical emissions of carbon dioxide, methane, and nitrous oxide during the years 1851–2021, which are consistent with the latest findings of the IPCC. We calculate the global mean surface temperature response to historical emissions of the three gases, including recent refinements which account for the short atmospheric lifetime of CH4. We report national contributions to global warming resulting from emissions of each gas, including a disaggregation to fossil and land use sectors. This dataset will be updated annually as national emissions datasets are updated.

Effect of Resolving Ocean Eddies on the Transient Response of Global Mean Surface Temperature to Abrupt 4xCO2 Forcing

AbstractThe magnitude of global mean surface temperature (GMST) response to increasing atmospheric CO2 concentrations is affected by the efficiency of ocean heat uptake, which in turn can be affected by oceanic mesoscale eddies. Using the Max Planck Institute ‐ Earth System Model (MPI‐ESM1.2), we find that resolving eddies leads to a −0.1°C cooler response of GMST to an abrupt CO2 quadrupling, which is related to a larger rate of heat uptake by an eddying ocean. This is consistent with changes in the energy budget of the whole climate system induced by increasing ocean resolution under the same radiative forcing and climate feedback. As a fast response, heat is taken up by the deep ocean, independent of resolution. The change in deep ocean heat uptake due to resolved eddies is an amplification in the magnitude of the responses of all heat processes, including eddy heat advection, mean heat advection, and diffusive processes.

DeepMIP: model intercomparison of early Eocene climatic optimum (EECO) large-scale climate features and comparison with proxy data

Abstract. We present results from an ensemble of eight climate models, each of which has carried out simulations of the early Eocene climate optimum (EECO, ∼ 50 million years ago). These simulations have been carried out in the framework of the Deep-Time Model Intercomparison Project (DeepMIP; http://www.deepmip.org, last access: 10 January 2021); thus, all models have been configured with the same paleogeographic and vegetation boundary conditions. The results indicate that these non-CO2 boundary conditions contribute between 3 and 5 ∘C to Eocene warmth. Compared with results from previous studies, the DeepMIP simulations generally show a reduced spread of the global mean surface temperature response across the ensemble for a given atmospheric CO2 concentration as well as an increased climate sensitivity on average. An energy balance analysis of the model ensemble indicates that global mean warming in the Eocene compared with the preindustrial period mostly arises from decreases in emissivity due to the elevated CO2 concentration (and associated water vapour and long-wave cloud feedbacks), whereas the reduction in the Eocene in terms of the meridional temperature gradient is primarily due to emissivity and albedo changes owing to the non-CO2 boundary conditions (i.e. the removal of the Antarctic ice sheet and changes in vegetation). Three of the models (the Community Earth System Model, CESM; the Geophysical Fluid Dynamics Laboratory, GFDL, model; and the Norwegian Earth System Model, NorESM) show results that are consistent with the proxies in terms of the global mean temperature, meridional SST gradient, and CO2, without prescribing changes to model parameters. In addition, many of the models agree well with the first-order spatial patterns in the SST proxies. However, at a more regional scale, the models lack skill. In particular, the modelled anomalies are substantially lower than those indicated by the proxies in the southwest Pacific; here, modelled continental surface air temperature anomalies are more consistent with surface air temperature proxies, implying a possible inconsistency between marine and terrestrial temperatures in either the proxies or models in this region. Our aim is that the documentation of the large-scale features and model–data comparison presented herein will pave the way to further studies that explore aspects of the model simulations in more detail, for example the ocean circulation, hydrological cycle, and modes of variability, and encourage sensitivity studies to aspects such as paleogeography, orbital configuration, and aerosols.

Estimating the Effect of Radiative Feedback Uncertainties on Climate Response to Changes in the Concentration of Stratospheric Aerosols

Using the two-box energy balance model (EBM), we explore the climate system response to radiative forcing generated by variations in the concentrations of stratospheric aerosols and estimate the effect of uncertainties in radiative feedbacks on changes in global mean surface temperature anomaly used as an indicator of the response of the climate system to external radiative perturbations. Radiative forcing generated by stratospheric sulfate aerosols from the second-largest volcanic eruption in the 20th century, the Mount Pinatubo eruption in June 1991, was chosen for this research. The global mean surface temperature response to a specified change in radiative forcing is estimated as a convolution of the derived impulse response function corresponding to EBM with a function that describes the temporal change in radiative forcing. The influence of radiative feedback uncertainties on changes in the global mean surface temperature is estimated using several “versions” of the EBM. The parameters for different “versions” were identified by applying a specific procedure for calibrating the two-box EBM parameters using the results of climate change simulations conducted with coupled atmosphere–ocean general circulation models from the Coupled Model Intercomparison Project phase 5 (CMIP5). Changes in the global mean surface temperature caused by stratospheric aerosol forcing are found to be highly sensitive not only to radiative feedbacks but also to climate system inertia defined by the effective heat capacity of the atmosphere–land–ocean mixed layer system, as well as to deep-ocean heat uptake. The results obtained have direct implications for a better understanding of how uncertainties in climate feedbacks, climate system inertia and deep-ocean heat uptake affect climate change modelling.

Statistical estimation of global surface temperature response to forcing under the assumption of temporal scaling

Abstract. Reliable quantification of the global mean surface temperature (GMST) response to radiative forcing is essential for assessing the risk of dangerous anthropogenic climate change. We present the statistical foundations for an observation-based approach using a stochastic linear response model that is consistent with the long-range temporal dependence observed in global temperature variability. We have incorporated the model in a latent Gaussian modeling framework, which allows for the use of integrated nested Laplace approximations (INLAs) to perform full Bayesian analysis. As examples of applications, we estimate the GMST response to forcing from historical data and compute temperature trajectories under the Representative Concentration Pathways (RCPs) for future greenhouse gas forcing. For historic runs in the Model Intercomparison Project Phase 5 (CMIP5) ensemble, we estimate response functions and demonstrate that one can infer the transient climate response (TCR) from the instrumental temperature record. We illustrate the effect of long-range dependence by comparing the results with those obtained from one-box and two-box energy balance models. The software developed to perform the given analyses is publicly available as the R package INLA.climate.

Local and remote mean and extreme temperature response to regional aerosol emissions reductions

Abstract. The climatic implications of regional aerosol and precursor emissions reductions implemented to protect human health are poorly understood. We investigate the mean and extreme temperature response to regional changes in aerosol emissions using three coupled chemistry–climate models: NOAA GFDL CM3, NCAR CESM1, and NASA GISS-E2. Our approach contrasts a long present-day control simulation from each model (up to 400 years with perpetual year 2000 or 2005 emissions) with 14 individual aerosol emissions perturbation simulations (160–240 years each). We perturb emissions of sulfur dioxide (SO2) and/or carbonaceous aerosol within six world regions and assess the statistical significance of mean and extreme temperature responses relative to internal variability determined by the control simulation and across the models. In all models, the global mean surface temperature response (perturbation minus control) to SO2 and/or carbonaceous aerosol is mostly positive (warming) and statistically significant and ranges from +0.17 K (Europe SO2) to −0.06 K (US BC). The warming response to SO2 reductions is strongest in the US and Europe perturbation simulations, both globally and regionally, with Arctic warming up to 1 K due to a removal of European anthropogenic SO2 emissions alone; however, even emissions from regions remote to the Arctic, such as SO2 from India, significantly warm the Arctic by up to 0.5 K. Arctic warming is the most robust response across each model and several aerosol emissions perturbations. The temperature response in the Northern Hemisphere midlatitudes is most sensitive to emissions perturbations within that region. In the tropics, however, the temperature response to emissions perturbations is roughly the same in magnitude as emissions perturbations either within or outside of the tropics. We find that climate sensitivity to regional aerosol perturbations ranges from 0.5 to 1.0 K (W m−2)−1 depending on the region and aerosol composition and is larger than the climate sensitivity to a doubling of CO2 in two of three models. We update previous estimates of regional temperature potential (RTP), a metric for estimating the regional temperature responses to a regional emissions perturbation that can facilitate assessment of climate impacts with integrated assessment models without requiring computationally demanding coupled climate model simulations. These calculations indicate a robust regional response to aerosol forcing within the Northern Hemisphere midlatitudes, regardless of where the aerosol forcing is located longitudinally. We show that regional aerosol perturbations can significantly increase extreme temperatures on the regional scale. Except in the Arctic in the summer, extreme temperature responses largely mirror mean temperature responses to regional aerosol perturbations through a shift of the temperature distributions and are mostly dominated by local rather than remote aerosol forcing.

Using Late Pleistocene sea surface temperature reconstructions to constrain future greenhouse warming

Future greenhouse warming projections conducted with coupled climate models still exhibit a substantial spread in response to a given anthropogenic greenhouse gas concentration scenario. In order to constrain this spread and to provide robust warming projections, our understanding of Earth's global-mean surface temperature response to radiative forcing (referred to as climate sensitivity) needs to be further refined. Here we estimate an averaged glacial/interglacial climate sensitivity using 25 transient Earth system model simulations of the Last Glacial Cycle and a global-mean sea surface temperature (SST) reconstruction derived from 64 globally-distributed paleo-proxies of SST. Our results document that Earth's averaged Late Pleistocene equilibrium climate sensitivity is in the order of ∼4.2 K per CO2 doubling. Using the Representative Concentration Pathway 8.5 for future greenhouse radiative forcing, this value translates into a global-mean surface warming of ∼5.0 K by the year 2100 relative to pre-industrial levels. This estimate is in excellent agreement with the ensemble-mean projection of climate simulations conducted as part of the Coupled Model Intercomparison Project Phase 5 (CMIP5). Our uncertainty analysis reveals further that the lack of robust reconstructions of glacial aerosol forcing is a key contributor to the overall uncertainty of paleo-based estimates of climate sensitivity.

The response of the ozone layer to quadrupled CO2 concentrations.

An accurate quantification of the stratospheric ozone feedback in climate change simulations requires knowledge of the ozone response to increased greenhouse gases. Here, we present an analysis of the ozone layer response to an abrupt quadrupling of CO2 concentrations in four chemistry-climate models. We show that increased CO2 levels lead to a decrease in ozone concentrations in the tropical lower stratosphere, and an increase over the high latitudes and throughout the upper stratosphere. This pattern is robust across all models examined here, although important inter-model differences in the magnitude of the response are found. As a result of the cancellation between upper and lower stratospheric ozone, the total column ozone response in the tropics is small, and appears to be model dependent. A substantial portion of the spread in the tropical column ozone is tied to inter-model spread in upwelling. The high latitude ozone response is strongly seasonally dependent, and shows increases peaking in late-winter and spring of each hemisphere, with prominent longitudinal asymmetries. The range of ozone responses to CO2 reported in this paper has the potential to induce significant radiative and dynamical effects on the simulated climate. Hence, these results highlight the need of using an ozone dataset consistent with CO2 forcing in models involved in climate sensitivity studies.

Global Mean Surface Temperature Response to Large‐Scale Patterns of Variability in Observations and CMIP5

AbstractGlobal mean surface temperature (GMST) fluctuates over decadal to multidecadal time scales. Patterns of internal variability are partly responsible, but the relationships can be conflated by anthropogenically forced signals. Here we adopt a physically based method of separating internal variability from forced responses to examine how trends in large‐scale patterns, specifically the Interdecadal Pacific Oscillation (IPO) and Atlantic Multidecadal Variability (AMV), influence GMST. After removing the forced responses, observed variability of GMST is close to the central estimates of Coupled Model Intercomparison Project phase 5 simulations, but models tend to underestimate IPO variability at time scales >10 years, and AMV at time scales >20 years. Correlations between GMST trends and these patterns are also underrepresented, most strongly at 10‐ and 35‐year time scales, for IPO and AMV, respectively. Strikingly, models that simulate stronger variability of IPO and AMV also exhibit stronger relationships between these patterns and GMST, predominately at the 10‐ and 35‐year time scales, respectively.

Nonlocal Effects Dominate the Global Mean Surface Temperature Response to the Biogeophysical Effects of Deforestation

AbstractDeforestation influences surface temperature locally (“local effects”), but also at neighboring or remote regions (“nonlocal effects”). Observations indicate that local effects induce a warming in most locations, while many climate models show a global mean cooling when simulating global deforestation. We show that a nonlocal cooling in models, which is excluded from observations, may strongly contribute to these conflicting results. For the MPI‐ESM, the globally averaged nonlocal cooling exceeds the globally averaged local warming by a factor of three, for global deforestation but also for realistic areal extents and spatial distributions of deforestation. Furthermore, the globally averaged nonlocal effects dominate the local effects in realistic scenarios across a range of climate models. We conclude that observations alone are not sufficient to capture the full biogeophysical effects, and climate models are needed to better understand and quantify the full effects of deforestation before they are considered in strategies for climate mitigation.

Low clouds link equilibrium climate sensitivity to hydrological sensitivity

Equilibrium climate sensitivity (ECS) and hydrological sensitivity describe the global mean surface temperature and precipitation responses to a doubling of atmospheric CO2. Despite their connection via the Earth’s energy budget, the physical linkage between these two metrics remains controversial. Here, using a global climate model with a perturbed mean hydrological cycle, we show that ECS and hydrological sensitivity per unit warming are anti-correlated owing to the low-cloud response to surface warming. When the amount of low clouds decreases, ECS is enhanced through reductions in the reflection of shortwave radiation. In contrast, hydrological sensitivity is suppressed through weakening of atmospheric longwave cooling, necessitating weakened condensational heating by precipitation. These compensating cloud effects are also robustly found in a multi-model ensemble, and further constrained using satellite observations. Our estimates, combined with an existing constraint to clear-sky shortwave absorption, suggest that hydrological sensitivity could be lower by 30% than raw estimates from global climate models. The connections between global mean temperature and precipitation responses to CO2 doubling (equilibrium climate and hydrological sensitivity) are driven through low-cloud responses to surface warming, according to MIROC5 perturbation experiments.

Quantification of temperature response to CO2 forcing in atmosphere–ocean general circulation models

The present study establishes a general formulation to represent the behavior and variation of an ensemble of complex climate models in terms of the global mean surface temperature response to atmospheric CO2 increase. The response parameters of this formulation provide a set of metrics that extends the conventional concept of climate sensitivity and quantifies transient temperature changes with sufficient simplicity and transparency to serve studies on climate change mitigation. Two commonly used metrics for transient and equilibrium climate sensitivity are analytically derived from the formulation, such that conventional estimates of equilibrium climate sensitivity based on standard numerical experiments for quadrupling CO2 increase are properly scaled down to the reference level of doubling CO2. The characteristics and variations of a specific ensemble of complex climate models can be simulated with a statistical model built using the principal components (PCs) of the response parameters. This approach is applied to the probabilistic assessment of temperature changes as well as to the diagnosis of the base ensemble. In current complex climate models, the ratio of transient-to-equilibrium sensitivity decreases with an increase of equilibrium sensitivity, as identified in variations associated with two specific PCs that characterize coherence between transient temperature response and properties of heat uptake by the ocean.

Uncertainty in crossing time of 2°C warming threshold over China

The 2°C warming target has been used widely in global and regional climate change research. Previous studies have shown large uncertainties in the time when surface air temperature (SAT) change over China will reach 2°C relative to the pre-industrial era. To understand the uncertainties, we analyzed the projected SAT in the twenty-first century using 40 state-of-the-art climate models under two Representative Concentration Pathways (RCP4.5 and RCP8.5) from the Coupled Model Intercomparison Project Phase 5. The 2°C threshold-crossing time (TCT) of SAT averaged across China was around 2033 and 2029 for RCP4.5 and RCP8.5, respectively. Considering a ±1σ range of intermodel SAT change, the upper and lower bounds of the 2°C TCT could differ by about 25 years or even more. Uncertainty in the projected SAT and the warming rate around the TCT are the two main factors responsible for the TCT uncertainty. The former is determined by the climate sensitivity represented by the global mean surface temperature response. About 45 % of the intermodel variance of the projected 2°C TCT for averaged SAT over China can be explained by climate sensitivity across the models, which is contributed mainly by central and southern China. In a climate more sensitive to CO2 forcing, stronger greenhouse effect, less stratus cloud over the East Asian monsoon region, and less snow cover on the Tibetan Plateau result in increased downward longwave radiation, increased shortwave radiation, and decreased shortwave radiation reflected by the surface, respectively, all of which may advance the TCT.

A review of progress towards understanding the transient global mean surface temperature response to radiative perturbation

The correct understanding of the transient response to external radiative perturbation is important for the interpretation of observed climate change, the prediction of near-future climate change, and committed warming under climate stabilization scenarios, as well as the estimation of equilibrium climate sensitivity based on observation data. It has been known for some time that the radiative damping rate per unit of global mean surface temperature increase varies with time, and this inconstancy affects the transient response. Knowledge of the equilibrium response alone is insufficient, but understanding the transient response of the global mean surface temperature has made rapid progress. The recent progress accompanies the relatively new concept of the efficacies of ocean heat uptake and forcing. The ocean heat uptake efficacy associates the temperature response induced by ocean heat uptake with equilibrium temperature response, and the efficacy of forcing compares the temperature response caused by non-CO2 forcing with that by CO2 forcing. In this review article, recent studies on these efficacies are discussed, starting from the classical global feedback framework and basis of the transient response. An attempt is made to structure different studies that emphasize different aspects of the transient response and to stress the relevance of those individual studies. The implications on the definition and computation of forcing and on the estimation of the equilibrium response in climate models are also discussed. Along with these discussions, examples are provided with MIROC climate model multi-millennial simulations.

Global warming projections derived from an observation-based minimal model

Abstract. A simple conceptual model for the global mean surface temperature (GMST) response to CO2 emissions is presented and analysed. It consists of linear long-memory models for the GMST anomaly response ΔT to radiative forcing and the atmospheric CO2-concentration response ΔC to emission rate. The responses are connected by the standard logarithmic relation between CO2 concentration and its radiative forcing. The model depends on two sensitivity parameters, αT and αC, and two "inertia parameters," the memory exponents βT and βC. Based on observation data, and constrained by results from the Climate Model Intercomparison Project Phase 5 (CMIP5), the likely values and range of these parameters are estimated, and projections of future warming for the parameters in this range are computed for various idealised, but instructive, emission scenarios. It is concluded that delays in the initiation of an effective global emission reduction regime is the single most important factor that influences the magnitude of global warming over the next 2 centuries. The most important aspect of this study is the simplicity and transparency of the conceptual model, which makes it a useful tool for communicating the issue to non-climatologists, students, policy makers, and the general public.

Interactions between perturbations to different Earth system components simulated by a fully-coupled climate model

We examine the global mean surface temperature and carbon cycle responses to the A1B emissions scenario for a new 57 member perturbed-parameter ensemble of simulations generated using the fully coupled atmosphere-ocean-carbon cycle climate model HadCM3C. The model variants feature simultaneous perturbation to parameters that control atmosphere, ocean, land carbon cycle and sulphur cycle processes in this Earth system model, and is the first experiment of its kind. The experimental design, based on four earlier ensembles with parameters varied within each individual Earth system component, allows the effects of interactions between uncertainties in the different components to be explored. A large spread in response is obtained, with atmospheric CO2 at the end of the twenty-first century ranging from 615 to 1,100 ppm. On average though, the mean effect of the parameter perturbations is to significantly reduce the amount of atmospheric CO2 compared to that seen in the standard HadCM3C model. Global temperature change for 2090–2099 relative to the pre-industrial period ranges from 2.2 to 7.5 °C, with large temperature responses occurring when atmospheric model versions with high climate sensitivities are combined with carbon cycle components that emit large amounts of CO2 to the atmosphere under warming. A simple climate model, tuned to reproduce the responses of the separate Earth system component ensembles, is used to demonstrate that interactions between uncertainties in the different components play a significant role in determining the spread of responses in global mean surface temperature. This ensemble explores a wide range of interactions and response, and therefore provides a useful resource for the provision of regional climate projections and associated uncertainties.

Deep-ocean heat uptake and equilibrium climate response

We integrate the coupled climate model ECHAM5/MPIOM to equilibrium under atmospheric CO2 quadrupling. The equilibrium global-mean surface-temperature change is 10.8 K. The surface equilibrates within about 1,200 years, the deep ocean within 5,000 years. The impact of the deep ocean on the equilibrium surface-temperature response is illustrated by the difference between ECHAM5/MPIOM and ECHAM5 coupled with slab ocean model (ECHAM5/SOM). The equilibrium global-mean surface temperature response is 11.1 K in ECHAM5/SOM and is thus 0.3 K higher than in ECHAM5/MPIOM. ECHAM5/MPIOM shows less warming over the northern-hemisphere mid and high latitudes, but larger warming over the tropical ocean and especially over the southern-hemisphere high latitudes. ECHAM5/MPIOM shows similar polar amplification in both the Arctic and the Antarctic, in contrast to ECHAM5/SOM, which shows stronger polar amplification in the northern hemisphere. The southern polar warming in ECHAM5/MPIOM is greatly delayed by Antarctic deep-ocean warming due to convective and isopycnal mixing. The equilibrium ocean temperature warming under CO2 quadrupling is around 8.0 K and is near-uniform with depth. The global-mean steric sea-level rise is 5.8 m in equilibrium; of this, 2.3 m are due to the deep-ocean warming after the surface temperature has almost equilibrated. This result suggests that the surface temperature change is a poor predictor for steric sea-level change in the long term. The effective climate response method described in Gregory et al. (2004) is evaluated with our simulation, which shows that their method to estimate the equilibrium climate response is accurate to within 10 %.

Global temperature response to radiative forcing: Solar cycle versus volcanic eruptions

I show that the peak‐to‐peak amplitude of the global mean surface temperature response to the 11‐year cyclic total irradiance forcing is an order of magnitude less than the amplitude of a cyclic component roughly in phase with the solar forcing which has been observed in the temperature record in the period 1959–2004. If this cyclic temperature component were a response to the solar forcing, it would imply the existence of strong amplifying feedbacks which operate exclusively for solar forcing, such as top‐down mechanisms responding to the large variability in the ultraviolet part of the solar spectrum. I demonstrate, however, that the apparent cyclic component in the temperature record is dominated by the response to five major volcanic eruptions some of which incidentally took place a few years before solar minimum in four consecutive solar cycles, and hence that the correlation with the solar cycle is coincidental. A temperature rise of approximately 0.15 K over the 20th century ascribed to an increasing trend in solar forcing is more than offset by a cooling trend of about 0.3 K due to stratospheric aerosols from volcanic eruptions.

The Shape of Things to Come: Why Is Climate Change So Predictable?

Abstract The framework of feedback analysis is used to explore the controls on the shape of the probability distribution of global mean surface temperature response to climate forcing. It is shown that ocean heat uptake, which delays and damps the temperature rise, can be represented as a transient negative feedback. This transient negative feedback causes the transient climate change to have a narrower probability distribution than that of the equilibrium climate response (the climate sensitivity). In this sense, climate change is much more predictable than climate sensitivity. The width of the distribution grows gradually over time, a consequence of which is that the larger the climate change being contemplated, the greater the uncertainty is about when that change will be realized. Another consequence of this slow growth is that further efforts to constrain climate sensitivity will be of very limited value for climate projections on societally relevant time scales. Finally, it is demonstrated that the effect on climate predictability of reducing uncertainty in the atmospheric feedbacks is greater than the effect of reducing uncertainty in ocean feedbacks by the same proportion. However, at least at the global scale, the total impact of uncertainty in climate feedbacks is dwarfed by the impact of uncertainty in climate forcing, which in turn is contingent on choices made about future anthropogenic emissions.

Long run surface temperature dynamics of an A-OGCM: the HadCM3 4×CO2 forcing experiment revisited

The global mean surface temperature (GMST) response of HadCM3 to a 1,000 year 4×CO2 forcing is analysed using a transfer function methodology. We identify a third order transfer function as being an appropriate characterisation of the dynamic relationship between the radiative forcing input and GMST output of this Atmosphere-Ocean General Circulation Model (A-OGCM). From this transfer function the equilibrium climate sensitivity is estimated as 4.62 (3.92–11.88) K which is significantly higher than previously estimated for HadCM3. The response is also characterised by time constants of 4.5 (3.2–6.4), 140 (78–191) and 1,476 (564–11,737) years. The fact that the longest time constant element is significantly longer than the 1,000 year simulation run makes estimation of this element of the response problematic, highlighting the need for significantly longer model runs to express A-OGCM behaviour fully. The transfer function is interpreted in relation to a three box global energy balance model. It was found that this interpretation gave rise to three fractions of ocean heat capacity with effective depths of 63.0 (46.7–85.4), 1291.7 (787.3–2,955.3) and 2,358.0 (661.3–17,283.8) meters of seawater, associated with three discrete time constants of 4.6 (3.2–6.5), 107.7 (68.9–144.3) and 537.1 (196.2–1,243.1) years. Given this accounts for approximately 94% of the ocean heat capacity in HadCM3, it appears HadCM3 could be significantly more well mixed than previously thought when viewed on the millennial timescale.