This doctoral dissertation presents a study of the glacial history of the North American Cordillera using numerical ice sheet modelling calibrated against field evidence. This area, characterized by the steep topography of several mountain ranges separated by large inter-montane depressions, was once covered by a large-scale ice mass: the former Cordilleran ice sheet. Because of the irregular topography on which the ice sheet formed, geological studies have often had only local or regional relevance, thus leaving the Cordilleran ice sheet least understood among Pleistocene ice sheets in terms of its extent, volume, and dynamics. Here, I present numerical simulations that allow quantitative reconstructions of the former ice sheet evolution based on approximated physics of glacier flow. These simulations show that the geometry of the Last Glacial Maximum Cordilleran ice sheet was largely controlled by sharp contrasts in regional temperature, precipitation, and daily temperature variability associated with the presence of mountain ranges. However, this maximum stage appears short-lived and out of balance with contemporaneous climate. During most of the simulated last glacial cycle, the North American Cordillera is characterized by an intermediate state of glaciation including isolated glaciers and ice caps covering major mountain ranges, the largest of which is located over the Skeena Mountains. The numerically modelled Cordilleran ice sheet appears in constant imbalance with evolving climate conditions, while the complexity of this transient response transcends that encapsulated in two-dimensional, conceptual models of ice sheet growth and decay. This thesis demonstrates the potential of numerical ice sheet modelling to inform on ice sheet history and former climate conditions over a glacial cycle, given that ice sheet models can be calibrated against field constraints. / <p>At the time of the doctoral defense, the following paper was unpublished and had a status as follows: Paper 4: Mansucript.</p>

Glaciers and ice sheets covered extensive areas in the Northern Hemisphere during the last glacial period. Subsequently to the Last Glacial Maximum (LGM), they retreated rapidly and, except for Greenland and some other ice caps and glaciers, they vanished after the last glacial termination. This thesis examines the dynamics of deglacial environments by analysing the glacial geomorphological record with focus on the landforms created by glacial meltwater. The aims are (i) to evaluate the data available for mapping glacial meltwater features at the regional scale, and (ii) to demonstrate the potential of such features for regional ice retreat reconstructions in high-relief landscapes. Meltwater landforms such as ice-marginal meltwater channels, eskers, deltas and fossil glacial lake shorelines are used to infer former ice surface slope directions and successive positions of retreating ice margins. Evaluated high-resolution satellite imagery and digital elevation models reveal their potential to replace aerial photographs as the primary data for mapping glacial meltwater landforms. Following a methods study, reconstructions of the deglacial dynamics are carried out for central Transbaikalia, Siberia, Russia, and for the Cordilleran Ice Sheet (CIS) in central British Columbia, Canada, using regional geomorphological mapping surveys. Mapped glacial landforms in central Transbaikalia show evidence of a significant glaciation that possibly extended beyond the high mountain areas. Large glacial lakes were formed as advancing glaciers blocked rivers, and of these, Glacial Lake Vitim was the most prominent. Deglacial dynamics of the CIS reveals that the ice divide shifted to the Coast Mountains in north-central British Columbia and the eastern ice margin retreated towards the ice divide in late glacial time. This thesis demonstrates the potential to reconstruct ice retreat patterns and deglacial dynamics at regional scales by interpretation of the meltwater landform record. / At the time of the doctoral defense, the following papers were unpublished and had a status as follows: Paper 5: Submitted. Paper 6: Manuscript.

This dissertation presents results from a multi-year interdisciplinary study of the Late Quaternary geomorphic history of northwest Calvert Island, British Columbia, Canada. There is a considerable knowledge gap in the region pertaining to Cordilleran ice cover and extent as well as landscape response to a uniquely stable relative sea-level history. The objective of this study was to reconstruct this regional landscape response to deglaciation including post-LGM ice cover and extent, relative sea-level changes, coastal landform development, and climate and ecological variance. Methods used to inform this reconstruction included airborne lidar, aerial photography interpretation, sedimentary stratigraphy and detailed sedimentology of samples from shovel pits and lake cores, surficial geology and geomorphic mapping, palaeoecological examinations, and the development of a geochronology using radiocarbon and optical dating. To assist with landscape reconstruction, a new method was developed and used to differentiate littoral and aeolian sands in sediment samples that range in age from Mid to Late Holocene by using modern reference samples. The method utilized a standard optical microscope paired with freely available software (ImageJ) to characterize grain shape parameters. The method was tested on nearly 6,000 sand grains from samples of known and hypothesized depositional settings and was able to correctly identify the depositional setting for 76% of the samples. After testing, the method was used to differentiate littoral and aeolian sands in a number of shovel pit, exposure, and core sediment samples to give context to stratigraphic and geomorphic interpretations. A short-lived Late Pleistocene re-advance of Cordilleran ice occurred in the study area, with radiocarbon ages indicating ice advanced to, and then retreated from, the western edge of Calvert Island between 14.2 and 13.8 ka cal BP, respectively. Sedimentological and palaeoecological information that suggests a cold climate and advancing/retreating glacier as well as lidar remote sensing and field-based geomorphic mapping of moraines in the region provide evidence of the re-advance. After ice retreated from the area, a broad suite of geomorphic landforms developed, including flood plains, iv aeolian dunes, beaches, spits, marshes, and tombolos. Coastal reworking was extensive, with progradation rates greater than 1 m a-1 occurring in some locations during the Late Holocene. These data provide the first evidence of a re-advance of the retreating ice sheet margin on the central coast of British Columbia, contribute an important methodology to advance Quaternary reconstructions, and give a unique account of the geomorphic development of a Pacific Northwest coastline that experienced little relative sea-level change over the Late Pleistocene and Holocene. Results help fill a spatial and temporal gap in the landscape history of British Columbia and have implications for climate and sea-level reconstructions, early human migration patterns, and the palaeoenvironment of an understudied area of the Pacific Northwest coast of North America. / Graduate / 0368 / 0372 / 0426 / jeamer7@gmail.com

New relative sea-level (RSL) observations dating from the late Pleistocene and early Holocene, during and after the collapse of the Cordilleran ice-sheet (CIS), are provided for two regions in southern coastal British Columbia. They record the glacial isostatic adjustment (GIA) response of the Earth to the changing surface load of the waning CIS. The data provide a new RSL curve for Sechelt, on the mainland coast north of Vancouver, and extend and revise a previously constructed curve for Barkley Sound on the west coast of Vancouver Island. The observations create a new profile of RSL curves oriented southwest-northeast across Vancouver Island and the Strait of Georgia. A previously-defined profile of RSL curves is oriented northwest-southeast profile along the east coast of Vancouver Island. The two profiles intersect in the central Strait of Georgia. The new RSL curves sample different parts of the Cascadia Subduction Zone (CSZ) and provide constraints on the history of the CIS. The Juan de Fuca plate subducts beneath the North American plate in roughly the same southwest to northeast direction as the RSL profile. GIA modelling of the RSL observations along this profile may indicate spatial variations related to the structure of the Cascadia Subduction Zone (CSZ). The CIS flowed roughly from northeast to southwest over the regions of interest. RSL observations along this path indicate how sea-level change differed with distance from the edge of the ice-sheet towards its centre. The CIS model of James et al. (2009b) is refined to fit observed sea levels while applying glacial geological constraints to regional ice sheet advance and retreat. Sea level in Barkley Sound dropped from greater than 27 m elevation before 15 cal kyr BP to -46 m below present around 12 cal kyr BP. At Sechelt, sea level closely follows the same trend as in the central Strait of Georgia, dropping from over 150 m before 14 cal kyr BP and falling past present levels after 12.4 cal kyr BP to a poorly constrained lowstand between 12 and 9 cal kyr BP. The initial crustal uplift rate near Sechelt was at least 85 mm/yr, comparable to that of the central Strait of Georgia. The sea-level observations are best fit with predictions employing an Earth model with a 60-km effective lithosphere thickness and asthenospheric viscosity and thickness of 4 × 1019 Pa s and 380 km, respectively. The transition zone and lower mantle viscosities are based on the VM2 Earth model (Peltier 2002). Sea level in Barkley Sound fell quickly (15-30 mm/yr), and observed sea level is best fit with the same asthenospheric viscosity, but with a thinner 30-km thick lithosphere, consistent with the regional tectonic structure. Revisions to the ice model are consistent with radiocarbon constraints on ice sheet history and provide good agreement with the observed sea-level history for the study regions as well as RSL histories previously described for the Strait of Georgia and southern Vancouver Island. / Graduate / 0372

In the late Pleistocene, most of British Columbia and northern Washington was covered by the Cordilleran ice sheet. The weight of the ice sheet caused up to several hundred metres of depression of the Earth’s crust. This caused relative sea level to be higher in southwestern British Columbia despite lower global eustatic sea level. After deglaciation, postglacial rebound of the crust caused sea level to quickly drop to below present levels. The rate of sea-level fall is used here to determine the rheology of the mantle in southwestern British Columbia. The first section of this study deals with determination of the postglacial sea-level history in the Victoria area. Constraints on sea-level position come from isolation basin cores collected in 2000 and 2001, as well as from previously published data from the past 45 years. The position of sea-level is well constrained at elevations greater than -4 m, and there are only loose constraints below that. The highstand position in the Victoria area is between 75-80 m. Sea level fell rapidly from the highstand position to below 0 m between 14.3 and 13.2 thousand calendar years before present (cal kyr BP). The magnitude of the lowstand position was between -11 and -40 m. Though there are few constraints on the lowstand position, analysis of the crustal response favours larger lowstand. Well constrained sea-level histories from Victoria, central Strait of Georgia and northern Strait of Georgia are used to model the rheology of the mantle in southwestern British Columbia. A new ice sheet model for the southwestern Cordillera was developed as older models systematically underpredicted the magnitude of sea level in late glacial times. Radiocarbon dates are compiled to provide constraints on ice sheet advance and retreat. The Cordillera ice sheet reached maximum extent between 17 and 15.4 cal kyr BP. After 15.4 cal kyr, the ice sheet retreated, and by 13.7 cal kyr BP Puget Sound, Juan de Fuca Strait and Strait of Georgia were ice free. By 10.7 cal kyr BP, ice was restricted to mountain glaciers at levels similar to present. With the new ice model, and using an Earth model with a 60 km lithosphere, asthenosphere with variable viscosity and thickness, and transitional and lower mantle viscosity based on the VM2 Earth model, predicted sea level matches the observed sea level constraints in southwestern British Columbia. Nearly identical predicted sea-level curves are found using asthenosphere thicknesses between 140-380 km with viscosity values between 3x10^18 and 4x10^19 Pa s. Predicted sea level is almost completely insensitive to the mantle below the asthenosphere. Modeled present day postglacial uplift rates are less than 0.5 mm yr^-1. Despite the tight fit of the predicted sea level to observed late-glacial sea level observations, the modelling was not able to fit the early Holocene rise of sea level to present levels in the central and northern Strait of Georgia.

Northeastern British Columbia was occupied by the Cordilleran (CIS) and the Laurentide (LIS) ice sheets, however, the timing and extent remains contentious. The late Quaternary and Holocene history of this area is examined by exploring geomorphic, stratigraphic, geochemical and geochronologic components of glacial, deglacial, paraglacial and non-glacial landsystems. New tools, such as GIS, LiDAR, and new geochronologic methods, such as optical dating are used to understand the Quaternary geology and geomorphology of the region. Bedrock topography represents the base of the Quaternary section and modelling shows that paleovalleys, common in this region, host extensive Neogene sedimentary records. Stratigraphies from the Murray and Pine valleys indicate glaciation prior to the Mid-Wisconsinan (MIS 3) and during the Late Wiconsinan (MIS 2). Glacial landforms record Late Wisconsinan ice-sheet coalescence and reflect the complex interaction of the LIS and CIS margins. During deglaciation, the LIS and CIS separated and glacial Lake Peace (GLP) formed. Shoreline features enable reconstruction of lake and ice configurations. Four phases of GLP are preserved. Optical ages from Phase II indicate GLP occupied the area some time between ca. 16 – 14 ka yrs ago. The apparent tilt on the shorelines provides a measure of isostatic adjustments and suggests asynchronous retreat of first the LIS, then the CIS. The transition from paraglacial to boreal conditions was driven by climate change and is recorded by vegetation sucession and cessation of paraglacial processes. Optical ages from stabilized dunes and radiocarbon ages from organics date the transition between 12 – 11.5 ka yrs ago with full boreal conditions established by 10 ka yrs ago. The Holocene is dominated by erosional processes, however some systems are aggrading. A case study on a floodplain demonstrates that resistivity (Ohmmapper) surveys provide a grain-size proxy to suppliant GPR studies, which is essential for geophysical fluvial architectural analysis. In the study, the discrepancy between planform style (classic meander model) and subsurface geophysical surveys (indicative of vertical accretion associated with braided and wandering fluvial styles) reiterates cautions that planform may not always be a functions of depositional process and one may not be used to predict the other. / Graduate / 0372 / 0373 / 0368