The Response of Arctic Sea Ice to Global Change

The sea ice-covered polar oceans have received wider attention recently for two reasons. Firstly, the global conveyor belt circulation of the ocean is believed to be forced in the North and South Atlantic through deep water formation, which to a large degree is controlled by the variations of the sea ice margin and especially by the sea ice export to lower latitudes. Secondly, CO₂ response experiments with coupled climate models show an enhanced warming in polar regions for increased concentrations of atmospheric greenhouse gases. Whether this large response in high latitudes is due to real physical feedback processes or to unrealistic simplifications of the sea ice model component remains to be determined. Coupled climate models generally use thermodynamic sea ice models or sea ice models with oversimplified dynamics schemes. Realistic dynamic-thermodynamic sea ice models are presently implemented only at a few modeling centers. Sensitivity experiments with thermodynamic and dynamic-thermodynamic sea ice models show that the more sophisticated models are less sensitive to perturbations of the atmospheric and oceanic boundary conditions. Because of the importance of the role of sea ice in mediating between atmosphere and ocean an improved representation of sea ice in global climate models is required. This paper discusses present sea ice modeling as well as the sensitivity of the sea ice cover to changes in the atmospheric boundary conditions. These numerical experiments indicate that the sea ice follows a smooth response function: sea ice thickness and export change by 2% of the mean value per 1 Wm^-2 change of the radiative forcing.     

Ice Age mystery: a proposed theory for the cause of long term climate change

Summary ?Analysis of available data shows that the duration of the glacial/interglacial cycle is determined by the time for the ocean to go through one major temperature cycle. At the start of an interglacial period, clear skies with consequent release of CO₂ from the ocean, warms the atmosphere, which in turn eventually warms the ocean to its maximum. Cloudy skies then cause the climate (land and air temperature) to cool and the CO₂ to be reabsorbed to start glaciation preliminaries. The albedo feedback effect of the glacial ice, a relatively warm ocean, which produces enhanced cloud cover, and the increased solubility of CO₂ in cold seawater ensure a long period of glaciation. Glacial periods end when pack ice spreads out on the ocean cooling it until reduced cloud cover once again allows the Sun’s heat, unreflected by cloud cover, to melt the ice and release CO₂ back into the atmosphere. Received May 22, 2002; accepted June 20, 2002     

Vegetation feedback under future global warming

It has been well documented that vegetation plays an important role in the climate system. However, vegetation is typically kept constant when climate models are used to project anthropogenic climate change under a range of emission scenarios in the Intergovernmental Panel on Climate Change (IPCC) Special Report on Emission Scenarios. Here, an atmospheric general circulation model, and an asynchronously coupled system of an atmospheric and an equilibrium terrestrial biosphere model are forced by monthly sea surface temperature and sea ice extent for the periods 2051?C2060 and 2090?C2098 as projected with 17 atmosphere?Cocean general circulation models participating in the IPCC Fourth Assessment Report, and by appropriate atmospheric carbon dioxide concentrations under the A2 emission scenario. The effects of vegetation feedback under future global warming are then investigated. It is found that the simulated composition and distribution of vegetation during 2051?C2060 (2090?C2098) differ greatly from the present, and global vegetation tends to become denser as expressed by a 21% (36%) increase in global mean leaf area index, which is most pronounced at the middle and high northern latitudes. Vegetation feedback has little effect on globally averaged surface temperature. On a regional scale, however, it induces statistically significant changes in surface temperature, in particular over most parts of continental Eurasia east of about 60°E where annual surface temperature is expected to increase by 0.1?C1.0?K, with an average of about 0.4?K for each future period. These changes can mostly be explained by changes in surface albedo resulting from vegetation changes in the context of future global warming.     

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.     

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.     

The effects of rerouting aircraft around the arctic circle on arctic and global climate

Climate data suggest greater warming over the Arctic than lower latitudes, and the most abundant direct source of black carbon and other climate-relevant pollutants over the Arctic is cross-polar flights by international aviation. A relevant question is whether rerouting cross-polar flights to circumnavigate the Arctic Circle reduces or enhances such warming. To study this issue, a model accounting for subgrid exhaust plumes from each individual commercial flight worldwide was used with 2006 global aircraft emission inventories that treated cross-polar flights and flights rerouted around the Arctic Circle (66.56083 °N), respectively. Rerouting increased fuel use by 0.056?% in the global average, mostly right outside the Arctic Circle, but most of the associated black carbon and other emissions were removed faster because they were now over latitudes of greater precipitation and lesser stability. Rerouting also reduced fuel use and emissions within the Arctic Circle by 83?% and delayed pollutant transport to the Arctic. The Arctic reduction in pollutants, particularly of black carbon, decreased Arctic and global temperature and increased Arctic sea ice over 22?years. Although the slight increase in total CO₂ emissions due to rerouting may dampen the benefit of rerouting over more decades, rerouting or even partial rerouting (allowing cross-polar flights during polar night only) may delay the elimination of Arctic sea ice, which will otherwise likely occur within the next 2–3 decades due to global warming in general. Rerouting may increase worldwide fuel plus operational costs by only ~$99 million/yr, 47–55 times less than an estimated 2025 U.S.-alone cost savings due to the global warming reduction from rerouting.     

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

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.     

The Pangean ice age: studies with a coupled climate-ice sheet model

 Application of an ice sheet model developed for the Pleistocene to the extensive Carboniferous glaciation on Gondwana yields an ice sheet which has several features consistent with observations. While complete deglaciation is not achieved without CO₂ changes, the Milankovich-induced fluctuations in ice sheet volume are comparable to Pleistocene glacial/ interglacial signals. This result is shown to hold for a large fraction of physically reasonable parameter space. The model also exhibits multiple equilibria and sharp bifurcations, as infinitesimal changes in the solar constant or precipitation can lead to a qualitatively different climate. The success of the model in predicting ice location in an environment quite different from the Pleistocene provides additional support for the robustness of the basic model physics and suggests that the model can be applied with some confidence to other pre-Pleistocene glaciations. Received: 30 June 1998 / Accepted: 5 January 1999     

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     

Ice and climate modeling: An editorial essay

The growth and decay of ice sheets are driven by forces affecting the seasonal cycles of snowfall and snowmelt. The external forces are likely to be variations in the earth's orbit which cause differences in the solar radiation received. Radiational control of snowmelt is modulated by the seasonal cycles of snow albedo and cloud cover. The effects of orbital changes can be magnified by feedbacks involving atmospheric CO₂ content, ocean temperatures and desert areas. Climate modeling of the causes of the Pleistocene ice ages involves modeling the interactions of all components of the climate system; snow, sea ice, glacier ice, the ocean, the atmosphere, and the solid earth. Such modeling is also necessary for interpreting oxygen isotope records from ice and ocean as paleoclimatic evidence.     

Heat budget of the upper Arctic Ocean under a warming climate

The heat budget of the upper Arctic Ocean is examined in an ensemble of coupled climate models under idealised increasing CO₂ scenarios. All of the experiments show a strong amplification of surface air temperatures but a smaller increase in sea surface temperature than the rest of the world as heat is lost to the atmosphere as the sea-ice cover is reduced. We carry out a heat budget analysis of the Arctic Ocean in an ensemble of model runs to understand the changes that occur as the Arctic becomes ice free in summer. We find that as sea-ice retreats heat is lost from the ocean surface to the atmosphere contributing to the amplification of Arctic surface temperatures. Furthermore, heat is mixed upwards into the mixed layer as a result of increased upper ocean mixing and there is increased advection of heat into the Arctic as the ice edge retreats. Heat lost from the upper Arctic Ocean to the atmosphere is therefore replenished by mixing of warmer water from below and by increased advection of warm water from lower latitudes. The ocean is therefore able to contribute more to Arctic amplification.     

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.     

Relative contribution of feedback processes to Arctic amplification of temperature change in MIROC GCM

The finding that surface warming over the Arctic exceeds that over the rest of the world under global warming is a robust feature among general circulation models (GCMs). While various mechanisms have been proposed, quantifying their relative contributions is an important task in order to understand model behavior. Here we apply a recently proposed feedback analysis technique to an atmosphere–ocean GCM under two and four times CO₂ concentrations which approximately lead to seasonally and annually sea ice-free climates. The contribution of feedbacks to Arctic temperature change is investigated. The surface warming in the Arctic is contributed by albedo, water vapour and large-scale condensation feedbacks and reduced by the evaporative cooling feedback. The surface warming contrast between the Arctic and the global averages (AA) is maintained by albedo and evaporative cooling feedbacks. The latter contributes to AA predominantly by cooling the low latitudes more than the Arctic. Latent heat transport into the Arctic increases and hence evaporative cooling plus large-scale condensation feedback contributes positively to AA. On the other hand, dry-static energy transport into the Arctic decreases and hence dynamical heating feedback contributes negatively to AA. An important contribution is thus made via changes in hydrological cycle and not via the ‘dry’ heat transport process. A larger response near the surface than aloft in the Arctic is maintained by the albedo, water vapour, and dynamical heating feedbacks, in which the albedo and water vapour feedbacks contribute through warming the surface more than aloft, and the dynamical heating feedback contributes by cooling aloft more than the surface. In our experiments, ocean and sea ice dynamics play a secondary role. It is shown that a different level of CO₂ increase introduces a latitudinal and seasonal difference into the feedbacks.     

The steric component of sea level rise associated with enhanced greenhouse warming: a model study

Climate change due to enhanced greenhouse warming has been calculated using the coupled GFDL general circulation model of the atmosphere and ocean. The results of the model for a sustained increase of atmospheric carbon dioxide of 1% per year over a century indicate a marked warming of the upper ocean. Results of the model are used to study the rise in sea level caused by increase in ocean temperatures and associated changes in ocean circulation. Neglecting possible contributions due to changes in the volume of polar ice sheets and mountain glaciers, the model predicts an average rise in sea level of approximately 15 ± 5 cm by the time atmospheric carbon dioxide doubles. Heating anomalies are greatest in subpolar latitudes. This effect leads to a weakening of the ocean thermohaline circulation. Changes in thermohaline circulation redistribute heat within the ocean from high latitudes toward the equator, and cause a more uniform sea level rise than would occur otherwise.     

Future changes in tropical cyclone genesis in fully dynamic ocean- and mixed layer ocean-coupled climate models: a low-resolution model study

The global and regional projected changes in tropical cyclone (TC) genesis due to increased CO₂ concentrations has been investigated through a large-scale TC genesis parameter (convective seasonal genesis parameter, ConvGP) in two perturbed physics ensembles. The ensembles are based on the third generation Hadley Centre atmosphere?Cocean general circulation model with the first ensemble using a coupled fully dynamic ocean (HadCM3) and the second coupled to a simplified mixed layer thermodynamic ocean (HadSM3) both consisting of 17 members. In each ensemble, parameters are identically perturbed to provide a wide range of climate sensitivity whilst retaining a credible present-day climate simulation. It is found, by comparing the ConvGP climatology from reanalysis data with the best track genesis, that it is possible to reproduce the observed genesis distribution. Future changes in the spatial ConvGP distribution are explored with respect to each tropical ocean basin. Whilst there is a similarity in the gross pattern of the ensemble-mean projected ConvGP change between HadCM3 and HadSM3, there is a non-trivial difference in the tropical Pacific Ocean, arising from different patterns of tropical Pacific sea surface temperature change. This indicates that ocean representation can be important for regional scale projections. The quantitative contribution of individual constituent parameters (i.e. vorticity parameter, shear parameter and convective potential) to the projected ConvGP change is estimated. It is found that all three large-scale parameters generally contribute constructively, but with different magnitude, in the regions where a large doubled CO₂ response is found.     

A new tool for evaluating the physics of coupled atmosphere�Cocean variability in nature and in general circulation models

Intermediate models of the coupled tropical atmosphere?Cocean system have been used to illuminate the physics of interannual climate phenomenon such as El Ni?o Southern Oscillation (ENSO) in the tropical Pacific and to explore how the tropics might respond to a forcing such as changing insolation (Milankovitch) or atmospheric carbon dioxide. Importantly, most of the intermediate models are constructed as anomaly models: models that evolve on a prescribed climatological mean state, which is typically prescribed and done so on a rather ad hoc basis. Here we show how the observed climatological mean state fields [ocean currents and upwelling, sea surface temperature (SST) and atmospheric surface winds] can be incorporated into a linearized intermediate model of the tropical coupled atmosphere?Cocean system: called Linear Ocean?CAtmosphere Model (LOAM), it is a linearized version of the Zebiak and Cane model. With realistic, seasonally varying mean state fields, we find that the essential physics of the ENSO mode is very similar to that in the original model and to that in the observations and that the observed mean fields support an ENSO mode that is stable to perturbations. Thus, our results provide further evidence that ENSO is generated and maintained by stochastic (uncoupled) perturbations. The method that we have outlined can be used to assimilate any set of ocean and atmosphere climatological data into the linearized atmosphere?Cocean model. In a companion paper, we apply this same method to incorporate mean field output from two global climate models into the linearised model. We use the latter to diagnose the physics of the leading coupled mode (ENSO) that is supported by the climate models, and to illuminate why the structure and variance in the ENSO mode changes in the models when they are forced by early Holocene and Last Glacial Maximum boundary conditions.     

Parallel climate model (PCM) control and transient simulations

The Department of Energy (DOE) supported Parallel Climate Model (PCM) makes use of the NCAR Community Climate Model (CCM3) and Land Surface Model (LSM) for the atmospheric and land surface components, respectively, the DOE Los Alamos National Laboratory Parallel Ocean Program (POP) for the ocean component, and the Naval Postgraduate School sea-ice model. The PCM executes on several distributed and shared memory computer systems. The coupling method is similar to that used in the NCAR Climate System Model (CSM) in that a flux coupler ties the components together, with interpolations between the different grids of the component models. Flux adjustments are not used in the PCM. The ocean component has 2/3° average horizontal grid spacing with 32 vertical levels and a free surface that allows calculation of sea level changes. Near the equator, the grid spacing is approximately 1/2° in latitude to better capture the ocean equatorial dynamics. The North Pole is rotated over northern North America thus producing resolution smaller than 2/3° in the North Atlantic where the sinking part of the world conveyor circulation largely takes place. Because this ocean model component does not have a computational point at the North Pole, the Arctic Ocean circulation systems are more realistic and similar to the observed. The elastic viscous plastic sea ice model has a grid spacing of 27?km to represent small-scale features such as ice transport through the Canadian Archipelago and the East Greenland current region. Results from a 300?year present-day coupled climate control simulation are presented, as well as for a transient 1% per year compound CO₂ increase experiment which shows a global warming of 1.27?°C for a 10?year average at the doubling point of CO₂ and 2.89?°C at the quadrupling point. There is a gradual warming beyond the doubling and quadrupling points with CO₂ held constant. Globally averaged sea level rise at the time of CO₂ doubling is approximately 7?cm and at the time of quadrupling it is 23?cm. Some of the regional sea level changes are larger and reflect the adjustments in the temperature, salinity, internal ocean dynamics, surface heat flux, and wind stress on the ocean. A 0.5% per year CO₂ increase experiment also was performed showing a global warming of 1.5?°C around the time of CO₂ doubling and a similar warming pattern to the 1% CO₂ per year increase experiment. El Niño and La Niña events in the tropical Pacific show approximately the observed frequency distribution and amplitude, which leads to near observed levels of variability on interannual time scales.     

CO2 feedbacks and the 100 K year cycle in climate

Summary In an earlier paper (Lindzen, 1986), it was shown that allowing CO₂ to vary with snow/sea ice position could lead to a greatly enhanced response in glaciation to 100 K year orbital forcing—even when 20 K year forcing was much stronger. In that model, snow/sea ice position (SSIP) and glaciation were different: the former was the forcing for the latter. However, SSIP and glaciation were not decorrelated. Observations (Berner et al., 1979; Lorius et al., 1985; Neftel et al., 1982) suggest that CO₂ may be independently related to both SSIP and glaciation. In the present paper, we allow (in a highly simplified manner) such independent dependence, and show how it alters the earlier results. Briefly, the dependence of CO₂ on glaciation can contribute to and even cause a highly enhanced response to the 100 K year component of the forcing. However, the CO₂ dependence on SSIP is, on the whole, more effective in this regard. Thus, we expect time series of CO₂ to show variation on the faster time scales than does glaciation.With 5 Figures     

An examination of climate sensitivity for idealised climate change experiments in an intermediate general circulation model

Radiative forcing and climate sensitivity have been widely used as concepts to understand climate change. This work performs climate change experiments with an intermediate general circulation model (IGCM) to examine the robustness of the radiative forcing concept for carbon dioxide and solar constant changes. This IGCM has been specifically developed as a computationally fast model, but one that allows an interaction between physical processes and large-scale dynamics; the model allows many long integrations to be performed relatively quickly. It employs a fast and accurate radiative transfer scheme, as well as simple convection and surface schemes, and a slab ocean, to model the effects of climate change mechanisms on the atmospheric temperatures and dynamics with a reasonable degree of complexity. The climatology of the IGCM run at T-21 resolution with 22 levels is compared to European Centre for Medium Range Weather Forecasting Reanalysis data. The response of the model to changes in carbon dioxide and solar output are examined when these changes are applied globally and when constrained geographically (e.g. over land only). The CO₂ experiments have a roughly 17% higher climate sensitivity than the solar experiments. It is also found that a forcing at high latitudes causes a 40% higher climate sensitivity than a forcing only applied at low latitudes. It is found that, despite differences in the model feedbacks, climate sensitivity is roughly constant over a range of distributions of CO₂ and solar forcings. Hence, in the IGCM at least, the radiative forcing concept is capable of predicting global surface temperature changes to within 30%, for the perturbations described here. It is concluded that radiative forcing remains a useful tool for assessing the natural and anthropogenic impact of climate change mechanisms on surface temperature.