Long-term climate changes due to increased CO2 concentration in the coupled atmosphere-ocean general circulation model ECHAM3/LSG

 The long-term adjustment processes of atmosphere and ocean in response to gradually increased atmospheric CO₂ concentration have been analysed in two 850-year integrations with a coupled atmosphere-ocean general circulation model (AOGCM). In these experiments the CO₂ concentration has been increased to double and four times the initial concentration, respectively, and is kept fixed thereafter. Three characteristic time scales have been identified: a very fast response associated with processes dominated by the atmospheric adjustment, an intermediate time scale of a few decades connected with processes in the upper ocean, and adjustment processes with time scales of centuries and longer due to the inertia of the deep ocean. The latter in particular is responsible for a still ongoing adjustment of the atmosphere-ocean system at the end of the integrations after 850 years. After 60 years, at the time of CO₂ doubling, the global mean near-surface air temperature rises by 1.4 K. In spite of the constant CO₂ concentration during the following centuries the warming continues to 2.6 K after 850 years. The behaviour of the quadrupling run is similar: global mean near-surface air temperature increases by 3.8 K at the time of CO₂ quadrupling and by 4.8 K at the end of the simulation. The thermohaline circulation undergoes remarkable changes. Temporarily, the North Atlantic overturning circulation weakens by up to 30% in the CO₂ doubling experiment and up to 50% in the CO₂ quadrupling experiment. After reaching the minimum the North Atlantic overturning slowly recovers in both experiments. Received: 23 August 1999 / Accepted: 27 April 2000     

On the sensitivity of global warming experiments to the parametrisation of sub-grid scale ocean mixing

An ocean general circulation model coupled to an energy-moisture balance atmosphere model is used to investigate the sensitivity of global warming experiments to the parametrisation of sub-grid scale ocean mixing. The climate sensitivity of the coupled model using three different parametrisations of sub-grid scale mixing is 3°C for a doubling of CO₂ (6°C for a quadrupling of CO₂). This suggests that the ocean has only a weak feedback on global mean surface air temperature although significant regional differences, notably at high latitudes, exist with different sub-grid scale parametrisations. In the experiment using the Gent and McWilliams parametrisation for mixing associated with mesoscale eddies, an enhancement of the surface response in the Southern Ocean is found. This enhancement is largely due to the existence of more realistic sea-ice in the climatological control integration and the subsequent enhanced ice-albedo feedback upon warming. In accordance with earlier analyses, the Gent and McWilliams scheme decreases the global efficiency of ocean heat uptake. During the transient phase of all experiments, the North Atlantic overturning initially weakened but ultimately recovered, surpassing its former strength. This suggests that in the region around the North Atlantic the ocean acts as a negative feedback on local warming during the transient phase but a positive feedback at equilibrium. During the transient phase of the experiments with a more sophisticated and realistic parametrisation of sub-grid scale mixing, warmed Atlantic water was found to penetrate at depth into the Arctic, consistent with recent observations in the region. Received: 14 October 1998 / Accepted: 27 April 1999     

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 CO₂ 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 CO₂, 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, CO₂, 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, CO₂ 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.     

Sea‐ice dynamics and CO2 sensitivity in a global climate model

Abstract

Present‐day results and CO₂ sensitivity are described for two versions of a global climate model (genesis) with and without sea‐ice dynamics. Sea‐ice dynamics is modelled using the cavitating‐fluid method of Flato and Hibler (1990, 1992). The atmospheric general circulation model originated from the NCAR Community Climate Model version 1, but is heavily modified to include new treatments of clouds, penetrative convection, planetary boundary‐layer mixing, solar radiation, the diurnal cycle and the semi‐Lagrangian transport of water vapour. The surface models include an explicit model of vegetation (similar to BATS and SiB), multilayer models of soil, snow and sea ice, and a slab ocean mixed layer.

When sea‐ice dynamics is turned off, the CO₂‐induced warming increases drastically around ～60–80°S in winter and spring. This is due to the much greater (and unrealistic) compactness of the Antarctic ice cover without dynamics, which is reduced considerably when CO₂ is doubled and exposes more open ocean to the atmosphere. With dynamics, the winter ice is already quite dispersed for 1 × CO₂ so that its compactness does not decrease as much when CO₂ is doubled.     

Climate Change to the End of the Millennium

Anthropogenic climate change will continue long after anthropogenic CO₂ emissions cease. Atmospheric CO₂, global warming and ocean circulation will approach equilibrium on the millennial timescale, whereas thermal expansion of the ocean, ice sheet melt and their contributions to sea level rise are unlikely to be complete. Atmospheric CO₂ in year 3000 depends non-linearly on the total amount of CO₂ emitted and is very likely to exceed the present level of ∼380 ppmv. CO₂ is doubled for ∼2500 GtC emitted, quadrupled if all ∼5000 GtC of conventional fossil fuel resources are emitted, and increases by a factor of ∼32 if a further 20,000 GtC of exotic fossil fuel resources are emitted. Global warming in year 3000 will also depend on climate sensitivity to doubling CO₂, which is most probably ∼3 ^∘C but highly uncertain. Thermal expansion will contribute 0.5–2 m to millennial sea level rise for each doubling of CO₂. The Greenland ice sheet could melt completely within the millennium under > 8×CO₂, adding a further ∼7 m to sea level. The rate of melt depends on the magnitude of forcing above a regional warming threshold of 1–3 ^∘C. The West Antarctic ice sheet could be threatened by 4–10 ^∘C local warming, and its potential contribution to millennial sea level rise exceeds current maximum estimates of ∼1 m. The fate of the ocean thermohaline circulation may depend on the rate as well as the magnitude of forcing.     

Multi-Model Projection of July–August Climate Extreme Changes over China under CO_2 Doubling. Part II: Temperature

This is the second part of the authors’ analysis on the output of 24 coupled climate models from the Twentieth-Century Climate in Coupled Models (20C3M) experiment and 1% per year CO 2 increase experiment (to doubling) (1pctto2x) of phase 3 of the Coupled Model Inter-comparison Project (CMIP3). The study focuses on the potential changes of July–August temperature extremes over China. The pattern correlation coefficients of the simulated temperature with the observations are 0.6–0.9, which are higher than the results for precipitation. However, most models have cold bias compared to observation, with a larger cold bias over western China (>5°C) than over eastern China (<2°C). The multi-model ensemble (MME) exhibits a significant increase of temperature under the 1pctto2x scenario. The amplitude of the MME warming shows a northwest–southeast decreasing gradient. The warming spread among the models (～1°C– 2°C) is less than MME warming (～2°C–4°C), indicating a relatively robust temperature change under CO 2 doubling. Further analysis of Geophysical Fluid Dynamics Laboratory coupled climate model version 2.1 (GFDL-CM2.1) simulations suggests that the warming pattern may be related to heat transport by summer monsoons. The contrast of cloud effects also has contributions. The different vertical structures of warming over northwestern China and southeastern China may be attributed to the different natures of vertical circulations. The deep, moist convection over southeastern China is an effective mechanism for "transporting" the warming upward, leading to more upper-level warming. In northwestern China, the warming is more surface-orientated, possibly due to the shallow, dry convection.     

Estimated bioresource sensitivity to climatic change in Alberta,Canada

An index of dry-matter productivity is used to assess the sensitivity of bioresources in Alberta, western Canada, to changes in the thermal and precipitation regimes, particularly to climatic warming. Results suggest that warming would improve productivity in northern Alberta, but reduce it, because of the associated increased moisture stress, in the relatively warm, dry southeastern part of the province. Estimated productivity generally increased with precipitation regardless of location or temperature. Warming induced by CO₂ quadrupling would apparently give a net increase of about 18% in bioresource productivity for the province, and CO₂ doubling would lead to a 16% increase. However, the bioclimate would be changed much more than this might suggest. For CO₂ quadrupling, Alberta would acquire a thermal regime similar to that of present-day Nebraska, some 2,000 km to the southeast. Also, the increase in productivity of plants due to the effects of CO₂ on photosynthesis might be much larger than the climate-related effects. Climate impact assessment in Canada is made especially challenging by the shortness of the period of instrumental record, the relatively high degree of sensitivity to climatic change, and the sparseness of the station network, particularly in the most sensitive areas.     

Climate Change on Newfoundland and Labrador Shelves: Results From a Regional Downscaled Ocean and Sea-Ice Model Under an A1B Forcing Scenario 2011–2069

Climate change may affect ocean and ice conditions in coastal oceans and thus have significant impacts on coastal infrastructure, marine navigation, and marine ecosystems. In this study a three-dimensional ice–ocean model is developed to examine likely changes of ocean and ice conditions over the Newfoundland and Labrador Shelves in response to climate change. The model is configured with a horizontal grid of approximately 7?km and a vertical grid of 46 levels and is run from 1979 to 2069. The projection period is 2011 to 2069 under a median emission scenario A1B used by the Intergovernmental Panel on Climate Change. For the projection period, the surface atmospheric forcing fields used are from the Canadian Regional Climate Model over the North Atlantic. The open boundary conditions come from the Canadian Global Climate Model, Version 3 (CGCM3), adjusted for the 1981–2010 mean of the Simple Ocean Data Assimilation model output. The simulated fields over the 1981–2010 period have patterns consistent with observations. Over the Newfoundland and Labrador Shelves during the projection period, the model shows general trends of warming, freshening, and decreasing ice. From 2011 to 2069, the model projects that under A1B sea surface temperature will increase by 1.4°C; bottom temperature will increase by 1.6°C; sea surface salinity will decrease by 0.7; bottom salinity will decrease by 0.3; and sea-ice extent will decrease by 70%. The sea level will rise by 0.11?m at the St. John's tide-gauge station because of oceanographic change, and the freshwater transport of the Labrador Current will double as a result of freshening. The regional ice–ocean model reproduces more realistic present climate conditions and projects considerably different future climate conditions than CGCM3.     

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.     

Sensitivity of the Humboldt Current system to global warming: a downscaling experiment of the IPSL-CM4 model

The impact of climate warming on the seasonal variability of the Humboldt Current system ocean dynamics is investigated. The IPSL-CM4 large scale ocean circulation resulting from two contrasted climate scenarios, the so-called Preindustrial and quadrupling CO₂, are downscaled using an eddy-resolving regional ocean circulation model. The intense surface heating by the atmosphere in the quadrupling CO₂ scenario leads to a strong increase of the surface density stratification, a thinner coastal jet, an enhanced Peru–Chile undercurrent, and an intensification of nearshore turbulence. Upwelling rates respond quasi-linearly to the change in wind stress associated with anthropogenic forcing, and show a moderate decrease in summer off Peru and a strong increase off Chile. Results from sensitivity experiments show that a 50% wind stress increase does not compensate for the surface warming resulting from heat flux forcing and that the associated mesoscale turbulence increase is a robust feature.     

Transient responses of the CSIRO climate model to two different rates of CO2 increase

 The responses of the CSIRO coupled atmosphere-ocean-sea ice model to two greenhouse gas induced warming scenarios are described and compared to a control run with the current CO₂ level. In one scenario, denoted IS92a, the atmospheric CO₂ increases such that it reaches doubling after 128 years. In the other, the CO₂ increases at 1% per year compounding (doubling after 70 y). As the CO₂ increases in both scenarios, the top-of-atmosphere outgoing longwave radiation increases giving enhanced cooling of the coupled system, while the outgoing short wave radiation decreases contributing to a warming of the system. The latter overcompensates the former leading to a global mean net radiative heat gain. The distribution of this heat gain produces the well-known interhemispheric asymmetry in warming, despite a decrease in the sea ice around Antarctica in this model. It is found that the volume mean temperature response over the southern ocean is greater than that over the northern hemispheric oceans, and a maximum warming takes place at the subsurface rather at the surface of the ocean in the southern mid-to-high latitude region. The enhanced high-latitude freshening associated with the strengthened hydrological cycle significantly affects the latitudinal distribution of warming and other responses. It enhances the warming immediately equatorward of the deep water formation regions while produces a reduced warming, even a cooling, in these regions. In both runs, there is a decrease in the large-scale oceanic currents which have a significant thermohaline-driven component. The reduction in these currents reduces the poleward transport of salt out of the tropical and subtropical regions of these oceans. This and the enhanced evaporation contribute to considerable increases in surface salinity in the tropical and subtropical regions. In IS92a, the warming rate before doubling is smaller than that in 1% scenario, but the cumulative effects of the two experiments at the time of doubling are similar. Nevertheless, significant contrasts exist. For example, at the time of doubling in IS92a, the warming of the upper ocean is greater because a more developed temperature-albedo feedback occurs. In addition, a longer time is allowed for heat anomalies to spread downward, and so the effective heat penetration depth is greater than that in the 1% scenario. Thus the oceanic response is influenced by the CO₂ increase scenario used. Received: 2 September 1997 / Accepted: 21 January 1998     

Deep-ocean heat uptake and equilibrium climate response

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

Response of the Subtropical Gyre Circulation in the North Pacific Ocean to CO2 Quadrupling

ABSTRACT

This study investigates the response of the subtropical gyre circulation in the North Pacific Ocean to quadrupled CO₂ using the Community Earth System Model, version 1 (CESM1). In particular, an overriding technique is applied to isolate and quantify the effects of wind stress and thermal warming caused by CO₂ emissions. Results show that, in response to the increase in CO₂, the total mass transport in the subtropical gyre is reduced by approximately 11%. This reduction results mainly from negative anomalies of the wind stress curl over the subtropical region, with a smaller contribution from the thermal warming effect. Furthermore, a detailed analysis finds that the change in the subtropical gyre is baroclinic in nature [i.e., the gyre appears to be spin-up in the upper ocean (above 300?m) but spin-down in the lower thermocline (from 300 to 1500?m)]. This reversal between the upper ocean and lower thermocline is a result of the thermal warming effect, which intensifies ocean stratification, hindering the transfer of momentum from the upper layers to the lower layers and leading to an acceleration of the gyre in the upper ocean but a deceleration in the lower thermocline. Another feature of the response of the subtropical gyre to quadrupled CO₂ is the respective poleward and equatorward movements of its northern and southern boundaries, which is a result of the change in the zero wind stress curl lines.     

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.     

Towards quantifying uncertainty in transient climate change

Ensembles of coupled atmosphere–ocean global circulation model simulations are required to make probabilistic predictions of future climate change. “Perturbed physics” ensembles provide a new approach in which modelling uncertainties are sampled systematically by perturbing uncertain parameters. The aim is to provide a basis for probabilistic predictions in which the impact of prior assumptions and observational constraints can be clearly distinguished. Here we report on the first perturbed physics coupled atmosphere–ocean model ensemble in which poorly constrained atmosphere, land and sea-ice component parameters are varied in the third version of the Hadley Centre model (the variation of ocean parameters will be the subject of future study). Flux adjustments are employed, both to reduce regional sea surface temperature (SST) and salinity biases and also to admit the use of combinations of model parameter values which give non-zero values for the global radiation balance. This improves the extent to which the ensemble provides a credible basis for the quantification of uncertainties in climate change, especially at a regional level. However, this particular implementation of flux-adjustments leads to a weakening of the Atlantic overturning circulation, resulting in the development of biases in SST and sea ice in the North Atlantic and Arctic Oceans. Nevertheless, model versions are produced which are of similar quality to the unperturbed and un-flux-adjusted version. The ensemble is used to simulate pre-industrial conditions and a simple scenario of a 1% per year compounded increase in CO₂. The range of transient climate response (the 20 year averaged global warming at the time of CO₂ doubling) is 1.5–2.6°C, similar to that found in multi-model studies. Measures of global and large scale climate change from the coupled models show simple relationships with associated measures computed from atmosphere-mixed-layer-ocean climate change experiments, suggesting that recent advances in computing the probability density function of climate change under equilibrium conditions using the perturbed physics approach may be extended to the transient case.     

The transport of CO2-induced warming into the ocean: an analysis of simulations by the OSU coupled atmosphere-ocean general circulation model

The OSU global coupled atmosphere-ocean general circulation model has been used to investigate a 2xCO₂-induced climate change. A previous analysis of the simulated 2xCO₂–1xCO₂ temperature differences showed that the CO₂-induced warming penetrated into the ocean and thereby caused a delay in the equilibration of the climate system with an estimatede-folding time of 50–75 years. The objective of the present study is to determine by what pathways and through which physical processes the simulated ocean general circulation produces the penetration of the CO₂-induced warming into the ocean.A global-mean oceanic heat budget analysis shows that the ocean gains heat at a rate of 3 W/m² due to the CO₂ doubling, and that this heat penetrates downward into the ocean predominantly through the reduction in the convective overturning. A zonal-mean oceanic heat budget analysis shows that the surface warming increases from the tropics toward the midlatitudes of both hemispheres and gradually penetrated into the deeper ocean, with a greater penetration in the subtropics and midlatitudes than in the equatorial region. The zonal-mean heat budget analysis also shows that the CO₂-induced warming of the ocean occurs predominantly through the down-ward transport of heat, with the meridional heat flux being only of secondary importance. In the tropics the penetration of the CO₂-induced heating is minimized by the upwelling of cold water. In the subtropics the heating is transported down-ward more readily by the downwelling existing there. In the high latitudes the suppressed convection plays the dominant role in the downward penetration of the CO₂-induced heating. The latter result should be considered as tentative, however, as the ocean component of the coupled model employed a prescribed surface salinity field and did not include the mechanism of brine rejection when sea water freezes into sea ice.     

An Overview of BCC Climate System Model Development and Application for Climate Change Studies

This paper reviews recent progress in the development of the Beijing Climate Center Climate System Model(BCC-CSM) and its four component models(atmosphere,land surface,ocean,and sea ice).Two recent versions are described:BCC-CSM1.1 with coarse resolution(approximately 2.8125°×2.8125°) and BCC-CSM1.1(m) with moderate resolution(approximately 1.125°×1.125°).Both versions are fully coupled climate-carbon cycle models that simulate the global terrestrial and oceanic carbon cycles and include dynamic vegetation.Both models well simulate the concentration and temporal evolution of atmospheric CO_2 during the 20th century with anthropogenic CO2 emissions prescribed.Simulations using these two versions of the BCC-CSM model have been contributed to the Coupled Model Intercomparison Project phase five(CMIP5) in support of the Intergovernmental Panel on Climate Change(IPCC) Fifth Assessment Report(AR5).These simulations are available for use by both national and international communities for investigating global climate change and for future climate projections.Simulations of the 20th century climate using BCC-CSMl.l and BCC-CSMl.l(m) are presented and validated,with particular focus on the spatial pattern and seasonal evolution of precipitation and surface air temperature on global and continental scales.Simulations of climate during the last millennium and projections of climate change during the next century are also presented and discussed.Both BCC-CSMl.l and BCC-CSMl.l(m) perform well when compared with other CMIP5 models.Preliminary analyses indicate that the higher resolution in BCC-CSM1.1(m) improves the simulation of mean climate relative to BCC-CSMl.l,particularly on regional scales.     

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%).     

EMIC Intercomparison Project (EMIP–CO2): comparative analysis of EMIC simulations of climate, and of equilibrium and transient responses to atmospheric CO2 doubling

An intercomparison of eight EMICs (Earth system Models of Intermediate Complexity) is carried out to investigate the variation and scatter in the results of simulating (1) the climate characteristics at the prescribed 280 ppm atmosphere CO₂ concentration, and (2) the equilibrium and transient responses to CO₂ doubling in the atmosphere. The results of the first part of this intercomparison suggest that EMICs are in reasonable agreement with the present-day observational data. The dispersion of the EMIC results by and large falls within the range of results of General Circulation Models (GCMs), which took part in the Atmospheric Model Intercomparison Project (AMIP) and Coupled Model Intercomparison Project, phase 1 (CMIP1). Probable reasons for the observed discrepancies among the EMIC simulations of climate characteristics are analysed. A scenario with gradual increase in CO₂ concentration in the atmosphere (1% per year compounded) during the first 70 years followed by a stabilisation at the 560 ppm level during a period longer than 1,500 years is chosen for the second part of this intercomparison. It appears that the EMIC results for the equilibrium and transient responses to CO₂ doubling are within the range of the corresponding results of GCMs, which participated in the atmosphere-slab ocean model intercomparison project and Coupled Model Intercomparison Project, phase 2 (CMIP2). In particular EMICs show similar temperature and precipitation changes with comparable magnitudes and scatter across the models as found in the GCMs. The largest scatter in the simulated response of precipitation to CO₂ change occurs in the subtropics. Significant differences also appear in the magnitude of sea ice cover reduction. Each of the EMICs participating in the intercomparison exhibits a reduction of the strength of the thermohaline circulation in the North Atlantic under CO₂ doubling, with the maximum decrease occurring between 100 and 300 years after the beginning of the transient experiment. After this transient reduction, whose minimum notably varies from model to model, the strength of the thermohaline circulation increases again in each model, slowly rising back to a new equilibrium.     

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.