Self-potential investigations of a gravel bar in a restored river corridor

Self-potentials (SP) are sensitive to water fluxes and concentration gradients in both saturated and unsaturated geological media, but quantitative interpretations of SP field data may often be hindered by the superposition of different source contributions and time-varying electrode potentials. Self-potential mapping and close to two months of SP monitoring on a gravel bar were performed to investigate the origins of SP signals at a restored river section of the Thur River in northeastern Switzerland. The SP mapping and subsequent inversion of the data indicate that the SP sources are mainly located in the upper few meters in regions of soil cover rather than bare gravel. Wavelet analyses of the time-series indicate a strong, but non-linear influence of water table and water content variations, as well as rainfall intensity on the recorded SP signals. Modeling of the SP response with respect to an increase in the water table elevation and precipitation indicate that the distribution of soil properties in the vadose zone has a very strong influence. We conclude that the observed SP responses on the gravel bar are more complicated than previously proposed semi-empiric relationships between SP signals and hydraulic head or the thickness of the vadose zone. We suggest that future SP monitoring in restored river corridors should either focus on quantifying vadose zone processes by installing vertical profiles of closely spaced SP electrodes or by installing the electrodes within the river to avoid signals arising from vadose zone processes and time-varying electrochemical conditions in the vicinity of the electrodes.

Full-waveform inversion of cross-hole ground-penetrating radar data to characterize a gravel aquifer close to the Thur River, Switzerland

Cross-hole radar tomography is a useful tool for mapping shallow subsurface electrical properties viz. dielectric permittivity and electrical conductivity. Common practice is to invert cross-hole radar data with ray-based tomographic algorithms using first arrival traveltimes and first cycle amplitudes. However, the resolution of conventional standard ray-based inversion schemes for cross-hole ground-penetrating radar (GPR) is limited because only a fraction of the information contained in the radar data is used. The resolution can be improved significantly by using a full-waveform inversion that considers the entire waveform, or significant parts thereof. A recently developed 2D time-domain vectorial full-waveform crosshole radar inversion code has been modified in the present study by allowing optimized acquisition setups that reduce the acquisition time and computational costs significantly. This is achieved by minimizing the number of transmitter points and maximizing the number of receiver positions. The improved algorithm was employed to invert cross-hole GPR data acquired within a gravel aquifer (4-10 m depth) in the Thur valley, Switzerland. The simulated traces of the final model obtained by the full-waveform inversion fit the observed traces very well in the lower part of the section and reasonably well in the upper part of the section. Compared to the ray-based inversion, the results from the full-waveform inversion show significantly higher resolution images. At either side, 2.5 m distance away from the cross-hole plane, borehole logs were acquired. There is a good correspondence between the conductivity tomograms and the natural gamma logs at the boundary of the gravel layer and the underlying lacustrine clay deposits. Using existing petrophysical models, the inversion results and neutron-neutron logs are converted to porosity. Without any additional calibration, the values obtained for the converted neutron-neutron logs and permittivity results are very close and similar vertical variations can be observed. The full-waveform inversion provides in both cases additional information about the subsurface. Due to the presence of the water table and associated refracted/reflected waves, the upper traces are not well fitted and the upper 2 m in the permittivity and conductivity tomograms are not reliably reconstructed because the unsaturated zone is not incorporated into the inversion domain.

Pollution and pollutant transport in the geosphere - An introduction to the symposium

Environmental research in earth sciences is focused on the geosphere, i.e. (1) waters and sediments of rivers, lakes and oceans, and (2) soils and underlying shallow rock formations,both water-unsaturated and -saturated. The subsurface is studied down to greater depths at sites where waste repositories or tunnels are planned and mining activities exist. In recent years, earth scientists have become more and more involved in pollution problems related to their classical field of interest, e.g. groundwater, ore deposits, or petroleum and non-metal natural deposits (gravel, clay, cement precursors). Major pollutants include chemical substances, radioactive isotopes and microorganisms. Mechanisms which govern the transport of pollutants are of physical, chemical (dissolution, precipitation, adsorption), or microbiological (transformation) nature. Land-use planning must reflect a sustainable development and sound scientific criteria. Today's environmental pollution requires working teams with an interdisciplinary background in earth sciences, hydrology, chemistry, biology, physics as well as engineering. This symposium brought together for the first time in Switzerland earth and soil scientists, physicists and chemists, to present and discuss environmental issues concerning the geosphere.

Conservation of a threatened European tree frog (Hyla arborea) metapopulation

Summary: Amphibians are among the most vulnerable animals of the world. One third of all species are currently threatened with extinction. Habitat loss is the major menace to pond- and stream-breeding species in the old world. In highly urbanized landscape like the Swiss Plateau, most species suffer from habitat reduction and fragmentation. Among all indigenous species, the European tree frog (Hyla arborea L., 1758) is one of the most endangered. It experienced an alarming decline during the last century and its regional long-term persistence is not guaranteed. We developed a monitoring framework based on calling male counts which included multiple visits to each wetland during the reproduction period in order to precisely determine its distribution on the Lemanic coast. Our results indicate that visiting populations 3 limes under suitable climatic conditions (temperature >20°C) provides reliable presence/absence data. Based on our monitoring data, we analyzed the species requirements regarding its breeding habitat. It appeared that anthropogenic activities had paradoxical effects on the species. On one hand, urbanization, traffic and intensive agriculture had a strong detrimental effect on tree frog distribution. On the other hand, large tree frog populations were frequently associated with gravel pits and military training grounds. Our results allowed us to create a habitat suitability map taking into account detrimental landscape elements around ponds (>1100m away from urban areas and >500m away from first class roads). In parallel, we developed a metapopulation model of the European tree frog in order to identify the critical threats to the long term persistence of the species. Our results indicated that suitable pond density is at the low end of the species requirements. Pond creation must therefore be considered an essential complementary approach to pond conservation and restoration. Our model also provided a mapping solution permitting the location of the must suitable area for pond creation from a metapopulation perspective. As many other amphibians, the European tree frog is not only exposed to an aquatic habitat (breeding and larval period), but also to a terrestrial stage (summer and overwintering habitats). Unfortunately, animals in their terrestrial phase are less conspicuous and, as a consequence, their terrestrial needs are relatively unknown. Using a recent tracking method (the Harmonic Direction Finder), we followed post-breeding frogs and identified favored terrestrial habitats, thus providing another practical conservation tool. We conclude that only the combination of multiple spatially explicit approaches (landscape-scale habitat suitability, metapopulation dynamics and terrestrial needs) is likely to provide wildlife managers with effective tools for the conservation of highly endangered amphibians. Résumé: Les amphibiens font partie des animaux les plus vulnérables du monde. Un tiers des espèces est actuellement menacé d'extinction. Dans l'ancien monde, la disparition des habitats constitue la principale menace pour les grenouilles, crapauds, tritons et salamandres. Dans les paysages fortement urbanisés comme le Plateau Suisse, la plupart des espèces souffrent d'une réduction et d'une fragmentation de leurs habitats. Parmi toutes les espèces indigènes, la rainette verte (Hyla arborea L., 1758) est l'une des plus menacée. Sa distribution a régressé de manière alarmante durant le siècle passé et sa survie régionale à long terme n'est pas assurée. Nous avons développé une méthode de suivi des populations se basant sur le comptage des mâles chanteurs durant la période de reproduction. Cette méthode requiert plusieurs visites à chaque plan d'eau de manière à déterminer précisément la distribution de l'espèce. Nos résultats démontrent que 3 visites par population dans des conditions climatiques favorable (température >20°C) permettent d'obtenir des données de présence/ absence valables. Sur la base de nos comptages sur la Côte lémanique, nous avons analysé les exigences de l'espèce concernant ses sites de reproduction. Il est apparu que les activités humaines avaient un effet paradoxal sur l'espèce. D'une part, l'urbanisation, le trafic routier et l'intensification de l'agriculture ont un effet fortement préjudiciable, tandis que d'autre part les plus grandes populations sont souvent associées à des gravières et autres places d'armes. Nos résultats ont permis de créer une carte de qualité d'habitat prenant en compte les éléments paysagers préjudiciables à la rainette (situé à plus de 1100m de zones urbaines et à plus de 500m de routes de première classe). En parallèle, nous avons développé un modèle métapopulationnel (incluant l'ensemble des populations) de manière à identifier les menaces prépondérantes sur la survie à long terme de l'espèce. Nos résultats ont permis de déterminer que la densité actuelle de plans d'eau adéquats est à la limite inférieure des exigences de l'espèce. La création d'étangs doit donc être considérée comme une approche indispensable et complémentaire à la protection et à la restauration des sites existants. Notre modèle a également fourni des résultats cartographiables permettant l'identification des sites les plus appropriés dans une perspective métapopulationnelle. Comme de nombreux autres amphibiens, la rainette verte est exposée à un habitat aquatique (reproduction et développement larvaire) ainsi qu'à un habitat terrestre (été et hiver). Les animaux étant particulièrement cryptiques dans cette seconde phase, leurs besoins terrestres sont relativement mal connus. Nous avons donc développé une nouvelle méthode de télémétrie basée sur le goniomètre harmonique. Cette méthode nous a permis de suivre des rainettes dans leurs migrations jusqu'à leurs habitats d'été et d'établir ainsi des recommandations pratiques pour la conservation de la rainette. Nous concluons que la combinaison de multiples approches spatialement explicites (qualité d'habitat, dynamique de métapopulation et habitats terrestres) est seule à même de produire des outils efficaces pour la conservation des espèces menacées d'amphibiens.

Remediation of a marine shore tailings deposit and the importance of water-rock interaction on element cycling in the coastal aquifer

We present the study of the geochemical processes associated with the first successful remediation of a marine shore tailings deposit in a coastal desert environment (Bahia de Ite, in the Atacama Desert of Peru). The remediation approach implemented a wetland on top of the oxidized tailings. The site is characterized by a high hydrauliz gradient produced by agricultural irrigation on upstream gravel terraces that pushed river water (similar to 500 mg/L SO(4)) toward the sea and through the tailings deposit. The geochemical and isotopic (delta(2)H(water) and delta(18)O(water), delta(34)S(sulfate) , delta(18)O(sulfate)) approach applied here revealed that evaporite horizons (anhydrite and halite) in the gravel terraces are the source of increased concentrations of SO(4), Cl, and Na up to similar to 1500 mg/L in the springs at the base of the gravel terraces. Deeper groundwater interacting with underlying marine sequences increased the concentrations of SO(4), Cl, and Na up to 6000 mg/L and increased the alkalinity up to 923 mg/L CaCO(3) eq. in the coastal aquifer. These waters infiltrated into the tailings deposit at the shelf-tailings interface. Nonremediated tailings had a low-pH oxidation zone (pH 1-4) with significant accumulations of efflorescent salts (10-20 cm thick) at the surface because of upward capillary transport of metal cations in the arid climate. Remediated tailings were characterized by neutral pH and reducing conditions (pH similar to 7, Eh similar to 100 mV). As a result, most bivalent metals such as Cu, Zn, and Ni had very low concentrations (around 0.01 mg/L or below detection limit) because of reduction and sorption processes. In contrast, these reducing conditions increased the mobility of iron from two sources in this system: (1) The originally Fe(III)-rich oxidation zone, where Fe(II) was reduced during the remediation process and formed an Fe(II) plume, and (2) reductive dissolution of Fe(III) oxides present in the original shelf lithology formed an Fe-Mn plume at 10-m depth. These two Fe-rich plumes were pushed toward the shoreline where more oxidizing and higher pH conditions triggered the precipitation of Fe(HI)hydroxide coatings on silicates. These coatings acted as a filter for the arsenic, which naturally infiltrated with the river water (similar to 500 mu g/L As natural background) into the tailings deposit.

Genetic and maternal contributions to offspring viability under different environmental conditions in various salmonids

Summary : Due to anthropogenic impacts and natural fluctuations, fish usually have to cope with constantly changing and often hostile environments. Whereas adult fish have various possibilities to counteract unfavourable environmental conditions, embryos have much fewer options. Besides by their developing immune system, they are protected by the egg envelopes and several immune substances provided by their mothers. In addition to this, they may also adjust their hatching timing in reaction to various risks. However, individuals may vary in their defensive potential. This variation may be either based on their genetics and/or on differential maternal investments and may be dependent on the experienced stress. Nevertheless, in fish, the impact of such parental contributions on embryo and/or juvenile viability is still poorly investigated. The main objective of this thesis was to investigate the importance of paternal (i.e. genetic) and maternal (i.e. genetic + egg investment) contributions to offspring viability under different environmental conditions and at different life stages. In order to investigate this, we used gametes of various salmonids for in vitro fertilisation experiments based on full-factorial breeding designs. The individual studies are summarised in the following chapters: In the first chapter, we tested the effectiveness of the embryonic immune system in Lake whitefish (Coregonus palaea). Namely, we investigated paternal and maternal contributions to the embryos' tolerance to different kinds of pathogen exposure. Additionally, we tested whether an early sub-léthal exposure has a positive or a negative effect on an embryo's susceptibility to later pathogen exposures with the same pathogen. We found that pre-challenged embryos were more susceptible to future challenges. Moreover, pathogen susceptibility was dependent on maternal investments and/or the embryos' own genetics, depending on the challenge level. Chapter 2 summarises a similar study with brown trout (Salmo trutta). In addition to the previously described investigations, we analysed if genetic effects on offspring viability are mediated either by parental MHC genotypes or relatedness based on neutral microsatellite markers, and we tested if males signal their genetic quality either by their body size or their melanin-based skin colouration. We found that embryo survival was lower at higher stress levels and dependent on the embryos' genetics. Addirionally, parents with similar and/or, very common MHC genotypes had higher offspring viabilities. Finally, darker males produced more viable offspring. In the first two chapters we investigated the embryos' defensive potential based on their immune system, i.e. their pathogen tolerance. In chapter 3 we investigate whether hatching timing of Lake whitefìsh (C. palaea) is dependent on parental contributions and/or on pathogen pressure, and whether there are parental-environmental interactions. We found that whitefish embryos hatch earlier under increasing pathogen pressure. Moreover, hatching timing was affected by embryo genetics and/or maternally provided resources, but the magnitude of the effect was dependent on the pathogen. pressure. We also found a significant paternal-environmental interaction, indicating that the hatching efficiency of a certain sib group is dependent on the pathogen environment. Chapter 4 describes an analogous study with brown trout (S. trutta), with similar findings. In the former chapters, we only looked at offspring performance during the embryonic period, and only under semi-natural conditions. In chapter 5 we now test the performance and viability of embryonic and juvenile brown trout (S. trutta) under natural conditions. To measure embryo viability, we put them in brood boxes, buried them in the gravel of a natural river, and analysed survival after several months. To investigate juvenile survival and performance, wé reared embryos under different stress levels in the laboratory and subsequently released the resulting hatchlings in to a closed river section. Juvenile size and survival was then determined one year later. Additionally, we investigated if sires differ in their genetic quality, determined by embryo and juvenile survival as well as juvenile size, and if they signal their quality by either body size or melanin-based body darkness. We found hat juvenile size was dependent on genetic effects and on maternal investment, whereas this was neither the case for embryo nor for juvenile survival. Additionally, we found that offspring of darker males grew larger, and larger juveniles had also an increased survival. Finally, we found acarry-over effect of the early non-lethal challenge: exposing embryos to higher stress levels resulted in smaller juveniles. To evaluate the long-term performance of differently treated groups, mark-recapture studies are inevitable. For this purpose, effective mass-marking techniques are essential. In chapter 6 we tested the suitability of the fluorescent pigment spray marking method for the mass marking of European graylings (Thymallus thymallus), with very promising results. Our in vitro fertilisation studies on whitefish may reveal new insights on potential genetic benefits of mate choice, but the mating system of whitefish under natural conditions is still poorly investigated. In order to study this, we installed underwater cameras at the spawning place of a Coregonus suidteri population, recorded the whole mating period and subsequently analysed the recordings. Confirmations of previous findings as well as exciting new observations are listed and discussed in chapter 7. Dus aux impacts anthropogéniques et aux fluctuations naturelles, les poissons doivent faire face à des environnements en perpétuel changement. Ces changements font que les poissons doivent s'adapter à de nouvelles situations, souvent hostiles pour eux. Les adultes ont différentes possibilités d'échapper à un environnement peu favorable, ce n'est par contre pas le cas des embryons. Les embryons sont protégés d'une part par leur système immunitaire en développement, d'autre part, par la coquille de l'eeuf et différentes substances immunitaires fournies par leur mère. De plus, ils sont capables d'influencer leur propre date d'éclosion en réponse à différents facteurs de stress. Malgré tout, les individus varient dans leur capacité à se défendre. Cette variation peut être basé sur des facteurs génétiques et/ou sur des facteurs maternels, et est dépendante du stress subi. Néanmoins, chez les poissons, l'impact de telles contributions parentales sur la survie d'embryons et/ou juvéniles est peu étudié. L'objectif principal de cette thèse a été d'approfondir les connaissances sur l'importance de la contribution paternelle (c.a.d. génétique) et maternelle (c.a.d. génétique + investissement dans l'oeuf) sur la survie des jeunes dans différentes conditions expérimentales et stades de vie. Pour faire ces analyses, nous avons utilisé des gamètes de divers salmonidés issus de croisements 'full-factorial'. Les différentes expériences sont résumées dans les chapitres suivants: Dans le premier chapitre, nous avons testé l'efficacité du système immunitaire des embryons chez les corégones (Coregonus palea). Plus précisément nous avons étudié la contribution paternelle et maternelle à la tolérance des embryons à différents niveaux de stress pathogène. Nous avons aussi testé, si une première exposition non létale à un pathogène avait un effet positif ou négatif sur la susceptibilité d'un embryon a une deuxième exposition au même pathogène. Nous avons trouvé que des embryons qui avaient été exposés une première fois étaient plus sensibles au pathogène par la suite. Mais aussi que la sensibilité au pathogène était dépendante de l'investissement de la mère et/ou des gènes de l'embryon, dépendamment du niveau de stress. Le deuxième chapitre résume une étude similaire avec des truites (Salmo truffa). Nous avons examiné, si la survie des jeunes variait sous différentes intensités de stress, et si la variance observée était due aux gènes des parents. Nous avons aussi analysé si les effets génétiques sur la survie des juvéniles étaient dus au MHC (Major Histocompatibility Complex) ou au degré de parenté des parents. De plus, nous avons analysé si les mâles signalaient leur qualité génétique par la taille du corps ou par leur coloration noire, due à la mélanine. On a trouvé que la survie des embryons était plus basse quand le niveau de stress était plus haut mais que la variation restait dépendante de la génétique des embryons. De plus, les parents avec des MHC similaires et/ou communs avaient des embryons avec une meilleure survie. Par contre, des parents avec un degré de parenté plus haut produisent des embryons avec une survie plus mauvaise. Finalement nous avons montré que les mâles plus foncés ont des embryons qui survivent mieux, mais que la taille des mâles n'a pas d'influence sur la survie de ces mêmes embryons. Dans les deux premiers chapitres, nous avons étudié le potentiel de défense des embryons basé sur leur système immunitaire, c.a.d. leur tolérance aux pathogènes. Dans le troisième chapitre, nous nous intéressons à la date d'éclosion des corégones (C. palea), pour voir si elle est influencée par les parents ou par la pression des pathogènes, et si il y a une interaction entre ces deux facteurs. Nous avons trouvé que les jeunes naissent plus rapidement lorsque la pression en pathogènes augmente. La date d'éclosion est influencée par la génétique des embryons et/ou l'investissement des parents, mais c'est la magnitude des effets qui est dépendante de la pression du pathogène. Nous avons aussi trouvé une interaction entre l'effet paternel et l'environnement, ce qui indique que la rapidité d'éclosion de certains croisements est dépendante des pathogènes dans l'environnement. Le chapitre 4 décrit une étude analogue avec de truites (S. truffa), avec des résultats sitzimilaires. Dans les précédents chapitres nous nous sommes uniquement concentrés sur les performances des jeunes durant leur stade embryonnaire, et seulement dans des conditions semi naturelles. Dans le chapitre 5 nous testons la performance et la viabilité des embryons et de juvéniles de truites (S. truffa) dans des conditions naturelles. Nous avons trouvé que la taille des juvéniles était dépendante d'effets génétiques et de l'investissement maternel, mais ceci n'était ni les cas pour les survie des embryons et des juvéniles. De plus, nous avons trouvé que les jeunes des mâles plus foncés devenaient plus grands et que les grands ont un meilleur taux de survie. Finalement nous avons trouvé un 'carry-over effect' d'une première exposition non létale à un pathogène: exposer des embryons à des plus hauts niveaux de stress donnait des juvéniles plus petits. Pour évaluer la performance à long terme de groupes traités dé manières différentes, une méthode de marquage-recapture est inévitable. Pour cette raison, des techniques de marquage en masse sont nécessaires. Dans le chapitre 6, nous avons testé l'efficacité de la technique `fluorescent pigment spray marking' pour le marquage en masse de l'Ombre commun (Thymallus thymallus), avec des résultats très prometteurs. Les études de fertilisations in vitro avec les corégones nous donnent une idée du potentiel bénéfice génétique que représente la sélection d'un bon partenaire, même si le système d'accouplement des corégones en milieu naturel reste peu connu. Pour combler cette lacune, nous avons installé des caméras sous-marines autour de la frayère d'une population de corégones (C. suidteri), nous avons enregistré toute la période de reproduction et nous avons analysé les données par la suite. Ainsi, nous avons été capables de confirmer bien des résultats trouvés précédemment, mais aussi de faire de nouvelles observations. Ces résultats sont reportés dans le septième chapitre, où elles sont comparées avec des observations antérieures.

Premières nidifications du Guêpier d'Europe (Merops apiaster) dans le Canton de Vaud et problématique de conservation du site

The first breeding records of the Bee-eater in canton de Vaud was recorded in a temporarily disused gravel-pit in the Venoge valley (W Switzerland). In 1996 the reproduction of a single pair failed after the eclosion of the chicks during a long period of uninterrupted rain. In 1997, probably the same pair returned and raised 3 fledglings. In 1998, 4 pairs bred raising 20 feldglings. The site is particularly rich with insects and birds: 12 bird species of the red data list breed in the site and surrounding area. However, it is threatened by state project: the gravel-pit is planned to be filled in the coming years. Its future management in relation to the preservation of the Bee-eater and the other bird species on the red data list is presented and discussed.

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.

Topographic forcing of flow partition and flow structures at river bifurcations

This paper presents the predicted flow dynamics from the application of a Reynolds-averaged NavierStokes model to a series of bifurcation geometries with morphologies measured during previous flume experiments. The topography of the bifurcations consists of either plane or bedform-dominated beds which may or may not possess discordance between the two bifurcation distributaries. Numerical predictions are compared with experimental results to assess the ability of the numerical model to reproduce the division of flow into the bifurcation distributaries. The hydrodynamic model predicts: (1) diverting fluxes in the upstream channel which direct water into the distributaries; (2) super-elevation of the free surface induced at the bifurcation edge by pressure differences; and (3) counter-rotating secondary circulation cells which develop upstream of the apex of the bifurcation and move into the downstream channels, with water converging at the surface and diverging at the bed. When bedforms are not present, weak transversal fluxes characterize the upstream channel for almost its entire length, associated with clearly distinguishable secondary circulation cells, although these may be under-estimated by the turbulence model used in the solution. In the bedform dominated case, the same hydrodynamic conditions were not observed, with the bifurcation influence restricted and depth scale secondary circulation cells not forming. The results also demonstrate the dominant effect bed discordance has upon flow division between the two distributaries. Finally, results indicate that in bedform dominated rivers. Consequently, we suggest that sand-bed river bifurcations are more likely to have an influence that extends much further upstream and have a greater impact upon water distribution. This may contribute to observed morphological differences between sand-bedded and gravel-bedded braided river networks. Copyright (C) 2012 John Wiley & Sons, Ltd.

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.

Zonation for 3D aquifer characterization based on joint inversions of multimethod crosshole geophysical data

Predictive groundwater modeling requires accurate information about aquifer characteristics. Geophysical imaging is a powerful tool for delineating aquifer properties at an appropriate scale and resolution, but it suffers from problems of ambiguity. One way to overcome such limitations is to adopt a simultaneous multitechnique inversion strategy. We have developed a methodology for aquifer characterization based on structural joint inversion of multiple geophysical data sets followed by clustering to form zones and subsequent inversion for zonal parameters. Joint inversions based on cross-gradient structural constraints require less restrictive assumptions than, say, applying predefined petro-physical relationships and generally yield superior results. This approach has, for the first time, been applied to three geophysical data types in three dimensions. A classification scheme using maximum likelihood estimation is used to determine the parameters of a Gaussian mixture model that defines zonal geometries from joint-inversion tomograms. The resulting zones are used to estimate representative geophysical parameters of each zone, which are then used for field-scale petrophysical analysis. A synthetic study demonstrated how joint inversion of seismic and radar traveltimes and electrical resistance tomography (ERT) data greatly reduces misclassification of zones (down from 21.3% to 3.7%) and improves the accuracy of retrieved zonal parameters (from 1.8% to 0.3%) compared to individual inversions. We applied our scheme to a data set collected in northeastern Switzerland to delineate lithologic subunits within a gravel aquifer. The inversion models resolve three principal subhorizontal units along with some important 3D heterogeneity. Petro-physical analysis of the zonal parameters indicated approximately 30% variation in porosity within the gravel aquifer and an increasing fraction of finer sediments with depth.

Constraining 3-D electrical resistance tomography with GPR reflection data forimproved aquifer characterization

Surface-based ground penetrating radar (GPR) and electrical resistance tomography (ERT) are common tools for aquifer characterization, because both methods provide data that are sensitive to hydrogeologically relevant quantities. To retrieve bulk subsurface properties at high resolution, we suggest incorporating structural information derived from GPR reflection data when inverting surface ERT data. This reduces resolution limitations, which might hinder quantitative interpretations. Surface-based GPR reflection and ERT data have been recorded on an exposed gravel bar within a restored section of a previously channelized river in northeastern Switzerland to characterize an underlying gravel aquifer. The GPR reflection data acquired over an area of 240×40 m map the aquifer's thickness and two internal sub-horizontal regions with different depositional patterns. The interface between these two regions and the boundary of the aquifer with then underlying clay are incorporated in an unstructured ERT mesh. Subsequent inversions are performed without applying smoothness constraints across these boundaries. Inversion models obtained by using these structural constraints contain subtle resistivity variations within the aquifer that are hardly visible in standard inversion models as a result of strong vertical smearing in the latter. In the upper aquifer region, with high GPR coherency and horizontal layering, the resistivity is moderately high (N300 Ωm). We suggest that this region consists of sediments that were rearranged during more than a century of channelized flow. In the lower low coherency region, the GPR image reveals fluvial features (e.g., foresets) and generally more heterogeneous deposits. In this region, the resistivity is lower (~200 Ωm), which we attribute to increased amounts of fines in some of the well-sorted fluvial deposits. We also find elongated conductive anomalies that correspond to the location of river embankments that were removed in 2002.