Laser altimetry and satellite gravity surveys indicate that the St Elias Icefields are currently losing mass and are among the largest non-polar sea level contributors in the world. However, a poor understanding of glacier dynamics in the region is a major hurdle in evaluating regional variations in ice motion and the relationship between changing surface conditions and ice flux. This study combines in-situ dGPS measurements and advanced Radarsat-2 (RS-2) processing techniques to determine daily and seasonal ice velocities for the Kaskawulsh Glacier from summer 2009 to summer 2011. Three permanent dGPS stations were installed along the centreline of the glacier in 2009, with an additional permanent station on the South Arm in 2010. The Precise Point Positioning (PPP) method is used to process the dGPS data using high accuracy orbital reconstruction. RS-2 imagery was acquired on a 24-day cycle from January to March 2010, and from October to March 2010-2011 in a combination of ultra-fine and fine beam modes. Seasonal velocity regimes are readily identifiable in the dGPS results, with distinct variations in both horizontal velocity and vertical motion. The Spring Regime consists of an annual peak in horizontal velocity that corresponds closely with the onset of the melt season and progresses up-glacier, following the onset of melt at each station. The Summer Regime sees variable horizontal velocity and vertical uplift, superimposed on a long-term decline in motion. The Fall Regime sees a gradual slowing at all stations with little variation in horizontal velocity or vertical position. Rapid but short accelerations lasting up to 10 days were seen in the Winter regimes in both 2010 and 2011, occurring at various times throughout each regime. These events initiated at the Upper Station and progress down-glacier at propagation speeds up to 16,380 m day-1 and were accompanied by vertical uplift lasting for similar periods. Three velocity maps, one from the winter of 2010 and two from the fall/winter of 2011, produced from speckle tracking were validated by comparison with dGPS velocity, surface flow direction, and bedrock areas of zero motion, with an average velocity error of 2.0% and average difference in orientation of 4.3º.

Laser altimetry and satellite gravity surveys indicate that the St Elias Icefields are currently losing mass and are among the largest non-polar sea level contributors in the world. However, a poor understanding of glacier dynamics in the region is a major hurdle in evaluating regional variations in ice motion and the relationship between changing surface conditions and ice flux. This study combines in-situ dGPS measurements and advanced Radarsat-2 (RS-2) processing techniques to determine daily and seasonal ice velocities for the Kaskawulsh Glacier from summer 2009 to summer 2011. Three permanent dGPS stations were installed along the centreline of the glacier in 2009, with an additional permanent station on the South Arm in 2010. The Precise Point Positioning (PPP) method is used to process the dGPS data using high accuracy orbital reconstruction. RS-2 imagery was acquired on a 24-day cycle from January to March 2010, and from October to March 2010-2011 in a combination of ultra-fine and fine beam modes. Seasonal velocity regimes are readily identifiable in the dGPS results, with distinct variations in both horizontal velocity and vertical motion. The Spring Regime consists of an annual peak in horizontal velocity that corresponds closely with the onset of the melt season and progresses up-glacier, following the onset of melt at each station. The Summer Regime sees variable horizontal velocity and vertical uplift, superimposed on a long-term decline in motion. The Fall Regime sees a gradual slowing at all stations with little variation in horizontal velocity or vertical position. Rapid but short accelerations lasting up to 10 days were seen in the Winter regimes in both 2010 and 2011, occurring at various times throughout each regime. These events initiated at the Upper Station and progress down-glacier at propagation speeds up to 16,380 m day-1 and were accompanied by vertical uplift lasting for similar periods. Three velocity maps, one from the winter of 2010 and two from the fall/winter of 2011, produced from speckle tracking were validated by comparison with dGPS velocity, surface flow direction, and bedrock areas of zero motion, with an average velocity error of 2.0% and average difference in orientation of 4.3º.

A consequence of climate change is rising global sea levels, predicted to bring increased socio-economic and environmental impacts to coastal communities. The Greenland Ice Sheet has become a prominent contributor to rising sea levels, a consequence of the Arctic warming at twice the rate of the global average. Mass loss from the ice sheet is separated between changes in surface mass balance and ice discharge at marine terminating outlet glaciers, with the later dominating mass loss over the past fifty years. While advances in ice sheet modelling have provided greater clarity on Greenland’s future mass loss, there remains inefficiencies in modelling the response of outlet glaciers in Greenland’s fjords. This thesis aims to provide greater insight into behaviour of Kangerlussuaq Glacier, SE Greenland, by employing a 2D flowline model to understand the processes governing ice dynamics and to explore how the glacier may respond to a warming climate. Results indicate that the presence of a winter ice mélange is the principle dictator of Kangerlussuaq Glacier’s behaviour and likely protects against further retreat towards a reverse sloped section of bedrock. However, if such a retreat does materialise, then large overdeepenings in Kangerlussuaq Fjord raise the spectre of uncontrollable retreat and excessive mass loss.

<p>NR 20140805</p>

Despite their relatively small ice volume, mountain glaciers contributed nearly one third of global sea level rise since 2000, with one of the largest total mass loss rates (73 ± 17 Gt a-1) occurring in the Yukon-Alaska region. However, there is uncertainty surrounding how ice dynamics are being affected by such losses and whether glacier flow instabilities, such as surges, are changing in a warming climate. The St. Elias Mountains contain a major cluster of surge-type glaciers, yet a detailed analysis of their characteristics, including surge frequency, morphology, magnitude, and propensity over time has not been undertaken on a regional basis. This thesis presents a review of surging behaviour and an updated surge event inventory in the St. Elias Mountains, and quantifies the processes influencing both surging and non-surging glacier velocity variability using a variety of remote sensing and field measurements. An updated inventory of surge-type glaciers and observed surge events (1874-2023), compiled from existing inventories, recently published articles, and velocity analysis, is used to analyze the characteristics of surge-type glaciers and velocity patterns during surge events. The modern (1985-2023) trends in annual, winter and summer velocities of selected surge-type glaciers is then used to classify dynamic instability events into 4 categories. While 231 glaciers were classified as surge-type, only 42 were observed to have experienced rapid velocity events over the period 1985-2023, through either direct measurements or remote sensing observations. For glaciers with observed rapid velocity events, these predominantly fall into two categories: Alaskan-style surges with short active and quiescent phases, and glacier pulses, which are velocity accelerations that are limited in both magnitude and extent. An unnamed former tributary to Kluane Glacier underwent a dramatic surge from 2013-18. Using a combination of air photos, remote sensing and field observations, the characteristics and changes of ‘Little Kluane Glacier’ were reconstructed from the 1940s until 2021. While only the single full surge of 2013-18 was identified, it is likely that a partial surge of just the upper north arm occurred between 1963 and 1972. Repeat Digital Elevation Models (DEMs) and velocity profiles show that the recent surge initiated from the upper north arm accumulation area in 2013, which developed into a full surge of the main trunk from 2017-18. Terminus positions show long-term retreat from 1949-2017, followed by rapid advance of >2 km from May to September 2018, with surface velocities reaching a peak of ~3600 m a-1 in summer 2018 over the lower ablation area. This was likely enhanced by the drainage of supraglacial lakes and streams to the glacier bed through crevassing as the surge progressed. Changes in surface topography caused by initial mass movement, the resulting reorganization of the supraglacial hydrological system, and ponding of surface water, may drive a partial surge into a full surge, and therefore exert a direct control on glacier dynamics. In May 2016, Kaskawulsh Glacier underwent a dramatic proglacial hydrologic reorganization instigated by the rapid drainage of proglacial Slims Lake: as a result, water which previously drained north into Ä’äy Chú, (Slims River) toward Lhú’áán Män (Kluane Lake), was redirected south into Kaskawulsh River, eventually flowing into the Gulf of Alaska. A long-term (up to ∼120 year) record of terminus retreat, thinning and surface velocities from in-situ and remote sensing observations is used to determine the impact of this reorganization on glacier dynamics. After an initial deceleration during the late 1990s, terminus velocities increased at a rate of 3 m a-2 from 2000-12, while the area of proglacial Slims Lake increased simultaneously. The rapid drainage of the lake substantially altered the velocity profile of the adjacent glacier, decreasing annual velocities by 48% within 3 km of the terminus between 2015 and 2021, at an average rate of ∼12.5 m a-2. A key cause of the rapid drop in glacier motion was a reduction in flotation of the lower part of the terminus after lake drainage. This has important implications for glacier dynamics and provides one of the first assessments of the impacts of a rapid proglacial lake drainage event on local terminus velocities. The results of this study provide an examination of factors controlling glacier dynamics, as well as the characteristics of rapid glacier velocity events, in the St. Elias Mountains. This provides insights into the behaviour of mountain glaciers, how they are changing in a warming climate, controls on glacier surging, and the hazards they may pose for downstream communities, which are particularly vulnerable to disturbances.

This work aims to increase our understanding of the nature of large scale features of sea ice from a dynamics point of view.Sea ice plays an important part in the exchange of heat and humidity between sea and air and thus is an important component of the climate system. Its physical presence also directly impacts the various forms of life such as diatoms, polar bears and humans alike.The dynamics of sea ice affect both weather and climate, through the large scale drift in the Arctic from the Siberian coast towards Fram Strait, through creation of cracks in the ice called leads or polynyas, and through ridging and other mechanical deformations of ice floes.In this work, we have focused on modelling of ridged ice for a number of reasons. Direct observations of the internal ice state is very difficult to perform and in general, observations of sea ice are either sparse or of limited information density. Ridged ice can be seen as the memory of high ice stress events, giving us a view on these highly dynamic events. Ridging is of major importance for the ice thickness distribution, as the thickest ice can only be formed through mechanical processes. Further, ridged ice is of direct interest for anyone conducting shipping through seasonal or perennial ice covered seas as it can form impenetrable barriers or in extreme even cases crush a ship caught within the ice pack. To this end, a multi-category sea ice model, the HELsinki Multi category Ice model (HELMI), was implemented into the Rossby Centre Ocean model (RCO). HELMI has explicit formulations for ridged and rafted ice, as well as sub-grid scale ice thickness distribution (a feature shared with other multi category models) and an ice strength based on energetics. These features give RCO better representation of sub-grid scale physics and gives us the possibility to study the deformed ice in detail. In paper I we look at the change in behaviour in the Arctic as the ice becomes more mobile, leading to a slight increase in modelled ridged ice volume in the central Arctic, despite a general trend of a decreasing ice cover.Paper II takes us to the Baltic Sea and the possibilities of modelling ridge ice concentration with a statistical model.In Paper III we investigate how the diminishing ice cover in future scenarios affects the biological activity in the Baltic Sea.Finally Paper IV investigates how the ice stress and the internal ice force can be interpreted in terms of ice compression on the ship scale. / <p>At the time of the doctoral defence the following paper was unpublished and had a status as follows: Paper 4: Manuscript</p>

Lake Vostok, East Antarctica, represents an extensive water surface at the base of the ice sheet. Snow, ice and atmospheric pressure loads applied anywhere within the lake area produce a hydrostatic response, involving deformations of the ice surface, ice–water interface and particle horizons. A modelling scheme is developed to derive height changes of these surfaces for a given load pattern. It is applied to a series of load scenarios, and predictions based on load fields derived from a regional climate model are compared to observational datasets. Our results show that surface height changes due to snow-buildup anomalies are damped over the lake area, reducing the spatial standard deviation by one-third. The response to air pressure variations, in turn, adds surface height variability. Atmospheric pressure loads may produce height changes of up to 4 cm at daily resolution, but decay rapidly with integration time. The hydrostatic load response has no significant impact neither on ICESat laser campaign biases determined over the lake area nor on vertical particle movements derived from GNSS observations.

<p>This thesis contributes to the understanding of glacier response to climate change by ice dynamical studies on Storglaciären, Sweden, and Bonnevie-Svendsenbreen, Kibergbreen and Plogbreen in Dronning Maud Land, Antarctica. Ice surface velocities, ice geometry and temperature information is fed through a force budget model to calculate ice mass outflux of these glacial systems via three-dimensional stress distributions for a flux-gate. </p><p>Field data were collected through repeated DGPS and GPR observations on Storglaciären between July 2000 to September 2001 and on Kibergbreen and Plobreen during the SWEDARP 2002/03 expedition to Antarctica. The work was strongly supported by remotely-sensed information.</p><p>The results from Storglaciären show a strength in the force budget model to discern both spatial and temporal variability in ice dynamical patterns. It highlights the influence of seasonality and bedrock topography upon glacier flow. A modeling experiment on Bonnevie-Svendsenbreen suggested that ice temperature increases substantially under conditions of high stress (≥0.4 MPa) due to strain-heating. This provides a positive feedback loop, increasing ice deformation, as long as it overcomes the advection of cool ice from the surface. These results explain, to some extent, the mechanism behind fast flowing ice streams. Mass flux caclulations from Bonnevie-Svendsenbreen suggest that the outflux given from force budget calculations can be used as a gauge for influx assuming steady state conditions. Plogbreen receives an influx of 0.48±0.1 km³ a^-1 and expedites a discharge volume of 0.55±0.05 km³ a^-1. This indicative negative mass balance is explained by a falling trend in upstream accumulation and the recent rise in global sea level, as it is likely to induce glacier acceleration due to a reduction in resistive forces at the site of the gate. This result is comparable with other Antarctic studies reporting negative mass balances, e.g. from WAIS, as caused by changes in the global atmospheric circulation pattern.</p>

This thesis contributes to the understanding of glacier response to climate change by ice dynamical studies on Storglaciären, Sweden, and Bonnevie-Svendsenbreen, Kibergbreen and Plogbreen in Dronning Maud Land, Antarctica. Ice surface velocities, ice geometry and temperature information is fed through a force budget model to calculate ice mass outflux of these glacial systems via three-dimensional stress distributions for a flux-gate. Field data were collected through repeated DGPS and GPR observations on Storglaciären between July 2000 to September 2001 and on Kibergbreen and Plobreen during the SWEDARP 2002/03 expedition to Antarctica. The work was strongly supported by remotely-sensed information. The results from Storglaciären show a strength in the force budget model to discern both spatial and temporal variability in ice dynamical patterns. It highlights the influence of seasonality and bedrock topography upon glacier flow. A modeling experiment on Bonnevie-Svendsenbreen suggested that ice temperature increases substantially under conditions of high stress (≥0.4 MPa) due to strain-heating. This provides a positive feedback loop, increasing ice deformation, as long as it overcomes the advection of cool ice from the surface. These results explain, to some extent, the mechanism behind fast flowing ice streams. Mass flux caclulations from Bonnevie-Svendsenbreen suggest that the outflux given from force budget calculations can be used as a gauge for influx assuming steady state conditions. Plogbreen receives an influx of 0.48±0.1 km3 a-1 and expedites a discharge volume of 0.55±0.05 km3 a-1. This indicative negative mass balance is explained by a falling trend in upstream accumulation and the recent rise in global sea level, as it is likely to induce glacier acceleration due to a reduction in resistive forces at the site of the gate. This result is comparable with other Antarctic studies reporting negative mass balances, e.g. from WAIS, as caused by changes in the global atmospheric circulation pattern.

Sea ice plays a key role in the global climate system. Indeed, through the albedo effect it reflects significant solar radiation away from the oceans, while it also plays a key role in the momentum and heat transfer between the atmosphere and ocean by acting as an insulating layer between the two. Furthermore, as more sea ice melts due to climate change, additional fresh water is released into the upper oceans, affecting the global circulation of the ocean as a whole. While there has been significant effort in recent decades, the ability to simulate sea ice has lagged behind other components of the climate system and most Earth System Models fail to capture the observed losses of Arctic sea ice, which is largely attributed to our inability to resolve sea ice dynamics. The most widely accepted model for sea ice dynamics is the Viscous- Plastic (VP) rheology, which leads to a very non-linear set of partial differential equations that are known to be intrinsically difficult to solve numerically. This work builds on recent advances in solving these equations with a Jacobian-Free Newton- Krylov (JFNK) solver. We present an improved JFNK solver, where a fully second order discretization is achieved via the Crank Nicolson scheme and consistency is improved via a novel approach to the rheology term. More importantly, we present a significant improvement to the Jacobian approximation used in the Newton iterations, and partially form the action of the matrix by expressing the linear and nearly linear terms in closed form and approximating the remaining highly non-linear term with a second order approximation of its Gateaux derivative. This is in contrast with the previous approach which used a first order approximation for the Gateaux derivative of the whole functional. Numerical tests on synthetic equations confirm the theoretical convergence rate and demonstrate the drastic improvements seen by using a second order approximation in the Gateaux derivative. To produce a fast and efficient solver for VP sea ice dynamics, the improved JFNK solver is then coupled with a non- oscillatory, central differencing scheme for transporting sea ice as well as a novel method for tracking the ice domain using a level set method. Two idealized test cases are then presented and simulation results discussed, demonstrating the solver’s ability to efficiently produce Viscous-Plastic, physically motivated solutions. / Graduate