No description available.

No description available.

Au-delà de leur rôle d’icône du changement climatique, les glaciers de montagne sont une composante essentielle de notre planète. Ils sont, de plus, de véritables « climat-mètres » naturels. Malgré leur faible superficie (0.5% des terres émergées), les glaciers de montagne contribuent à hauteur de 30% à la hausse du niveau des mers. Dans certaines régions, ils constituent de véritables enjeux quant à l’eau potable, l’agriculture, la production hydroélectrique ou les aléas glaciaires. Peu sont en revanche instrumentés (<0.0025%) et leurs fluctuations à l’échelle de régions entières sont mal connues.Grâce au développement de capteurs satellitaires à haute résolution spatiale (métrique à décamétrique), le développement de méthodes automatisées permet aujourd’hui d’augmenter considérablement le nombre de glaciers observés. Après avoir dressé un état de l’art des méthodes existantes et identifié les verrous méthodologiques, nous avons développé deux méthodes en particulier.La première se base sur la détection automatique de l’altitude de la limite glace/neige (i.e. ligne de neige) à la surface du glacier, à partir d’images satellites optiques. Cette altitude, lorsque mesurée à la fin de l’été, est un marqueur du changement de masse à la surface du glacier ayant eu lieu au cours de l’année (appelé bilan de masse de surface). Cette approche nous a permis d’estimer le bilan de masse de surface annuel de 239 glaciers dans les Alpes européennes et de 82 glaciers dans les Andes tropicales pour la période 2000-2016 et 2000-2018, respectivement. La perte moyenne annuelle observée est de -0.74 et de -1.29 m équivalent eau par an pour les deux régions respectivement. A notre connaissance, cette approche a permis d’établir le premier jeu de données de bilans de masse de surface annuels pour des glaciers individuels à échelle régionale à partir d’images satellites optiques. Une dépendance du bilan de masse de surface moyen par glacier à des caractères morpho-topographiques (e.g. pente, altitude médiane …) a été observée, où plus les glaciers sont pentus et hauts en altitude, moins leur perte de masse est importante. Une comparaison avec des mesures in situ dans les Alpes Européennes révèle une surestimation de la perte de masse par ces dernières si on les extrapole spatialement, notamment à cause de la faible représentation de glaciers à forte pente (>20°) dans les mesures in situ. Notre étude sur les Alpes Européennes a de plus permis d’identifier une variabilité interannuelle hétérogène sur cette région, en partie expliquée par des contextes climatiques différents grâce à l’utilisation de données issues de ré-analyses.Le développement d’une autre méthode a permis, à partir de l’analyse de carte d’albédo issues du capteur MODIS, de caractériser le bilan de masse de surface annuel et estival de 30 glaciers dans les Alpes françaises. Cette étude ouvre la porte à l’utilisation de cette méthode pour l’analyse du bilan annuel et saisonnier à l’échelle régionale.Ce travail a permis, à travers des applications dans différentes régions englacées, de développer et valider des méthodes capables, à partir d’images satellites optiques, d’estimer le bilan de masse de surface annuel et saisonnier de glaciers de montagne à l’échelle de régions entières. Ces estimations peuvent ensuite être utilisées pour : (1) étudier l’impact du climat local sur les glaciers de montagne ; (2) d’investiguer de possibles conditions météorologiques favorisant les fluctuations observées ; (3) calibrer et valider les modèles glacio-hydrologiques utilisés pour estimer les contributions actuelles et futures des glaciers de montagne au fonctionnement hydrologique des bassins versants et à l'élévation du niveau des mers. / Beyond their iconic role of climate change, mountain glaciers can be considered as Earth’ essential component and natural “climate-meter”. Despite their small spatial coverage (0.5% of emerged land), mountain glaciers contribute as high as 30% of the observed sea-level rise. In some regions, they are considered as essential issues because of their importance in terms of potable water, agriculture, hydroelectricity or natural hazards. A small share is however monitored in situ (<0.0025%) and their fluctuations at regional scale are poorly known.Thanks to the development of high spatial resolution satellite sensors (metric to decametric), new methods are today available to significantly increase the number of monitored glaciers. After a state of the art of the existing methods and an identification of the limitations, we focused our attention on the development of two methods.The first one is based on the automatic detection of the snow/ice interface altitude (i.e. snowline) at the glacier surface from optical satellite images. This altitude, when estimated at the end of summer, is a proxy of the annual glacier-wide mass change at the glacier surface (called surface mass balance, SMB). Using this approach, we estimated the annual SMBs of 239 glaciers in the European Alps and 82 glaciers in the tropical Andes for the period 2000-2016 and 2000-2018, respectively. The mean mass loss are -0.74 and -1.29 m water equivalent per year for the two regions, respectively. This approach allowed to derive the first dataset of annual SMBs for individual glaciers at regional scale from optical remote sensing. We found significant relationships between the computed SMBs and the glacier morpho-topographic features (e.g. slope, median altitude, …), with steeper and higher glaciers, experiencing less mass losses. Comparison with in situ monitored SMBs revealed an overestimation of mass losses from in situ estimates, due to a low representativeness of steep glaciers (>20°) in the in situ datasets. Our study also revealed heterogeneous inter-annual variability across the European Alps, partially explained by the climatic context of the studied sub-regions, thanks to the analysis of climate reanalysis data.We developed a second method to derive the annual and summer SMBs from albedo maps, computed from MODIS images. With an application on 30 glaciers in the French Alps, this work opened the way toward a regional application of this method, in order to estimate both annual and summer SMBs.By performing regional applications on different glacierized regions, we developed and validated methods capable of deriving the annual and summer SMBs of individual mountain glaciers at regional scale, from optical remote sensing data. These data could then be used to (1) assess the impact of peculiar climatic conditions onto mountain glaciers; (2) investigate possible meteorological conditions driving the documented glacier fluctuations; (3) calibrate and validate glacio-hydrological models used to estimate the current and future contributions of mountain glaciers to the hydrological functioning of mountain catchments and to sea level rise.

Since the mid-20th century, global average temperatures have dramatically risen mostly due to the increasing amount of greenhouse gas emissions in the atmosphere. The effects of this recent global warming are already evident and could be exacerbated in the near future if no real action is taken. Recent ice loss in West Antarctica, monitored by satellite measurements and other techniques, gives cause for concern in such a warming world. A major part of this loss has been driven by warm water masses penetrating underneath the ice shelves in this region. This has led to a flow acceleration of the inland outlet glaciers and a greater discharge of ice to the ocean. The actual resulting contribution of West Antarctica to sea-level rise is estimated to be around 0.2 mm per year between 1992 and 2011, i.e. about one third of the ice-sheet contribution (Antarctica and Greenland), and is expected to increase in the near future.<p><p>In this thesis, we first clearly demonstrate that modeling grounding-line (the boundary between grounded and floating ice) migration depends on both the numerical approach and the physical approximation of the ice-sheet model used. Ice-sheet models prescribing the ice flux at the grounding line and using appropriate physical level and numerical approach converge to the same steady-state grounding-line position irrespective of the grid size used. However, the transient behavior of those models is less accurate than other models and leads to an overestimated grounding-line discharge. Therefore, they need to be used with particular attention on short time scales. Furthermore, the non-inclusion of vertical shear stress in those models increases the effective viscosity and gives steady-state grounding-line positions further downstream when compared to full-Stokes models.<p><p>The second major finding of this thesis is the high control of geometry (glacier width and bedrock topography) on Thwaites Glacier, one of the fastest-flowing outlet glaciers in West Antarctica. A flowline finite-difference Shallow-Shelf Approximation (SSA) model is applied to the glacier and shows that ice-flow convergence (through width parameterization) slows down the grounding-line retreat when compared to simulations where the width is constant. A new buttressing parameterization is also tested on the glacier and permits a better understanding of this effect. Finally, the three-dimensional version of the model above is applied to Thwaites Glacier and highlights the strong control of lateral variations in bedrock topography on grounding-line migration./Depuis le milieu du 20e siècle, les températures moyennes globales ont fortement augmenté principalement à cause de l'augmentation des émissions de gaz à effet de serre d'origine humaine. Les effets de ce réchauffement global récent sont déjà détectables et pourraient s'accentuer dans un futur proche si aucune mesure réelle n'est prise. La perte récente de glace en Antarctique de l'Ouest, enregistrée par mesures satellites et d'autres techniques, est préoccupante dans un monde qui se réchauffe. Une grande partie de cette perte de glace est due à la pénétration de masses d'eau chaude sous les plateformes de glace flottante dans cette région. Cela engendre une accélération de l'écoulement des glaciers émissaires et une plus grande décharge de glace vers l'océan. Ainsi, la contribution récente à la hausse du niveau de la mer de l'Antarctique de l'Ouest s'élève à environ 0.2 mm par an entre 1992 et 2011, c'est-à-dire près du tiers de la contribution des calottes glaciaires (Antarctique et Groenland). On estime que cette contribution va continuer à augmenter dans le futur proche.<p>Dans cette thèse, nous démontrons clairement que la modélisation de la migration de la ligne d'ancrage (frontière entre glaces posée et flottante) dépend de l'approche numérique et de l'approximation physique du modèle cryosphérique utilisé. Les modèles cryosphériques qui prescrivent le flux glaciaire à la ligne d'ancrage et qui utilisent un niveau de physique et une approche numérique appropriés convergent vers la même position d'équilibre de la ligne d'ancrage quelle que soit la taille de maille utilisée. Cependant, le comportement transitoire de ces modèles est moins précis que d'autres modèles et mène à une surestimation du flux à la ligne d'ancrage. Dès lors, ces modèles doivent être utilisés avec précaution sur de courtes périodes de temps. De plus, la non inclusion des contraintes verticales de cisaillement dans ces modèles augmente la viscosité effective et donne des positions d'équilibre de la ligne d'ancrage plus en aval en comparaison avec les modèles « full-Stokes ».<p>La seconde découverte majeure de cette thèse est le contrôle important exercé par la géométrie (largeur du glacier et topographie du lit rocheux) sur Thwaites Glacier, l'un des glaciers émissaires les plus rapides en Antarctique de l'Ouest. Un modèle « Shallow-Shelf Approximation » (SSA) résolvant les différences finies le long d'une ligne d'écoulement est appliqué au glacier et montre que la convergence de l'écoulement glaciaire (au travers de la paramétrisation de la largeur) ralentit le retrait de la ligne d'ancrage comparé aux simulations où la largeur est constante. Une nouvelle paramétrisation de l'effet arc-boutant est testée sur le glacier et permet de mieux comprendre cet effet. Finalement, la version en trois dimensions du modèle ci-dessus est appliquée à Thwaites Glacier et met en évidence le contrôle important des variations latérales de l'altitude du lit rocheux sur la migration de la ligne d'ancrage. / Doctorat en Sciences / info:eu-repo/semantics/nonPublished

Glaciers respond to local climate changes making them important indicators of regional climate change. The North Cascades region of Washington is the most glaciated region in the lower-48 states with approximately 25% of all glaciers and 40% of the total ice-covered area. While there are many on-going investigations of specific glaciers, little research has addressed the entire glacier cover of the region. A reference inventory of glaciers was derived from a comparison of two different inventories dating to about 1958. The different inventories agree within 93% of total number of glaciers and 94% of total ice-covered area. To quantify glacier change over the past century aerial photographs, topographic maps, and geologic maps were used. In ~1900 total area was about 533.89 ± 22.77 km2 and by 2009 the area was reduced by -56% ± 3% to 236.20 ± 12.60 km2. Most of that change occurred in the first half of the 20th century, between 1900 and 1958, -245.59 ± 25.97 km2 (-46% ± 5%) was lost, followed by a period of stability/growth in mid-century (-1% ± 3% from 1958-1990) then decline since the 1990s (-9% ± 3% from 1990-2009). The century-scale loss is associated with increasing regional temperatures warming in winter and summer; precipitation shows no trend. On a decadal time scale winter precipitation and winter and summer temperatures are important factors correlated with area loss. Topographically, smaller glaciers at lower elevations with steeper slopes and higher mean insolation exhibited greater loss than higher, gentler more shaded glaciers.

In this dissertation, we have examined mineral assemblages and geomorphologic features in the Sisyphi Planum region of Mars, as well as examined the mineral assemblage of palagonite in Iceland. Chapter 2 is focused on the mineral assemblages detected on possible glaciovolcanic edifices in the Sisyphi Planum region of Mars. Minerals were identified utilizing visible/near-infrared orbital spectra from the Compact Reconnaissance Imaging Spectrometer for Mars (CRISM). Analysis of eleven CRISM images located on the volcanic edifices revealed three distinct spectral classes in the region which are interpreted to be: gypsum-dominated, smectite-zeolite- iron oxide-dominated (possibly palagonite), and polyhydrated sulfate-dominated material. The possible palagonite detections on the volcanic edifices, the geomorphology of the region, and the analogous terrestrial mineralogy of subglacial eruptions strongly suggests the formation of these minerals during subglacial eruptions or associated hydrothermal systems. This implies that thick water ice sheets were present in this region in the late Noachian or early Hesperian, and that the subglacial hydrothermal systems could have supported habitable environments with excellent biosignature preservation potential. Chapter 3 is focused on evaluating the variability of the composition and crystallinity of palagonite on Earth in order to inform efforts to identify it on Mars. We hypothesized that variability in palagonite composition and crystallinity could occur due to differences in environmental conditions during formation. Palagonite samples were collected in Iceland at subglacial volcanic sites around Reykjavík in the Western Volcanic Zone, on the southern coast in the Eastern Volcanic Zone, and from the Herðubreið tuya and Askja volcano in the Northern Volcanic Zone. Visible/near-infrared reflectance spectroscopy, thermal-infrared emission spectroscopy, and quantitative XRD were used to assess the bulk mineralogy, crystallinity, and clay composition of all samples. Results show the sampled palagonites contain partially devitrified glass, unaltered glass, and secondary minerals including clay minerals, poorly crystalline ferric oxides, and zeolites. However, one sample (SCoast01) shows a vastly different mineral assemblage in all sample techniques, including well-crystalline Fe/Mg-clays as opposed to the poorly-crystalline Al-clays observed in our other samples. Based on previous studies of subaqueous palagonites and the location this sample was collected from, we hypothesize that the SCoast01 sample was formed in a submarine environment rather than subglacial. This suggests that it may be possible to differentiate submarine vs. subglacial palagonite on Earth based on composition and from remote sensing observations on Mars. Chapter 4 is a geomorphologic study of the Sisyphi Planum region of Mars where we identified and classified the tops of the Sisyphi Montes as well as geomorphologically mapped the Sisyphi Planum region. Here, we address an overarching question: What is the relationship between the Sisyphi Montes and the ice in this region? To do this, we identified 106 edifices in the region and classified them into five categories: 1) flat topped, 2) rounded tops, 3) sharp peaks, 4) cratered peaks, and 5) height less than 300 meters – a “catch-all” category for all features below the specified height, which exhibit less distinctive morphologies in MOLA topography. While many of the edifices could be sub-glacial in origin, we find that the only morphologic class that exhibits uniquely subglacial morphologies are the flat-topped edifices. These edifices are similar to terrestrial tuyas, which form when a subglacial volcano breaches an ice sheet and erupts a plateau of sub-aerial lavas. Based on the geomorphologic map and topographic data, we have shown that flat-topped edifices are all located outside of regions that we map as the Mantled Unit, which we infer to be related to the Dorsa Argentina Formation. The combination of the flat topped edifices and their location outside of the mapped ice-related regions strongly suggests that the ice in the region was once more extensive than what is currently observed. While this has been proposed in the past, it has not been documented how far the ice sheet could have extended. Here we show that the ice must have extended to at least as far as the flat topped edifices in the region. The combination of these chapters using both mineralogy and morphology suggest that the Sisyphi Planum region of Mars was subglacial in origin. <br><br>

The internal structure of ice sheets reflects the history of flow and deformation experienced by the ice mass. Flow and deformation are controlled by processes occurring within the ice mass and at its boundaries, including surface accumulation or ablation, ice rheology, basal topography, basal sliding, and basal melting or freezing. The internal structure and basal environment of ice sheets is studied with ice-penetrating radar. Recently, radar observations in Greenland and Antarctica have imaged large englacial structures rising from near the bed that deform the overlying stratigraphy into anticlines, synclines, and overturned folds. The mechanisms that may produce these structures include basal freeze-on, travelling slippery patches at the ice base, and rheological contrasts within the ice column. In this thesis, I explore the setting and mechanisms that produce large basal stratigraphic structures inside ice sheets. First, I use radar data to map subglacial hydrologic networks that deliver meltwater uphill towards freeze-on structures in East Antarctica. Next, I use a thermomechanical flowline model to demonstrate that trains of alternating slippery and sticky patches can form underneath ice sheets and travel downstream over time. The disturbances to the ice flow field produced by these travelling patches produce stratigraphic folds resembling the observations. I then examine the overturned folds produced by a single travelling sticky patch using a kinematic flowline model. This model is used to interpret stratigraphic measurements in terms of the dynamic properties of basal slip. Finally, I use a simple local one-dimensional model to estimate the thickness of basal freeze-on that can be produced based on the supply of available meltwater, the thermal boundary conditions, ice sheet geometry, and the ice flow regime.

The dynamic response of a faster-flowing Greenland Ice Sheet to climate change is modulated by feedbacks between ice flow and surface meltwater delivery to the basal environment. While supraglacial melt processes have been thoroughly examined and are well constrained, the response of the englacial and subglacial environment to these seasonal perturbations still represent the least-studied, understood, and parameterised processes of glacier dynamics due to a paucity of direct observation. To better understand these processes in the wake of a changing climate, novel in-situ geophysical experiments were undertaken on Store Glacier in west Greenland to quantify rates of englacial deformation and basal melting. The records produced from these experiments yield new insights into the thermodynamic setting of a major outlet glacier, and the physical mechanisms underlying and resulting from fast glacier motion. The deployment of autonomous phase-sensitive radio-echo sounders (ApRES) $\SI{30}{\kilo\metre}$ from the calving terminus of Store Glacier between 2014 and 2016 revealed high rates of both englacial deformation and basal melting, driven primarily by the dynamic response of the basal hydrological system to seasonal surface meltwater influxes. Thermodynamic modelling of this process revealed a convergence of large-scale basal hydrological pathways that aggregated large amounts of water towards the field site. The warm, turbulent water routed from surface melt contained and dissipated enough energy at the ice-bed interface to explain the observed high melt rates. Simultaneously, changes in the local strain field, involving seasonal variations in the morphology of internal layers, were found to be the result of far-field perturbations in downstream ice flow which propagated tens of kilometres upglacier through longitudinal stress coupling. When observed in multiple dimensions, the layer structure revealed complex internal reflection geometries, demonstrating ApRES as not just a monitor of depth changes in ice thickness, but also as an imaging instrument capable of characterising the subsurface environment within and beneath ice sheets. Altogether, the observations and analyses comprising this thesis provide new and significant insight and understanding into the structural, thermal, and mechanical processes tied to Store Glacier and its fast, complex, and dynamic ice flow.

To increase our understanding of how ice sheets and glaciers interact with the climate system, numerical models have become an indispensable tool. However, the complexity of these systems and the natural limitation in computational power is reflected in the simplifications of the represented processes and the spatial and temporal resolution of the models. Whether the effect of these limitations is acceptable or not, can be assessed by theoretical considerations and by validating the output of the models against real world data. Equally important is to verify if the numerical implementation and computational method accurately represent the mathematical description of the processes intended to be simulated. This thesis concerns a set of numerical models used in the field of glaciology, how these are applied and how they relate to other study areas in the same field. The dynamical flow of glaciers, which can be described by a set of non-linear partial differential equations called the Full Stokes equations, is simulated using the finite element method. To reduce the computational cost of the method significantly, it is common to lower the order of the used elements. This results in a loss of stability of the method, but can be remedied by the use of stabilization methods. By numerically studying different stabilization methods and evaluating their suitability, this work contributes to constraining the values of stabilization parameters to be used in ice sheet simulations. Erroneous choices of parameters can lead to oscillations of surface velocities, which affects the long term behavior of the free-surface ice and as a result can have a negative impact on the accuracy of the simulated mass balance of ice sheets. The amount of basal sliding is an important component that affects the overall dynamics of the ice. A part of this thesis considers different implementations of the basal impenetrability condition that accompanies basal sliding, and shows that methods used in literature can lead to a difference in velocity of 1% to 5% between the considered methods. The subglacial hydrological system directly influences the glacier's ability to slide and therefore affects the velocity distribution of the ice. The topology and dominant mode of the hydrological system on the ice sheet scale is, however, ill constrained. A third contribution of this thesis is, using the theory of R-channels to implement a simple numerical model of subglacial water flow, to show the sensitivity of subglacial channels to transient processes and that this limits their possible extent. This insight adds to a cross-disciplinary discussion between the different sub-fields of theoretical, field and paleo-glaciology regarding the characteristics of ice sheet subglacial hydrological systems. In the study, we conclude by emphasizing areas of importance where the sub-fields have yet to unify: the spatial extent of channelized subglacial drainage, to what degree specific processes are connected to geomorphic activity and the differences in spatial and temporal scales. As a whole, the thesis emphasizes the importance of verification of numerical models but also acknowledges the natural limitations of these to represent complex systems. Focusing on keeping numerical ice sheet and glacier models as transparent as possible will benefit end users and facilitate accurate interpretations of the numerical output so it confidently can be used for scientific purposes. / <p>At the time of the doctoral defense, the following papers were unpublished and had a status as follows: Paper 3: Manuscript. Paper 4: Manuscript.</p> / Greenland Analogue Project

Ice cores are a valuable component with regards to paleoclimate reconstruction due to the ability to use stable water isotopic concentrations in ice as a proxy for paleo-temperature records. It is therefore important to understand the processes and conditions under which isotopic concentrations can be altered after ice has formed. Historically, sublimation has been considered to only have a trivial impact on the isotopic record in glacial ice due to the low diffusivity of solid ice (~10-15 m2 s-1). Recent publications have shown that diffusion of impurities through ice can occur at much faster rates than the diffusivity of solid ice would imply, and have proposed that networks of unfrozen liquid (premelt) between ice grains may expedite the diffusion process. However, the application of this mode of diffusion to isotopic concentrations in ice under non-equilibrium conditions has been largely unexplored. Here I model changes in isotopic concentrations in ice using a two-dimensional diffusion mechanism, which incorporates premelt, coupled with a sublimation flux at the surface. Model results show an increase in δD at the ice surface and in near-surface ice. Concentrations exponentially decrease from the surface value to the initial concentration at depth. These results are consistent with recent experimental results.