The role of glacial dynamics in the development of ice divides and the Horseshoe Intersection Zone of the northeastern Labrador Sector of the Laurentide Ice Sheet

The Labrador Sector of the Laurentide Ice Sheet is characterized by a complex network of ice divides and an extensive landform record outlining two broad and opposing ice flows that are separated by a narrow Horseshoe Intersection Zone (HIZ). This geomorphic system gave rise to contrasting reconstructions, which reflect uncertainties on the temporal evolution of ice divides, the nature and age of the main landform systems and the overall pattern of ice retreat during the last deglaciation. Here, we address these issues through systematic mapping of glacial landforms and ice-flow indicators in a large area of northeastern Quebec and Labrador. The application of cosmogenic (10Be and 26Al) dating to esker-fed glaciomarine deltas and different rock surfaces brings new constraints on the chronological framework and insights on the subglacial thermal regime. Our results outline four main landform assemblages, two of which relating to a major system of opposite ice flows that delineate the position of the eastern Ancestral Labrador ice divide. Our reconstruction shows that the HIZ occurs several tens of km to the west of this divide and that these two features are genetically distinct. The HIZ relates to the lateglacial development of ice streams showing a massive convergent ice flow into Ungava Bay. The upstream extent of this ice flow system shows a strong coupling with the outline of Ungava Bay drainage divide, suggesting that the configuration of the HIZ could be the expression of a topographic control on the late-glacial ice sheet dynamics. Landforms also indicate a significant shift in the deglaciation mode, which evolved from an overall warmbased to an areally-confined cold-based ice retreat near the center of the former ice mass - a change in dynamics that also played a role in the configuration of the HIZ. Overall, areas indicative of cold-based ice conditions are limited to highly elevated terrains, thus ruling out models arguing for an extensive cold-based ice cover over Ungava Bay throughout the last deglaciation.