Point defects and texture of Pd1-xInx sputtered intermetallic thin films

Pd1-xInx thin films (0.4 < x < 0.56) were prepared by radio frequency sputtering from a multi-zone target. The properties of these Hume-Rothery alloys were studied by X-ray diffractometry, electron probe microanalysis and scanning tunneling microscopy. The diffraction spectra were analyzed to obtain the intensity ratio of the (100) superlattice line to the (200) normal line, together with the variations of the lattice constant. The results ape explained quantitatively by a model based on point defects, i.e. Pd vacancies in In-rich films and Pd antisite atoms in Pd-rich films. In-rich films grow preferentially in the [100] direction while Pd-rich films grow preferentially in the [110] direction. The grains in indium-rich sputtered films appear to be enclosed in an atomically thick, indium-rich layer. The role of texture and the influence of point defects on electrical resistivity is also reported. (C) 1996 Elsevier Science Limited.

Rate and processes of river network rearrangement during incipient faulting: The case of the Cahabon river, Guatemala

Deeply incised river networks are generally regarded as robust features that are not easily modified by erosion or tectonics. Although the reorganization of deeply incised drainage systems has been documented, the corresponding importance with regard to the overall landscape evolution of mountain ranges and the factors that permit such reorganizations are poorly understood. To address this problem, we have explored the rapid drainage reorganization that affected the Cahabon River in Guatemala during the Quaternary. Sediment-provenance analysis, field mapping, and electrical resistivity tomography (ERT) imaging are used to reconstruct the geometry of the valley before the river was captured. Dating of the abandoned valley sediments by the Be-10-Al-26 burial method and geomagnetic polarity analysis allow us to determine the age of the capture events and then to quantify several processes, such as the rate of tectonic deformation of the paleovalley, the rate of propagation of post-capture drainage reversal, and the rate at which canyons that formed at the capture sites have propagated along the paleovalley. Transtensional faulting started 1 to 3 million years ago, produced ground tilting and ground faulting along the Cahabon River, and thus generated differential uplift rate of 0.3 +/- 0.1 up to 0.7 +/- 0.4 mm . y(-1) along the river's course. The river responded to faulting by incising the areas of relative uplift and depositing a few tens of meters of sediment above the areas of relative subsidence. Then, the river experienced two captures and one avulsion between 700 ky and 100 ky. The captures breached high-standing ridges that separate the Cahabon River from its captors. Captures occurred at specific points where ridges are made permeable by fault damage zones and/or soluble rocks. Groundwater flow from the Cahabon River down to its captors likely increased the erosive power of the captors thus promoting focused erosion of the ridges. Valley-fill formation and capture occurred in close temporal succession, suggesting a genetic link between the two. We suggest that the aquifers accumulated within the valley-fills, increased the head along the subterraneous system connecting the Cahabon River to its captors, and promoted their development. Upon capture, the breached valley experienced widespread drainage reversal toward the capture sites. We attribute the generalized reversal to combined effects of groundwater sapping in the valley-fill, axial drainage obstruction by lateral fans, and tectonic tilting. Drainage reversal increased the size of the captured areas by a factor of 4 to 6. At the capture sites, 500 m deep canyons have been incised into the bedrock and are propagating upstream at a rate of 3 to 11 mm . y(-1) deepening at a rate of 0.7 to 1 5 mm . y(-1). At this rate, 1 to 2 million years will be necessary for headward erosion to completely erase the topographic expression of the paleovalley. It is concluded that the rapid reorganization of this drainage system was made possible by the way the river adjusted to the new tectonic strain field, which involved transient sedimentation along the river's course. If the river had escaped its early reorganization and had been given the time necessary to reach a new dynamic equilibrium, then the transient conditions that promoted capture would have vanished and its vulnerability to capture would have been strongly reduced.

A filtering method to correct time-lapse 3D ERT data and improve imaging of natural aquifer dynamics

We have developed a processing methodology that allows crosshole ERT (electrical resistivity tomography) monitoring data to be used to derive temporal fluctuations of groundwater electrical resistivity and thereby characterize the dynamics of groundwater in a gravel aquifer as it is infiltrated by river water. Temporal variations of the raw ERT apparent-resistivity data were mainly sensitive to the resistivity (salinity), temperature and height of the groundwater, with the relative contributions of these effects depending on the time and the electrode configuration. To resolve the changes in groundwater resistivity, we first expressed fluctuations of temperature-detrended apparent-resistivity data as linear superpositions of (i) time series of riverwater-resistivity variations convolved with suitable filter functions and (ii) linear and quadratic representations of river-water-height variations multiplied by appropriate sensitivity factors; river-water height was determined to be a reliable proxy for groundwater height. Individual filter functions and sensitivity factors were obtained for each electrode configuration via deconvolution using a one month calibration period and then the predicted contributions related to changes in water height were removed prior to inversion of the temperature-detrended apparent-resistivity data. Applications of the filter functions and sensitivity factors accurately predicted the apparent-resistivity variations (the correlation coefficient was 0.98). Furthermore, the filtered ERT monitoring data and resultant time-lapse resistivity models correlated closely with independently measured groundwater electrical resistivity monitoring data and only weakly with the groundwater-height fluctuations. The inversion results based on the filtered ERT data also showed significantly less inversion artefacts than the raw data inversions. We observed resistivity increases of up to 10% and the arrival time peaks in the time-lapse resistivity models matched those in the groundwater resistivity monitoring data.

Time-lapse and probabilistic inversion strategies for plane-wave electromagnetic data.

L'utilisation efficace des systèmes géothermaux, la séquestration du CO2 pour limiter le changement climatique et la prévention de l'intrusion d'eau salée dans les aquifères costaux ne sont que quelques exemples qui démontrent notre besoin en technologies nouvelles pour suivre l'évolution des processus souterrains à partir de la surface. Un défi majeur est d'assurer la caractérisation et l'optimisation des performances de ces technologies à différentes échelles spatiales et temporelles. Les méthodes électromagnétiques (EM) d'ondes planes sont sensibles à la conductivité électrique du sous-sol et, par conséquent, à la conductivité électrique des fluides saturant la roche, à la présence de fractures connectées, à la température et aux matériaux géologiques. Ces méthodes sont régies par des équations valides sur de larges gammes de fréquences, permettant détudier de manières analogues des processus allant de quelques mètres sous la surface jusqu'à plusieurs kilomètres de profondeur. Néanmoins, ces méthodes sont soumises à une perte de résolution avec la profondeur à cause des propriétés diffusives du champ électromagnétique. Pour cette raison, l'estimation des modèles du sous-sol par ces méthodes doit prendre en compte des informations a priori afin de contraindre les modèles autant que possible et de permettre la quantification des incertitudes de ces modèles de façon appropriée. Dans la présente thèse, je développe des approches permettant la caractérisation statique et dynamique du sous-sol à l'aide d'ondes EM planes. Dans une première partie, je présente une approche déterministe permettant de réaliser des inversions répétées dans le temps (time-lapse) de données d'ondes EM planes en deux dimensions. Cette stratégie est basée sur l'incorporation dans l'algorithme d'informations a priori en fonction des changements du modèle de conductivité électrique attendus. Ceci est réalisé en intégrant une régularisation stochastique et des contraintes flexibles par rapport à la gamme des changements attendus en utilisant les multiplicateurs de Lagrange. J'utilise des normes différentes de la norme l2 pour contraindre la structure du modèle et obtenir des transitions abruptes entre les régions du model qui subissent des changements dans le temps et celles qui n'en subissent pas. Aussi, j'incorpore une stratégie afin d'éliminer les erreurs systématiques de données time-lapse. Ce travail a mis en évidence l'amélioration de la caractérisation des changements temporels par rapport aux approches classiques qui réalisent des inversions indépendantes à chaque pas de temps et comparent les modèles. Dans la seconde partie de cette thèse, j'adopte un formalisme bayésien et je teste la possibilité de quantifier les incertitudes sur les paramètres du modèle dans l'inversion d'ondes EM planes. Pour ce faire, je présente une stratégie d'inversion probabiliste basée sur des pixels à deux dimensions pour des inversions de données d'ondes EM planes et de tomographies de résistivité électrique (ERT) séparées et jointes. Je compare les incertitudes des paramètres du modèle en considérant différents types d'information a priori sur la structure du modèle et différentes fonctions de vraisemblance pour décrire les erreurs sur les données. Les résultats indiquent que la régularisation du modèle est nécessaire lorsqu'on a à faire à un large nombre de paramètres car cela permet d'accélérer la convergence des chaînes et d'obtenir des modèles plus réalistes. Cependent, ces contraintes mènent à des incertitudes d'estimations plus faibles, ce qui implique des distributions a posteriori qui ne contiennent pas le vrai modèledans les régions ou` la méthode présente une sensibilité limitée. Cette situation peut être améliorée en combinant des méthodes d'ondes EM planes avec d'autres méthodes complémentaires telles que l'ERT. De plus, je montre que le poids de régularisation des paramètres et l'écart-type des erreurs sur les données peuvent être retrouvés par une inversion probabiliste. Finalement, j'évalue la possibilité de caractériser une distribution tridimensionnelle d'un panache de traceur salin injecté dans le sous-sol en réalisant une inversion probabiliste time-lapse tridimensionnelle d'ondes EM planes. Etant donné que les inversions probabilistes sont très coûteuses en temps de calcul lorsque l'espace des paramètres présente une grande dimension, je propose une stratégie de réduction du modèle ou` les coefficients de décomposition des moments de Legendre du panache de traceur injecté ainsi que sa position sont estimés. Pour ce faire, un modèle de résistivité de base est nécessaire. Il peut être obtenu avant l'expérience time-lapse. Un test synthétique montre que la méthodologie marche bien quand le modèle de résistivité de base est caractérisé correctement. Cette méthodologie est aussi appliquée à un test de trac¸age par injection d'une solution saline et d'acides réalisé dans un système géothermal en Australie, puis comparée à une inversion time-lapse tridimensionnelle réalisée selon une approche déterministe. L'inversion probabiliste permet de mieux contraindre le panache du traceur salin gr^ace à la grande quantité d'informations a priori incluse dans l'algorithme. Néanmoins, les changements de conductivités nécessaires pour expliquer les changements observés dans les données sont plus grands que ce qu'expliquent notre connaissance actuelle des phénomenès physiques. Ce problème peut être lié à la qualité limitée du modèle de résistivité de base utilisé, indiquant ainsi que des efforts plus grands devront être fournis dans le futur pour obtenir des modèles de base de bonne qualité avant de réaliser des expériences dynamiques. Les études décrites dans cette thèse montrent que les méthodes d'ondes EM planes sont très utiles pour caractériser et suivre les variations temporelles du sous-sol sur de larges échelles. Les présentes approches améliorent l'évaluation des modèles obtenus, autant en termes d'incorporation d'informations a priori, qu'en termes de quantification d'incertitudes a posteriori. De plus, les stratégies développées peuvent être appliquées à d'autres méthodes géophysiques, et offrent une grande flexibilité pour l'incorporation d'informations additionnelles lorsqu'elles sont disponibles. -- The efficient use of geothermal systems, the sequestration of CO2 to mitigate climate change, and the prevention of seawater intrusion in coastal aquifers are only some examples that demonstrate the need for novel technologies to monitor subsurface processes from the surface. A main challenge is to assure optimal performance of such technologies at different temporal and spatial scales. Plane-wave electromagnetic (EM) methods are sensitive to subsurface electrical conductivity and consequently to fluid conductivity, fracture connectivity, temperature, and rock mineralogy. These methods have governing equations that are the same over a large range of frequencies, thus allowing to study in an analogous manner processes on scales ranging from few meters close to the surface down to several hundreds of kilometers depth. Unfortunately, they suffer from a significant resolution loss with depth due to the diffusive nature of the electromagnetic fields. Therefore, estimations of subsurface models that use these methods should incorporate a priori information to better constrain the models, and provide appropriate measures of model uncertainty. During my thesis, I have developed approaches to improve the static and dynamic characterization of the subsurface with plane-wave EM methods. In the first part of this thesis, I present a two-dimensional deterministic approach to perform time-lapse inversion of plane-wave EM data. The strategy is based on the incorporation of prior information into the inversion algorithm regarding the expected temporal changes in electrical conductivity. This is done by incorporating a flexible stochastic regularization and constraints regarding the expected ranges of the changes by using Lagrange multipliers. I use non-l2 norms to penalize the model update in order to obtain sharp transitions between regions that experience temporal changes and regions that do not. I also incorporate a time-lapse differencing strategy to remove systematic errors in the time-lapse inversion. This work presents improvements in the characterization of temporal changes with respect to the classical approach of performing separate inversions and computing differences between the models. In the second part of this thesis, I adopt a Bayesian framework and use Markov chain Monte Carlo (MCMC) simulations to quantify model parameter uncertainty in plane-wave EM inversion. For this purpose, I present a two-dimensional pixel-based probabilistic inversion strategy for separate and joint inversions of plane-wave EM and electrical resistivity tomography (ERT) data. I compare the uncertainties of the model parameters when considering different types of prior information on the model structure and different likelihood functions to describe the data errors. The results indicate that model regularization is necessary when dealing with a large number of model parameters because it helps to accelerate the convergence of the chains and leads to more realistic models. These constraints also lead to smaller uncertainty estimates, which imply posterior distributions that do not include the true underlying model in regions where the method has limited sensitivity. This situation can be improved by combining planewave EM methods with complimentary geophysical methods such as ERT. In addition, I show that an appropriate regularization weight and the standard deviation of the data errors can be retrieved by the MCMC inversion. Finally, I evaluate the possibility of characterizing the three-dimensional distribution of an injected water plume by performing three-dimensional time-lapse MCMC inversion of planewave EM data. Since MCMC inversion involves a significant computational burden in high parameter dimensions, I propose a model reduction strategy where the coefficients of a Legendre moment decomposition of the injected water plume and its location are estimated. For this purpose, a base resistivity model is needed which is obtained prior to the time-lapse experiment. A synthetic test shows that the methodology works well when the base resistivity model is correctly characterized. The methodology is also applied to an injection experiment performed in a geothermal system in Australia, and compared to a three-dimensional time-lapse inversion performed within a deterministic framework. The MCMC inversion better constrains the water plumes due to the larger amount of prior information that is included in the algorithm. The conductivity changes needed to explain the time-lapse data are much larger than what is physically possible based on present day understandings. This issue may be related to the base resistivity model used, therefore indicating that more efforts should be given to obtain high-quality base models prior to dynamic experiments. The studies described herein give clear evidence that plane-wave EM methods are useful to characterize and monitor the subsurface at a wide range of scales. The presented approaches contribute to an improved appraisal of the obtained models, both in terms of the incorporation of prior information in the algorithms and the posterior uncertainty quantification. In addition, the developed strategies can be applied to other geophysical methods, and offer great flexibility to incorporate additional information when available.

Localization and characterization of an active fault in an urbanized area in central Guatemala by means of geoelectrical imaging

The Polochic and Motagua faults define the active plate boundary between the North American and Caribbean plates in central Guatemala. A splay of the Polochic Fault traverses the rapidly growing city of San Miguel Uspantan that is periodically affected by destructive earthquakes. This fault splay was located using a 2D electrical resistivity tomography (ERT) survey that also characterized the fault damage zone and evaluated the thickness and nature of recent deposits upon which most of the city is built. ERT images show the fault as a similar to 50 m wide, near-vertical low-resistivity anomaly, bounded within a few meters by high resistivity anomalies. Forward modeling reproduces the key aspects of the observed electrical resistivity data with remarkable fidelity thus defining the overall location, geometry, and internal structure of the fault zone as well as the affected lithologies. Our results indicate that the city is constructed on a similar to 20 m thick surficial layer consisting of poorly consolidated, highly porous, water-logged pumice. This soft layer is likely to amplify seismic waves and to liquefy upon moderate to strong ground shaking. The electrical conductivity as well as the major element chemistry of the groundwater provides evidence to suggest that the local aquifer might, at least in part, be fed by water rising along the fault. Therefore, the potential threat posed by this fault splay may not be limited to its seismic activity per se, but could be compounded its potential propensity to enhance seismic site effects by injecting water into the soft surficial sediments. The results of this study provide the basis for a rigorous analysis of seismic hazard and sustainable development of San Miguel Uspantan and illustrate the potential of ERT surveying for paleoseismic studies.

Screen printed La2/3Sr1/3MnO3 thick films on alumina substrates

We report here on the preparation of La2/3Sr1/3MnO3 magnetoresistive thick films on polycrystalline Al2O3 substrates by using the screen printing technique. It is shown that films can be obtained using high temperature sintering. While there is a reacted layer, this improves adhesion and is not too troublesome if the films are made thick enough. It is shown that PbO-B2O3-SiO2 glass additives allow sintering at lower temperatures and can be used to improve the mechanical stress of the films. However, it is found that glass concentrations large enough to significantly improve the film adherence result in a weak low field magnetoresistance probably because grains are coated with high resistivity material. Strategies to overcome these difficulties are discussed.

Gas-sensing properties of sprayed films of (CdO)x(ZnO)1-x mixed oxide

A novel NO2 sensor based on (CdO)x(ZnO)1-x mixed-oxide thin films deposited by the spray pyrolysis technique is developed. The sensor response to 3-ppm NO2 is studied in the range 50°C-350°C for three different film compositions. The device is also tested for other harmful gases, such as CO (300 ppm) and CH4 (3000 ppm). The sensor response to these reducing gases is different at different temperatures varying from the response typical for the p-type semiconductor to that typical for the n-type semiconductor. Satisfactory response to NO2 and dynamic behavior at 230°C, as well as low resistivity, are observed for the mixed-oxide film with 30% Cd. The response to interfering gas is poor at working temperature (230°C). On the basis of this study, a possible sensing mechanism is proposed.

Enhanced electron-electron correlations in nanometric epitaxial SrRuO3 epitaxial films

Epitaxial and fully strained SrRuO3 thin films have been grown on SrTiO3(100). At initial stages the growth mode is three-dimensional- (3D-)like, leading to a finger-shaped structure aligned with the substrate steps and that eventually evolves into a 2D step-flow growth. We study the impact that the defect structure associated with this unique growth mode transition has on the electronic properties of the films. Detailed analysis of the transport properties of nanometric films reveals that microstructural disorder promotes a shortening of the carrier mean free path. Remarkably enough, at low temperatures, this results in a reinforcement of quantum corrections to the conductivity as predicted by recent models of disordered, strongly correlated electronic systems. This finding may provide a simple explanation for the commonly observed¿in conducting oxides-resistivity minima at low temperature. Simultaneously, the ferromagnetic transition occurring at about 140 K, becomes broader as film thickness decreases down to nanometric range. The relevance of these results for the understanding of the electronic properties of disordered electronic systems and for the technological applications of SrRuO3¿and other ferromagnetic and metallic oxides¿is stressed.

Influence of RF power on the properties of sputtered ZnO:Al thin films

Transparent conducting, aluminium doped zinc oxide thin films (ZnO:Al) were deposited by radio frequency (RF) magnetron sputtering. The RF power was varied from 60 to 350Wwhereas the substrate temperature was kept at 160 °C. The structural, electrical and optical properties of the as-deposited films were found to be influenced by the deposition power. The X-ray diffraction analysis showed that all the films have a strong preferred orientation along the [001] direction. The crystallite size was varied from 14 to 36 nm, however no significant change was observed in the case of lattice constant. The optical band gap varied in the range 3.44-3.58 eV. The lowest resistivity of 1.2×10 -3Vcm was shown by the films deposited at 250 W. The mobility of the films was found to increase with the deposition power.

Transparent conducting thin films by the co-sputtering of ZnO-ITO targets

Transparent and conductive Zn-In-Sn-O (ZITO) amorphous thin films have been deposited at room temperature by the rf magnetron co-sputtering of ITO and ZnO targets. Co-sputtering gives the possibility to deposit multicomponent oxide thin films with different compositions by varying the power to one of the targets. In order to make ZITO films with different Zn content, a constant rf power of 50 W was used for the ITO target, where as the rf power to ZnO target was varied from 25 W to 150 W. The as deposited films showed an increase in Zn content ratio from 17 to 67 % as the power to ZnO target was increased from 25 to 150 W. The structural, electrical and optical properties of the as deposited films are reported. The films showed an average transmittance over 80% in the visible wavelength range. The electrical resistivity and optical band gap of the ZITO films were found to depend on the Zn content in the film. The ZITO films deposited at room temperature with lower Zn content ratios showed better optical transmission and electrical properties compared to ITO film.

Electrical and optical properties of Zn-In-Sn-O transparent conducting thin films

Indium tin oxide (ITO) is one of the widely used transparent conductive oxides (TCO) for application as transparent electrode in thin film silicon solar cells or thin film transistors owing to its low resistivity and high transparency. Nevertheless, indium is a scarce and expensive element and ITO films require high deposition temperature to achieve good electrical and optical properties. On the other hand, although not competing as ITO, doped Zinc Oxide (ZnO) is a promising and cheaper alternative. Therefore, our strategy has been to deposit ITO and ZnO multicomponent thin films at room temperature by radiofrequency (RF) magnetron co-sputtering in order to achieve TCOs with reduced indium content. Thin films of the quaternary system Zn-In-Sn-O (ZITO) with improved electrical and optical properties have been achieved. The samples were deposited by applying different RF powers to ZnO target while keeping a constant RF power to ITO target. This led to ZITO films with zinc content ratio varying between 0 and 67%. The optical, electrical and morphological properties have been thoroughly studied. The film composition was analysed by X-ray Photoelectron Spectroscopy. The films with 17% zinc content ratio showed the lowest resistivity (6.6 × 10 - 4 Ω cm) and the highest transmittance (above 80% in the visible range). Though X-ray Diffraction studies showed amorphous nature for the films, using High Resolution Transmission Electron Microscopy we found that the microstructure of the films consisted of nanometric crystals embedded in a compact amorphous matrix. The effect of post deposition annealing on the films in both reducing and oxidizing atmospheres were studied. The changes were found to strongly depend on the zinc content ratio in the films.

Magnetic field effect on quantum corrections to the low-temperature conductivity in mettalic perovskite oxides

The transport and magnetotransport properties of the metallic and ferromagnetic SrRuO3 (SRO) and the metallic and paramagnetic LaNiO3 (LNO) epitaxial thin films have been investigated in fields up to 55 T at temperatures down to 1.8 K . At low temperatures both samples display a well-defined resistivity minimum. We argue that this behavior is due to the increasing relevance of quantum corrections to the conductivity (QCC) as temperature is lowered; this effect being particularly relevant in these oxides due to their short mean free path. However, it is not straightforward to discriminate between contributions of weak localization and renormalization of electron-electron interactions to the QCC through temperature dependence alone. We have taken advantage of the distinct effect of a magnetic field on both mechanisms to demonstrate that in ferromagnetic SRO the weak-localization contribution is suppressed by the large internal field leaving only renormalized electron-electron interactions, whereas in the nonmagnetic LNO thin films the weak-localization term is relevant.

Magneto-optical Kerr effect in laser-patterned La2/3Sr1/3MnO3 epitaxial thin films

In this study, we have performed magneto-optical Kerr effect (MOKE) measurement on epitaxial La2/3Sr1/3MnO3 thin films containing artificial interfaces created by laser-patterning the SrTiO3 substrate. The observed increase of the resistivity and of the high-field magnetoresistance when measuring the films across the interface arrays are related to the reduction of the magnetization of the interfaces with respect to the rest of the film. As observed by the local MOKE probe, the structural disorder in the manganite film induced by the underlying patterned substrate leads to a large spin disorder responsible for a strong high-field susceptibility of the resistance.

Self-interference of charge carriers in ferromagnetic SrRuO3

We report a systematic study of the low-temperature electrical conductivity in a series of SrRuO3 epitaxial thin films. At relatively high temperature the films display the conventional metallic behavior. However, a well-defined resistivity minimum appears at low temperature. This temperature dependence can be well described in a weak localization scenario: the resistivity minimum arising from the competition of electronic self-interference effects and the normal metallic character. By appropriate selection of the film growth conditions, we have been able to modify the mean-free path of itinerant carriers and thus to tune the relative strength of the quantum effects. We show that data can be quantitatively described by available theoretical models.

Detection and Characterization of Preferential Flow Paths in the Downstream Area of a Hazardous Landfill

We have used surface-based electrical resistivity tomography to detect and characterize preferential hydraulic pathways in the immediate downstream area of an abandoned, hazardous landfill. The landfill occupies the void left by a former gravel pit and its base is close to the groundwater table and lacking an engineered barrier. As such, this site is remarkably typical of many small- to medium-sized waste deposits throughout the densely populated and heavily industrialized foreland on both sides of the Alpine arc. Outflows of pollutants lastingly contaminated local drinking water supplies and necessitated a partial remediation in the form of a synthetic cover barrier, which is meant to prevent meteoric water from percolating through the waste before reaching the groundwater table. Any future additional isolation of the landfill in the form of lateral barriers thus requires adequate knowledge of potential preferential hydraulic pathways for outflowing contaminants. Our results, inferred from a suite of tomographically inverted surfaced-based electrical resistivity profiles oriented roughly perpendicular to the local hydraulic gradient, indicate that potential contaminant outflows would predominantly occur along an unexploited lateral extension of the original gravel deposit. This finds its expression as a distinct and laterally continuous high-resistivity anomaly in the resistivity tomograms. This interpretation is ground-truthed through a litholog from a nearby well. Since the probed glacio-fluvial deposits are largely devoid of mineralogical clay, the geometry of hydraulic and electrical pathways across the pore space of a given lithological unit can be assumed to be identical, which allows for an order-of-magnitude estimation of the overall permeability structure. These estimates indicate that the permeability of the imaged extension of the gravel body is at least two to three orders-of-magnitude higher than that of its finer-grained embedding matrix. This corroborates the preeminent role of the high-resistivity anomaly as a potential preferential flow path.