No description available.

Glacial sediments are found across formerly glaciated regions across the world and host a variety of important resources, ranging from groundwater to hydrocarbons, aggregate material, and mineral deposits. In southern Ontario, Canada, thick successions (up to 200 m) of Quaternary glacial sediments are truncated by large valleys (>30km long, 2 to >8.5 km wide, and up to 200m deep) that formed subglacially and have characteristic morphology and infill stratigraphy. These valleys are interpreted as (a new class of) tunnel valleys and strongly affect groundwater resources and flow systems at local and regional scales. The overall context of the valleys is evaluated through an introduction to the study area, objectives, and background information on subglacial systems and geologic history of south-central Ontario (Chapter 1). Interpretation of valley genesis in Simcoe County is provided through an integrated, multi-faceted approach, involving: description of the morphology and sediment infill succession within the valleys from surficial mapping, sedimentological logging of continuously-cored boreholes, and geophysical surveys (Chapter 2); delineation and characterization of seismic architecture from high-resolution lake-based sub-bottom profiles in one of the valleys (Chapter 3); detailed site-scale field description of the internal characteristics of the regional Late Wisconsin till sheet in various subglacial settings (Niagara Escarpment, uplands, lowlands; Chapter 4); comparison of the characteristics of the subglacial bed within the study area to adjacent regions in southern Ontario (Chapter 5); and a synthesis of the major findings from all the different components of this investigation and suggestions for future work to shed further light on several questions that arise from this study (Chapter 6). Together, key data from these studies of tunnel valleys and related deposits – a near-continuous till sheet on the surface of uplands and along the flanks and floors of the tunnel valleys, multi-stage drumlinization of the till sheet following development of the tunnel valleys, variations in internal facies and physical properties within the till sheet in different subglacial settings, localized distribution of coarse-grained tunnel valley in-fill sediments, and gradational upward transitions from tunnel valley in-fills to fossiliferous proglacial lacustrine sediments – indicate multiple phases of subglacial meltwater, and direct subglacial, erosion and deformation contributed to the development of the valleys over a protracted time period during the Late Wisconsin. Landform and sediment associations within the valleys in Simcoe County and surrounding parts of the bed of the former Laurentide ice sheet in south-central Ontario, are inconsistent with previous conceptualizations involving the presence of large (>1000 km2) subglacial lakes and the storage and discharge of regional-scale subglacial meltwater sheetfloods followed by ice stagnation. This study provides new data and insight to help refine reconstructions and better understand the evolution of past ice dynamics and subglacial processes, evaluate competing theories of regional landscape evolution, and provide new conceptual and (hydro)stratigraphic frameworks for future hydrogeological investigations related to groundwater exploration and use. / Thesis / Doctor of Philosophy (PhD)

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

It is now accepted that the last British-Irish Ice Sheet (BIIS) was highly dynamic and drained by numerous fast flowing ice streams. This dynamic nature combined with its maritime location made the BIIS sensitive to the rapid climate change that characterised the Last Glacial Interglacial Transition. Gaining an understanding of the behaviour of the BIIS at this time is important to explore the nature of forcing between ice sheets and climate. This thesis presents new chronological data relating to the deglaciation of the northwest sector of the BIIS (NW-BIIS) from onshore dating of moraines using cosmogenic exposure dating. This improved chronological framework is supported by offshore data in the form of a newly constructed Ice Rafted Detritus (IRD) record from the offshore sediment core MD95-2007. These data suggest that deglaciation commenced sometime after 18 ka and that the NW-BIIS was located close to the present day shoreline by 16 ka. Further provenance analysis of the IRD using U-Pb dating of detrital minerals demonstrates that during the Last Glacial-Interglacial Transition MD95-2007 was being supplied distal IRD from a source(s) to the west. The absence of diagnostic Scottish material suggests that after retreat to the coastline at 16 ka calving margins were not re-established during Greenland Interstadial 1. By combining these results with existing data relating to the deglaciation of the NW-BIIS it is possible to summarise the deglaciation history of the NW-BIIS from the continental shelf to mountainous source regions and compare this to numerical models of BIIS behaviour during this time. With a better understanding of the chronology of NW-BIIS retreat it is possible to relate the timing of initial deglaciation to possible forcing factors and gain a better understanding of the response of a marine based sector of an ice sheet to rapid climate change.

Glacial isostatic adjustment (GIA) is the process by which the solid Earth responds to past and present-day changes in glaciers, ice caps, and ice sheets. This thesis focuses on vertical crustal motion of the Earth caused by GIA, which is influenced by several factors including lithosphere thickness, mantle viscosity profile, and changes to the thickness and extent of surface ice. The viscosity of the mantle beneath Antarctica is a poorly constrained quantity due to the rarity of relative sea-level and heat flow observations. Other methods for obtaining a better-constrained mantle viscosity model must be investigated to obtain more accurate GIA model predictions. The first section of this study uses seismic wave tomography to determine mantle viscosity. By calculating the deviation of the P- and S-wave velocities relative to a reference Earth model (PREM), the viscosity can be determined. For Antarctica mantle viscosities obtained from S20A (Ekstrom and Dziewonski, 1998) seismic tomography in the asthenosphere range from 1016 Pa∙s to 1023 Pa∙s, with smaller viscosities beneath West Antarctica and higher viscosities beneath East Antarctica. This agrees with viscosity expectations based on findings from the Basin and Range area of North America, which is an analogue to the West Antarctic Rift System. Section two compares bedrock elevations in Antarctica to crustal thicknesses, to infer mantle temperatures and draw conclusions about mantle viscosity. Data from CRUST 2.0 (Bassin et al., 2000), BEDMAP (Lythe and Vaughan, 2001) and specific studies of crustal thickness in Antarctica were examined. It was found that the regions of Antarctica that are expected to have low viscosities agree with the hot mantle trend found by Hyndman (2010) while the regions expected to have high viscosity are in better agreement with the trend for cold mantle. Bevis et al. (2009) described new GPS observations of crustal uplift in Antarctica and compared the results to GIA model predictions, including IJ05 (Ivins and James, 2005). Here, we have generated IJ05 predictions for a three layered mantle (viscosities ranging over more than four orders of magnitude) and compared them to the GPS observations using a χ2 measure of goodness-of-fit. The IJ05 predictions that agree best with the Bevis et al. observations have a χ2 of 16, less than the null hypothesis value of 42. These large values for the best-fit model indicate the need for model revisions and/or that uncertainties are too optimistic. Equally important, the mantle viscosities of the best-fit models are much higher than expected for West Antarctica. The smallest χ2 values are found for an asthenosphere viscosity of 1021 Pa•s, transition zone viscosity of 1023 Pa∙s and lower mantle viscosity of 2 x 1023 Pa∙s, whereas the expected viscosity of the asthenosphere beneath West Antarctica is probably less than 1020 Pa∙s. This suggests that revisions to the IJ05 ice sheet history are required. Simulated annealing was performed on the ice sheet history and it was found that changes to the recent ice load history have the strongest effect on GIA predictions. / Graduate

La Côte Nord de l’Estuaire et du Golfe du St. Laurent (Québec, Canada) est caractérisée par une série de complexes deltaïques mis en place en contexte de chute de Niveau Marin Relatif (NMR) forcée par le rebond glacio-isostatique, lors du retrait des marges de l’Inlandsis Laurentidien. L’étude de trois complexes deltaïques montre qu’un motif d’évolution morphostratigraphique contrôlé par le retrait des marges glaciaires prévalait au premier ordre. Le premier stade de sédimentation se caractérise par la mise en place de cônes d’épandage juxtaglaciaires sous-aquatiques. Leur localisation au front de la marge glaciaire fait que la répartition spatiale des corps sédimentaires montre localement un empilement en rétrogradation. Dès l’émergence d’un front glaciaire continental, des deltas proglaciaires se développent en contexte de chute de NMR, formant des lobes dont l’altitude décroît vers le bassin. Ces derniers sont initialement associés à un système fluviatile en tresse alimenté en sédiments glaciogéniques par les marges glaciaires en retrait. Malgré des taux de chute de NMR de plusieurs cm/an, aucune incision fluviatile n’est observée à ce stade et la dynamique de transfert est prédominante du fait des forts taux d’apports sédimentaires. Plus tard, à la suite du retrait des marges glaciaires hors des bassins versants, le remaniement paraglaciaire se développe du fait de la réduction drastique des apports en eaux et sédiments. Le système fluviatile, devenu méandriforme, s’encaisse dans les anciens dépôts deltaïques maintenant inactifs et les bordures de deltas sont remaniées par les processus littoraux (plages soulevées, terrasses marines). Cette étude révèle que la majorité du volume de ces complexes deltaïques (10-20 km3) est mis en place en quelques milliers d’années seulement, immédiatement après la déglaciation ; le remaniement paraglaciaire n’ayant contribué à l’accrétion deltaïque que très marginalement. La modélisation numérique (Dionisos) valide les différents forçages identifiés par l’analyse morphosédimentaire. Une synthèse des complexes deltaïques à l’échelle de toute la Côte Nord du St. Laurent a permis de catégoriser 21 complexes deltaïques en 4 scénarios d’évolution morphosédimentaire, directement liés à la dynamique de retrait de la marge glaciaire. Leur décryptage offre une clef de lecture originale pour l’historique du retrait des marges glaciaires de l’Inlandsis Laurentidien sur la période 12-7.5 ka cal BP. / The North Shore of the St. Lawrence Estuary and Gulf is characterized by the presence of deltaic complexes that were emplaced under falling Relative Sea Level (RSL) forced by the glacio-isostatic rebound, immediately after the retreat of the Laurentide Ice Sheet (LIS) margin. The study of three deltaic complexes reveals that a common morphostratigraphic evolution forced by the retreat of the LIS prevailed for the edification of these structures, reflecting the retreat of the LIS margin. The emplacement of subaqueous outwash fans beyond the retreating or stillstanding glaciomarine margin constitutes the first stage of this evolution. As tied to the ice-margin position, these fans are characterized in places by a backstepping stacking pattern, in spite of the forced regressive setting. From the emergence of a continental ice front, the proglacial deltaic system develops and forms lobes staged accordingly to the RSL fall. These deltaic systems actively prograde at that time because they are fed in glaciogenics by the retreating LIS margin through braided meltwater streams. In spite of the RSL fall reaching several centimeters per years, no fluvial entrenchment occurs mainly owing to the significant amount of sediment supply. Later, when the LIS margin retreats from the drainage basins of feeding rivers, fluvial systems experience a drastic drop in sediment supply that forced the interruption of the deltaic progradation and the onset of paraglacial reworking. The paraglacial reworking consists in the entrenchment of a meandering fluvial system within former deltaic deposits and shows the prevalence of shallow-marine processes (waves, tides) at the delta rim (raised beaches, marine terraces). This study reveals that the bulk of the deltaic volume (c.a. 10-20 km3) for each complex was emplaced in only a few thousands of years following the LIS margin retreat when the latter was still located in the drainage basin. The paraglacial reworking had a minor influence on the deltaic accretion. A forward stratigraphic model (Dionisos) is used to validate the variety of forcing as understood from the sedimentary analysis. A synthesis including 21 deltaic complexes of the St. Lawrence North Shore allowed the establishment of a fourfold categorization. This scheme of deltaic evolution was used in order to refine the position of the LIS margin retreating upland for a period ranging from 12 to 7.5 ka cal BP.

The recent discovery of the 31-km-wide Hiawatha impact crater has raised unanswered questions about its age, impactor and highly unusual organic carbon component. Previous research suggests a fractionated iron meteorite impactor, a probable maximum 3–2.4 Ma impact age and a possible Younger Dryas impact age. The first objective in this study has been to investigate a possible link between the Cape York meteorites and the Hiawatha impact crater by comparing the chromium isotopic signature in chromite from a Cape York meteorite with the chromium isotopic signature in potential chromite from the Hiawatha impactor. The second objective has been to investigate a possible Hiawatha signature in the Younger Dryas deposits from Baffin Bay. The third objective has been to study the organic carbon component in impactites derived from the Hiawatha impact crater. Heavy mineral grains were separated from glaciofluvial sediment which contains Hiawatha impactite grains. Not a single chromite grain was found and the possible link to the Cape York meteorites could not be tested. The petrographic examination of Younger Dryas marine deposits resulted in absence of impact-related Hiawatha grains. A petrological investigation revealed that organic carbon was likely found in five of six variably shocked impactites derived from the Hiawatha impact crater. The character of the organic carbon varies between the samples and also within individual samples. Vitrinite reflectance measurements of the organic carbon in two impactites yielded low reflectance values compared to charcoalification experiments of wood. Organic particles with different reflectance in the same sample suggest that the particles had different impact histories prior to settling and becoming a rock. Diagnostic conifer cellular texture was found in at least one of the samples. The character of the organic particles in the impactites supports the suggestion in a previous study that the sources of the Hiawatha organic carbon component are unmetamorphosed surficial deposits containing dead conifer tree trunks and fine-grained layered clay and organic matter.  In this study it is concluded that the apparent absence of chromite in the examined glaciofluvial sediment sample corroborates the significance of previous research which suggests that the Hiawatha impactor was an iron meteorite. The apparent absence of impact related grains in the Younger Dryas deposits suggests that although a Younger Dryas age for the Hiawatha impact crater is less likely now, the possibility remains open. The organic carbon with diagnostic conifer cellular texture in the Hiawatha impactites corroborates the conclusion in a previous study that the Hiawatha impact-related organic carbon component stems from local, thermally degraded conifer trees with a probable age of ca. 3–2.4 Ma. It is also concluded that the relatively low reflectance values of the organic carbon in the Hiawatha impactites seem to be related to the short duration of the high-temperature excursion during the hypervelocity impact event.

<p>Understanding climate-driven changes in global land-based ice volume is a critical component in our capability to predict how global sea level will rise as a consequence of the current human-driven climate change. At the last glacial maximum (LGM, which peaked around 20 ka), ephemeral ice sheets covered vast regions of the northern hemisphere while both the Greenland and Antarctic ice sheets were more extensive than at present. As global temperatures rose at the transition into the Holocene, driving the LGM deglaciation, eustatic sea level rose by approximately 125 m. The east Antarctic ice sheet (EAIS) is the largest ice sheet on Earth today, holding an ice volume equivalent to ca. 53 m rise in global sea level. Considering current trends in global climate, specifically rapidly increasing atmospheric CO₂ levels and global temperature, it is important to improve our understanding of how the EAIS will respond to global warming so that we can make better predictions of future sea level changes to guide community adaptation and planning efforts. Numerical ice sheet models which inform projections of future ice volume changes, and can, therefore, yield projections of sea level rise, rely on empirical data to test their ability to accurately represent former and present ice configurations. However, there is a general lack of data on the paleoglaciology of the EAIS along the western Dronning Maud Land (DML) margin. In order to address this situation, the paleoglaciology of western DML forms the focus of the work presented in this thesis.</p><p><b> </b></p><p>Together with collaborators within the MAGIC-DML consortium (Mapping, Measuring and Modelling Antarctic Geomorphology and Ice Change in Dronning Maud Land) that provides the funding for this MS project, the author has performed geomorphological mapping across western DML; an area of approximately 200,000 km². The results of the mapping presented in this thesis will provide the basis for a detailed glacial reconstruction of the region. The geomorphological mapping was completed almost entirely by remote sensing using very high-resolution (sub-meter in the panchromatic) WordView-2 and WorldView-3 (WV) satellite imagery, combined with ground validation studies during field work. Compared to Landsat products, the improved spatial resolution provided by WV imagery has fundamentally changed the scale and detail at which remote sensing based geomorphological mapping can be completed. The mapping presented here is focused on the glacial geomorphology of mountain summits and flanks that protrude through the ice sheet’s surface (nunataks). In our study area of western DML these nunatak surfaces make up <0.2 % of the total surface area, and the landforms mapped here are generally smaller than can be identified from Landsat products (30 m spatial resolution). The detail achieved in our mapping, across such a vast, remote area that presents numerous obstacles to accessibility highlights the benefits of utilizing the new VHR WV data. As such an evaluation of the WV data, as applied to geomorphological mapping is presented here together with our mapping of the glacial geomorphology of western DML. The results of which provides evidence of ice having overridden sites at all elevations across the entire study area; from the highest elevation inland nunataks that form the coast-parallel escarpment, to low-elevation emerging nunataks close to the coast. Hence from our studies of the glacial geomorphology of this region we can ascertain that, at some point in the glacial history of western DML, ice covered all of the mountain summits that are exposed today, indicating an ice sheet surface lowering of up to 700 m in some places.</p>

L’Ordovicien (485–444 Ma) est une période géologique caractérisée par laconcomitance d’une glaciation majeure et de l’une des 5 plus grandes extinctions de masse del’histoire de la Terre. Cette thèse avait pour objectif d’améliorer la compréhension de l’évolutiondu climat à cette époque à l’aide de la modélisation numérique, ain de fournir une imagecohérente de la glaciation. Nous avons d’abord démontré que la coniguration continentaleordovicienne induit une dynamique océanique particulière, à l’origine d’une instabilité climatiquepermettant un refroidissement brutal du climat global sans variation importante de laconcentration atmosphérique en CO2 (pCO2). Dans un deuxième temps, un modèle innovantcouplé climat-calotte a permis de produire la première simulation de la mise en place de la glaciationsupportée par les données géologiques, sous un scénario cohérent de baisse de la pCO2.Les résultats indiquent que les premières glaces continentales se seraient mises en place dèsl’Ordovicien Moyen (465 Ma), quelque 20 millions d’années plus tôt qu’initialement envisagé.Dans ce scénario, le franchissement de l’instabilité climatique ordovicienne marque la miseen place du maximum glaciaire au cours de l’Ordovicien terminal Hirnantien (445–444 Ma).Des expériences réalisées avec un modèle de végétation primitive montrent que le développementdes plantes non-vasculaires a pu constituer le mécanisme à l’origine de la chute de lapCO2, via une intensiication de l’altération des surfaces continentales. Enin, les interactionsentre climat et biosphère marine ont été envisagées selon 2 axes complémentaires. (i) De nouvellescontraintes ont été fournies pour comprendre la paléobiogéographie des communautésmarines, par la publication de cartes de la circulation océanique de surface modélisée sousdiférentes pCO2 au cours de l’Ordovicien Inférieur, Moyen et Supérieur. (ii) Les relationsentre variations climatiques et état redox de l’océan ont été étudiées avec un modèle d’océanrécent bénéiciant d’un module de biogéochimie marine (MITgcm). Les simulations suggèrentdes anoxies partielles (durant le Katien) ou globales (durant le Silurien inférieur) au cours dela transition Ordovicien–Silurien. Elles démontrent également que l’extinction de l’Ordovicienterminal ne serait pas liée à un évènement d’anoxie. / The Ordovician (485–444 Ma) is a geological period characterized by theconcomitance of a major glaciation and one of the “Big Five” mass extinction events thatpunctuated the Earth’s history. This dissertation aimed at developing a better understandingof the climatic evolution at that time through numerical modeling, in order to providea consistent picture of the glaciation. First, it was shown that the Ordovician continentalconiguration leads to a particular ocean dynamics, which induces in turn the development ofa climatic instability that allows global climate to cool suddenly in response to subtle changesin the atmospheric partial pressure of CO2 (pCO2). Secondly, an innovative climate-ice sheetcoupled model produced the irst simulation of the glaciation that is supported by geologicaldata, in the context of a decrease in pCO2. Results show that glacial onset may have occurredas early as the Mid Ordovician (465 Ma), i.e., some 20 million years earlier than thoughtinitially. In this scenario, the climatic instability is reached during the latest Ordovician andaccounts for the onset of the Hirnantian glacial maximum (445–444 Ma). Experiments conductedwith a non-vascular vegetation model reveal that the origination and expansion of theirst land plants signiicantly intensiied continental weathering during the Ordovician andpotentially drove the drop in atmospheric CO2. Finally, the interactions between climate andthe marine biosphere were investigated based on 2 complementary axes. (i) News constraintson the paleobiogeography of marine living communities were brought through the publicationof maps showing the ocean surface circulation modeled at various pCO2 levels during theEarly, Middle and Late Ordovician. (ii) The relationships between climatic variations andthe redox state of the ocean were studied using a recent ocean model with biogeochemical capabilities(MITgcm). The simulations suggest partial and global oceanic anoxic events duringthe Katian and the early Silurian respectively. They also show that anoxia is probably notresponsible for the latest Ordovician mass extinction event.

Cette thèse est consacrée à l’étude de l’interaction cryosphère-climat depuis le milieu du Pliocène jusqu’au quaternaire pendant l’installation pérenne de la calotte groenlandaise. Nous étudions d’abord les causes du développement et de la disparition de l’importante mais courte glaciation qui a eu lieu pendant le stade isotopique marin M2 (MIS M23.264-¬3.312 Ma). Ensuite, dans le cadre du programme international sur la modélisation du Pliocène (PLIOMIP2), nous étudions le climat de la période chaude du Plaisancien moyen(MPWP, 3.3-3.0Ma). Enfin, la troisième période étudiée est la transition Plio-Pléistocène transition (PPT, 3.0-2.5Ma), que nous avons étudiée grâce à un couplage asynchrone entre un modèle de climat et un modèle de calotte. A travers ces différentes périodes, nous avons amélioré la connaissance des relations entrepCO2, tectonique et climat pendant la transition d’un monde chaud et riche en CO2 vers le monde bien plus froid et à faible pCO2 des glaciations quaternaires. Ce résultat montre l’importance de mieux comprendre les relations entre dynamique océanique, pCO2 et climat. / This thesis is devoted tounderstanding the interaction betweencryosphere and climate from the mid Plioceneto the early Quaternary during the onset ofNorthern Hemisphere Glaciation (NHG).Firstly, we investigate the causes for thedevelopment and decay of the large but shortliving glaciation that occurred during MarineIsotope Stage 2 (M2, 3.264-¬3.312 Ma);Secondly, in the framework of the internationalPliocene Model Intercomparison Project(PLIOMIP2), we study the climate of Mid-Piacenzian Warm Period (MPWP, 3.3-3.0Ma).Thirdly, we explore the Plio-PleistoceneTransition (PPT, 3.0-2.5Ma) with anappropriate asynchronously coupled climatecryosphere model. Through these differentperiods, we provide a better understanding ofthe relationship between pCO2, tectonics andclimat during the transition from a warm andhigh-CO2 world to the cold and low-CO2Quaternary glaciations. This work also pointsout the necessity to further study the linkbetween ocean dynamics, carbon cycle andclimate.