Athabasca River lowlands area

A 110-km-long, NNE-SSW aligned corridor of highly elongate subglacial bedforms (megaflutings), parallels the axis of the Tawatinaw and Athabasca Rivers, north of Rochester (study area 1, Fig. 2). Described by Richard (Reference Richard1979, Reference Richard1986, Reference Richard1987) as “drumlins, drumlinoids, flutings, grooves and furrows,” these landforms record vigorous ice flow along the northernmost CAIS (Evans et al., Reference Evans, Clark and Rea2008; cf. Margold et al., Reference Margold, Stokes and Clark2015a, Reference Margold, Stokes, Clark and Kleman2015b, Reference Margold, Stokes and Clark2018). In detail (Fig. 3), they are discontinuous features, locally superimposed upon a number of largely elongate and smoothed mounds that extend obliquely to the main megafluting direction. The megaflutings display highly variable spacing, lengths, and widths. Plan forms also vary, with some comprising wider stoss- and tapering lee-ends but most displaying tapering stoss- and lee-ends typical of spindle drumlins (cf. Gluckert, Reference Gluckert1973; Shaw, Reference Shaw1983). Some features also display multiple flow-parallel ridges and/or summit grooves, and in places densely spaced lineations give the appearance of ridges and grooves of similar amplitudes. Large blocks (rafts) occur either in areas between ridges or embedded in the stoss-ends; in at least one case, elongate mounds (likely rafts) within partially grooved drumlins lie downflow of elongate depressions (Fig. 3, insets). North of the Athabasca townsite, megaflutings reported by Shaw et al. (Reference Shaw, Faragini, Kvill and Rains2000) clearly initiate at a scarp, in this case the margin of the modern Athabasca River valley. This relationship indicates that the megaflutings are genetically related to the scarp, which predates these streamlined landforms and hence potentially represents an incipient palaeo-valley margin within the vicinity of the ancestral Athabasca River (Fig. 3).

Figure 3. (color online) Annotated LiDAR digital elevation model (top) and geomorphology map (bottom) of the megafluting field located near the Athabasca townsite, showing the locations of exposures studied as well as that of Shaw et al. (Reference Shaw, Faragini, Kvill and Rains2000).

Other examples of megaflutings initiating at scarps (either along valley margins or the edges of uplands) occur around the Athabasca townsite (Fig. 3) and on the south shore of Calling Lake, 25 km northeast (Fig. 4). Indeed, such scarp-fluting associations are very common in the region, particularly in areas where bedrock occurs at, or close to, the modern land surface (Atkinson et al., Reference Atkinson, Pawley, Andriashek, Hartman, Utting and Atkinson2020), and hence we propose the term “initiator scarp” to describe such occurrences. Many megaflutings do not initiate directly at the scarp but rather at various distances down-ice, where they often appear to have conical or spindle-shaped stoss-ends and taper southwards (stoss-and-lee flutings/drumlins). Also apparent are flat-floored grooves separating flat-topped ridges, the latter locally exhibiting grooved summits. Northeast of Calling Lake, the width and amplitude of all of the megafluting ridges decrease over a distance of 25 km, where they merge into tracts of prairie mounds and/or doughnuts and other low-amplitude hummocks (Figs. 3 and 4). Ridge widths also tend to alternately widen and narrow, and many groove floors undulate due to the occurrence of chains of straight-sided depressions (Fig. 4). At south Calling Lake, the initiator scarp has been heavily subglacially moulded to form a band of short drumlins with surface ridge and groove patterns and interspersed flutings (Fig. 4). Drumlins generally taper out/flatten over 3 km, with overall streamlining (fluting) disappearing ca. 8 km downflow of the scarp. Many drumlins appear rectilinear at their stoss-ends, hence their tendency to flatten in long profile rather than taper in plan form, potentially indicative of raft cores.

Figure 4. (color online) Annotated LiDAR digital elevation model of megafluting fields, drumlins, and initiator scarps near Calling Lake. The Google Earth image (inset) shows an enlarged area of the western flutings and illustrates examples of the chains of straight-sided depressions in the flat-floored grooves.

The complex sedimentology of the megaflutings was reported by Shaw et al. (Reference Shaw, Faragini, Kvill and Rains2000) based on exposures through the Schiller site, north of Athabasca town (Figs. 3 and 5). The main lithofacies were described as: (1) a mélange, comprising sand and gravel intrabeds within a matrix of fine-grained diamicton; and (2) interbedded deposits, characterised by diamicton and plane-bedded and cross-laminated sorted sediments. Glacitectonic structures included widespread folding, localised but densely spaced normal faults, and diapiric or injection features. Clast macrofabrics from the diamictons revealed a wide range of orientations with weak to moderately strong S₁ eigenvalues (0.44–0.71; see Fig. 5), although two broad alignments are visible at N-S and W-E. Overall, Shaw et al. (Reference Shaw, Faragini, Kvill and Rains2000) considered a diapiric mélange rather than shear zone mélange (sensu Orange, Reference Orange1990) origin for these sediments and consequently dismissed megafluting production by subglacial deformation, preferring instead an erosional origin due to subglacial megaflooding. Additionally, the initiation of the megaflutings at an upflow-facing bedrock escarpment was explained as the product of longitudinal vortices being generated at this escarpment to form high-velocity, eroding streaks in subglacial sheet flows. The widespread deformation was interpreted as a product of the gradients generated between highly pressurised groundwater and low-pressure subglacial cavities.

Figure 5. Sketch of sedimentary structures and clast macrofabric data (including one palaeocurrent sample from the sands) at Schiller excavation site (modified after Shaw et al., Reference Shaw, Faragini, Kvill and Rains2000).

New exposures through a fluting ridge south of the Athabasca townsite (Fig. 3) provide further ridge-parallel and ridge-transverse sections displaying complexly deformed diamictons and stratified sediments and mélange. This tripartite sequence (Fig. 6) comprises a lower pseudo-stratified, clay-rich diamicton (LF1); a middle gravel-rich, sandy diamicton containing irregular, deformed lenses of silt/clay diamicton (LF2); and an upper massive, matrix-supported diamicton (Dmm; LF3).

Figure 6. (color online) Details of the south Athabasca townsite section in the megafluting complex. (a and b) The main section face orientated parallel to the fluting long axis; (c) the exposure orientated oblique to fluting long axis with associated clast macrofabrics. Dmm–massive, matrix-supported diamicton. MLA–mean lineation azimuth.

In detail, LF1 characterises as a Type I/II mélange (sensu Cowan, Reference Cowan1985; cf. Evans, Reference Evans2018), wherein deformed intraclasts of stratified sand and gravel within a clay-rich diamicton have been highly folded and thrust. Also apparent are irregular, highly deformed rafts of lower Cretaceous sandstone, the closest available outcrops of which occur ~100 km to the north (Pelican Formation; Prior et al., Reference Prior, Hathway, Glombick, Pana, Banks, Hay, Schneider, Grobe, Elgr and Weiss2013). The numerous intraclasts/rafts give the impression of pseudo-stratification in a diamictic mélange. The geometry of the thrusts and asymmetrical folds record a sense of shear towards the S/SE, consistent with an applied stress direction from the N/NW (Fig. 6b). Towards the top of LF1, forming a ≤0.30-m-thick deformed contact with overlying LF2, is a zone of more attenuated lenses and recumbent rootless folds more similar to a Type IV mélange (sensu Cowan, Reference Cowan1985; cf. Evans, Reference Evans2018). The diamicton of LF2 has a coarser sand/gravel matrix and contains deformed intraclasts of poorly sorted cobble to pebble gravel as well as clayey/silty diamicton similar to the matrix of underlying LF1 (Type II/III mélange). A clast macrofabric from LF2 displays a moderately strong (S₁ = 0.65; Fig. 6c), WNW-ESE clast alignment towards 120°. This mélange (LF2) is unconformably overlain by LF3 (Dmm), the contact and outcrop of which is well exposed in upper sections along the sides of a track cut oblique to the main face. The Dmm is ≤1 m thick, compact, and locally fissile, with a moderately strong (S₁ = 0.58) clast macrofabric orientated 353° (Fig. 6c). The two general alignments of the LF2 and LF3 macrofabrics replicate most of the data presented by Shaw et al. (Reference Shaw, Faragini, Kvill and Rains2000), with broadly N-S alignments matching the NNE-SSW orientation of the megaflutings, and the broadly W-E alignments representing fluting-transverse signatures (cf. Figs. 3, 4, and 6c).

Contrary to the genesis proposed by Shaw et al. (Reference Shaw, Faragini, Kvill and Rains2000), we advocate a raft-ploughing origin for the Athabasca megaflutings. At sites bordering the Athabasca River valley in particular (Figs. 3, 5, and 6), this is consistent with both: (1) the landforms, specifically the close association of initiator scarps, rafts, elongate hill-hole pairs, and stoss-and-lee flutings with rectilinear stoss bumps; and (2) the internal sedimentology, comprising mélanges of fragmented bedrock and glacigenic sediments representative of a vertical continuum from Type A to Type B glacitectonites (sensu Benn and Evans, Reference Benn and Evans1996; Evans et al., Reference Evans, Phillips, Hiemstra and Auton2006; Evans, Reference Evans2018) capped by a more homogenised Dmm or subglacial traction till. Deformation of this sequence was initially driven by ice advancing from the northwest, leading to the construction of a series of ice flow-transverse ridges (likely thrust ridges) that were subsequently overridden to form the NE-SW aligned streamlined hills. These ridges were then ploughed obliquely by erodent rafts detached from the initiator scarps, giving rise to increasingly intense glacitectonite production and homogenisation to create the streamlined till carapace beneath the southward-flowing ice stream. A clear pattern of downflow dissipation of glacial bedforms represents raft fragmentation and the propagation of lee-side till deformation. Also important are the chains of straight-sided depressions in the flat-floored grooves (Fig. 4). These features resemble chatter marks observed at smaller scales across subglacially eroded bedrock (e.g., striations and grooves; Laverdiere et al., Reference Laverdiere, Guimont and Pharand1979, Reference Laverdiere, Guimont and Dionne1985). Consequently, we speculate that these features represent soft-bedded mega-chattermarks created by the stick-slip motion of rafts being dragged across a deformable substrate to create megafluting, a concept that requires further testing.

Otter–Orloff Lakes area

Located west of Calling Lake (Fig. 2b), this area comprises substantial hill-hole pairs (sensu Aber et al., Reference Aber, Croot and Fenton1989) as well as numerous smaller hills with rectilinear edges separated by water-filled depressions (Fig. 7), all developed in an area of thick Quaternary sediment (≤45 m). The most obvious hill-hole pairs occur around Otter–Orloff Lakes, where a series of rectangular NNW-SSE aligned, water-filled depressions lie to the north of a 37 km² complex of thrust masses (hills) derived from the depressions. Diagnostic of a thrust origin are arcuate inset subsidiary ridges or crenulations on the summits of the main hills, as well as linear hill margins and internal structural lineaments marked by channels, both aligning with the predominantly straight margins of the source depressions; the straight margins represent wrench faults initiated by stick-slip displacement. These parallel fault structures have also locally influenced/controlled the morphology of the adjacent thrust masses, resulting in the linear nature of the thrust hill margins. Paraxial ridges also occur at the lateral edges of some thrust masses (Fig. 7). Glacial overriding of these features is indicated by more hummocky to smoothed rather than crenulate surfaces and the continuation of eskers from the holes and onto the proximal slopes, where they connect with meltwater channels incised through the hills (cf. Moran et al., Reference Moran, Clayton, Hooke, Fenton and Andriashek1980; Fig. 7). The displacement direction of the hills is towards the SSE, parallel to, but partially crosscutting, the westernmost megaflutings of the Athabasca Ice Stream (the up-ice component of the CAIS). This indicates that the features were displaced during the early shutdown of the western edge of the CAIS. Elongate assemblages of prairie doughnut chains and eskers on the distal slopes and down-ice flow of the main thrust mass record the release of pressurised groundwater due to glacitectonic stress/ice sheet loading; this was linked to subglacial drainage networks at the junction with the active part of the CAIS immediately to the east (cf. Moran et al., Reference Moran, Clayton, Hooke, Fenton and Andriashek1980; Evans et al., Reference Evans, Atkinson and Phillips2020).

Figure 7. (color online) Annotated LiDAR digital elevation model (top) and geomorphology map (bottom) of the Otter–Orloff Lakes and Calling Lake area, showing the major landforms.

To the west of the Otter-Orloff Lakes hill-hole pairs, a similar but less prominent thrust moraine arc forms the easternmost extent of a 14-km-wide corridor of streamlined ridges and small hills. This moraine arc records a late-stage readvance from the WNW (hereafter the Otter Lake readvance), which appears to have partially modified the earlier thrust mass excavated from Otter Lake by the shutdown of the Athabasca Ice Stream. This is manifest by overprinted surface crenulations indicative of thrusting and folding, which are aligned E-W (oldest) on the front of the thrust mass and SW-NE (youngest) on the back (Fig. 7). The streamlining created by the readvance comprises subtle fluting that has modified underlying arcuate to crenulated ridges and chaotically distributed hills. Directly west and northwest of Otter Lake, the ridges and hills are less modified, appearing as rectangular, flat-topped uplands (≤0.4 km wide × 1.4 km long), in some cases displaying crenulated summits interspersed with rectilinear depressions. The linear edges of the hills parallel those of a prominent escarpment located immediately to the west, giving the impression that they have been horizontally detached and then extensionally displaced along a shallow décollement zone towards the western side of the Otter Lake depression (Fig. 7). We hereafter adopt the term “en échelon hill-hole complex” (EHHC) for this landform assemblage and relate the preservation of individual thrust masses here to late-stage production, in this case the result of either: (1) the basal submarginal freeze-on of the ice lobe responsible for the Otter Lake readvance; or (2) their detachment due to transient variations in groundwater pressure during the rapid ice advance, possibly also explaining the presence of doughnut chains and eskers (cf. Moran et al., Reference Moran, Clayton, Hooke, Fenton and Andriashek1980).

Jenkins Lake-Grosmont area

Remarkable features on the escarpment flanking the western edge of the Athabasca River lowlands (Figs. 2b and 8) relate to a late-stage, curvilinear, NW-SE pattern of ice streaming recorded in a strongly fluted terrain that feeds into the more dominant N-S flowing CAIS (Fig. 8). A final WNW-ESE linear flow crosscuts the northernmost part of this curvilinear pattern but similarly terminates at the footprint of the CAIS. The suture zone between these easterly flowing ice streams and the CAIS is demarcated by a shear margin moraine (sensu Dyke and Morris, Reference Dyke and Morris1988; Stokes and Clark, Reference Stokes and Clark2002) some 20 km long and grading into a northward-widening, elongate zone of thrust masses and associated up-ice source depressions. This shear margin moraine separates the present-day Athabasca River valley and the glacifluvial landforms and large depressions containing Island and Baptiste Lakes. Located immediately to the west of this heavily pitted glacifluvial complex and its esker ridges is a prominent rectilinear scarp. To the west of Island and Majors Lakes, right-angled re-entrants that cut back into the scarp are floored by shallow depressions, one of which is now occupied by Ghost Lake. Further south, Baptiste Lake appears to occupy a similar depression but also contains islands and peninsulas composed of crenulated thrust masses; indeed, a larger thrust mass occurs at the southern end of Baptiste Lake, suggesting that its twin lake basins are two inset hill-hole pairs (Fig. 8). A further hill-hole pair occurs below the right-angled re-entrants, where the highly linear Majors Lake has been excavated to form a 5-km-long and 0.75-km-wide groove, ending at an upstanding mound of the same width. This groove also displays a clear paraxial ridge comprising four slightly offset segments, likely due to the progressive narrowing of the thrust mass while it continued to excavate the groove (Fig. 8, inset). We classify this type of feature a “paraxial ridge and groove association” (PRGA).

Figure 8. (color online) The glacial geomorphology of the Laurence Lake-Grosmont area to the northwest of Athabasca. Main panel is annotated LiDAR digital elevation model of the major landforms of the area. Bottom left inset shows the details of the Majors Lake paraxial ridge and groove association (PRGA) and top right inset shows detailed mapping of the glacitectonic landforms in the Jenkins Lake–Grosmont area.

The most obvious hill-hole pairs and associated PRGAs lie to the north of Grosmont, where narrow thrust masses have been displaced southeastward, carving deep grooves into the uplands and surrounding scarp. The most spectacular of these is the Jenkins Lake hill-hole pair, where a 6.5 km² crumpled thrust mass originating from Laurence Lake has been displaced over 11 km to form a 45-m-deep mega-groove (Atkinson et al., Reference Atkinson, Utting and Pawley2018; Fig. 8), effectively creating a hill-groove pair. Along the groove edge are dislocated compressive ridges curving downflow and truncated by small flutings or paraxial ridges. Emerging from the up-ice and distal ends of the thrust mass are multiple eskers that coalesce into a single ridge immediately to the southwest and then wind 15 km southwards into the Baptiste Lake depression. This is reminiscent of the thrust mass–esker relationships identified by Bluemle and Clayton (Reference Bluemle and Clayton1984), due to transient variations in porewater pressures associated with thrusting, with eskers forming during water escape and depressurisation of the hydrogeologic system, leading to the cessation of thrust mass displacement.

Further groove-shaped hill-hole pairs north of Jenkins Lake indicate that thrusting was initiated when ice flow shifted slightly from SE to ESE in an onset zone of a tributary of the CAIS. Also apparent here is the tendency for the thrust masses to have partially fragmented downflow, producing partially streamlined linear rubble trains (see Lac la Biche-Beaver Lake area below). The compressive, constructional nature and excellent preservation of the Laurence-Jenkins Lakes hill-hole pair and other smaller examples to the north indicate that they were late-stage features, likely developed beneath a thinning ice margin. This contrasts with the more slablike thrust masses (rafts) that were detached from the rectilinear source scarps, the re-entrants along which are thought to be controlled by the intersection of fracture sets within the underlying bedrock (cf. Ozoray, Reference Ozoray1972; Misra et al., Reference Misra, Slaney, Graham and Harris1991; Pana and Waters, Reference Pana and Waters2016). These bedrock structures were exploited during subsequent erosion of the source scarp, forming the prominent, linear “sidewalls” and “backwalls” (Fig. 8) left behind after raft removal. Two genetic scenarios are possible here. First, rectilinear re-entrants within down-ice flow-facing scarps are similar to headwalls of translational mass movement features, ubiquitous in Alberta, and hence the detachment of incipient rafts could have been initiated by shallow slope failure before ice advance. Second, and not mutually exclusively, we suggest a two-stage raft-detachment process comprising: (1) basal freeze-on during ice stream thinning and shutdown; followed by (2) dislocation and lateral displacement of thin rafts across near-surface décollements, likely aided by variations in porewater pressure beneath these areas of the ice stream when its fast flow switched back on.

Lac la Biche-Beaver Lake area

The glacial landforms in the Lac la Biche and Beaver Lake area (Figs. 2a and 9) record vigorous flow within the Lac la Biche and Seibert Ice Streams (cf. Andriashek and Fenton, Reference Andriashek and Fenton1989), which form the northernmost branches of the onset region of the Buffalo-James Lobe Ice Stream system (Ross et al., Reference Ross, Campbell, Parent and Adams2009; Margold et al., Reference Margold, Stokes and Clark2015a, Reference Margold, Stokes, Clark and Kleman2015b). Regional orientations of subglacial lineations and fracture patterns in the thrust masses and rafts record spatiotemporal switching in the relative dominance of competing ice streams in the Lac la Biche area, previously reported by Andriashek and Fenton (Reference Andriashek and Fenton1989) and modified here based upon newly recognised landform patterns. Specifically, an early NW-SE flow of the Lac la Biche Ice Stream was replaced by the N-S flow of the Seibert Ice Stream as its margin migrated westwards. Later narrowing of the Seibert Ice Stream allowed the re-establishment of a vigorous but similarly narrow southeasterly flowing Lac la Biche Ice Stream, as defined by the limit of its most recent streamlining directly south of Beaver Lake. Between the most recent, highly streamlined parts of these ice stream trunks, immediately south and east of Lac la Biche and around Beaver Lake, a variety of glacitectonically displaced rafts lie only relatively short distances from their source depressions and appear variably fragmented downflow (Fig. 9b and d). In contrast, the highly streamlined trunks of the Lac la Biche and Seibert Ice Streams contain more intensely smoothed source depressions, thrust masses, and linear block trains or rubble stripes (Evans et al., Reference Evans, Atkinson and Phillips2020; cf. aligned rubble of Atkinson et al., Reference Atkinson, Utting and Pawley2018), some of which have drumlinoid plan forms and hence are termed “rubble drumlinoids” (Fig. 9a and c). In combination, the landforms in the two areas of highly and less streamlined ice stream footprints can be used to demonstrate the spectrum of ice-marginal to potentially subglacial raft detachment, substrate moulding, reworking, and fragmentation of earlier detached rafts, and their role as downflow dispersing erodents to create MSGLs.

Figure 9. (color online) Annotated LiDAR digital elevation model of the glacial landforms in the area around Lac la Biche and Beaver Lake. Inset boxes show: (a) streamlined inner zone of the Lac la Biche Ice Stream trunk; (b) western boundary zone of the Seibert Ice Stream; (c) streamlined rubble of the inner zone of the Lac la Biche Ice Stream trunk; (d) the Beaver Lake en échelon hill-hole complex (EHHC).

Early-stage raft-detachment patterns are evident in the complex outline of Beaver Lake, which is a product of a series of arcuate, en échelon, flow-transverse, and largely flat-topped ridges (Fig. 9d). Although curvilinear in plan form, the cliffs of the ridges are composed of rectilinear segments, with each ridge giving the appearance that they all fit together like a jigsaw puzzle; that is, the land surface has been fractured by extension and transported southwards over a shallow décollement zone. The surfaces of the ridges exhibit faint N-S orientated flutings and are mantled by geometric ridge networks, with both landform types being visible but not directly traceable from one ridge to another and onto the land surfaces surrounding the lake. Hence, these features are interpreted as an EHHC, developed after ice streaming had created flutings and geometric ridge networks. These subglacial landforms are diagnostic of surging activity and are recognised on palaeo–ice stream beds across the region (cf. Sharp, Reference Sharp1985a, Reference Sharp1985b; Evans and Rea, Reference Evans and Rea1999, Reference Evans, Rea and Evans2003; Evans et al., Reference Evans, Lemmen and Rea1999, Reference Evans, Clark and Rea2008, Reference Evans, Storrar and Rea2016, Reference Evans, Atkinson and Phillips2020; Ó Cofaigh et al., Reference Cofaigh, Evans and Smith2010). The subsequent fracturing of this subglacial surface to form an EHHC most likely represents the freeze-on along the narrowing outer western boundary of the Seibert ice stream after a phase of surging.

Directly east of Lac la Biche and Beaver Lake, the western boundary of the Seibert Ice Stream can be further subdivided based on the intensity of streamlining; highly elongate subglacial bedforms in the east are separated by an abrupt boundary from thrust masses and rubble terrain and rubble stripes in the west (Fig. 9b). The types of thrust mass in this area include arcuate, flat-topped ridges with a source scarp in the north (an EHHC), numerous examples of further-travelled and more highly fragmented rafts and their source scarps (e.g., immediately east of Beaver Lake), hill-hole pairs forming the sub-basins of southern Beaver Lake and south Elinor Lake, and closely spaced rubble stripes to the south of Beaver and Elinor Lakes. Importantly, these rubble stripes resemble pre-existing thrust masses that have become overridden and fragmented subglacially as they were displaced variable distances southwards as a result of continued ice advance from the NW/NNW (see Fig. 9b). The linearity of the margins of the rubble stripes likely relates to the original morphology of the thrust masses and their straight-sided source hollows (see Fig. 7). It is possible that the flow-parallel arrangement of the rubble stripes is indicative of the subglacial dispersal of the blocks.

In some cases, the narrow thrust masses have been fragmented along their displacement track to form aligned chains of transverse ridges reminiscent of the “ladder” morphology identified for some ribbed moraine by Dunlop and Clark (Reference Dunlop and Clark2006). Importantly, these ladder morphologies and rubble stripes change abruptly on their eastern margin into highly streamlined forms (flutings, spindle drumlins, stoss-and-lee drumlins), thereby demarcating two zones on the Seibert Ice Stream footprint comprising: (1) a relatively narrow but vigorous flow zone represented by the highly streamlined forms to the east; and 2) a boundary zone to the west, in which less subglacially streamlined features such as EHHCs, fragmented rafts, hill-hole pairs, ladder morphologies, and rubble stripes record freeze-on and displacement of rafts over a relatively shallow décollement due to late-stage ice stream narrowing and shutdown.

The relatively highly streamlined inner zone of the Lac la Biche Ice Stream trunk (Fig. 9a) displays a number of thrust mass source depressions, identifiable by: 1) their rectilinear (fracture-controlled) boundaries; and 2) rectangular, parallel sub-basins defined by flutings that likely originated as paraxial ridges. In some cases (e.g., Burnt Lake), the streamlined and fluted thrust masses occur immediately down-ice flow of the depressions, thereby constituting hill-hole pairs. However, more common are highly fragmented thrust masses (e.g., Missawawi and Kinosiu Lakes), manifest either as single, partially streamlined blocks (small cupola hills) and stoss-and-lee drumlins/flutings located at various distances from their source depressions or as flutings with hummocky long profiles, hereby termed “Type 1 hogsback flutings”; the latter appear to be the product of streamlined rubble stripes. Features also widespread in this inner zone of intensive ice stream streamlining are blocks or thrust masses that appear to have been: (1) ploughed through deformable subglacial material to create a prow and extended limbs or paraxial ridges; and (2) lodged, such that deformable subglacial material has been advected downflow to form two extended limbs, akin to the horned crag-and-tails of Jansson and Kleman (Reference Jansson and Kleman1999). A quarry within the thrust mass of a horned crag-and-tail near the Lac la Biche townsite reveals a stratigraphy composed of stacked glacitectonites (Fig. 10). Here, an attenuated, 40-m-long, highly fissile mudstone raft with small-scale tectonic laminae and slickensides occurs within a sand and gravel glacitectonite, the uppermost contact of which grades abruptly into a massive, matrix-supported diamicton representative of till production by cannibalisation (cf. Evans et al., Reference Evans, Phillips, Hiemstra and Auton2006; Evans, Reference Evans2018).

Figure 10. (color online) Details of a quarry exposure in a horned crag-and-tail near the Lac la Biche townsite.

Whitefish Lake-Goodfish Lake area

The three lake basins containing Whitefish, Goodfish, and Garner Lakes are defined by prominent rectilinear outlines, the largest and most continuous of which is a NNW-SSE trending scarp that stretches more than 30 km (Fig. 11). They are located in an ice stream confluence zone (Andriashek and Fenton, Reference Andriashek and Fenton1989) on the southern edge of the Lac la Biche fluting field (fast-flow trunk zone of Ice Stream 2/2A; Ó Cofaigh et al., Reference Cofaigh, Evans and Smith2010; Norris et al., Reference Norris, Evans and Cofaigh2018). The scarp and its fronting basins constitute a series of depressions that appear to be the source of the glacitectonically displaced masses located between 7.5 and 9 km to the southeast. The rectilinear shapes of the main and subsidiary depressions are likely controlled by bedrock fracture sets (joints) orientated NW-SE (ice flow parallel) and NE-SW (ice flow transverse; Ozoray, Reference Ozoray1972; Misra et al., Reference Misra, Slaney, Graham and Harris1991; Pana and Waters, Reference Pana and Waters2016). The displaced masses, although occupying a broadly linear assemblage, lie at various distances from their source depressions, and many individual mounds appear to relate to subsidiary basins in the major lakes; this effectively creates assemblages of closely spaced, often narrow hill-hole-pairs, collectively representing an EHHC. Some of the narrowest rafts have been fragmented along their displacement tracks to form ladder morphologies (Fig. 11). These features can be seen to merge into linear block trains that are variably streamlined to form Type 1 hogsback flutings and rubble stripes (cf. Evans et al., Reference Evans, Atkinson and Phillips2020).

Figure 11. (color online) Annotated LiDAR digital elevation model of the glacitectonic landforms in the Whitefish Lake-Goodfish Lake area, on the southern edge of the Lac la Biche fluting field. EHHC–en échelon hill-hole complex.

The overall thrust mass displacement and streamlining direction, as indicated by the various structural lineations, smoothed thrust blocks, rubble stripes, and flutings, is towards the southeast. However, the size and coherence of rafts changes rapidly within a ≤9.5-km-wide fragmentation zone comprising densely spaced, locally crenulate, and parallel ridges and hollows aligned approximately NE-SW (Fig. 11). The geometry of the ridges and intervening hollows is consistent with the fragmentation of the rafts as a result of ice flow-parallel extension occurring during transport, likely along pre-existing, ice flow-transverse joints. These relationships indicate that individual slablike rafts were detached and transported by sliding upon a subhorizontal décollement, with the relative intensity of fragmentation increasing towards the east, consistent with the increasing transport distance from the source hollows. At the far northern and eastern edge of this zone of fragmented rafts, superimposed lineaments (e.g., rubble stripes) appear to be aligned with Lac la Biche Ice Stream megaflutings (Andriashek and Fenton, Reference Andriashek and Fenton1989; Fig. 11). This indicates that the south-southeasterly flow of the Lac la Biche Ice Stream postdates the southeastern flow signature recorded in the EHHC and that intensive thrust mass fragmentation was at least partially related to the lateral drag of the Lac la Biche Ice Stream as its flow direction turned from southeasterly to south-southeasterly, effectively creating a “pull-apart” zone in the fractured substrate (Fig. 11; cf. Andriashek and Fenton, Reference Andriashek and Fenton1989). This is well illustrated by EHHC (locally developed into large-scale ladder morphology) in the area of Mann Lakes, where thrust masses have been displaced S-SE, parallel to another prominent NNW-SSE aligned scarp (see section on Ashmont-St. Paul area).

A large composite thrust moraine to the west of Whitefish Lake occurs immediately south of a complex area of bedrock hollows and fragmented rafts, located close to the southern margin of the Lac la Biche Ice Stream (Fig. 11). Its surface morphology indicates a complex internal structure, with a series of sinuous ridges representing the crests of folds and/or individual thrust-bound slices within an imbricate thrust stack. The ridge traces are crenulated by a set of open to moderately tight, NW-SE trending folds (Fig. 11). This moraine likely originated from the shallow depressions/hollows immediately to the north, and its smoothed morphology suggests that it has been subsequently overridden by ice. These characteristics indicate that the moraine is a cupola hill, constructed by ice push from the northwest/north-northwest, compatible with the orientation of the subglacial bedforms of the Lac la Biche Ice Stream to the north and east (Fig. 11). Furthermore, this cupola hill is located southwest of the Lac la Biche Ice Stream, on the northern side of an area of fragmented rafts centred upon Whitefish and Goodfish Lakes. This area was likely overlain by static or very slow-moving ice, which allowed the preservation of earlier landforms associated with raft detachment, transport, and fragmentation. The cupola hill therefore formed subglacially at a compressive margin between the Lac la Biche Ice Stream to the north and east and static ice to the south. This is unusual, in that cupola hills are typically considered to form proglacially, due to the need for space to accommodate the growing imbricate thrust stack. Consequently, the Whitefish Lake cupola hill is unlikely to have been constructed beneath thick ice. Rather, increased flow rates associated with the widening of the ice stream will have thinned the ice mass, with southerly directed thrusting effecting both the overlying ice as well as the underlying sediments and bedrock (see cross section in Fig. 11). It is possible that static ice to the south and west of the ice stream was frozen to its bed and therefore acted as an aquitard. This prohibited the pressurised subglacial meltwater, typically equated with fast ice flow, from moving beneath the ice stream into the area underlying the immobile ice. The resulting increase in subglacial water pressure on the up-ice side of the static ice mass is thought to have promoted the detachment of the thrust slices, which were transported southwards and subsequently stacked into the evolving cupola hill. Corroboration of the large-scale glacitectonic displacement of bedrock in this area is provided by Andriashek and Fenton (Reference Andriashek and Fenton1989), who reported the superimposition of more than 80 m of Belly River Formation sandstone over Empress Group sediments near Whitefish Lake. Pressurisation of the groundwater system in this area is evidenced by the widespread occurrence of prairie doughnuts and eskers immediately southwest of the composite thrust moraine (cf. Evans et al., Reference Evans, Atkinson and Phillips2020; see Fig. 11).

Ashmont-St. Paul area

Southeast of Whitefish–Goodfish Lakes is an extensive area of intensely fragmented EHHC, including numerous straight-sided, elongate, and curvilinear lakes such as Lower Mann and Upper Mann Lakes and the area west of Vincent Lake (Fig. 12; cf. Andriashek and Fenton, Reference Andriashek and Fenton1989). The western limit of these lakes is marked by a NNW-SSE escarpment, which likely marks a major fracture set or fault within the underlying bedrock (Ozoray, Reference Ozoray1972; Misra et al., Reference Misra, Slaney, Graham and Harris1991; Pana and Waters, Reference Pana and Waters2016). A series of similarly trending, highly elongate linear ridges and intervening hollows occurs in a ≤3.6-km-wide zone that parallels megaflutings along the trunk of the adjacent Lac la Biche Ice Stream (Fig. 12). This highly lineated zone forms the abrupt easterly edge to the fragmented EHHC and occupies a position relative to the ice stream that is typically associated with the development of lateral shear margin moraines (Dyke and Morris, Reference Dyke and Morris1988; Stokes and Clark, Reference Stokes and Clark2002). In the Ashmont-St. Paul area, this ice stream marginal zone comprises a complex assemblage of individual landform components, such as ladder-type morphologies or narrow zones of ribbed moraine, elongate rafts with transverse surface ridges (hereafter termed “ridged spindles”), and their more attenuated forms, here referred to as “Type 2 hogsback flutings.” We classify this style and pattern of raft evolution here as “incipient ice stream shear margin moraine,” thereby providing an explanation for the accumulation of linear subglacial debris mounds at the suture zones of ice stream beds (cf. Dyke and Morris, Reference Dyke and Morris1988; Stokes and Clark, Reference Stokes and Clark2002). This linear pattern of ice stream marginal landforms crosscuts the fragmented EHHC, indicating that in this area, the Lac la Biche Ice Stream postdated and migrated westwards over the footprint of an earlier WNW-ESE ice flow that was responsible for the detachment, transport, and fragmentation of the bedrock rafts from the elongate basins now occupied by Lower and Upper Mann Lakes.

Figure 12. Landforms in the Ashmont-St. Paul area, at the junction of the onset zone of the westerly tributary to the Lac la Biche Ice Stream: (a) annotated LiDAR digital elevation model and major glacial geomorphology of the area, with locations of enlarged areas outlined. Black lines represent the crests of major glacitectonic thrust masses; green lines depict major flutings and associated streamlined glacial bedforms; brown outlines demarcate major landform assemblages (Lottie Lake-Saddle Lake composite thrust terrain in the west and overridden and partially cannibalised megafluting in the east); and red lines demarcate discrete streamlined zones on the palaeo–ice stream beds defined by landform styles. Circled cross marks the location of section depicted in b. Detailed areas include the strike-slip raft complexes (SSRCs; incipient shear margin moraine) at western edge of the Lac la Biche palaeo–ice stream (A), hairpin-shaped landforms or murdlins amongst streamlined landforms (B), megablocks/rafts grading into linear rubble stripes (C), a downflow-tapering rubble band located east of Saddle Lake (D), and a fragmented and streamlined raft and groove (hill-groove pair) derived from a source quarry in an overridden megascale glacial lineations (MSGL) (E). (b) Section through raft in the tapering rubble band located immediately east of Lottie Lake, with upper panels showing sedimentologic and structural details and lower panel showing the nature of the rafts and streamlining on an enlargement of the LiDAR imagery. EHHC, en échelon hill-hole complex. (For interpretation of the references to color in this figure legend, the reader is referred to the web version of this article.)

The landforms in the areas to the northwest, south, and southeast of Ashmont are dominated by WNW-ESE aligned megaflutings (Fig. 12a), defining the northern part of a much larger subglacial assemblage of an ice stream onset zone that joins the Lac la Biche trunk to form Ice Stream 2/2A of Ó Cofaigh et al. (Reference Cofaigh, Evans and Smith2010) and Norris et al. (Reference Norris, Evans and Cofaigh2018). This flow set can, however, be traced into the area of fragmented EHHC centred upon the Upper and Lower Mann Lakes area, even though the smaller-scale flutings appear to have been locally subtly rotated on the surfaces of the fragmented rafts. This indicates that the Upper and Lower Mann Lakes area also forms part of the ice stream onset zone. At the centre of the onset zone, in the area south of Saddle Lake, is a large glacitectonic composite ridge (Fig. 12a). To the north and west of Saddle Lake, comparable landforms are increasingly fragmented, with individual ridges becoming more hummocky and crenulate in plan form. To the north of Lottie Lake, the ridges turn abruptly from their N-S alignment to verge with the WNW-ESE flutings south of Ashmont, where the ridges become increasingly fragmented, forming individual blocks/rafts that were evidently dragged ESE by ice stream flow to form linear rubble stripes (Fig. 12a, panel C, and 12b). This relationship, noted initially by Andriashek and Fenton (Reference Andriashek and Fenton1989), provides another example of the pull apart of glacitectonically detached blocks during transport. Similar fragmentation and downflow modification of thrust masses, detached rafts, and large plucked bedrock blocks into linear and/or downflow tapering rubble stripes are apparent throughout the onset zone, such that raft density and size can be used to differentiate ice flow parallel zones in the subglacial footprint (Fig. 12a; cf. Andriashek and Fenton, Reference Andriashek and Fenton1989). A particularly good example of a downflow tapering rubble stripe lies east of Saddle Lake and likely originated as a thrust mass that was displaced from the source depression now occupied by the lake (Fig. 12a, panel D, and 12b). In all examples of linear rubble stripes, trailing paraxial ridges are often well developed in association with megablocks/rafts, with smaller blocks often being almost completely buried in their associated prow and paraxial ridges to form hairpin-shaped landforms akin to the murdlins of Stalker (Reference Stalker1973; Figs. 1 and 12a, panel B). Raft detachment and down-ice transport and variable fragmentation also occur where pre-existing subglacial bedforms have been cannibalised, an excellent example occurring in the southern part of the onset zone. Here, three parallel-aligned and offset MSGLs recording earlier ice flow towards the SSW have been locally modified by the removal of rafts, as defined by straight-sided source quarries. A fragmented and streamlined raft (1.3-km-long and 0.5-km-wide rubble drumlinoid) derived from one of these quarries is visible 5.5 km downflow and is clearly associated with a 1.5-km-long proximal furrow indicative of its ploughing through the substrate (Fig. 12a, panel E), constituting a hill-groove pair. Linear rubble stripes, downflow tapering rubble bands (drumlinoids), ploughed rafts with paraxial ridges, and megaflutings have also been superimposed on this MSGL complex (Fig. 12a).

The rafted bedrock origins of the blocks identified throughout the onset zone are verified by roadside excavations, best illustrated by an exposure through a megablock in the tapering rubble band located immediately east of Lottie Lake (Fig. 12). Here, a vertically stacked sequence of displaced blocks of sandstone bedrock and diamicton (till), in places separated by thin thrust slices of mudstone, records displacement from the west. Diamicton with a clear vertical fissility wraps around the front of the uppermost sandstone raft to form a surface carapace. This sequence and its structural architecture indicate that the sandstone rafts have been ploughed through till, and hence the exposure represents a cross-section through a ploughing raft. As this area lies over ≤130 m of Quaternary sediment (Andriashek and Fenton, Reference Andriashek and Fenton1989), the fragmented rafts here must be relatively far-travelled erratics.

Goodridge

A ploughed raft (270 m wide and ≥20 m high) with prow and trailing paraxial ridges either side of a linear depression or groove occurs near Goodridge, south of the Beaver River (Figs. 2 and 13a). The groove represents the plough mark left behind by the displaced raft, and hence we term this association a “hill-groove pair,” similar to a PRGA. It lies amongst streamlined thrust masses and N-S orientated flutings within the footprint of the northern tributary of the onset zone of Ice Stream 2/2A of Ó Cofaigh et al. (Reference Cofaigh, Evans and Smith2010) and Norris et al. (Reference Norris, Evans and Cofaigh2018) (cf. Andriashek and Fenton, Reference Andriashek and Fenton1989). Two exposures reveal the highly deformed nature of the materials in one of the paraxial ridges and in the parent thrust mass (Fig. 13), comprising a sandy diamicton containing small lenses (intraclasts) of sand, as well as larger lenticular blocks of glacitectonised mudstone and sandstone (Type III and IV mélanges; cf. Cowan, Reference Cowan1985; Evans, Reference Evans2018). Structures include truncated recumbent folds, up-ice dipping low-angle thrusts and brittle detachments, ductile shears defining subhorizontal Y-type and down-ice dipping R-type Riedel shears, lenticular asymmetrical intraclasts, and asymmetrical extensional crenulation cleavage and SC shear fabrics (Fig. 13b), consistent with brittle-ductile shearing during the formation of this glacitectonic mélange. These mélanges are stacked within sequences comprising bedrock glacitectonite and sandy diamictic glacitectonite overlain by a more homogenous upper zone of clast-poor, sandy diamicton (Fig. 13b). The geometry of the thrusts/shears and other kinematic indicators within the highly deformed lower part of the sequence record an easterly directed sense of shear (see Fig. 13b) consistent with the lateral displacement of material from beneath the raft as it ploughed through a pre-existing glacitectonite substrate. Similar brittle-ductile shear-related structures in a quarry section farther to the north (Fig. 13b) record a more southerly directed sense of shear associated with the construction of the thrust mass, which was later dissected by the ploughing of the raft (Fig. 13a). This section also reveals that the mélange contains a number of semicoherent sandstone blocks. Together, these sediment–landform associations clearly demonstrate the linkage between the ploughing of a longitudinal groove by a gradually fragmenting sandstone raft, the base of which was displacing material laterally to form the paraxial ridges while concomitantly being brecciated and then homogenized into a sandy diamicton (cf. Atkinson et al., Reference Atkinson, Utting and Pawley2018). As the local Quaternary deposits are relatively thick in this area, the raft must represent a relatively far-travelled erratic that was locally constructing paraxial ridges in heterogeneous tills and glacitectonites that themselves were being contemporaneously generated.

Figure 13. (color online) Hill-groove pair near Goodridge: (a) annotated LiDAR digital elevation model of the glacial landforms in the area; (b) ground overview (i) of landform looking north, with exposures in the parent thrust mass (ii and iii) and in one of the paraxial ridges (iv).

Drayton Valley

Numerous exposures through Quaternary glacial deposits and associated bedrock rafts occur in a small tributary of the North Saskatchewan River at Drayton Valley (study area 2; Fig. 2). The site is in a pronounced glacitectonic landscape (Fenton et al., Reference Fenton, Trudell, Pawlowicz, Jones, Moran and Nikols1986; Tsui et al., Reference Tsui, Cruden and Thomson1989) within the footprint of the southerly flowing HPIS (Evans et al., Reference Evans, Clark and Rea2008, Reference Evans, Young and Cofaigh2014; Ó Cofaigh et al., Reference Cofaigh, Evans and Smith2010; Atkinson et al., Reference Atkinson, Utting and Pawley2014; Utting et al., Reference Utting, Atkinson, Pawley and Livingstone2016). Five main cliff sections contain important information on the processes occurring during the transport, disruption, and accretion of rafts (Figs. 14–18).

Figure 14. (color online) Sedimentologic and stratigraphic details of section 1 at Drayton Valley, including clast macrofabric data. Dmm–massive, matrix-supported diamicton; Dml–laminated, matrix-supported diamicton; Dmf–fissile, matrix-supported diamicton; MLA–mean lineation azimuth.

Figure 15. (color online) Sedimentologic and stratigraphic details of section 2 at Drayton Valley, including clast macrofabric and structural data. MLA–mean lineation azimuth.

Figure 16. (color online) Sedimentologic and stratigraphic details of section 3 at Drayton Valley, showing details of fragmented rafts of Paskapoo Sandstone and Empress Group sands and enclosing diamictons, together with clast macrofabric data and their locations.

Figure 17. (color online) Sedimentologic and stratigraphic details of section 4 at Drayton Valley, including clast macrofabric and structural data.

Figure 18. (color online) Sedimentologic and stratigraphic details of section 5 at Drayton Valley, including clast macrofabric and structural data.

Section 1 (Fig. 14) is regarded by Evans (Reference Evans2018) as an example of the vertical development from bedrock raft, to glacitectonite, to subglacial till. A raft of middle to late Palaeocene Paskapoo Formation mudstone and sandstone (Hamblin, Reference Hamblin2004; Lyster and Andriashek, Reference Lyster and Andriashek2012) here displays internal shear zone development in the form of a mélange within the upper layers of the mudstone. This is overprinted by a Dms-derived glacitectonite that has been deformed and attenuated laterally over the raft. The Dms is a coarsening-upward sequence comprising interbedded sand/silt/clay laminae, diamictons, and minor matrix-supported gravelly lenses and displays only localised soft-sediment deformation, thrusting, and dislocated recumbent folds. It appears to have originated as a subaqueous infill of a depression into which laminated sands, silts, and clays and gravity mass-flow diamictons were deposited, potentially in the low points between thrust masses before they were glacially overrun. Above this, vertical homogenisation is reflected in the zone of pseudo-laminated to fissile diamicton (Dml/Dmf) that passes upwards into Dmm with sand-filled convexo-planar lenses, correlative to the uppermost sliding/deforming bed mosaic tills in the Drayton Valley area. Clast macrofabrics from the subtill diamictons and glacitectonite shear zone (F1–F7; Fig. 14) display a W-E or WNW-ESE alignment, compatible with the WNW-ESE stress directions reflected in the early phases of glacitectonic disturbance at Drayton Valley.

In section 2 (Fig. 15), the original lithostratigraphic relationship between the pale-yellow sandstone of the Paskapoo Formation and the unconformably overlying quartzite-rich gravel of the Plio-Pleistocene Empress Group (Whitaker and Christiansen, Reference Whitaker and Christiansen1972) is preserved. Primary bedding and cross bedding within the sandstone indicate that the bedrock is the right way up. Additionally, the planar, erosive contact at the base of the Empress deposits is undeformed, indicating that the poorly cemented sandstone and unlithified quartzite-rich gravels within the raft were detached, transported, and accreted as a coherent slablike block with negligible internal disruption. The upper boundary of the raft with the overlying diamicton is marked by a narrow zone of thrusting containing highly attenuated lenses of sand and sandstone. This brittle-ductile shear zone can be traced laterally towards the east, where it forms the upper bounding surface of a second raft of Paskapoo Formation, within which the bedding is deformed by a large-scale, recumbent to gently inclined fold, the upper limb being truncated by the brittle-ductile shear zone. Bedding within the hinge of this fold has been offset by a number of low-angle brittle faults. The base of this second raft is a gently dipping fault that separates the Paskapoo Formation in its hanging wall from the structurally underlying stratified/laminated diamicton. The layering within this diamicton is defined by thin stringers of sand and is locally truncated at the thrust that forms the base of the raft. Locally this thrust is marked by a thin (≤10-cm-thick) lens of highly fragmented mudstone glacitectonite.

In section 3 (Fig. 16), the rafts are lenticular to irregular in form and internally deformed by tight to isoclinal, recumbent to very gently inclined, mesoscale disharmonic folds. These rootless folds possess highly attenuated limbs with a marked thickening of the deformed bedrock and sediments around the hinges of these complex, highly ductile structures. Although deformed, the sequence preserved within the largest raft is similar to that observed elsewhere and comprises mudstone glacitectonite (sensu Evans, Reference Evans2018; Figs. 15, 16a and b, and 18) at the base, separated from the overlying tectonised sandstone by a narrow shear zone or thrust, which is in turn overlain by quartzite-rich sandy gravel. The southwestern end of this raft is highly attenuated and thins rapidly to form a highly elongate “tail” that extends for several meters, where it is coplanar with the lamination within the host diamicton. The latter is defined by thin stringers of sand and is locally observed wrapping around the rafts. Deformation within the rafts has resulted in the variable overprinting/transposition of the primary bedding within the sandstone. The sandstone also locally contains highly irregular patches or augens of diamicton with flame-like stringers and ribbons. The boundary between the Paskapoo Formation and Empress Group within the larger raft has been strongly modified, with a similar complex, highly irregular to flame-like boundary occurring between the quartzite-rich sandy gravel and the host diamicton. These relationships indicate that, with increasing deformation, the bedrock within the raft had begun to lose its integrity and was being incorporated into the host diamicton-rich glacitectonite. This conclusion is supported by the fact that the sandstone within these highly deformed rafts is poorly indurated. The combination of the highly ductile, disharmonic nature of the folding, coupled with the complex, flame-like nature of the lithological boundaries, suggests that the raft and host sediments became water saturated during glacitectonism, which reduced cohesive strength and facilitated their breakup.

Section 4 (Fig. 17) displays a strongly modified part of the boundary between the Paskapoo Formation and the Empress Group, where highly irregular bodies of massive to weakly bedded sandstone extend upwards into the quartzite-rich gravels. Smaller (≤1-m-long) lenticular to irregular blocks of sandstone, possibly representing eroded intraclasts, occur near the base of the Empress gravels. The highly complex form of this boundary is clearly secondary in nature and suggests that the poorly lithified sandstone has been mobilised and injected upwards into the gravels. In a more deformed raft, the original stratigraphic relationship between the Paskapoo Formation and the Empress Group has been removed as a result of glacitectonism. The upper boundary of this lenticular raft (ca. 2.5 m thick and thinning abruptly towards the southeast) is a sharp, faulted contact separating the sandstone within its footwall from the structurally overlying clay-rich diamicton. The base of the raft is a complex, gently NW-dipping brittle-ductile shear zone (0.5–1.0 m thick) composed of attenuated and thrust-repeated lenses of tectonised sandstone, diamicton, and quartzite-rich gravel.

Section 5 (Fig. 18) displays thrust repetition of the Paskapoo Formation, Empress Group, and host glacigenic sediments. It is uncertain whether this raft is composed of a single thrust-repeated, W/NW-dipping slab of Paskapoo Formation and Empress Group or several, individual thrust-stacked blocks that preserve the original lithostratigraphic relationships, albeit strongly modified by deformation. Although they are repeated by thrusting, the sedimentary structures within the sandstone suggest that the bedrock blocks are the right way up. The structurally lowest, westerly dipping thrust/shear zone observed towards the eastern end of the section separates the lowest unit of sandstone from a unit of mudstone glacitectonite. The highly fissile mudstone has been broken into angular, tabular fragments (≤5–10 cm long) within a dark grey, clay-rich matrix derived from disaggregated mudstone. The moderately to gently NW-dipping thrust forming the upper boundary of the thrust stack clearly truncates bedding within the sandstone as well as the modified lithostratigraphic boundary between the displaced raft of Paskapoo Formation and Empress Group. The overlying diamicton contains lenses of sand and gravel that are wrapped by a variably developed stratification/lamination defined by thin sand stringers. The diamicton is in turn overlain by a mudstone glacitectonite, lithologically similar to that at the base of the thrust sequence. These fissile mudstones represent the overbank/lacustrine deposits within the Paskapoo Formation and, where present, will have constituted weak horizons within the bedrock that could be exploited as décollement surfaces during raft detachment, a conclusion supported by the close association between the mudstone glacitectonites and the large-scale thrusts that bound the rafts.

The diamictons that host the variably deformed rafts in sections 1–5 range from heterogeneous, pseudo-stratified, or pseudo-laminated/banded (Dms, Dml) to homogeneous, matrix-supported, and massive with zones of distinct fissility (Dmm/Dmf). The homogeneous diamictons also contain lenses of stratified sediments. The characteristics of the heterogeneous diamictons vary considerably due to the range of material that has been incorporated into them, including intraclasts, boudins, stringers. and wisps of sand, gravel, and mudstone, often appearing as densely spaced glacitectonic foliation that pinches and swells across the section and becomes more attenuated where it wraps around rafts (Figs. 14 and 16). This glacitectonic foliation resulted from the mixing and ingestion of different parent materials (bedrock, sand, gravel) within an actively deforming mélange.

The characteristics of the more matrix-dominated, locally fissile diamictons (Dmm) vary only in colour, likely reflecting a variable influence of their matrix constituents, with grey diamictons originating from local mudstones and dark brown diamictons from the Paskapoo Formation sandstone. This variation is best illustrated by section 3 (Fig. 16), where five lithofacies can be identified:

1. 1. LF 1 is a Dmm/Dml with attenuated intraclasts (boudins and stringers) of both mudstone and sand. This is separated from overlying LF 2 by a lens of interbedded and upward-coarsening sand and matrix-supported gravel with discontinuous diamicton beds. These lenses possess a planar top and convex base, resulting in a distinctive channel-like (convexo-planar) cross section;

2. 2. LF 2 is a Dmm containing large, internally heavily deformed sand and gravel intraclasts (rafts) as well as lenses similar to those at the junction with LF 1;

3. 3. LF 3 is a Dmm/Dml that grades upwards into a Dms composed of discontinuous and contorted beds that merge laterally into the attenuated limb of a large composite raft of Paskapoo Formation and Empress Group;

4. 4. LF 4 lies above the raft and like LF 3 comprises a sequence of Dmm/Dml grading upwards to Dms capped by more-fragmented rafts of Paskapoo Formation and Empress Group materials. This raft pinches out towards the southwestern end of the section, where it forms a thin, highly attenuated tail separating LF 4 and LF 5; and

5. 5. LF 5 is a Dml/Dmm in which pseudo-lamination is developed only in a thin zone directly above the raft tail and grades rapidly upwards to a fissile structure before becoming massive.

The sedimentary characteristics of these lithofacies are typical of subglacial traction tills and glacitectonites, with the channel-like form of the convexo-planar lenses of bedded sediment considered diagnostic of subglacial canal infills formed within accreting deforming layers (Eyles et al., Reference Eyles, Sladen and Gilroy1982; Clark and Walder, Reference Clark and Walder1994; Evans et al., Reference Evans, Owen and Roberts1995, Reference Evans, Phillips, Hiemstra and Auton2006; Evans, Reference Evans2018). Stress directions are evident in clast macrofabrics, with a clear WNW-ESE a-axis alignment in LF 1 (F1) being replaced by a NW-SE alignment in LF 3 (F2). This alignment is also apparent in the Dml and overlying fissile Dmm of LF 5 (F3 and F4), which changes to a southerly dip in the top of LF 5 (F5). These fabric alignments are compatible with the surface flutings in the area (cf. Atkinson et al., Reference Atkinson, Utting and Pawley2014) but importantly, clasts dip in the down-ice flow direction in the capping till (LF 5), contrary to the normal alignment of passive strain markers. As the till has been plastered over a significant bedrock raft assemblage (LFs 3 and 4), these macrofabrics likely reflect a perturbation in the stress regime due to the presence of this raft (obstacle) within the deforming bed. This is illustrated also by macrofabrics in the diamictons that cap multiple rafts in section 4 (Fig. 17). The southerly-directed stress of the last ice flow imprint at Drayton Valley is recorded also in the emplacement of a mudstone glacitectonite raft over a diamicton in section 5 (Fig. 18), in which clasts were subject to alignment by the imposed stress of the overriding raft.

Southeasterly ice flow is recorded also in clast macrofabrics from diamictons at the base and top of a sandstone raft in section 2 (Fig. 15). This is reflected in the fabric from the capping diamicton, which characterises a basal shear zone comprising thrust interdigitated Empress Group, truncated by multiple convexo-planar lenses of bedded sediments. Such characteristics are all typical of a subglacial till produced by alternating sliding-bed and deformation processes (Evans et al., Reference Evans, Phillips, Hiemstra and Auton2006; Piotrowski et al., Reference Piotrowski, Larsen, Menzies and Wysota2006; Evans, Reference Evans2018).

In summary, the Drayton Valley sections reveal a progression from: (1) internally relatively undeformed, composite rafts comprising Paskapoo Formation overlain by Empress Group gravels and sands; to (2) structurally complex rafts of bedrock and quartzite-rich gravels repeated by large-scale thrusting and folding; to c) intensely deformed and disrupted rafts. The complexity of the Quaternary sequence exposed within these sections qualifies it as a glacitectonic mélange (sensu Cowan, Reference Cowan1985; cf. Evans, Reference Evans2018). The geometry and sense of displacement of the thrusts, as well as the asymmetry of the folds and asymmetrical shape of the displaced blocks within this mélange, record a relatively consistent E/SE-directed sense of shear, compatible with the regional pattern recorded by the streamlined landforms in the area.

Warburg

A cluster of depositional crag-and-tail or stoss-and-lee megaflutings are clearly visible within the fluted and drumlinised terrain around Warburg, south of Strawberry Creek (Fig. 19a). That the crag or stoss-ends are bedrock rafts is evident in their morphology, which predominantly features straight-sided and largely flat-topped mounds associated with either relatively short, downflow paraxial ridges or single flat-topped flutings, as well as the exposure of rafted Paskapoo Formation sandstone within the quarried crags of the megaflutings (Fig. 19a, S1 and S2 sections). The megaflutings dissipate westwards (downflow), where they widen and flatten, presumably as a result of the reduced influence of the stoss rafts with distance (cf. Boulton, Reference Boulton1976; Rose, Reference Rose1989; Benn, Reference Benn1994; Evans et al., Reference Evans, Nelson and Webb2010, Reference Evans, Roberts, Hiemstra, Nye, Wright and Steer2018; Eyles et al., Reference Eyles, Boyce and Putkinen2015). Very short flutings or downflow paraxial ridges likely reflect very late stage emplacement of the rafts, explaining also their relatively modest streamlining (cf. Evans and Rea, Reference Evans, Rea and Evans2003; see Evans [2018] for modern, minor fluting analogs). Together with several broad hummocky arcs of streamlined mounds (likely thrust masses), they form the outermost 20 km of the subglacial bedform imprint of a deglacial fast-flow lobe that moved west-southwest from the Edmonton area towards the trunk zone of the High Plains Ice Stream near Drayton Valley (Atkinson et al., Reference Atkinson, Utting and Pawley2014).

Figure 19. (color online) Depositional crag-and-tail or stoss-and-lee megaflutings near Warburg: (a) annotated LiDAR digital elevation model of the glacial landforms in the area showing sections at S and N. (b) Top panel shows road cut through a less streamlined raft at location “S,” including clast macrofabric from the capping till. Bottom panels show photograph and sketch of quarry exposure in the margin of a megafluting located at “N,” including clast macrofabric from the sandy glacitectonite. Dmm–massive, matrix-supported diamicton. MLA–mean lineation azimuth.

Road cuts through two of the less streamlined rafts located to the south-southeast of Warburg reveal significant bedrock outcrops comprising subhorizontal to very gently dipping, well-bedded Paskapoo Formation sandstone and minor mudstone strata. The largest of the exposures (location S2; Fig. 19a) displays at least 3 m of bedded Paskapoo Formation that is characterised by variably coherent and incoherent beds of sandstone and mudstone. The more coherent sandstone beds represent more indurated/highly cemented layers, within which subvertical fractures represent pre-existing joints that create a blocky appearance (Fig. 19b, upper). A thin (≤1 m) unit of heavily brecciated mudstone separates the upper sandstone beds from a capping diamicton, which is massive and clay matrix supported but also contains numerous stringers and wisps of sand and is fissile, likely the product of attenuation and pulverization of sandstone intraclasts (Fig. 19b, upper). These characteristics are indicative of a subglacial traction-till origin (sensu Evans et al., Reference Evans, Phillips, Hiemstra and Auton2006; Evans, Reference Evans2018), verified by a northeasterly to east-northeasterly dipping clast macrofabric, indicative of imposed stress from that direction and compatible with the alignment of the flutings.

A quarry exposure along the margin of a megafluting located north of Warburg (location N; Fig. 19a) provides a representative stratigraphy of the processes operating in the lee side of a raft. This comprises a lower heterogeneous, sandy diamicton, containing lenses and stringers of sand and attenuated sandy intraclasts (Type III/IV mélange), separated by a thin shear zone from an overlying grey/brown, clast-poor and massive, matrix-supported diamicton with pseudo-laminated, attenuated clay stringers (Fig. 19b, lower). The shear zone at this location comprises a 0.20-m-thick fissile unit of mixed sand and silt/clay pseudo-laminae and sandy diamicton (Type IV mélange) that passes laterally to wrap around a lenticular raft of Paskapoo Formation, itself separated into two parts by a shallow thrust. A thinner raft occurs also at the base of the sequence and is separated from the overlying raft by a thrust-stacked sequence of gravel lenses and heterogeneous diamicton. The numerous thrusts, overturned folds, and fissile structures throughout the sequence, in addition to a clast macrofabric from the lower diamicton, indicate that the deposits are glacitectonites capped by till, all emplaced by ice flow from the northeast. A northeasterly dipping fabric on the flank of an ENE-WSW aligned megafluting is not unusual for deformation-generated bedforms of this type, being reminiscent of herringbone fabrics documented in minor flutings (e.g., Rose, Reference Rose1989, Reference Rose1992; Benn, Reference Benn1994; Evans et al., Reference Evans, Nelson and Webb2010). The attenuated sandstone rafts likely represent fragments of the stoss raft, transported downflow and gradually disaggregated to form bedrock glacitectonite.

Wandering River

Evidence of the initial emplacement of ploughed rafts, indicative of lodgement in subglacial deforming materials before lee-side streamlining, is well illustrated by a cluster of rafts and lee-side prows in the Wandering River area (Fig. 20). These features lie on the outer edge of the southward-orientated CAIS footprint but are associated with more subtle, superimposed WNW-ESE aligned flutings that record a late-stage/deglacial ice flow pattern. This association is manifest in the orientation of the sediment prows at the front of the megablocks (0.15–0.52 km²), which clearly show the direction in which they were pushed and crumpled based upon the multiple, arcuate ridge crests. A relatively short transport distance and lack of intensive ice overriding is evidenced by both the presence of the prows and the lack of surface smoothing of the rafts. The shallow, straight-sided depressions located immediately up-ice of the rafts are significantly (2.0–3.3 times) longer than the flow-parallel (i.e., compressed) block width and, as such, may potentially represent the source of both the bedrock raft and sediments forming the prow. The weak, ridge-like summits of the blocks could be an indication of either compression (i.e., thrusting and/or tight folding) within the raft or the expression of bedding exposed by the back rotation/tilting of the raft as it was accreted onto the up-ice side of the developing frontal prow. Consequently, these landforms constitute clones of hill-hole pairs, in that they comprise a relatively coherent raft and less coherent sediment prow lying downflow of its source depression.

Figure 20. (color online) Annotated LiDAR digital elevation model of the ploughed megablocks (rafts) near Wandering River.