Ces travaux de thèse se focalisent sur la modélisation numérique de la réponse thermique des glaciers froids aux variations climatiques. La prise en compte du forçage climatique par la modification des conditions de surface des glaciers est l'un des points clefs de cette étude. On montre l'importance des processus de fusion/regel en zone d'accumulation froide ainsi que de l'épaisseur de la couverture de neige et névé autour de la ligne d'équilibre. Les modèles développés dans cette thèse s'appliquent à des cas concrets sur trois glaciers du massif du Mont Blanc (Alpes françaises). Ces études permettent d'analyser des situations à risque et de reconstituer le climat passé à haute altitude. Le jeu de données très complet acquis sur le site du Col du Dôme (4250 m) comprend des mesures profondes de températures et de densités ainsi que des mesures d'accumulation et de vitesses de surface. Il a permis de valider et développer un modèle à trois dimensions de régime thermique couplé à un modèle d'écoulement et forcé par des données climatiques. Ce modèle permet d'étudier en détail la réponse thermique des zones d'accumulation froides aux changements climatiques. Il peut aussi être utilisé pour la prédiction du réchauffement futur des glaciers suspendus présentant un risque comme le glacier de Taconnaz dont une analyse préliminaire est effectuée ici. Le réchauffement basal jusqu'au point de fusion pourrait conduire à la déstabilisation de ce type de glacier. Le site du glacier de Tête Rousse ( environ 3200 m) a permis d'établir clairement l'influence de la couverture de neige et de névé sur le régime thermique des glaciers autour de la ligne d'équilibre. Les résultats obtenus à partir du couplage d'un modèle de neige (CROCUS) avec un modèle de régime thermique sont en excellent accord avec les observations de températures réalisées dans le glacier. On explique ainsi les causes de la structure thermique particulière du glacier de Tête Rousse responsable de la formation d'une poche présentant un risque pour la commune à l'aval. Dans le futur, les simulations indiquent que le glacier va continuer de se refroidir, limitant ainsi le risque de formation de poche d'eau. Dans le dernier chapitre, nous avons développé une nouvelle méthode inverse pour reconstituer les températures du passé à partir de profils de températures observées dans la glace. Elle est basée sur l'inversion simultanée de plusieurs profils de températures provenant de différents sites de forage ayant subi le même forçage climatique. La méthode permet de s'affranchir de l'impact des processus de fonte/regel en surface pour reconstruire les variations passées des températures de l'air. Les résultats obtenus sur le Col du Dôme montrent que le réchauffement climatique sur le site à 4250 m est similaire à celui observé à basse altitude. Ces résultats tendent à montrer que le réchauffement climatique n'est pas amplifié avec l'altitude.

Published also as Thesis--Tübingen. / Bibliography: p. [150]-160.

Ces travaux de thèse se focalisent sur la modélisation numérique de la réponse thermique des glaciers froids aux variations climatiques. La prise en compte du forçage climatique par la modification des conditions de surface des glaciers est l'un des points clefs de cette étude. On montre l'importance des processus de fusion/regel en zone d'accumulation froide ainsi que de l'épaisseur de la couverture de neige et névé autour de la ligne d'équilibre. Les modèles développés dans cette thèse s'appliquent à des cas concrets sur trois glaciers du massif du Mont Blanc (Alpes françaises). Ces études permettent d'analyser des situations à risque et de reconstituer le climat passé à haute altitude. Le jeu de données très complet acquis sur le site du Col du Dôme (4250 m) comprend des mesures profondes de températures et de densités ainsi que des mesures d'accumulation et de vitesses de surface. Il a permis de valider et développer un modèle à trois dimensions de régime thermique couplé à un modèle d'écoulement et forcé par des données climatiques. Ce modèle permet d'étudier en détail la réponse thermique des zones d'accumulation froides aux changements climatiques. Il peut aussi être utilisé pour la prédiction du réchauffement futur des glaciers suspendus présentant un risque comme le glacier de Taconnaz dont une analyse préliminaire est effectuée ici. Le réchauffement basal jusqu'au point de fusion pourrait conduire à la déstabilisation de ce type de glacier. Le site du glacier de Tête Rousse ( environ 3200 m) a permis d'établir clairement l'influence de la couverture de neige et de névé sur le régime thermique des glaciers autour de la ligne d'équilibre. Les résultats obtenus à partir du couplage d'un modèle de neige (CROCUS) avec un modèle de régime thermique sont en excellent accord avec les observations de températures réalisées dans le glacier. On explique ainsi les causes de la structure thermique particulière du glacier de Tête Rousse responsable de la formation d'une poche présentant un risque pour la commune à l'aval. Dans le futur, les simulations indiquent que le glacier va continuer de se refroidir, limitant ainsi le risque de formation de poche d'eau. Dans le dernier chapitre, nous avons développé une nouvelle méthode inverse pour reconstituer les températures du passé à partir de profils de températures observées dans la glace. Elle est basée sur l'inversion simultanée de plusieurs profils de températures provenant de différents sites de forage ayant subi le même forçage climatique. La méthode permet de s'affranchir de l'impact des processus de fonte/regel en surface pour reconstruire les variations passées des températures de l'air. Les résultats obtenus sur le Col du Dôme montrent que le réchauffement climatique sur le site à 4250 m est similaire à celui observé à basse altitude. Ces résultats tendent à montrer que le réchauffement climatique n'est pas amplifié avec l'altitude. / This thesis focuses on the numerical modeling of the thermal response of glaciers to climate changes. The change in surface conditions related to climate forcing is one of the key point of this study. We point out the crucial role of the meltwater refreezing in cold accumulation zone and the thickness of the snow and firn cover around the equilibrium line. The models are developed from the study of specific cases of three glaciers in Mont Blanc range ( French Alps). These studies allow us to analyze the risk related to hanging glaciers and to reconstruct the past temperaures at high altitude from observede nglacial temperatures. The very extensive dataset obtained at the site of the Col du Dôme (4250 m) includes measurements of deep temperatures and density profiles as well as measurements of snow accumulation and surface velocities. These observations allow us to develop and validate a three-dimensional thermal regime model coupled with a flow model and forced by meteorological data. This model is used to perform a thorough study of the thermal response of cold accumulation zone to climate change. It can also be used for the prediction of future glacial hazard related to hanging glaciers warming. A preliminary analysis has been carried out on the hanging glacier of Taconnaz. The warming could have a major impact on the stability of this kind of glaciers frozen to their beds if the melting point is reached. The study performed on Tête Rousse glacier (about 3200 m) shows clearly the influence of the snowpack and firn cover on the thermal regime of glaciers around the equilibrium line. The results obtained from coupling of a snow model (CROCUS) with a model of thermal regime reveal a very good agreement with the observed ice temperatures. This study provides insights into the thermal processes responsible for water storage inside a small almost static glacier, which can lead to catastrophic outburst floods. In the future, according to atmospheric temperature increase scenarios for the coming century, most of the glacier will become cold which will reduce the risk of water filled cavity formation. Finally, a new inverse method has been developed to reconstruct temperatures in the past from observed englacial temperatures. The method is based on the simultaneous inversion of several temperature profiles coming from different drilling sites with the same climate forcing. The method overcomes the impact of refreezing meltwater to reconstruct past air temperatures variations. This is similar to the observed regional low altitude trend. The results obtained on the Col du Dome shows that climate warming on the site at 4250 m is similar to the observed regional low altitude trend in the northwestern Alps, suggesting that air temperature trends are not altitude dependent.

La région de l'archipel canadien, située en Arctique, connaît actuellement d'importants changements climatiques, se traduisant notamment par une augmentation des températures, une réduction de l'étendue de la banquise marine et du couvert nival terrestre ou encore une perte de masse significative des calottes glaciaires disséminées sur les îles de l'archipel. Parmi ces calottes glaciaires, la calotte Barnes, située en Terre de Baffin, ne fait pas exception comme le montrent les observations satellitaires qui témoignent d'une importante perte de masse ainsi que d'une régression de ses marges, sur les dernières décennies. Bien que les calottes glaciaires de l'archipel canadien ne représentent que quelques dizaines de centimètres d'élévation potentielle du niveau des mers, leur perte de masse est une composante non négligeable de l'augmentation actuelle du niveau des mers. Les projections climatiques laissent à penser que cette contribution pourrait rester significative dans les décennies à venir. Cependant, afin d'estimer les évolutions futures de ces calottes glaciaires et leur impact sur le climat ou le niveau des mers, ils est nécessaire de caractériser les processus physiques tels que les modifications du bilan de masse de surface. Cette connaissance est actuellement très limitée du fait notamment du sous-échantillonnage des régions arctiques en terme de stations météorologiques permanentes. Une autre particularité de certaines calottes de l'archipel canadien, et de la calotte Barnes en particulier, est de présenter un processus d'accumulation de type glace surimposée, ce phénomène étant à prendre en compte dans l'étude des processus de surface. Pour palier au manque de données, l'approche retenue a été d'utiliser des données de télédétection, qui offrent l'avantage d'une couverture spatiale globale ainsi qu'une bonne répétitivité temporelle. En particulier les données acquises dans le domaine des micro-ondes passives est d'un grand intérêt pour l'étude de surfaces enneigées. En complément de ces données, la modélisation du manteau neigeux, tant d'un point de vue des processus physiques que de l'émission électromagnétique permet d'avoir accès à une compréhension fine des processus de surface tels que l'accumulation de la neige, la fonte, les transferts d'énergie et de matière à la surface, etc. Ces différents termes sont regroupés sous la notion de bilan de masse de surface. L'ensemble du travail présenté dans ce manuscrit a donc consisté à développer des outils permettant d'améliorer la connaissance des processus de surface des calottes glaciaires du type de celles que l'on rencontre dans l'archipel canadien, l'ensemble du développement méthodologique ayant été réalisé sur la calotte Barnes à l'aide du schéma de surface SURFEX-CROCUS pour la modélisation physique et du modèle DMRT-ML pour la partie électromagnétique. Les résultats ont tout d'abord permis de mettre en évidence une augmentation significative de la durée de fonte de surface sur la calotte Barnes (augmentation de plus de 30% sur la période 1979-2010), mais aussi sur la calotte Penny, elle aussi située en Terre de Baffin et qui présente la même tendance (augmentation de l'ordre de 50% sur la même période). Ensuite, l'application d'une chaîne de modélisation physique contrainte par diverses données de télédétection a permis de modéliser de manière réaliste le bilan de masse de surface de la dernière décennie, qui est de +6,8 cm/an en moyenne sur la zone sommitale de la calotte, qui est une zone d'accumulation. Enfin, des tests de sensibilité climatique sur ce bilan de masse ont permis de mettre en évidence un seuil à partir duquel cette calotte voit disparaître sa zone d'accumulation. Les modélisations effectuées suggèrent que ce seuil a de fortes chances d'être atteint très prochainement, pour une augmentation de température moyenne inférieure à 1°C, ce qui aurait pour conséquence une accélération de la perte de masse de la calotte. / Significant climate change is curently monitored in the Arctic, and especially in the region of the canadian arctic archipellago. This climate warming leads to recession of sea-ice extent and seasonnal snow cover, and also to large mass loss of the archipellago's ice caps. One of the most southern ice cap, the Barnes Ice Cap, located on the Baffin Island, is no exception to significant mass loss and margins recession as satellite observations exhibited over the last decades. Despite the relative low sea level potential of the small ice caps located in the canadian arctic achipellago in regards to major ice sheets, Antarctica and Greenland, their contribution to the current sea level rise is significant. Climate projections show that this contribution could accelerate significant over the next decades. However, to estimate the future evolution of these ice caps and their impact on climate or sea level rise, a better characterisation of the surface processes such as the evolution of the surface mass balance is needed. This knowledge is currently very limited, mainly due to the sparse covering of automatic weather stations or in-situ measurements over the Arctic. Furthermore, several ice caps, among with the Barnes Ice Cap, present a superimposed ice accumulation area which particularities have to be taken into account in the surface processes studies. Given the lack of in-situ data, the approach choosen in this work is to use remote sensing data, that have the advantage to offer a good spatial and temporal coverage. In particular, passive microwave data are very suitable for snowy surfaces studies. To complement these data, physical and electromagnetic snowpack modeling provide a fine characterisation of surface processes such as snow accumulation. The whole work presented in this manuscript thus consisted in developping specific tools to improve the understanding of surface processes of small arctic ice caps. This methodological development was performed and applied on the Barnes Ice Cap using the surface scheme SURFEX-CROCUS and the electromagnetic model DMRT-ML. First results highlight a significant increase in surface melt duration over the past 3 decades on the Barnes Ice Cap (increase of more than 30% over 1979-2010 period). A similar trend is also monitored over the Penny Ice Cap, located in the south part of the Baffin Island (increase of more than 50% over the same period). Then, the surface mass balance over the last decade was modeled by using a physical based modeling chain constrained by remote sensing data. The results give a mean net accumulation of +6,8 cm/an on the summit area of the ice cap. Finaly, sensitivity tests, performed to investigate the climatic sensitivity of the surface mass balance, highlight a threshold effect that may lead to a complete disapearence of the accumulation area of the Barnes Ice Cap. With a temperature increase less than 1°C, modeling results suggest it is likely that the threshold will be reached rapidly leading to an increase in mass loss from the ice cap.

La région de Charlevoix, dans le sud-est du Québec, est caractérisée par un paysage spectaculaire formé par la jonction des Hautes-Laurentides, de l’astroblème de Charlevoix et de l'estuaire du Saint-Laurent. Cette étude, basée sur deux campagnes de terrain et sur l'interprétation de photos aériennes numériques de haute résolution 3D et de données LiDAR, documente un géosystème glaciaire dominé par le développement de lacs proglaciaires lors de la déglaciation en raison d'une topographie très irrégulière et de la présence de grands systèmes morainiques. Les résultats fournissent une chronologie relative pour le développement et le drainage d’une série de lacs proglaciaires pendant chaque étape de la déglaciation. Dans la région, les Hautes-Laurentides sont caractérisées par un haut plateau (700 - 1174 m) et des vallées profondément encaissées (Malbaie, Sainte-Anne, à Mars). L’astroblème de Charlevoix est caractérisé par un graben circulaire (vallées du Gouffre et Malbaie) et une remontée centrale (mont des Éboulements et collines de Saint-Hilarion). Quatre systèmes morainiques sont observés dans la région: les systèmes Rochette et Brûlée-Sainte-Anne, probablement développés pendant la Dryas ancien (~ 14 cal. ka BP.), et les systèmes de Mars-Batiscan et de Saint-Narcisse, datant du Dryas récent (12,9 - 11,7 cal. ka BP). La présence de deltas, de rythmites, de dépôts littoraux, de terrasses lacustres, de brèches de drainage et de dépôts de drainage a été utilisée pour identifier 91 lacs proglaciaires. Cinq systèmes de lacs proglaciaires ont été identifiés sur la base de leur emplacement, leur chemin de drainage, leur type de barrages, leur complexe morainique associé et leur stratigraphie: 1- Système des lacs se drainant la vallée du Saint-Laurent : des lacs étant contenus par le courant glaciaire du Saint-Laurent ou par les moraines de Rochette; 2- Système des vallées latérales de l’astroblème : des lacs se drainant dans le graben de l’astroblème, contenus par les langues glaciaires du Gouffre et Malbaie; 3- Système du secteur de Sainte-Aimée-des-lacs : des lacs situés dans le graben de l’astroblème, contenus par des crêtes morainiques associées au complexe de Saint-Narcisse; 4- Système du Secteur Sainte-Anne : des lacs se drainant par les vallées des rivières Sainte-Anne et Montmorency, contenus par les moraines du complexe de Brûlé-Sainte-Anne ou par des glaciers de vallées; 5- Système des Laurentides Nord : des lacs se drainant vers le nord (Saguenay, rivière Noire), contenus d’abord par l’Inlandsis laurentidien et ensuite par les moraines des complexes de Saint-Narcisse et de Mars-Batiscan. / The Charlevoix region, in southeastern Québec, is characterized by a dramatic landscape formed by the junction of the Laurentian Highlands, the Charlevoix Astrobleme and the St. Lawrence Estuary. This study, based on two extensive field campaigns and the interpretation of high-resolution 3D digital air photos and LiDAR data, documents a glacial landsystem dominated by the development of proglacial lakes during deglaciation due to a very irregular topography and the presence of major moraine systems. Our results provide a relative chronology for the development and drainage of series of proglacial lakes formed during each successive stage of deglaciation. In the region, the Laurentian Highlands are characterized by a high plateau (700-1174 m) and deeply incised river valleys (Malbaie, Saint-Anne, à Mars). The Charlevoix Astrobleme is characterized by a circular graben (Malbaie and du Gouffre river valleys) and a central uplift (Mount of Éboulements and Saint-Hilarion hills). Four moraine systems are observed in the region: the Rochette and Brûlée-Sainte-Anne systems, probably developed during the Older Dryas (~ 14 cal. ka BP), and the Saint-Narcisse and Mars-Batiscan systems dating from the Younger Dryas (12.9 to 11.7 cal. ka BP). Deltas, rythmites, coastal deposits, lacustrine terraces, breaches and drainage deposits were used to identify 91 proglacial lakes. Five types of proglacial lakes were identified based on their location, drainage path, type of dams, associated moraine complex and stratigraphy: 1- Lakes draining towards the St. Lawrence valley : constrained by glacial ice in the St. Lawrence Estuary or by the Rochette moraine, and draining later into the Goldthwait Sea; 2- Astrobleme side valley lakes : constrained by du Gouffre and Malbaie glacial tongues and draining into the graben of the astrobleme; 3- Sainte-Aimée-des-Lacs lakes : constrained by moraine ridges associated with Saint-Narcisse moraine systems; 4- Sainte-Anne Sector lakes : constrained by Brûlé-Sainte-Anne moraines or by valley glaciers lakes and draining through the Sainte-Anne and Montmorency river valleys; 5- Northern Laurentian lakes: constrained first by the Laurentide Ice Sheet and then by the moraines of the Saint-Narcisse and Mars-Batiscan moraine systems and draining northward towards Lake St-Jean and Saguenay River.

Résumé : Le contexte récent d’accélération de la fonte des glaciers et calottes glaciaires (GCG) de l’archipel arctique canadien, jumelé aux difficultés de suivi des GCG de cette région, rendent essentiels le développement et l’utilisation de nouvelles approches innovatrices de suivi. Le potentiel de la stéréo radargrammétrie (SRG) RADARSAT-2 est ici caractérisé pour l’extraction d’élévations et le calcul de changements d’élévation et de bilans de masse (historiques et récents) sur les calottes glaciaires Barnes et Penny (Nunavut, Canada). Par la méthode semi-automatisée de recherche de corrélation à partir de couples stéréoscopiques RADARSAT-2 de 2013 (mode wide ultra-fin; résolution spatiale de 3 m; taille d’image de 50 km x 50 km), une précision verticale de ~7 m (LE68) est mesurée sur la terre ferme, et cette valeur de précision est possiblement légèrement supérieure sur la calotte Barnes, étant donné la variabilité de profondeur de pénétration. Par captage 3D, une précision altimétrique de ~3-4 m (LE68) est mesurée par différents photo-interprètes à partir de couples RADARSAT de 2012 en zone d’ablation de la calotte Penny. Sur la calotte Barnes, les changements d’élévation mesurés par rapport aux premiers modèles numériques de terrain disponibles permettent de mesurer un bilan de masse spécifique historique (1960-2013) de -0,49 ± 0,20 m w.e./année, pour un bilan de masse total de -2,9 Gt/année. Entre 2005 et 2013, le bilan de masse spécifique de cette calotte augmente significativement à -1,20 ± 0,86 m w.e./année, pour un bilan de masse total de -7 Gt/année. En zone d’ablation de la calotte Penny, un changement d’élévation annuel moyen de -0,59 m/année est mesuré entre 1958 et 2012. Parallèlement, plusieurs aspects méthodologiques et techniques sont discutés et analysés. Des profondeurs de pénétration nulles (bande C) sont mesurées à partir des images acquises sur la calotte Barnes à la toute fin de la saison d’ablation (fin septembre/début octobre), alors que cette profondeur augmente à ~2,5-3 m pour des images acquises à la fin octobre/début novembre (période de gel). Nos résultats suggèrent aussi que le modèle de fonction rationnelle, lorsqu’utilisé avec des images RADARSAT-2 en mode wide ultra-fin, permet d’obtenir des précisions plus constantes que le modèle hybride de Toutin. De par son indépendance des conditions météorologiques, son utilisation possible sans point de contrôle et sa simplicité de traitement, la SRG RADARSAT-2 s’avère donc être une excellente alternative aux technologies actuelles pour le suivi de GCG situés dans des régions affectées par des contraintes opérationnelles importantes. / Abstract : Given the recent melt acceleration of the Canadian arctic archipelago’s ice caps and the monitoring difficulties of this remote region, the development of new innovative monitoring tools has become essential. Here, the potential of the RADARSAT-2 stereo radargrammetry (SRG) is characterized for elevations extraction, as well as for elevation changes/mass balances calculations (historical and recent) on Barnes and Penny ice caps (Nunavut, Canada). Using the semi-automatic approach of correlation search from RADARSAT-2 stereoscopic couples of 2013 (wide ultra-fine mode; spatial resolution of 3 m; coverage of 50 km x 50 km), a vertical precision of ~7 m (LE68) is measured on ice-free terrain and this precision is possibly slighty worse on the ice cap because of the penetration depth’s variability. On the other hand, the 3D vision extraction approach reveals an altimetric precision of ~3-4 m (LE68) on the ablation area of the Penny Ice Cap. On the Barnes Ice Cap, elevation changes calculated relative to the oldest digital elevation models available allows to calculate an historical specific mass balance (1960-2013) of -0,49 ± 0,20 m w.e./year, resulting in a total annual mass balance of -2,9 Gt/year. Between 2005 and 2013, the specific mass balance of this ice cap increases to -1,20 ± 0,86 m w.e./year, which equals to a total annual mass balance f -7 Gt/year. On Penny Ice Cap’s ablation area, an average elevation change of -0,59 m/year is measured between 1958 and 2012. As also suggested in the literature, the recent melt acceleration is highly linked to warmer summer temperatures. Methodological and technical aspects are also presented and analyzed. No penetration depth (C band) is perceived on elevations derived from late ablation season images (late September/beginning of October), while a penetration of ~2,5-3 m is measured from images acquired in late October/beginning of November (freeze period). Our results also suggest the superiority and better consistency of the rational function model for geometrical correction of wide ultra-fine mode RADARSAT-2 images, compared to the hybrid Toutin’s model. Because of its all-weather functionality, its possible use without any ground control point and the simplicity and facility of its treatment, the RADARSAT-2 SRG represents a really good technology for glacier monitoring in regions affected by serious operational constraints.

Les glaciers tropicaux sont des indicateurs climatiques très sensibles. Pour relier les fluctuations de leur bilan de masse aux changements climatiques, nous avons réalisé un bilan d'énergie à la surface du Glacier Zongo (16°S), en Bolivie, entre mars 96 et août 98. Cette thèse présente, d'une part les méthodes employées pour obtenir le cycle annuel du bilan d'énergie ponctuel, et d'autre part l'intérêt de cette étude pour la compréhension du fonctionnement du glacier et sa relation avec la climatologie locale. Dans la première partie, nous nous attachons à décrire les conditions climatiques régnant sur ce glacier, à partir des enregistrements des stations météorologiques installées sur le Glacier Zongo. Les glaciers tropicaux sont soumis à des températures et des radiations solaires incidentes peu variables au cours de l'année, alors que le débit du torrent émissaire est très saisonnier. Pour comprendre quels sont les flux d'énergie responsables de la saisonnalité de la fusion en surface du glacier, nous avons dressé le bilan d'énergie en un point du glacier, dans la zone d'ablation, à 5150 m d'altitude. Ce bilan d'énergie traduit un métabolisme propre aux glaciers tropicaux. L'albédo est le paramètre essentiel qui contrôle la quantité d'énergie disponible pour la fusion en surface du glacier. La particularité fondamentale de ces glaciers réside en la très forte consommation d'énergie sous forme de chaleur latente (ce qui se traduit par une sublimation élevée). En outre, le flux de chaleur latente présente une saisonnalité marquée à l'origine du régime hydrologique du torrent qui s'échappe du glacier. En effet, en saison sèche, la sublimation est maximale et il reste peu d'énergie pour la fusion d'où le faible débit du torrent émissaire. La situation opposée se produit en saison humide. À l'échelle de l'ensemble de la période étudiée, le bilan de masse nettement déficitaire de l'année El Nino 1997-98 s'explique surtout par un déficit de précipitations, responsable d'une chute de l'albédo moyen annuel. Le manteau neigeux réduit lors de la saison d'accumulation 1997-98 a disparu plus vite que celui de l'année précédente, laissant affleurer une surface de glace sale, d'albédo faible, pendant une période plus longue. De plus, la chaleur sensible plus élevée et la sublimation réduite ont aussi contribué à augmenter la fusion superficielle, diminuant ainsi le bilan de masse du glacier.

Sous l'effet de charges en surface, la surface terrestre s'affaisse pour retrouver un état d'équilibre isostatique. Inversement, quand ces charges viennent à diparaître, la surface remonte pour retrouver sa configuration initiale. Dans le cas d'une calotte glaciaire, la taille mais aussi la vitesse à laquelle peut évoluer la masse de glace est à l'origine d'une réponse hautement spécifique et particulièrement de comportements transitoires associés. Après avoir introduit la notion d'isostasie, nous décrivons quelles peuvent être les conséquences à la fois de la taille et de la vitesse d'évolution de la charge sur la réponse terrestre. Cela permet de clarifier la notion de "contexte glaciaire" en matière de charge à la surface de la Terre. Ensuite sont présentés les nombreux impacts que cette réponse isostatique peut avoir en retour sur la dynamique glaciaire, ainsi que les façons d'en rendre compte. Cette première partie aura montré tout l'intérêt de developper un modèle de Terre élaboré dont la physique à la base s'avère plus rigoureuse que dans les paramétrisations usuelles utilisées jusqu'à maintenant dans les modèles de glace. Les caractéristiques majeures ainsi que les équations à la base de ce modèle sont alors présentées en détail et révèlent le degré de sophistication de cette nouvelle approche. La richesse de données issues de mesures sur les zones de rebond postglaciaire comme la Fennoscandie sont un excellent moyen de valider pareil modèle. En comparant les résultats du rebond simulé avec ces données il est possible d'estimer le réalisme du modèle. Bien que l'ensemble de ces données ne soit pas simultanément reproduit de manière parfaite, le modèle montre cependant un comportement suffisamment réaliste pour pouvoir par la suite être couplé à un modèle de glace. Les résultats de cette incorporation sont alors comparés à ceux obtenus avec les paramétrisations usuelles au cours d'une simulation de la calotte Antarctique au cours du dernier cycle glaciaire. Ces résultats confirment le relativement bon comportement de notre modèle complet de Terre par rapport aux autres approches, ce qui tendrait à discréditer l'emploi de certaines de ces paramétrisations de la réponse terrestre dans les modèles évolutifs de calottes glaciaires.

Lors de la dernière déglaciation au Wisconsinien supérieur, la région de la côte est de la baie d'Hudson a subi une transgression marine (mer de Tyrrell) synchrone au recul des glaces continentales de l'Inlandsis laurentidien vers le centre du Québec nordique. Le relèvement glacio-isostatique postglaciaire a provoqué l'émersion d'importantes masses sédimentaires déposées dans cette mer. L'étude morphologique et stratigraphique de ces dépôts a permis de reconstituer les étapes de la déglaciation ainsi que les différentes phases du retrait de la mer. Dans la région de la rivière Biscarat, des indices d'une stabilisation temporaire de la marge glaciaire ont été retrouvés dans les collines côtières, à environ 7 km de la côte. Ce système est la continuité vers le nord de la ± Nastapoka drift belt ¿ dont l'âge a été évalué à 8 ka BP. Une autre pause ou ralentissement du front a été observée à l'est de la limite de l'invasion marine à l'intérieur des terres.

La basse vallée de la rivière du Gouffre est remblayée par des dépôts quaternaires appartenant au Wisconsinien et à l'épisode marin goldthwaitien. Le lit rocheux préglaciaire correspond à une faille orientée NO-SE. Le refroidissement climatique qui a sévi sur la région au Fini-glaciaire a permis une réavancée du front glaciaire associé à la récurrence deSaint-Narcisse, phase durant laquelle s'effectua la mise en place de matériel glacio-marin. La mer de Goldthwait a par la suite submergé la vallée jusqu'à la cote de 170 m. L'analyse des coupes stratigraphiques et les datations 14C obtenues à partir d'échantillons prélevés dans les sédiments fini-glaciaires et holocènes ont permis de définir la chronologie des événements depuis plus de 11 000 ans BP. À cette époque, un lobe de glace flottante s'est avancée jusqu'à l'embouchure de la baie Saint-Paul et vraisemblablement à l'intérieur de l'estuaire moyen. La photo-interprétation a permis de circonscrire le plus ancien glissement de terrain qui est survenu dans la vallée, lequel date l'abaissement du niveau marin au niveau de la terrasse de 20-35 m. Une discontinuité liée aux argiles basculées, non liquéfiées, est observable sous les dépôts de la terrasse de 8-10 m et au niveau du haut estran. L'existence d'un niveau marin plus bas que l'actuel vers 4600 ans BP a aussi été suggérée dans la baie. La courbe d'émersion post-glaciaire proposée pour ce secteur rend compte d'un déglacement qui semble s'être produit de l'est vers l'ouest de la vallée.