Andreas Schmittner

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.

Global warming preceded by increasing carbon dioxide concentrations during the last deglaciation

The covariation of carbon dioxide (CO2) concentration and temperature in Antarctic ice-core records suggests a close link between CO2 and climate during the Pleistocene ice ages. The role and relative importance of CO2 in producing these climate changes remains unclear, however, in part because the ice-core deuterium record reflects local rather than global temperature. Here we construct a record of global surface temperature from 80 proxy records and show that temperature is correlated with and generally lags CO2 during the last (that is, the most recent) deglaciation. Differences between the respective temperature changes of the Northern Hemisphere and Southern Hemisphere parallel variations in the strength of the Atlantic meridional overturning circulation recorded in marine sediments. These observations, together with transient global climate model simulations, support the conclusion that an antiphased hemispheric temperature response to ocean circulation changes superimposed on globally in-phase warming driven by increasing CO2 concentrations is an explanation for much of the temperature change at the end of the most recent ice age.

Influence of CO2 emission rates on the stability of the thermohaline circulation

Present estimates of the future oceanic uptake of anthropogenic CO2 and calculations of CO2-emission scenarios are based on the assumption that the natural carbon cycle is in steady state. But it iswell known from palaeoclimate records and modelling studies that the climate system has more than one equilibrium state, and that perturbations can trigger transitions between them. Anticipated future changes in todays climate system due to human activities have the potential to weaken the thermohaline circulation of the North Atlantic Ocean, which would greatly modify estimates of future oceanic CO2 uptake. Here we use a simple coupled atmosphere–ocean climate model to show that the Atlantic thermohaline circulation is not only sensitive to the final atmospheric CO2 concentration attained, but also depends on the rate of change of the CO2 concentration in the atmosphere. A modelled increase to 750 parts per million by volume (p.p.m.v.) CO2 within 100 years (corresponding approximately to a continuation of todays growth rate) leads to a permanent shut-down of the thermohaline circulation. If the final atmospheric concentration of 750 p.p.m.v. CO2 is attained more slowly, the thermohaline circulation simply slows down. The reason for this rate-sensitive response of the climate system lies with the transfer of buoyancy in the form of heat and fresh water from the uppermost layers of the ocean into the deep waters below. This sensitivity of the simulated thermohaline circulation to the rate ofchange of atmospheric CO2 concentration has potentially important implications for the choice of future CO2-emission scenarios.

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.

Decline of the marine ecosystem caused by a reduction in the Atlantic overturning circulation

Reorganizations of the Atlantic meridional overturning circulation were associated with large and abrupt climatic changes in the North Atlantic region during the last glacial period. Projections with climate models suggest that similar reorganizations may also occur in response to anthropogenic global warming. Here I use ensemble simulations with a coupled climate–ecosytem model of intermediate complexity to investigate the possible consequences of such disturbances to the marine ecosystem. In the simulations, a disruption of the Atlantic meridional overturning circulation leads to a collapse of the North Atlantic plankton stocks to less than half of their initial biomass, owing to rapid shoaling of winter mixed layers and their associated separation from the deep ocean nutrient reservoir. Globally integrated export production declines by more than 20 per cent owing to reduced upwelling of nutrient-rich deep water and gradual depletion of upper ocean nutrient concentrations. These model results are consistent with the available high-resolution palaeorecord, and suggest that global ocean productivity is sensitive to changes in the Atlantic meridional overturning circulation.

Model projections of the North Atlantic thermohaline circulation for the 21st century assessed by observations

Most climate models predict a weakening of the North Atlantic thermohaline circulation for the 21st century when forced by increasing levels of greenhouse gas concentrations. The model spread, however, is rather large, even when the forcing scenario is identical, indicating a large uncertainty in the response to forcing. In order to reduce the model uncertainties a weighting procedure is applied considering the skill of each model in simulating hydrographic properties and observation-based circulation estimates. This procedure yields a “best estimate” for the evolution of the North Atlantic THC during the 21st century by taking into account a measure of model quality. Using 28 projections from 9 different coupled global climate models of a scenario of future CO2 increase (SRESA1B) performed for the upcoming fourth assessment report of the Intergovernmental Panel on Climate Change, the analysis predicts a gradual weakening of the North Atlantic THC by 25(±25)% until 2100.

Simulated 21st century's increase in oceanic suboxia by CO2-enhanced biotic carbon export

The primary impacts of anthropogenic CO2 emissions on marine biogeochemical cycles predicted so far include ocean acidification, global warming induced shifts in biogeographical provinces, and a possible negative feedback on atmospheric CO2 levels by CO2‐fertilized biological production. Here we report a new potentially significant impact on the oxygen‐minimum zones of the tropical oceans. Using a model of global climate, ocean circulation, and biogeochemical cycling, we extrapolate mesocosm‐derived experimental findings of a pCO2‐sensitive increase in biotic carbon‐to‐nitrogen drawdown to the global ocean. For a simulation run from the onset of the industrial revolution until A.D. 2100 under a “business‐as‐usual” scenario for anthropogenic CO2 emissions, our model predicts a negative feedback on atmospheric CO2 levels, which amounts to 34 Gt C by the end of this century. While this represents a small alteration of the anthropogenic perturbation of the carbon cycle, the model results reveal a dramatic 50% increase in the suboxic water volume by the end of this century in response to the respiration of excess organic carbon formed at higher CO2 levels. This is a significant expansion of the marine “dead zones” with severe implications not only for all higher life forms but also for oxygen‐sensitive nutrient recycling and, hence, for oceanic nutrient inventories.

Ice-shelf collapse from subsurface warming as a trigger for Heinrich events

Episodic iceberg-discharge events from the Hudson Strait Ice Stream (HSIS) of the Laurentide Ice Sheet, referred to as Heinrich events, are commonly attributed to internal ice-sheet instabilities, but their systematic occurrence at the culmination of a large reduction in the Atlantic meridional overturning circulation (AMOC) indicates a climate control. We report Mg/Ca data on benthic foraminifera from an intermediate-depth site in the northwest Atlantic and results from a climate-model simulation that reveal basin-wide subsurface warming at the same time as large reductions in the AMOC, with temperature increasing by approximately 2 °C over a 1–2 kyr interval prior to a Heinrich event. In simulations with an ocean model coupled to a thermodynamically active ice shelf, the increase in subsurface temperature increases basal melt rate under an ice shelf fronting the HSIS by a factor of approximately 6. By analogy with recent observations in Antarctica, the resulting ice-shelf loss and attendant HSIS acceleration would produce a Heinrich event.

Climate Sensitivity Estimated from Temperature Reconstructions of the Last Glacial Maximum

Last Glacial Maximum temperature reconstructions and model simulations can constrain the equilibrium climate sensitivity. Assessing the impact of future anthropogenic carbon emissions is currently impeded by uncertainties in our knowledge of equilibrium climate sensitivity to atmospheric carbon dioxide doubling. Previous studies suggest 3 kelvin (K) as the best estimate, 2 to 4.5 K as the 66% probability range, and nonzero probabilities for much higher values, the latter implying a small chance of high-impact climate changes that would be difficult to avoid. Here, combining extensive sea and land surface temperature reconstructions from the Last Glacial Maximum with climate model simulations, we estimate a lower median (2.3 K) and reduced uncertainty (1.7 to 2.6 K as the 66% probability range, which can be widened using alternate assumptions or data subsets). Assuming that paleoclimatic constraints apply to the future, as predicted by our model, these results imply a lower probability of imminent extreme climatic change than previously thought.

Glacial greenhouse-gas fluctuations controlled by ocean circulation changes

Earth’s climate and the concentrations of the atmospheric greenhouse gases carbon dioxide (CO2) and nitrous oxide (N2O) varied strongly on millennial timescales during past glacial periods. Large and rapid warming events in Greenland and the North Atlantic were followed by more gradual cooling, and are highly correlated with fluctuations of N2O as recorded in ice cores. Antarctic temperature variations, on the other hand, were smaller and more gradual, showed warming during the Greenland cold phase and cooling while the North Atlantic was warm, and were highly correlated with fluctuations in CO2. Abrupt changes in the Atlantic meridional overturning circulation (AMOC) have often been invoked to explain the physical characteristics of these Dansgaard–Oeschger climate oscillations, but the mechanisms for the greenhouse-gas variations and their linkage to the AMOC have remained unclear. Here we present simulations with a coupled model of glacial climate and biogeochemical cycles, forced only with changes in the AMOC. The model simultaneously reproduces characteristic features of the Dansgaard–Oeschger temperature, as well as CO2 and N2O fluctuations. Despite significant changes in the land carbon inventory, CO2 variations on millennial timescales are dominated by slow changes in the deep ocean inventory of biologically sequestered carbon and are correlated with Antarctic temperature and Southern Ocean stratification. In contrast, N2O co-varies more rapidly with Greenland temperatures owing to fast adjustments of the thermocline oxygen budget. These results suggest that ocean circulation changes were the primary mechanism that drove glacial CO2 and N2O fluctuations on millennial timescales.