The impact of 21st century Arctic sea ice decline on the climate and mass balance of the Greenland ice sheet and Svalbard’s glaciers and ice caps

The Arctic is a region particularly susceptible to rapid climate change. General circulation models (GCMs) suggest a polar amplification of any global warming signal by a factor of about 1.5 due, in part, to sea ice feedbacks. The dramatic recent decline in multi-year sea ice cover lies outside the standard deviation of the CMIP3 ensemble GCM predictions. Sea ice acts as a barrier between cold air and warmer oceans during winter, as well as inhibiting evaporation from the ocean surface water during the summer. An ice free Arctic would likely have an altered hydrological cycle with more evaporation from the ocean surface leading to changes in precipitation distribution and amount. Using the U.K. Met Office Regional Climate Model (RCM), HadRM3, the atmospheric effects of the observed and projected reduction in Arctic sea ice are investigated. The RCM is driven by the atmospheric GCM HadAM3. Both models are forced with sea surface temperature and sea ice for the period 2061-2090 from the CMIP3 HadGEM1 experiments. Here we use an RCM at 50km resolution over the Arctic and 25km over Svalbard, which captures well the present-day pattern of precipitation and provides a detailed picture of the projected changes in the behaviour of the oceanic-atmosphere moisture fluxes and how they affect precipitation. These experiments show that the projected 21stCentury sea ice decline alone causes large impacts to the surface mass balance (SMB) on Svalbard. However Greenland’s SMB is not significantly affected by sea ice decline alone, but responds with a strongly negative shift in SMB when changes to SST are incorporated into the experiments. This is the first study to characterise the impact of changes in future sea ice to Arctic terrestrial cryosphere mass balance.

A Possible Change in Mass Balance of Greenland and Antarctic Ice Sheets in the Coming Century

A high-resolution GCM is found to simulate precipitation and surface energy balance of high latitudes with high accuracy. This opens new possibilities to investigate the future mass balance of polar glaciers and its effect on sea level. The surface mass balance of the Greenland and the Antarctic ice sheets is simulated using the ECHAM3 GCM with TI06 horizontal resolution. With this model, two 5-year integrations for the present and doubled carbon dioxide conditions based on the boundary conditions provided by the ECHAM1/T21 transient experiment have been conducted. A comparison of the two experiments over Greenland and Antarctica shows to what extent the effect of climate change on the mass balance on the two largest glaciers of the world can differ. On Greenland one sees a slight decrease in accumulation and a substantial increase in melt, while on Antarctica a large increase in accumulation without melt is projected. Translating the mass balances into terms of sea-level equivalent. the Greenland discharge causes a sea level rise of 1.1 mm yr−1, while the accumulation on Antarctica tends to lower it by 0.9 mm yr−1. The change in the combined mass balance of the two continents is almost zero. The sea level change of the next century can be affected more effectively by the thermal expansion of seawater and the mass balance of smaller glaciers outside of Greenland and Antarctica.

The Greenland Ice Sheet's surface mass balance in a seasonally sea-ice free Arctic

General circulation models predict a rapid decrease in sea ice extent with concurrent increases in near surface air temperature and precipitation in the Arctic over the 21st century. This has led to suggestions that some Arctic land ice masses may experience an increase in accumulation due to enhanced evaporation from a seasonally sea ice free Arctic Ocean. To investigate the impact of this phenomenon on Greenland ice sheet climate and surface mass balance (SMB) a regional climate model, HadRM3, was used to force an insolation-temperature melt SMB model. A set of experiments designed to investigate the role of sea ice independently from sea surface temperature (SST) forcing are described. In the warmer and wetter SI + SST simulation Greenland experiences a 23% increase in winter SMB but 65% reduced summer SMB, resulting in a net decrease in the annual value. This study shows that sea ice decline contributes to the increased winter balance, causing 25% of the increase in winter accumulation; this is largest in eastern Greenland as the result of increased evaporation in the Greenland Sea. These results indicate that the seasonal cycle of Greenland's SMB will increase dramatically as global temperatures increase, with the largest changes in temperature and precipitation occurring in winter. This demonstrates that the accurate prediction of changes in sea ice cover is important for predicting Greenland SMB and ice sheet evolution.

Stability of an ice sheet on an elastic bed

The stability of stationary flow of a two-dimensional ice sheet is studied when the ice obeys a power flow law (Glen's flow law). The mass accumulation rate at the top is assumed to depend on elevation and span and the bed supporting the ice sheet consists of an elastic layer lying on a rigid surface. The normal perturbation of the free surface of the ice sheet is a singular eigenvalue problem. The singularity of the perturbation at the front of the ice sheet is considered using matched asymptotic expansions, and the eigenvalue problem is seen to reduce to that with fixed ice front. Numerical solution of the perturbation eigenvalue problem shows that the dependence of accumulation rate on elevation permits the existence of unstable solutions when the equilibrium line is higher than the bed at the ice divide. Alternatively, when the equilibrium line is lower than the bed, there are only stable solutions. Softening of the bed, expressed through a decrease of its elastic modulus, has a stabilising effect on the ice sheet.

Glaciodynamics of the central sector of the last British–Irish ice sheet in Northern England

The central sector of the last British–Irish Ice Sheet (BIIS) was characterised by considerable complexity, both in terms of its glacial stratigraphy and geomorphological signature. This complexity is reflected by the large number and long history of papers that have attempted to decipher the glaciodynamic history of the region. Despite significant advances in our understanding, reconstructions remain hotly debated and relatively local, thereby hindering attempts to piece together BIIS dynamics. This paper seeks to address these issues by reviewing geomorphological mapping evidence of palimpsest flow signatures and providing an up-to-date stratigraphy of the region. Reconciling geomorphological and sedimentological evidence with relative and absolute dating constraints has allowed us to develop a new six-stage glacial model of ice-flow history and behaviour in the central sector of the last BIIS, with three major phases of glacial advance. This includes: I. Eastwards ice flow through prominent topographic corridors of the north Pennines; II. Cessation of the Stainmore ice flow pathway and northwards migration of the North Irish Sea Basin ice divide; III. Stagnation and retreat of the Tyne Gap Ice Stream; IV. Blackhall Wood–Gosforth Oscillation; V. Deglaciation of the Solway Lowlands; and VI. Scottish Re-advance and subsequent final retreat of ice out of the central sector of the last BIIS. The ice sheet was characterised by considerable dynamism, with flow switches, initiation (and termination) of ice streams, draw-down of ice into marine ice streams, repeated ice-marginal fluctuations and the production of large volumes of meltwater, locally impounded to form ice-dammed glacial lakes. Significantly, we tie this reconstruction to work carried out and models developed for the entire ice sheet. This therefore situates research in the central sector within contemporary understanding of how the last BIIS evolved over time.

Probabilistic parameterisation of the surface mass balance--elevation feedback in regional climate model simulations of the Greenland ice sheet

We present a new parameterisation that relates surface mass balance (SMB: the sum of surface accumulation and surface ablation) to changes in surface elevation of the Greenland ice sheet (GrIS) for the MAR (Modèle Atmosphérique Régional: Fettweis, 2007) regional climate model. The motivation is to dynamically adjust SMB as the GrIS evolves, allowing us to force ice sheet models with SMB simulated by MAR while incorporating the SMB–elevation feedback, without the substantial technical challenges of coupling ice sheet and climate models. This also allows us to assess the effect of elevation feedback uncertainty on the GrIS contribution to sea level, using multiple global climate and ice sheet models, without the need for additional, expensive MAR simulations. We estimate this relationship separately below and above the equilibrium line altitude (ELA, separating negative and positive SMB) and for regions north and south of 77� N, from a set of MAR simulations in which we alter the ice sheet surface elevation. These give four “SMB lapse rates”, gradients that relate SMB changes to elevation changes. We assess uncertainties within a Bayesian framework, estimating probability distributions for each gradient from which we present best estimates and credibility intervals (CI) that bound 95% of the probability. Below the ELA our gradient estimates are mostly positive, because SMB usually increases with elevation: 0.56 (95% CI: −0.22 to 1.33) kgm−3 a−1 for the north, and 1.91 (1.03 to 2.61) kgm−3 a−1 for the south. Above the ELA, the gradients are much smaller in magnitude: 0.09 (−0.03 to 0.23) kgm−3 a−1 in the north, and 0.07 (−0.07 to 0.59) kgm−3 a−1 in the south, because SMB can either increase or decrease in response to increased elevation. Our statistically founded approach allows us to make probabilistic assessments for the effect of elevation feedback uncertainty on sea level projections (Edwards et al., 2014).

Effect of uncertainty in surface mass balance--elevation feedback on projections of the future sea level contribution of the Greenland ice sheet

We apply a new parameterisation of the Greenland ice sheet (GrIS) feedback between surface mass balance (SMB: the sum of surface accumulation and surface ablation) and surface elevation in the MAR regional climate model (Edwards et al., 2014) to projections of future climate change using five ice sheet models (ISMs). The MAR (Modèle Atmosphérique Régional: Fettweis, 2007) climate projections are for 2000–2199, forced by the ECHAM5 and HadCM3 global climate models (GCMs) under the SRES A1B emissions scenario. The additional sea level contribution due to the SMB– elevation feedback averaged over five ISM projections for ECHAM5 and three for HadCM3 is 4.3% (best estimate; 95% credibility interval 1.8–6.9 %) at 2100, and 9.6% (best estimate; 95% credibility interval 3.6–16.0 %) at 2200. In all results the elevation feedback is significantly positive, amplifying the GrIS sea level contribution relative to the MAR projections in which the ice sheet topography is fixed: the lower bounds of our 95% credibility intervals (CIs) for sea level contributions are larger than the “no feedback” case for all ISMs and GCMs. Our method is novel in sea level projections because we propagate three types of modelling uncertainty – GCM and ISM structural uncertainties, and elevation feedback parameterisation uncertainty – along the causal chain, from SRES scenario to sea level, within a coherent experimental design and statistical framework. The relative contributions to uncertainty depend on the timescale of interest. At 2100, the GCM uncertainty is largest, but by 2200 both the ISM and parameterisation uncertainties are larger. We also perform a perturbed parameter ensemble with one ISM to estimate the shape of the projected sea level probability distribution; our results indicate that the probability density is slightly skewed towards higher sea level contributions.

Ice-sheet configuration in the CMIP5/PMIP3 Last Glacial Maximum experiments

We describe the creation of a data set describing changes related to the presence of ice sheets, including ice-sheet extent and height, ice-shelf extent, and the distribution and elevation of ice-free land at the Last Glacial Maximum (LGM), which were used in LGM experiments conducted as part of the fifth phase of the Coupled Modelling Intercomparison Project (CMIP5) and the third phase of the Palaeoclimate Modelling Intercomparison Project (PMIP3). The CMIP5/PMIP3 data sets were created from reconstructions made by three different groups, which were all obtained using a model-inversion approach but differ in the assumptions used in the modelling and in the type of data used as constraints. The ice-sheet extent in the Northern Hemisphere (NH) does not vary substantially between the three individual data sources. The difference in the topography of the NH ice sheets is also moderate, and smaller than the differences between these reconstructions (and the resultant composite reconstruction) and ice-sheet reconstructions used in previous generations of PMIP. Only two of the individual reconstructions provide information for Antarctica. The discrepancy between these two reconstructions is larger than the difference for the NH ice sheets, although still less than the difference between the composite reconstruction and previous PMIP ice-sheet reconstructions. Although largely confined to the ice-covered regions, differences between the climate response to the individual LGM reconstructions extend over the North Atlantic Ocean and Northern Hemisphere continents, partly through atmospheric stationary waves. Differences between the climate response to the CMIP5/PMIP3 composite and any individual ice-sheet reconstruction are smaller than those between the CMIP5/PMIP3 composite and the ice sheet used in the last phase of PMIP (PMIP2).

Impact of oceanic processes on the carbon cycle during the last termination

During the last termination (from ~18 000 years ago to ~9000 years ago), the climate significantly warmed and the ice sheets melted. Simultaneously, atmospheric CO2 increased from ~190 ppm to ~260 ppm. Although this CO2 rise plays an important role in the deglacial warming, the reasons for its evolution are difficult to explain. Only box models have been used to run transient simulations of this carbon cycle transition, but by forcing the model with data constrained scenarios of the evolution of temperature, sea level, sea ice, NADW formation, Southern Ocean vertical mixing and biological carbon pump. More complex models (including GCMs) have investigated some of these mechanisms but they have only been used to try and explain LGM versus present day steady-state climates. In this study we use a coupled climate-carbon model of intermediate complexity to explore the role of three oceanic processes in transient simulations: the sinking of brines, stratification-dependent diffusion and iron fertilization. Carbonate compensation is accounted for in these simulations. We show that neither iron fertilization nor the sinking of brines alone can account for the evolution of CO2, and that only the combination of the sinking of brines and interactive diffusion can simultaneously simulate the increase in deep Southern Ocean δ13C. The scenario that agrees best with the data takes into account all mechanisms and favours a rapid cessation of the sinking of brines around 18 000 years ago, when the Antarctic ice sheet extent was at its maximum. In this scenario, we make the hypothesis that sea ice formation was then shifted to the open ocean where the salty water is quickly mixed with fresher water, which prevents deep sinking of salty water and therefore breaks down the deep stratification and releases carbon from the abyss. Based on this scenario, it is possible to simulate both the amplitude and timing of the long-term CO2 increase during the last termination in agreement with ice core data. The atmospheric δ13C appears to be highly sensitive to changes in the terrestrial biosphere, underlining the need to better constrain the vegetation evolution during the termination.

Twentieth-century global-mean sea-level rise: is the whole greater than the sum of the parts?

Confidence in projections of global-mean sea level rise (GMSLR) depends on an ability to account for GMSLR during the twentieth century. There are contributions from ocean thermal expansion, mass loss from glaciers and ice sheets, groundwater extraction, and reservoir impoundment. Progress has been made toward solving the “enigma” of twentieth-century GMSLR, which is that the observed GMSLR has previously been found to exceed the sum of estimated contributions, especially for the earlier decades. The authors propose the following: thermal expansion simulated by climate models may previously have been underestimated because of their not including volcanic forcing in their control state; the rate of glacier mass loss was larger than previously estimated and was not smaller in the first half than in the second half of the century; the Greenland ice sheet could have made a positive contribution throughout the century; and groundwater depletion and reservoir impoundment, which are of opposite sign, may have been approximately equal in magnitude. It is possible to reconstruct the time series of GMSLR from the quantified contributions, apart from a constant residual term, which is small enough to be explained as a long-term contribution from the Antarctic ice sheet. The reconstructions account for the observation that the rate of GMSLR was not much larger during the last 50 years than during the twentieth century as a whole, despite the increasing anthropogenic forcing. Semiempirical methods for projecting GMSLR depend on the existence of a relationship between global climate change and the rate of GMSLR, but the implication of the authors' closure of the budget is that such a relationship is weak or absent during the twentieth century.

The Dubawnt Lake palaeo-ice stream: evidence for dynamic ics sheet behaviour on the Canadian Shield and insights regarding the controls on ice-stream location and vigour

We report evidence for a major ice stream that operated over the northwestern Canadian Shield in the Keewatin Sector of the Laurentide Ice Sheet during the last deglaciation 9000-8200 (uncalibrated) yr BP. It is reconstructed at 450 km in length, 140 km in width, and had an estimated catchment area of 190000 km. Mapping from satellite imagery reveals a suite of bedforms ('flow-set') characterized by a highly convergent onset zone, abrupt lateral margins, and where flow was presumed to have been fastest, a remarkably coherent pattern of mega-scale glacial lineations with lengths approaching 13 km and elongation ratios in excess of 40:1. Spatial variations in bedform elongation within the flow-set match the expected velocity field of a terrestrial ice stream. The flow pattern does not appear to be steered by topography and its location on the hard bedrock of the Canadian Shield is surprising. A soft sedimentary basin may have influenced ice-stream activity by lubricating the bed over the downstream crystalline bedrock, but it is unlikely that it operated over a pervasively deforming till layer. The location of the ice stream challenges the view that they only arise in deep bedrock troughs or over thick deposits of 'soft' fine-grained sediments. We speculate that fast ice flow may have been triggered when a steep ice sheet surface gradient with high driving stresses contacted a proglacial lake. An increase in velocity through calving could have propagated fast ice flow upstream (in the vicinity of the Keewatin Ice Divide) through a series of thermomechanical feedback mechanisms. It exerted a considerable impact on the Laurentide Ice Sheet, forcing the demise of one of the last major ice centres.

Abrupt climate change: ice sheets and oceans in action

Large temperature variations on land, in the air, and at the ocean surface, and highly variable flux of ice-rafted debris (IRD) delivered to the North Atlantic Ocean show that rapid climate fluctuations took place during the last glacial period. These quasi-periodic, high-amplitude climate variations followed a sequence of events recognized as a rapid warming, followed by a phase of gradual cooling, and terminating with more rapid cooling and increased flux of IRD to the north Atlantic Ocean. Each cycle lasted ˜1500 years, and was followed by an almost identical sequence. These cycles are referred to as Dansgaard/Oechger cycles (D/O cycles), and approximately every fourth cycle culminated in a more pronounced cooling with a massive discharge of IRD into the north Atlantic Ocean over an interval of ˜500 years. These massive discharges of IRD are known as Heinrich layers. “Heinrich events” are thus characterized as a rapid transfer of IRD from a “source,” the bed of the Laurentide Ice Sheet (LIS), to a “sink,” the North Atlantic.

Laurentide ice streaming on the Canadian Shield: A conflict with the soft-bedded ice stream paradigm?

Elucidating the controls on the location and vigor of ice streams is crucial to understanding the processes that lead to fast disintegration of ice flows and ice sheets. In the former North American Laurentide ice sheet, ice stream occurrence appears to have been governed by topographic troughs or areas of soft-sediment geology. This paper reports robust evidence of a major paleo-ice stream over the northwestern Canadian Shield, an area previously assumed to be incompatible with fast ice flow because of the low relief and relatively hard bedrock. A coherent pattern of subglacial bedforms (drumlins and megascalle glacial lineations) demarcates the ice stream flow set, which exhibits a convergent onset zone, a narrow main trunk with abrupt lateral margins, and a lobate terminus. Variations in bedform elongation ratio within the flow set match theoretical expectations of ice velocity. In the center of the ice stream, extremely parallel megascalle glacial lineations tens of kilometers long with elongation ratios in excess of 40:1 attest to a single episode of rapid ice flow. We conclude that while bed properties are likely to be influential in determining the occurrence and vigor of ice streams, contrary to established views, widespread soft-bed geology is not an essential requirement for those ice streams without topographic control. We speculate that the ice stream acted as a release valve on ice-sheet mass balance and was initiated by the presence of a proglacial lake that destabilized the ice-sheet margin and propagated fast ice flow through a series of thermomechanical feedbacks involving ice flow and temperature.