Since 1996, super-critical CO$_2$ has been injected at a rate of $\sim$0.85~Mt~yr$^{-1}$ into a pristine, saline aquifer at the Sleipner carbon capture and storage project. A suite of time-lapse, three-dimensional seismic reflection surveys have been acquired over the injection site. This suite includes a pre-injection survey acquired in 1994 and seven post-injection surveys acquired between 1999 and 2010. Nine consistently bright reflections within the reservoir, mapped on all post-injection surveys, are interpreted to be thin layers of CO$_2$ trapped beneath mudstone horizons. The areal extents of these CO$_2$ layers are observed to either increase or remain constant with time. However, volume flux of CO$_2$ into these layers has proven difficult to measure accurately. In addition, the complex planform of the shallowest layer, Layer 9, has proven challenging to explain using reservoir simulations. In this dissertation, the spatial distribution of CO$_2$ in Layer~9 is measured in three dimensions using a combination of seismic reflection amplitudes and changes in two-way travel time between time-lapse seismic reflection surveys. The CO$_2$ volume in this layer is shown to be growing at an increasing rate through time. To investigate CO$_2$ flow within Layer~9, a numerical gravity current model that accounts for topographic gradients is developed. This vertically-integrated model is computationally efficient, allowing it to be inverted to find reservoir properties that minimise differences between measured and modelled CO$_2$ distributions. The best-fitting reservoir permeability agrees with measured values from nearby wells. Rapid northward migration of CO$_2$ in Layer~9 is explained by a high permeability channel, inferred from spectral decomposition of the seismic reflection surveys. This numerical model is found to be capable of forecasting CO$_2$ flow by comparing models calibrated on early seismic reflection surveys to observed CO$_2$ distributions from later surveys. Numerical and analytical models are then used to assess the effect of the proximity of an impermeable base on the flow of a buoyant fluid, motivated by the variable thickness of the uppermost reservoir. Spatial gradients in the confinement of the reservoir are found to direct the flow of CO$_2$ when the current is of comparable thickness to the reservoir. Finally, CO$_2$ volume in the second shallowest layer, Layer~8, is measured using structural analysis and numerical modelling. CO$_2$ in Layer~8 is estimated to have reached the spill point of its structural trap by 2010. CO$_2$ flux into the upper two layers is now $\sim$40\% of total CO$_2$ flux injected at the base of the reservoir, and is increasing with time. This estimate is supported by observations of decreasing areal growth rate of the lower layers. The uppermost layers are therefore expected to contribute significantly to the total reservoir storage capacity in the future. CO$_2$ flow within Layer~9 beyond 2010 is forecast to be predominantly directed towards a topographic dome located $\sim$3~km north of the injection point. This dissertation shows that advances in determining the spatial distribution and flow of CO$_2$ in the sub-surface can be made by a combination of careful seismic interpretation and numerical flow modelling.

Seit mehr als 100 Jahren wird Oberflächenwasser mit verschiedenen Verfahren versickert, wobei sich dessen Qualität stark verbessert. Sie finden zunehmend weltweit Anwendung, um eine zuverlässige Versorgung mit sauberem Trinkwasser zu gewährleisten. Solche Verfahren erfordern nur einen geringen technischen Aufwand. Während bei der Versicherung mit Schluckbrunnen und bei künstlichen Infiltrationsbecken das Wasser noch gepumpt und oft auch vorgereinigt werden muss, reduziert die Nutzung eines natürlichen Gewässers den Aufwand nochmals. Am Beispiel eines Untersuchungsgebiets am Tegeler See in Berlin werden die hydraulischen Prozesse modelliert. Die regionale, instationäre und 3-dimensionale Modellierung eines langen Zeitraums zeigt, dass die bisher verwendeten linearen Ansätze zur Beschreibung der Durchlässigkeit der Kolmationsschicht sowohl die infiltrierten Wassermengen als auch die Infiltrationsprozesse nur unzureichend wiedergeben. Grundwasserpiegelschwankungen werden stärker als bisher angenommen gedämpft. Als Folge dieser Wasserspiegelschwankungen wird die Bodenluft in der ungesättigten Bodenzone ausgetauscht und Sauerstoff eingetragen. Auf diese Weise erhöht sich die Durchlässigkeit der Kolmationsschicht um mindestens eine Größenordnung. Auf Grundlage dieser Ergebnisse wurden eine instationäre Wasserbilanz aufgestellt und die Uferfiltratanteile bestimmt. Die genaue lokale Lage eines Grundwasserstauers wird mit Hilfe eines neu entwickelten inversen Modellkonzeptes bestimmt. Er wird mit räumlich verteilten Pilot Points und Überparameterisierung unter Nebenbedingungen kalibriert. Die Zielgröße sind Piezometerhöhenschwankungen. Kreuzvalidierung bestätigt, dass die Lage des Grundwasserstauers präzise bestimmt wird. Diese hat einen entscheidenden Einfluss auf das räumliche Strömungsfeld und liefert neue Anhaltspunkte zur Genese des Tegeler Sees. Ein Transportmodell liefert lediglich mit den hydraulisch kalibrierten Parametern eine genaue Anpassung der Durchbruchskurven von Chlorid, 18O und Temperatur. Dies kann als Vorhersage angesehen werden und zeigt die Güte des Modells. Es konnte gezeigt werden, dass der Sauerstoffeintrag die Redoxprozesse in der Infiltrationszone dominiert. In Kombination mit der Lage des Grundwasserstauers können so auch zuvor widersprüchliche und lokale geochemische Verhältnisse schlüssig interpretiert werden. Es wird eine Theorie entwickelt, die erklärt, wie die vertikale Schichtung von oxidiertem und reduziertem Wasser zu starker Verockerung der Brunnen führt. / Several kinds of managed aquifer recharge techniques provide very good purification of surface water since more than 100 years. In order to maintain a reliable supply of clean water, they are becoming increasingly popular all over the world. These methods require low technical effort. At Aquifer Storage and Recovery and ponded infiltration the recharged amounts are technically controlled. The infiltration water has to be pumped and often pretreated. At bank filtration this is dispensable, the approach, of using existing surface water bodies is even more consequent. Exemplarily, at a test site at Lake Tegel, Berlin, Germany, the hydraulic processes are modelled. By means of 3D long term regional and transient hydraulic modelling it was detected that the existing approaches for determining the leakance induce large errors in the water balance and describe the infiltration zone insufficiently. The leakance could be identified to be triggered by the groundwater table, causing air exchange and intrusion of atmospheric oxygen, which reduces clogging by altered redox conditions by at least one order of magnitude. This causes that changes of the groundwater table are mitigated much more than previously assumed. Taking these findings into account, a transient water balance is determined and bank filtration ratios are quantified. A new inverse modelling concept has been developed and applied to a 3D short term local and transient hydraulic model. It comprises spatially distributed pilot points and overparameterisation constrained by regularisation and calibration to head differences. Significance of the results is demonstrated by cross validation. With this approach the spatial distribution of an aquitard have been identified with high precision. The highly transient and heterogeneous flow conditions are specified and a new viewpoint on the geologic formation of Lake Tegel is obtained. The good fit of modelled and observed breakthrough curves of 18O, chloride and temperature by just using transferred parameters obtained with the previous hydraulic methods, show the very good model performance and predictive capabilities. The intrusion of atmospheric oxygen into the unsaturated zone is identified to be the principal redox determining factor during infiltration. Previously inconsistent and also local geochemical conditions are identified to be determined by interaction of infiltration processes with the spatial extent of the aquitard. A theory for chemical clogging of abstraction wells is developed, identifying the strong vertical redox zoning as principal factor of influence.

La modélisation inverse de la dispersion atmosphérique consiste à reconstruire les caractéristiques d’une source (quantité de polluants rejetée, position) à partir de mesures de concentration dans l’air, en utilisant un modèle direct de dispersion et un algorithme d’inversion. Nous avons utilisé dans cette étude deux modèles directs de dispersion atmosphérique SIRANE (Soulhac, 2000; Soulhac et al., 2011) et SIRANERISK (Cierco et Soulhac, 2009a; Lamaison et al., 2011a, 2011b). Il s’agit de deux modèles opérationnels de « réseau des rues », basés sur le calcul du bilan de masse à différents niveaux du réseau. Leur concept permet de décrire correctement les différents phénomènes physiques de dispersion et de transport de la pollution atmosphérique dans des réseaux urbains complexes. L’étude de validation de ces deux modèles directs de dispersion a été effectuée après avoir évalué la fiabilité des paramétrages adoptés pour simuler les échanges verticaux entre la canopée et l'atmosphère, les transferts aux intersections de rues et la canalisation de l’écoulement à l’intérieur du réseau de rues. Pour cela, nous avons utilisé des mesures en soufflerie effectuées dans plusieurs configurations académiques. Nous avons développé au cours de cette thèse un système de modélisation inverse de dispersion atmosphérique (nommé ReWind) qui consiste à déterminer les caractéristiques d’une source de polluant (débit, position) à partir des concentrations mesurées, en résolvant numériquement le système matriciel linéaire qui relie le vecteur des débits au vecteur des concentrations. La fiabilité des résultats et l’optimisation des temps de calcul d’inversion sont assurées par le couplage de plusieurs méthodes mathématiques de résolution et d’optimisation, bien adaptées pour traiter le cas des problèmes mal posés. L’étude de sensibilité de cet algorithme d’inversion à certains paramètres d’entrée (comme les conditions météorologiques, les positions des récepteurs,…) a été effectuée en utilisant des observations synthétiques (fictives) fournies par le modèle direct de dispersion atmosphérique. La spécificité des travaux entrepris dans le cadre de ce travail a consisté à appliquer ReWind dans des configurations complexes de quartier urbain, et à utiliser toute la variabilité turbulente des mesures expérimentales obtenues en soufflerie pour qualifier ses performances à reconstruire les paramètres sources dans des conditions représentatives de situations de crise en milieu urbain ou industriel. L’application de l’approche inverse en utilisant des signaux instantanés de concentration mesurés en soufflerie plutôt que des valeurs moyennes, a montré que le modèle ReWind fournit des résultats d’inversion qui sont globalement satisfaisants et particulièrement encourageants en termes de reproduction de la quantité de masse totale de polluant rejetée dans l’atmosphère. Cependant, l’algorithme présente quelques difficultés pour estimer à la fois le débit et la position de la source dans certains cas. En effet, les résultats de l’inversion sont assez influencés par le critère de recherche (d’optimisation), le nombre de récepteurs impactés par le panache, la qualité des observations et la fiabilité du modèle direct de dispersion atmosphérique. / The aim of this study is to develop an inverse atmospheric dispersion model for crisis management in urban areas and industrial sites. The inverse modes allows for the reconstruction of the characteristics of a pollutant source (emission rate, position) from concentration measurements, by combining a direct dispersion model and an inversion algorithm, and assuming as known both site topography and meteorological conditions. The direct models used in these study, named SIRANE and SIRANERISK, are both operational "street network" models. These are based on the decomposition of the urban atmosphere into two sub-domains: the urban boundary layer and the urban canopy, represented as a series of interconnected boxes. Parametric laws govern the mass exchanges between the boxes under the assumption that the pollutant dispersion within the canopy can be fully simulated by modelling three main bulk transfer phenomena: channelling along street axes, transfers at street intersections and vertical exchange between a street canyon and the overlying atmosphere. The first part of this study is devoted to a detailed validation of these direct models in order to test the parameterisations implemented in them. This is achieved by comparing their outputs with wind tunnel experiments of the dispersion of steady and unsteady pollutant releases in idealised urban geometries. In the second part we use these models and experiments to test the performances of an inversion algorithm, named REWind. The specificity of this work is twofold. The first concerns the application of the inversion algorithm - using as input data instantaneous concentration signals registered at fixed receptors and not only time-averaged or ensemble averaged concentrations. - in urban like geometries, using an operational urban dispersion model as direct model. The application of the inverse approach by using instantaneous concentration signals rather than the averaged concentrations showed that the ReWind model generally provides reliable estimates of the total pollutant mass discharged at the source. However, the algorithm has some difficulties in estimating both emission rate and position of the source. We also show that the performances of the inversion algorithm are significantly influenced by the cost function used to the optimization, the number of receptors and the parameterizations adopted in the direct atmospheric dispersion model.

Les eaux souterraines constituent un des principaux facteurs de contrôle influençant l'activité des glissements de terrain. L'étude hydrogéologique constitue donc une étape incontournable dans la compréhension de leur fonctionnement, en vue de prédire leur déplacement ou d'adapter les solutions d'assainissement. Cependant, les aquifères qui se développent au sein des glissements de terrain sont généralement complexes, et la compréhension de l'hydrogéologie des grands glissements reste encore aujourd'hui une problématique majeure. Dans un contexte climatique extrême de type tropical humide, les glissements du cirque de Salazie (île de La Réunion) constituent des cas d'étude particulièrement porteurs pour l'apport de nouvelles connaissances sur l'hydrogéologie et le fonctionnement des mouvements de terrain de grande ampleur. Pour caractériser ces glissements, une étude pluridisciplinaire associant géologie, hydrochimie, modélisation numérique (globale et déterministe), est mise en œuvre. Sur le secteur nord du cirque de Salazie, l'étude géologique réalisée permet de reconstruire la géométrie de la base d'une instabilité de très grande ampleur, dont le volume est estimé à 370 Mm3. Le glissement de Grand Ilet (175 Mm3) en est la partie la plus active. Les formations volcano-détritiques (brèches) constituant le glissement sont des dépôts d'une ancienne avalanche de débris (> 350 Ma), réactivés localement. Cette étude met également en évidence un régime de déformation extensif associée à l'activité du glissement, favorisant le développement de zones décomprimées subverticales. Au sein des brèches constituant la zone instable, un aquifère continu est identifié. Les écoulements souterrains y sont contrôlés par la géométrie de la base du glissement. Par ailleurs, la nature du couvert végétal en surface, le climat, les propriétés des sols et l'épaisseur de la zone non saturée font que seuls les épisodes pluvieux intenses dont le cumul dépasse 80 mm/ évènement génèrent une recharge de l'aquifère des brèches. Lors de cette recharge, les zones décomprimées favorisent des transferts rapides (< 1 jour) à travers la zone non-saturée, malgré des épaisseurs pouvant atteindre localement 160 mètres. Une infiltration plus lente percole à travers la matrice poreuse des brèches. Pour l'étude de la dynamique des grands glissements dans le cirque de Salazie, l'application de l'outil de modélisation inverse met en évidence que les variations de vitesses de déplacement sont directement reliée aux processus de recharge et de tarissement de la nappe des brèches. Ces modèles s'avèrent particulièrement performants pour prédire les déplacements des grands glissements et comprendre leur fonctionnement. / The landslide activity is commonly controlled by the variation of hydraulic head inside the instable mass. Thus, the hydrogeological study of landslides is an essential step to predict landslide dynamic, and for the remediation choices. However, the aquifers developed in landslide are generally complexes. Actually, the comprehension of landslide hydrogeological functioning is a major problematic, especially for the deep-seated landslides.Under humid tropical climate, the Salazie landslides (Reunion Island) allow to implement an interesting study to characterized deep-seated landslide hydrogeology and functioning. In this study, we performed a multidisciplinary approach, combining geology, hydrochemistry and numerical modeling (global and deterministic);The geological study allows the construction of the deep seated landslide geological model (Grand Ilet sector). The Grand Ilet landslide, corresponding to 175 Mm3, is the most active part of deep seated instability whose the total volume is estimated at 370 Mm3. The volcano-detritic lithologies (i.e basic breccia) constitute the main geological formation of the unstable mass. This breccia is a present-day reactivation of an old destabilization (> 350 Ma). Furthermore, the extensive deformation generated by the landslide activity allows the formation of decompressed zones.Inside breccias, a continuous aquifer is identified. Here, the groundwater flows are controlled by the geometry of the landslide base. Natural land cover, soils properties, unsaturated zone thickness and warm temperature limit the groundwater recharge. Only the intense rainfall episode (80 mm/event) can generated the recharge of landslide aquifer. During these events, rapid transfers circulating inside the decompressed zones have a significant effect on hydraulic charge variations. A slow component is infiltrated inside the porous medium of breccias.The inverse modeling methods with bimodal Gaussian-Exponential transfer function is applied to study the Salazie landslides dynamic. The results show that landslide speed variations are directly controlled by groundwater level variations during the hydrologic cycle (recharge and recession stage). Consequently, the inverse model is a powerful tool for predicting deep-seated landslide movements and for studying how they function.

In the medtech industry models of the cardiiovascular systems and simulations are valuable tools for the development of new products ad therapies. The simulator Aplysia has been developed over several decade and is able to replicate a wide range of phenomena involved in the physiology and pathophysiology of breathing and circulation. Aplysia is also able to simulate the hemodynamics phenomena starting from a set of patient model parameters enhancing the idea of a "digital twin", i.e. a patient-specific representative simulation. Having a good starting estimate of the patient model parameters is a crucial aspect to start the simulation. A first estimate can be given by looking at patient monitoring data but medical expertise is required. The goal of this thesis is to address the parameter estimation task by developing machine learning and deep learning model to give an estimate of the patient model parameter starting from a set of time-varying data that we will refers as state variables. Those state variables are descriptive of a specific patient and for our project we will generate them through Aplysia starting from the simulation presets already available in the framework. Those presets simulates different physiologies, from healthy cases to different cardiovascular diseases. The thesis propose a comparison between a machine learning pipeline and more complex deep learning architecture to simultaneously predicting all the model parameters. This task is referred as Multi Target Regression (MTR) so the performances will be assessed in terms of MTR performance metrics. The results shows that a gradient boosting regressor with a regressor-stacking approach achieve overall good performances, still it shows some lack of performances on some target model parameters. The deep learning architectures did not produced any valuable results because of the amount of our data: to deploy deep architectures such as ResNet or more complex Convolutional Neural Network (CNN) we need more simulations then the one that were done for this thesis work. / Simulatorn Aplysia har under flera decennier utvecklats för forskning och FoU inom området kardiovaskulära systemmodeller och simuleringar och kan idag replikera ett brett spektrum av fenomen involverade i andningens och cirkulationens fysiologi och patofysiologi. Aplysia kan också simulera hemodynamiska fenomen med utgångspunkt från en uppsättning patientmodellparametrar och detta förstärker idén om en digital tvilling", det vill säga en patientspecifik representativ simulering. Att ha en bra startuppskattning av patientmodellens parametrar är en avgörande aspekt för att starta simuleringen. En första uppskattning kan ges genom att titta på patientövervakningsdata men medicinsk expertis krävs för tolkningen av sådana data. Målet med denna mastersuppsats är att addressera parameteruppskattningsuppgiften genom att utveckla maskininlärnings-och djupinlärningsmodeller för att erhålla en uppskattning av patientmodellparametrar utgående från en uppsättning tidsvarierande data som vi kommer att referera till som tillståndsvariabler. Dessa tillståndsvariabler är beskrivande för en specifik patient och för vårt projekt kommer vi att generera dem med hjälp av Aplysia med utgångspunkt från de modellförinställningar som redan finns tillgängliga i ramverket. Dessa förinställningar simulerar olika fysiologier, från friska fall till olika hjärt-kärlsjukdomar. Uppsatsen presenterar en jämförelse mellan en maskininlärningspipeline och en mer komplex djupinlärningsarkitektur för att samtidigt förutsäga alla modellparametrar. Denna uppgift bygger på MTR så resulterande prestanda kommer att bedömas i termer av MTR prestationsmått. Resultaten visar att en gradientförstärkande regressor med en regressor-stacking-metod uppnår överlag goda resultat, ändå visar den en viss brist på prestanda på vissa målmodellparametrar. Deep learning-arkitekturerna gav inga värdefulla resultat på grund av den begränsade mängden av data vi kunde generera. För att träna djupa arkitekturer som ResNet eller mer komplexa CNN behöver vi fler simuleringar än den som gjordes för detta examensarbete.

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 confronte aux problèmes de non unicité des solutions du problème inverse. Ces contraintes limitent l’interprétation des modèles a un aspect qualitatif de la distribution des contrastes de résistivité modélises en 2D ou 3D, résultant du choix des paramètres d’inversion et de l’association de paramètres du milieu non identifiables a 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 empilements 2D de billes de verre sphériques (de l’ordre de la 100èn de μm) dans un réservoir en plexiglas transparent. Deux réseaux verticaux de 21 électrodes sont disposes sur les bords latéraux du réservoir pour effectuer le suivi ERT du traceur sale (NaCl dissout) a 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 dispose face a un panneau lumineux permettant de réalise 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 a 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 appliques (L1 ou L2) sur les données et les paramètres, et au facteur d’amortissement _. Des valeurs trop élevées de _ et 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 de α et du facteur d’acceptante donne des résultats plus proches l’analyse vidéo, mais produit des effets de bosses à l’avant et à l’arrière du front.La vitesse interstitiel 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 impose 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 ou 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 testes ici n’ont pas d’influences significatives sur la dispersion qui est dominée par le débit impose a 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éalises 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ées 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 paramètres de transport u, α (ou D), et Rf . / Electrical resistivity tomography is a common geophysical method often used to detect and follow plumepollutants in aquifers. However ERT is an integrative method whose reliability of the models is faced tothe non-unicity of the inverse problem solutions. These constraints limit the interpretation to a qualitativeview of the resistivity contrasts modelled in 2D or 3D, resulting of the chosen inverse parameters and thecombination of several hydrodynamic paramaters related to the poral network.The purpose of this thesis was to test the abilities of the ERT imaging to quantify solute transport parametersin miscible displacement occurred in groundwater and the sensitivity of inverse parameters most affectingthe modelled dispersion front.Laboratory experiments are conducted on glass beads poured into a transparent plexiglas container. Twovertical lines of 21 stainless steel electrodes are fixed on the lateral sides of the container to perform the ERTmonitoring, of the NaCl dissolved in the tracer, from a sequence of 210 quadripole measurements acquiredin transverse dipole-dipole each 5 minutes. A light panel is placed behind the experimental device and avideo follow up of the dyed part of the tracer is acquired from the other side.Video analysis reveal a faster propagation of the dye in contact of the vertical edge, which is negligible forthe lowest flow rates imposed by the pump. In contrast, ERT mesurements are strongly disturbed by theresistant edges of the plexiglas container which affect the resulting models. Normalisation of the apparentresisitivity measurements acquired at the experimental stage and by the Archie’s law strongly tones downthese resistive artefacts.ERT modelling is here particularly sensitive to the grid mesh, the norm (L1 or L2 ) applied on data andparameters, and the damping factor λ. High values of λ and the cutoff factor tend to smooth the resistivityconstrasts in the area of the mixing front and increase the weight of the horizontal edge effects on D andα. While results from inverse modelling constraint by low λ and cutoff factors are much closer to the videoanalysis but with enhanced side slope effects at the rear and the front of the mixing area.The interstitial velocity u is independant of the chosen inverse parameters. For both methods u is alwaysinferior to the flow rate provided by the pump, whose the gap is expressed as a retardation factor Rf . Thisretardation is due to adsorption of Na+ on the beads surfaces, which contributes to delay the dispersionfront followed by ERT. The retardation expressed by the video analysis can be due to the important sizeof the molecule of the dye which is easily slowed down in lower permeability areas, added to an eventualaffinity to the solid surface.The ranges of fluid conductivity contrasts and stability of the interfaces tested here have no influences onthe dispersion which is dominated by the flow velocity u. Estimations of the dispersion coefficient D asfunction of the Péclet is consistent with the theoretical curve of Bachmat (1968) and Bijeljic & al (2004).Field experiments are first conducted in 2D on homogeneous unsaturated sand which is considered as aninfinite half-space. However, data normalisation is much more complicated since the saturation state andthe initial fluid conductivities need to be estimated to calibrate the Archie’s law. Because of the 3D tracerinflitration, the RMS error of 2D-ERT models highlights that the inversion process is not enough constraintby data which does not allow to quantify the transport parameters. 3D experiments were then adaptedto detect and follow plumes of saline tracers injected in the centre of the electrode device. From 3D ERTmeasurements we are able to produce reliable models in order to estimate such transport parameters as themean flow velocity, and transverse and longitudinal dispersivities.