Les propriétés hydrauliques d'un till du Nord du Québec en condition saturée et non saturée ont été étudiées en laboratoire à l'aide de l'essai de conductivité hydraulique et l'essai de succion qui permet d'obtenir la courbe de rétention d'eau. Les essais ont été réalisés avec des échantillons compactés sous différentes conditions de teneur en eau et de masse volumique sèche, de manière à évaluer l'influence de ces conditions sur les propriétés hydrauliques. Les études précédentes ont montré qu'une structure agrégée se développe dans le till lorsqu'il est compacté de côté sec de l'optimum Proctor, c'est-à-dire à un degré de saturation de compactage inférieur au degré de saturation optimum (Sr0 < Sr-opt). Le till étudié dans le cadre de cette étude diffère de ceux utilisés dans les études précédentes par son très faible pourcentage de particules argileuses (< 2 μm) de 0,7%. L'objectif principal de cette étude est d'évaluer l'influence des conditions de compactage sur la conductivité hydraulique et la courbe de rétention d'eau et de déterminer si une structure agrégée se développe lorsque le till est compacté du côté sec de l'optimum. Les résultats expérimentaux typiques de la conductivité hydraulique saturée obtenus pour les échantillons compactés du côté humide de l'optimum (Sr0 > Sr-opt) avec un indice des vides de 0,328 varient de 2,6 x 10-7 à 8,6 x 10-7 m/s. Les conductivités hydrauliques typiques obtenues pour les échantillons compactés du côté sec de l'optimum (Sr0 < Sr-opt) avec le même indice des vides étaient très similaires, variant entre 2,4 x 10-7 et 9,1 x 10-7 m/s. Les courbes de rétention d'eau mesurées en laboratoire sont très semblables entre elles, présentant toutes une structure de sol homogène, peu importe si l'échantillon a été compacté du côté sec ou du côté humide. Il apparait que la conductivité hydraulique et la courbe de rétention d'eau de ce till sont très peu influencées par le degré de saturation de compactage. Ceci peut être expliqué par le très faible pourcentage de particules argileuses, qui prévient la formation d'agrégations de particules dans la structure du till lorsqu'il est compacté du côté sec de l'optimum Proctor. Plusieurs modèles d'estimation de la conductivité hydraulique sont utilisés. Cette étude propose trois nouveaux modèles d'estimation de la conductivité hydraulique saturée pour les tills du Nord du Québec. Le premier modèle met en relation la conductivité hydraulique saturée avec le diamètre des grains de la fraction fine (d50 FF) et le deuxième modèle est une équation modifiée de Kozeny-Carman et Hazen basée sur la porosité et une taille de particule effective (d10). Finalement, un modèle permettant d'évaluer l'effet de l'agrégation avec l'équation de Kozeny-Carman modifiée dans les tills compactés est proposé.

Les tills compactés sont largement utilisés comme noyau des barrages en remblai dans le Nord québécois. La mise en place en différentes couches successives laisse présager un rapport d’anisotropie (rK) de la perméabilité différent de l’unité. Différents montages expérimentaux sont disponibles dans la littérature pour évaluer l’anisotropie de la conductivité hydraulique sur plusieurs types de sols. Toutefois, ces méthodes présentent deux principaux désavantages : i) le rapport d’anisotropie est évalué sur deux échantillons distincts et ii) la procédure expérimentale peut potentiellement remanier l’échantillon par des manipulations hasardeuses. Dans l’objectif d’améliorer la caractérisation du rapport d’anisotropie des tills compactés, un nouveau perméamètre cubique a été conçu. Cette nouvelle cellule permet de quantifier le paramètre rK dans le domaine saturé ou quasi-saturé, et ce sans démanteler la cellule. Successivement, une percolation verticale et une horizontale sont effectuées sur un même échantillon. Une procédure spéciale a été développée afin de limiter l’écoulement dans les pierres poreuses inactives, i.e. parallèles à l’écoulement. Une solution de gélatine bovine qui est solide à une température inférieure à 30°C est utilisée afin de sceller ces pierres inactives. De l’eau chaude à 60°C permet de chasser cette gélatine sans démanteler la cellule. Un nouvel essai de perméabilité peut ainsi être effectué sur le même échantillon, mais selon une autre direction. Le nouveau perméamètre et la procédure ont permis d’évaluer l’effet de la mise en place et de la teneur en pierre sur rK. Pour une analyse dans le domaine quasi-saturé, le rapport d’anisotropie varie de 3 à 7 pour une mise en place du côté sec de la courbe Proctor. Pour le côté humide, le rapport d’anisotropie est près de l’unité. Cependant, la procédure développée pour des conditions saturées devrait être utilisée pour confirmer ces résultats. En conclusion, un outil intéressant a été conçu pour quantifier le rapport d’anisotropie des tills compactés, une problématique importante pour l’analyse des barrages en remblai. Mots clés: perméabilité, anisotropie, perméamètre cubique, tills compactés, barrage en remblai / Compacted tills are widely use as core of embankment dams in northern Quebec. It is expected that placement method in successive layers may result in a permeability anisotropy different from unity. Various experimental techniques are available to assess the permeability anisotropy for different types of materials. However, these methods have two principal disadvantages: i) the anisotropy rapport is evaluated on separate samples and ii) the experimental procedure could potentially disturb the sample by inappropriate manipulations. In order to improve the characterization of the anisotropy rapport of compacted tills, a new square cell was built. This new cell is able to quantify the anisotropy rapport in both states, saturated and almost saturated conditions. The permeability under vertical and horizontal flow is assessed on a single sample and a new experimental procedure has been developed to minimize preferential seepage paths along the porous stones. Injection of bovine gelatin is used to temporarily seal the porous stones perpendicular to the flow path (unused). Continuous hot water flow is used to extract the gelatin once the permeability testing is done. This new device and procedure have been successfully tested to evaluate the effect of placing method and stone content on the permeability anisotropy of compacted tills. For an analysis in quasi-saturated state, the anisotropy rapport is greater than 3 or 7 for a dry compaction. On the wet side, the anisotropy rapport is close to one. However, the procedure developed in the saturated condition should be used to confirm these results. In conclusion, this project has helped to build a useful tool to quantify the anisotropy rapport of tills, an important issue in the evaluation of embankment dams. Key words: permeability, anisotropy, square permeability cell, compacted tills, embankment dam.

No description available.

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.

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.

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.