Understanding how debris-covered glaciers respond to climate is necessary in order to evaluate future water resources and glacier flood hazard potential, and to make sense of the glacier chronology in mountain regions, In order achieve this, it is necessary to improve the current understanding of how surface debris affects glacier ablation rate, and to develop methods by which the ablation of debris-covered glaciers can be predicted under various climatic scenarios. This thesis develops a numerical surface energy balance model that uses simple meteorological data to calculate melt beneath a debris layer of given thickness and thermal characteristics. Field data from three contrasting sites demonstrate that the assumptions made within the model concerning the thermal properties of supraglacial debris are valid during most ablation conditions and that model performance is considerably better than previous models. Model results indicate that the effect of debris on melt rate is highly dependent on meteorological conditions. Under colder climates, thin debris can accelerate ice melt by extending the ablation period at both diurnal and seasonal scales. However, in milder mid- summer conditions, even a very thin debris cover inhibits melt rate compared to that of exposed ice. The new melt model is applied to produce the first quantified ablation gradients for debris- covered glaciers, and to model the evolution of ice surfaces under a debris layer of variable thickness. Modelled ablation gradients are qualitatively similar to hypothetical ones outlined previously, and quantitatively similar to those measured in the field. The ablation gradients are used to explore the factors affecting the response of debris-covered glaciers to climate change. Beneath a debris layer of variable thickness, the melt model produced ablation topography, as observed in the field, which underwent topographic inversion over time in response to debris redistribution. Debris thickness variability was found to cause calculated ablation rate to increase compared to that calculated using a mean debris thickness by one to two orders of magnitude, suggesting that melt calculations made on the basis of spatially averaged debris thickness may be inaccurate.

An improved understanding of the processes controlling the dynamics of the Greenland Ice Sheet is needed to enable more accurate determination of the response of the ice sheet to projected climate change. Meltwater produced on the ice sheet surface can penetrate to the bed and cause ice motion to speed up through enhanced basal sliding. However, the importance of coupled hydro-dynamics both to current ice sheet motion and future stability over the coming century is unclear. This thesis presents observations from the south-west Greenland Ice Sheet which improve our understanding of coupled hydro-dynamics. It commences with an investigation of the response of ice motion to exceptional meltwater forcing during summer 2012. Simultaneous field observations of ice motion (by GPS) and proglacial discharge show that, despite two extreme melt events during July 2012 and summer ice sheet runoff 3.9 s.d. above the 1958– 2011 mean which resulted in faster summer motion, net annual motion was slower than in the average melt year of 2009. This suggests that surface melt-induced acceleration of land-terminating regions of the ice sheet will remain insignificant even under extreme melting scenarios. The thesis then examines spatial variability in ice motion, in relation to an inferred subglacial drainage axis, using GPS and satellite radar observations from a land-terminating margin up to 20 km inland where ice is 800 m thick. Whilst spatial variability in subglacial drainage system configuration is found to control ice motion at short timescales, the proportional contribution of summer motion to annual motion is almost invariant. The structure of the subglacial drainage system does not therefore appear to significantly influence spatial variations in net summer speedup. Lastly, observations are made by applying feature tracking to 30 years of optical satellite imagery in a ~170 by 50 km area along the ice sheet margin (where ice reaches ~850 m thick) to examine whether coupled hydrology-dynamics affects inter-annual ice motion. Hydro-dynamic coupling resulted in net ice motion slowdown during a period of clear climate warming. Further increases in meltwater production may therefore reduce ice sheet motion. The thesis concludes that at land-terminating margins of the Greenland Ice Sheet, (1) larger annual meltwater volumes do not result in faster annual ice motion; (2) the detailed structure of the subglacial drainage network appears unimportant to the role of summer motion in determining annual motion; and (3) atmospheric warming over several decades has been accompanied by a slowdown in ice motion. As such, hydro-dynamic coupling is unlikely to form a significant positive feedback between surface melting and ice motion in response to projected climate warming. The wider relevance of these findings to tidewater systems requires further investigation.

Glacial-interglacial cycles drive changes in the discharge and sediment flux from the headwaters of glaciated basins, which are recorded by proglacial fluvial sediments. Linking braided river stratigraphy to the Quaternary climate record could indicate the control of climate-driven variations in discharge and sediment flux on fluvial processes, and the magnitude and frequency of past climate variations. New Zealand is a key location for investigating terrestrial records of Southern Hemisphere climate change. The Late Quaternary braided river deposits on the Canterbury Plains, South Island, New Zealand have formed over the last 400 ka. The coastal cliff marking the southeastern margin of the Canterbury Plains provides excellent exposure of fluvial sediments deposited during the last glacial period, from ~40 ka until the end of the Last Glacial Maximum (LGM) (~18 ka).Deposition at the modern coastline of the Canterbury Plains is interpreted in the context of the climate event stratigraphy for New Zealand, which requires a precise geochronology. This thesis demonstrates the first successful application of optically stimulated luminescence (OSL) dating to glaciofluvial sediments on South Island: a technique that was previously considered unsuitable for this region. Ages produced for the coastal stratigraphy range from 36.7 ± 2.9 to 18.2 ± 1.3 ka, indicating that deposition took place during the last glacial, with little or no postglacial aggradation. Three adjacent catchments on the southern Canterbury Plains - the Rakaia, Ashburton and Rangitata - have undergone glacial-drainage capture during the period represented by the coastal stratigraphy. During glacials, transfluent ice reversed the flow direction in several key tributaries, resulting in dramatic variations in effective drainage area of the Ashburton and Rangitata, and variations in discharge and sediment flux which are recorded in the stratigraphy of these catchments. The magnitude, timing and duration of drainage capture were quantified using the Plummer and Phillips (2003) glacier model. The Ashburton catchment increased to 160% of the modern effective drainage area when temperature change relative to modern conditions exceeded -6°C during the LGM. Meanwhile, the effective drainage area of the Rangitata decreased to 63%, and the Rakaia to 93%, reducing discharge from these catchments. Furthermore, glaciation dramatically affects the seasonality of the annual hydrograph. At four coastal sites, the fluvial stratigraphy was surveyed to investigate possible variations in depositional architecture, due to both climate variations, and glacial-drainage capture in the Ashburton and Rangitata. Unexpectedly, little vertical variation in depositional architecture was found, indicating that the deposits created by the braided rivers represent sediment transport during a similar set of flow (and by inference, climate) conditions. Laterally extensive erosional surfaces separating storeys of one or two flow depths in thickness, in combination with the OSL geochronology, suggest that the gravel-bed braided river stratigraphy primarily records a response to climate variations within glacial maxima, rather than on the scale of the glacial-interglacial cycle.

The aim of this study is to understand the extent, depth, magnitude and significance of subglacial sediment deformation. It will examine the role of this deformation in controlling glacier dynamics and landform generation in glaciers in general, and polythermal glaciers in particular. A detailed multi-dimensional approach is used to study recently exposed glacigenic sediments on the forefields of three polyglaciers in the Tarfala Valley, northern Sweden. Overridden fluted moraines and diamicton plains occur in each forefield. These palimpsest landforms consist of multiple subglacial traction tills. Flutes have quasi-regular geometry and about half of those studied have no initiating boulder. It is suggested here that flute formation by forced-mechanisms was superimposed on flute formation related to a topographically-induced flow instability. In each forefield the depth of the deforming-bed averaged between 0.2m and 0.6m thickness. Detailed clast fabric data suggest the diamicton plain is composed of thin layers of traction tills that accreted over time as the zone of deformation moved upwards. Laboratory shear box tests show that subglacial deformation required elevated pore-water pressures, which suggests deforming-bed conditions and flute formation were restricted to the temperate zones of polythermal glaciers. Magnetic fabrics suggest strain magnitudes were moderate (≤10), rather than the very high strain magnitudes (>10²) required by the deforming-bed model. The application of the micro-structural mapping technique demonstrates that subglacial deformation was multi-phase, heterogeneous, and partitioned into the softer and more easily deformed parts of the matrix. Consequently, deformation is controlled by variations in sediment granulometry and pore-water pressure, and is likely to have been spatially and temporally variable, a finding that supports the ice-bed mosaic model. The strain magnitudes and deforming-bed thickness suggest that soft-bed deformation did not exert a major control on glacier dynamics during the Little Ice Age advance.

The Cairngorm Mountains contain an outstanding assemblage of glacial landforms from both the deglaciation of the last British Irish Ice Sheet and the Younger Dryas readvance. Glaciers are recognised as sensitive indicators of past and present climate change and, thus, these landforms provide information about past climate and glacier-climate interaction that can be used to contextualise the present climate change. Previous interpretations have left doubt over the extent and style of the Younger Dryas readvance. In addition, the pattern and timing of deglaciation in the southern Cairngorms and, particularly, how local and external ice masses interacted is unclear. New geomorphological mapping from aerial images and fieldwork has been compiled in a GIS for a 600km2 area of the Cairngorm Mountains. This has allowed a complex pattern of ice-dammed lakes and local and regionally sourced ice margins to be reconstructed during the retreat of the last British Irish Ice Sheet. The mapping has been combined with new cosmogenic surface exposure ages taken from areas of hummocky moraine previously subject to differing age interpretations. The effect of moraine denudation on apparent 10Be ages has been checked by inverse modelling of the 10Be concentration vs. boulder height. The results indicate more extensive Younger Dryas glaciation, with glacier reconstructions and equilibrium-line altitudes (ELAs) comparable with the surrounding areas. Reconstruction of both valley and plateau-fed glaciers are presented, with modelling of local topoclimatic factors, such as radiation, avalanche and snow drifting, combined with precipitation gradients, explaining most of the variation within the glacier ELAs. The geomorphological evidence and palaeoclimatic inferences are important, alongside a growing number of palaeoglaciological studies, in acting as evaluation areas for current numerical models of ice-sheet growth and decay.

The reliance on freshwater released by mountain glaciers and ice caps demands that the effects of climate change on these thermally-sensitive systems are evaluated thoroughly. Coupling climate variability to processes of mass and energy exchange at the glacier scale is challenged, however, by a lack of climate data at an appropriately fine spatial resolution. The thesis addresses this challenge through attempting to reconcile this scale mismatch: glacier boundary-layer observations of meteorology and ablation at Vestari Hagafellsjökull, Iceland, and Storglaciären, Sweden, are related to synoptic-scale meteorological variability recorded in gridded, reanalysis data. Specific attention is directed toward synoptic controls on: i) near-surface air temperature lapse rates; ii) stationarity of temperature-index melt model parameters; and iii) glacier-surface ablation. A synoptic weather-typing procedure, which groups days of similar reanalysis meteorology into weather categories , forms the basis of the analytical approach adopted to achieve these aims. Lapse rates at Vestari Hagafellsjökull were found to be shallowest during weather categories characterised by warm, cloud-free weather that encouraged katabatic drainage; steep lapse rates were encountered in weather categories associated with strong synoptic winds. Quantitatively, 26% to 38% of the daily lapse-rate variability could be explained by weather-category and regression-based models utilizing the reanalysis data: a level of skill sufficient to effect appreciable improvements in the accuracy of air temperatures extrapolated vertically over Vestari Hagafellsjökull. Weather categories also highlighted the dynamic nature of the temperature-ablation relationship. Notably, the sensitivity of ablation to changes in air temperature was observed to be non-stationary between weather categories, highlighting vulnerabilities of temperature-index models. An innovative solution to this limitation is suggested: the relationship between temperature and ablation can be varied as a function of weather-category membership. This flexibility leads to an overall improvement in the simulation of daily ablation compared to traditional temperature-index formulations (up to a 14% improvement in the amount of variance explained), without the need for additional meteorological data recorded in-situ. It is concluded that weather categories are highly appropriate for evaluating synoptic controls on glacier meteorology and surface energetics; significant improvements in the parameterization of boundary-layer meteorology and ablation rates are realised through their application.

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.

Les vallées tunnels sont les structures de drainage d'eau de fonte les plus imposantes de l'environnement sous-glaciaire. L'inaccessibilité des milieux sous-glaciaires actuels limitent nos connaissances des mécanismes impliqués dans leurs formations, des paramètres contrôlant leur morphologie et de leurs influences sur la dynamique glaciaire. Ce travail présente une nouvelle approche expérimentale visant à mieux contraindre la formation et le fonctionnement des systèmes de vallées tunnels. Cette approche repose sur le développement d'un nouveau dispositif expérimental simulant la circulation d'eau pressurisée au sein d'un substrat poreux et perméable sous une couverture visqueuse. Les résultats des expériences menées avec ce dispositif ont permis de déterminer des relations étroites entre les paramètres du substrat et les modalités de l'écoulement d'eau sur la formation et la morphologie des vallées tunnels. Les résultats issus des expériences démontrent que ce dispositif permet de recréer des systèmes de vallées tunnels. L'étude de ces vallées expérimentales suggèrent que la topographie du substrat et la production d'eau de fonte joue un rôle primordial sur la genèse des vallées tunnels et sur leurs morphologies. Deux morphotypes de vallées tunnels ont pu être identifiés avec des morphologies et des mécanismes de formation indépendants. L'analyse de la dynamique de la calotte sus-jacente a permis de mettre en évidence un lien étroit entre le développement des vallées tunnels et la dynamique des "ice streams". L'évolution de la capacité de drainage des vallées tunnels semble contrôler la dynamique glaciaire en régulant le flux de glace transitant dans les "ice streams". / Tunnel valleys are major components of the subglacial meltwater drainage system. The inaccessibility of modern subglacial environments reduces our knowledge on the mechanisms involved in tunnel valleys formation, the parameters controlling their morphology and their influence on ice-sheet dynamics. This work presents a new experimental approach aiming to better assess the processes of tunnel valleys development. This approach relies on the development of a new experimental device simulating a pressurized water flow within a porous and permeable substratum underneath a viscous layer simulating the ice-sheet. The main results of the experiments conducted with this device have demonstrated the influence of both substratum properties and meltwater drainage on tunnel valleys formation and morphology. Using the device, we first manage to reproduce tunnel valley systems experimentally. Analyses conducted on these valleys experimental valleys suggest that the substratum topography and meltwater production play a key role on tunnel valleys genesis and morphology. Two tunnel valleys morphotypes have been identified, each one being charaterized by a unique morphology and mechanism of formation. Monitoring of the experimental ice sheet during tunnel valley formation shows close relationship between tunnel valleys development and "ice streams" dynamics. The evolution of tunnel valley drainage capacity seems to have a strong influence on ice sheet stability by regulating ice flux within "ice stream corridors".

It is now accepted that the last British-Irish Ice Sheet (BIIS) was highly dynamic and drained by numerous fast flowing ice streams. This dynamic nature combined with its maritime location made the BIIS sensitive to the rapid climate change that characterised the Last Glacial Interglacial Transition. Gaining an understanding of the behaviour of the BIIS at this time is important to explore the nature of forcing between ice sheets and climate. This thesis presents new chronological data relating to the deglaciation of the northwest sector of the BIIS (NW-BIIS) from onshore dating of moraines using cosmogenic exposure dating. This improved chronological framework is supported by offshore data in the form of a newly constructed Ice Rafted Detritus (IRD) record from the offshore sediment core MD95-2007. These data suggest that deglaciation commenced sometime after 18 ka and that the NW-BIIS was located close to the present day shoreline by 16 ka. Further provenance analysis of the IRD using U-Pb dating of detrital minerals demonstrates that during the Last Glacial-Interglacial Transition MD95-2007 was being supplied distal IRD from a source(s) to the west. The absence of diagnostic Scottish material suggests that after retreat to the coastline at 16 ka calving margins were not re-established during Greenland Interstadial 1. By combining these results with existing data relating to the deglaciation of the NW-BIIS it is possible to summarise the deglaciation history of the NW-BIIS from the continental shelf to mountainous source regions and compare this to numerical models of BIIS behaviour during this time. With a better understanding of the chronology of NW-BIIS retreat it is possible to relate the timing of initial deglaciation to possible forcing factors and gain a better understanding of the response of a marine based sector of an ice sheet to rapid climate change.

La Côte Nord de l’Estuaire et du Golfe du St. Laurent (Québec, Canada) est caractérisée par une série de complexes deltaïques mis en place en contexte de chute de Niveau Marin Relatif (NMR) forcée par le rebond glacio-isostatique, lors du retrait des marges de l’Inlandsis Laurentidien. L’étude de trois complexes deltaïques montre qu’un motif d’évolution morphostratigraphique contrôlé par le retrait des marges glaciaires prévalait au premier ordre. Le premier stade de sédimentation se caractérise par la mise en place de cônes d’épandage juxtaglaciaires sous-aquatiques. Leur localisation au front de la marge glaciaire fait que la répartition spatiale des corps sédimentaires montre localement un empilement en rétrogradation. Dès l’émergence d’un front glaciaire continental, des deltas proglaciaires se développent en contexte de chute de NMR, formant des lobes dont l’altitude décroît vers le bassin. Ces derniers sont initialement associés à un système fluviatile en tresse alimenté en sédiments glaciogéniques par les marges glaciaires en retrait. Malgré des taux de chute de NMR de plusieurs cm/an, aucune incision fluviatile n’est observée à ce stade et la dynamique de transfert est prédominante du fait des forts taux d’apports sédimentaires. Plus tard, à la suite du retrait des marges glaciaires hors des bassins versants, le remaniement paraglaciaire se développe du fait de la réduction drastique des apports en eaux et sédiments. Le système fluviatile, devenu méandriforme, s’encaisse dans les anciens dépôts deltaïques maintenant inactifs et les bordures de deltas sont remaniées par les processus littoraux (plages soulevées, terrasses marines). Cette étude révèle que la majorité du volume de ces complexes deltaïques (10-20 km3) est mis en place en quelques milliers d’années seulement, immédiatement après la déglaciation ; le remaniement paraglaciaire n’ayant contribué à l’accrétion deltaïque que très marginalement. La modélisation numérique (Dionisos) valide les différents forçages identifiés par l’analyse morphosédimentaire. Une synthèse des complexes deltaïques à l’échelle de toute la Côte Nord du St. Laurent a permis de catégoriser 21 complexes deltaïques en 4 scénarios d’évolution morphosédimentaire, directement liés à la dynamique de retrait de la marge glaciaire. Leur décryptage offre une clef de lecture originale pour l’historique du retrait des marges glaciaires de l’Inlandsis Laurentidien sur la période 12-7.5 ka cal BP. / The North Shore of the St. Lawrence Estuary and Gulf is characterized by the presence of deltaic complexes that were emplaced under falling Relative Sea Level (RSL) forced by the glacio-isostatic rebound, immediately after the retreat of the Laurentide Ice Sheet (LIS) margin. The study of three deltaic complexes reveals that a common morphostratigraphic evolution forced by the retreat of the LIS prevailed for the edification of these structures, reflecting the retreat of the LIS margin. The emplacement of subaqueous outwash fans beyond the retreating or stillstanding glaciomarine margin constitutes the first stage of this evolution. As tied to the ice-margin position, these fans are characterized in places by a backstepping stacking pattern, in spite of the forced regressive setting. From the emergence of a continental ice front, the proglacial deltaic system develops and forms lobes staged accordingly to the RSL fall. These deltaic systems actively prograde at that time because they are fed in glaciogenics by the retreating LIS margin through braided meltwater streams. In spite of the RSL fall reaching several centimeters per years, no fluvial entrenchment occurs mainly owing to the significant amount of sediment supply. Later, when the LIS margin retreats from the drainage basins of feeding rivers, fluvial systems experience a drastic drop in sediment supply that forced the interruption of the deltaic progradation and the onset of paraglacial reworking. The paraglacial reworking consists in the entrenchment of a meandering fluvial system within former deltaic deposits and shows the prevalence of shallow-marine processes (waves, tides) at the delta rim (raised beaches, marine terraces). This study reveals that the bulk of the deltaic volume (c.a. 10-20 km3) for each complex was emplaced in only a few thousands of years following the LIS margin retreat when the latter was still located in the drainage basin. The paraglacial reworking had a minor influence on the deltaic accretion. A forward stratigraphic model (Dionisos) is used to validate the variety of forcing as understood from the sedimentary analysis. A synthesis including 21 deltaic complexes of the St. Lawrence North Shore allowed the establishment of a fourfold categorization. This scheme of deltaic evolution was used in order to refine the position of the LIS margin retreating upland for a period ranging from 12 to 7.5 ka cal BP.