La modélisation numérique des écoulements de glace est indispensable pour prédire l’évolution des calottes polaires suite au réchauffement climatique. De récentes études ont souligné l’importance des modèles d’écoulement dits d’ordre supérieur voir même de Stokes au lieu de la traditionnelle approximation de couche mince dont les hypothèses ne sont pas valables dans certaines zones critiques mais à l’étendue limitée. Cependant, ces modèles d’ordre supérieur sont difficiles à utiliser à l’échelle d’un continent en raison de leurs temps de calculs prohibitifs. Ce travail de thèse propose une nouvelle technique qui permet de réduire les temps de calculs tout en maximisant la précision des modèles. Plusieurs modèles d’écoulement de glace de complexité variables ont été mis en place dans ISSM (Ice Sheet System Model), un code élément fini massivement parallèle développé par le Jet Propulsion Laboratory. L’analyse et la comparaison des différents modèles, à la fois sur des cas théoriques et réels, montrent que l’utilisation des modéles les plus complets est principalement nécessaire au voisinage de la zone d’échouage, transition entre les parties flottantes et posées de la glace, mais aussi que des modèles plus simples peuvent être utilisés sur la majeure partie des glaciers. Coupler différents modèles présente donc un avantage significatif en terme de temps de calcul mais aussi d’amélioration de la physique utilisées dans les modèles. Plusieurs méthodes de couplage de modèles existent et sont présentées dans ce manuscrit. Une nouvelle technique, dite de tuilage, particulièrement adaptée au couplage de modèles d’écoulement de glace est décrite ici : son principe repose sur la superposition et le raccordement de plusieurs modèles mécaniques. Une analyse mathématique est effectuée afin de définir les conditions d’utilisation de cette méthode de tuilage. Le traitement du couplage entre un modèle de Stokes et des modèles simplifiés, pour lesquels le calcul des vitesses horizontales et verticales est découplé, est ensuite présenté. Cette technique a été mise en place dans ISSM afin de pouvoir créer des modèles hybrides combinant plusieurs modèles d’écoulement de complexité variable. Après avoir été validée sur des cas synthétiques, cette technique est utilisée sur des glaciers réels comme Pine Island Glacier, dans l’Antarctique de l’Ouest, afin d’illustrer sa pertinence. Les modèles hybrides ont le potentiel d’améliorer la précision des résultats en combinant différents modèles mécaniques, utilisés chacun dans les zones où leurs approximations sont valides, tout en réduisant les temps de calcul et en étant compatibles avec les ressources informatiques actuelles. / Ice flow numerical models are essential for predicting the evolution of ice sheets in a warming climate. Recent research emphasizes the need for higher-order and even full-Stokes flow models instead of the traditional Shallow-Ice Approximation whose assumptions are not valid in certain critical but spatially limited areas. These higher-order models are however computationally intensive and difficult to use at the continental scale. The purpose of this work, therefore, is to develop a new technique that reduces the computational cost of ice flow models while maximizing their accuracy. To this end, several ice flow models of varying order of complexity have been implemented in the Ice Sheet System Model, a massively parallelized finite element software developed at the Jet Propulsion Laboratory. Analysis and comparison of model results on both synthetic and real geometries shows that sophisticated models are only needed in the grounding line area, transition between grounded and floating ice, whereas simpler models yield accurate results in most of the model domain. There is therefore a strong need for coupling such models in order to balance computational cost and physical accuracy. Several techniques and frameworks dedicated to model coupling already exist and are investigated. A new technique adapted to the specificities of ice flow models is developed: the Tiling method, a multi-model computation strategy based on the superposition and linking of different numerical models. A mathematical analysis of a mixed Tiling formulation is first performed to define the conditions of application. The treatment of the junction between full-Stokes and simpler models that decouple horizontal and vertical equation is then elaborated in order to rigorously combine all velocity components. This method is finally implemented in the Ice Sheet System Model to design hybrid models that combine several ice flow approximations of varying order of complexity. Following a validation on synthetic geometries, this method is applied to real cases, such as Pine Island Glacier, in West Antarctica, to illustrate its relevance. Hybrid models have the potential to significantly improve physical accuracy by combining models in their domain of validity, while preserving the computational cost and being compatible with the actual computational resources.

This thesis is concerned with layered reflections observed in the Antarctic and Greenland ice sheets during radio-echo sounding. It describes in detail layering seen with 60 and 300 MHz equipment _in the Antarctic ice sheet during three field seasons betwe~n 1974 and 1979, and discusses the effects of glaciological and equipment factors, particularly radio pulse length. The effects of changing pulse length and accumulation rate on layer echo separation are analysed, and the importance of closely spaced groups of reflectors assessed. Reflection coefficient vs depth profiles are presented for layering observed at nearly twenty sites in the Antarctic and Greenland, and are used to demonstrate the existence of two separate reflection mechanisms: changes in (a) ice density, and (b) loss tangent. Supporting density and conductivity data from ice cores are summarised. It is concluded that most layer echoes arise from layers of ice containing acidic impurities of volcanic origin~ Layer reflectivity variations_are observed over distances of tens of kilometres. Short-period fading is also observed, from which estimates of layer reflector roughness are made using the autocorrelation function and variance of the received power. The roughness estimates are shown to be consistent with a depositional origin for the surfaces. Relevant theory is summarised and a procedure developed for the remote estimation of elevated acidity levels in ice from radio-echo sounding. Such estimates are shown to compare well with direct ice core measurements. The method is used to present estimated elevated acidity profiles for the Antarctic (to ~100 kaBP) and Greenland (to ~30 kaBP), which are interpreted in terms of variations in the input of volcanogenic acid impurities to the ice sheets. The use of layering as isochronous horizons in ice flow studies is discussed in the light of new measurements, with particular emphasis on t he zone close to bedrock.

There is increasing recognition that the hydrology of the Greenland Ice Sheet plays an important role in the dynamics and therefore mass balance of the ice sheet. Understanding the hydrology of the ice sheet and being able to predict its future behaviour is therefore a key aspect of glaciological research. To date, the ice sheet’s hydrology has tended to be inferred from the analysis of surface velocity measurements, or modelled in a theoretical, idealised way. This study focuses on the development of a high spatial (100 m) and temporal (1 hour) resolution, physically based, time-dependent hydrological model which is applied to the ~2,300 km2 Paakitsoq region, West Greenland, and is driven, calibrated, and evaluated using measured data. The model consists of three components. First, net runoff is calculated across the ice sheet from a distributed, surface energy- balance melt model coupled to a subsurface model, which calculates changes in temperature, density and water content in the snow, firn and upper-ice layers, and hence refreezing. The model is calibrated by adjusting key parameter values to minimize the error between modelled output and surface height and albedo measurements from the three Greenland Climate Network (GC-Net) stations, JAR 1, JAR 2 and Swiss Camp. Model performance is evaluated in two ways by comparing: i) modelled snow and ice distribution with that derived from Landsat-7 ETM+ satellite imagery using Normalised Difference Snow Index (NDSI) classification and supervised image thresholding; and ii) modelled albedo with that retrieved from the Moderate- resolution Imaging Spectroradiometer (MODIS) sensor MOD10A1 product. Second, a surface routing / lake filling model takes the time-series of calculated net runoff over the ice sheet and calculates flow paths and water velocities over the snow / ice covered surface, routing the water into ‘open’ moulins or into topographic depressions which can fill to form supraglacial lakes. This model component is calibrated against field measurements of a filling lake in the study area made during June 2011. Supraglacial lakes are able to drain by a simulated hydrofracture mechanism if they reach a critical volume. Once water is at the ice / bed interface, discharge and hydraulic head within subglacial drainage pathways are modelled using the third model component. This consists of an adaptation of a component (EXTRAN) of the U.S. Environmental Protection Agency Storm Water Management Model (SWMM), modified to allow for enlargement and closure of ice-walled conduits. The model is used to identify how the subglacial hydrological system evolves in space and time in response to varying surface water inputs due to melt and lake drainage events, driven ultimately by climate data. A key output from the model is the spatially and temporally varying water pressures which are of interest in helping to explain patterns of surface velocity and uplift found by others, and will ultimately be of interest for driving ice dynamics models.

Active microwave sensors are known to provide valuable information regarding snow and ice surfaces in the polar regions, where darkness and cloud cover prevail. Here, data collected in the Arctic by a Ku-band microwave radar altimeter, designed and constructed in the UK, are analysed. The two main components of this study comprise data gathered in the East Greenland Sea marginal ice zone and over two Svalbard ice caps. A systematic treatment is made of the electromagnetic properties of snow and ice at 13.81 GHz, and the differences between various polar surface media are highlighted. Theoretical and empirical models are presented which enable calculation of the relevant dielectric and scattering properties of snow and ice layers. Parametric studies are undertaken to give insight into the range of scattering conditions likely to be encountered by a radar altimeter in the regions investigated. Examples of altimetric data and results of their analysis are presented, demonstrating the effects of different ice types and terrain upon incident altimeter pulses. Waveforms are characterised by their shape, and certain forms are linked with particular physical properties of the surface. To this a variety of supporting information is added in order to verify and validate interpretations of these results. Algorithms are proposed which enable geophysical information to be derived from altimetric data.

The growth of the Weddell Sea ice cover in 1980 is examined, using Nimbus-7 satellite Scanning Multi-channel Microwave Radiometer (SMMR) data in combination with data (positional, air temperature and atmospheric pressure) from 4 Nimbus-tracked drifting buoys. Ice concentrations are retrieved from the SMMR data by applying a cluster analysis algorithm developed by J. Comiso of NASA. Analyses of computed differential kinematic parameters (DKP)s of the buoy array offer insight into the complex mesoscale behaviour of the underlying Weddell Gyre. High frequency divergence, convergence and deformation events isolated in the DKP results, and driven largely by the regular passage of cyclones, are related to changes in ice concentration observed in the SMMR data. The profound role of the Antarctic Peninsula in influencing both atmospheric and oceanic circulation (and thus ice formation, drift and eventual decay) in the region is evaluated. Possible relationships between buoy drift in the inner pack and ice edge advance are examined, yielding information on the relative importance of ice growth in open water within the pack and that at the ice edge. After an introductory chapter, Chapter 2 describes the physical setting of the Weddell Sea. Comparisons are drawn both with other sectors of the Southern Ocean and the Arctic, emphasizing the uniqueness of the region not only in terms of its climate and oceanography but also its sea ice cover. Chapter 3 traces the evolution of passive microwave remote sensing from space as a tool for monitoring Antarctic sea ice extent and concentration; the relative merits and disadvantages of these techniques are evaluated. Chapter 4 concentrates on the use of SMMR data. Detailed comparisons are made of algorithms available for the extraction of ice concentrations from the raw brightness temperature data. The choice of algorithm used is justified. Chapter 5 is largely concerned with the analysis of the buoy data, and the kinematic behaviour of the array as a unit. These results are combined with the SMMR data in Chapter 6 to identify distinct dynamic zones and meridional advective sectors, and to compare the behaviour of the inner pack with that of the unconstrained ice edge. The evolution of a high concentration core within the unique perennial sea ice zone hugging the east coast of the Peninsula, which persists throughout the period of study is unusual enough to merit a separate sub-section. Conclusions are drawn in chapter 7.

been used to investigate the form and topography of the Antarctic ice sheet and to relate these to the physical processes of ice flow and basal conditions. Topographic roughness typically increases towards the thin ice of coastal reg ions as surface undulation wav el eng ths decrease and amplitudes increase. Temperature and velocity variations also have significant effects. The coastal zone is punctuated by embayments of severe topography immediately inland of outlet glaciers. This topographic variability has been summarized in a statistical model for the purposes of simulating satellite radar altimeter waveforms. Consideration of the relationship between bedrock and surface profiles has shown that ice temperature is a major influence on the response of the surface to bedrock irregularities. Regional subglacial water layers may al so have an important effect on surface topography. A re-analysis of models of longitudinal stress grad i ents suggests that er ystal fabrics favouring faster flow develop with distance from ice divides and that the relative depth of the zone of maximum shear fluctuates in response to topographic and glaciological constraints. Driving stress patterns have been associated with characteristic glaciological regimes and have suggested a qualitative difference between outlet glaciers and ice streams. The transition to high velocity flow in outlet glaciers has been shown to be triggered ab ruptly in response to subglac ial fjord heads. The dependence of fast flow on subglacial topography indicates a significant stabilizing effect on discharge from ice sheets and suggests that surge behaviour is unlikely within existing ice sheet outlets. The onset of basal sliding at the head of subglac ial fjords suggests a mechanism for the production of overdeepened fjords and steep headwalls through concentrated erosion. This may help in the reconstruction of the dynamics of former ice sheets. Some West Antarctic ice streams do not exhibit this rapid transition in behaviour.

Glacier facies represent distinct regions of a glacier surface characterized by near surface structure and density that develop as a function of spatial variations in surface melt and accumulation. In post freeze-up (autumn) synthetic aperture radar (SAR) satellite imagery, the glacier ice zone and dry snow zone have a relatively low backscatter due to the greater penetration of the radar signal into the surface. Conversely, the saturation and percolation zones are identifiable based on their high backscatter due to the presence of ice lenses and pipes acting as efficient scatterers. In this study, EnviSat ASAR imagery is used to monitor the progression of facies zones across Devon Ice Cap (DIC) from 2004 to 2011. This data is validated against in situ surface temperatures, mass balance data, and ground penetrating radar surveys from the northwest sector of DIC. Based on calibrated (sigma nought) EnviSat ASAR backscatter values, imagery from autumn 2004 to 2011 shows the disappearance of the ‘pseudo’ dry snow zone at high elevations, the migration of the glacier and superimposed ice zones to higher elevations, and reduction in area of the saturation/percolation zone. In 2011, the glacier and superimposed ice zone were at their largest extent, occupying 92% of the ice cap, leaving the saturation/percolation zone at 8% of the total area. This is indicative of anomalously high summer melt and strongly negative mass balance conditions on DIC, which results in the infilling of pore space in the exposed firn and consequent densification of the ice cap at higher elevations.

There has been a marked increase in melt season length over the past two decades on glaciers and ice caps within Canada’s Queen Elizabeth Islands (QEI). Prior to the year ~2000 land ice was in a state of slightly negative mass balance (-11 +/- 11.5 Gt yr⁻¹ over 1958-1995), but recent GRACE measurements suggest that mass losses averaged -33 +/- 5 Gt yr⁻¹ between 2003-2015. These losses have primarily been attributed to meltwater runoff, making the QEI one of the largest recent contributors to sea level rise outside of the ice sheets. Despite these losses, there is a lack of information concerning how a warming climate is affecting small (<1 km²) ice bodies, which are considered sensitive indicators of change due to their short response time. In this study, historical and contemporary aerial photographs, high resolution optical satellite imagery, and ground penetrating radar (GPR) surveys are used to determine area, thickness, mass and volume changes of Adams Icefield within Expedition Fiord, Axel Heiberg Island, Nunavut, over the past seven decades (1948/59-2019). Area changes are determined from a comparison of air photos acquired in 1948/59 with satellite images acquired since 1979. Contemporary (2001, 2012, 2019) digital elevation models (DEMs) were either collected or created from stereo satellite images, and via aerial photo surveys using Structure from Motion photogrammetry. DEM of Difference maps calculated from these DEMs provide volume and mass changes. Results illustrate a steady reduction in glacier area, thickness, and volume prior to the year ~2000, followed by a rapid increase in losses over the past two decades. As a result, Adams Icefield is now rapidly declining and is likely to completely disappear early in the twenty-second century.

The Coast Mountains in British Columbia contain over 10,000 km2 of glacial ice. While these glaciers have lost significant mass since the Little Ice Age (LIA; around 300 years before present), the melting rate has significantly increased over the past decade, likely due to the effects of climate change. The purpose of this study was to develop an approach to quantifying the isostatic response to LIA glacier change and investigate how it can further our understanding of the Earth’s rheology through GIA modelling. The Coast Mountains in southwestern British Columbia were chosen due to their significant ice mass loss since the LIA, their location in a tectonically active region, which includes a volcanic arc, and the presence of information of vertical land motion. The GIA models in this study use a wide range of Earth rheological parameters that are then constrained through comparison to observations of vertical land motion in the region. The study used available Global Navigation Satellite System (GNSS) vertical velocity data as the observable from seven GNSS sites in southwestern BC, using a combination of Western Canada Deformation Array (WCDA) and British Columbia Active Control System (BCACS) GNSS stations. Raw data were analyzed using the GIPSY 6.4 software following the Precise Point Positioning processing strategy. Two ice load histories were developed based on gridded estimates of present-day ice thicknesses in the region in order to simulate the change in the surface loading as the glacial ice mass fluctuates over time. Ice Load A used a simple uniform thickness change profile over 3 time-steps based on extrapolated modern melt rates. Ice Load B is more complex and utilized a published profile of glacier change through time basing the magnitude of volume changes on the volume-area scaling relationship with a range of coefficient values. This allowed for a range of ice change magnitudes to be tested. The Earth models used were spherically symmetric Preliminary Reference Earth Models (PREM). Their viscosity structure is based on VM5a for the transition zone and lower mantle, but with variable lithospheric thickness and asthenospheric viscosity. The goodness of fit for the modeled velocities were compared to the observed velocities using a normalized RMS (NRMS) statistic. Ice Load A models had a best fitting lithospheric thickness of 50 km and an asthenospheric viscosity of 2×1019 Pa s. For all variations of Ice Load B, the best fitting model parameters had lithospheric thicknesses ranging from 45 km to 55 km and asthenospheric viscosities between 6×1018 Pa s and 3×1019 Pa s. Corrected GNSS vertical velocity observations were tested to check the effects of interseismic vertical signal and assumed residual GIA from the Cordilleran Ice Sheet. However, the corrections did not improve the NRMS fit. Overall, the asthenospheric viscosity results from this study overlap with all the ranges found in the previous studies while lithospheric thicknesses agree with some past studies. The results of this study generally align with previous work and the current understanding of the Coast Mountains region and can inform a future round of sea-level projections for the region as ice mass loss continues in the Coast Mountains. This study serves to further refine constraints on Earth rheology and can be used to guide future work on GIA in the region. / Graduate

No description available.