Simulation of the interannual variability of the wind-driven Arctic sea-ice cover during 1958–1998

A thermodynamic-dynamic sea-ice model based on a granular material rheology developed by Tremblay and Mysak is used to study the interannual variability of the Arctic sea-ice cover during the 41-year period 1958–98. Monthly wind stress forcing derived from the National Centers for Environmental Prediction (NCEP) Reanalysis data is used to produce the year-to-year variations in the sea-ice circulation and thickness. We focus on analyzing the variability of the sea-ice volume in the Arctic Basin and the subsequent changes in sea-ice export into the Greenland Sea via Fram Strait. The relative contributions of the Fram Strait sea-ice thickness and velocity anomalies to the sea-ice export anomalies are first investigated, and the former is shown to be particularly important during several large export events. The sea-ice export anomalies for these events are next linked to prior sea-ice volume anomalies in the Arctic Basin. The origin and evolution of the sea-ice volume anomalies are then related to the sea-ice circulation and atmospheric forcing patterns in the Arctic. Large sea-ice export anomalies are generally preceded by large volume anomalies formed along the East Siberian coast due to anomalous winds which occur when the Arctic High is centered closer than usual to this coastal area. When the center of this High relocates over the Beaufort Sea and the Icelandic Low extends far into the Arctic Basin, the ice volume anomalies are transported to the Fram Strait region via the Transpolar Drift Stream. Finally, the link between the sea-ice export and the North Atlantic Oscillation (NAO) index is briefly discussed. The overall results from this study show that the Arctic Basin and its ice volume anomalies must be considered in order to fully understand the export through Fram Strait. Received: 27 January 1999 / Accepted: 8 July 1999     

On the origin and evolution of sea-ice anomalies in the Beaufort-Chukchi Sea

 The origin and space-time evolution of Beaufort-Chukchi Sea ice anomalies are studied using data and a recently developed dynamic-thermodynamic sea-ice model. First, the relative importance of anomalies of river runoff, atmospheric temperature and wind in creating anomalous sea-ice conditions in the Beaufort-Chukchi Sea is investigated. The results indicate that wind anomalies are the dominant factor responsible for creating interannual variability in the Beaufort-Chukchi Sea ice cover. Temperature anomalies appear to play a major role for longer time scale fluctuations, whereas the effects of runoff anomalies are small. The sea-ice model is then used to track the position of a positive sea-ice anomaly as it is transported by the Beaufort Gyre toward the Transpolar Drift Stream and then exported out of the Arctic Basin into the Greenland Sea via Fram Strait. The model integration shows that sea-ice anomalies originating in the western Beaufort Sea can survive a few seasonal cycles as they propogate through the Arctic Basin and can account for a notable amount of anomalous ice export into the Greenland Sea. These anomalies, however, represent a small contribution to the fresh water budget in this area when compared with sea-ice fluctuations generated by interannually varying local winds. Received: 1 May 1997/Accepted: 22 October 1997     

A mechanism of decadal variability of the sea-ice volume in the Northern Hemisphere

A long-term simulation performed with a coarse-resolution, global, atmosphere-ocean-sea-ice model displays strong decadal variability of the sea-ice volume in the Northern Hemisphere with a significant peak at about 15-18 years. This model results from the coupling of ECBILT, a spectral T21, 3-level quasi-geostrophic atmospheric model, and CLIO, a sea-ice-ocean general circulation model. First, the mechanism underlying the variability of ice volume in the model was studied by performing correlation analyses between the simulated variables. In a second step, a series of additional sensitivity experiments was performed in order to illustrate the role of specific physical processes. This has allowed us to identify a feedback loop in the ice-ocean system, which proceeds as follows: an increase in Arctic sea-ice volume induces an increase in the salinity there. This salinity anomaly is transported to the Greenland Sea where it promotes convective activity. This warms up the surface oceanic layer and the atmosphere in winter and induces a decrease of the ice volume, completing half a cycle. The changes in ice volume are driven by a geopotential height pattern characterised by centres of action of opposite signs over Greenland and the Barents-Kara-Central Arctic area. Thermodynamic feedback between the ice and the atmosphere appear also to be very important for the persistence of the oscillation. The dynamical response of the atmosphere to sea-ice and temperature anomalies at surface plays a smaller role.     

Precursor Role of Winter Sea-Ice in the Labrador Sea for Following-Spring Precipitation over Southeastern North America and Western Europe

The role of winter sea-ice in the Labrador Sea as a precursor for precipitation anomalies over southeastern North America and Western Europe in the following spring is investigated. In general terms, as the sea ice increases, the precipitation also increases. In more detail, however, analyses indicate that both the winter sea-ice and the sea surface temperature(SST)anomalies related to increases in winter sea-ice in the Labrador Sea can persist into the following spring. These features play a forcing role in the spring atmosphere, which may be the physical mechanism behind the observational relationship between the winter sea-ice and spring precipitation anomalies. The oceanic forcings in spring include Arctic sea-ice anomalies and SST anomalies in the tropical Pacific and high-latitude North Atlantic. Multi-model Coupled Model Intercomparison Project Phase 5 and Atmospheric Model Intercomparison Project simulation results show that the atmospheric circulation response to the combination of sea-ice and SST is similar to that observed, which suggests that the oceanic forcings are indeed the physical reason for the enhanced spring precipitation. Sensitivity experiments conducted using an atmospheric general circulation model indicate that the increases in precipitation over southeastern North America are mainly attributable to the effect of the SST anomalies, while the increases over Western Europe are mainly due to the sea-ice anomalies. Although model simulations reveal that the SST anomalies play the primary role in the precipitation anomalies over southeastern North America, the observational statistical analyses indicate that the area of sea-ice in the Labrador Sea seems to be the precursor that best predicts the spring precipitation anomaly.     

Impact of a projected future Arctic Sea Ice reduction on extratropical storminess and the NAO

The impact of a reduced Arctic sea ice cover on wintertime extratropical storminess is investigated by conducting atmospheric general circulation model (AGCM) experiments. The AGCM ECHAM5 is forced by the present and a projected future seasonal cycle of Arctic sea ice. In the experiment with projected sea-ice concentrations significant reductions in storminess were found during December and January in both midlatitudes and towards the Arctic. However, a substantially larger reduction in extratropical storminess was found in March, despite a smaller change in surface energy fluxes in March than in the other winter months. The projected decrease in storminess is also related to the negative phase of the North Atlantic Oscillation (NAO). The March response is consistent with a forcing from transient and quasi-stationary eddies associated with negative NAO events. The greater sensitivity to sea-ice anomalies in late winter sets this study apart from earlier ones.     

Sea-Ice Concentration Multivariate Assimilation for the Canadian East Coast in a Coupled Sea Ice–Ocean Model

The present study focuses on the impact of ocean state (i.e., salinity and temperature) updates on the sea-ice analysis and short-term forecast in an assimilative sea ice–ocean coupled system. A relatively simple sea-ice assimilation scheme was applied to the sea ice–ocean coupled North Atlantic Nucleus for European Modelling of the Ocean (NEMO) system with a focus on the Canadian East Coast. In this assimilation scheme the ocean state was updated directly based on the correlations between the model's sea-ice concentration and the upper ocean salinity and temperature. These correlations were based on a limited time ensemble generated by applying random perturbations to the atmospheric forcing fields. High deviations in the sea-ice conditions were found along the ice edge, implying that the sea-ice edge position is sensitive to small atmospheric forcing variations. Assimilation runs with and without ocean state updates (i.e., sea-ice concentration nudging) were conducted and compared for the winter of 2002. Both continuous and intermittent assimilation schemes were examined. In a continuous sea-ice assimilation experiment, the ocean direct update is unnecessary. When the sea-ice updates are introduced intermittently the ocean state has to be altered to accommodate them, or they will be rapidly diminished by the model's dynamics. The correlations between sea-ice concentration and ocean salinity and temperature based on the first 15 days of January were used for corrections during the entire winter season when, in addition to thermodynamic processes, dynamic processes are responsible for, and even dominate, sea-ice evolution on the Labrador and Newfoundland shelves. This was an adequate choice as was demonstrated by the results of the study which showed that the experiments with ocean state adjustments generated more accurate short-term sea-ice forecasts.     

Linking the northern hemisphere sea-ice reduction trend and the quasi-decadal arctic sea-ice oscillation

The nature of the reduction trend and quasi-decadal oscillation in Northern Hemisphere sea-ice extent is investigated. The trend and oscillation that seem to be two separate phenomena have been found in data. This study examines a hypothesis that the Arctic sea-ice reduction trend in the last three decades amplified the quasi-decadal Arctic sea-ice oscillation (ASIO) due to a positive ice/ocean-albedo feedback, based on data analysis and a conceptual model proposed by Ikeda et al. The theoretical, conceptual model predicts that the quasi-decadal oscillation is amplified by the thinning sea-ice, leading to the ASIO, which is driven by the strong positive feedback between the atmosphere and ice-ocean systems. Such oscillation is predicted to be out-of-phase between the Arctic Basin and the Nordic Seas with a phase difference of 3/4, with the Nordic Seas leading the Arctic. The wavelet analysis of the sea ice data reveals that the quasi-decadal ASIO occurred actively since the 1970s following the trend starting in the 1960s (i.e., as sea-ice became thinner and thinner), as the atmosphere experienced quasi-decadal oscillations during the last century. The wavelet analysis also confirms the prediction of such out-of-phase feature between these two basins, which varied from 0.62 in 1960 to 0.25 in 1995. Furthermore, a coupled ice-ocean general circulation model (GCM) was used to simulate two scenarios, one without the greenhouse gas warming and the other having realistic atmospheric forcing along with the warming that leads to sea-ice reduction trend. The quasi-decadal ASIO is excited in the latter case compared to the no-warming case. The wavelet analyses of the simulated ice volume were also conducted to derive decadal ASIO and similar phase relationship between the Arctic Ocean and the Nordic Seas. An independent data source was used to confirm such decadal oscillation in the upper layer (or freshwater) thickness, which is consistent with the model simulation. A modified feedback loop for the sea-ice trend and ASIO was proposed based on the previous one by Mysak and Venegas and the ice/albedo and cloud/albedo feedabcks, which are responsible for the sea ice reduction trend.     

Influence of Low-frequency Solar Forcing on the East Asian Winter Monsoon Based on HadCM3 and Observations

In this study, we investigate the influence of low-frequency solar forcing on the East Asian winter monsoon(EAWM)by analyzing a four-member ensemble of 600-year simulations performed with Had CM3(Hadley Centre Coupled Model,version 3). We find that the EAWM is strengthened when total solar irradiance(TSI) increases on the multidecadal time scale. The model results indicate that positive TSI anomalies can result in the weakening of Atlantic meridional overturning circulation, causing negative sea surface temperature(SST) anomalies in the North Atlantic. Especially for the subtropical North Atlantic, the negative SST anomalies can excite an anomalous Rossby wave train that moves from the subtropical North Atlantic to the Greenland Sea and finally to Siberia. In this process, the positive sea-ice feedback over the Greenland Sea further enhances the Rossby wave. The wave train can reach the Siberian region, and strengthen the Siberian high. As a result, low-level East Asian winter circulation is strengthened and the surface air temperature in East Asia decreases. Overall,when solar forcing is stronger on the multidecadal time scale, the EAWM is typically stronger than normal. Finally, a similar linkage can be observed between the EAWM and solar forcing during the period 1850–1970.     

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.     

The role of sea ice in the temperature-precipitation feedback of glacial cycles

The response of the hydrological cycle to climate variability and change is a critical open question, where model reliability is still unsatisfactory, yet upon which past climate history can shed some light. Sea ice is a key player in the climate system and in the hydrological cycle, due to its strong albedo effect and its insulating effect on local evaporation and air-sea heat flux. Using an atmospheric general circulation model with specified sea surface temperature and sea-ice distribution, the role of sea ice in the hydrological cycle is investigated under last glacial maximum (LGM) and present day conditions, and by studying its contribution to the “temperature-precipitation feedback”. By conducting a set of sensitivity experiments in which the albedo and thickness of the sea ice are varied, the various effects of sea ice in the hydrological cycle are isolated. It is demonstrated that for a cold LGM like state, a warmer climate (as a result of reduced sea-ice cover) leads to an increase in snow precipitation over the ice sheets. The insulating effect of the sea ice on the hydrological cycle is found to be larger than the albedo effect. These two effects interact in a nonlinear way and their total effect is not equal to summing their separate contribution.     

Modelling the sea ice-ocean seasonal cycle in Hudson Bay, Foxe Basin and Hudson Strait, Canada

The seasonal cycle of water masses and sea ice in the Hudson Bay marine system is examined using a three-dimensional coastal ice-ocean model, with 10 km horizontal resolution and realistic tidal, atmospheric, hydrologic and oceanic forcing. The model includes a level 2.5 turbulent kinetic energy equation, multi-category elastic-viscous-plastic sea-ice rheology, and two layer sea ice with a single snow layer. Results from a two-year long model simulation between August 1996 and July 1998 are analyzed and compared with various observations. The results demonstrate a consistent seasonal cycle in atmosphere-ocean exchanges and the formation and circulation of water masses and sea ice. The model reproduces the summer and winter surface mixed layers, the general cyclonic circulation including the strong coastal current in eastern Hudson Bay, and the inflow of oceanic waters into Hudson Bay. The maximum sea-ice growth rates are found in western Foxe Basin, and in a relatively large and persistent polynya in northwestern Hudson Bay. Sea-ice advection and ridging are more important than local thermodynamic growth in the regions of maximum sea-ice cover concentration and thickness that are found in eastern Foxe Basin and southern Hudson Bay. The estimate of freshwater transport to the Labrador Sea confirms a broad maximum during wintertime that is associated with the previous summers freshwater moving through Hudson Strait from southern Hudson Bay. Tidally driven mixing is shown to have a strong effect on the modeled ice-ocean circulation.     

Late Pleistocene variations in Antarctic sea ice I: effect of orbital insolation changes

The effect of orbitally induced insolation changes on Antarctic sea-ice cover are examined by means of a dynamic-thermodynamic seaice model. Results are compared with modified CLIMAP 18 000 B.P. sea-ice reconstructions. Calculations suggest that changes in insolation receipt had only a minor influence on Pleistocene sea-ice distributions. The small response can be explained by a number of factors: albedo effects reduce the insolation perturbation at the surface; some of the shortwave radiation entering the ocean contributes to bottom ablation rather than lateral melting; the radiation perturbation at the upper surface of the ice must go to warming the surface to the melting point before melting ensues; and, finally, the relatively high heat capacity of open water dampens the surface temperature response to altered seasonal insolation perturbations.     

A 20-year coupled ocean-sea ice-atmosphere variability mode in the North Atlantic in an AOGCM

In order to understand potential predictability of the ocean and climate at the decadal time scales, it is crucial to improve our understanding of internal variability at this time scale. Here, we describe a 20-year mode of variability found in the North Atlantic in a 1,000-year pre-industrial simulation of the IPSL-CM5A-LR climate model. This mode involves the propagation of near-surface temperature and salinity anomalies along the southern branch of the subpolar gyre, leading to anomalous sea-ice melting in the Nordic Seas, which then forces sea-level pressure anomalies through anomalous surface atmospheric temperatures. The wind stress associated to this atmospheric structure influences the strength of the East Greenland Current across the Denmark Strait, which, in turn, induces near-surface temperature and salinity anomalies of opposite sign at the entrance of the Labrador Sea. This starts the second half of the cycle after approximatively 10 years. The time scale of the cycle is thus essentially set by advection of tracers along the southern branch of the subpolar gyre, and by the time needed for anomalous East Greenland Current to accumulate heat and freshwater anomalies at the entrance of the Labrador Sea. The Atlantic meridional overturning circulation (AMOC) does not play a dominant role in the mode that is confined in the subpolar North Atlantic, but it also has a 20-year preferred timescale. This is due to the influence of the propagating salinity anomalies on the oceanic deep convection. The existence of this preferred timescale has important implications in terms of potential predictability of the North Atlantic climate in the model, although its realism remains questionable and is discussed.     

A data-model intercomparison study of Arctic sea-ice variability

The dynamic-thermodynamic granular rheology sea-ice model of Tremblay and Mysak is validated against 40 years of observed sea-ice concentration (SIC) data. Subsequently, the mechanisms responsible for producing SIC anomalies in the model are evaluated by studying the coupled variance (using the singular value decomposition method, SVD) between the simulated SIC anomalies and the ice speed and air temperature anomalies. To execute this validation, a 49-year (1949-97) simulation (including a 9-year spin-up period) of the Arctic and peripheral sea-ice cover using daily varying winds and monthly mean air temperatures is produced. In general, the simulated SIC variations for 1958-97 in the East Siberian, Chukchi and Beaufort seas are in agreement with observations, while larger discrepancies occur in the Laptev and Kara seas. Moreover, the sensitivity of the model to southerly wind anomalies in creating summer SIC anomalies compares well with the observed sensitivity; however, the model's sensitivity to summer air temperature anomalies is weaker than observed. The summer SIC anomalies over an entire sea are not influenced by variations in the level of river runoff. Results from the SVD analysis show that the main source of variability in the peripheral seas is associated with the variation in the strength of the Arctic High; in the East Siberian and Laptev seas, the strengthening and weakening of the Transpolar Drift Stream also play an important role. Over the entire Arctic domain, surface air temperature anomalies are negatively correlated with sea-ice anomalies. Finally, the observed downward trend in total sea-ice cover in the last two decades as well as record minima in the East Siberian Sea are well reproduced in the simulation.     

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     

Solar irradiance forcing of centennial climate variability during the Holocene

Centennial climate variability during the Holocene has been simulated in two 10,000 year experiments using the intermediate-complexity ECBilt model. ECBilt contains a dynamic atmosphere, a global 3-D ocean model and a thermodynamic sea-ice model. One experiment uses orbital forcing and solar irradiance forcing, which is based on the Stuiver et al. residual ¹⁴C record spliced into the Lean et al. reconstruction. The other experiment uses orbital forcing alone. A glacier model is coupled off-line to the climate model. A time scale analysis shows that the response in atmospheric parameters to the irradiance forcing can be characterised as the direct response of a system with a large thermal inertia. This is evident in parameters like surface air temperature, monsoon precipitation and glacier length, which show a stronger response for longer time scales. The oceanic response, on the other hand, is strongly modified by internal feedback processes. The solar irradiance forcing excites a (damped) mode of the thermohaline circulation (THC) in the North Atlantic Ocean, similar to the loop-oscillator modes associated with random-noise freshwater forcing. This results in a significant peak (at time scales 200–250 year) in the THC spectrum which is absent in the reference run. The THC response diminishes the sea surface temperature response at high latitudes, while it gives rise to a signal in the sea surface salinity. A comparison of the model results with observations shows a number of encouraging similarities.     

On the role of high latitude ice,snow, and vegetation feedbacks in the climatic response to external forcing changes

A seasonal energy balance climate model containing a detailed treatment of surface and planetary albedo, and in which seasonally varying land snow and sea ice amounts are simulated in terms of a number of explicit physical processes, is used to investigate the role of high latitude ice, snow, and vegetation feedback processes. Feedback processes are quantified by computing changes in radiative forcing and feedback factors associated with individual processes. Global sea ice albedo feedback is 5–8 times stronger than global land snowcover albedo feedback for a 2% solar constant increase or decrease, with Southern Hemisphere cryosphere feedback being 2–5 times stronger than Northern Hemisphere cryosphere feedback.In the absence of changes in ice extent, changes in ice thickness in response to an increase in solar constant are associated with an increase in summer surface melting which is exactly balanced by increased basal winter freezing, and a reduction in the upward ocean-air flux in summer which is exactly balanced by an increased flux in winter, with no change in the annual mean ocean-air flux. Changes in the mean annual ocean-air heat flux require changes in mean annual ice extent, and are constrained to equal the change in meridional oceanic heat flux convergence in equilibrium. Feedback between ice extent and the meridional oceanic heat flux obtained by scaling the oceanic heat diffusion coefficient by the ice-free fraction regulates the feedback between ice extent and mean annual air-sea heat fluxes in polar regions, and has a modest effect on model-simulated high latitude temperature change.Accounting for the partial masking effect of vegetation on snow-covered land reduces the Northern Hemisphere mean temperature response to a 2% solar constant decrease or increase by 20% and 10%, respectively, even though the radiative forcing change caused by land snowcover changes is about 3 times larger in the absence of vegetational masking. Two parameterizations of the tundra fraction are tested: one based on mean annual land air temperature, and the other based on July land air temperature. The enhancement of the mean Northern Hemisphere temperature response to solar constant changes when the forest-tundra ecotone is allowed to shift with climate is only 1/3 to 1/2 that obtained by Otterman et al. (1984) when the mean annual parameterization is used here, and only 1/4 to 1/3 as large using the July parameterization.The parameterized temperature dependence of ice and snow albedo is found to enhance the global mean temperature response to a 2% solar constant increase by only 0.04 °C, in sharp contrast to the results of Washington and Meehl (1986) obtained with a mean annual model. However, there are significant differences in the method used here and in Washington and Meehl to estimate the importance of this feedback process. When their approach is used in a mean annual version of the present model, closer agreement to their results is obtained.     

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.     

Multiple sea-ice states and abrupt MOC transitions in a general circulation ocean model

Sea ice has been suggested, based on simple models, to play an important role in past glacial–interglacial oscillations via the so-called “sea-ice switch” mechanism. An important requirement for this mechanism is that multiple sea-ice extents exist under the same land ice configuration. This hypothesis of multiple sea-ice extents is tested with a state-of-the-art ocean general circulation model coupled to an atmospheric energy–moisture-balance model. The model includes a dynamic-thermodynamic sea-ice module, has a realistic ocean configuration and bathymetry, and is forced by annual mean forcing. Several runs with two different land ice distributions represent present-day and cold-climate conditions. In each case the ocean model is initiated with both ice-free and fully ice-covered states. We find that the present-day runs converge approximately to the same sea-ice state for the northern hemisphere while for the southern hemisphere a difference in sea-ice extent of about three degrees in latitude between the different runs is observed. The cold climate runs lead to meridional sea-ice extents that are different by up to four degrees in latitude in both hemispheres. While approaching the final states, the model exhibits abrupt transitions from extended sea-ice states and weak meridional overturning circulation, to less extended sea ice and stronger meridional overturning circulation, and vice versa. These transitions are linked to temperature changes in the North Atlantic high-latitude deep water. Such abrupt changes may be associated with Dansgaard–Oeschger events, as proposed by previous studies. Although multiple sea ice states have been observed, the difference between these states is not large enough to provide a strong support for the sea-ice-switch mechanism.     

Properties of adjoint sea ice sensitivities to atmospheric forcing and implications for the causes of the long term trend of Arctic sea ice

Based on adjoint sensitivities of the coupled Massachusetts Institute of Technology ocean–sea ice circulation model, the potential influence of thermodynamic atmospheric forcing on the interannual variability of the September sea ice area (AREA) and volume (VOLUME) in the Arctic is investigated for the three periods 1980–1989, 1990–1999 and 2000–2009. Sensitivities suggest that only large forcing anomalies prior to the spring melting onset in May can influence the September sea ice characteristics while even small changes in the atmospheric variables during subsequent months can significantly influence September sea ice state. Specifically, AREA close to the ice edge in the Arctic seas is highly sensitive to thermodynamic atmospheric forcing changes from June to July. In contrast, VOLUME is highly sensitive to atmospheric temperature changes occurring during the same period over the central parts of the Arctic Ocean. A comparison of the sea ice conditions and sensitivities during three different periods reveals that, due to the strong decline of sea ice concentration and sea ice thickness, sea ice area became substantially more sensitive to the same amplitude thermodynamic atmospheric forcing anomalies during 2000–2009 relative to the earlier periods. To obtain a quantitative estimate of changes that can be expected from existing atmospheric trends, adjoint sensitivities are multiplied by monthly temperature differences between 1980s and two following decades. Strongest contributions of surface atmospheric temperature differences to AREA and VOLUME changes are observed during May and September. The strongest contribution from the downward long-wave heat flux to AREA changes occurs in September and to VOLUME changes in July–August. About 62 % of the AREA decrease simulated by the model can be explained by summing all contributions to the thermodynamic atmospheric forcing. The changing sea ice state (sensitivity) is found to enhance the decline and accounts for about one third of the explained reduction. For the VOLUME decrease, the explained fraction of the decrease is only about 37 %.