Le changement climatique et les événements climatiques sévères tels que les périodes de sécheresse/humidification prolongées sont à l'origine du phénomène de retrait-gonflement dans les sols argileux. Ce phénomène est affecté par les interactions sol-végétation-atmosphère (SVA) et peut causer d’importants dommages structurels aux constructions légères telles que les bâtiments résidentiels. L’objectif de ce travail de recherche est de modéliser le comportement in situ du retrait-gonflement des sols gonflants dans un contexte SVA en se basent sur des outils numériques. Une méthode d'interaction sol-atmosphère est initialement présentée accompagnée d’un modèle couplé hydro-thermique du sol. Cette approche a été principalement mise en place afin de déterminer les conditions aux limites temporelles à la surface du sol en se basent sur la notion du bilan de masse et d'énergie pour déterminer a posteriori, les modifications spatio-temporelles de la succion du sol, de la teneur en eau et de la température. Cette approche a été validée à l'aide des observations in situ des sites instrumentés. Par la suite, l’influence de l’absorption d’eau par les végétations a été intégrée dans le terme source de l’écoulement de l’eau dans un milieu non saturé, à l’aide d’un modèle d’absorption d’eau de racine existant.Les variations temporelles de succion ont été postérieurement reliées au comportement volumique du sol en appliquant une approche simple développée à partir des résultats expérimentaux des essais de séchage/humidification réalisés dans la littérature. Les indices associés dans le plan indice des vides-log succion, ainsi que les paramètres complémentaires du modèle linéaire ont été corrélés aux paramètres géotechniques de base. L'approche proposée a été ultérieurement validée avec des données in situ fournies par la surveillance d’un site expérimental. Le site expérimental de Roaillan a été instrumenté afin de surveiller les modifications physiques du sol ainsi que le comportement structurel du bâtiment. Les comparaisons entre les résultats de la modélisation et les observations in situ de la succion du sol, la teneur en eau, la température et les mouvements du sol dans le temps ont montré une performance acceptable du modèle. L’approche a ensuite été appliquée pour étudier l’influence des projections climatiques futures (2050) sur les variables physiques et les mouvements du sol sur ce site. Trois scénarios RCP relatifs aux changements climatiques ont été examinés dans cette étude, qui ont révélé des différents comportements possibles à court terme et à long terme. Finalement, l'approche développée a été appliquée au territoire français en le divisant en six régions climatiques. Différents paramètres de sol ont été attribués à chacune de ces régions climatiques afin de définir les conditions de référence. En conséquence, l’influence de différents facteurs externes sur les mouvements du sol a été analysée sur une période donnée. Enfin, l’étude suggère les mesures adéquates à prendre pour minimiser l’amplitude du phénomène de retrait et de gonflement dans un contexte SVA. / Climate change and severe climatic events such as long drought/rehydration periods are at the origin of the shrinkage and swelling phenomenon in expansive soils. This phenomenon is affected by Soil-Vegetation-Atmosphere (SVA) interactions and can cause severe structural damage to lightly loaded constructions such as residential buildings. The objective of this re-search work is to simulate the in-situ behavior of the shrinkage-swelling in expansive soils in a SVA context using numerical tools. A soil-atmosphere interaction method is primarily presented along with a coupled hydro-thermal soil model. This approach was established in order to determine primarily, the natural time variable boundary conditions at the considered soil surface based on the mass and energy balance concept, and secondly to determine the spatial-temporal changes of the soil suction, water content and temperature. This approach was validated using in situ observations of monitored sites. Thereafter, the influence of the water uptake by vegetation was incorporated in the source term of the unsaturated water flow theory, using an existing root water uptake model. Subsequently, the temporal variations of the soil suction were related to the volume change behavior using a simple approach developed based on the experimental results of drying/wetting tests performed in the literature. The associated volumetric indices in the void ratio-log suction plan, along with the complementary parameters of the linear model were correlated with basic geotechnical parameters. The proposed approach was validated with in situ data provided from an experimental site. The Roaillan experimental site was instrumented in order to monitor the soil’s physical changes along with the structural behavior of the building. Comparisons between the simulated and observed soil suction, soil water content, temperature and soil movements in time and depth showed an acceptable performance of the predictions. The approach was then extended to study the influence of future climate projections (2050) on the soil’s physical variables and movements. Three RCP climate change scenarios were considered in this analysis which revealed different possible behavior in both short term and long term. Finally, the developed approach was applied to the French territory by dividing it to six different climatic regions. Different soil parameters were attributed to each of these climatic regions in order to set the reference condition. Thereafter, the influence of different external factors was analyzed on the soil movements over a chosen period. The study finally suggests the adequate actions to take for minimizing the amplitude of the shrinkage and swelling phenome-non in a SVA context.

Les marais de l'Ouest ont été gagnés sur les sédiments fluvio-marins par poldérisation débutée dès le moyen âge. Les aménagements hydrauliques avaient, initialement, des objectifs sanitaires et d'élevage. A partir des endiguements, ils ont divisé les territoires en marais mouillés et marais desséchés. Depuis les années 1970 le développement des cultures intensives de céréales a nécessité un rabaissement complémentaire de la nappe par drainage en partie pour augmenter la désalinisation de surface et limiter l’engorgement des sols. Néanmoins, ces territoires restent caractérisés par des nappes proches de la surface. En conséquence les profils hydriques vont être gouvernés par les conditions météorologiques y compris l'évapotranspiration, la pluviométrie, mais également par les remontées capillaires issues d'une nappe salée. D'autre par la nature argileuse des sols et ses propriétés de retrait conditionnent énormément le fonctionnement hydrodynamique et les évolutions de structure.La première partie du travail a été de suivre les évolutions des profils hydriques et de salinité en parcelles non drainées (prairies) et en parcelles drainées (Maïs, Blé, Tournesol). Ces suivis ont été complétés par les mesures des niveaux de nappes et par des mesures tensiométriques via des bougies poreuses implantées à différentes profondeurs. L'objectif final a été de calculer et modéliser les profils de réserve utile (RU) et de RU "résiduelle" utilisable par les plantes. Dans ces systèmes alimentés par les remontées capillaires, le réseau racinaire puise l'eau dans la zone non saturée de surface (vadose) puis dans la zone saturée sous-jacente. La teneur en eau caractéristique de l'interface zone non saturée - zone saturée a été déterminée par analogie entre les chemins d'état de la matrice argileuse le long de sa courbe de retrait et du sol le long de sa courbe de compaction. Les profils de RU résiduelle utilisable par les plantes ont été calculés à partir des profils de teneurs en eau puis comparés aux profils de RU obtenus via les données de station météorologique. Ces profils de RU résiduelle ont pu être écrits sous forme d'équation polynomiale du second degré puis modélisés. Dans cet environnement alimenté par les remontées capillaires, ces profils de RU résiduelles peuvent être modélisés à partir d'un paramètre facilement mesurable en surface qui prend en compte la structure du sol et les conditions météorologique : soit la teneur en eau à 10 cm de profondeur. Cette modélisation reste suffisamment réaliste pour être utilisée comme un outil prédictif face à la pédodiversité et/ou les rendements de culture.A ce travail s'ajoute deux études préliminaires : - les mesures des conductivités thermiques effectives de ces sols par la méthode du fil chaud et leurs modélisations dans les systèmes biphasés : eau - argile et air - argile, mais également pour les systèmes triphasés non saturés : eau - air - argile. Les perspectives sont la modélisation des transferts thermiques et hydriques dans le sol à partir de la surface, - et l'élaboration d'un protocole d'imprégnation-polymérisation des sols argileux humides par des résines de type HEMA. Cette imprégnation permet d'envisager la confection de lames minces dans le matériau argileux induré avec conservation de sa structure initiale humide. Les perspectives sont la pétrographie quantitative à l'interface racine - sol le long de profils verticaux dans les environnements argileux à degrés de saturation et structure évolutives. / The coastal marshlands are territories generally reclaimed on primary fluvio-marine sediments. They result from hydraulic managements and/or polderization which may date from the Middle Ages. Historically these hydraulic managements were built for goals of wholesomeness, breeding and farming. They isolate two territories: the dried marshes and the wet marshes. For the intensive cereal crops the slow drying caused by land reclamation was recently improved by the drainage, in part for increase the depth of desalinization and decrease waterlogging. Nevertheless, these territories remain characterized by shallow ground water of initial salt water. Consequently, the hydric profiles are governed by the meteoric conditions including the Evapotranspiration, the rainfall, but also the capillarity rises from the salt groundwater. Moreover, the clay dominated nature of the soils and their drastic shrinkage properties govern the hydrodynamic functioning and the soil structure behavior.The first part of the work was the monitoring of the water content and salinity profiles in drained cereal crops and in undrained grasslands. These measurements have been completed by the ground water level and tensiometric monitoring. The final goal was the calculation and modeling of the available water capacity (AWC) and plant available water (PAW) profiles. In these systems mainly supplied by the capillarity rises, the root network gets water in the subsurface vadose zone and then in the deeper saturated groundwater zone. The water content characteristic of the interface between the vadose and saturated zone was determined by comparison between the clay material state paths along its shrinkage curve and along its compaction curve. The PAW profiles were calculated from the water content profiles and then compared to the AWC profiles. The PAW profiles have been equated as polynomial second degree equations. In these shallow groundwater environments the PAW profiles have been modeled taking into account an easy measurable surface parameter which includes the soil structure behavior and the meteoric conditions: i.e. the water content measured at 10 cm depth. The PAW modelling remains sufficiently realistic to be used as a tool for farming management. Two preliminary studies were added to this work: - the measurement of effective thermal conductivity of the clayey soils by the transient hot wire method, and the modeling of the effective thermal conductivity of biphasic air-clay and water-clay media, but also triphasic unsaturated air-water-clay media. The prospect is the modeling of thermal and hydric transfer from the surface to the depth. - and the elaboration of a protocol of impregnation - hardening for wet clay dominated soils by HEMA resins. This impregnation allows the making of thin sections in these clay materials with conservation of their initial wet structures. The prospective is the quantitative petrography at the root - clay matrix interface along vertical profiles in clayey soils at different degrees of saturation and different structures.