Ocean dynamics determine the response of oceanic CO2 uptake to climate change

The increase of atmospheric CO₂ concentrations due to anthropogenic activities is substantially damped by the ocean, whose CO₂ uptake is determined by the state of the ocean, which in turn is influenced by climate change. We investigate the mechanisms of the ocean’s carbon uptake within the feedback loop of atmospheric CO₂ concentration, climate change and atmosphere/ocean CO₂ flux. We evaluate two transient simulations from 1860 until 2100, performed with a version of the Max Planck Institute Earth System Model (MPI-ESM) with the carbon cycle included. In both experiments observed anthropogenic CO₂ emissions were prescribed until 2000, followed by the emissions according to the IPCC Scenario A2. In one simulation the radiative forcing of changing atmospheric CO₂ is taken into account (coupled), in the other it is suppressed (uncoupled). In both simulations, the oceanic carbon uptake increases from 1 GT C/year in 1960 to 4.5 GT C/year in 2070. Afterwards, this trend weakens in the coupled simulation, leading to a reduced uptake rate of 10% in 2100 compared to the uncoupled simulation. This includes a partial offset due to higher atmospheric CO₂ concentrations in the coupled simulation owing to reduced carbon uptake by the terrestrial biosphere. The difference of the oceanic carbon uptake between both simulations is primarily due to partial pressure difference and secondary to solubility changes. These contributions are widely offset by changes of gas transfer velocity due to sea ice melting and wind changes. The major differences appear in the Southern Ocean (?45%) and in the North Atlantic (?30%), related to reduced vertical mixing and North Atlantic meridional overturning circulation, respectively. In the polar areas, sea ice melting induces additional CO₂ uptake (+20%).     

Can ocean iron fertilization mitigate ocean acidification?

Ocean iron fertilization has been proposed as a method to mitigate anthropogenic climate change, and there is continued commercial interest in using iron fertilization to generate carbon credits. It has been further speculated that ocean iron fertilization could help mitigate ocean acidification. Here, using a global ocean carbon cycle model, we performed idealized ocean iron fertilization simulations to place an upper bound on the effect of iron fertilization on atmospheric CO₂ and ocean acidification. Under the IPCC A2 CO₂ emission scenario, at year 2100 the model simulates an atmospheric CO₂ concentration of 965 ppm with the mean surface ocean pH 0.44 units less than its pre-industrial value of 8.18. A globally sustained ocean iron fertilization could not diminish CO₂ concentrations below 833 ppm or reduce the mean surface ocean pH change to less than 0.38 units. This maximum of 0.06 unit mitigation in surface pH change by the end of this century is achieved at the cost of storing more anthropogenic CO₂ in the ocean interior, furthering acidifying the deep-ocean. If the amount of net carbon storage in the deep ocean by iron fertilization produces an equivalent amount of emission credits, ocean iron fertilization further acidifies the deep ocean without conferring any chemical benefit to the surface ocean.     

The effect of atmospheric CO2 and ice sheet topography on LGM climate

The role of reduced atmospheric CO₂ concentration and ice sheet topography plus its associated land albedo on the LGM climate is investigated using a coupled atmosphere-ocean-sea ice climate system model. The surface cooling induced by the reduced CO₂ concentration is larger than that by the ice sheet topography plus other factors by about 30% for the surface air temperature and by about 100% for the sea surface temperature. A large inter-hemispheric asymmetry in surface cooling with a larger cooling in the Northern Hemisphere is found for both cases. This asymmetric inter-hemispheric temperature response is consistent in the ice sheet topography case with earlier studies using an atmospheric model coupled with a mixed-layer ocean representation, but contrasts with these results in the reduced CO₂ case. The incorporation of ocean dynamics presumably leads to a larger snow and sea ice feedback as a result of the reduction in northward ocean heat transport, mainly as a consequence of the decrease in the North Atlantic overturning circulation by the substantial freshening of the North Atlantic convection regions. A reversed case is found in the Southern Ocean. Overall, the reduction in atmospheric CO₂ concentration accounts for about 60% of the total LGM climate change.     

Wind-driven changes in Southern Ocean residual circulation, ocean carbon reservoirs and atmospheric CO2

The effect of idealized wind-driven circulation changes in the Southern Ocean on atmospheric CO₂ and the ocean carbon inventory is investigated using a suite of coarse-resolution, global coupled ocean circulation and biogeochemistry experiments with parameterized eddy activity and only modest changes in surface buoyancy forcing, each experiment integrated for 5,000 years. A positive correlation is obtained between the meridional overturning or residual circulation in the Southern Ocean and atmospheric CO₂: stronger or northward-shifted westerly winds in the Southern Hemisphere result in increased residual circulation, greater upwelling of carbon-rich deep waters and oceanic outgassing, which increases atmospheric pCO₂ by ～20 μatm; weaker or southward-shifted winds lead to the opposing result. The ocean carbon inventory in our model varies through contrasting changes in the saturated, disequilibrium and biogenic (soft-tissue and carbonate) reservoirs, each varying by O(10–100) PgC, all of which contribute to the net anomaly in atmospheric CO₂. Increased residual overturning deepens the global pycnocline, warming the upper ocean and decreasing the saturated carbon reservoir. Increased upwelling of carbon- and nutrient-rich deep waters and inefficient biological activity results in subduction of unutilized nutrients into the ocean interior, decreasing the biogenic carbon reservoir of intermediate and mode waters ventilating the Northern Hemisphere, and making the disequilibrium carbon reservoir more positive in the mode waters due to the reduced residence time at the surface. Wind-induced changes in the model carbon inventory are dominated by the response of the global pycnocline, although there is an additional abyssal response when the peak westerly winds change their latitude, altering their proximity to Drake Passage and changing the depth extent of the southward return flow of the overturning: a northward shift of the westerly winds isolates dense isopycnals, allowing biogenic carbon to accumulate in the deep ocean of the Southern Hemisphere, while a southward shift shoals dense isopycnals that outcrop in the Southern Ocean and reduces the biogenic carbon store in the deep ocean.     

Simulation of sea ice and ocean variability in the arctic during 1955–2002 with an intermediate complexity model

Abstract

A number of recent sea‐ice and ocean changes in the Arctic and subarctic regions are simulated using the global University of Victoria (UVic) Earth System Climate Model version 2.6. This is an intermediate complexity model which includes a three‐dimensional ocean model (MOM 2.2), an energy‐moisture balance model for the atmosphere with heat and moisture transport, and a dynamic‐thermodynamic sea‐ice model with elastic‐viscous‐plastic rheology. The model is first spun up for 1800 years with monthly wind stress forcing derived from the National Centers for Environmental Prediction (NCEP) climatology winds and a pre‐industrial atmospheric CO₂ concentration of 280 ppm. After a second spin‐up for the period 1800–1947 with daily climatology winds‐tress forcing, and a linearly increasing atmospheric CO₂ concentration, the model is run with interannually varying wind stresses for the period 1948–2002 with an average forcing interval of 2.5 days and an exponentially increasing atmospheric CO₂ concentration varying from 315 to 365 ppm. However, the analysis of the model output is only carried out for the years 1955–2002.

The simulated maximum and minimum sea‐ice areas for the Arctic are within 6% of the observed climatologies for the years 1978–2001. The model output also shows a small downward trend in sea‐ice extent, which, however, is smaller than has been observed during the past few decades. In addition, the model simulates a decrease in sea‐ice thickness in the SCICEX (SCientific ICe EXpeditions) measurement area in the central Arctic that is consistent with, but smaller than, that observed from submarine sonar profiling data.

The observed variability and magnitude of the export of sea ice through Fram Strait is quite well captured in the simulation. The change in correlation between the North Atlantic Oscillation (NAO) index and the sea‐ice export around 1977 as found in a data study by Hilmer and Jung (2000) is also reproduced. Within the Arctic basin the model simulates well the patterns and the timing of the two major regimes of wind‐forced sea‐ice drift circulation (cyclonic and anticyclonic) as found earlier by Proshutinsky and Johnson (1997). The influence of variations in the Fram Strait ice export on the strength of the North Atlantic thermohaline circulation and surface air temperature are also determined. In particular, it is shown that 3–4 years after a large ice export, the maximum meridional overturning streamfunction decreases by more than 10%.

The temperature and salinity increase at depths of 200–300 m, as observed in the eastern Arctic by Morison et al. (1998), between the USS Pargo cruise in 1993 and the Environmental Working Group (EWG) Joint USRussian Arctic Atlas climatology for the years 1948–87, are just visible in the model simulation. The increases are more noticeable, however, when the ocean model data are averaged over the pentade 1995–2000 and compared with model data averaged over the pentade 1955–60. The fact that these, and some of the other modelled changes, are smaller than the observed changes can likely be attributed to the relatively coarse resolution of the UVic Earth System Climate Model (3.6°E‐W and 1.8°N‐S). Nevertheless, the fact that the model captures qualitatively many of the recent sea‐ice and ocean changes in the Arctic suggests that it can be successfully used to investigate other Arctic‐North Atlantic Ocean climate interactions during past and future eras.     

Snow cover validation and sensitivity to CO2 in the UVic ESCM

Abstract

The University of Victoria's (UVic) Earth System Climate model is used to conduct equilibrium atmospheric CO₂ sensitivity experiments over the range 200–1600 ppm in order to explore changes in northern hemisphere snow cover and feedbacks on terrestrial surface air temperature (SAT). Simulations of warmer climates predict a retreat of snow cover over northern continents, in a northeasterly direction. The decline in northern hemisphere global snow mass is estimated to reach 33% at 600 ppm and 54% at 1200 ppm. In the most northerly regions, annual mean snow depth increases for simulations with CO₂ levels higher than present day. The shift in the latitude of maximum snowfall is estimated to be inversely proportional to the CO₂ concentration. The northern hemisphere net shortwave radiation changes are found to be greater over land than over the ocean, suggesting a stronger albedo feedback from changes in terrestrial snow cover than from changes in sea ice. Results also reveal high sensitivity of the snow mass balance under low CO₂ conditions. The amplification feedback (defined as the zonal SAT anomaly caused by doubling CO₂ divided by the equatorial anomaly) is greatest for scenarios with less than 300 ppm, reaching 1.9 at the pole for 250 ppm. The stronger feedback is attributed to the significant albedo changes over land areas. The simulation with 200 ppm triggers continuous accumulation of snow ('glaciation') in regions which, according to paleo‐reconstructions, were covered by ice during the last glacial cycle (the Canadian Arctic, Scandinavia, and the Taymir Peninsula).     

Response of Southern Ocean circulation to global warming may enhance basal ice shelf melting around Antarctica

We investigate the large-scale oceanic features determining the future ice shelf–ocean interaction by analyzing global warming experiments in a coarse resolution climate model with a comprehensive ocean component. Heat and freshwater fluxes from basal ice shelf melting (ISM) are parameterized following Beckmann and Goosse [Ocean Model 5(2):157–170, 2003]. Melting sensitivities to the oceanic temperature outside of the ice shelf cavities are varied from linear to quadratic (Holland et al. in J Clim 21, 2008). In 1% per year CO₂-increase experiments the total freshwater flux from ISM triples to 0.09 Sv in the linear case and more than quadruples to 0.15 Sv in the quadratic case after 140 years at which 4 × 280 ppm = 1,120 ppm was reached. Due to the long response time of subsurface temperature anomalies, ISM thereafter increases drastically, if CO₂ concentrations are kept constant at 1,120 ppm. Varying strength of the Antarctic circumpolar current (ACC) is crucial for ISM increase, because southward advection of heat dominates the warming along the Antarctic coast. On centennial timescales the ACC accelerates due to deep ocean warming north of the current, caused by mixing of heat along isopycnals in the Southern Ocean (SO) outcropping regions. In contrast to previous studies we find an initial weakening of the ACC during the first 150 years of warming. This purely baroclinic effect is due to a freshening in the SO which is consistent with present observations. Comparison with simulations with diagnosed ISM but without its influence on the ocean circulation reveal a number of ISM-related feedbacks, of which a negative ISM-feedback, due to the ISM-related local oceanic cooling, is the dominant one.     

Benchmarking Climate-Carbon Model Simulations against Forest FACE Data

Terrestrial carbon fluxes are an important factor in regulating concentrations of atmospheric carbon dioxide (CO₂). In this study, we use a coupled climate model with interactive biogeochemistry to benchmark the simulation of net primary productivity (NPP) and its response to elevated atmospheric CO₂. Short-term field experiments such as Free-Air Carbon Dioxide Enrichment (FACE) studies have examined this phenomenon but it is difficult to infer trends from only a few years of field data. Here, we employ the University of Victoria's Earth System Climate Model (UVic ESCM) version 2.8 to compare simulated changes in NPP due to an elevated atmospheric CO₂ concentration of 550 ppm to observed increases in NPP of 23% ±2% from four temperate forest FACE studies between 1997 and 2002. We further compare two scenarios: elevated CO₂ with climate change, and elevated CO₂ without climate change, the latter being consistent with FACE methodology. In the climate change scenario global terrestrial and forest-only NPP increased by 24.5% and 27.9%, respectively, while these increases were 21.0% and 17.2%, respectively, in the latitude band most representative of the location of the FACE studies. In the scenario without climate change, terrestrial and forest-only NPP increased instead by 28.3% and 30.6%, respectively, while these increases were 24.3% and 14.4%, respectively, in the FACE latitudes. This suggests that the model may underestimate temperate forest NPP increases when compared to results from temperate forest FACE studies and highlights the need for both increased experimental study of other forest biomes and further model development.     

Simulation of glacial-interglacial cycles by simple relaxation models: consistency with observational results

The model of Paillard and Parrenin (Earth Planet Sci Lett 227:263–271, 2004) was modified to obtain a closer fit to δ¹⁸O and CO₂ time series for the last 800 kyr. The model performance can be improved if its CO₂ sensitivity to I₆₅ insolation is eliminated and if different response times are assumed for ablation/accumulation of ice. Correlations between simulated and experimental time series for CO₂ and ice volume V increase from 0.59 and 0.63 to 0.79 and 0.88, respectively. According to these models, terminations are produced by I₆₅ amplification through CO₂-T and T-CO₂ feedbacks, in synergy with an extra CO₂ contribution from the deep ocean. This contribution is strongly dependent on ice-sheet extent and ice volume (or alternatively, CO₂ concentration, which is a good proxy of Antarctic temperature) but is insensitive to Southern Ocean (SO) insolation on 21 February (I₆₀). Change of deep SO state may be the “order parameter” for nonlinear deglacial changes. According to these models, 100 kyr periodicity of glacial cycles arises from the characteristic time of Antarctic ice sheet advance to the continental slope.     

On the response of the oceanic wind-driven circulation to atmospheric CO2 increase

A global, flux-corrected climate model is employed to predict the surface wind stress and associated wind-driven oceanic circulation for climate states corresponding to a doubling and quadrupling of the atmospheric CO₂ concentration in a simple 1% per year CO₂ increase scenario. The model indicates that in response to CO₂ increase, the position of zero wind stress curl in the mid-latitudes of the Southern Hemisphere shifts poleward. In addition, the wind stress intensifies significantly in the mid-latitudes of the Southern Hemisphere. As a result, the rate of water circulation in the subpolar meridional overturning cell in the Southern Ocean increases by about 6 Sv (1 Sv=10⁶ m³ s⁻¹) for doubled CO₂ and by 12 Sv for quadrupled CO₂, implying an increase of deep water upwelling south of the circumpolar flow and an increase of Ekman pumping north of it. In addition, the changes in the wind stress and wind stress curl translate into changes in the horizontal mass transport, leading to a poleward expansion of the subtropical gyres in both hemispheres, and to strengthening of the Antarctic Circumpolar Current. Finally, the intensified near-surface winds over the Southern Ocean result in a substantial increase of mechanical energy supply to the ocean general circulation.     

The role of land surface dynamics in glacial inception: a study with the UVic Earth System Model

The first results of the UVic Earth System Model coupled to a land surface scheme and a dynamic global vegetation model are presented in this study. In the first part the present day climate simulation is discussed and compared to observations. We then compare a simulation of an ice age inception (forced with 116 ka BP orbital parameters and an atmospheric CO₂ concentration of 240 ppm) with a preindustrial run (present day orbital parameters, atmospheric [CO₂] = 280 ppm). Emphasis is placed on the vegetations response to the combined changes in solar radiation and atmospheric CO₂ level. A southward shift of the northern treeline as well as a global decrease in vegetation carbon is observed in the ice age inception run. In tropical regions, up to 88% of broadleaf trees are replaced by shrubs and C₄ grasses. These changes in vegetation cover have a remarkable effect on the global climate: land related feedbacks double the atmospheric cooling during the ice age inception as well as the reduction of the meridional overturning in the North Atlantic. The introduction of vegetation related feedbacks also increases the surface area with perennial snow significantly.     

Evaluation of ocean and climate models using present‐day observations and forcing

Abstract

The most common method used to evaluate climate models involves spinning them up under perpetual present‐day forcing and comparing the model results with present‐day observations. This approach clearly ignores any potential long‐term memory of the model ocean to past climatic conditions. Here we examine the validity of this approach through the 6000‐year integration of a coupled atmosphere–ocean–sea‐ice model. The coupled model is initially spun‐up with atmospheric CO₂ concentrations and orbital parameters applicable for 6KBP. The model is then integrated forward in time to 2100. Results from this transient coupled model simulation are compared with the results from two additional simulations, in which the model is spun up with perpetual 1850 (preindustrial) and 1998 (present‐day) atmospheric CO₂ concentrations and orbital parameters. This comparison leads to substantial differences between the equilibrium climatologies and the transient simulation, even at 1850 (in weakly ventilated regions), prior to any significant changes in atmospheric CO₂. When compared to the present‐day equilibrium climatology, differences are very large: the global mean surface air and sea surface temperatures are ,0.5°C and ,0.4°C colder, respectively, deep ocean temperatures are substantially cooler, Southern Hemisphere sea‐ice cover is 38% larger, and the North Atlantic conveyor 16% weaker in the transient case. These differences are due to the long timescale memory of the deep ocean to climatic conditions which prevailed throughout the late Holocene, as well as to its large thermal inertia. It is also demonstrated that a ‘cold start’ global warming simulation (one that starts from a 1998 equilibrium climatology) underestimates the global temperature increase at 2100 by ,10%. Our results question the accuracy of current techniques for climate model evaluation and underline the importance of using paleoclimatic simulations in parallel with present‐day simulations in this evaluation process.     

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.     

Estimates of anthropogenic CO<Subscript>2</Subscript> uptake in a global ocean model

A global ocean general circulation model (L30T63) is employed to study the uptake and distribution of anthropogenic CO₂ in the ocean. A subgrid-scale mixing scheme called GM90 is used in the model. There are two main GM90 parameters including isopycnal diffusivity and skew (thickness) diffusivity. Sensitivities of the ocean circulation and the redistribution of dissolved anthropogenic CO₂ to these two parameters are examined. Two runs estimate the global oceanic anthropogenic CO₂ uptake to be 1.64 and 1.73 Pg C yr^-1 for the 1990s, and that the global ocean contained 86.8 and 92.7 Pg C of anthropogenic CO₂ at the end of 1994, respectively. Both the total inventory and uptake from our model are smaller than the data-based estimates. In this presentation, the vertical distributions of anthropogenic CO₂ at three meridional sections are discussed and compared with the available data-based estimates. The inventory in the individual basins is also calculated. Use of large isopycnal diffusivity can generally improve the simulated results, including the exchange flux, the vertical distribution patterns, inventory, storage, etc. In terms of comparison of the vertical distributions and column inventory, we find that the total inventory in the Pacific Ocean obtained from our model is in good agreement with the data-based estimate, but a large difference exists in the Atlantic Ocean, particularly in the South Atlantic. The main reasons are weak vertical mixing and that our model generates small exchange fluxes of anthropogenic CO₂ in the Southern Ocean. Improvement in the simulation of the vertical transport and sea ice in the Southern Ocean is important in future work.     

Modelling the concentration of atmospheric CO2 during the Younger Dryas climate event

 The Younger Dryas (YD, dated between 12.7–11.6 ky BP in the GRIP ice core, Central Greenland) is a distinct cold period in the North Atlantic region during the last deglaciation. A popular, but controversial hypothesis to explain the cooling is a reduction of the Atlantic thermohaline circulation (THC) and associated northward heat flux as triggered by glacial meltwater. Recently, a CH₄-based synchronization of GRIP δ¹⁸O and Byrd CO₂ records (West Antarctica) indicated that the concentration of atmospheric CO₂ (CO^atm ₂) rose steadily during the YD, suggesting a minor influence of the THC on CO^atm ₂ at that time. Here we show that the CO^atm ₂ change in a zonally averaged, circulation-biogeochemistry ocean model when THC is collapsed by freshwater flux anomaly is consistent with the Byrd record. Cooling in the North Atlantic has a small effect on CO^atm ₂ in this model, because it is spatially limited and compensated by far-field changes such as a warming in the Southern Ocean. The modelled Southern Ocean warming is in agreement with the anti-phase evolution of isotopic temperature records from GRIP (Northern Hemisphere) and from Byrd and Vostok (East Antarctica) during the YD. δ¹³C depletion and PO₄ enrichment are predicted at depth in the North Atlantic, but not in the Southern Ocean. This could explain a part of the controversy about the intensity of the THC during the YD. Potential weaknesses in our interpretation of the Byrd CO₂ record in terms of THC changes are discussed. Received: 27 May 1998 / Accepted: 5 November 1998     

Declining Temporal Effectiveness of Carbon Sequestration: Implications for Compliance with the United National Framework Convention on Climate Change

Carbon sequestration is increasingly being promoted as a potential response to the risks of unrestrained emissions of CO₂, either in place of or as a complement to reductions in the use of fossil fuels. However, the potential role of carbon sequestration as an (at-least partial) substitute for reductions in fossil fuel use can be properly evaluated only in the context of a long-term acceptable limit (or range of limits) to the increase in atmospheric CO₂ concentration, taking into account the response of the entire carbon cycle to artificial sequestration. Under highly stringent emission-reduction scenarios for non-CO₂ greenhouse gases, 450 ppmv CO₂ is the equivalent, in terms of radiative forcing of climate,to a doubling of the pre-industrial concentration of CO₂. It is argued in this paper that compliance with the United Nations Framework Convention on Climate Change (henceforth, the UNFCCC) implies that atmospheric CO₂ concentration should be limited, or quickly returned to, a concentration somewhere below 450 ppmv. A quasi-one-dimensional coupled climate-carbon cycle model is used to assess the response of the carbon cycle to idealized carbon sequestration scenarios. The impact on atmospheric CO₂ concentration of sequestering a given amount of CO₂ that would otherwise be emitted to the atmosphere, either in deep geological formations or in the deep ocean, rapidly decreases over time. This occurs as a result of a reduction in the rate of absorption of atmospheric CO₂ by the natural carbon sinks (the terrestrial biosphere and oceans) in response to the slower buildup of atmospheric CO₂ resulting from carbon sequestration. For 100 years of continuous carbon sequestration, the sequestration fraction (defined as the reduction in atmospheric CO₂ divided by the cumulative sequestration) decreases to 14% 1000 years after the beginning of sequestration in geological formations with no leakage, and to 6% 1000 years after the beginning of sequestration in the deep oceans. The difference (8% of cumulative sequestration) is due to an eflux from the ocean to the atmosphere of some of the carbon injected into the deep ocean.The coupled climate-carbon cycle model is also used to assess the amount of sequestration needed to limit or return the atmospheric CO₂ concentration to 350–400 ppmv after phasing out all use of fossil fuels by no later than 2100. Under such circumstances, sequestration of 1–2 Gt C/yr by the latter part of this century could limit the peak CO₂ concentration to 420–460 ppmv, depending on how rapidly use of fossilfuels is terminated and the strength of positive climate-carbon cycle feedbacks. To draw down the atmospheric CO₂ concentration requires creating negative emissions through sequestration of CO₂ released as a byproduct of the production of gaseous fuels from biomass primary energy. Even if fossil fuel emissions fall to zero by 2100, it will be difficult to create a large enough negative emission using biomass energy to return atmospheric CO₂ to 350 ppmv within 100 years of its peak. However, building up soil carbon could help in returning CO₂ to 350 ppmv within 100 years of its peak. In any case, a 100-year period of climate corresponding to the equivalent of a doubled-CO₂ concentration would occur before temperatures decreased. Nevertheless, returning the atmospheric CO₂concentration to 350 ppmv would reduce longterm sea level rise due to thermal expansion and might be sufficient to prevent the irreversible total melting of the Greenland ice sheet, collapse of the West Antarctic ice sheet, and abrupt changes in ocean circulation that might otherwise occur given a prolonged doubled-CO₂ climate. Recovery of coral reef ecosystems, if not already driven to extinction, could begin.     

Role of the ocean in controlling atmospheric CO2 concentration in the course of global glaciations

Responses of ocean circulation and ocean carbon cycle in the course of a global glaciation from the present Earth conditions are investigated by using a coupled climate-biogeochemical model. We investigate steady states of the climate system under colder conditions induced by a reduction of solar constant from the present condition. A globally ice-covered solution is obtained under the solar constant of 92.2% of the present value. We found that because almost all of sea water reaches the frozen point, the ocean stratification is maintained not by temperature but by salinity just before the global glaciation (at the solar constant of 92.3%). It is demonstrated that the ocean circulation is driven not by the surface cooling but by the surface freshwater forcing associated with formation and melting of sea ice. As a result, the deep ocean is ventilated exclusively by deep water formation in southern high latitudes where sea ice production takes place much more massively than northern high latitudes. We also found that atmospheric CO₂ concentration decreases through the ocean carbon cycle. This reduction is explained primarily by an increase of solubility of CO₂ due to a decrease of sea surface temperature, whereas the export production weakens by 30% just before the global glaciation. In order to investigate the conditions for the atmospheric CO₂ reduction to cause global glaciations, we also conduct a series of simulations in which the total amount of carbon in the atmosphere?Cocean system is reduced from the present condition. Under the present solar constant, the results show that the global glaciation takes place when the total carbon decreases to be 70% of the present-day value. Just before the glaciation, weathering rate becomes very small (almost 10% of the present value) and the organic carbon burial declines due to weakened biological productivity. Therefore, outgoing carbon flux from the atmosphere?Cocean system significantly decreases. This suggests the atmosphere?Cocean system has strong negative feedback loops against decline of the total carbon content. The results obtained here imply that some processes outside the atmosphere?Cocean feedback loops may be required to cause global glaciations.     

Modelling the ocean response to secular surface temperature transitions

The ocean response to surface temperature transients is simulated with the use of the Hamburg large-scale geostrophic (LSG) ocean general circulation model (OGCM). The transition, from the present to a climate corresponding to a doubling of the atmospheric CO₂ content, is compared with the reversed transition. For the Atlantic, the time scale for the deep ocean to adjust to the temperature changes was similar for both transitions. In the Pacific, the time scale is shorter for the present to warm transition than for the reverse case, a result of increased production of Antarctic bottom water (AABW) during the warm climate. While the transition from cold to warm climate shows no secular variability, the reversed transition generates considerable variability on time scales of 300–400 years. For the warm climate, oscillations with periods of 45 years are found in the Southern Ocean. Results of principal oscillation pattern (POP) analysis indicate that these oscillations are due to interaction between convection in the Southern Ocean and advected salinity anomalies in the Antarctic Circumpolar Current (ACC) and the Southern Pacific Ocean. Received: 19 September 1995 / Accepted: 15 March 1996     

On the multi-century Southern Hemisphere response to changes in atmospheric CO2-concentration in a Global Climate Model

Summary The transient response of the Southern Hemisphere to climate change is examined using an intermediate complexity climate model. Unlike previous studies, the Southern Ocean response on the centennial to multi-centennial time-scale is assessed in some detail. It is shown that changes in atmospheric CO₂-concentrations lead to an increase in the strength of the Antarctic Circumpolar Current (ACC) by ∼20 Sv by 2750 for an atmospheric CO₂-concentration of 750 ppm. This increase is predominantly the result of an enhanced steric height gradient. The increase in the strength of the ACC induces changes in its steering around topographic features. This change in ACC pathway causes increased surface flow of colder waters into some regions (reducing the rate of warming) and increased surface flow of warmer waters into others (increasing the rate of warming). This meridional shifting of the ACC causes changes in atmospheric temperature in the Southern Hemisphere to be nonuniform. It is also shown that the strength and location of the Antarctic Bottom Water (AABW) overturning cell is affected by increased atmospheric CO₂. For a CO₂-concentration scenario increasing gradually to 750 ppm, AABW production initially decreases, then recovers and eventually increases. New production zones form, which extend AABW production all the way from the Weddell Sea eastward into the Ross Sea. These new production zones are the result of increased areas of atmosphere-ocean interactions, due to decreased sea-ice coverage, although the overturned waters are now warmer and fresher due to climate change. A new production zone of Antarctic Intermediate water is also established in the Southeast Pacific Ocean, poleward of its present-day location.     

General circulation model CO2 sensitivity experiments: Snow-sea ice albedo parameterizations and globally averaged surface air temperature

Two experiments are performed with the NCAR Community Climate Model (CCM) coupled to a swamp ocean with annually averaged solar forcing. A swamp ocean model is one in which the ocean temperature is computed from a surface energy balance. Both experiments are run with present (1 × CO₂) and doubled (2 × CO₂) amounts of atmospheric carbon dioxide (CO₂). The first tests the sensitivity of the model to a snow and sea-ice-albedo formulation which facilitates relatively greater ice melt. The second assesses the model response when the basic state of the model in the control run is colder due to a 2% decrease in solar constant. Both are compared to a previous experiment with the same model using a different snow and sea-ice-albedo formulation and the present value of the solar constant. It is found that the globally averaged surface air temperature increase due to a doubling of CO₂ is highly dependent on (1) the type of snow-sea-ice-albedo formulation used such that the parameterization which better facilitates relatively greater ice melt exhibits a greater sensitivity to increased CO₂, and (2) the basic state of the control run such that the colder the basic state, the greater the warming due to increased CO₂.A portion of this study is supported by the U.S. Department of Energy as part of its Carbon Dioxide Research Program.The National Center for Atmospheric Research is sponsored by the National Science Foundation