The geology, hydrogeology, and groundwater geochemistry are described for an interface aquifer in Eastern Ontario exhibiting anomalously high proportions of iodine (I) as iodide (I-) and dissolved methane (CH4). The studied area is unique in that it shows a significant marine influence, attributed to the most recent Champlain Sea incursion 10 – 12 ka BP, which has implications for I and CH4 enrichment. I and CH4 in groundwater are found in high proportions in reducing fossil seawaters, which are typically observed in depressions in the bedrock surface that are overlain by thick layers of glaciomarine muds. I is released via microbial decomposition of marine phytoplankton into mud porewaters, where it is then leached to underlying groundwaters. 129I and 14C isotopic signatures of I and C compounds highlight the importance of allochthonous I and C sources in the Champlain Sea basin derived from glacial abrasion of the surrounding terrain and imported via glacial meltwater. CH4 is microbial in origin, with marine phytoplankton from the Champlain Sea incursion and ancient terrestrial organic matter from an Early Wisconsinian interstadial period (60 – 75 ka BP) as the dominant substrates. A thermogenic CH4 component was observed for areas underlain by the Billings shale unit. Both I and CH4 originate at least partially from the same marine phytoplankton source within the muds, demonstrate similar controls on enrichment, and have a Spearman’s rank coefficient of 0.62, indicating that the correlation between I and CH4 in groundwater in the studied area is significant.

Sediment accumulated in lakes stores valuable information about past environments and paleoclimatological conditions. Cores previously obtained from Saint Albans Bay, located in the Northeast Arm of Lake Champlain, VT record the transition from the Champlain Sea to Lake Champlain. Belrose (2015) documented the presence of a peat horizon separating the sediments of the Champlain Sea from those of Lake Champlain. Initially, this layer was thought to comprise the transition from the marine environment of the Champlain Sea to a freshwater wetland. However, based on the results from this study, the transition between marine and freshwater conditions is thought to be represented by an erosional unconformity, indicative of a lowstand at the end of the Champlain Sea period. For this study, five additional cores were collected from Saint Albans Bay along a transect following the long axis of the bay moving into progressively deeper water. These cores better constrain the spatial extent, thickness and age variability of the peat layer within the bay and allow us to better understand the environmental conditions that preceded the period of peat deposition. In each of the cores there is evidence of sediment reworking in the uppermost Champlain Sea sediments, indicated by the presence of coarse-grained sediment, which is suggestive of a lowstand at the end of the Champlain Sea period before the inception of Lake Champlain. This coarse-grained layer is immediately overlain by a thick peat horizon. The widespread occurrence of the peat layer points to a large wetland that occupied the entire inner portion of Saint Albans Bay, and lake level ~ 9 m lower than at present during the Early Holocene. Based on radiocarbon dating, this paleo-wetland existed in Saint Albans Bay from ~ 9,600-8,400 yr BP. The development of this wetland complex is time transgressive, reflecting rapidly increasing lake level during the Early Holocene. This hypothesis is supported by the basal peat radiocarbon dates, as well as by the composition of plant macrofossils recovered from the peat horizons. The shift from peat deposition to fine-grained, low organic content lacustrine sedimentation is believed to have occurred at ~8.6-8.4 ka and is likely the result of continued isostatically driven lake level rise coupled with a changing climate. Although it was not its primary focus, this study also seeks to address the variations in sediment composition in the Lake Champlain sections of the cores. Evidence from the Lake Champlain record in Saint Albans Bay indicates that there were notable fluctuations in sedimentation, which were likely linked to both climatic variations and a change in the morphology of the bay. The rebound in productivity from ~8-5 ka is likely the result of warmer conditions during the Hypsithermal period. An increase in terrigenous sedimentation during this same time suggests a change in the morphology of the bay in which the Mill River delta migrated towards the inner bay. Initially, the cooler conditions of the Neoglacial are reflected in Saint Albans Bay by a decrease in organic matter content from ~5-3 ka. During the latter part of the Neoglacial (~3-1 ka), increases in organic matter content and detrital input point to enhanced productivity in response to increased precipitation and runoff from the watershed. The most recently deposited sediments in Saint Albans Bay bear out the legacy of anthropogenic nutrient enrichment of the bay in the form of increased algal productivity.

The City of Ottawa is situated in an area known as the Champlain Sea, 17,000 years before present (BP) the entire area was covered with sea water. This area deposited marine clays which are known to be highly sensitive. The City of Ottawa needs to expand land use to allow for the expansion of infrastructure and housing to support its growth. This study is intended to assist the City of Ottawa’s geotechnical engineers in their decision-making by identifying future sensitive areas prone to landslides due to rainfall based on future climate model data. The project incorporates rainfall intensities from downscaled climate model data in the Transient Rainfall Infiltration and Grid-based Regional Slope-Stability (TRIGRS) model to investigate areas sensitive to landslides, then within a GIS platform, the future landslide susceptibility maps were created based on Factor of Safety (FS) values showing the areas prone to landslides. The data input for the model includes climate model data, topography, hydrogeology, geology and geophysical data obtained from a previous study. These data were prepared using ArcGIS software and converted into ascii format for TRIGRS model. The model was calibrated using historical rainfall intensities and validated by comparing to historical landslide areas. Sensitivity analysis were performed to ranges of geotechnical properties found within sensitive marine clays in the area to find the values best to create the ideal scenario, normal scenario and worst-case model scenario for the prediction. Rainfall intensities from projected climate data Intensities Duration Frequency (IDF) of 10 years and 50 years returning period and rainfall intensities of 12 hr, 24 hr, and 48 hr were selected for the model. Results from simulations find the projected climate rainfall intensity do not have impact or has minimal impact to slope stability in sensitive marine clay areas in Ottawa directly. However, higher rainfall runoff is expected from projected rainfall RCP8.5 than the RCP4.5. The infiltration rate remains constant throughout each simulation, which is the same value as the hydraulic conductivity. The time when the slope becomes unstable varies depending on initial water levels. Results from the ideal and normal scenario show no areas prone to slope failure after 48 hours of rainfall duration. However, the factor of safety decreases as the rainfall duration increases and is expected to decrease with longer rainfall durations. The worst-case scenario shows some areas prone to slope failure (FS < 1) with 2% probability of slope failure at 48 hours of rainfall duration. The distribution of these unstable areas are located along the Ottawa River, Rideau River, Carp River, Mississippi River and valleys along their tributaries, the majority of the area prone to slope instability from rainfall are in the east part of the City of Ottawa. While there are many uncertainties and limitations which contribute to the model results, this study is useful to engineers and planners in initial implementation of mitigation strategies to mitigate the damages and cost from landslides events. The susceptibility maps can also assist in decision making for planners in developing into these areas.

Les connaissances sur la morphologie de la Terre sont essentielles à la compréhension d’une variété de processus géomorphologiques et hydrologiques. Des avancées récentes dans le domaine de la télédétection ont significativement fait progresser notre habilité à se représenter la surface de la Terre. Parmi celles-ci, les données d’élévation LiDAR permettent la production de modèles numériques d’altitude (MNA) à haute résolution sur de grands territoires. Le LiDAR est une avancée technologique majeure permettant aux scientifiques de visualiser en détail la morphologie de la Terre et de représenter des reliefs peu prononcés, et ce, même sous la canopée des arbres. Une telle avancée technologique appelle au développement de nouvelles approches innovantes afin d’en réaliser le potentiel scientifique. Dans ce contexte, le présent travail vise à développer deux approches de cartographie numérique utilisant des données d’élévation LiDAR et servant à l’évaluation de la composition du sous-sol. La première approche à être développée utilise la localisation de crêtes de plage identifiées sur des MNA LiDAR afin de modéliser l’étendue maximale de la mer de Champlain, une large paléo-mer régionalement importante. Cette approche nous a permis de cartographier avec précision les 65 000 km2 autrefois inondés par la mer. Ce modèle sert à l’évaluation de la distribution des sédiments marins et littoraux dans les basses-terres du Saint-Laurent. La seconde approche utilise la relation entre des échantillons de matériel parental des sols (MPS) et des attributs topographiques dérivés de données LiDAR afin de cartographier à haute résolution et à une échelle régionale le MPS sur le Bouclier canadien. Pour ce faire, nous utilisons une approche novatrice combinant l’analyse d’image orientée-objet (AIOO) avec une classification par arbre décisionnel. Cette approche nous a permis de produire une carte du MPS à haute résolution sur plus de 185 km2 dans un environnement hétérogène de post-glaciation. Les connaissances issues de la production de ces deux modèles ont permis de conceptualiser la composition du sous-sol dans les régions limitrophes entre les basses-terres du Saint-Laurent et le Bouclier canadien. Ce modèle fournit aux chercheurs et aux gestionnaires de ressources des connaissances détaillées sur la géomorphologie de cette région et contribue à l’amélioration de notre capacité à saisir les services écosystémiques et à prédire les aléas environnementaux liés aux processus du sous-sol. / Knowledge of the earth’s morphology is essential to the understanding of many geomorphic and hydrologic processes. Recent advancements in the field of remote sensing have significantly improved our ability to assess the earth’s surface. From these, LiDAR elevation data permits the production of high-resolution digital elevation models (DEMs) over large areas. LiDAR is a major technological advance as it allows geoscientists to visualize the earth’s morphology in high detail, even allowing us to resolve low-relief landforms in forested areas where the surface is obstructed by vegetation cover. Such a technological advance calls for the development of new and novel approaches to realize the scientific potential of this new spatial data. In this context, the present work aims to develop two digital mapping approaches that use LiDAR elevation data for assessing the earth’s subsurface composition. The first approach to be developed uses the location of low-relief beach ridges observed on LiDAR-derived DEMs to map the extent of a large and regionally important paleo-sea, the Champlain Sea. This approach allowed us to accurately map the 65,000 km2 area once inundated by sea water. The model serves to the assessment of the distribution of marine and littoral sediments in the St. Lawrence Lowlands. The second approach uses the relationship between field-acquired samples of soil parent material (SPM) and LiDAR-derived topographic attributes to map SPM at high-resolution and at a regional scale on the Canadian Shield. To do so, we used a novel approach that combined object-based image analysis (OBIA) with a classification tree algorithm. This approach allowed us to produce a fine-resolution 185 km2 map of SPM in a heterogeneous post-glaciation Precambrian Shield setting. The knowledge obtained from producing these two models allowed us to conceptualize the subsurface composition at the limit between the St. Lawrence Lowlands and the Canadian Shield. This insight provides researchers and resource managers with a more detailed understanding of the geomorphology of this area and contributes to improve our capacity to grasp ecosystem services and predict environmental hazards related to subsurface processes.