La minéralisation aurifère de la propriété Amaruq (Nunavut, CA) est recouverte par divers sédiments glaciaires ayant formé des traînées de débris minéralisés sur des surfaces importantes. Ce projet a pour objectif de caractériser la signature minéralogique et géochimique du gisement aurifère dans le till de surface afin de développer des méthodes de prospection pour les travaux d’exploration. 61 échantillons de till ont été récoltés le long de quatre transects pour déterminer les variations de composition dans la direction principale de l’écoulement glaciaire (NNO). Quatre profils verticaux d’ostioles, un par transect, ont été échantillonnés pour documenter la variation de la signature géochimique avec la profondeur. La signature chimique de la scheelite dans le gisement est comparée à celle retrouvée dans le till et permet de tester son utilisation comme outil de prospection glacio-sédimentaire. Trois des quatre profils d’ostioles présentent des concentrations stables des éléments traces avec la profondeur, indiquant une bonne homogénéisation du till par les processus cryogéniques. Les ostioles sont donc un milieu d’échantillonnage pertinent en région de pergélisol. Les analyses en composantes principales (ACP) de la géochimie des échantillons d’une campagne de till fournie par Mines Agnico Eagle ltée (AEM), et les données des transects, décrivent la variance de la roche encaissante mafique/ultramafique et de la minéralisation aurifére (CP1 et CP2 respectivement). Le krigeage des scores de la CP1 et CP2 indiquent des trains de dispersion vers le NNO, en aval glaciaire des zones minéralisées. Les comptages des minéraux indicateurs (MI : or et scheelite) augmentent de façon importante à environ 1.4 km en aval de la minéralisation dans le secteur est de la propriété, dans un till distal formant des formes drumlinoides. Dans le secteur ouest, recouvert d’un till proximal formant des crêtes morainiques, les comptages des MI et les scores de la CP1 augmentent directement en aval de la minéralisation. Les différences de distance de transport entre ces deux secteurs reflètent un changement de la dynamique glaciaire et démontrent l’importance du contexte local dans les travaux de prospection glaciaire. La comparaison de la signature géochimique de la scheelite entre le gisement et et des échantillons de till du transect de Whale Tail, grâce à des diagrammes discriminants et spectres des terres rares, montrent d’importantes similitudes, démontrant sa potentielle utilisation en prospection glacio-sédimentaire pour des gisements de type orogéniques. / On the Amaruq property (Nunavut, CA), the gold mineralization is covered by a patchwork of glacial sediments, including trains of mineralized debris dispersed over large surfaces. This project aims to define the deposit mineralogical and multi-elemental signatures in surface sediments to guide exploration in the area. A total of 61 samples were collected from frost boils along four NNW transects, parallel to the major ice flow direction, in order to assess the spatial variability of the deposit signature using till matrix geochemistry and indicator minerals. Four profiles in frost boils were sampled, one in each transect, in order to document trace elements depth variation. The scheelite geochemical signature from the deposit is compared to that in till from the Whale Tail transect to test its applicability in drift prospecting. Three out of four frost boil profiles exposed constant trace elements concentrations at depth, reflecting a good till homogenization by cryoturbation. Therefore, frost boils are a suitable medium of sampling for exploration purpose in permafrost terrain. Principal component analysis (PCA) of a previous till survey by Agnico Eagle Mines ldt (AEM) and transect samples define the ultramafic/mafic host rock of the mineralization and a pathfinder element suite characteristic of lode gold deposits (PC1 and PC2, respectively). Kriging of PC1 and PC2 defines a dispersal train oriented NNW, down ice of known mineralization. In the eastern part of the property, covered mostly by a distal till forming drumlinoid ridges, indicator mineral counts (IM: gold and scheelite) increase to a maximum approximately 1.4 km down ice from the outcropping zone. Western transects, largely covered by a proximal till forming morainal ridges, show a sharp increase of counts directly down-ice of mineralization. The difference in glacial transport between these two sectors represents a local change of ice dynamic and shows the importance of understanding the local context in exploration when interpreting glacial dispersal patterns. The comparison of the trace element signature between the deposit and the till scheelite using discriminant diagrams and REE abundance shows similarities, supporting its usefullness in drift prospecting for orogenic gold deposits.

No description available.

L'érosion de contact est un processus d'érosion interne pouvant être retrouvé dans les barrages en remblai et les digues qui peut mener à la rupture de l'ouvrage. Elle se définit comme le détachement et l'entraînement au travers d'un sol grossier de particules d'un sol fin, sous l'effet d'une sollicitation hydraulique parallèle à l'interface entre les deux matériaux. Cette étude vise à comprendre et à identifier les processus contrôlant l'érosion de contact avec un sol fin étalé (till de Romaine 3) lorsque combiné à différents matériaux grossiers. Pour étudier le processus d'érosion, un dispositif expérimental original a été développé. De plus, dans le but de bien comprendre les processus en jeu avec le sol fin étalé, des essais ont également été effectués avec des sols uniformes afin de servir de comparatif. Il s'avère que l'étalement granulométrique du sol fin a un impact important sur l'initiation et la progression de l'érosion de contact alors que différents mécanismes comme le pavage et le colmatage peuvent se développer et limiter l'érosion, ou même l'arrêter. En se basant les résultats expérimentaux de la présente étude ainsi que ceux recensés dans la littérature, un modèle d'évaluation du seuil de l'érosion est proposé. De plus, une approche expérimentale est proposée afin de faciliter l'étude de l'érosion de contact des sols fins pulvérulents à granulométrie étalée. Finalement, une conceptualisation de l'évolution de l'érosion de contact en fonction de la dimension des particules du sol fin et du sol grossier ainsi que de la sollicitation hydraulique est proposée. Avec cette dernière, la susceptibilité du till de Romaine 3 à l'érosion de contact est définie. / Contact erosion is an internal erosion process that can be found in embankment dams and dikes which can lead to failure. It is defined as the pullout and dragging of fines from a base material through an adjacent coarser material under seepage parallel to the interface. This study intends to understand and to identify the mechanisms that control contact erosion with a well-graded base soil (Romaine 3 till) combined with different filter materials. To study the erosion process, a new experimental setup has been developed. Also, in order to properly assess the mechanisms involved with a well-graded base soil, comparative tests have also been performed with poorly-graded base soils. It was found that the grading of the base material as an important impact on the initiation and progression of contact erosion as different mechanisms such as paving and clogging can develop to limit or even stop the erosion. With the experimental results of this study as well as those found in the literature, a model for contact erosion threshold evaluation is proposed. Additionally, an experimental approach is proposed in order to simplify the study of contact erosion for well-graded base soils. Finally, a conceptualization of the erosion process (initiation and evolution) depending on the grain size of both filter and base soils and on the hydraulic solicitation is proposed. The Romaine 3 till susceptibility to contact erosion is assessed with latter conceptualization.

L'évolution des calottes polaires est régie à la fois par une dynamique d'écoulement complexe et par des mécanismes tel le glissement à la base, la température de la glace ou le bilan de masse en surface. De plus, de nombreuses boucles de rétroactions sont constatées entre les différents phénomènes impliquées. Tout ceci rend la modélisation de cette évolution complexe. Malgré tout, un certain nombre de modèles ont été développés dans cette optique. Ceux-ci font tous intervenir des paramètres influents qui dans certains cas sont peu ou pas connus. Ils nécessitent donc d'être correctement spécifiés. L'assimilation de données peut permettre une meilleure estimation de ces paramètres grâce à l'utilisation d'observations qui sont peu nombreuses en glaciologie. Dans cette thèse, nous nous intéressons à la mise en place de systèmes d'assimilation performants pour deux problèmes inverses concernant l'évolution des calottes polaires. Pour mieux nous concentrer sur ce point, nous avons travaillé avec un modèle d'évolution de calotte simplifié (appelé Winnie) qui, cependant, représente bien la plupart des processus complexes de la dynamique de la glace, et permet de travailler à différentes échelles de temps. Dans un premier temps, nous mettons en place une approche 4D-Var pour la reconstruction de l'évolution d'un paramètre climatique influant sur l'évolution d'une calotte sur une échelle de temps typique de 20 000 ans. Elle nécessite notamment l'écriture du code adjoint du modèle. Dans un second temps, nous nous intéressons au problème du spin-up. Ce problème de calibration du modèle pour des simulations à échelle de temps courtes (pas plus de 100 ans) consiste plus particulièrement en la reconstruction conjointe de l'état initial, de la topographie du socle rocheux et des paramètres de glissement basal. Nous développons ici une approche filtre de Kalman d'ensemble pour résoudre ce problème.

Les glaces de l'Antarctique constituent un matériau d'archivage unique pour garder la trace de l'atmosphère du passé. La plus ancienne archive glaciaire disponible (800 ka) ne permet malheureusement pas d'étudier une transition climatique majeure qui s'est déroulé entre 1.2 Ma et 0.9 Ma (transition du mi-Pleistocène). Localiser un site ayant permis la conservation d'une glace aussi ancienne est donc un enjeu majeur au sein de la communauté des carottes de glace. Plusieurs régions autour du Dôme C présentent de bonnes caractéristiques d'un point de vue thermique et mécanique. Le but de ces travaux est d'identifier précisément des sites propices à un archivage de long terme, en utilisant différents modèles successifs. D'abord un modèle thermique 1D, qui permet de calculer le flux géothermique local, et d'affirmer que les hauteurs du socle sont exemptes de fusion depuis au moins 800 ka. Ensuite un modèle 3D d'écoulement, qui permet la description des trajectoires des particules de glace depuis la surface jusqu'au socle. Pour minimiser les perturbations stratigraphiques, il conviendra de placer le site de forage proche de la crête, et pour assurer une résolution suffisante le site ne devra pas se trouver à l'aplomb du sommet du relief basal, mais sur son flanc. Enfin, on utilise un modèle de datation 1D contraints par les isochrones, pour connaître la distribution des âges basaux et résolutions compatibles avec les isochrones observées, et estimer les risques de succès ou d'échec du forage. On propose deux sites de forage distants de 37 km de Dôme C, pour lesquels la présence de glace exploitable de 1.5 Ma et d'une résolution suffisante est hautement probable. / Antarctic ice is an exceptional archive of ancient atmospheric air. The oldest available ice archive (800 ka) is not old enough to study a main climatic transition that occurred between 1.2 Ma and 0.9 Ma (mid-Pleistocene transition). Locating a future old-ice drill site is now one of the main goal of the ice-core community. In the region of Dome C, several old-ice targets have required characteristics, from a thermal and mechanical point of view, to hosts very old ice. The goal of this work is to precisely locate suitable coring sites, using several glaciological models. First a 1D heat model to compute the local geothermal flux, that make us confident in the fact that bedrock heights have been free of basal melting for at least the last 800 ka. It is very likely that such sites are able to host a 1.5 Ma ice archive. Second, a 3D ice flow model, to describe the ice trajectories from the ice surface to the bed. Drilling less than 10 km from the surface ridge mimimizes the risks of stratigraphic disturbance, and drilling on the flanks of the bed relief - and not on its top - is mandatory to ensure a sufficient age resolution. Third, a 1D dating model constrained by isochrones let us describe the distribution of basal ages and age resolutions compatible with observed isochrones, and assess the risk of success or failure of the drilling. With currently available data and modelling tools, two drilling sites are proposed that are located 37 km from Dome C, at which 1.5 Ma-old ice is very likely to be archived in required conditions and with sufficient resolution.

Sur la Cordillère Royale en Bolivie, nous avons reconstitué les bilans de masse de 21 glaciers représentatifs, par photogrammétrie numérique. Les MNT's ont été construit à partir des clichés aériens de l'IGM-Bolivie de 1956, 1963, 1975, 1983, 1997 et 2006, avec une précision décimétrique ou métrique.<br />L'analyse des résultats sur le glacier du Zongo montre qu'un tiers de la surface du glacier (zone d'ablation) contrôle 80% du bilan de masse annuel. La série du Zongo constitue la série de bilan de masse glaciaire la plus longue et la plus précise des Andes tropicales.<br />Nous avons analysé les bilans de masse des 21 glaciers depuis 1956, lorsqu'on s'affranchit des tendances de ces glaciers, les variations temporelles du bilan de masse montrent un signal commun d'origine climatique.<br />Nous avons montré ainsi que 88% de la variance des bilans cumulés sont expliqués par l'exposition et l'altitude. Cette corrélation et les restitutions photogrammétriques ont permis de calculer les bilans de masse et la perte de surfaces de 376 glaciers au cours de la période 1963-2006 (43% de perte de leur volume et 48% de perte de leur surface depuis 1975).<br />Enfin, notre étude évalue la contribution des glaciers au système d'alimentation en eau potable de la ville de La Paz. En moyenne sur la période 1975-2006, la fonte glaciaire a contribué à hauteur de 18% et 15% sur la période récente 1997-2006 (12% en saison humide et 27% en saison sèche). Nous avons enfin calculé les débits attendus en cas de disparition totale des glaciers (12% en moins à l'année, 9% en moins en saison humide et 25% en moins en saison sèche).

Les glaciers intègrent naturellement des informations hydrologiques et énergétiques dans des régions climatiques variées du globe et sont sensibles à des changements minimes de bilan d’énergie de surface. Parmi les paramètres mesurés in-situ traduisant les variations climatiques, le bilan de masse de surface est la variable directement reliée aux conditions atmosphériques qui contrôlent les processus d’accumulation et d’ablation. Un grand nombre de modèles, de complexités diverses, existent pour simuler l’évolution du bilan de masse. Cependant, des questions perdurent quant à l’utilisation de l’approche la plus appropriée, notamment sur de longues périodes de temps. Basés sur de longues séries de mesures réalisées sur quatre glaciers des Alpes françaises (issues du service d’observation GLACIOCLIM), ces travaux de thèse ont pour but d’étudier la sensibilité du bilan de masse aux variables atmosphériques et topographiques, afin d’identifier celles essentielles pour contraindre les modèles.En s’appuyant tout d’abord sur une approche empirique, nos résultats indiquent que pour chaque glacier, la variabilité temporelle de l’ablation s’explique principalement par la variabilité de la température tandis que la variabilité spatiale de la fonte de la neige et de la glace est largement gouvernée par l’insolation. Néanmoins, à l’échelle du glacier, considérer l’insolation en plus de la température ne conduit pas à une amélioration de la performance des modèles empiriques et un modèle type degré-jour apparaît suffisant pour simuler son bilan de masse. L’utilisation du modèle à base physique SURFEX/ISBA-Crocus a permis d’étudier la sensibilité du bilan de masse estival à chacune des variables du bilan d’énergie de surface. Le modèle s’avère très performant pour la simulation des bilans spatialisés, mais nécessite un forçage météorologique précis et distribué spatialement. En particulier, le vent, souvent complexe à représenter dans les modèles atmosphériques,apparaît comme un facteur prépondérant pour modéliser correctement les processus de fonte. Avec des forçages météorologiques bien contraints, ces deux approches (empirique et physique) se révèlent très performantes quant à la modélisation du bilan de masse estival. Toutefois, nos résultats montrent que les bilans annuels restent très conditionnés par les bilans hivernaux du fait de la rétroaction de l’albédo de surface. Une attention particulière a donc été portée sur la simulation de l’accumulation.A partir des données de l’observatoire GLACIOCLIM, notre étude n’a pas permis d’établir de relation significative entre la variabilité spatio-temporelle des bilans hivernaux et les variables morpho-topographiques.Pour cette raison, une campagne spécifique basée sur des mesures LiDAR a été mise en œuvre au Col du Midi (Massif du Mont-Blanc) afin d’affiner le réseau de mesures et d’améliorer la compréhension de la variabilité spatio-temporelle de la couverture neigeuse. Nos résultats montrent la complexité d’obtenir des mesures d’accumulation par différence de MNT, du fait des vitesses d’émergence/submergence. / Glaciers naturally record hydrological and energetic information into varied climatic regions of the world, and are sensitive to small changes in their surface energy balance. The glacier surface mass balance is directly linked to the atmospheric variables which control accumulation and ablation processes. Many models, with different complexities, enable simulations of mass balance evolution, but some questions arise about the best approach to use, especially over long time periods. Based on long measurement series on four glaciers in the French Alps (GLACIOCLIM observatory), this thesis aims at studying the mass balance sensitivity to climatic variations, to identify relevant meteorological variables for modelling. Based on an empirical approach performed on each glacier, temperature was found to be the main driver of temporal ice/snow ablation variability while solar radiations strongly influence the spatial distribution of summer mass balance. However, at glacier scale, to include solar radiation in melt models does not improvethe performances and a classical degree-day model is sufficient to simulate glacier-wide mass balance.The snowpack model SURFEX-ISBA-Crocus was then used to study summer mass balance sensitivity to all surface energy balance fluxes. Results indicate good performances of the model to simulate mass balance at each stake, but it requires accurate meteorological forcing. In particular wind, generally difficult to assessby atmospheric models, appears as a key factor for an accurate ablation modelling. Our work emphasizes that both physical and empirical approaches are very efficient for ablation modeling when forced with accurate meteorological data. Yet, annual mass balances remain very sensitive to wintermass balance due to the surface albedo feedback. For this reason, simulating accumulation processes remainsas important as ablation modelling. Using GLACIOCLIM data measurements, our study failed to find relevant relationship between winter mass balance and topographical variables such as slope, aspect, elevation, due to insufficient spatio-temporal resolution. To fill this gap of measurements, we set up a dedicated field campaign to improve our understanding of the spatio-temporal variability of the snow depth distribution on Alpineglaciers. This year-round campaign was performed using LiDAR acquisitions at Col du Midi (Mont Blancmassif). Results underline the complexity of accurate accumulation measurements from DEM differences due to submergence/emergence velocities.

Les forages polaires profonds contiennent des enregistrements des conditions climatiques du passé et de l'air piégé qui témoignent des compositions atmosphériques du passé, notamment des gaz à effet de serre. Cette archive nous permet de décrypter le rôle des gaz à effet de serre dans les variations climatiques pendant huit cycles glaciaire-interglaciaires, soit l'équivalent de plus de 800 000 ans. Les carottes de glace, comme toute archive paléoclimatique, sont caractérisées par des incertitudes liées aux processus qui traduisent les variables climatiques en proxy, ainsi que par des incertitudes dues aux chronologies de la glace et des bulles d'air piégées. Nous développons un cadre méthodologique, basé sur la modélisation inverse dite Bayesienne et l'évaluation de fonctions complexes de densité de probabilité, pour traiter les incertitudes liées aux enregistrements paléoclimatiques des carottes de glace de manière précise. Nous proposons deux études dans ce cadre. Pour la première étude, nous identifions les probabilités de localisation des points de changement de pente de l'enregistrement du CO2 dans la carotte de WAIS Divide et d'un stack d'enregistrements de paléotempérature a partir de cinq carottes Antarctiques avec des fonctions linéaires par morceaux. Nous identifions aussi les probabilités pour chaque enregistrement individuel de température. Cela nous permet d'examiner les changements de pente à l'échelle millénaire dans chacune des séries, et de calculer les déphasages entre les changements cohérents. Nous trouvons que le déphasage entre la température en Antarctique et le CO2 à probablement varié (en restant inferieur, generalement, à 500 ans) lors de la déglaciation. L'âge des changements de temperature varie probablement entre les sites de carottage aussi. Ce résultat indique que les mécanismes qui reliaient la température en Antarctique et le CO2 lors de la déglaciation pouvaient être differents temporellement et spatialement. Dans la deuxième étude nous développons une méthode Bayesienne pour la synchronisation des carottes de glace dans le modèle inverse chronologique IceChrono. Nos simulations indiquent que cette méthode est capable de synchroniser des séries de CH4 avec précision, tout en prenant en compte des observations chronologiques externes et de l'information à priori sur les caractéristiques glaciologiques aux sites de forage. La méthode est continue et objective, apportant de la précision à la synchronisation des carottes de glace. / Deep polar ice cores contain records of both past climate and trapped air that reflects past atmospheric compositions, notably of greenhouse gases. This record allows us to investigate the role of greenhouse gases in climate variations over eight glacial-interglacial cycles. The ice core record, like all paleoclimate records, contains uncertainties associated both with the relationships between proxies and climate variables, and with the chronologies of the records contained in the ice and trapped air bubbles. In this thesis, we develop a framework, based on Bayesian inverse modeling and the evaluation of complex probability densities, to accurately treat uncertainty in the ice core paleoclimate record. Using this framework, we develop two studies, the first about Antarctic Temperature and CO2 during the last deglaciation, and the second developing a Bayesian synchronization method for ice cores. In the first study, we use inverse modeling to identify the probabilities of piecewise linear fits to CO2 and a stack of Antarctic Temperature records from five ice cores, along with the individual temperature records from each core, over the last deglacial warming, known as Termination 1. Using the nodes, or change points in the piecewise linear fits accepted during the stochastic sampling of the posterior probability density, we discuss the timings of millenial-scale changes in trend in the series, and calculate the phasings between coherent changes. We find that the phasing between Antarctic Temperature and CO2 likely varied, though the response times remain within a range of ~500 years from synchrony, both between events during the deglaciation and accross the individual ice core records. This result indicates both regional-scale complexity and modulations or variations in the mechanisms linking Antarctic temperature and CO2 accross the deglaciation. In the second study, we develop a Bayesian method to synchronize ice cores using corresponding time series in the IceChrono inverse chronological model. Tests show that this method is able to accurately synchronize CH4 series, and is capable of including external chronological observations and prior information about the glaciological characteristics at the coring site. The method is continuous and objective, bringing a new degree of accuracy and precision to the use of synchronization in ice core chronologies.

Un carottier à triple paroi a permis la récupération d'échantillons de calibre PQ3 et ZQ3, dans le till glaciaire des fondations du barrage Jules-Allard, dans le sud du Québec. Des essais cross-hole ont permis d'obtenir des modules de Young variant entre 1300 et 2600 MPa. Les essais triaxiaux consolidé non drainés de chargement monotones ont permis l'obtention de modules de Young Emax variant entre 290 et 6000 MPa, notamment en fonction de la fraction de gravier de l'échantillon. Les modules de Young des échantillons de calibre ZQ3 étaient similaires à ceux obtenus avec les essais cross-hole et leurs coefficients de Poisson étaient régulièrement nuls. Lors d'essais triaxiaux de chargement cycliques, les échantillons ont subi une contraction, avant de dilater à l'atteinte d'un point de transformation de phase, lorsque les contraintes ont croisé le cheminement de contraintes d'un échantillon chargé de façon monotone, à la même contrainte effective de consolidation. Mots clés: transformation de phase, comportement cyclique, till glaciaire, module de Young, vitesse de propagation d'onde de cisaillement, coefficient de Poisson. / A triple tube wireline core barrel allowed to retrieve PQ3 and ZQ3 calibre samples from the foundation glacial till of the Jules-Allard dam, South of Quebec. Cross-hole testing allowed to obtain Young's modulus ranging from 1300 to 2600 MPa. Consolidated-undrained triaxial monotonic testing, yielded Young's modulus Emax ranging between 290 MPa and 6000 MPa, notably as a function of the amount of gravel particles present in the samples. ZQ3 samples Young modulus were similar to those obtained using the cross-hole testing. In addition, the Poisson's ratio of these samples were regularly zero. During cyclic triaxial testing, contraction was followed by dilation when the samples reached a phase transformation point and the stress path intersected the stress path of a monotonically loaded sample having the same effective consolidation stress. keywords: phase transformation, cyclic behaviour, glacial till, Young's modulus, shear wave velocity, Poisson's ratio.

La dernière période glaciaire a été ponctuée de variations climatiques abruptes, appelées événements d'Heinrich, qui sont caractérisées par des épisodes de vêlages massifs d'icebergs et d'apports en eaux de fonte dans l'Atlantique Nord. Les événements d'Heinrich (EH) ont fait l'objet de nombreuses publications scientifiques portant sur la paléocéanographie de la zone Nord Atlantique. Cependant, peu d'études se sont intéressées aux modifications sédimentaires intervenant entre la source et la zone de dépôt. L'objectif de cette recherche est de caractériser l'enregistrement sédimentaire des événements d'Heinrich Hl et H2 le long de la marge du Labrador. Pour cela, une approche multi-proxy à haute résolution (propriétés physiques, géochimiques et datations au ¹⁴C) a été réalisée sur des carottes prélevées le long de la marge du Labrador en 2005 et 2006 à partir de l'embouchure du détroit d'Hudson jusqu'à l'éperon d'Hamilton. En combinant les différentes analyses, il a été possible d'illustrer l'évolution des dépôts de H1-H2 en fonction de la distance à la source. Les résultats ont permis de déterminer la répartition spatiale de Hl et H2 le long de la marge et les différents processus de déposition qui y sont associés. Les brusques variations de concentrations géochimiques, la faible susceptibilité magnétique et la présence de carbonates inorganiques des dépôts étudiés indiquent une provenance de l'érosion des roches carbonatées du plancher du détroit d'Hudson ainsi qu'une dilution des minéraux ferrimagnétiques par les carbonates (Andrews et Tedesco, 1992). La présence d'une couche noire à la base de H2 composée de sédiments qui proviendrait du détroit de Cumberland (Jennings et al., 1996) indique un transport sédimentaire précurseur à la déposition des débris carbonates du détroit d'Hudson. Les concentrations d'IRD ont permis de détecter des fluctuations d'apports en eaux de fonte et en icebergs dans Hl et H2. La variabilité des concentrations géochimiques de ces couches d'Heinrich s'explique par des apports sédimentaires de sources différentes ou par des réactions chimiques dues aux échanges eau-particules selon la localisation le long de la marge.