Sulphate (SO4), predominantly derived from sulphur (S)-bearing glacial sediments distributed widely across the Canadian Interior Plains, contributes to high groundwater salinity and can be detrimental to riparian and dry-land ecosystems, agricultural production, and water use. While previous researchers investigated SO4 distribution and dynamics in shallow groundwater at local scales (<1500 km2), we examine SO4 occurrence in groundwater at larger scales, and to depths of ∼150 m, considering variations in geology, glacial history, climate, and geochemical and hydrogeological settings in the Canadian province of Alberta. Sulphate concentrations in groundwater vary considerably, with 15 % of 139,130 samples above the 500 mg/L Canadian drinking water aesthetic objective. Analysis of δ34SSO4 and δ18OSO4 from 179 wells throughout Alberta shows SO4 is dominantly derived from oxidation of S-bearing material. At this large scale, marine shale bedrock subcrops, glacial ice flow directions, and climate control SO4 distribution. Groundwater contains less SO4 (<200 mg/L) in western Alberta compared to the east due to a lack of S-bearing bedrock, higher groundwater recharge rates, thinner glacial sediments, and increased groundwater circulation. Focusing on a region in south-central Alberta (Red Deer area), hydrogeological characteristics relating to SO4 formation, recharge, permeability, hydraulic gradient, and travel time are included in a principal component analysis to identify six groundwater groups at varying stages of evolution and SO4 concentrations depending on their hydrogeological setting. Low SO4 concentrations are associated with areas of high recharge, where SO4 has been leached from the sediments, or with more evolved bedrock groundwater where SO4 has been reduced or was recharged in pre-glacial times. High SO4 concentrations are found in areas of lower recharge and permeability, where flushing occurs more slowly. Combining the Red Deer area and province-wide analyses, we identify seven regions across Alberta with characteristic distribution and controls on SO4 occurrence in shallow groundwater.

Maps showing the thickness of sediments above the bedrock (depth to bedrock, or DTB) are important for many geoscience studies and are necessary for many hydrogeological, engineering, mining, and forestry applications. However, it can be difficult to accurately estimate DTB in areas with varied topography, like lowland and mountainous terrain, because traditional methods of predicting bedrock elevation often underestimate or overestimate the elevation in rugged or incised terrain. Here, we describe a machine learning spatial prediction approach that uses information from traditional digital elevation model derived estimates of terrain morphometry and satellite imagery, augmented with spatial feature engineering techniques to predict DTB across Alberta, Canada. First, compiled measurements of DTB from borehole lithologs were used to train a natural language model to predict bedrock depth across all available lithologs, significantly increasing the dataset size. The combined data were then used for DTB modelling employing several algorithms (XGBoost, Random forests, and Cubist) and spatial feature engineering techniques, using a combination of geographic coordinates, proximity measures, neighbouring points, and spatially lagged DTB estimates. Finally, the results were contrasted with DTB predictions based on modelled relationships with the auxiliary variables, as well as conventional spatial interpolations using inverse-distance weighting and ordinary kriging methods. The results show that the use of spatially lagged variables to incorporate information from the spatial structure of the training data significantly improves predictive performance compared to using auxiliary predictors and/or geographic coordinates alone. Furthermore, unlike some of the other tested methods such as using neighbouring point locations directly as features, spatially lagged variables did not generate spurious spatial artifacts in the predicted raster maps. The proposed method is demonstrated to produce reliable results in several distinct physiographic sub-regions with contrasting terrain types, as well as at the provincial scale, indicating its broad suitability for DTB mapping in general.

The shift in climate that occurred between the Last Glacial Maximum (LGM) and the Early Holocene (ca. 18–12 kyr BP) displayed rates of temperature increase similar to present-day warming trends. The most rapid recorded changes in temperature occurred during the abrupt climate oscillations known as the Bølling–Allerød interstadial (14.7–12.9 kyr BP) and the Younger Dryas stadial (12.9–11.7 kyr BP). Reconstructing ice sheet dynamics during these climate oscillations provides the opportunity to assess long-term ice sheet evolution in reaction to a rapidly changing climate. Here, we use glacial geomorphological inversion methods (flowsets) to reconstruct the ice flow dynamics and the marginal retreat pattern of the southwestern sector of the Laurentide Ice Sheet (SWLIS). We combine our reconstruction with a recently compiled regional deglaciation chronology to depict ice flow dynamics that encompass the time period from pre-LGM to the Early Holocene. Our reconstruction portrays three macroscale reorganizations in the orientation and dynamics of ice streaming followed by regional deglaciation associated with rapid warming during the Bølling–Allerød interstadial. Initial westward flow is documented, likely associated with an early set of ice streams that formed during the advance to the LGM. During the LGM ice streaming displays a dominant north to south orientation. Ice sheet thinning at ∼15 ka is associated with a macroscale reorganization in ice stream flow, with a complex of ice streams recording south-eastward flow. A second macroscale reorganization in ice flow is then observed at ∼14 ka, in which southwestern ice flow is restricted to the Hay, Peace, Athabasca, and Churchill river lowlands. Rates of ice sheet retreat then slowed considerably during the Younger Dryas stadial; at this time, the ice margin was situated north of the Canadian Shield boundary and ice flow continued to be sourced from the northeast. Resulting from these changes in ice sheet dynamics, we recognize a three-part pattern of deglacial landform zonation within the SWLIS characterized by active ice margin recession, stagnation and downwasting punctuated by local surging (terrestrial ice sheet collapse): the outer deglacial zone contains large recessional moraines aligned with the direction of active ice margin retreat; the intermediate deglacial zone contains large regions of hummocky and stagnation terrain, in some areas crosscut by the signature of local surges, reflecting punctuated stagnation and downwasting; and the inner deglacial zone contains inset recessional moraines demarcating progressive regional ice margin retreat. We attribute these macroscale changes in ice flow geometry and associated deglacial behaviour to external climatic controls during the Bølling–Allerød and Younger Dryas but also recognize the role of internal (glaciological, lithological and topographic) controls in SWLIS dynamics.

Basal gravel and sand mantling the bedrock floors of buried valleys throughout the Canadian Interior Plains, and conformably overlying proglacial lacustrine sediment, comprise the Empress Group. While previously conceptualized as stratigraphically equivalent deposits of preglacial rivers prior to the first and most extensive continental and montane glaciations, subsequent stratigraphic studies indicated that buried valley basal gravel must have been deposited between, or during, progressively more extensive continental glaciations and could not be stratigraphically equivalent throughout the buried valley network. However, in the general absence of formation-rank stratigraphic description of basal gravel units that might better inform the geologic history of the deposits, most workers simply consider Empress Group sediments time-transgressive. In this paper, we examine basal gravel at provincial and regional scales to understand its genesis and geologic history. At the provincial scale, we map basal gravel in three dimensions using a novel machine learning approach. At the regional scale, we formally define basal gravel formations at either end of the largest buried valley system in Alberta, which informs its glacial history and physiographic development and shows the importance of formation-rank stratigraphic description. Our results indicate that the buried valley network across Alberta is palimpsest in genesis and basal gravel units within it are chronostratigraphically intercalated between tills. We advocate that the Empress Group definition be extended across Alberta with modifications to improve its clarity and utility, and formally define the Old Fort, Unchaga, Ipiatik, and Winefred formations as part of the Empress Group.

Abstract Landforms and sediments on the palaeo–ice stream beds of central Alberta record glacitectonic raft production and subsequent progressive disaggregation and moulding, associated substrate ploughing, and grooving. We identify a subglacial temporal or developmental hierarchy that begins with incipient rafts, including en échelon hill-hole complexes, hill-hole pairs, and strike-slip raft complexes, all of which display patterns typical of transcurrent fault activation and pull apart. Many display jigsaw puzzle–style fragmentation, indicative of substrate displacement along shallow décollement zones and potentially related to patchy ice stream freeze-on. Their gradual fragmentation and smoothing produces ice flow-transverse ridges (ribbed moraine), hill-groove pairs, and paraxial ridge and groove associations. Initiator scarp and megafluting associations are indicative of raft dislodgement and groove ploughing, leading to the formation of murdlins, crag-and-tails, stoss-and-lee type flutings and drumlins, and Type 1 hogsback flutings. Downflow modification of rafts creates linear block trains (rubble stripes), stoss-and-lee type megaflutings, horned crag-and-tails, rubble drumlinoids, and murdlins, diagnostic of an immature palaeo–ice stream footprint. Lateral ice stream margin migration ingests disaggregated thrust masses to form ridged spindles, ladder-type morphologies, and narrow zones of ribbed terrain and Type 2 hogsback flutings, an assemblage diagnostic of ice stream shear margin moraine formation.

Conceptual ice stream land systems derived from geomorphological and sedimentological observations provide constraints on ice–meltwater–till–bedrock interactions on palaeo-ice stream beds. Within these land systems, the spatial distribution and formation processes of ribbed bedforms remain unclear. We explore the conditions under which these bedforms may develop and their spatial organization with (i) an experimental model that reproduces the dynamics of ice streams and subglacial land systems and (ii) an analysis of the distribution of ribbed bedforms on selected examples of palaeo-ice stream beds of the Laurentide Ice Sheet. We find that a specific kind of ribbed bedform can develop subglacially through soft-bed deformation, where the ice flow undergoes lateral or longitudinal velocity gradients and the ice–bed interface is unlubricated; oblique ribbed bedforms develop beneath lateral shear margins, whereas transverse ribbed bedforms develop below frontal lobes. We infer that (i) ribbed bedforms strike orthogonally to the compressing axis of the horizontal strain ellipse of the ice surface and (ii) their development reveals distinctive types of subglacial drainage patterns: linked cavities below lateral shear margins and efficient meltwater channels below frontal lobes. These ribbed bedforms may act as convenient geomorphic markers to reconstruct lateral and frontal margins, constrain ice flow dynamics, and infer meltwater drainage characteristics of palaeo-ice streams.

Over the three last decades, great efforts have been undertaken by the glaciological community to characterize the behaviour of ice streams and better constrain the dynamics of ice sheets. Studies of modern ice stream beds reveal crucial information on ice-meltwater-till-bedrock interactions, but are restricted to punctual observations limiting the understanding of ice stream dynamics as a whole. Consequently, theoretical ice stream landsystems derived from geomorphological and sedimentological observations were developed to provide wider constraints on those interactions on palaeo-ice stream beds. Within these landsystems, the spatial distribution and formation processes of subglacial periodic bedforms transverse to the ice flow direction – ribbed bedforms – remain unclear. The purpose of this study is (i) to explore the conditions under which these ribbed bedforms develop and (ii) to constrain their spatial organisation along ice stream beds. We performed physical experiments with silicon putty (to simulate the ice), water (to simulate the meltwater) and sand (to simulate a soft sedimentary bed) to model the dynamics of ice streams and produce analog subglacial landsystems. We compare the results of these experiments with the distribution of ribbed bedforms on selected examples of palaeo-ice stream beds of the Laurentide Ice Sheet. Based on this comparison, we can draw several conclusions regarding the significance of ribbed bedforms in ice stream contexts: * Ribbed bedforms tend to form where the ice flow undergoes high velocity gradients and the ice-bed interface is unlubricated. Where the ribs initiate, we hypothesize that high driving stresses generate high basal shear stresses, accommodated through bed deformation of the active uppermost part of the bed. * Ribbed bedforms can develop subglacially from a flat sediment surface beneath shear margins (i.e., lateral ribbed bedforms) and stagnant lobes (i.e., submarginal ribbed bedforms) of ice streams, while they do not develop beneath surging lobes. * The orientation of ribbed bedforms reflects the local stress state along the ice-bed interface, with transverse bedforms formed by compression beneath ice lobes and oblique bedforms formed by transgression below lateral shear margins. * The development of ribbed bedforms where the ice-bed interface is unlubricated reveals distinctive types of discontinuous basal drainage systems below shear and lobe margins: linked-cavities and efficient meltwater channels respectively. Ribbed bedforms could thus constitute convenient geomorphic markers for the reconstruction of palaeo-ice stream margins, palaeo-ice flow dynamics and palaeo-meltwater drainage characteristics.

Abstract. Conceptual ice stream landsystems derived from geomorphological and sedimentological observations provide constraints on ice-meltwater-till-bedrock interactions on palaeo-ice stream beds. Within these landsystems, the spatial distribution and formation processes of ribbed bedforms remain unclear. We explore the conditions under which these bedforms develop and their spatial organisation with (i) an experimental model that reproduces the dynamics of ice streams and subglacial landsystems and (ii) an analysis of the distribution of ribbed bedforms on selected examples of paleo-ice stream beds of the Laurentide Ice Sheet. We find that a specific kind of ribbed bedforms can develop subglacially from a flat bed beneath shear margins (i.e., lateral ribbed bedforms) and lobes (i.e., submarginal ribbed bedforms) of ice streams. These bedforms initiate where the ice flow undergoes high velocity gradients and the ice-bed interface is unlubricated. We suggest that (i) their orientation reflects the local stress state along the ice-bed interface and (ii) their development reveals distinctive types of subglacial drainage patterns below these two kinds of margins: linked-cavities and efficient meltwater channels respectively. These ribbed bedforms are thus convenient geomorphic markers to reconstruct palaeo-ice stream margins and constrain palaeo-ice flow dynamics and palaeo-meltwater drainage characteristics.

Abstract. Conceptual ice stream landsystems derived from geomorphological and sedimentological observations provide constraints on ice-meltwater-till-bedrock interactions on palaeo-ice stream beds. Within these landsystems, the spatial distribution and formation processes of ribbed bedforms remain unclear. We explore the conditions under which these bedforms develop and their spatial organisation with (i) an experimental model that reproduces the dynamics of ice streams and subglacial landsystems and (ii) an analysis of the distribution of ribbed bedforms on selected examples of paleo-ice stream beds of the Laurentide Ice Sheet. We find that a specific kind of ribbed bedforms can develop subglacially from a flat bed beneath shear margins (i.e., lateral ribbed bedforms) and lobes (i.e., submarginal ribbed bedforms) of ice streams. These bedforms initiate where the ice flow undergoes high velocity gradients and the ice-bed interface is unlubricated. We suggest that (i) their orientation reflects the local stress state along the ice-bed interface and (ii) their development reveals distinctive types of subglacial drainage patterns below these two kinds of margins: linked-cavities and efficient meltwater channels respectively. These ribbed bedforms are thus convenient geomorphic markers to reconstruct palaeo-ice stream margins and constrain palaeo-ice flow dynamics and palaeo-meltwater drainage characteristics.

The Neutral Hills Uplands of southern Alberta, Canada is an area of complex and varied glacial landforms dominated by glacitectonic compressional structures but also containing expansive areas of hummocky terrain and kame and kettle topography. It lies between the strongly streamlined trunks of the former Central Alberta (CAIS) and Maskwa palaeo-ice streams of the SW Laurentide Ice Sheet (LIS) and hence comprises an inter-ice stream regional moraine zone, constructed at around 15.5 cal ka BP. This study aimed to compile a regional map of the glacial geomorphology of central southeast Alberta in order to decipher the landform-sediment signatures of overprinted ice stream margins in terrestrial continental environments, and to refine the palaeoglaciological reconstructions for the southwest LIS. Detailed mapping from LiDAR and aerial imagery identifies distinctive glacial landsystems diagnostic of the partial overprinting of cross-cutting ice stream trunks and fast flow lobes. Widespread evidence of surge-diagnostic features indicates that the ice streams experienced repeated flow instabilities, consistent with the broader scenario of a highly dynamic and unstable SW LIS, characterised by markedly transitory and cross-cutting palaeo-ice streams. The inter-ice stream moraine zone is characterised by spectacular glacitectonic compression of bedrock, cupola hill construction and mega raft displacement but also displays evidence of multi-phase stagnant ice melt-out, where partially overprinted surge lobes advanced into large areas of buried glacier ice. Contemporaneous ice melting led to the widespread development of glacier karst and the production of eskers at a range of scales, the largest of which record deranged drainage patterns indicative of ice-walled channel sedimentation controlled by the regional bedrock slope towards the northeast. These process-form regimes have created a significant local relief that is a product of not only glacitectonic compression of bedrock but also the creation and melting of a melange of ice and bedrock/sediment blocks of variable ice volume, which are representative of former buried snout ice with a glacier karst system that was repeatedly proglacially thrust due to surging. Widespread evidence for subglacial channel cutting is likely strongly linked to the transitory, surging and cross-cutting nature of the palaeo-ice streams in the region, whereby ice streams switched on and surged in response to the build-up, migration and marginal outbursts of subglacial water reservoirs. In addition to the reduced basal friction caused by the low permeability of the Cretaceous bedrock, pressurized groundwater and potentially also shallow biogenic gas deposits were likely important to the process-form regimes of surging lobes of soft-bedded ice streams in a region where ice flow was against an adverse bed slope; a scenario that gave rise to a variety of enigmatic landforms such as doughnuts, doughnut chains, apparent blow-out features and possible till eskers, as well as glacitectonic mega-rafts. (C) 2019 Elsevier B.V. All rights reserved.

The North American Ice Sheet Complex (NAISC; consisting of the Laurentide, Cordilleran and Innuitian ice sheets) was the largest ice mass to repeatedly grow and decay in the Northern Hemisphere during the Quaternary. Understanding its pattern of retreat following the Last Glacial Maximum is critical for studying many facets of the Late Quaternary, including ice sheet behaviour, the evolution of Holocene landscapes, sea level, atmospheric circulation, and the peopling of the Americas. Currently, the most up-to-date and authoritative margin chronology for the entire ice sheet complex is featured in two publications (Geological Survey of Canada Open File 1574 [Dyke et al., 2003]; 'Quaternary Glaciations - Extent and Chronology, Part II' [Dyke, 2004]). These often-cited datasets track ice margin recession in 36 time slices spanning 18 ka to 1 ka (all ages in uncalibrated radiocarbon years) using a combination of geomorphology, stratigraphy and radiocarbon dating. However, by virtue of being over 15 years old, the ice margin chronology requires updating to reflect new work and important revisions. This paper updates the aforementioned 36 ice margin maps to reflect new data from regional studies. We also update the original radiocarbon dataset from the 2003/2004 papers with 1541 new ages to reflect work up to and including 2018. A major revision is made to the 18 ka ice margin, where Banks and Eglinton islands (once considered to be glacial refugia) are now shown to be fully glaciated. Our updated 18 ka ice sheet increased in areal extent from 17.81 to 18.37 million km(2), which is an increase of 3.1% in spatial coverage of the NAISC at that time. Elsewhere, we also summarize, region-by-region, significant changes to the deglaciation sequence. This paper integrates new information provided by regional experts and radiocarbon data into the deglaciation sequence while maintaining consistency with the original ice margin positions of Dyke et al. (2003) and Dyke (2004) where new information is lacking; this is a pragmatic solution to satisfy the needs of a Quaternary research community that requires up-to-date knowledge of the pattern of ice margin recession of what was once the world's largest ice mass. The 36 updated isochrones are available in PDF and shapefile format, together with a spreadsheet of the expanded radiocarbon dataset (n = 5195 ages) and estimates of uncertainty for each interval. (C) 2020 Elsevier Ltd. All rights reserved.

The extent of proglacial lakes following the initial separation of the southwest Laurentide Ice Sheet from the Cordilleran Ice Sheet and its eastward retreat from the Canadian Rocky Mountains has been reconstructed across regions of Alberta at a range of scales. However, to date, no studies have integrated all available geological information to produce a province-wide deglacial reconstruction that considers the evolution of proglacial lakes as components of the ice-marginal system. In this paper, we utilize a geologically constrained shoreline projection method with a high resolution digital elevation model to reconstruct the evolution of the ice-marginal system along the southwest LIS during the last deglaciation. This method provides new details on the configuration, volume, drainage history and routing of similar to 240 proglacial lakes as they migrated across Alberta and establishes a succession of paleogeographic reconstructions that can place geological evidence of regional and local ice-flow reorganizations into a spatiotemporal context. These reconstructions demonstrate that although the evolution of proglacial lakes was largely driven by the topography of the emerging landscape and the configuration of the ice margin, positive feedbacks in the ice-marginal system, particularly where margins transitioned from terrestrial to subaqueous settings played a major role in deglacial ice dynamics. Narrow, ribbon-shaped lakes that paralleled the ice margin induced relatively minor changes in style and rate of deglaciation, whereas the evolution of progressively larger, longer-lived lakes extending obliquely to the ice margin promoted surging and subsequent rapid retreat. Crown Copyright (C) 2019 Published by Elsevier Ltd. All rights reserved.

We examined seven landslide dams and their changes over time in the Peace River region of Canada. These landslides had subchannel rupture surfaces in glacial and glaciolacustrine sediments. We assessed the stability of the dams using 6 separate, morphometric-based stability indices (with a total of 10 stability thresholds).

Ice stream bed topography contains key evidence for the ways ice streams interact with, and are potentially controlled by, their beds. Here we present the first application of two‐dimensional Fourier analysis to 22 marine and terrestrial topographies from 5 regions in Antarctica and Canada, with and without mega‐scale glacial lineations (MSGLs). We find that the topography of MSGL‐rich ice stream sedimentary beds is characterized by multiple, periodic wavelengths between 300 and 1200 m and amplitudes from decimeters to a few meters. This periodic topography is consistent with the idea that instability is a key element to the formation of MSGL bedforms. Dominant wavelengths vary among locations and, on one paleo ice stream bed, increase along the direction of ice flow by 1.7 ± 0.52% km−1. We suggest that these changes are likely to reflect pattern evolution via downstream wavelength coarsening, even under potentially steady ice stream geometry and flow conditions. The amplitude of MSGLs is smaller than that of other fluvial and glacial topographies but within the same order of magnitude. However, MSGLs are a striking component of ice stream beds because the topographic amplitude of features not aligned with ice flow is reduced by an order of magnitude relative to those oriented with the flow direction. This study represents the first attempt to automatically derive the spectral signatures of MSGLs. It highlights the plausibility of identifying these landform assemblages using automated techniques and provides a benchmark for numerical models of ice stream flow and subglacial landscape evolution.

This paper presents the results of a detailed multidisciplinary study of the deformed bedrock and overlying Quaternary sediments exposed at the Mud Buttes in southern Alberta, Canada. This large, arcuate cupola hill is composed of intensely folded and thrust sandstones, siltstones and mudstones of the Cretaceous Belly River Group. Glacitectonism responsible for the development of this internally complex landform occurred at the margin of the newly defined Prospect Valley lobe of the Laurentide Ice Sheet. Analysis of the deformation structures reveals that construction of this landform occurred in response to at least two phases of south-directed ice sheet advance separated by a period of retreat. The first phase led to the formation of a forward propagating imbricate thrust stack leading to polyphase deformation of the Belly River Group. D1 thrusting led to the detachment of thrust-bound slices of bedrock which were accreted to the base of the developing imbricate stack. This process resulted in the structurally higher and older thrust-slices being progressively “back-rotated” (tilted), accompanied by D2 thrusting and folding. Further thrusting during D3 was restricted to the core of the Mud Buttes as the deforming sequence accommodated further compression imposed by the advancing ice. Minor oscillations of the ice margin led to localised brittle-ductile shearing (D4) of the bedrock immediately adjacent to the ice contact part of the thrust stack. The second phase of ice advance led to the accretion of a relatively simple thrust and folded sequence seen the northern side of Mud Buttes. The resulting composite thrust moraine was subsequently overridden by ice advancing from the NNW to form a dome-like cupola-hill. This readvance of the Prospect Valley lobe led to the formation of a thin carapace of Quaternary sediments mantling the Mud Buttes which include glacitectonite, till and an organic-rich clay-silt (?palaeosol).

ABSTRACTThis paper presents a reconstruction of the geometry, dynamics and flow pattern of the Laurentide and Cordilleran ice sheets along the Albertan portion of the eastern Canadian Rocky Mountains during the last glaciation. Sediment–landform associations relating to the evolution of these ice sheets document their initial convergence and mutual deflection across west‐central Alberta. The continued advance of Laurentide ice locally displaced Cordilleran ice westward into the Front Ranges, and deflected trunk ice emerging from the Athabasca River valley south‐eastwards along the Foothills. These convergent flow patterns remained major features in the regional geometry of both ice sheets throughout the Last Glacial Maximum. The onset of deglaciation was characterized by the northward retreat of the Laurentide Ice Sheet, which progressively unblocked valleys along the mountain front. The configuration of buttress removal initiated the Lateglacial expansion and sequential drawdown of the eastern Cordilleran Ice Sheet along fast‐flowing outlet glaciers. Subsequent regional retreat of Laurentide and Cordilleran ice was associated with the development of large proglacial lakes which amplified deglaciation by destabilizing the margins of both ice sheets, triggering extensive readvances due to periodic changes in water level associated with the opening and closing of outlets by ice margin fluctuations.