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.

Sea surface temperature climate change initiative: alternative image classification algorithms for sea-ice affected oceans

We present a Bayesian image classification scheme for discriminating cloud, clear and sea-ice observations at high latitudes to improve identification of areas of clear-sky over ice-free ocean for SST retrieval. We validate the image classification against a manually classified dataset using Advanced Along Track Scanning Radiometer (AATSR) data. A three way classification scheme using a near-infrared textural feature improves classifier accuracy by 9.9 % over the nadir only version of the cloud clearing used in the ATSR Reprocessing for Climate (ARC) project in high latitude regions. The three way classification gives similar numbers of cloud and ice scenes misclassified as clear but significantly more clear-sky cases are correctly identified (89.9 % compared with 65 % for ARC). We also demonstrate the poetential of a Bayesian image classifier including information from the 0.6 micron channel to be used in sea-ice extent and ice surface temperature retrieval with 77.7 % of ice scenes correctly identified and an overall classifier accuracy of 96 %.

Ice-stream initiation, duration and thinning on James Ross Island, northern Antarctic Peninsula

Predicting the future response of the Antarctic Ice Sheet to climate change requires an understanding of the ice streams that dominate its dynamics. Here we use cosmogenic isotope exposure-age dating (26Al, 10Be and 36Cl) of erratic boulders on ice-free land on James Ross Island, north-eastern Antarctic Peninsula, to define the evolution of Last Glacial Maximum (LGM) ice in the adjacent Prince Gustav Channel. These data include ice-sheet extent, thickness and dynamical behaviour. Prior to ∼18 ka, the LGM Antarctic Peninsula Ice Sheet extended to the continental shelf-edge and transported erratic boulders onto high-elevation mesas on James Ross Island. After ∼18 ka there was a period of rapid ice-sheet surface-lowering, coincident with the initiation of the Prince Gustav Ice Stream. This timing coincided with rapid increases in atmospheric temperature and eustatic sea-level rise around the Antarctic Peninsula. Collectively, these data provide evidence for a transition from a thick, cold-based LGM Antarctic Peninsula Ice Sheet to a thinner, partially warm-based ice sheet during deglaciation.

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.

Landscape evolution and ice-sheet behaviour in a semi-arid polar environment: James Ross Island, NE Antarctic Peninsula

This study of landscape evolution presents both new modern and palaeo process-landform data, and analyses the behaviour of the Antarctic Peninsula Ice Sheet through the Last Glacial Maximum (LGM), the Holocene and to the present day. Six sediment-landform assemblages are described and interpreted for Ulu Peninsula, James Ross Island, NE Antarctic Peninsula: (1) the Glacier Ice and Snow Assemblage; (2) the Glacigenic Assemblage, which relates to LGM sediments and comprises both erratic-poor and erratic-rich drift, deposited by cold-based and wet-based ice and ice streams respectively; (3) the Boulder Train Assemblage, deposited during a Mid-Holocene glacier readvance; (4) the Ice-cored Moraine Assemblage, found in front of small cirque glaciers; (5) the Paraglacial Assemblage including scree, pebble-boulder lags, and littoral and fluvial processes; and (6) the Periglacial Assemblage including rock glaciers, protalus ramparts, blockfields, solifluction lobes and extensive patterned ground. The interplay between glacial, paraglacial and periglacial processes in this semi-arid polar environment is important in understanding polygenetic landforms. Crucially, cold-based ice was capable of sediment and landform genesis and modification. This landsystem model can aid the interpretation of past environments, but also provides new data to aid the reconstruction of the last ice sheet to overrun James Ross Island.

Dynamic Devensian ice flow in NE England: a sedimentological reconstruction

This study reconstructs the depositional environments that accompanied both ice advance and ice retreat of the last British–Irish Ice Sheet in NE England during the Last Glacial Maximum, and proposes three regional ice-flow phases. The Late Devensian (29–22 cal. ka BP) Tyne Gap Ice Stream initially deposited the Blackhall Till Formation during shelf-edge glaciation (Phase I). This subglacial traction till comprises several related facies, including stratified and laminated diamictons, tectonites, and sand and gravel beds deposited both in subglacial canals and in proglacial streams. Eventually, stagnation of the Tyne Gap Ice Stream led to ice-marginal sedimentation in County Durham (Phase II). During the Dimlington Stadial (21 cal. ka BP), the North Sea Lobe advanced towards the coastline of N Norfolk. This resulted initially in sandur deposition (widespread, tabular sand and gravel; the Peterlee Sand and Gravel Formation; Phase II) and ultimately in deposition of the Horden Till Formation (Phase III), a massive subglacial till. As the North Sea Lobe overrode previous formations, it thrusted and stacked sediments in County Durham, and dammed proglacial lakes between the east-coast ice, the Pennine uplands and the remaining Pennine ice. The North Sea Lobe retreated after Heinrich Event 1 (16 ka). This study highlights the complexity of ice flow during the Late Devensian glaciation of NE England, with changing environmental and oceanic conditions forcing a mobile and sensitive ice sheet.

Antarctic Peninsula ice sheet evolution during the Cenozoic Era

The Antarctic Peninsula region is currently undergoing rapid environmental change, resulting in the thinning, acceleration and recession of glaciers and the sequential collapse of ice shelves. It is important to view these changes in the context of long-term palaeoenvironmental complexity and to understand the key processes controlling ice sheet growth and recession. In addition, numerical ice sheet models require detailed geological data for tuning and testing. Therefore, this paper systematically and holistically reviews published geological evidence for Antarctic Peninsula Ice Sheet variability for each key locality throughout the Cenozoic, and brings together the prevailing consensus of the extent, character and behaviour of the glaciations of the Antarctic Peninsula region. Major contributions include a downloadable database of 186 terrestrial and marine calibrated dates; an original reconstruction of the LGM ice sheet; and a new series of isochrones detailing ice sheet retreat following the LGM. Glaciation of Antarctica was initiated around the Eocene/Oligocene transition in East Antarctica. Palaeogene records of Antarctic Peninsula glaciation are primarily restricted to King George Island, where glacigenic sediments provide a record of early East Antarctic glaciations, but with modification of far-travelled erratics by local South Shetland Island ice caps. Evidence for Neogene glaciation is derived primarily from King George Island and James Ross Island, where glaciovolcanic strata indicate that ice thicknesses reached 500–850 m during glacials. This suggests that the Antarctic Peninsula Ice Sheet draped, rather than drowned, the topography. Marine geophysical investigations indicate multiple ice sheet advances during this time. Seismic profiling of continental shelf-slope deposits indicates up to ten large advances of the Antarctic Peninsula Ice Sheet during the Early Pleistocene, when the ice sheet was dominated by 40 kyr cycles. Glacials became more pronounced, reaching the continental shelf edge, and of longer duration during the Middle Pleistocene. During the Late Pleistocene, repeated glacials reached the shelf edge, but ice shelves inhibited iceberg rafting. The Last Glacial Maximum (LGM) occurred at 18 ka BP, after which transitional glaciomarine sediments on the continental shelf indicate ice-sheet retreat. The continental shelf contains large bathymetric troughs, which were repeatedly occupied by large ice streams during Pleistocene glaciations. Retreat after the LGM was episodic in the Weddell Sea, with multiple readvances and changes in ice-flow direction, but rapid in the Bellingshausen Sea. The late Holocene Epoch was characterised by repeated fluctuations in palaeoenvironmental conditions, with associated glacial readvances. However, this has been subsumed by rapid warming and ice-shelf collapse during the twentieth century.

Accelerating shrinkage of Patagonian glaciers from the Little Ice Age (∼AD 1870) to 2011

We used Little Ice Age (LIA) trimlines and moraines to assess changes in South American glaciers over the last ∼140 years. We determined the extent and length of 640 glaciers during the LIA (∼ AD 1870) and 626 glaciers (the remainder having entirely disappeared) in 1986, 2001 and 2011. The calculated reduction in glacierized area between the LIA and 2011 is 4131 km2 (15.4%), with 660 km2 (14.2%) being lost from the Northern Patagonia Icefield (NPI), 1643 km2 (11.4%) from the Southern Patagonia Icefield (SPI) and 306 km2 (14.4%) from Cordillera Darwin. Latitude, size and terminal environment (calving or land-terminating) exert the greatest control on rates of shrinkage. Small, northerly, land-terminating glaciers shrank fastest. Annual rates of area loss increased dramatically after 2001 for mountain glaciers north of 52° S and the large icefields, with the NPI and SPI now shrinking at 9.4 km2 a–1 (0.23% a–1) and 20.5 km2 a–1 (0.15% a–1) respectively. The shrinkage of glaciers between 52° S and 54° S accelerated after 1986, and rates of shrinkage from 1986 to 2011 remained steady. Icefield outlet glaciers, isolated glaciers and ice caps south of 54° S shrank faster from 1986 to 2001 than they did from 2001 to 2011.