Modeling the eco-hydrologic response of a Mediterranean type ecosystem to the combined impacts of projected climate change and altered fire frequencies

Global Climate Models (GCMs) project moderate warming along with increases in atmospheric CO₂ for California Mediterranean type ecosystems (MTEs). In water-limited ecosystems, vegetation acts as an important control on streamflow and responds to soil moisture availability. Fires are also key disturbances in semi-arid environments, and few studies have explored the potential interactions among changes in climate, vegetation dynamics, hydrology, elevated atmospheric CO2 concentrations and fire. We model ecosystem productivity, evapotranspiration, and summer streamflow under a range of temperature and precipitation scenarios using RHESSys, a spatially distributed model of carbon–water interactions. We examine the direct impacts of temperature and precipitation on vegetation productivity and impacts associated with higher water-use efficiency under elevated atmospheric CO2. Results suggest that for most climate scenarios, biomass in chaparral-dominated systems is likely to increase, leading to reductions in summer streamflow. However, within the range of GCM predictions, there are some scenarios in which vegetation may decrease, leading to higher summer streamflows. Changes due to increases in fire frequency will also impact summer streamflow but these will be small relative to changes due to vegetation productivity. Results suggest that monitoring vegetation responses to a changing climate should be a focus of climate change assessment for California MTEs.     

The potential of vegetation feedback to alleviate climate aridity over the United States associated with a 2×CO2 climate condition

This study examines the potential impact of vegetation feedback on changes in summer climate aridity over the contiguous United States (US) due to the doubling of atmospheric CO₂ concentration using a set of 100-year-long climate simulations made by a global climate model interactively coupled with a dynamic vegetation model. The Thornthwaite moisture index (I _m ), which quantifies climate aridity on the basis of atmospheric water supply (i.e., precipitation) and atmospheric water demand (i.e., potential evapotranspiration, PET), is used to measure climate aridity. Warmer atmosphere and drier surface resulting from increased CO₂ concentration increase climate aridity over most of the contiguous US. This phenomenon is due to larger increments in PET than in precipitation, regardless of the presence or absence of vegetation feedback. Compared to simulations without active dynamic vegetation feedback, the presence of vegetation feedback significantly alleviates the increase in aridity. This vegetation-feedback effect is most noticeable in the subhumid regions such as southern, midwestern and northwestern US, primarily by the increasing vegetation greenness. In these regions, the greening in response to warmer temperatures enhances moisture transfer from soil to atmosphere by evapotranspiration (ET). The increased ET and subsequent moistening over land areas result in weaker surface warming (1–2?K) and PET (3–10?mm?month^?1), and greater precipitation (4–10?mm?month^?1). Collectively, they result in moderate increases in I _m . Our results suggest that moistening by enhanced vegetation feedback may prevent aridification under climatic warming especially in areas vulnerable to climate change, with consequent implications for mitigation strategies.     

SEASONAL CHANGES OF SOIL MOISTURE INDUCED BY INCREASED CO2 AS SIMULATED BY THE ATMOSPHERIC GENERAL CIRCULATION AND MIXED LAYER OCEAN MODEL

The atmospheric general circulation model coupled to the mixed layer ocean model has been used to simulate thechanges of the global soil moisture.Comparing the simulated results with observations,it is shown that the model is ca-pable of doing sensitive experiments about the carbon dioxide change.The 2×CO_2/1×CO_2 comparison shows that there are the obvious changes of the soil moisture in the global forfour seasons.There are the wet soil moisture in the lower latitudes of both hemispheres and dry soil moisture in the mid-dle latitudes of both hemispheres for four seasons.The dry soil moisture in summer and wet in other seasons are foundin the northern higher latitudes.The analyses of the physical feedbacks responsible for the CO_2-induced changes of soil moisture show that the bud-gets of the surface water and heat are the important factors.     

Influence of dynamic vegetation on climate change arising from increasing CO<Subscript>2</Subscript>

A dynamic global vegetation model (DGVM) is coupled to an atmospheric general circulation model (AGCM) to investigate the influence of vegetation dynamics on climate change under conditions of global warming. The model results are largely in agreement with observations and the results of previous studies in terms of the present climate, present potential vegetation, present net primary productivity (NPP), and pre-industrial carbon budgets. The equilibrium state of climate properties are compared among pre-industrial, doubled, and quadrupled atmospheric CO₂ values using DGVM–AGCM and current AGCM with fixed vegetation to evaluate the influence of dynamic vegetation change. We also separated the contributions of temperature, precipitation and CO₂ fertilization on vegetation change. The results reveal an amplification of global warming climate sensitivity by 10% due to the inclusion of dynamic vegetation. The total effects of elevated CO₂ and climate change also lead to an increase in NPP and vegetation coverage globally. The reduction of albedo associated with this greening results in enhanced global warming. Our separation analysis indicates that temperature alters vegetation at high latitudes such as Siberia or Alaska, where there is a switch from tundra to forest. On the other hand, CO₂ fertilization provides the largest contribution to greening in arid/semi-arid region. Precipitation change did not cause any drastic vegetation shift.     

Sensitivity of FORGRO to climatic change scenarios: A case study on Betula pubescens, Fagus sylvatica and Quercus robur in The Netherlands

The impacts of the climate change predictions of four general circulation models (GFDL, GISS, OSU and UKMO) on net primary production (NPP) ofBetula pubescens, Fagus sylvatica and Quercus robur in The Netherlands were analysed using the process-based model FORGRO. FORGRO is a model suitable to simulate growth of managed mono-species stands. For the GCMs mentioned, both transient and equilibrium 2 × CO₂ scenarios of temperature and precipitation change were evaluated and compared with responses under current climate. It was found that the NPP increases in the transient scenarios, but remains the same or declines in the 2 × CO₂ scenarios. This is because respiration increases more with rising temperature than photosynthesis. During the transient scenarios this effect gradually increases, while in the 2 × CO₂ scenario this effect is operating over the entire simulation period.If water limitation is taken into account, then the NPP of the reference scenario is reduced. In both the transient and 2 × CO₂ scenarios mis water limitation is annulated, resulting in a stronger response of NPP compared to the situation without water limitation. This enhancement of the response is most pronounced in the transient scenario due to the gradual effect of temperature on respiration.Similar results were obtained with a version of FORGRO in which the photosynthesis module of HYBRID (PGEN) is incorporated, although the response in FORGRO-PGEN is usually higher than that of FORGRO. This is because the response of photosynthesis to CO₂ rises with increasing temperature as defined in the PGEN-model, but not according to FORGRO.     

Modeling responses of the meadow steppe dominated by Leymus chinensis to climate change

Grassland is one of the most widespread vegetation types worldwide and plays a significant role in regional climate and global carbon cycling. Understanding the sensitivity of Chinese grassland ecosystems to climate change and elevated atmospheric CO₂ and the effect of these changes on the grassland ecosystems is a key issue in global carbon cycling. China encompasses vast grassland areas of 354 million ha of 17 major grassland types, according to a national grassland survey. In this study, a process-based terrestrial model the CENTURY model was used to simulate potential changes in net primary productivity (NPP) and soil organic carbon (SOC) of the Leymus chinensis meadow steppe (LCMS) under different scenarios of climatic change and elevated atmospheric CO₂. The LCMS sensitivities, its potential responses to climate change, and the change in capacity of carbon stock and sequestration in the future are evaluated. The results showed that the LCMS NPP and SOC are sensitive to climatic change and elevated CO₂. In the next 100 years, with doubled CO₂ concentration, if temperature increases from 2.7-3.9˚C and precipitation increases by 10% NPP and SOC will increase by 7-21% and 5-6% respectively. However, if temperature increases by 7.5-7.8˚C and precipitation increases by only 10% NPP and SOC would decrease by 24% and 8% respectively. Therefore, changes in the NPP and SOC of the meadow steppe are attributed mainly to the amount of temperature and precipitation change and the atmospheric CO₂ concentration in the future.     

Long-term effects of anthropogenic CO<Subscript>2</Subscript> emissions simulated with a complex earth system model

A new complex earth system model consisting of an atmospheric general circulation model, an ocean general circulation model, a three-dimensional ice sheet model, a marine biogeochemistry model, and a dynamic vegetation model was used to study the long-term response to anthropogenic carbon emissions. The prescribed emissions follow estimates of past emissions for the period 1751–2000 and standard IPCC emission scenarios up to the year 2100. After 2100, an exponential decrease of the emissions was assumed. For each of the scenarios, a small ensemble of simulations was carried out. The North Atlantic overturning collapsed in the high emission scenario (A2) simulations. In the low emission scenario (B1), only a temporary weakening of the deep water formation in the North Atlantic is predicted. The moderate emission scenario (A1B) brings the system close to its bifurcation point, with three out of five runs leading to a collapsed North Atlantic overturning circulation. The atmospheric moisture transport predominantly contributes to the collapse of the deep water formation. In the simulations with collapsed deep water formation in the North Atlantic a substantial cooling over parts of the North Atlantic is simulated. Anthropogenic climate change substantially reduces the ability of land and ocean to sequester anthropogenic carbon. The simulated effect of a collapse of the deep water formation in the North Atlantic on the atmospheric CO₂ concentration turned out to be relatively small. The volume of the Greenland ice sheet is reduced, but its contribution to global mean sea level is almost counterbalanced by the growth of the Antarctic ice sheet due to enhanced snowfall. The modifications of the high latitude freshwater input due to the simulated changes in mass balance of the ice sheet are one order of magnitude smaller than the changes due to atmospheric moisture transport. After the year 3000, the global mean surface temperature is predicted to be almost constant due to the compensating effects of decreasing atmospheric CO₂ concentrations due to oceanic uptake and delayed response to increasing atmospheric CO₂ concentrations before.     

Effects of climate change and elevated atmospheric CO2 on soil organic carbon: a response equation

Climate change, such as warming and precipitation change, as well as elevated CO₂ can affect soil organic carbon (SOC) dynamics and cause changes in soil carbon sequestration. In this study, we introduced a response equation, relating the relative change of SOC to the relative changes of annual average temperature, annual precipitation, and atmospheric CO₂ concentration, as well as their inter-products. Using Nelson Farm as a case study, based on simulations of CENTURY model and multiple regressions, we examined the response equation for three vegetation covers (i.e., soybean, corn, and grass) and scenarios with different soil erosion rates and initial SOC contents. The response equation fit the simulation results very well with high adjusted coefficients of determination (R ²) (0.982 to 0.990). The results showed that the SOC was negatively related to the annual average temperature, positively related to the annual precipitation, and positively related to the elevated CO₂ for all the vegetation covers (p?<?0.001). The SOC was also significantly impacted by the interaction effects between elevated CO₂ and warming or precipitation change (p?<?0.001). The general form of the response equations for the different vegetation covers, soil erosion rates, and initial SOC contents was the same although the parameters varied with the different conditions. Based on the response equation, ??cutoff surfaces?? were defined to clearly quantify the synthesis effects of any possible combination of climate change and elevated CO₂ on the SOC, and the SOC sequestration potential was assessed under climate change and elevated CO₂ for different vegetations. Compared with the empirical models in the literature, this response equation provides a simple yet but robust method to represent the relationship between the SOC relative change vs. the relative changes of atmospheric temperature, precipitation, and atmospheric CO₂ concentration.     

Hydrological changes in the climate system from leaf responses to increasing CO2

Vegetation is a major component of the climate system because of its controls on the energy and water balance over land. This functioning changes because of the physiological response of leaves to increased CO₂. A climate model is used to compare these changes with the climate changes from radiative forcing by greenhouse gases. For this purpose, we use the Community Earth System Model coupled to a slab ocean. Ensemble integrations are done for current and doubled CO₂. The consequent reduction of transpiration and net increase of surface radiative heating from reduction in cloudiness increases the temperature over land by a significant fraction of that directly from the radiative warming by CO₂. Large-scale atmospheric circulation adjustments result. In particular, over the tropics, a low-level westerly wind anomaly develops associated with reduced geopotential height over land, enhancing moisture transport and convergence, and precipitation increases over the western Amazon, the Congo basin, South Africa, and Indonesia, while over mid-latitudes, land precipitation decreases from reduced evapotranspiration. On average, land precipitation is enhanced by 0.03 mm day^?1 (about 19 % of the CO₂ radiative forcing induced increase). This increase of land precipitation with decreased ET is an apparent negative feedback, i.e., less ET makes more precipitation. Global precipitation is slightly reduced. Runoff increases associated with both the increased land precipitation and reduced evapotranspiration. Examining the consistency of the variations among ensemble members shows that vegetation feedbacks on precipitation are more robust over the tropics and in mid to high latitudes than over the subtropics where vegetation is sparse and the internal climate variability has a larger influence.     

用全球大气混合层海洋环流模式模拟二氧化碳增加对土壤湿度季节变化的影响

本文利用全球三维大气耦合混合层海洋环流模式模拟大气中二氧化碳浓度增加对土壤湿度的影响。敏感试验(2×CO_2)与控制试验(1×CO_2)对照表明,当大气中二氧化碳浓度增加时,全球土壤湿度在各季发生明显变化。其中两半球低纬度地区在冬季土壤温度变温,两半球中纬度地区则在各季土壤湿度变干,北半球高纬度地区土壤湿度在夏季变干,其余各季变温。分析大气中二氧化碳浓度增加造成土壤温度全球变化的可能物理机制表明,地面水循环和热量循环是重要的因素。     

Hydrological effects of the increased CO2 and climate change in the Upper Mississippi River Basin using a modified SWAT

Increased atmospheric CO₂ concentration and climate change may significantly impact the hydrological and meteorological processes of a watershed system. Quantifying and understanding hydrological responses to elevated ambient CO₂ and climate change is, therefore, critical for formulating adaptive strategies for an appropriate management of water resources. In this study, the Soil and Water Assessment Tool (SWAT) model was applied to assess the effects of increased CO₂ concentration and climate change in the Upper Mississippi River Basin (UMRB). The standard SWAT model was modified to represent more mechanistic vegetation type specific responses of stomatal conductance reduction and leaf area increase to elevated CO₂ based on physiological studies. For estimating the historical impacts of increased CO₂ in the recent past decades, the incremental (i.e., dynamic) rises of CO₂ concentration at a monthly time-scale were also introduced into the model. Our study results indicated that about 1–4% of the streamflow in the UMRB during 1986 through 2008 could be attributed to the elevated CO₂ concentration. In addition to evaluating a range of future climate sensitivity scenarios, the climate projections by four General Circulation Models (GCMs) under different greenhouse gas emission scenarios were used to predict the hydrological effects in the late twenty-first century (2071–2100). Our simulations demonstrated that the water yield would increase in spring and substantially decrease in summer, while soil moisture would rise in spring and decline in summer. Such an uneven distribution of water with higher variability compared to the baseline level (1961–1990) may cause an increased risk of both flooding and drought events in the basin.     

The Effect of Spatial Scale of Climatic Change Scenarios on Simulated Maize, Winter Wheat, and Rice Production in the Southeastern United States

We use the CERES family of crop models to assess the effect of different spatial scales of climate change scenarios on the simulated yield changes of maize (Zea mays L.), winter wheat (Triticum aestivum L.),and rice (Oryza sativa L.) in the Southeastern United States. The climate change scenarios were produced with the control and doubled CO₂ runs of a high resolution regional climate model anda coarse resolution general circulation model, which provided the initial and lateral boundary conditions for the regional model. Three different cases were considered for each scenario: climate change alone, climate change plus elevated CO₂, and the latter with adaptations. On the state level,for most cases, significant differences in the climate changed yields for corn were found, the coarse scale scenario usually producing larger modeled yield decreases or smaller increases. For wheat, however, which suffered large decreases in yields for all cases, very little contrast in yield based on scale of scenario was found. Scenario scale resulted in significantly different rice yields, but mainly because of low variability in yields. For maize the primary climate variable that explained the contrast in the yields calculated from the two scenarios is the precipitation during grain fill leading to different water stress levels. Temperature during vernalization explains some contrasts in winter wheat yields. With adaptation, the contrasts in the yields of all crops produced by the scenarios were reduced but not entirely removed. Our results indicate that spatial resolution of climate change scenarios can be an important uncertainty in climate change impact assessments, depending on the crop and management conditions.     

Climate Change in Northern Africa: The Past is Not the Future

By using a climate system model of intermediate complexity, we have simulated long-term natural climate changes occurring over the last 9000 years. The paleo-simulations in which the model is driven by orbital forcing only, i.e., by changes in insolation caused by changes in the Earth's orbit, are compared with sensitivity simulations in which various scenarios of increasing atmospheric CO₂ concentration are prescribed. Focussing on climate and vegetation change in northern Africa, we recapture the strong greening of the Sahara in the early and mid-Holocene (some 9000–6000 years ago), and we show that some expansion of grasslandinto the Sahara is theoretically possible, if the atmospheric CO₂ concentration increases well above pre-industrial values and if vegetation growth is not disturbed. Depending on the rate of CO₂ increase, vegetation migration into the Sahara can be rapid, up to 1/10th of the Saharan area per decade, but could not exceed a coverage of 45%. In ourmodel, vegetation expansion into today's Sahara is triggered by an increase in summer precipitation which is amplified by a positive feedback between vegetation and precipitation. This is valid for simulations with orbital forcing and greenhouse-gas forcing. However, we argue that the mid-Holocene climate optimum some 9000 to 6000 years ago with its marked reduction of deserts in northern Africa is not a direct analogue for future greenhouse-gas induced climate change, as previously hypothesized. Not only does the global pattern of climate change differ between the mid-Holocene model experiments and the greenhouse-gas sensitivity experiments, but the relative role of mechanisms which lead to a reduction of the Sahara also changes. Moreover, the amplitude of simulated vegetation cover changes in northern Africa is less than is estimated for mid-Holocene climate.     

Climatic change and the broad-scale distribution of terrestrial ecosystem complexes

The broad-scale distribution of terrestrial ecosystem complexes is determined in large part by climate and can be altered by climatic change due to natural causes or due to human activities such as those leading to increasing atmospheric CO₂ concentration. Classifications that recognize the dependence of natural vegetation on climate provide one means of constructing maps to display the impact of climatic change on the geography of major vegetation zones. A world map of the Holdridge Life-Zone Classification, developed from approximately 8,000 meteorological records, is compared with a Holdridge Map with average temperature increments simulated by a. model of climate under elevated atmospheric CO₂ concentration. The largest changes are indicated at high latitudes, where the simulated temperature increase is largest and the temperature intervals defining life zones are smallest. Boreal Forest Zones are replaced by either Cool Temperate Forest or Cool Temperate Steppe, depending on average precipitation. Changes in the tropics are smaller; however, in some regions, Subtropical Moist Forest is replaced by Tropical Dry Forest.Research supported by the National Science Foundation's Ecosystem Studies Program under Interagency Agreement Nos. DEB81-15316 and DEB83-15185.     

Climate change impacts on Laurentian Great Lakes levels

Scenarios of water supplies reflecting CO₂-induced climatic change are used to determine potential impacts on levels of the Laurentian Great Lakes and likely water management policy implications. The water supplies are based on conceptual models that link climate change scenarios from general circulation models to estimates of basin runoff, overlake precipitation, and lake evaporation. The water supply components are used in conjunction with operational regulation plans and hydraulic routing models of outlet and connecting channel flows to estimate water levels on Lakes Superior, Michigan, Huron, St. Clair, Erie, and Ontario. Three steady-state climate change scenarios, corresponding to modeling a doubling of atmospheric CO₂, are compared to a steady-state simulation obtained with historical data representing an unchanged atmosphere. One transient climate change scenario, representing a modeled transition from present conditions to doubled CO₂ concentrations, is compared to a transient simulation with historical data. The environmental, socioeconomic, and policy implications of the climate change effects modeled herein suggest that new paradigms in water management will be required to address the prospective increased allocation conflicts between users of the Great Lakes.GLERL Contribution No. 645.     

Effects of increased CO2 levels on monsoons

Increased atmospheric carbon dioxide concentration provided warmer atmospheric temperature and higher atmospheric water vapor content, but not necessarily more precipitation. A set of experiments performed with a state-of-the-art coupled general circulation model forced with increased atmospheric CO₂ concentration (2, 4 and 16 times the present-day mean value) were analyzed and compared with a control experiment to evaluate the effect of increased CO₂ levels on monsoons. Generally, the monsoon precipitation responses to CO₂ forcing are largest if extreme concentrations of carbon dioxide are used, but they are not necessarly proportional to the forcing applied. In fact, despite a common response in terms of an atmospheric water vapor increase to the atmospheric warming, two out of the six monsoons studied simulate less or equal summer mean precipitation in the 16×CO₂ experiment compared to the intermediate sensitivity experiments. The precipitation differences between CO₂ sensitivity experiments and CTRL have been investigated specifying the contribution of thermodynamic and purely dynamic processes. As a general rule, the differences depending on the atmospheric moisture content changes (thermodynamic component) are large and positive, and they tend to be damped by the dynamic component associated with the changes in the vertical velocity. However, differences are observed among monsoons in terms of the role played by other terms (like moisture advection and evaporation) in shaping the precipitation changes in warmer climates. The precipitation increase, even if weak, occurs despite a weakening of the mean circulation in the monsoon regions (??precipitation-wind paradox??). In particular, the tropical east-west Walker circulation is reduced, as found from velocity potential analysis. The meridional component of the monsoon circulation is changed as well, with larger (smaller) meridional (vertical) scales.     

Climate Change Impacts for the Conterminous USA: An Integrated Assessment

As carbon dioxide and other greenhouse gases accumulate in the atmosphere and contribute to rising global temperatures, it is important to examine how derivative changes in climate may affect natural and managed ecosystems. In this series of papers, we study the impacts of climate change on agriculture, water resources and natural ecosystems in the conterminous United States using twelve scenarios derived from General Circulation Model (GCM) projections to drive biophysical impact models. These scenarios are described in this paper. The scenarios are first put into the context of recent work on climate-change by the IPCC for the 21st century and span two levels of global-mean temperature change and three sets of spatial patterns of change derived from GCM results. In addition, the effect of either the presence or absence of a CO₂ fertilization effect on vegetation is examined by using two levels of atmospheric CO₂ concentration as a proxy variable. Results from three GCM experiments were used to produce different regional patterns of climate change. The three regional patterns for the conterminous United States range from: an increase in temperature above the global-mean level along with a significant decline in precipitation; temperature increases in line with the global-mean with an average increase in precipitation; and, with a sulfate aerosol effect added to in the same model, temperature increases that are lower than the global-mean. The resulting set of scenarios span a wide range of potential climate changes and allows examination of the relative importance of global-mean temperature change, regional climate patterns, aerosol cooling, and CO₂ fertilization effects.     

1981-2008年中国陆地植被NPP对气候变化响应的敏感性试验

在验证CENTURY模型对中国陆地植被净初级生产力(Net Primary Productivity,NPP)模拟能力的基础上,利用该模型探讨了1981-2008年中国陆地植被NPP的年际变异和变化趋势对CO₂浓度、温度和降水变化的响应。结果表明,中国陆地植被NPP对不同气候因子的响应程度存在明显不同。其中,CO₂浓度变化对植被NPP年际变异的影响不显著,但能够引起中国大部分地区植被NPP趋势系数增大;温度对中国中高纬度地区植被NPP的年际变化影响显著,但就全国范围而言,植被NPP年际变异对温度变化的响应程度总体低于对降水变化的响应程度;降水变化是对中国植被NPP变化趋势起主导作用的气候因子。此外,综合考虑温度和降水变化的影响发现,植被NPP变化趋势的响应特征类似于降水单独变化时植被NPP变化趋势的响应特征。     

Recent climatic change in Australia: Implications for a CO2-warmed earth

A significant change in mean precipitation occurred over much of Australia between 1913–45 and 1946–78. This is described on a seasonal basis and related to possible changes in the atmospheric circulation. It now appears that during this time mean surface temperatures in the mid southern latitude zone increased by up to 1 °C. This temperature change could be at least partly due to an increase in atmospheric CO₂ concentrations from about 260 ppmv in the early nineteenth century. In any case the observed temperature increase is similar to the predicted future effects of a 50% increase in atmospheric CO₂ concentrations. Thus the climatic change which occurred earlier this century is at least a good analogy for the effects of a CO₂-induced global warming which is expected to occur over a similar time interval in the future. This allows the construction of more detailed and quantitative climate scenarios. The most noteworthy conclusion is that marked changes in the seasonally of precipitation should be anticipated, with seasonal changes in some areas being of the order of 50% or more for a doubling of CO₂ content. The results are in general consistent with earlier more qualitative scenarios for Australia.     

CO2 threshold for millennial-scale oscillations in the climate system: implications for global warming scenarios

We present several equilibrium runs under varying atmospheric CO₂ concentrations using the University of Victoria Earth System Climate Model (UVic ESCM). The model shows two very different responses: for CO₂ concentrations of 400 ppm or lower, the system evolves into an equilibrium state. For CO₂ concentrations of 440 ppm or higher, the system starts oscillating between a state with vigorous deep water formation in the Southern Ocean and a state with no deep water formation in the Southern Ocean. The flushing events result in a rapid increase in atmospheric temperatures, degassing of CO₂ and therefore an increase in atmospheric CO₂ concentrations, and a reduction of sea ice cover in the Southern Ocean. They also cool the deep ocean worldwide. After the flush, the deep ocean warms slowly again and CO₂ is taken up by the ocean until the stratification becomes unstable again at high latitudes thousands of years later. The existence of a threshold in CO₂ concentration which places the UVic ESCM in either an oscillating or non-oscillating state makes our results intriguing. If the UVic ESCM captures a mechanism that is present and important in the real climate system, the consequences would comprise a rapid increase in atmospheric carbon dioxide concentrations of several tens of ppm, an increase in global surface temperature of the order of 1–2°C, local temperature changes of the order of 6°C and a profound change in ocean stratification, deep water temperature and sea ice cover.