La compréhension des mécanismes de transport de solutés dans un milieu poreux variablement saturé tel que les sols à nappe drainée par tuyaux enterrés, est essentielle afin de répondre aux interrogations concernant l'impact de ces systèmes sur la qualité des eaux de surface. L'approche choisie repose sur des expériences de traçage réalisées en régime permanent sur modèle physique de sol au laboratoire, permettant un contrôle des conditions initiales et aux limites. La première partie de la thèse est dédiée aux études monodimensionnelles, afin de caractériser la zone non saturée du sol. Nous avons utilisé une colonne de sol (1*0,15 m²) et un modèle physique de grande taille (1*2*0,5 m3), remplis tous deux d'un sol sablo-limoneux. Les traceurs utilisés sont le H218O et le Cl-. La deuxième partie du travail traite le cas des écoulements bidimensionnels, réalisées uniquement sur le modèle physique. Deux scénarios ont été choisis : un soluté apporté en surface et un soluté résidant dans le sol. Seul le Cl- a alors été utilisé comme traceur. Le transport de solutés a été étudié par un suivi in situ (sondes de conductivité électrique) ainsi qu'à l'exutoire du système. L'analyse des courbes d'élution obtenues pour les études monodimensionnelles a permis l'identification des mécanismes de transport et le calcul des paramètres de transport pertinents (par résolution inverse). Pour les expériences bidimensionnelles, il a été observé une hétérogénéité des champs de vitesses, avec une zone proche du drain où le lessivage est rapide et une zone proche de l'interdrain où le soluté progresse à faible vitesse. Dans le cas du soluté résidant dans le sol, les courbes d'élution (concentration et conductivité électrique) au drain montrent un palier stable puis une diminution de forme « exponentielle » jusqu'à atteindre les concentrations du niveau de base. La migration du soluté dans le cas d'un pulse apporté en surface du sol est singulier. La courbe d'élution en sortie du drain présente deux modes : (i) un premier mode sous forme d'un pic, (ii) et le deuxième sous forme d'un palier constant de concentration. Ces résultats expérimentaux ont été confrontés à la modélisation mécaniste (code HYDRUS-2D), utilisant les paramètres de transport obtenus dans les études monodimensionnelles. L'utilisation d'une deuxième méthode, le suivi de particules (approche lagrangienne), nous a permis d'obtenir le cheminement de l'eau et le temps de transfert des particules dans le sol.

La compréhension des mécanismes de transport de solutés dans un milieu poreux variablement saturé tel que les sols à nappe drainée par tuyaux enterrés, est essentielle afin de répondre aux interrogations concernant l'impact de ces systèmes sur la qualité des eaux de surface.<br />L'approche choisie repose sur des expériences de traçage réalisées en régime permanent sur modèle physique de sol au laboratoire, permettant un contrôle des conditions initiales et aux limites. La première partie de la thèse est dédiée aux études monodimensionnelles, afin de caractériser la zone non saturée du sol. Nous avons utilisé une colonne de sol (1*0,15 m²) et un modèle physique de grande taille (1*2*0,5 m3), remplis tous deux d'un sol sablo-limoneux. Les traceurs utilisés sont le H218O et le Cl-. La deuxième partie du travail traite le cas des écoulements bidimensionnels, réalisées uniquement sur le modèle physique. Deux scénarios ont été choisis : un soluté apporté en surface et un soluté résidant dans le sol. Seul le Cl- a alors été utilisé comme traceur. Le transport de solutés a été étudié par un suivi in situ (sondes de conductivité électrique) ainsi qu'à l'exutoire du système.<br />L'analyse des courbes d'élution obtenues pour les études monodimensionnelles a permis l'identification des mécanismes de transport et le calcul des paramètres de transport pertinents (par résolution inverse). Pour les expériences bidimensionnelles, il a été observé une hétérogénéité des champs de vitesses, avec une zone proche du drain où le lessivage est rapide et une zone proche de l'interdrain où le soluté progresse à faible vitesse. Dans le cas du soluté résidant dans le sol, les courbes d'élution (concentration et conductivité électrique) au drain montrent un palier stable puis une diminution de forme « exponentielle » jusqu'à atteindre les concentrations du niveau de base. La migration du soluté dans le cas d'un pulse apporté en surface du sol est singulier. La courbe d'élution en sortie du drain présente deux modes : (i) un premier mode sous forme d'un pic, (ii) et le deuxième sous forme d'un palier constant de concentration. Ces résultats expérimentaux ont été confrontés à la modélisation mécaniste (code HYDRUS-2D), utilisant les paramètres de transport obtenus dans les études monodimensionnelles. L'utilisation d'une deuxième méthode, le suivi de particules (approche lagrangienne), nous a permis d'obtenir le cheminement de l'eau et le temps de transfert des particules dans le sol.

La détermination des modalités de transfert de l'eau dans un aquifère et la connaissance des volumes d'eau qu'il contient , fait l'objet de nombreuses études en hydrologie. Pour obtenir ces informations, l'hydrogéologue cherche à représenter le transfert de l'eau par une fonction simple qui lui permettrait de prévoir la réponse du milieu souterrain aux stimuli d'entrée , tels que les pluies ou les pollutions, de connaître le temps de séjour moyen de l'eau dans la matrice et par suite de calculer les volumes d'eau contenus dans celle-ci . Dans cette étude, il est proposé de considérer un aquifère comme une "boite noire" dont on ne connait que les fonctions d'entrée et de sortie , et de lui appliquer les méthodes de l'analyse-système qui se sont développées en électronique. A partir de deux hypothèses , vérifiées sur un modèle physique de laboratoire, il est montré qu'un aquifère dont l'entrée est représentée par la surface du bassin hydrologique , peut être assimilée à un système linéaire dont la réponse impulsionnelle est une fonction exponentielle unlque. L'argument de cette exponentielle définit le temps de séjour moyen de l'eau dans le système et le produit du débit moyen des écoulements par ce temps moyen permet de calculer le volume d'eau contenu dans l'aquifère. L'application de cette méthode sur un bassin hydrologique situé à St Genis, Hautes Alpes, s'est avérée satisfaisante . Elle permet de définir une réponse impulsionnelle constituée de deux fonctions exponentielles représentatives de deux systèmes qui composent l'aquifère du bassin et d'en calculer les volumes utiles . L'utilisation des données isotopiques du milieu, nécessairement pondérées par les débits, apparaît comme un supplément d'information indispensable à la bonne détermination de la réponse impulsionnelle et à l'estimation des réserves de l'aquifère.