Validity of the ice shelf water plume concept beneath Filchner-Ronne Ice Shelf

In winter, brine rejection from sea ice formation and export in the Weddell Sea, offshore of Filchner-Ronne Ice Shelf (FRIS), leads to the formation of High Salinity Shelf Water (HSSW). This dense water mass enters the cavity beneath FRIS by sinking southward down the sloping continental shelf towards the grounding line. Melting occurs when the HSSW encounters the ice shelf, and the meltwater released cools and freshens the HSSW to form a water mass known as Ice Shelf Water (ISW). If this ISW rises, the ‘ice pump’ is initiated (Lewis and Perkin, 1986), whereby the ascending ISW becomes supercooled and deposits marine ice at shallower locations due to the pressure increase in the in-situ freezing temperature. Sandh¨ager et al. (2004) were able to infer the thickness patterns of marine ice deposits at the base of FRIS (figure 1), so the primary aim of this work is to try to understand the ocean flows that determine these patterns. The plume model we use to investigate ISW flow is described fully by Holland and Feltham (accepted) so only a relatively brief outline is presented here. The plume is simulated by combining a parameterisation of ice shelf basal interaction and a multiplesize- class frazil dynamics model with an unsteady, depth-averaged reduced-gravity plume model. In the model an active region of ISW evolves above and within an expanse of stagnant ambient fluid, which is considered to be ice-free and has fixed profiles of temperature and salinity. The two main assumptions of the model are that there is a well-mixed layer underneath the ice shelf and that the ambient fluid outside the plume is stagnant with fixed properties. The topography of the ice shelf that the plume flows beneath is set to the FRIS ice shelf draft calculated by Sandh¨ager et al. (2004) masked with the grounding line from the Antarctic Digital Database (ADD Consortium, 2002). To initiate the plumes, we assume that the intrusion of dense HSSW initially causes melting at the points on the grounding line where the glaciological tributaries feeding FRIS go afloat.

Development of frazil modelling in ice shelf water plumes

Abstract Preliminary results are presented from a modelling study directed at the spatial variation of frazil ice formation and its effects on flow underneath large ice shelves. The chosen plume and frazil models are briefly introduced, and results from two simplified cases are outlined. It is found that growth and melting dominate the frazil model in the short term. Secondary nucleation converts larger crystals into several nuclei due to crystal collisions (microattrition) and fluid shear and therefore governs the ice crystal dynamics after the initial supercooling has been quenched. Frazil formation is found to have a significant depth-dependence in an idealised study of an Ice Shelf Water plume. Finally, plans for more extensive and realistic studies are discussed.

Method of determining a geophysical-scale sea ice rheology from laboratory measurements

We present a methodology that allows a sea ice rheology, suitable for use in a General Circulation Model (GCM), to be determined from laboratory and tank experiments on sea ice when combined with a kinematic model of deformation. The laboratory experiments determine a material rheology for sea ice, and would investigate a nonlinear friction law of the form τ ∝ σ n⅔, instead of the more familiar Amonton's law, τ = μσn (τ is the shear stress, μ is the coefficient of friction and σ n is the normal stress). The modelling approach considers a representative region R containing ice floes (or floe aggregates), separated by flaws. The deformation of R is imposed and the motion of the floes determined using a kinematic model, which will be motivated from SAR observations. Deformation of the flaws is inferred from the floe motion and stress determined from the material rheology. The stress over R is then determined from the area-weighted contribution from flaws and floes

Recovery mechanisms of Arctic summer sea ice

We examine the recovery of Arctic sea ice from prescribed ice-free summer conditions in simulations of 21st century climate in an atmosphere–ocean general circulation model. We find that ice extent recovers typically within two years. The excess oceanic heat that had built up during the ice-free summer is rapidly returned to the atmosphere during the following autumn and winter, and then leaves the Arctic partly through increased longwave emission at the top of the atmosphere and partly through reduced atmospheric heat advection from lower latitudes. Oceanic heat transport does not contribute significantly to the loss of the excess heat. Our results suggest that anomalous loss of Arctic sea ice during a single summer is reversible, as the ice–albedo feedback is alleviated by large-scale recovery mechanisms. Hence, hysteretic threshold behavior (or a “tipping point”) is unlikely to occur during the decline of Arctic summer sea-ice cover in the 21st century.

Assimilation of sea-ice concentration in a global climate model — physical and statistical aspects

We investigate the initialization of Northern-hemisphere sea ice in the global climate model ECHAM5/MPI-OM by assimilating sea-ice concentration data. The analysis updates for concentration are given by Newtonian relaxation, and we discuss different ways of specifying the analysis updates for mean thickness. Because the conservation of mean ice thickness or actual ice thickness in the analysis updates leads to poor assimilation performance, we introduce a proportional dependence between concentration and mean thickness analysis updates. Assimilation with these proportional mean-thickness analysis updates significantly reduces assimilation error both in identical-twin experiments and when assimilating sea-ice observations, reducing the concentration error by a factor of four to six, and the thickness error by a factor of two. To understand the physical aspects of assimilation errors, we construct a simple prognostic model of the sea-ice thermodynamics, and analyse its response to the assimilation. We find that the strong dependence of thermodynamic ice growth on ice concentration necessitates an adjustment of mean ice thickness in the analysis update. To understand the statistical aspects of assimilation errors, we study the model background error covariance between ice concentration and ice thickness. We find that the spatial structure of covariances is best represented by the proportional mean-thickness analysis updates. Both physical and statistical evidence supports the experimental finding that proportional mean-thickness updates are superior to the other two methods considered and enable us to assimilate sea ice in a global climate model using simple Newtonian relaxation.

Assimilation of sea-ice concentration in a global climate model — physical and statistical aspects

We investigate the initialisation of Northern Hemisphere sea ice in the global climate model ECHAM5/MPI-OM by assimilating sea-ice concentration data. The analysis updates for concentration are given by Newtonian relaxation, and we discuss different ways of specifying the analysis updates for mean thickness. Because the conservation of mean ice thickness or actual ice thickness in the analysis updates leads to poor assimilation performance, we introduce a proportional dependence between concentration and mean thickness analysis updates. Assimilation with these proportional mean-thickness analysis updates leads to good assimilation performance for sea-ice concentration and thickness, both in identical-twin experiments and when assimilating sea-ice observations. The simulation of other Arctic surface fields in the coupled model is, however, not significantly improved by the assimilation. To understand the physical aspects of assimilation errors, we construct a simple prognostic model of the sea-ice thermodynamics, and analyse its response to the assimilation. We find that an adjustment of mean ice thickness in the analysis update is essential to arrive at plausible state estimates. To understand the statistical aspects of assimilation errors, we study the model background error covariance between ice concentration and ice thickness. We find that the spatial structure of covariances is best represented by the proportional mean-thickness analysis updates. Both physical and statistical evidence supports the experimental finding that assimilation with proportional mean-thickness updates outperforms the other two methods considered. The method described here is very simple to implement, and gives results that are sufficiently good to be used for initialising sea ice in a global climate model for seasonal to decadal predictions.

The transient versus the equilibrium response of sea ice to global warming

To examine the long-term stability of Arctic and Antarctic sea ice, idealized simulations are carried out with the climate model ECHAM5/MPIOM. Atmospheric CO2 concentration is increased over 2000 years from pre-industrial levels to quadrupling, is then kept constant for 5940 years, is afterwards decreased over 2000 years to pre-industrial levels, and finally kept constant for 3940 years.Despite these very slow changes, the sea-ice response significantly lags behind the CO2 concentration change. This lag, which is caused by the ocean’s thermal inertia, implies that the sea-ice equilibrium response to increasing CO2 concentration is substantially underestimated by transient simulations. The sea-ice response to CO2 concentration change is not truly hysteretic and in principle reversible.We find no lag in the evolution of Arctic sea ice relative to changes in annual-mean northern-hemisphere surface temperature. The summer sea-ice cover changes linearly with respect to both CO2 concentration and temper...

Predictability of large interannual Arctic sea-ice anomalies

In projections of twenty-first century climate, Arctic sea ice declines and at the same time exhibits strong interannual anomalies. Here, we investigate the potential to predict these strong sea-ice anomalies under a perfect-model assumption, using the Max-Planck-Institute Earth System Model in the same setup as in the Coupled Model Intercomparison Project Phase 5 (CMIP5). We study two cases of strong negative sea-ice anomalies: a 5-year-long anomaly for present-day conditions, and a 10-year-long anomaly for conditions projected for the middle of the twenty-first century. We treat these anomalies in the CMIP5 projections as the truth, and use exactly the same model configuration for predictions of this synthetic truth. We start ensemble predictions at different times during the anomalies, considering lagged-perfect and sea-ice-assimilated initial conditions. We find that the onset and amplitude of the interannual anomalies are not predictable. However, the further deepening of the anomaly can be predicted for typically 1 year lead time if predictions start after the onset but before the maximal amplitude of the anomaly. The magnitude of an extremely low summer sea-ice minimum is hard to predict: the skill of the prediction ensemble is not better than a damped-persistence forecast for lead times of more than a few months, and is not better than a climatology forecast for lead times of two or more years. Predictions of the present-day anomaly are more skillful than predictions of the mid-century anomaly. Predictions using sea-ice-assimilated initial conditions are competitive with those using lagged-perfect initial conditions for lead times of a year or less, but yield degraded skill for longer lead times. The results presented here suggest that there is limited prospect of predicting the large interannual sea-ice anomalies expected to occur throughout the twenty-first century.

Seasonal to interannual Arctic sea-ice predictability in current GCMs

We establish the first inter-model comparison of seasonal to interannual predictability of present-day Arctic climate by performing coordinated sets of idealized ensemble predictions with four state-of-the-art global climate models. For Arctic sea-ice extent and volume, there is potential predictive skill for lead times of up to three years, and potential prediction errors have similar growth rates and magnitudes across the models. Spatial patterns of potential prediction errors differ substantially between the models, but some features are robust. Sea-ice concentration errors are largest in the marginal ice zone, and in winter they are almost zero away from the ice edge. Sea-ice thickness errors are amplified along the coasts of the Arctic Ocean, an effect that is dominated by sea-ice advection. These results give an upper bound on the ability of current global climate models to predict important aspects of Arctic climate.

Changes of the equilibrium line altitude since the Little Ice Age in the Nepalese Himalayas

Changes of the equilibrium-line altitude (ELA) since the end of the Little Ice Age (LIA) in eastern Nepal have been studied using glacier inventory data. The toe-to-headwall altitude ratios (THARs) for individual glaciers were calculated for 1992, and used to estimate the ELA in 1959 and at the end of the LIA. THAR for debris-free glaciers is found to be smaller than for debris-covered glaciers. The ELAs for debris-covered glaciers are higher than those for debris-free glaciers in eastern Nepal. There is considerable variation in the reconstructed change in ELA (ΔELA) between glaciers within specific regions and between regions. This is not related to climate gradients, but results from differences in glacier aspect: southeast- and south-facing glaciers show larger ΔELAs in eastern Nepal than north- or west-facing glaciers. The data suggest that the rate of ELA rise may have accelerated in the last few decades. The limited number of climate records from Nepal, and analyses using a simple ELA–climate model, suggest that the higher rate of the ΔELA between 1959 and 1992 is a result of increased warming that occurred after the 1970s at higher altitudes in Nepal.

The impact of sea ice dynamics on the Arctic climate system

Five paired global climate model experiments, one with an ice pack that only responds thermodynamically (TI) and one including sea-ice dynamics (DI), were used to investigate the sensitivity of Arctic climates to sea-ice motion. The sequence of experiments includes situations in which the Arctic was both considerably colder (Glacial Inception, ca 115,000 years ago) and considerably warmer (3 × CO2) than today. Sea-ice motion produces cooler anomalies year-round than simulations without ice dynamics, resulting in reduced Arctic warming in warm scenarios and increased Arctic cooling in cold scenarios. These changes reflect changes in atmospheric circulation patterns: the DI simulations favor outflow of Arctic air and sea ice into the North Atlantic by promoting cyclonic circulation centered over northern Eurasia, whereas the TI simulations favor southerly inflow of much warmer air from the North Atlantic by promoting cyclonic circulation centered over Greenland. The differences between the paired simulations are sufficiently large to produce different vegetation cover over >19% of the land area north of 55°N, resulting in changes in land-surface characteristics large enough to have an additional impact on climate. Comparison of the DI and TI experiments for the mid-Holocene (6000 years ago) with paleovegetation reconstructions suggests the incorporation of sea-ice dynamics yields a more realistic simulation of high-latitude climates. The spatial pattern of sea-ice anomalies in the warmer-than-modern DI experiments strongly resembles the observed Arctic Ocean sea-ice dipole structure in recent decades, consistent with the idea that greenhouse warming is already impacting the high-northern latitudes.

Radiative forcing effect of ice-age dust

During glacial periods, dust deposition rates and inferred atmospheric concentrations were globally much higher than present. According to recent model results, the large enhancement of atmospheric dust content at the last glacial maximum (LGM) can be explained only if increases in the potential dust source areas are taken into account. Such increases are to be expected, due to effects of low precipitation and low atmospheric (CO2) on plant growth. Here the modelled three-dimensional dust fields from Mahowald et al. and modelled seasonally varying surface-albedo fields derived in a parallel manner, are used to quantify the mean radiative forcing due to modern (non-anthropogenic) and LGM dust. The effect of mineralogical provenance on the radiative properties of the dust is taken into account, as is the range of optical properties associated with uncertainties about the mixing state of the dust particles. The high-latitude (poleward of 45°) mean change in forcing (LGM minus modern) is estimated to be small (–0.9 to +0.2 W m–2), especially when compared to nearly –20 W m–2 due to reflection from the extended ice sheets. Although the net effect of dust over ice sheets is a positive forcing (warming), much of the simulated high-latitude dust was not over the ice sheets, but over unglaciated regions close to the expanded dust source region in central Asia. In the tropics the change in forcing is estimated to be overall negative, and of similarly large magnitude (–2.2 to –3.2 W m–2) to the radiative cooling effect of low atmospheric (CO2). Thus, the largest long-term climatic effect of the LGM dust is likely to have been a cooling of the tropics. Low tropical sea-surface temperatures, low atmospheric (CO2) and high atmospheric dust loading may be mutually reinforcing due to multiple positive feedbacks, including the negative radiative forcing effect of dust.