H. Damon Matthews

The UVic earth system climate model: Model description, climatology, and applications to past, present and future climates

Abstract A new earth system climate model of intermediate complexity has been developed and its climatology compared to observations. The UVic Earth System Climate Model consists of a three‐dimensional ocean general circulation model coupled to a thermodynamic/dynamic sea‐ice model, an energy‐moisture balance atmospheric model with dynamical feedbacks, and a thermomechanical land‐ice model. In order to keep the model computationally efficient a reduced complexity atmosphere model is used. Atmospheric heat and freshwater transports are parametrized through Fickian diffusion, and precipitation is assumed to occur when the relative humidity is greater than 85%. Moisture transport can also be accomplished through advection if desired. Precipitation over land is assumed to return instantaneously to the ocean via one of 33 observed river drainage basins. Ice and snow albedo feedbacks are included in the coupled model by locally increasing the prescribed latitudinal profile of the planetary albedo. The atmospheric model includes a parametrization of water vapour/planetary longwave feedbacks, although the radiative forcing associated with changes in atmospheric CO2 is prescribed as a modification of the planetary longwave radiative flux. A specified lapse rate is used to reduce the surface temperature over land where there is topography. The model uses prescribed present‐day winds in its climatology, although a dynamical wind feedback is included which exploits a latitudinally‐varying empirical relationship between atmospheric surface temperature and density. The ocean component of the coupled model is based on the Geophysical Fluid Dynamics Laboratory (GFDL) Modular Ocean Model 2.2, with a global resolution of 3.6° (zonal) by 1.8° (meridional) and 19 vertical levels, and includes an option for brine‐rejection parametrization. The sea‐ice component incorporates an elastic‐viscous‐plastic rheology to represent sea‐ice dynamics and various options for the representation of sea‐ice thermodynamics and thickness distribution. The systematic comparison of the coupled model with observations reveals good agreement, especially when moisture transport is accomplished through advection. Global warming simulations conducted using the model to explore the role of moisture advection reveal a climate sensitivity of 3.0°C for a doubling of CO2, in line with other more comprehensive coupled models. Moisture advection, together with the wind feedback, leads to a transient simulation in which the meridional overturning in the North Atlantic initially weakens, but is eventually re‐established to its initial strength once the radiative forcing is held fixed, as found in many coupled atmosphere General Circulation Models (GCMs). This is in contrast to experiments in which moisture transport is accomplished through diffusion whereby the overturning is reestablished to a strength that is greater than its initial condition. When applied to the climate of the Last Glacial Maximum (LGM), the model obtains tropical cooling (30°N‐30°S), relative to the present, of about 2.1°C over the ocean and 3.6°C over the land. These are generally cooler than CLIMAP estimates, but not as cool as some other reconstructions. This moderate cooling is consistent with alkenone reconstructions and a low to medium climate sensitivity to perturbations in radiative forcing. An amplification of the cooling occurs in the North Atlantic due to the weakening of North Atlantic Deep Water formation. Concurrent with this weakening is a shallowing of, and a more northward penetration of, Antarctic Bottom Water. Climate models are usually evaluated by spinning them up under perpetual present‐day forcing and comparing the model results with present‐day observations. Implicit in this approach is the assumption that the present‐day observations are in equilibrium with the present‐day radiative forcing. The comparison of a long transient integration (starting at 6 KBP), forced by changing radiative forcing (solar, CO2, orbital), with an equilibrium integration reveals substantial differences. Relative to the climatology from the present‐day equilibrium integration, the global mean surface air and sea surface temperatures (SSTs) are 0.74°C and 0.55°C colder, respectively. Deep ocean temperatures are substantially cooler and southern hemisphere sea‐ice cover is 22% greater, although the North Atlantic conveyor remains remarkably stable in all cases. The differences are due to the long timescale memory of the deep ocean to climatic conditions which prevailed throughout the late Holocene. It is also demonstrated that a global warming simulation that starts from an equilibrium present‐day climate (cold start) underestimates the global temperature increase at 2100 by 13% when compared to a transient simulation, under historical solar, CO2 and orbital forcing, that is also extended out to 2100. This is larger (13% compared to 9.8%) than the difference from an analogous transient experiment which does not include historical changes in solar forcing. These results suggest that those groups that do not account for solar forcing changes over the twentieth century may slightly underestimate (∼3% in our model) the projected warming by the year 2100.

The proportionality of global warming to cumulative carbon emissions

The global temperature response to increasing atmospheric CO2 is often quantified by metrics such as equilibrium climate sensitivity and transient climate response. These approaches, however, do not account for carbon cycle feedbacks and therefore do not fully represent the net response of the Earth system to anthropogenic CO2 emissions. Climate–carbon modelling experiments have shown that: (1) the warming per unit CO2 emitted does not depend on the background CO2 concentration; (2) the total allowable emissions for climate stabilization do not depend on the timing of those emissions; and (3) the temperature response to a pulse of CO2 is approximately constant on timescales of decades to centuries. Here we generalize these results and show that the carbon–climate response (CCR), defined as the ratio of temperature change to cumulative carbon emissions, is approximately independent of both the atmospheric CO2 concentration and its rate of change on these timescales. From observational constraints, we estimate CCR to be in the range 1.0–2.1 °C per trillion tonnes of carbon (Tt C) emitted (5th to 95th percentiles), consistent with twenty-first-century CCR values simulated by climate–carbon models. Uncertainty in land-use CO2 emissions and aerosol forcing, however, means that higher observationally constrained values cannot be excluded. The CCR, when evaluated from climate–carbon models under idealized conditions, represents a simple yet robust metric for comparing models, which aggregates both climate feedbacks and carbon cycle feedbacks. CCR is also likely to be a useful concept for climate change mitigation and policy; by combining the uncertainties associated with climate sensitivity, carbon sinks and climate–carbon feedbacks into a single quantity, the CCR allows CO2-induced global mean temperature change to be inferred directly from cumulative carbon emissions.

Future changes in climate, ocean circulation, ecosystems, and biogeochemical cycling simulated for a business‐as‐usual CO2 emission scenario until year 4000 AD

A new model of global climate, ocean circulation, ecosystems, and biogeochemical cycling, including a fully coupled carbon cycle, is presented and evaluated. The model is consistent with multiple observational data sets from the past 50 years as well as with the observed warming of global surface air and sea temperatures during the last 150 years. It is applied to a simulation of the coming two millennia following a business-as-usual scenario of anthropogenic CO2 emissions (SRES A2 until year 2100 and subsequent linear decrease to zero until year 2300, corresponding to a total release of 5100 GtC). Atmospheric CO2 increases to a peak of more than 2000 ppmv near year 2300 (that is an airborne fraction of 72% of the emissions) followed by a gradual decline to ∼1700 ppmv at year 4000 (airborne fraction of 56%). Forty-four percent of the additional atmospheric CO2 at year 4000 is due to positive carbon cycle–climate feedbacks. Global surface air warms by ∼10°C, sea ice melts back to 10% of its current area, and the circulation of the abyssal ocean collapses. Subsurface oxygen concentrations decrease, tripling the volume of suboxic water and quadrupling the global water column denitrification. We estimate 60 ppb increase in atmospheric N2O concentrations owing to doubling of its oceanic production, leading to a weak positive feedback and contributing about 0.24°C warming at year 4000. Global ocean primary production almost doubles by year 4000. Planktonic biomass increases at high latitudes and in the subtropics whereas it decreases at midlatitudes and in the tropics. In our model, which does not account for possible direct impacts of acidification on ocean biology, production of calcium carbonate in the surface ocean doubles, further increasing surface ocean and atmospheric pCO2. This represents a new positive feedback mechanism and leads to a strengthening of the positive interaction between climate change and the carbon cycle on a multicentennial to millennial timescale. Changes in ocean biology become important for the ocean carbon uptake after year 2600, and at year 4000 they account for 320 ppmv or 22% of the atmospheric CO2 increase since the preindustrial era.

Transient climate-carbon simulations of planetary geoengineering

Geoengineering (the intentional modification of Earths climate) has been proposed as a means of reducing CO2-induced climate warming while greenhouse gas emissions continue. Most proposals involve managing incoming solar radiation such that future greenhouse gas forcing is counteracted by reduced solar forcing. In this study, we assess the transient climate response to geoengineering under a business-as-usual CO2 emissions scenario by using an intermediate-complexity global climate model that includes an interactive carbon cycle. We find that the climate system responds quickly to artificially reduced insolation; hence, there may be little cost to delaying the deployment of geoengineering strategies until such a time as “dangerous” climate change is imminent. Spatial temperature patterns in the geoengineered simulation are comparable with preindustrial temperatures, although this is not true for precipitation. Carbon sinks in the model increase in response to geoengineering. Because geoengineering acts to mask climate warming, there is a direct CO2-driven increase in carbon uptake without an offsetting temperature-driven suppression of carbon sinks. However, this strengthening of carbon sinks, combined with the potential for rapid climate adjustment to changes in solar forcing, leads to serious consequences should geoengineering fail or be stopped abruptly. Such a scenario could lead to very rapid climate change, with warming rates up to 20 times greater than present-day rates. This warming rebound would be larger and more sustained should climate sensitivity prove to be higher than expected. Thus, employing geoengineering schemes with continued carbon emissions could lead to severe risks for the global climate system.

The long-term effect of increasing the albedo of urban areas

Solar reflective urban surfaces (white rooftops and light-colored pavements) can increase the albedo of an urban area by about 0.1. Increasing the albedo of urban and human settlement areas can in turn decrease atmospheric temperature and could potentially offset some of the anticipated temperature increase caused by global warming. We have simulated the long-term (decadal to centennial) effect of increasing urban surface albedos using a spatially explicit global climate model of intermediate complexity. We first carried out two sets of simulations in which we increased the albedo of all land areas between 20 and 45 latitude respectively. The results of these simulations indicate a long-term global cooling effect of 3 10 15 K for each 1 m 2 of a surface with an albedo increase of 0.01. This temperature reduction corresponds to an equivalent CO2 emission reduction of about 7 kg, based on recent estimates of the amount of global warming per unit CO2 emission. In a series of additional simulations, we increased the albedo of urban locations only, on the basis of two independent estimates of the spatial extent of urban areas. In these simulations, global cooling ranged from 0.01 to 0.07 K, which corresponds to a CO2 equivalent emission reduction of 25‐150 billion tonnes of CO2.

Terrestrial Carbon Cycle Dynamics under Recent and Future Climate Change

Abstract The behavior of the terrestrial carbon cycle under historical and future climate change is examined using the University of Victoria Earth System Climate Model, now coupled to a dynamic terrestrial vegetation and global carbon cycle model. When forced by historical emissions of CO2 from fossil fuels and land-use change, the coupled climate–carbon cycle model accurately reproduces historical atmospheric CO2 trends, as well as terrestrial and oceanic uptake for the past two decades. Under six twenty-first-century CO2 emissions scenarios, both terrestrial and oceanic carbon sinks continue to increase, though terrestrial uptake slows in the latter half of the century. Climate–carbon cycle feedbacks are isolated by comparing a coupled model run with a run where climate and the carbon cycle are uncoupled. The modeled positive feedback between the carbon cycle and climate is found to be relatively small, resulting in an increase in simulated CO2 of 60 ppmv at the year 2100. Including non-CO2 greenhouse ...

Sensitivity of ocean acidification to geoengineered climate stabilization

] Climate engineering has been proposed as a possibleresponse to anthropogenic climate change. While climateengineering may be able to stabilize temperatures, it isgenerally assumed that this will not prevent continued oceanacidification. However, due to the strong coupling betweenclimate and the carbon cycle, climate engineering couldindirectly affect ocean chemistry. We used a global Earth-system model to investigate how climate engineering mayaffectsurfaceoceanpHandthedegreeofaragonitesaturation.Climate engineering could significantly re-distribute carbonemissions among atmosphere, land and ocean reservoirs.This could slow pH decreases somewhat relative to the non-engineered case, but would not affect the level of aragonitesaturation due to opposing responses of pH and aragonitesaturation to temperature change. However, these effects aredependent on enhanced carbon accumulation in the landbiosphere; without this, climate engineering has little effecton pH, and leads to accelerated declines in aragonitesaturation.

Accounting for carbon cycle feedbacks in a comparison of the global warming effects of greenhouse gases

Greenhouse gases other than CO2 make a significant contribution to human-induced climate change, and multi-gas mitigation strategies are cheaper to implement than those which limit CO2 emissions alone. Most practical multi-gas mitigation strategies require metrics to relate the climate warming effects of CO2 and other greenhouse gases. Global warming potential (GWP), defined as the ratio of time-integrated radiative forcing of a particular gas to that of CO2 following a unit mass emission, is the metric used in the Kyoto Protocol, and we define mean global temperature change potential (MGTP) as an equivalent metric of the temperature response. Here we show that carbon–climate feedbacks inflate the GWPs and MGTPs of methane and nitrous oxide by ~ 20% in coupled carbon–climate model simulations of the response to a pulse of 50 × 1990 emissions, due to a warming-induced release of CO2 from the land biosphere and ocean. The magnitude of this effect is expected to be dependent on the model, but it is not captured at all by the analytical models usually used to calculate metrics such as GWP. We argue that the omission of carbon cycle dynamics has led to a low bias of uncertain but potentially substantial magnitude in metrics of the global warming effect of other greenhouse gases, and we suggest that the carbon–climate feedback should be considered when greenhouse gas metrics are calculated and applied.

Cumulative carbon as a policy framework for achieving climate stabilization.

The primary objective of the United Nations Framework Convention on Climate Change is to stabilize greenhouse gas concentrations at a level that will avoid dangerous climate impacts. However, greenhouse gas concentration stabilization is an awkward framework within which to assess dangerous climate change on account of the significant lag between a given concentration level and the eventual equilibrium temperature change. By contrast, recent research has shown that global temperature change can be well described by a given cumulative carbon emissions budget. Here, we propose that cumulative carbon emissions represent an alternative framework that is applicable both as a tool for climate mitigation as well as for the assessment of potential climate impacts. We show first that both atmospheric CO2 concentration at a given year and the associated temperature change are generally associated with a unique cumulative carbon emissions budget that is largely independent of the emissions scenario. The rate of global temperature change can therefore be related to first order to the rate of increase of cumulative carbon emissions. However, transient warming over the next century will also be strongly affected by emissions of shorter lived forcing agents such as aerosols and methane. Non-CO2 emissions therefore contribute to uncertainty in the cumulative carbon budget associated with near-term temperature targets, and may suggest the need for a mitigation approach that considers separately short- and long-lived gas emissions. By contrast, long-term temperature change remains primarily associated with total cumulative carbon emissions owing to the much longer atmospheric residence time of CO2 relative to other major climate forcing agents.

A climate sensitivity estimate using Bayesian fusion of instrumental observations and an Earth System model

[1] Current climate model projections are uncertain. This uncertainty is partly driven by the uncertainty in key model parameters such as climate sensitivity (CS), vertical ocean diffusivity (Kv), and strength of anthropogenic sulfate aerosol forcing. These parameters are commonly estimated using ensembles of model runs constrained by observations. Here we obtain a probability density function (pdf) of these parameters using the University of Victoria Earth System Climate Model (UVic ESCM) - an intermediate complexity model with a dynamic three-dimensional ocean. Specifically, we run an ensemble of UVic ESCM runs varying parameters that affect CS, ocean vertical diffusion, and the effects of anthropogenic sulfate aerosols. We use a statistical emulator that interpolates the UVic ESCM output to parameter settings where the model was not evaluated. We adopt a Bayesian approach to constrain the model output with instrumental surface temperature and ocean heat observations. Our approach accounts for the uncertainties in the properties of model-data residuals. We use a Markov chain Monte Carlo method to obtain a posterior pdf of these parameters. The mode of the climate sensitivity estimate is 2.8°C, with the corresponding 95% credible interval ranging from 1.8 to 4.9°C. These results are generally consistent with previous studies. The CS pdf is sensitive to the assumptions about the priors, to the effects of anthropogenic sulfate aerosols, and to the background vertical ocean diffusivity. Our method can be used with more complex climate models.