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.     

Time-dependent response of a zonally averaged ocean–atmosphere–sea ice model to Milankovitch forcing

An ocean–atmosphere–sea ice model is developed to explore the time-dependent response of climate to Milankovitch forcing for the time interval 5–3 Myr BP. The ocean component is a zonally averaged model of the circulation in five basins (Arctic, Atlantic, Indian, Pacific, and Southern Oceans). The atmospheric component is a one-dimensional (latitudinal) energy balance model, and the sea-ice component is a thermodynamic model. Two numerical experiments are conducted. The first experiment does not include sea ice and the Arctic Ocean; the second experiment does. Results from the two experiments are used to investigate (1) the response of annual mean surface air and ocean temperatures to Milankovitch forcing, and (2) the role of sea ice in this response. In both experiments, the response of air temperature is dominated by obliquity cycles at most latitudes. On the other hand, the response of ocean temperature varies with latitude and depth. Deep water formed between 45°N and 65°N in the Atlantic Ocean mainly responds to precession. In contrast, deep water formed south of 60°S responds to obliquity when sea ice is not included. Sea ice acts as a time-integrator of summer insolation changes such that annual mean sea-ice conditions mainly respond to obliquity. Thus, in the presence of sea ice, air temperature changes over the sea ice are amplified, and temperature changes in deep water of southern origin are suppressed since water below sea ice is kept near the freezing point.     

Review of simulations of climate variability and change with the GFDL R30 coupled climate model

A review is presented of the development and simulation characteristics of the most recent version of a global coupled model for climate variability and change studies at the Geophysical Fluid Dynamics Laboratory, as well as a review of the climate change experiments performed with the model. The atmospheric portion of the coupled model uses a spectral technique with rhomboidal 30 truncation, which corresponds to a transform grid with a resolution of approximately 3.75° longitude by 2.25° latitude. The ocean component has a resolution of approximately 1.875° longitude by 2.25° latitude. Relatively simple formulations of river routing, sea ice, and land surface processes are included. Two primary versions of the coupled model are described, differing in their initialization techniques and in the specification of sub-grid scale oceanic mixing of heat and salt. For each model a stable control integration of near millennial scale duration has been conducted, and the characteristics of both the time-mean and variability are described and compared to observations. A review is presented of a suite of climate change experiments conducted with these models using both idealized and realistic estimates of time-varying radiative forcing. Some experiments include estimates of forcing from past changes in volcanic aerosols and solar irradiance. The experiments performed are described, and some of the central findings are highlighted. In particular, the observed increase in global mean surface temperature is largely contained within the spread of simulated global mean temperatures from an ensemble of experiments using observationally-derived estimates of the changes in radiative forcing from increasing greenhouse gases and sulfate aerosols.     

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.     

Centennial-scale climate change from decadally-paced explosive volcanism: a coupled sea ice-ocean mechanism

Northern Hemisphere summer cooling through the Holocene is largely driven by the steady decrease in summer insolation tied to the precession of the equinoxes. However, centennial-scale climate departures, such as the Little Ice Age, must be caused by other forcings, most likely explosive volcanism and changes in solar irradiance. Stratospheric volcanic aerosols have the stronger forcing, but their short residence time likely precludes a lasting climate impact from a single eruption. Decadally paced explosive volcanism may produce a greater climate impact because the long response time of ocean surface waters allows for a cumulative decrease in sea-surface temperatures that exceeds that of any single eruption. Here we use a global climate model to evaluate the potential long-term climate impacts from four decadally paced large tropical eruptions. Direct forcing results in a rapid expansion of Arctic Ocean sea ice that persists throughout the eruption period. The expanded sea ice increases the flux of sea ice exported to the northern North Atlantic long enough that it reduces the convective warming of surface waters in the subpolar North Atlantic. In two of our four simulations the cooler surface waters being advected into the Arctic Ocean reduced the rate of basal sea-ice melt in the Atlantic sector of the Arctic Ocean, allowing sea ice to remain in an expanded state for?>?100 model years after volcanic aerosols were removed from the stratosphere. In these simulations the coupled sea ice-ocean mechanism maintains the strong positive feedbacks of an expanded Arctic Ocean sea ice cover, allowing the initial cooling related to the direct effect of volcanic aerosols to be perpetuated, potentially resulting in a centennial-scale or longer change of state in Arctic climate. The fact that the sea ice-ocean mechanism was not established in two of our four simulations suggests that a long-term sea ice response to volcanic forcing is sensitive to the stability of the seawater column, wind, and ocean currents in the North Atlantic during the eruptions.     

Simulated decadal oscillations of the Atlantic meridional overturning circulation in a cold climate state

The significance of the Atlantic meridional overturning circulation (MOC) for regional and hemispheric climate change requires a complete understanding using fully coupled climate models. Here we present a persistent, decadal oscillation in a coupled atmosphere–ocean general circulation model. While the present study is limited by the lack of comparisons with paleo-proxy records, the purpose is to reveal a new theoretically interesting solution found in the fully-coupled climate model. The model exhibits two multi-century-long stable states with one dominated by decadal MOC oscillations. The oscillations involve an interaction between anomalous advective transport of salt and surface density in the North Atlantic subpolar gyre. Their time scale is fundamentally determined by the advection. In addition, there is a link between the MOC oscillations and North Atlantic Oscillation (NAO)-like sea level pressure anomalies. The analysis suggests an interaction between the NAO and an anomalous subpolar gyre circulation in which sea ice near and south of the Labrador Sea plays an important role in generating a large local thermal anomaly and a meridional temperature gradient. The latter induces a positive feedback via synoptic eddy activity in the atmosphere. In addition, the oscillation only appears when the Nordic Sea is completely covered by sea ice in winter, and deep convection is active only near the Irminger Sea. Such conditions are provided by a substantially colder North Atlantic climate than today.     

The stability of the MOC as diagnosed from model projections for pre-industrial, present and future climates

The stability of the Atlantic meridional overturning circulation (MOC) is investigated for various climate scenario runs, using data from the CMIP3 archive of coupled atmosphere-ocean models. Apart from atmospheric feedbacks, the sign of the salt flux into the Atlantic basin that is carried by the MOC determines whether the MOC is in the single or multiple equilibria regime. This salt advection feedback is analyzed by diagnosing the freshwater and salt budgets for the combined Atlantic and Arctic basins. Consistent with the finding that almost all coupled climate models recover from hosing experiments, it is found that most models feature a negative salt advection feedback in their pre-industrial climate: freshwater perturbations are damped by this feedback, excluding the existence of a stable off-state for the MOC. All models feature enhanced evaporation over the Atlantic basin in future climates, but for a moderate increase in radiative forcing (B1 and 2 CO₂ scenarios), there is a decrease of the fresh water flux carried by the MOC into the Atlantic (the deficit is made up by increased fresh water transport by the gyre circulation). In this forcing regime the salt advection feedback becomes less negative: for three models from an ensemble of eight it is positive in a 2 CO₂ climate, while two models feature a positive feedback in the pre-industrial climate. For even warmer climates (A1B-equilibrium and 4 CO₂) the salt feedback becomes more negative (damping) again. It is shown that the decrease in northward fresh water transport at 34°S by the MOC (in B1-equilibrium and 2 CO₂) is due to a reduction of the inflow of intermediate waters relative to thermocline waters, associated with a robust shoaling of the MOC in future, warmer climates. In A1B and 4 CO₂ climates northward freshwater transport increases again. The MOC keeps shoaling, but both intermediate and thermocline water masses freshen.     

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.     

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 ocean’s role in setting the mean position of the Inter-Tropical Convergence Zone

Through study of observations and coupled climate simulations, it is argued that the mean position of the Inter-Tropical Convergence Zone (ITCZ) north of the equator is a consequence of a northwards heat transport across the equator by ocean circulation. Observations suggest that the hemispheric net radiative forcing of climate at the top of the atmosphere is almost perfectly symmetric about the equator, and so the total (atmosphere plus ocean) heat transport across the equator is small (order 0.2 PW northwards). Due to the Atlantic ocean’s meridional overturning circulation, however, the ocean carries significantly more heat northwards across the equator (order 0.4 PW) than does the coupled system. There are two primary consequences. First, atmospheric heat transport is southwards across the equator to compensate (0.2 PW southwards), resulting in the ITCZ being displaced north of the equator. Second, the atmosphere, and indeed the ocean, is slightly warmer (by perhaps 2 °C) in the northern hemisphere than in the southern hemisphere. This leads to the northern hemisphere emitting slightly more outgoing longwave radiation than the southern hemisphere by virtue of its relative warmth, supporting the small northward heat transport by the coupled system across the equator. To conclude, the coupled nature of the problem is illustrated through study of atmosphere–ocean–ice simulations in the idealized setting of an aquaplanet, resolving the key processes at work.     

The impact of dynamic sea-ice on the climatology and climate sensitivity of a GCM: a study of past, present, and future climates

 We assess two parametrisations of sea-ice in a coupled atmosphere–mixed layer ocean–sea-ice model. One parametrisation represents the thermodynamic properties of sea-ice formation alone (THERM), while the other also includes advection of the ice (DYN). The inclusion of some sea-ice dynamics improves the model's simulation of the present day sea-ice cover when compared to observations. Two climate change scenarios are used to investigate the effect of these different parametrisations on the model's climate sensitivity. The scenarios are the equilibrium response to a doubling of atmospheric CO₂ and the response to imposed glacial boundary conditions. DYN produces a smaller temperature response to a doubling of CO₂ than THERM. The temperature response of THERM is more similar to DYN in the glacial case than in the 2×CO₂ case which implies that the climate sensitivity of THERM and DYN varies with the nature of the forcing. The different responses can largely be explained by the different distribution of Southern Hemisphere sea-ice cover in the control simulations, with the inclusion of ice dynamics playing an important part in producing the differences. This emphasises the importance of realistically simulating the reference climatic state when attempting to simulate a climate change to a prescribed forcing. The simulated glacial sea-ice cover is consistent with the limited palaeodata in both THERM and DYN, but DYN simulates a more realistic present day sea-ice cover. We conclude that the inclusion of simple ice dynamics in our model increases our confidence in the simulation of the anomaly climate. Received: 24 May 2000 / Accepted: 25 October 2000     

Effects of tropical instability wave (TIW)-induced surface wind feedback in the tropical Pacific Ocean

Tropical instability waves (TIWs) arise from oceanic instability in the eastern tropical Pacific and Atlantic Oceans, having a clear atmospheric signature that results in coupled atmosphere–ocean interactions at TIW scales. In this study, the extent to which TIW-induced surface wind feedback influences the ocean is examined using an ocean general circulation model (OGCM). The TIW-induced wind stress (τ_TIW) part is diagnostically determined using an empirical τ_TIW model from sea surface temperature (SST) fields simulated in the OGCM. The interactively represented TIW wind tends to reduce TIW activity in the ocean and influence the mean state, with largest impacts during TIW active periods in fall and winter. In December, the interactive τ_TIW forcing induces a surface cooling (an order of ?0.1 to ?0.3 °C), an increased heat flux into the ocean, a shallower mixed layer and a weakening of the South Equatorial Current in the eastern equatorial Pacific. Additionally, the TIW wind effect yields a pronounced latitudinal asymmetry of sea level field across the equator, and a change to upper thermal structure, characterized by a surface cooling and a warming below in the thermocline, leading to a decreased temperature gradient between the mixed layer and the thermocline. Processes responsible for the τ_TIW–induced cooling effects are analyzed. Vertical mixing and meridional advection are the two terms in the SST budget that are dominantly affected by the TIW wind feedback: the cooling effect from the vertical mixing on SST is enhanced, with the maximum induced cooling in winter; the warming effect from the meridional advection is reduced in July–October, but enhanced in November–December. Additional experiments are performed to separate the relative roles the affected surface momentum and heat fluxes play in the cooling effect on SST. This ocean-only modeling work indicates that the effect of TIW-induced wind feedback is small but not negligible, and may need to be adequately taken into account in large-scale climate modeling.     

The formation and maintenance of the North Water Polynya

Abstract

This paper investigates the formation and maintenance of the North Water Polynya, Baffin Bay in winter using a multi‐category sea‐ice model coupled with the Princeton ocean model. Monthly climatological atmospheric data from the National Centers for Environmental Prediction/National Center for Atmospheric Research (NCEP/NCAR) reanalysis provides the forcing. An objectively‐analysed climatology provides the initial ocean temperature and salinity. Wind stress drives the ice in a cyclonic gyre around northern Baffin Bay. Localized regions of thin ice form where wind drives ice away from coastlines or fast ice. The regions of thin ice are characterized by enhanced ice growth, exceeding 1.2 m mo^?1. In the regions of thin ice, surface ocean heat flux is also enhanced and is between 30–60 W m^?2. Surface heat flux is, in part, attributable to convective mixing and entrainment driven by ice growth. The surface heat flux reflects advection of the warm West Greenland Current. Heat and salt balances show that horizontal advective exchange counterbalances surface fluxes of heat and salt.     

On the role of flux adjustments in an idealized coupled climate model

 The effect of employing flux adjustments on the climatic response of an idealized coupled model to an imposed radiative forcing is investigated with two coupled models, one of which employs flux adjustments. A linear reduction (to the planetary longwave flux) of 4 W/m² is applied over a 70 y period and held constant thereafter. Similar model responses are found (during the initial 70 y period) for global-scale diagnostics of hemispheric air temperature due to the nearly linear surface-air temperature response to the radiative forcing. Significant regional scale differences do exist, however. As the perturbation away from the present climate grows, basin-scale diagnostics (such as meridional overturning rates) begin to diverge between flux adjusted and non-flux adjusted models. Once the imposed radiative forcing is held constant, differences in global mean air temperature of up to 0.5 °C are found, with large regional-scale differences in air temperature and overturning rates within the North Atlantic and Southern Ocean. Two additional experiments with the flux adjusted model (beginning from points further along the control integration) suggest that the elimination of much of the coupling shock before the radiative forcing is applied leads to results slightly closer to the non-flux adjusted case, although large differences still persist. In particular a dipole structure indicating an enhanced warming within the Pacific sector of the Southern Ocean, and cooling within the Atlantic sector is not reproduced by the flux adjusted models. This disparity is intimately linked to the Southern Ocean overturning cell along with the flux adjustments employed as well as the drift arising from coupling shock. If a similar form of sensitivity exists in more realistic coupled models, our results suggest: (1) perturbation experiments should not be undertaken until after the coupled model control experiment is carried out for several hundred years (thereby minimizing the coupling shock); (2) care should be exercised in the interpretation of regional-scale results (over the ocean) in coupled models which employ flux adjustments; (3) care should also be taken in interpreting even global-scale diagnostics in flux adjusted models for large perturbations about the present climate. Received: 15 November 1996 / Accepted: 4 June 1997     

Steric sea level rise over the Mediterranean Sea: present climate and scenario simulations

A regional atmosphere–ocean coupled model has been used to estimate sea level rise in the Mediterranean basin under present and future conditions. A present climate simulation has been forced by ERA40 reanalysis covering the period 1958–2001. Moreover a simulation has been forced by the global coupled model ECHAM5-MPIOM under present climate conditions for the period 1951–2000. Two other 50-year simulations have been performed under the SRESA1B scenario for the twenty-first century and differ only in temperature and salinity profiles used to relax the ocean model in the Atlantic buffer zone. The present climate simulation has been verified in terms of temperature, salinity and sea level against observed data, showing good performances both in mean values and variability over the whole Mediterranean Sea and over different sub-basins. The future scenario simulations show that the steric sea level averaged over the entire basin rises of about 2 or 7?cm in 50?years depending on the Atlantic boundary conditions. The difference of about 1?°C and 0.5?psu in the upper layers of the Atlantic sea reflects mainly on the halosteric component that contributes negatively to the sea level rise, when fresher and colder boundary conditions are used in the Atlantic buffer zone, and positively in the other case. The impact of the boundary conditions is not uniform in the basin and is particularly strong in some easternmost regions.     

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.     

Sensitivity of Hudson Bay Sea ice and ocean climate to atmospheric temperature forcing

A regional sea-ice?Cocean model was used to investigate the response of sea ice and oceanic heat storage in the Hudson Bay system to a climate-warming scenario. Projections of air temperature (for the years 2041?C2070; effective CO₂ concentration of 707?C950?ppmv) obtained from the Canadian Regional Climate Model (CRCM 4.2.3), driven by the third-generation coupled global climate model (CGCM 3) for lateral atmospheric and land and ocean surface boundaries, were used to drive a single sensitivity experiment with the delta-change approach. The projected change in air temperature varies from 0.8°C (summer) to 10°C (winter), with a mean warming of 3.9°C. The hydrologic forcing in the warmer climate scenario was identical to the one used for the present climate simulation. Under this warmer climate scenario, the sea-ice season is reduced by 7?C9?weeks. The highest change in summer sea-surface temperature, up to 5°C, is found in southeastern Hudson Bay, along the Nunavik coast and in James Bay. In central Hudson Bay, sea-surface temperature increases by over 3°C. Analysis of the heat content stored in the water column revealed an accumulation of additional heat, exceeding 3?MJ?m^?3, trapped along the eastern shore of James and Hudson bays during winter. Despite the stratification due to meltwater and river runoff during summer, the shallow coastal regions demonstrate a higher capacity of heat storage. The maximum volume of dense water produced at the end of winter was halved under the climate-warming perturbation. The maximum volume of sea ice is reduced by 31% (592?km3) while the difference in the maximum cover is only 2.6% (32,350?km²). Overall, the depletion of sea-ice thickness in Hudson Bay follows a southeast?Cnorthwest gradient. Sea-ice thickness in Hudson Strait and Ungava Bay is 50% thinner than in present climate conditions during wintertime. The model indicates that the greatest changes in both sea-ice climate and heat content would occur in southeastern Hudson Bay, James Bay, and Hudson Strait.     

Sea‐ice interaction and the stability of the thermohaline circulation

Abstract

The role of sea‐ice in affecting the stability and long‐term variability of the oceanic thermohaline circulation (THC) is studied in this paper. The emphasis is placed on studying how sea‐ice might affect the stability and the long‐term variability of the THC through modulations of the surface heat and freshwater fluxes. A simple box model is analyzed to elucidate qualitatively the distinct physical meanings of these two processes. The analytical solution of this simple model indicates that, for the long timescales considered here, the thermal insulation stabilizes the THC while the freshwater feedback increases the effective inertia of the coupled ice‐ocean system. Sea‐ice insulation lessens the negative feedback between heat flux and the SST, and therefore, allows the SST to play a greater role in counteracting changes of the THC and high latitude salinity field. The freshwater feedback effectively links the surface heat flux to a freshwater reservoir, and thus, increases the effective inertia of the coupled ocean‐ice system. A two‐dimensional ocean model coupled with a thermodynamic sea‐ice model is used to estimate quantitatively the magnitudes of these two feedbacks. The numerical experiments involve the model's responses both to initial anomalies and to changes of forcing fields. For the free response cases (model responses to initial anomalies without changing the forcing fields), the model shows that the decay rate of an initial anomaly is greater when sea‐ice is included. For small perturbations the thermal insulation effect dominates over the freshwater feedback. The latter becomes increasingly more important for larger perturbations. In response to a change of external forcing, the presence of sea‐ice reduces the magnitude and the pace of the model's response. The numerical results are qualitatively consistent with the analytical solution of the box model.     

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

Radiative forcing and response of a GCM to ice age boundary conditions: cloud feedback and climate sensitivity

 A general circulation model is used to examine the effects of reduced atmospheric CO₂, insolation changes and an updated reconstruction of the continental ice sheets at the Last Glacial Maximum (LGM). A set of experiments is performed to estimate the radiative forcing from each of the boundary conditions. These calculations are used to estimate a total radiative forcing for the climate of the LGM. The response of the general circulation model to the forcing from each of the changed boundary conditions is then investigated. About two-thirds of the simulated glacial cooling is due to the presence of the continental ice sheets. The effect of the cloud feedback is substantially modified where there are large changes to surface albedo. Finally, the climate sensitivity is estimated based on the global mean LGM radiative forcing and temperature response, and is compared to the climate sensitivity calculated from equilibrium experiments with atmospheric CO₂ doubled from present day concentration. The calculations here using the model and palaeodata support a climate sensitivity of about 1 Wm^-2 K^-1 which is within the conventional range. Received: 8 February 1997 / Accepted: 4 June 1997