Rocky Mountain Foothills Research Articles

Structural geometry and kinematic evolution of the central Canadian Rocky Mountain Foothills fold-and-thrust belt: Complex kinematic relationships controlled by detachment utilization

Abstract Balanced regional cross sections based on surface, seismic, and subsurface data show that the thin-skinned fold-and-thrust belt in the Rocky Mountain Foothills of the Kakwa area of the central Canadian Rockies consists of a lower buried thrust belt developed in Paleozoic and Triassic strata and an upper exposed faulted fold belt of Jurassic to Cretaceous strata. Changes in fold wavelength, amplitude, and geometry with stratigraphic level indicate that multiple detachments were utilized in the upper faulted fold belt. Exposed folds are chevron or box shaped. Most appear to be detachment or fault propagation folds formed by fault-to-fold displacement transfer. Geometric and kinematic relationships in the upper faulted fold belt vary from thrust faults congruently folded by underlying folds (early fault) to folds abruptly truncated by thrust faults (late fault). In contrast, folding of thrust sheets in the buried thrust belt is consistent with in-sequence deformation for all faults except one. A sequential restoration of the balanced regional cross section shows that the variable kinematic relationships observed in the upper faulted fold belt can be explained by changes in the detachment level utilized by successive faults as they climbed out of the buried thrust belt. Chevron-folded thrust faults indicate a younger fault with associated fault displacement transfer folds formed in the footwall of an older fault. These folded thrusts formed by in-sequence faulting and utilization of successively higher detachment levels. Late faults that truncate preexisting folds required out-of-sequence (hinterlandward) utilization of a higher detachment level, and they illustrate another mechanism by which critical taper is maintained in a fold-and-thrust belt.

Elastic anisotropic finite-difference full-wave modeling and imaging of 2D tilted transversely isotropic (TTI) media

Imaging below the tilted transversely isotropic (TTI) media pose serious problems and distortions of the subsurface geological targets of interest for hydrocarbon exploration, which is very common in fold-thrust belts or active tectonic areas like the Rocky Mountain foothills of the Canada or foothills of the Himalaya. To model and image below the TTI media having symmetry axis orthogonal to dipping anisotropic layers, a 2D TTI thrust sheet is considered for elastic anisotropic finite-difference full-wave modeling and generation of synthetic seismic data using the robust staggered-grid scheme of wave propagation. Free-surface boundary conditions at the top and absorbing boundary conditions for the other three sides of the model have been imposed to reduce undesirable edge effects and suppress dispersions. The synthetic seismic data generated for the 2D TTI thrust model are migrated using the anisotropic Kirchhoff pre-stack time (PSTM) and depth migration (PSDM) technique followed by the migration velocity analysis (MVA) algorithm. The traveltime computations are made using robust eikonal equations with fourth-order elastic anisotropic finite-difference solutions to handle large tilts (ν) along with anisotropic ray-tracing to image below the TTI thrust sheet. The MVA algorithm adopted also estimates the anisotropic parameters ε and δ accurately within the acceptable error limit through successive iterations of inversion and parameter update, which results in sufficient flattening of the reflection events in the common image gathers (CIGs). This is a powerful and stable methodology to handle large tilt (ν = 60°) of the steeply dipping reflectors in complex geological structure like the TTI thrust model without much distortion of CIGs for which the symmetry axis of each TTI block is set orthogonal to the bottom of the reflectors having different tilts.

A seismic structural overview of Liard Basin, Northeast British Columbia, Canada

The Liard Basin is an important subbasin of the Western Canada Sedimentary Basin located in Northeast British Columbia along the eastern margin of the Canadian Cordillera. It contains significant potential unconventional gas resources but is largely underrepresented in public literature. Using available-for-purchase 2D seismic data, a regional structural interpretation of the basin was completed providing the first seismically controlled, high-level overview of the structural features of the basin and its surrounding area. The shape of the Liard Basin largely reflects the orientation of older Paleozoic and Proterozoic extensional structures that localized structures formed during Cretaceous-Tertiary compressive deformation. The eastern boundary of the basin is marked by the well-documented Bovie Structure. The Liard Anticline and the Liard River Anticline found near 60°N latitude are the only large structures located within the Liard Basin proper. Inversion of the herein named Liard Basin Boundary Structure, a west-side-down fault zone of Early Paleozoic age, localized the northwest boundary of the basin with the Liard Fold and Thrust Belt. A triangle zone bounds the Rocky Mountain Foothills and the Liard Basin to the southwest. Reflectors in the Proterozoic strata below the Liard Basin were deformed by compressive and then extensional structures prior to the deposition of Paleozoic strata. Proterozoic strata are involved in all of the major structures of the adjacent Liard Fold and Thrust Belt, the Rocky Mountain Foothills, and the Bovie Structure. These structures controlled the location of the Liard Basin.

Marinus Pilon.

Marinus Pilon.

Beryllium-10 dating of the Foothills Erratics Train in Alberta, Canada, indicates detachment of the Laurentide Ice Sheet from the Rocky Mountains at ~15 ka

AbstractThe Foothills Erratics Train consists of large quartzite blocks of Rocky Mountains origin deposited on the eastern slopes of the Rocky Mountain Foothills in Alberta between ~53.5°N and 49°N. The blocks were deposited in their present locations when the western margin of the Laurentide Ice Sheet (LIS) detached from the local ice masses of the Rocky Mountains, which initiated the opening of the southern end of the ice-free corridor between the Cordilleran Ice Sheet and the LIS. We use 10Be exposure dating to constrain the beginning of this decoupling. Based on a group of 12 samples well-clustered in time, we date the detachment of the western LIS margin from the Rocky Mountain front to ~14.9 ± 0.9 ka. This is ~1000 years later than previously assumed, but a lack of a latitudinal trend in the ages over a distance of ~500 km is consistent with the rapid opening of a long wedge of unglaciated terrain portrayed in existing ice-retreat reconstructions. A later separation of the western LIS margin from the mountain front implies higher ice margin–retreat rates in order to meet the Younger Dryas ice margin position near the boundary of the Canadian Shield ~2000 years later.

Terrestrial Bryophyte and Lichen Responses to Canopy Opening in Pine-Moss-Lichen Forests

Pinus contorta-dominated montane forests of western Canada with relatively dense tree canopies have ground layers with abundant bryophytes, especially the feather mosses (Pleurozium schreberi and Hylocomium splendens), while those with more open canopies are dominated by species of reindeer lichens, especially Cladonia arbuscula s.l. and C. rangiferina s.l. Woodland caribou (Rangifer tarandus caribou), which are a threatened species in Alberta, prefer open, Cladonia-dominated forests for their winter food supply. This study investigated if opening the canopy by thinning mature montane forests of the Canadian Rocky Mountain foothills would change the abundance of lichens and bryophytes. In 1997, forests were thinned by removing 20%, 40%, and 60% by volume. In 2016, 19 years after treatment, we re-surveyed a subset of these plots (n = 97) for lichen and bryophyte abundance and species richness by utilizing the amount of canopy opening at the plot level as our prime gradient. We then used ordination to determine the relationship of control plots to treatment plots. In uncut forest, the control plots were highly variable, but were mostly dominated by feather mosses, with little or no bare ground. Feather moss abundance was lower in treatment plots when compared to control plots, while cover of bare ground was greater. Overall, 19 years after treatment, we found that, in treatment plots, lichen abundance remained stable or slightly increased, feather mosses decreased markedly, and unoccupied space was double that of the control plots. We conclude that the canopy opening had little effect on understory and ground layer diversity, but considering species abundance (1) bryophytes have not recovered after canopy opening, (2) populations of reindeer lichens increased marginally, but have not colonized areas left bare from bryophyte dieback, and (3), after 19 years there, remains unoccupied areas of bare ground in plots with a reduced canopy cover. Our study demonstrated that, with canopy cover reduction resulting from forest thinning operations, the ground layer diversity is maintained, but recovery of ground layers in old-growth pine-dominated forests is not promoted. Therefore, timber harvest that partially opens the tree canopy is unlikely to benefit caribou by augmenting or accelerating winter food availability and habitat suitability for caribou.

Shear Wave Splitting Discloses Two Episodes of Collision‐Related Convergence in Western North America

AbstractSeismic anisotropy imposes first‐order constraints on the strain history of crust and upper mantle rocks. In this study, we analyze the mantle seismic anisotropy of the Western Canada Sedimentary Basin using a new shear wave spitting data set consisting of 1,333 teleseismic arrivals from 82 seismic stations. The resulting 332 high‐quality measurements yield a regional mean apparent splitting time (i.e., the magnitude of anisotropy) of 1.1 ± 0.3 s and an average fast orientation (i.e., the direction of anisotropy) of 54.6° ± 17.2°, which favor a two‐layer anisotropic model based on the 90° back azimuthal periodicity in both parameters. The northeast trending fast orientations dominate the lower layer at lithospheric depths and are approximately parallel to the present‐day absolute plate motions (APMs; i.e., <35°) due to the active asthenospheric flow. On the other hand, deviations from the APMs along the Canadian Rocky Mountain foothills could reflect disrupted mantle flow surrounding a southwestward migrating cratonic lithosphere. Also revealed are two elongated upper‐layer anisotropic anomalies in the lithosphere that are spatially correlated with Moho depths. Their characteristics suggest frozen‐in anisotropy imprinted along two convergent boundaries: (1) the Paleoproterozoic Snowbird Tectonic Zone that separates northeast (north) from northwest (south) fast directions and (2) the foothills of the Rocky Mountains that exhibit northeast trending orientations consistent with those of the APMs, maximum crustal stress, and electromagnetic anisotropy. Compressions associated with the Cordilleran orogenesis could be responsible for the spatial changes in the shear wave anisotropy from the foothills to the cratonic interior.

Historical land use and land cover for assessing the northern Colorado Front Range urban landscape

ABSTRACT We describe historical land-use and land-cover (LULC) maps for the northern Colorado urban Front Range. The Front Range urban landscape is diverse and interspersed with highly productive agriculture, as well as natural land cover types including evergreen forest in the Rocky Mountain foothills and Great Plains grassland. To understand the dynamics of urban growth, raster maps were created at a 1-meter resolution for each of four time steps, nominally 1937, 1957, 1977, and 1997. In total, 38 detailed LULC classes were identified using manual interpretation techniques, aerial photographs, historical maps, and other available information. The maps provide high-resolution spatial data for understanding the historical progression of urbanization and will allow further analysis of the effects of urban growth on social and ecological systems.

Transplanting Following Non-Native Plant Control in Rocky Mountain Foothills Fescue Grassland Restoration

Disturbed areas within national parks, resulting from historic and current land use activities, can harbor large and diverse populations of non-native plant species. These species must be controlled to prevent their spread and to restore native grassland. Research on alternative revegetation methods is urgently needed since seed based grassland restoration is often unsuccessful. The effectiveness of transplanting to restore foothills fescue grassland following implementation of non-native plant management was investigated at three disturbed sites in Waterton Lakes National Park, Alberta, Canada. Non-native plant abundance was reduced by cutting, glyphosate application, and steaming. Greenhouse grown seedlings of native fescue grassland species (13 forbs, five grasses, two shrubs) were outplanted into plots following management treatments, as three different sizes of container stock. The effects of prior non-native plant management and container size on transplant survival and growth were assessed over two growing seasons. Rhizomatous forbs and bunchgrasses with well-developed root systems had highest survival ( 50% over 14 months). Transplants with 10- to 15-cm rooting depth in cones and root trainers had significantly higher within-year survival than tray transplants with a shallow rooting depth. Transplant survival was improved by glyphosate application to control non-native plants prior to planting. Transplanting was effective for increasing native cover and species richness although high winter mortality reduced this effect. Key species for fescue grassland function were introduced at all sites and persisted for two growing seasons.

Controls on accommodation during the early-middle Campanian in southern Alberta, western Canada foreland system

The Western Canada Sedimentary Basin is a retroarc foreland system produced by thrust sheet loading in the Canadian Cordillera and resultant flexural displacement of the lithosphere east of the orogenic front. Theoretical models of flexural tectonism and dynamic loading have long been used to predict the effects of these controls on the creation/destruction of accommodation in retroarc foreland settings, however, few field studies have attempted to evaluate their relative contributions for a given point in geologic time. As key regional drivers of subsidence/uplift and implicitly, relative sea level rise and fall, flexural tectonism and dynamic loading impart unique signatures (subsidence geometries) which can be discerned via the spatial and temporal analysis of event-significant surfaces preserved in the geologic record. We use magnetostratigraphy to constrain the timing of a major marine incursion (transgressive phase of the second-order Campanian Claggett Cycle, which spanned ∼8 Myr in its entirety) at three locations in the southern Alberta plains where the onset of this transgression is marked by a widespread ravinement surface. A chronostratigraphic framework is presented which includes new empirical ages for stratigraphic units above and below this key surface, informally termed the Milk River ‘shoulder’. The framework provides a basis for determining the earliest extents of the Claggett (Pakowki/Lea Park) sea in southern Alberta, and for characterization of its subsequent expansion during the early-middle Campanian. Although commonly depicted as a simple, westward-expanding seaway, we observe that the initial expansion of the Claggett sea in southern Alberta was more complicated and proceeded in multiple directions through time, prior to maximum transgression and the ensuing second-order regression. We interpret the Claggett sea as restricted to the Alberta foredeep at ∼82.0 Ma, before inundating a paleotopographic high associated with a prominent basement feature in southernmost Alberta (Sweetgrass Arch) at ∼81.0 Ma, likely establishing a connection with the Claggett sea in northern Montana. As relative sea level continued to rise, the sea reached the Rocky Mountain foothills in the west while also advancing more than 100 km eastward into the forebulge region during peak transgression (∼80.0 Ma). With this history established, we ascertain the overriding cause for the Claggett transgression was dynamic subsidence because (1) the magnitude, direction, and timing of relative sea level rise and the associated expansion of the Claggett sea fits published dynamic subsidence models for the Campanian of North America, and (2) several hundred metres of marine sediments are preserved on the early-middle Campanian forebulge, requiring accommodation which cannot be explained by widely accepted flexural models or eustatic estimates.

Structural interpretation and restoration of Rocky Mountain Brazeau zone

The Rocky Mountain Foothills lie along the eastern margin of the Rocky Mountain fold-thrust belt. The area has been the focus of extensive research aimed at locating oil and gas fields with the potential to be used as CO2 storage traps. In this study, we use a seismic line from the Canadian Rockies to interpret the geologic structures along a cross-section parallel to the tectonic transport direction. We then compare our results with those of previous studies. The section was restored using the MOVE software (manufactured by Midland Valley Exploration Ltd.). The primary objectives of this work are: (1) to conduct a stratigraphic and structural interpretation of a 2D seismic profile; and (2) to conduct a cross-sectional restoration of the structures in order to validate the seismic interpretation in terms of CO2 storage candidates. Additional data sources include maps of the surface geology, which show that the age of horizons decrease from west to east, and stratigraphic and structural profiles derived from well logs. The results of our structural restoration indicate a detachment fault between the foreland and hinterland. This fault is responsible for the cutting and subsequent upwards and eastwards movement of a stratum located between the basement and the Late Devonian formation. Large thrust faults are responsible for the deformation of strata (through both folding and faulting) in the foreland basin. As a result of continuous eastward tectonic stress, the strata from Jurassic have deformed, forming a duplex system in the middle of the section and resulting in the uplift of the upper part of the section. Following surface erosion, this uplifted area became exposed during the Tertiary Period. The high shortening rate (53%) detected through structural restoration is consistent with the thin-skinned tectonic model.

Atmospheric research in the Rocky Mountain foothills

As an NSF-funded center approaches its 60th birthday, scientific and fiscal challenges lie ahead.

Chapter 1: Integrated, High-Resolution Allostratigraphic, Biostratigraphic and Carbon-Isotope Correlation of Coniacian Strata (Upper Cretaceous), Western Alberta and Northern Montana

ABSTRACT Lower to upper Coniacian rocks in the foredeep of the Western Canada Foreland Basin are dominated by mudstone and subordinate sandstone and were deposited on a very low-gradient, storm-dominated marine ramp. The rocks are organized into several scales of upward-coarsening, upward-shoaling succession, bounded by marine flooding surfaces. Abundant, publicly available wireline log data permit flooding surfaces to be traced for hundreds of kilometers in subsurface. Flooding surfaces can be considered to approximate time surfaces that allow the subsidence history of the basin to be reconstructed. Particularly widely traceable flooding surfaces were chosen, on pragmatic grounds, as the boundaries of 24 informal allomembers, most of which can be mapped along the foredeep for >750 km. Allomembers can also be traced westward into the fold-and-thrust belt to outcrop in the Rocky Mountain Foothills. Some flooding surfaces are mantled with intra- or extrabasinal pebbles that imply a phase of shallowing and, ...

Spatial Heterogeneity of the Forest Canopy Scales with the Heterogeneity of an Understory Shrub Based on Fractal Analysis

Spatial heterogeneity of vegetation is an important landscape characteristic, but is difficult to assess due to scale-dependence. Here we examine how spatial patterns in the forest canopy affect those of understory plants, using the shrub Canada buffaloberry (Shepherdia canadensis (L.) Nutt.) as a focal species. Evergreen and deciduous forest canopy and buffaloberry shrub presence were measured with line-intercept sampling along ten 2-km transects in the Rocky Mountain foothills of west-central Alberta, Canada. Relationships between overstory canopy and understory buffaloberry presence were assessed for scales ranging from 2 m to 502 m. Fractal dimensions of both canopy and buffaloberry were estimated and then related using box-counting methods to evaluate spatial heterogeneity based on patch distribution and abundance. Effects of canopy presence on buffaloberry were scale-dependent, with shrub presence negatively related to evergreen canopy cover and positively related to deciduous cover. The effect of evergreen canopy was significant at a local scale between 2 m and 42 m, while that of deciduous canopy was significant at a meso-scale between 150 m and 358 m. Fractal analysis indicated that buffaloberry heterogeneity positively scaled with evergreen canopy heterogeneity, but was unrelated to that of deciduous canopy. This study demonstrates that evergreen canopy cover is a determinant of buffaloberry heterogeneity, highlighting the importance of spatial scale and canopy composition in understanding canopy-understory relationships.

Chemical, Physical, and Biological Factors Shape Littoral Invertebrate Community Structure in Coal-Mining End-Pit Lakes.

Aquatic invertebrates form the base of the consumer food web in lakes. In coal-mining end-pit lakes, invertebrates are exposed to an environment with potentially challenging physical and chemical features. We hypothesized that the physical and chemical features of end-pit lakes reduce critical littoral habitat and thus reduce invertebrate diversity, thereby limiting the potential for these lakes to be naturalized. We used a multivariate approach using principle component analysis and redundancy analysis to study relationships between invertebrate community structure, habitat features, and water quality in five end-pit lakes and five natural lakes in the Rocky Mountain foothills of west-central Alberta, Canada. Results show a significantly different invertebrate community structure was present in end-pit lakes as compared with reference lakes in the same region, which could be accounted for by water hardness, conductivity, slope of the littoral zone, and phosphorus concentrations. Habitat diversity in end-pit lakes was also limited, cover provided by macrophytes was scarce, and basin slopes were significantly steeper in pit lakes. Although water chemistry is currently the strongest influencing factor on the invertebrate community, physical challenges of habitat homogeneity and steep slopes in the littoral zones were identified as major drivers of invertebrate community structure. The addition of floating wetlands to the littoral zone of existing pit lakes can add habitat complexity without the need for large-scale alterations to basing morphology, while impermeable capping of waste-rock and the inclusion of littoral habitat in the planning process of new pit lakes can improve the success of integrating new pit lakes into the landscape.

Forms and drivers of annual streamflow variability in the headwaters of Canadian Prairies during the 20th century

AbstractHeadwater streamflows in the Rocky Mountain foothills are the key to water availability in the Canadian Prairies. Headwater characteristics, however, have been and continue to be subject to major variability and change. Here, we identify various forms of change in the annual mean streamflow and timing of the annual peak and attempt to distinguish between the effects of multiple drivers using a generalized regression scheme. Our investigation shows that the Pacific Decadal Oscillation (PDO) is the main driver of significant monotonic trends and shifts in the central tendency of annual mean streamflow in major headwaters. In parallel, the cumulative effects of non‐PDO climatic drivers and human‐induced land use and land management are the main causes of significant variations in the timing of the annual peak. Additional analyses show that time sequences with significant trends in annual mean streamflow and timing of the annual peak coincide with those that show significant trends in the PDO or non‐PDO component of the air temperature, respectively. The natural streamflow characteristics are substantially perturbed by anthropogenic river flow regulation, depending on the form of change and/or the level of regulation. Evidence suggests that the general tendency of human regulation is to alleviate the severity of above‐ and below‐average streamflow conditions; however, it may also intensify the variability in natural streamflow characteristics during drier years and/or those with earlier annual peak timing. These are circumstances to which the regional water resource system is vulnerable. Our findings are important for the provision of effective regional water resource management in the Canadian Prairies and contribute to a better understanding of the complex interactions between natural and anthropogenic drivers in coupled human–water systems.

Reconstructing the confluence zone between Laurentide and Cordilleran ice sheets along the Rocky Mountain Foothills, south‐west Alberta

ABSTRACTLandform mapping near the Athabasca valley was fundamental to determining whether the Late Wisconsinan Laurentide and Cordilleran Ice Sheets (LIS, CIS) coalesced. In this paper we examine the detailed landform record using new LiDAR digital elevation and hillshade models of the area of confluence and eventual desuturing of these ice sheets. This work reveals an earlier more extensive Cordilleran advance before coalescence with the LIS. When the ice sheets coalesced, the flow pattern was dominated by ice flow along the mountain front, with Cordilleran ice flowing out of major trunk valleys but Laurentide ice flowing into the Foothills up smaller valleys. This flow pattern implies that when the ice sheets coalesced the CIS was already waning, or at least was not at its maximum. Deglaciation was interrupted by several re‐advances, probably related to the destabilizing effect of proglacial lakes along the ice margin.

Fracture abundance and strain in folded Cardium Formation, Red Deer River anticline, Alberta Foothills, Canada

The folded and thrusted Mesozoic clastic sequence of the Canadian Rocky Mountain foothills forms important hydrocarbon reservoirs; however, the generation of fractures and their timing relative to folding and thrusting have remained unclear. Understanding the distribution and evolution of natural fractures, especially their timing of formation relative to the evolution of the fold-and-thrust system, could potentially improve exploration and development outcomes in these otherwise tight, low-permeability reservoirs. We investigated the relation between folding and fracture formation in the Upper Cretaceous Cardium Formation by combining field structural observations and kinematic modeling of the fold-and-thrust belt evolution. We evaluated the effects of structural position and strain on the formation of predominantly opening-mode fractures by analyzing fracture intensity and aperture in the backlimb, the shallow forelimb close to the crest, and the steep forelimb, away from the crest of a reservoir-scale anticline, the Red Deer River anticline. Fracture strain determined from both outcrop and microscanlines varies little across the three structural domains of the fold, although somewhat lower fracture strain was measured in the forelimb close to the crest. These fracture strain measurements are qualitatively consistent with calculated horizontal strain in the tectonic-transport direction obtained from kinematic numerical models that simulate fold development associated with slip along the underlying Burnt Timber thrust. The kinematic models predict similar magnitudes of horizontal extension in both the back and forelimbs, and somewhat lower extension in the upper forelimb during early development of the Red Deer River anticline. Predicted early fracture formation during fold development is consistent with field structural observations of shear reactivation during later stages of folding. This combined kinematic modeling and field structural study demonstrates that deforming fold-and-thrust belts can undergo a complex evolution of bed-parallel extension in both space and time, resulting in spatially variable fracture formation in structurally complex subsurface reservoirs.

Record of orogenic cyclicity in the Alberta foreland basin, Canadian Cordillera

Jurassic–Cretaceous sedimentary rocks of the Alberta foreland basin are a key record of the evolution of the Canadian Cordillera. We test a recent model for cyclical development of Cordilleran orogenic systems using detrital zircon analysis of the major sandstone units deposited between 145 and 80 Ma exposed in the Rocky Mountain Foothills near Grande Cache, Alberta. The basin history is well constrained by decades of study, and the stratigraphy has been previously subdivided into tectonostratigraphic wedges. U-Pb data from 14 detrital zircon samples are included in this study. All the major magmatic provinces of North America are represented in each sample, with the relative proportions varying between samples. The samples are assigned to five groups with the aid of multidimensional scaling. Groups 1–3 are interpreted to record recycling from specific passive-margin units of western North America with varying input from the Cordilleran magmatic arc. Group 4 is interpreted to record recycling from sedimentary strata in the United States and dispersal by basin-axial fluvial systems. Group 5 is dominated by Mesozoic zircon grains interpreted to have originated in the Cordilleran magmatic arc. Detrital zircon age spectra do not form groups based on the tectonostratigraphic wedges from which they were sampled; rather, within each tectonostratigraphic wedge, they exhibit evolution from diverse age spectra to a less-diverse distribution of detrital zircon ages.

Coalescence of late Wisconsinan Cordilleran and Laurentide ice sheets east of the Rocky Mountain Foothills in the Dawson Creek region, northeast British Columbia, Canada

Geomorphic, stratigraphic and geochronological evidence from northeast British Columbia (Canada) indicates that, during the late Wisconsinan (approximately equivalent to marine oxygen isotope stage [MIS] 2), a major lobe of western-sourced ice coalesced with the northeastern-sourced Laurentide Ice Sheet (LIS). High-resolution digital elevation models reveal a continuous 75 km-long field of streamlined landforms that indicate the ice flow direction of a major northeast-flowing lobe of the Cordilleran Ice Sheet (CIS) or a montane glacier (>200 km wide) was deflected to a north-northwest trajectory as it coalesced with the retreating LIS. The streamlined landforms are composed of till containing clasts of eastern provenance that imply that the LIS reached its maximum extent before the western-sourced ice flow crossed the area. Since the LIS only reached this region in the late Wisconsinan, the CIS/montane ice responsible for the streamlined landforms must have occupied the area after the LIS withdrew. Stratigraphy from the Murray and Pine river valleys supports a late Wisconsinan age for the surface landforms and records two glacial events separated by a non-glacial interval that was dated to be of middle Wisconsinan (MIS 3) age.