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.

Using the surface profiles of modern ice masses to inform palaeo-glacier reconstructions

Morphometric study of modern ice masses is useful because many reconstructions of glaciers traditionally draw on their shape for guidance Here we analyse data derived from the surface profiles of 200 modern ice masses-valley glaciers icefields ice caps and ice sheets with length scales from 10º to 10³ km-from different parts of the world Four profile attributes are investigated relief span and two parameters C* and C that result from using Nye s (1952) theoretical parabola as a profile descriptor C* and C respectively measure each profile s aspect ratio and steepness and are found to decrease in size and variability with span This dependence quantifies the competing influences of unconstrained spreading behaviour of ice flow and bed topography on the profile shape of ice masses which becomes more parabolic as span Increases (with C* and C tending to low values of 2.5-3.3 m ½) The same data reveal coherent minimum bounds in C* and C for modern ice masses that we develop into two new methods of palaeo glacier reconstruction In the first method glacial limits are known from moraines and the bounds are used to constrain the lowest palaeo ice surface consistent with modern profiles We give an example of applying this method over a three-dimensional glacial landscape in Kamchatka In the second method we test the plausibility of existing reconstructions by comparing their C* and C against the modern minimum bounds Of the 86 published palaeo ice masses that we put to this test 88% are found to be plausible The search for other morphometric constraints will help us formalise glacier reconstructions and reduce their uncertainty and subjectiveness

The impact of a seasonally ice free Arctic Ocean on the temperature, precipitation and surface mass balance of Svalbard

The observed decline in summer sea ice extent since the 1970s is predicted to continue until the Arctic Ocean is seasonally ice free during the 21st Century. This will lead to a much perturbed Arctic climate with large changes in ocean surface energy flux. Svalbard, located on the present day sea ice edge, contains many low lying ice caps and glaciers and is expected to experience rapid warming over the 21st Century. The total sea level rise if all the land ice on Svalbard were to melt completely is 0.02 m. The purpose of this study is to quantify the impact of climate change on Svalbard’s surface mass balance (SMB) and to determine, in particular, what proportion of the projected changes in precipitation and SMB are a result of changes to the Arctic sea ice cover. To investigate this a regional climate model was forced with monthly mean climatologies of sea surface temperature (SST) and sea ice concentration for the periods 1961–1990 and 2061–2090 under two emission scenarios. In a novel forcing experiment, 20th Century SSTs and 21st Century sea ice were used to force one simulation to investigate the role of sea ice forcing. This experiment results in a 3.5 m water equivalent increase in Svalbard’s SMB compared to the present day. This is because over 50 % of the projected increase in winter precipitation over Svalbard under the A1B emissions scenario is due to an increase in lower atmosphere moisture content associated with evaporation from the ice free ocean. These results indicate that increases in precipitation due to sea ice decline may act to moderate mass loss from Svalbard’s glaciers due to future Arctic warming.

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.

Evolution of Melt Pond Volume on the Surface of the Greenland Ice Sheet

The presence of surface meltwater on ice caps and ice sheets is an important glaciological and climatological characteristic. We describe an algorithm for estimating the depth and hence volume of surface melt ponds using multispectral ASTER satellite imagery. The method relies on reasonable assumptions about the albedo of the bottom surface of the ponds and the optical attenuation characteristics of the ponded meltwater. We apply the technique to sequences of satellite imagery acquired over the western margin of the Greenland Ice Sheet to derive changes in melt pond extent and volume during the period 2001 - 2004. Results show large intra- and interannual changes in ponded water volumes, and large volumes of liquid water stored in extensive slush zones.

Ice loss in the Canadian Arctic Archipelago

Though much attention has been focused in recent years on the melting of ice from Greenland and Antarctica, nearly half of the ice volume currently being lost to the ocean is actually coming from other mountain glaciers and ice caps. Ice loss from a group of islands in northern Canada accounts for much of that volume. In a study published in April 2011 in the journal Nature, a team of researchers led by Alex Gardner of the University of Michigan found that land ice in both the northern and southern Canadian Arctic Archipelago has declined sharply. The maps above show ice loss from surface melting for the northern portion of the archipelago from 2004-2006 (left) and 2007-2009 (right). Blue indicates ice gain, and red indicates ice loss. In the six years studied, the Canadian Arctic Archipelago lost an average of approximately 61 gigatons of ice per year. (A gigaton is a billion tons of ice.) The research team also found the rate of ice loss was accelerating. From 2004 to 2006, the average mass loss was roughly 31 gigatons per year; from 2007 to 2009, the loss increased to 92 gigatons per year. Gardner and colleagues used three independent methods to assess ice mass, all of which showed the same trends. The team used a model to estimate the surface mass balance of ice and the amount of ice discharged. They also compiled and analyzed measurements from NASA's Ice, Cloud and Land Elevation Satellite (ICESat) to assess changes in the surface height of ice. Finally, they gathered observations from NASA's Gravity Recovery and Climate Experiment (GRACE) to determine changes in the gravity field in the region, an indicator of the amount of ice gained or lost. The Canadian Arctic Archipelago generally receives little precipitation, and the amount of snowfall changes little from year to year. But the rate of snow and ice melting varies considerably, so changes in ice mass come largely from changes in summertime melt. During the 2004 to 2009 study period, the Canadian Arctic Archipelago experienced four of its five warmest years since 1960, likely fueling the melting. Gardner notes that from 2001 to 2004, the sum of melting from all mountain glaciers and ice caps around the world (but not the Greenland and Antarctic ice sheets) contributed an estimated 1 millimeter per year to global sea level rise. Recent estimates suggest the Greenland and Antarctic ice sheets add another 1.3 millimeters per year to sea level. "This means 1 percent of the land ice volume-mountain glaciers and ice caps-account for about half of all ice loss to the world's oceans," Gardner said. "Most of the ice loss is coming from the Canadian Arctic Archipelago, Alaska, Patagonia, the Himalayas, and the smaller ice masses surrounding the main Greenland and Antarctic ice sheets."

Heavy minerals, grain size distribution, and ice-rafted debris in the Greenland Sea during 44 to 30 Ma

The widely accepted age estimate for the onset of glaciation in the Northern Hemisphere ranges between 2 and 15 million years ago (Ma). However, recent studies indicate the date for glacial onset may be significantly older. We report the presence of ice-rafted debris (IRD) in ~44 to 30 Ma sediments from the Greenland Sea, evidence for glaciation in the North Atlantic during the Middle Eocene to Early Oligocene. Detailed sedimentological evidence indicates that glaciers extended to sea level in the region, allowing icebergs to be produced. IRD may have been sourced from tidewater glaciers, small ice caps, and/or a continental ice sheet.

Surface velocity and mass balance of Livingston Island ice cap, Antarctica

The mass budget of the ice caps surrounding the Antarctica Peninsula and, in particular, the partitioning of its main components are poorly known. Here we approximate frontal ablation (i.e. the sum of mass losses by calving and submarine melt) and surface mass balance of the ice cap of Livingston Island, the second largest island in the South Shetland Islands archipelago, and analyse variations in surface velocity for the period 2007–2011. Velocities are obtained from feature tracking using 25 PALSAR-1 images, and used in conjunction with estimates of glacier ice thicknesses inferred from principles of glacier dynamics and ground-penetrating radar observations to estimate frontal ablation rates by a flux-gate approach. Glacier-wide surface mass-balance rates are approximated from in situ observations on two glaciers of the ice cap. Within the limitations of the large uncertainties mostly due to unknown ice thicknesses at the flux gates, we find that frontal ablation (−509 ± 263 Mt yr−1, equivalent to −0.73 ± 0.38 m w.e. yr−1 over the ice cap area of 697 km2) and surface ablation (−0.73 ± 0.10 m w.e. yr−1) contribute similar shares to total ablation (−1.46 ± 0.39 m w.e. yr−1). Total mass change (δM = −0.67 ± 0.40 m w.e. yr−1) is negative despite a slightly positive surface mass balance (0.06 ± 0.14 m w.e. yr−1). We find large interannual and, for some basins, pronounced seasonal variations in surface velocities at the flux gates, with higher velocities in summer than in winter. Associated variations in frontal ablation (of ~237 Mt yr−1; −0.34 m w.e. yr−1) highlight the importance of taking into account the seasonality in ice velocities when computing frontal ablation with a flux-gate approach.

Ice-core evidence of abrupt climate changes

Ice-core records show that climate changes in the past have been large, rapid, and synchronous over broad areas extending into low latitudes, with less variability over historical times. These ice-core records come from high mountain glaciers and the polar regions, including small ice caps and the large ice sheets of Greenland and Antarctica.

Mechanics of deformable glacier beds

My focus in this thesis is to contribute to a more thorough understanding of the mechanics of ice and deformable glacier beds. Glaciers flow under their own weight through a combination of deformation within the ice column and basal slip, which involves both sliding along and deformation within the bed. Deformable beds, which are made up of unfrozen sediment, are prevalent in nature and are often the primary contributors to ice flow wherever they are found. Their granular nature imbues them with unique mechanical properties that depend on the granular structure and hydrological properties of the bed. Despite their importance for understanding glacier flow and the response of glaciers to changing climate, the mechanics of deformable glacier beds are not well understood.

Our general approach to understanding the mechanics of bed deformation and their effect on glacier flow is to acquire synoptic observations of ice surface velocities and their changes over time and to use those observations to infer the mechanical properties of the bed. We focus on areas where changes in ice flow over time are due to known environmental forcings and where the processes of interest are largely isolated from other effects. To make this approach viable, we further develop observational methods that involve the use of mapping radar systems. Chapters 2 and 5 focus largely on the development of these methods and analysis of results from ice caps in central Iceland and an ice stream in West Antarctica. In Chapter 3, we use these observations to constrain numerical ice flow models in order to study the mechanics of the bed and the ice itself. We show that the bed in an Iceland ice cap deforms plastically and we derive an original mechanistic model of ice flow over plastically deforming beds that incorporates changes in bed strength caused by meltwater flux from the surface. Expanding on this work in Chapter 4, we develop a more detailed mechanistic model for till-covered beds that helps explain the mechanisms that cause some glaciers to surge quasi-periodically. In Antarctica, we observe and analyze the mechanisms that allow ocean tidal variations to modulate ice stream flow tens of kilometers inland. We find that the ice stream margins are significantly weakened immediately upstream of the area where ice begins to float and that this weakening likely allows changes in stress over the floating ice to propagate through the ice column.

Using Earth Deformation Caused by Surface Mass Loading to Constrain the Elastic Structure of the Crust and Mantle

Surface mass loads come in many different varieties, including the oceans, atmosphere, rivers, lakes, glaciers, ice caps, and snow fields. The loads migrate over Earth's surface on time scales that range from less than a day to many thousand years. The weights of the shifting loads exert normal forces on Earth's surface. Since the Earth is not perfectly rigid, the applied pressure deforms the shape of the solid Earth in a manner controlled by the material properties of Earth's interior. One of the most prominent types of surface mass loading, ocean tidal loading (OTL), comes from the periodic rise and fall in sea-surface height due to the gravitational influence of celestial objects, such as the moon and sun. Depending on geographic location, the surface displacements induced by OTL typically range from millimeters to several centimeters in amplitude, which may be inferred from Global Navigation and Satellite System (GNSS) measurements with sub-millimeter precision. Spatiotemporal characteristics of observed OTL-induced surface displacements may therefore be exploited to probe Earth structure. In this thesis, I present descriptions of contemporary observational and modeling techniques used to explore Earth's deformation response to OTL and other varieties of surface mass loading. With the aim to extract information about Earth's density and elastic structure from observations of the response to OTL, I investigate the sensitivity of OTL-induced surface displacements to perturbations in the material structure. As a case study, I compute and compare the observed and predicted OTL-induced surface displacements for a network of GNSS receivers across South America. The residuals in three distinct and dominant tidal bands are sub-millimeter in amplitude, indicating that modern ocean-tide and elastic-Earth models well predict the observed displacement response in that region. Nevertheless, the sub-millimeter residuals exhibit regional spatial coherency that cannot be explained entirely by random observational uncertainties and that suggests deficiencies in the forward-model assumptions. In particular, the discrepancies may reveal sensitivities to deviations from spherically symmetric, non-rotating, elastic, and isotropic (SNREI) Earth structure due to the presence of the South American craton.

A review of topographic controls on moraine distribution

Ice-marginal moraines are often used to reconstruct the dimensions of former ice masses, which are then used as proxies for palaeoclimate. This approach relies on the assumption that the distribution of moraines in the modern landscape is an accurate reflection of former ice margin positions during climatically controlled periods of ice margin stability. However, the validity of this assumption is open to question, as a number of additional, nonclimatic factors are known to influence moraine distribution. This review considers the role played by topography in this process, with specific focus on moraine formation, preservation, and ease of identification (topoclimatic controls are not considered). Published literature indicates that the importance of topography in regulating moraine distribution varies spatially, temporally, and as a function of the ice mass type responsible for moraine deposition. In particular, in the case of ice sheets and ice caps ( > 1000 km2), one potentially important topographic control on where in a landscape moraines are deposited is erosional feedback, whereby subglacial erosion causes ice masses to become less extensive over successive glacial cycles. For the marine-terminating outlets of such ice masses, fjord geometry also exerts a strong control on where moraines are deposited, promoting their deposition in proximity to valley narrowings, bends, bifurcations, where basins are shallow, and/or in the vicinity of topographic bumps. Moraines formed at the margins of ice sheets and ice caps are likely to be large and readily identifiable in the modern landscape. In the case of icefields and valley glaciers (10–1000 km2), erosional feedback may well play some role in regulating where moraines are deposited, but other factors, including variations in accumulation area topography and the propensity for moraines to form at topographic pinning points, are also likely to be important. This is particularly relevant where land-terminating glaciers extend into piedmont zones (unconfined plains, adjacent to mountain ranges) where large and readily identifiable moraines can be deposited. In the case of cirque glaciers (< 10 km2), erosional feedback is less important, but factors such as topographic controls on the accumulation of redistributed snow and ice and the availability of surface debris, regulate glacier dimensions and thereby determine where moraines are deposited. In such cases, moraines are likely to be small and particularly susceptible to post-depositional modification, sometimes making them difficult to identify in the modern landscape. Based on this review, we suggest that, despite often being difficult to identify, quantify, and mitigate, topographic controls on moraine distribution should be explicitly considered when reconstructing the dimensions of palaeoglaciers and that moraines should be judiciously chosen before being used as indirect proxies for palaeoclimate (i.e., palaeoclimatic inferences should only be drawn from moraines when topographic controls on moraine distribution are considered insignificant).

Modélisation de répartition d’espèces aviaires et de feux en forêt boréale du Québec dans un contexte de changement climatique

Les changements climatiques prennent une importance grandissante dans l’étude des phénomènes spatiaux à grande échelle. Plusieurs experts affirment que les changements climatiques seront un des principaux moteurs de changement écologique dans les prochaines décennies et que leurs conséquences seront inévitables. Ces changements se manifesteront sur le milieu physique par la fonte des calottes glaciaires, le dégel du pergélisol, l’instabilité des versants montagneux en zone de pergélisol, l’augmentation de l’intensité, de la sévérité et de la fréquence des événements climatiques extrêmes tels les feux de forêt. Les changements climatiques se manifesteront aussi sur le milieu biologique, tel la modification de la durée de la saison végétative, l’augmentation des espèces exotiques invasives et les changements dans la distribution en espèces vivantes. Deux aspects sont couverts par cette étude : 1) les changements dans la répartition spatiale de 39 espèces d’oiseaux et 2) les modifications dans les patrons spatiaux des feux, en forêt boréale québécoise, tous deux dans l’horizon climatique de 2100. Une approche de modélisation statistique démontre que la répartition spatiale des oiseaux de la forêt boréale est fortement liée à des variables bioclimatiques (R2adj = 0.53). Ces résultats permettent d’effectuer des modélisations bioclimatiques pour le gros-bec errant et la mésange à tête noire quivoient une augmentation de la limite nordique de distribution de l’espèce suivant l’intensité du réchauffement climatique. Finalement, une modélisation spatialement explicite par automate cellulaire permet de démontrer comment les changements climatiques induiront une augmentation dans la fréquence de feux de forêt et dans la superficie brûlée en forêt boréale du Québec.