Sensitivity of sea ice to wind-stress and radiative forcing since 1500: a model study of the Little Ice Age and beyond

Three different reconstructed wind-stress fields which take into account variations of the North Atlantic Oscillation, one general circulation model wind-stress field, and three radiative forcings (volcanic activity, insolation changes and greenhouse gas changes) are used with the UVic Earth System Climate Model to simulate the surface air temperature, the sea-ice cover, and the Atlantic meridional overturning circulation (AMOC) since 1500, a period which includes the Little Ice Age (LIA). The simulated Northern Hemisphere surface air temperature, used for model validation, agrees well with several temperature reconstructions. The simulated sea-ice cover in each hemisphere responds quite differently to the forcings. In the Northern Hemisphere, the simulated sea-ice area and volume during the LIA are larger than the present-day area and volume. The wind-driven changes in sea-ice area are about twice as large as those due to thermodynamic (i.e., radiative) forcing. For the sea-ice volume, changes due to wind forcing and thermodynamics are of similar magnitude. Before 1850, the simulations suggest that volcanic activity was mainly responsible for the thermodynamically produced area and volume changes, while after 1900 the slow greenhouse gas increase was the main driver of the sea-ice changes. Changes in insolation have a small effect on the sea ice throughout the integration period. The export of the thicker sea ice during the LIA has no significant effect on the maximum strength of the AMOC. A more important process in altering the maximum strength of the AMOC and the sea-ice thickness is the wind-driven northward ocean heat transport. In the Southern Hemisphere, there are no visible long-term trends in the simulated sea-ice area or volume since 1500. The wind-driven changes are roughly four times larger than those due to radiative forcing. Prior to 1800, all the radiative forcings could have contributed to the thermodynamically driven changes in area and volume. In the 1800s the volcanic forcing was dominant, and during the first part of the 1900s both the insolation changes and the greenhouse gas forcing are responsible for thermodynamically produced changes. Finally, in the latter part of the 1900s the greenhouse gas forcing is the dominant factor in determining the sea-ice changes in the Southern Hemisphere.

Evidence bearing on the mechanism of rapid deglaciation

Six northeast Atlantic cores contain planktonic foraminiferal records implying a very abrupt glacial/interglacial surface-ocean warming roughly coincident with the last deglaciation (isotopic termination II) at 127 000 yr B.P. These faunal composition curves have, however, been substantially altered by sediment mixing processes on the sea floor; they are translated downward in the core record and made to look steeper than they actually were. The reason for this abnormally large mixing impact is an interval of sediment with very low to negligible concentrations of all microfossils (surface ocean and bottom living). These low concentrations reflect a several-thousand-year interval of low productivity and little or no life in the overlying surface waters. We interpret this thorough suppression of productivity as a consequence of meltwater and icebergs flooding into the subpolar Atlantic gyre from the surrounding Northern Hemisphere ice sheets during deglaciation. The meltwater influx inhibited warm-season productivity by maintaining a well-stratified low-salinity surface layer; in winter, the low salinity layer froze, stopping nutrientrich deep waters from surfacing in normal cold-season convection. The earth's orbital configuration during this deglaciation created an unusually strong summer insolation maximum and winter insolation minimum in the Northern Hemisphere. Rapid melting and disintegration of the Northern Hemisphere ice sheets induced by strong summer insolation apparently created the meltwater influx; combined with very low winter insolation, the presence of this low-salinity meltwater layer led to unusually extensive sea-ice formation. The existence of a large region of winter sea ice across the subpolar North Atlantic during deglaciation implies a reduced supply of moisture in winter to the wasting Northern Hemisphere ice sheets. This includes the loss of winter moisture both locally from ice-covered northern waters and regionally from low-latitude winter storms no longer penetrating northward. The winter sea-ice cover thus acts as an amplifier providing positive feedback to the insolation-driven deglaciation 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.     

Simulating the Holocene climate evolution at northern high latitudes using a coupled atmosphere-sea ice-ocean-vegetation model

The response of the climate at high northern latitudes to slowly changing external forcings was studied in a 9,000-year long simulation with the coupled atmosphere-sea ice-ocean-vegetation model ECBilt-CLIO-VECODE. Only long-term changes in insolation and atmospheric CO₂ and CH₄ content were prescribed. The experiment reveals an early optimum (9–8 kyr BP) in most regions, followed by a 1–3°C decrease in mean annual temperatures, a reduction in summer precipitation and an expansion of sea-ice cover. These results are in general agreement with proxy data. Over the continents, the timing of the largest temperature response in summer coincides with the maximum insolation difference, while over the oceans, the maximum response is delayed by a few months due to the thermal inertia of the oceans, placing the strongest cooling in the winter half year. Sea ice is involved in two positive feedbacks (ice-albedo and sea-ice insulation) that lead regionally to an amplification of the thermal response in our model (7°C cooling in Canadian Arctic). In some areas, the tundra-taiga feedback results in intensified cooling during summer, most notably in northern North America. The simulated sea-ice expansion leads in the Nordic Seas to less deep convection and local weakening of the overturning circulation, producing a maximum winter temperature reduction of 7°C. The enhanced interaction between sea ice and deep convection is accompanied by increasing interannual variability, including two marked decadal-scale cooling events. Deep convection intensifies in the Labrador Sea, keeping the overall strength of the thermohaline circulation stable throughout the experiment.     

Investigating the variability of the Arctic sea ice thickness in response to a stochastic thermodynamic atmospheric forcing

The natural low frequency variability of the sea-ice thickness in the Arctic is investigated based on a 10 000 years simulation with a one-dimensional thermodynamic sea-ice model forced by random perturbations of the air surface temperature and solar radiation. The simulation results suggest that atmospheric random perturbations are integrated by the sea-ice. Moreover those perturbations occurring at the onset of ice melting force the largest ice thickness anomalies, which are successively amplified in summer by the albedo feedback and damped in winter by the feedback of the heat conduction through the ice. They also result in a global shift of the melting season which, in the mean annual cycle, leads to earlier melting as compared to the mean climatological cycle. The power spectrum of the ice anomalies suggests that the thickness of the perennial ice should vary preferentially on a time scale of approximately 20 years. The shape of the spectrum is consistent with that of a first order Markov process in which the characteristic time scale of the ice fluctuations would be the relaxation time scale associated with the linear feedback. The equivalent Markov model is constructed by linearizing the ice growth rate anomaly equations and allows us to derive an analytical expression of the feedback and of the forcing of the anomalies. The characteristic time scale depends explicitly on those model parameters involved in the atmosphere-ice interaction but also on the mean seasonal characteristics of the forcing and of the ice thickness. Received: 18 August 1999 / Accepted: 10 May 2000     

Sensitivity experiments performed with an energy balance atmosphere model coupled to an advection-diffusion ocean model

Summary In this paper a simple climate model is presented which is used to perform some sensitivity experiments. The atmospheric part is represented by a vertically and zonally averaged layer in which the surface air temperature, radiative fluxes at the surface and at the top of the atmosphere, the turbulent fluxes between atmosphere and surface and the snow cover are calculated. This atmospheric layer is coupled to a two-dimensional advection-diffusion ocean model in which the zonal overturning pattern is prescribed. The ocean model evaluates the temperature distribution, the amount of sea-ice and the meridional and vertical heat fluxes. The present-day climate simulated by the model compares reasonably well with observations of the seasonal and latitudinal distribution of temperature, radiation, surface alebdo, sea-ice and snow cover and meridional energy fluxes. Then, the sensitivity of the model-simulated present-day climate to perturbations in the incident solar radiation at the top of the atmosphere is investigated. The temperature response displays large latitudinal and seasonal variations, which is in qualitative agreement with results obtained with other climate models. It is found that the seasonal variation of sea-ice cover (and hence, the effective oceanic heat capacity) is one of the most important elements determining seasonal variations in climate sensitivity. Differences in sensitivity between the seasonal and annual mean version of the model are discussed. Finally, the equilibrium response to perturbations in some selected model variables is presented; these variables include meridional diffusion coefficients, drag coefficient, sea-ice thickness, atmospheric CO₂-concentration and cloud optical thickness.With 13 Figures     

Quantifying the AMOC feedbacks during a 2×CO2 stabilization experiment with land-ice melting

The response of the Atlantic Meridional Overturning Circulation (AMOC) to an increase in atmospheric CO₂ concentration is analyzed using the IPSL-CM4 coupled ocean–atmosphere model. Two simulations are integrated for 70 years with 1%/year increase in CO₂ concentration until 2×CO₂, and are then stabilized for further 430 years. The first simulation takes land-ice melting into account, via a simple parameterization, which results in a strong freshwater input of about 0.13 Sv at high latitudes in a warmer climate. During this scenario, the AMOC shuts down. A second simulation does not include this land-ice melting and herein, the AMOC recovers after 200 years. This behavior shows that this model is close to an AMOC shutdown threshold under global warming conditions, due to continuous input of land-ice melting. The analysis of the origin of density changes in the Northern Hemisphere convection sites allows an identification as to the origin of the changes in the AMOC. The processes that decrease the AMOC are the reduction of surface cooling due to the reduction in the air–sea temperature gradient as the atmosphere warms and the local freshening of convection sites that results from the increase in local freshwater forcing. Two processes also control the recovery of the AMOC: the northward advection of positive salinity anomalies from the tropics and the decrease in sea-ice transport through the Fram Strait toward the convection sites. The quantification of the AMOC related feedbacks shows that the salinity related processes contribute to a strong positive feedback, while feedback related to temperature processes is negative but remains small as there is a compensation between heat transport and surface heat flux in ocean–atmosphere coupled model. We conclude that in our model, AMOC feedbacks amplify land-ice melting perturbation by 2.5.     

Radiation Exchange Between Stratus Clouds and Polar Marine Surfaces

The radiative energy exchange between arctic sea-ice and stratiform clouds is studied by means of aircraft measurements and a two-stream radiation transfer model. The data have been obtained by flights of two identically instrumented aircraft during the Radiation and Eddy Flux Experiments REFLEX I in autumn 1991 and REFLEX II in winter 1993 over the arctic marginal ice zone of Fram Strait. The instrumental equipment comprised Eppley pyranometers and pyrgeometers, which measure the solar and terrestrial upwelling and downwelling hemispheric radiation flux densities, and a line-scan-camera on one aircraft to monitor the surface structure of the sea-ice. An empirical parametrization of the albedo of partly ice-covered ocean surfaces is obtained from the data, which describes the albedo increasing linearly with the concentration of the snow-covered sea-ice and with the cosine of the sun zenith angle at sun elevations below 10°. Cloud optical parameters, such as single scattering albedo, asymmetry factor and shortwave and longwave height-dependent extinction coefficient are determined by adjusting modeled radiation flux densities to observations. We found significant influence of the multiple reflection of shortwave radiation between the ice surface and the cloud base on the radiation regime. Consistent with the data, a radiation transfer model shows that stratus clouds of 400 m thickness with common cloud parameters may double the global radiation at the surface of sea-ice compared to open water values. The total cloud-surface-albedo under these circumstances is 30% larger over sea-ice than over water. Parametrizations of the global and reflected radiation above and below stratus clouds are proposed on the basis of the measurements and modeling. The upwelling and downwelling longwave emission of stratus clouds with thicknesses of more than 500 m can be satisfactorily estimated by Stefan's law with an emissivity of nearly 1 and when the maximum air temperature within the cloud is used.     

Actual and insolation-weighted Northern Hemisphere snow cover and sea-ice between 1973–2002

Actual and insolation-weighted Northern Hemisphere snow cover and sea ice are binned by latitude bands for the years 1973–2002. Antarctic sea-ice is also analyzed for the years 1980–2002. The use of insolation weighting provides an improved estimate of the radiative feedbacks of snow cover and sea-ice into the atmosphere. One conclusion of our assessment is that while a decrease in both areal and insolation-weighted values have occurred, the data does not show a monotonic decrease of either Arctic sea-ice or Northern Hemisphere snow cover. If Arctic perennial sea-ice is decreasing since the total reduction in areal coverage is relatively small, a large portion of it is being replenished each year such that its radiative feedback to the atmosphere is muted. Antarctic sea-ice areal cover shows no significant long-term trend, while there is a slight decrease in the insolation-weighted values for the period 1980–2002. From the early 1990s to 2001, there was a slight increase in both values. The comparison of general circulation model simulations of changes over the last several decades to observed changes in insolation-weighted sea-ice and snow cover should be a priority research topic.     

Sensitivity of Arctic warming to sea surface temperature distribution over melted sea-ice region in atmospheric general circulation model experiments

Substantial reduction in Arctic sea ice in recent decades has intensified air-sea interaction over the Arctic Ocean and has altered atmospheric states in the Arctic and surrounding high-latitude regions. This study has found that the atmospheric responses related to Arctic sea-ice melt in the cold season (October–March) depend on sea-ice fraction and are very sensitive to in situ sea surface temperature (SST) from a series of atmospheric general circulation model (AGCM) simulations in which multiple combinations of SSTs and sea-ice concentrations are prescribed in the Arctic Ocean. It has been found that the amplitude of surface warming over the melted sea-ice region is controlled by concurrent in situ SST even if these simulations are forced by the same sea-ice concentration. Much of the sensitivity of surface warming to in situ SST are related with large changes in surface heat fluxes such as the outgoing long-wave flux in early winter (October–December) and the sensible and latent heat fluxes for the entire cold season. Vertical extension of surface warming and moistening is sensitive to these changes as well; the associated condensational heating modulates a static stability in the lower troposphere. This study also indicates that changes in SST fields in AGCM simulations must be implemented with extra care, especially in the melted sea-ice region in the Arctic. The statistical method introduced in this study for adjusting SSTs in conjunction with a given sea-ice change can help to model the atmospheric response to sea-ice loss more accurately.     

Simulation of climate phase lags in response to precession and obliquity forcing and the role of vegetation

Long (130,000 years) transient simulations with a coupled model of intermediate complexity (CLIMBER-2) have been performed. The main objective of this study is to examine leads and lags in the response to the climate system to separate obliquity and precession-induced insolation changes. Focus is on the role of internal feedbacks in the coupled atmosphere/ocean/sea-ice/vegetation system. No interactive ice sheets were used. The results show that leads and lags occur in response to the African/Asian monsoon, temperatures at high latitudes and the Atlantic thermohaline circulation. For the monsoon, leads and lags of the monthly precipitation with respect to the precession parameter were found, which are strongly modified by vegetation. In contrast, no lag was observed for the annual precipitation. At high latitudes during late winter/early spring a vegetation-induced lag with respect to the precession parameter was found in surface air temperatures. Again, no annual lag was detected. The lag in the monthly surface air temperatures induces a lag in the annual overturning in the Atlantic Ocean by changing the strength of the deep convection. The lag is several thousand years. The obliquity-related forcing does not give rise to lags in the climate system. We conclude that lags in monthly climatic variables, which are due to vegetation feedbacks, can result in an annual lag when a climatic process (like deep water formation) acts as a filter for certain months.     

Impact of Earth’s orbit and freshwater fluxes on Holocene climate mean seasonal cycle and ENSO characteristics

We use a state-of-the-art 3-dimensional coupled model to investigate the relative impact of long term variations in the Holocene insolation forcing and of a freshwater release in the North Atlantic. We show that insolation has a greater effect on seasonality and La Ni?a events and is the major driver of sea surface temperature changes. In contrast, the variations in precipitation reflect changes in El Ni?o events. The impact of ice-sheet melting may have offset the impact of insolation on El Ni?o Southern Oscillation variability at the beginning of the Holocene. These simulations provide a coherent framework to refine the interpretation of proxy data and show that changes in seasonality may bias the projection of relationships established between proxy indicators and climate variations in the east Pacific from present day records.     

Adjustment and feedbacks in a global coupled ocean-atmosphere model

 We report the analysis of two 20-year simulations performed with the low resolution version of the IPSL coupled ocean-atmosphere model, with no flux correction at the air-sea interface. The simulated climate is characterized by a global sea surface temperature warming of about 4 °C in 20 years, driven by a net heat gain at the top of the atmosphere. Despite this drift, the circulation is quite realistic both in the ocean and the atmosphere. Several distinct periods are analyzed. The first corresponds to an adjustment during which the heat gain weakens both at the top of the atmosphere and at the ocean surface, and the tropical circulation is slightly modified. Then, the surface warming is enhanced by an increase of the greenhouse feedback. We show that the mechanisms involved in the model share common features with sensitivity experiments to greenhouse gases or to SST warming. At the top of the atmosphere, most of the longwave trapping in the atmosphere is driven by the tropical circulation. At the surface, the reduction of longwave cooling is a direct response to increased temperature and moisture content at low levels in the atmospheric model. During the last part of the simulation, a regulation occurs from evaporation at the surface and longwave cooling at TOA. Most of the model drift is attributed to a too large heating by solar radiation in middle and high latitudes. The reduction of the north–south temperature gradient, and the related changes in the meridional equator-to-pole ocean heat transport lead to a warming of equatorial and subtropical regions. This is also well demonstrated by the difference between the two simulations which differ only in the parametrization of sea-ice. When the sea-ice cover is not restored to climatology the model does not maintain sea-ice at high latitudes. The climate warms more rapidly and the water vapor and clouds feedback occurs earlier. Received: 24 May 1996 / Accepted: 29 November 1996     

Responses of Mean and Extreme Precipitation to Different Climate Forcing Under Radiative-Convective Equilibrium

Understanding the responses of mean and extreme precipitation to climate change is of great importance.Previous studies have mainly focused on the responses to prescribed sea surface warming or warming due to increases of CO2.This study uses a cloud-resolving model under the idealization of radiative-convective equilibrium to examine the responses of mean and extreme precipitation to a variety of climate forcings,including changes in prescribed sea surface temperature,CO2,solar insolation,surface albedo,stratospheric volcanic aerosols,and several tropospheric aerosols.The different responses of mean precipitation are understood by examining the changes in the surface energy budget.It is found that the cancellation between shortwave scattering and longwave radiation leads to a small dependence of the mean precipitation response on forcings.The responses of extreme precipitation are decomposed into three components(thermodynamic,dynamic,and precipitation efficiency).The thermodynamic components for all climate forcings are similar.The dynamic components and the precipitation-efficiency components,which have large spreads among the cases,are negatively correlated,leading to a small dependence of the extreme precipitation response on the forcings.     

Warm winds from the Pacific caused extensive Arctic sea-ice melt in summer 2007

During summer 2007 the Arctic sea-ice shrank to the lowest extent ever observed. The role of the atmospheric energy transport in this extreme melt event is explored using the state-of-the-art ERA-Interim reanalysis data. We find that in summer 2007 there was an anomalous atmospheric flow of warm and humid air into the region that suffered severe melt. This anomaly was larger than during any other year in the data (1989?C2008). Convergence of the atmospheric energy transport over this area led to positive anomalies of the downward longwave radiation and turbulent fluxes. In the region that experienced unusual ice melt, the net anomaly of the surface fluxes provided enough extra energy to melt roughly one meter of ice during the melting season. When the ocean successively became ice-free, the surface-albedo decreased causing additional absorption of shortwave radiation, despite the fact that the downwelling solar radiation was smaller than average. We argue that the positive anomalies of net downward longwave radiation and turbulent fluxes played a key role in initiating the 2007 extreme ice melt, whereas the shortwave-radiation changes acted as an amplifying feedback mechanism in response to the melt.     

‘Modelling the Arctic Boundary Layer: An Evaluation of Six Arcmip Regional-Scale Models using Data from the Sheba Project’

A primary climate change signal in the central Arctic is the melting of sea ice. This is dependent on the interplay between the atmosphere and the sea ice, which is critically dependent on the exchange of momentum, heat and moisture at the surface. In assessing the realism of climate change scenarios it is vital to know the quality by which these exchanges are modelled in climate simulations. Six state-of-the-art regional-climate models are run for one year in the western Arctic, on a common domain that encompasses the Surface Heat Budget of the Arctic Ocean (SHEBA) experiment ice-drift track. Surface variables, surface fluxes and the vertical structure of the lower troposphere are evaluated using data from the SHEBA experiment. All the models are driven by the same lateral boundary conditions, sea-ice fraction and sea and sea-ice surface temperatures. Surface pressure, near-surface air temperature, specific humidity and wind speed agree well with observations, with a falling degree of accuracy in that order. Wind speeds have systematic biases in some models, by as much as a few metres per second. The surface radiation fluxes are also surprisingly accurate, given the complexity of the problem. The turbulent momentum flux is acceptable, on average, in most models, but the turbulent heat fluxes are, however, mostly unreliable. Their correlation with observed fluxes is, in principle, insignificant, and they accumulate over a year to values an order of magnitude larger than observed. Typical instantaneous errors are easily of the same order of magnitude as the observed net atmospheric heat flux. In the light of the sensitivity of the atmosphere–ice interaction to errors in these fluxes, the ice-melt in climate change scenarios must be viewed with considerable caution.     

Non-linear climate feedback analysis in an atmospheric general circulation model

 A method is described for evaluating the ‘partial derivatives’ of globally averaged top-of-atmosphere (TOA) radiation changes with respect to basic climate model physical parameters. This method is used to analyse feedbacks in the Australian Bureau of Meteorology Research Centre general circulation model. The parameters considered are surface temperature, water vapour, lapse rate and cloud cover. The climate forcing which produces the changes is a globally uniform sea surface temperature (SST) perturbation. The first and second order differentials of model parameters with respect to the forcing (i.e. SST changes) are estimated from quadratic least square fitting. Except for total cloud cover, variables are found to be strong functions of global SST. Strongly non-linear variations of lapse rate and high cloud amount and height appear to relate to the non-linear response in penetrative convection. Globally averaged TOA radiation differentials with respect to model parameters are also evaluated. With the exception of total cloud contributions, a high correlation is generally found to exist, on the global mean level, between TOA radiation and the respective parameter perturbations. The largest non-linear terms contributing to radiative changes are those due to lapse rate and high cloud. The contributions of linear and non-linear terms to the overall radiative response from a 4 K SST perturbation are assessed. Significant non-linear responses are found to be associated with lapse rate, water vapour and cloud changes. Although the exact magnitude of these responses is likely to be a function of the particular model as well as the imposed SST perturbation pattern, the present experiments flag these as processes which cannot properly be understood from linear theory in the evaluation of climate change sensitivity. Received: 16 January 1997/Accepted: 9 May 1997     

A diagnostic study of future evaporation changes projected in CMIP5 climate models

It has been pointed out that climatological-mean precipitation-evaporation difference (P–E) should increase under global warming mainly through the increasing saturation level of moisture. This study focuses on evaporation changes under global warming and their dependency on the direct warming effect, on the basis of future projections from the Coupled Model Intercomparison Project Phase 5 (CMIP5). Over most of the tropical, subtropical and midlatitude regions, the direct contribution from surface temperature increase is found to dominate the projected increase in evaporation. This contribution is nevertheless offset partially, especially over the oceans, by contributions from weakening surface winds and increasing near-surface relative humidity. Greater warming of surface air than of the sea surface also acts to reduce surface evaporation, by reducing both the exchange coefficient and humidity contrast at the surface. Though generally of secondary importance, this contribution is the dominant factor over the subpolar oceans. Over the polar oceans, the effect of sea-ice retreat dominantly contributes to the evaporation increase in winter, whereas the reduced exchange coefficient and surface humidity contrast coupled with the sea-ice retreat account for most of the response during summertime. Over the continents, changes in the surface exchange coefficient, reflecting changes in soil moisture and vegetation among other factors, are important to modulate the direct effects of the warming and the generally reduced surface air relative humidity.     

Climate Change Scenarios for the Hudson Bay Region: An Intermodel Comparison

General circulation models (GCMs) are unanimous in projecting warmer temperatures in an enhanced CO₂ atmosphere, with amplification of this warming in higher latitudes. The Hudson Bay region, which is located in the Arctic and subarctic regions of Canada, should therefore be strongly influenced by global warming. In this study, we compare the response of Hudson Bay to a transient warming scenario provided by six-coupled atmosphere-ocean models. Our analysis focuses on surface temperature, precipitation, sea-ice coverage, and permafrost distribution. The results show that warming is expected to peak in winter over the ocean, because of a northward retreat of the sea-ice cover. Also, a secondary warming peak is observed in summer over land in the Canadian and Australian-coupled GCMs, which is associated with both a reduction in soil moisture conditions and changes in permafrost distribution. In addition, a relationship is identified between the retreat of the sea-ice cover and an enhancement of precipitation over both land and oceanic surfaces. The response of the sea-ice cover and permafrost layer to global warming varies considerably among models and thus large differences are observed in the projected regional increase in temperature and precipitation. In view of the important feedbacks that a retreat of the sea-ice cover and the distribution of permafrost are likely to play in the doubled and tripled CO₂ climates of Hudson Bay, a good representation of these two parameters is necessary to provide realistic climate change scenarios. The use of higher resolution regional climate model is recommended to develop scenarios of climate change for the Hudson Bay region.     

An investigation into the mechanisms of changes in mid-latitude mean sea level pressure as greenhouse gases are increased

When greenhouse gases are increased in coupled GCM experiments there is both a direct effect and an indirect effect due to changes in the surface conditions. In this study we carry out experiments with a perpetual winter atmosphere only model in order to investigate the influence of changes to the surface conditions (sea surface temperatures, sea-ice and snow amount) on the Northern Hemisphere winter mid-latitude mean sea level pressure response. The surface conditions for the perpetual winter model experiments are prescribed from time averages of the HadCM2 control and greenhouse gas experiments. Forcing the perpetual winter model with the HadCM2 greenhouse gas surface conditions produces a negative mean sea level pressure (MSLP) response across both Northern Hemisphere ocean basins, as was found in the coupled model HadCM2 experiment. Additional PW model experiments show that the sea surface temperature forcing from the HadCM2 greenhouse gas experiment dominates the snow and soil moisture content forcings. The sea-ice forcing from the HadCM2 greenhouse gas experiment reduces MSLP at high latitudes. In the north Pacific region MSLP decreases when the global mean warming is applied to the sea surface temperature forcing field at all open sea points. In the north Atlantic region the increased tropics to mid-latitude meridional sea surface temperature gradient is required for MSLP to decrease. These experiments show that the MSLP response in the Northern Hemisphere mid-latitude storm track regions is sensitive to the non-local sea surface temperature anomaly pattern.