La caractérisation des transferts hydrologiques entre le domaine continental souterrain et la zone côtière est essentielle afin de mieux contraindre les bilans hydrologiques et de matière dans ces zones. Ce travail porte sur l'étude du fonctionnement hydrologique de la zone côtière sud de la lagune de Venise, et plus précisément sur les interactions entre les eaux de surface et les eaux souterraines des aquifères superficiels côtiers, où la charge hydraulique et le niveau du sol sont en dessous du niveau marin. <br />Les éléments majeurs et les isotopes stables de l'eau sont utilisés pour caractériser les processus de mélange au sein de l'aquifère de surface (0-30m). Nos résultats mettent en évidence la restriction de la salinisation de l'aquifère à une étroite bande côtière (<1km), en bon accord avec les résultats des travaux de modélisation hydrodynamique. En amont, le système de canaux constitue une barrière efficace de l'intrusion saline, malgré la subsidence de la zone, et le système d'irrigation/drainage permet un rafraîchissement permanent.<br />Les apports d'eau souterraine dans la lagune sont quantifiés en utilisant les traceurs des isotopes du radium et du radon. Leurs activités mesurées dans la lagune sont en excès par rapport à la mer et nécessitent un apport d'eau souterraine enrichie en ces traceurs, ce qui est confirmé par une approche de modélisation géochimique. Un bilan de masse à l'état stationnaire permet de calculer un flux variant entre 8.1x10^5 et 1.5x10^7 m3.j-1, dépendant fortement de la composition du pôle d'eau souterraine choisi. Ce flux, principalement dû à la recirculation d'eau lagunaire au sein des sédiments, est 2 à 10 fois plus élevé que le flux de surface, et représente potentiellement une source significative d'éléments dans les eaux côtières.

Les milieux poreux fracturés sont des milieux composés d'une roche présentant des zones de fracturation. Géologiquement, ces milieux sont caractérisés par la présence de fractures sur plusieurs échelles avec des propriétés hétérogènes créant des zones fortement perméables à comparer de la roche environnante. Hydrauliquement, ces milieux sont caractérisés par des temps de réponses courts correspondant aux sollicitations des structures fortement perméables et des temps de réponses longs correspondant aux sollicitations des structures faiblement perméables. Ce type de milieux est impliqué dans de nombreux enjeux sociétaux tels que l'étude de sites contaminés, le stockage d'éléments et l'exploitation de ressources. Le principal défi de la modélisation des milieux poreux fracturés réside en la représentation des hétérogénéités géométriques et physiques caractérisant ces milieux. Une représentation exacte du milieu naturel étant impossible, il s'agit de déterminer quelles sont les propriétés caractéristiques de ces milieux, c'est-à-dire les propriétés principales "responsables" de leur comportement. Mes travaux de thèse ont consisté à étudier l'impact des propriétés géométriques et physiques des fractures et de la matrice de l'échelle locale de la fracture et du bloc matriciel à l'échelle globale du réseau de fractures. Une première partie du travail correspond à la mise au point de méthodes d'évaluation des effets des structures et de leurs propriétés sur les échanges et une seconde partie en l'exploitation de ces méthodes. Au final, je propose un modèle double-porosité discret prenant en compte les propriétés du milieu identifiées comme caractéristiques de son comportement.

L'inversion des formes d'ondes (FWI) est une méthode d'imagerie sismique quantitative multiparamètre actuellement en plein essor. La prise en compte de l'ensemble des phénomènes de propagation enregistrés, en particulier des ondes de surface, en font une technique prometteuse pour des applications de la subsurface sur des problématiques géotechniques ou environnementales. Cependant, l'adaptation au contexte de la subsurface d'une méthode sismique validée numériquement ou sur des applications d'imagerie profonde est une tâche difficile. Nous proposons ici une approche par modélisation physique à échelle réduite à l'aide d'un laboratoire de Mesures laser-Ultrasonore Sans Contact (MUSC) mis en place dans le cadre de cette thèse. Les différents éléments du laboratoire ont été caractérisés pour la modélisation sismique à échelle réduite et une attention particulière a été portée à la simulation de la source. Les données acquises ont été validées par confrontation avec des données obtenues par modélisation numérique à l'aide des codes fournis par Géoazur dans le cadre du projet ANR SEISCOPE. Dans un deuxième temps, le potentiel de la FWI pour des applications de la subsurface a été évalué par l'approche conjointe de la modélisation numérique et de la modélisation expérimentale. Nous avons ainsi pu montrer par cette double approche que la FWI permet d'obtenir une information sur la morphologie de cavités souterraines grâce à la prise en compte des ondes de surface. Par ailleurs, sur les différents modèles multicouches étudiés, l'information sur les zones les plus superficielles fournie par les ondes de surface dispersives a pu être exploitée par la FWI.

La connaissance des processus de mobilisation et de migration de l’uranium dans les milieux souterrains est un enjeu majeur pour évaluer les risques environnementaux associés aux activités anciennes ou actuelles d'installations nucléaires en vue d'anticiper et gérer tout impact environnemental. La mobilisation et la migration sont habituellement décrites comme fortement contrôlées par les réactions d'adsorption, la présence de ligands, notamment sous forme colloïdale, étant susceptible de l'augmenter considérablement. La capacité des entités colloïdales à transporter l'uranium sur de longues distances va dépendre de leur origine, de leur composition, de leur structure, de leur stabilité et de leur réactivité, ainsi que de leur capacité à complexer l'uranium de manière forte et irréversible. Cependant, les faibles concentrations colloïdales dans les eaux souterraines s’avèrent parfois être un verrou analytique à la détection et à la caractérisation de la mobilisation et du transport colloïdaux de polluant(s) dans les eaux souterraines. Les protocoles d'extraction des phases colloïdales depuis les sols ou les méthodes destinées à concentrer les colloïdes depuis un échantillon d'eau souterraine ne sont pas toujours optimisés pour à la fois assurer un rendement d'extraction satisfaisant et respecter l'intégrité physiques et chimiques des entités colloïdales.Dans ce contexte, deux axes de recherche ont été définis dans ce travail de thèse afin 1) d’appréhender l’impact de l’utilisation de méthodes de préconcentration sur la phase colloïdale et 2) d’améliorer notre compréhension de la dynamique de mobilisation de l’uranium au sein d’un sol podzolisé et ainsi améliorer les modèles prédictifs du transport colloïdal de l’uranium dans les eaux souterraines.Les méthodes de séparation membranaires permettent de préconcentrer la phase colloïdale des eaux souterraines tout en assurant un recouvrement supérieur à 80%. Elles engendrent néanmoins des modifications des paramètres intrinsèques des entités colloïdales avec en particulier un décalage vers les grandes tailles de la distribution en taille de ces entités ainsi qu’une modification de la distribution des espèces entre phases dissoute et colloïdale. Ces modifications peuvent mener à une mauvaise appréciation non seulement des concentrations colloïdales mais aussi des mécanismes de complexation/sorption impliquant la surface des entités colloïdales et de leur stabilité dans les eaux souterraines. La matière organique constitutive des eaux du sol d’étude a été identifiée comme étant le vecteur principal de la remobilisation de l’uranium dans les eaux de sol, sans toutefois avoir un contrôle total de ce processus. Elle est constituée majoritairement d’acides fulviques et humiques de petite taille qui tendent à s’agglomérer au cours du temps, lors de la lixiviation du sol. La composition moléculaire de la matière organique varie de manière temporelle à la fois à l’échelle de l’horizon superficiel mais aussi le long du profil de sol. Divers processus (réduction, dénitrification…) interviennent durant la migration de la matière organique vers les horizons plus profonds et pourraient modifier ainsi sa nature. A l’échelle de l’horizon superficiel, les mécanismes de mobilisation diffèrent selon l’élément considéré. L’uranium en particulier est mobilisé rapidement mais de manière limitée (<2% du sol considéré), tandis que la matière organique est générée en continu et de manière pseudo-linéaire. / Knowledge of the processes of mobilization and migration of uranium in underground environments is a major challenge to assess the environmental risks associated with ancient or current activities of nuclear facilities in order to anticipate and manage any environmental impact. Mobilization and migration are usually described as strongly controlled by adsorption reactions, the presence of ligands, especially in colloidal form, being able to increase considerably. The ability of colloidal entities to transport uranium over long distances will depend on their origin, composition, structure, stability and reactivity, as well as their ability to complex uranium in a significant and irreversible manner. However, low colloidal concentrations in groundwater sometimes prove to be an analytical lock on the detection and characterization of colloidal mobilization and transport of pollutant(s) in groundwater. Protocols for extracting colloidal phases from soils or methods for concentrating colloids from a groundwater sample are not always optimized to both ensure satisfactory extraction yield and respect physical and chemical integrity colloidal entities.In this context, two lines of research have been defined in this thesis work to 1) to understand the impact of the use of preconcentration methods on the colloidal phase and 2) to improve our understanding of mobilization dynamics of uranium within a podzolized soil and thus improve predictive models of colloidal uranium transport in groundwater.The membrane separation methods make it possible to preconcentrate the colloidal phase of the groundwater while ensuring a recovery greater than 80%. Nevertheless, they cause modifications of the intrinsic parameters of the colloidal entities with, in particular, a shift towards the large sizes of the size distribution of these entities as well as a modification of the species distribution between dissolved and colloidal phases. These modifications may lead to an inaccurate appreciation not only of colloidal concentrations but also of complexation/sorption mechanisms involving the surface of colloidal entities and their stability in groundwater. The constituent organic matter of the study soil waters has been identified as the main vector for the remobilization of uranium in the soil water, without having a complete control of this process. It consists mainly of small fulvic and humic acids that tend to agglomerate over time, during leaching of the soil. The molecular composition of organic matter varies temporally both at the level of the superficial horizon and also along the soil profile. Various processes (reduction, denitrification ...) occur during the migration of organic matter to deeper horizons and could thus change its nature. At the level of the superficial horizon, mobilization mechanisms differ according to the element considered. Uranium in particular is mobilized quickly but in a limited way (<2% of the considered soil), while organic matter is generated continuously and in a pseudo-linear way.

La caractérisation de l'agencement spatial des propriétés hydrauliques est essentielle pour prédire les écoulements et le transport des solutés dans les milieux hétérogènes. Les méthodes de tomographie hydraulique, principalement développées pour estimer les propriétés des milieux poreux, n'ont qu'une faible résolution spatiale qui ne reflète pas la vraie hétérogénéité des distributions de fractures des milieux fracturés. Le principal objectif de cette thèse est de développer une nouvelle méthode d'inversion spécifique pour imager les propriétés hydrauliques et de transport des milieux fracturés à l'échelle du site. Pour atteindre ces objectifs, des expériences in situ ainsi qu'une nouvelle approche de modélisation inverse sont proposées, notamment en utilisant la température comme marqueur des écoulements. Nous proposons tout d'abord la tomographie d'écoulement basée sur des tests séquentiels de débimétrie entre puits, comme une nouvelle approche pour caractériser la connectivité des fractures ainsi que leur transmissivité. À partir de simulations numériques reproduisant des cas d'études synthétiques, nous montrons que l'approche par tomographie réduit significativement l'incertitude sur les paramètres estimés, et fournit une caractérisation détaillée du réseau de fracture sans requérir à l'utilisation d'obturateurs hydrauliques. Nous montrons ensuite comment les mesures de température peuvent être utilisées pour quantifier les écoulements dans les milieux fracturés. Le grand intérêt d'utiliser la température est d'obtenir facilement et de façon continue en puits des profils de température. En utilisant un modèle numérique d'écoulement et de transfert de chaleur à l'échelle du puits, une méthode d'inversion pour estimer les vitesses d'écoulement dans le puits à partir des données de température est proposée. Nous couplons ensuite les deux approches présentées précédemment dans une nouvelle approche expérimentale consistant en des enregistrements séquentiels de température dans un puits dans des conditions de pompage entre puits. L'application de cette approche de tomographie en température sur le site de Stang er Brune montre des résultats encourageants pour l'identification du réseau global de connectivité et des zones d'écoulement principales. Enfin, nous discutons de l'intérêt d'utiliser la chaleur comme traceur par rapport à l'utilisation de traceurs classiques. Nous montrons que réaliser des tests de traçage thermiques en milieu fracturé fournit des contraintes supplémentaires importantes sur les propriétés de transport du milieu. / The accurate characterization of distribution of hydraulic properties and connectivity distribution is essential to predict flow and transport in fractured media. Classical approaches were developed for homogeneous aquifers and result in smooth tomograms that often do not match true heterogeneity distribution of fractured media. The main goal of this thesis is to develop new inverse approaches specifically for imaging hydraulic and transport properties in fractured media at the field-scale. To attain this objective new in situ measurement methods as well as new inverse modelling frameworks are proposed.We first propose flow tomography (i.e., sequential cross-borehole flowmeter tests) as a new approach for characterizing fracture connectivity and transmissivities. Based on a discrete fracture network approach, we present a general method to invert flow tomography data. From synthetic case studies, we show that the tomographic approach reduces significantly the uncertainty on the parameter estimation. Flow tomography approach provides detailed characterization on fracture networks without the necessity of using packers. We then study the contribution of temperature measurements for quantifying flow in fractured media. The advantage of using temperature data is that temperature profiles can be obtained more easily and continuously in space, compared to flowmeter profiles. Using a numerical model of flow and heat transfer at the borehole scale, a method to invert temperature measurements to derive borehole flow velocities was proposed. We then couple the two previously proposed approaches in a new experimental approach which we call temperature tomography. This experiment consists of sequential borehole temperature logging under cross-borehole flow conditions. The full inverse framework is then presented to interpret temperature tomography experiments. Application of the temperature tomography approach to Stang er Brune field site showed encouraging results for the identification of general connectivity patterns and transmissivities of the main flowpaths. Finally, we explore the interest of using push-pull thermal tracer tests. Through field experiments and numerical modelling, we demonstrate that conducting push-pull heat tracer tests provide important constraints on the effective transport behavior.

Les aquifères fracturés hétérogènes se développent dans les roches cristallines (schiste, granite,...). Ces milieux présentent une large gamme de porosité de celle contenue dans les fissures et fractures qui constituent la zone de transfert principale de l'eau, à une porosité plus faible constituée de l'ensemble des microfissures et micropores de la matrice rocheuse. Une meilleure caractérisation de ces aquifères, appliquée plus particulièrement à la problématique de recherche de nouvelles ressources en eau en Bretagne, est réalisée au travers de deux axes principaux : (1) de suivis hydrochimiques, d'analyses isotopiques, d'expérimentations en laboratoire et (2) d'analyse des chlorofluorocarbones – CFC – (gaz anthropiques) en vue de la détermination des temps de résidence. <br />L'hydrochimie de ces aquifères est fortement liée à l'hétérogénéité physique des roches avec une variabilité spatiale beaucoup plus grande que la variabilité temporelle dans ces bassins-versants. Dans ces milieux, le cycle de l'azote et du soufre sont intimement liés au travers de la dénitrification par oxydation des sulfures. La potentialité d'une dénitrification par oxydation du fer dans les roches cristallines a été mise en évidence au travers d'expérimentations en batch. L'interprétation des données concernant les isotopes du soufre et de l'oxygène des sulfates montrent une cohérence entre la localisation spatiale et le processus mis en jeu quel que soit le site étudié. L'ensemble de ces résultats met en avant l'existence dans la zone profonde de ces aquifères, de solutés d'origine ancienne conservés dans la microporosité et qui contribuent à la chimie des eaux circulant dans les fractures connectées à cette microporosité.<br />Ces aquifères peuvent être décomposés en trois compartiments principaux que sont les zones de recharge, la zone altérée et la zone profonde. Dans chacun de ces compartiments les temps de résidence sont homogènes à l'échelle de la Bretagne : eaux récentes dans les zones de recharge, eaux anciennes (>50 ans) dans la zone profonde et un temps de résidence moyen de 18 ± 2 ans dans la zone altérée. Ces données indiquent que la gestion de ces aquifères doit prendre en compte l'échelle du temps de résidence de l'eau dans la zone altérée, réservoir de soutien des cours d'eau (supérieur à la dizaine d'années).

La tomographie de résistivité électrique (ERT) est une méthode courante géophysique de terrain, souvent utilisée pour détecter et l'évolution suivre les panaches de polluants en zone saturée. L'ERT est cependant une méthode intégratrice dont la fiabilité des modèles est confrontée aux problèmes de non-unicité des solutions du problème inverse. Ces contraintes limitent l'interprétation des modèles à un aspect qualitatif de la distribution des contrastes de résistivité modélisés en 2D ou 3D, résultant du choix des paramètres d'inversion et de l'association de paramètres du milieu non identifiables à l'échelle du volume poreux. Cette thèse propose de tester la faisabilité de la méthode pour quantifier les paramètres de transport de polluants et de solutés miscibles au contact des eaux souterraines, ainsi que la sensibilité des paramètres d'inversion les plus influents sur la modélisation. Les tests expérimentaux sont réalisés en laboratoire sur des empillements 2D de billes de verre sphériques (de l'ordre de la 100ène de μm) dans un réservoir en plexiglas transparent. Deux réseaux verticaux de 21 électrodes sont disposés sur les bords latéraux du réservoir pour effectuer le suivi ERT du traceur salé (NaCl dissout) à partir de 210 points de mesures en dipôle-dipôle transverse acquis toutes les 5 minutes afin d'optimiser la résolution temporelle. Le dispositif est également disposé face à un panneau lumineux permettant de réalisé un suivi vidéo simultané du colorant. L'analyse vidéo révèle une propagation plus rapide du colorant sur les bords latéraux qui reste somme toute négligeable pour les débits à la pompe les plus faibles. En revanche les mesures ERT sont fortement perturbées par les effets résistant de la cellule plexiglas qui se répercutent sur les modèles. La normalisation des mesures de résistivité apparente à partir d'une série acquise à l'état initial permet de les atténuer fortement. La modélisation est particulièrement sensible au choix du maillage, aux normes appliqués (L1 ou L2) sur les données et les paramètres, et au facteur d'amortissement. Des valeurs trop élevées du facteur d'acceptance tendent à lisser les contrastes au niveau du front de dispersion et augmentent l'impact des effets des bords horizontaux sur D et . A l'inverse, une modélisation contrainte par de faibles valeurs du facteur d'acceptance donnent des résultats plus proches de l'analyse vidéo, mais produit des effets de bosses à l'avant et à l'arrière du front. La vitesse interstitielle u est indépendante du choix des paramètres d'inversion pour l'ERT. Pour les deux méthodes u est toujours inférieure au débit imposé par la pompe, dont le décalage est exprimé par le facteur retard Rf . Les effets de retard résultent de l'adsorption du Na+ sur les surfaces des billes de verre chargées négativement qui retarde le front de dispersion du suivi ERT. Pour le suivi vidéo, la taille importante de la molécule du colorant favorise son piégeage dans les zones où la perméabilité est plus faible, en plus d'une éventuelle affinité avec la surface solide. Les contrastes de conductivité et la stabilité de l'interface créent par la différence de densité entre les fluides testés ici n'ont pas d'influences significatives sur la dispersion qui est dominée par le débit imposé à la pompe. Les estimations du coefficient de dispersion D en fonction du nombre de Péclet sont cohérentes avec la courbe théorique de Bachmat (1968). Cependant la dispersivité augmente pour les vitesses d'écoulement les plus élevées. Les premières expérimentations de terrain réalisées en 2D sur des sables de Fontainebleau présentent l'avantage de s'affranchir des effets de bords inhérents au laboratoire. En revanche la recalibration des données normalisée par la loi de Archie est plus complexe puisqu'il est nécessaire de tenir compte de l'état de saturation de la résistivité des fluides initialement présents. De plus l'erreur importante sur les modèles ne permet pas de déduire une estimation fiable des paramatères de transport u, (ou D), et Rf .

Depuis les années 50, les contaminations des sols et des eaux souterraines par les solvants chlorés se sont multipliées. Ces composés organiques volatils sont encore très fréquemment employés dans l'industrie et la plupart d'entre eux sont considérés comme des produits potentiellement cancérigènes pour l'homme. Ce travail est particulièrement focalisé sur une stratégie innovante de gestion de sites pollués : l'Atténuation Naturelle. Plusieurs mécanismes régissent cette démarche et parmi eux, c'est la déchloration réductrice anaérobie directe qui semble être le processus de dégradation le plus significatif. L'objectif de ce travail a été d'estimer si l'étude de la déchloration réductrice à plusieurs échelles d'observation permettait de mieux caractériser l'Atténuation Naturelle sur un site pilote. Le site d'étude est un complexe industriel encore en activité, où une source de pollution uniquement constituée de PCE (tétrachloroéthylène) a été mise en évidence. De plus, diverses zones polluées par des hydrocarbures d'huiles minérales ont été également caractérisées. Les campagnes de suivi physico-chimique ont permis de montrer qu'à l'aval hydraulique de la source, des sous-produits de biodégradation du PCE, comme le cis-1,2-DCE (cis-1,2-Dichloroéthylène) et le CV (Chlorure de Vinyle), ont été détectés très localement sur le site. La modélisation approximative selon la solution de Domenico et al. (1987) a démontré que le panache de composés organochlorés dissous variait latéralement. Afin de confirmer l'existence de la déchloration réductrice directe sur plusieurs points du site, plusieurs réacteurs fermés ont été préparés à partir de sous-sol et d'eaux souterraines. Les résultats ont indiqué que la biodégradation complète du PCE jusqu'à l'Ethylène est très localisée sur le site. L'activité bactérienne déchloratrice semble être fortement limitée par la quantité de matière organique dissoute, ce qui laisse suggérer un manque de donneurs d'électrons à l'échelle réelle. Les expérimentations effectuées dans une colonne de laboratoire ont fourni des informations intéressantes quant à l'influence des sulfates et de la quantité de matière organique sur la biodégradation du PCE. Les ajustements des constantes de dégradation du premier ordre sous Phreeqc-2.0 ont montré que la réduction des sulfates et la conversion complète du PCE jusqu'au cis-1,2-DCE pouvaient être simultanées dans certaines conditions physico-chimiques. Les différentes simulations exécutées avec le modèle d'écoulement, de transport et de biodégradation, selon une campagne précise, ont prouvé que le panache de pollution organochlorée est de faible largeur ; les écoulements étant rapides dans cette région, ce panache se déplace assez rapidement. La modélisation de la biodégradation s'est révélée délicate en raison des difficultés d'ajustement des concentrations en solvants chlorés mesurées, qui sont très hétérogènes spatialement et temporellement. Les constantes de dégradation estimées à l'échelle réelle semblent être relativement faibles par rapport aux données mesurées en laboratoire, du fait de l'hétérogénéité des conditions physico-chimiques modélisées à l'échelle réelle.

Les ressources en eau se définissent autant par la qualité de l'eau disponible que par sa quantité, et leur évolution future, sous l'influence des pressions anthropiques (pollutions, prélèvements, etc.) ou du changement climatique, est un enjeu majeur pour les sociétés. Ce mémoire fait la synthèse des travaux que j'ai menés dans ce cadre, en commençant par ceux qui portent sur la compréhension et la modélisation du fonctionnement des hydrosystèmes dans leur composante " naturelle ", si tant est que l'on puisse l'isoler en conditions récentes. Il est question des modèles hydrologiques TOPMODEL et CLSM, des incertitudes que révèlent leur confrontation à des observations issues de sites contrastés, et à certains processus couplés à l'hydrologie, comme la dénitrification dans les zones humides, l'hydraulique des cours d'eau et leur température. Une deuxième partie présente ensuite comment je me suis appuyée sur les modèles ainsi développés et évalués pour évaluer l'impact des changement globaux sur les hydrosystèmes, en insistant sur ceux du changement climatique. J'aborde enfin les perspectives de ces travaux, ciblées sur les influences réciproques entre surfaces continentales et climat, avec une attention particulière envers les stress hydriques à l'évapotranspiration et le rôle des eaux souterraines.

Pour d'évidentes raisons agronomiques et environnementales, être en mesure de concevoir et mettre en place des systèmes de culture dans lesquels les plantes accèdent aux ressources de manière optimale revêt une importance cruciale pour tous les intervenants impliqués dans la production agricole. Les techniques d'intensification telles que la mise en place de cultures d'inter-rang et l'agro-foresterie visent à accroître la productivité globale des terres tout en assurant la durabilité des agro-écosystèmes, via une optimisation de l'utilisation des ressources environnementales (lumière, eau et nutriments) par les plantes et une préservation des cycles géochimiques. En théorie, les moyens d'atteindre ces objectifs sont nombreux mais en pratique, les interactions souterraines sont complexes et difficiles à mesurer, de sorte que les progrès réalisés dans la conception d'agro-écosystèmes améliorés et durables demeurent modestes. Dans ce contexte, l'objectif de ce travail a été d'évaluer, au travers de mesures racinaires detaillées en rhizotron et au champ, les effets des cultures d'inter-rang sur la croissance des jeunes hévéas. La dynamique des interactions souterraines a été étudiée, tant au niveau de la racine individuelle qu'à celui du système racinaire entier, sur la base, notamment, de descriptions numériques détaillées. Une telle approche a permis de proposer, en outre, une voie novatrice pour l'analyse de la dynamique racinaire à l'echelle du systeme racinaire entier, à savoir l'analyse des trajectoires de croissance. Dans le cas de l'association maïs-hévéa, les expérimentations en rhizotron ont permis de mettre en evidence que les interactions souterraines entre ces deux plantes peuvent induire des modifications de la croissance de leurs racines, à la fois à l'échelle de la racine individuelle et à celle du système racinaire entier. Toutefois, une telle coordination des dynamiques racinaires des plantes associées n'a pas pu être confirmée dans le cas des traitements manioc-hévéa et arachide- hévéa. Les expérimentations au terrain ont fourni, de manière assez prévisible, une image complexe des interactions souterraines entre hévéa et cultures d'inter-rang. Toutefois, un premier résultat obtenu par le biais de la mise en place de 'pièges à racines' dans un traitement niébé-hévéa, a été de montrer que ces deux plantes n'avaient pas un comportement compétitif marqué l'une vis-à-vis de l'autre. Il est également apparu que les hévéas paraissent 'investir' dans des racines 'coûteuses', car de faible longueur spécifique, probablement pour assurer une certaine durabilité de ces organes, tandis que les cultures d'inter-rang favorisent l'allocation des assimilâts vers des racines de longueur spécifique élevée, de construction moins 'coûteuses', probablement en réponse a un impératif de croissance plus rapide (suggéré par les taux d'élongation racinaire mesurés au cours des expérimentations en rhizotron). Enfin, excepté le cas du manioc, l'introduction de cultures d'inter-rang telles que le maïs et l'arachide n'a pas eu d'impact significatif sur le développement des jeunes hévéas, comme en attestent l'évolution de leur circonférence, hauteur et développement foliaire. Ce résultat de terrain est compatible avec les résultats des expérimentations en rhizotron qui n'ont démontré aucun effet inhibiteur des cultures d'inter-rang sur le développement de la partie aérienne des hévéas. Bien que les travaux présentés dans ce rapport, ne permettent pas, à eux seuls de conclure de manière définitive sur la façon dont les espèces cultivées en association peuvent se compléter mutuellement sur le plan fonctionnel, ils apportent des éléments de réponse préliminaires à cette question complexe ainsi que des méthodes permettant de les obtenir. Au total, ce travail représente donc une contribution à la conception des agro-écosystèmes durables qui deviennent de plus en plus indispensables dans le contexte d'une demande mondiale croissante en produits alimentaires et matières premières. En outre, certains des résultats obtenus dans le cadre de cette thèse ouvrent des perspectives pour des recherches plus approfondies, avec une finalité agronomique appliquée