Binding of ²³⁹Pu and ⁹⁰Sr to organic colloids in soil solutions: evidence from a field experiment.

Colloidal transport has been shown to enhance the migration of plutonium in groundwater downstream from contaminated sites, but little is known about the adsorption of ⁹⁰Sr and plutonium onto colloids in the soil solution of natural soils. We sampled soil solutions using suction cups, and separated colloids using ultrafiltration to determine the distribution of ²³⁹Pu and ⁹⁰Sr between the truly dissolved fraction and the colloidal fraction of the solutions of three Alpine soils contaminated only by global fallout from the nuclear weapon tests. Plutonium was essentially found in the colloidal fraction (>80%) and probably associated with organic matter. A significant amount of colloidal ⁹⁰Sr was detected in organic-rich soil solutions. Our results suggest that binding to organic colloids in the soil solutions plays a key role with respect to the mobility of plutonium in natural alpine soils and, to a lesser extent, to the mobility of ⁹⁰Sr.

Early Cretaceous life, climate and anoxia

Early Cretaceous life and the environment were strongly influenced by the accelerated break up of Pangaea, which was associated with the formation of a multitude of rift basins, intensified spreading, and important volcanic activity on land and in the sea. These processes likely interacted with greenhouse conditions, and Early Cretaceous climate oscillated between "normal" greenhouse, predominantly arid conditions, and intensified greenhouse, predominantly humid conditions. Arid conditions were important during the latest Jurassic and early Berriasian, the late Barremian, and partly also during the late Aptian. Humid conditions were particularly intense and widespread during shorter episodes of environmental change (EECs): the Valanginian Weissert, the latest Hauterivian Faraoni, the latest Barremian earliest Aptian Taxy, the early Aptian Selli, the early late Aptian Fallot and the late Aptian-early Albian Paquier episodes. Arid conditions were associated with evaporation, low biogeochemical weathering rates, low nutrient fluxes, and partly stratified oceans, leading to oxygen depletion and enhanced preservation of laminated, organic-rich mud (LOM). Humid conditions enabled elevated biogeochemical weathering rates and nutrient fluxes, important runoff and the buildup of freshwater lids in proximal basins, intensified oceanic and atmospheric circulation, widespread upwelling and phosphogenesis, important primary productivity and enhanced preservation of LOM in expanded oxygen-minimum zones. The transition of arid to humid climates may have been associated with the net transfer of water to the continent owing to the infill of dried-out groundwater reservoirs in internally drained inland basins. This resulted in shorter-term sea-level fall, which was followed by sea-level rise. These sea-level changes and the influx of freshwater into the ocean may have influenced oxygen-isotope signatures. Climate change preceding and during the Early Cretaceous EECs may have been rapid, but in general, the EECs had a "pre"-history, during which the stage was set for environmental change. Negative feedback on the climate through increased marine LOM preservation was unlikely, because of the low overall organic-carbon accumulation rates during these episodes. Life and climate co-evolved during the Early Cretaceous. Arid conditions may have affected continental life, such as across the Tithonian/Berriasian boundary. Humid conditions and the corresponding tendency to develop dys- to anaerobic conditions in deeper ocean waters led to phases of accelerated extinction in oceans, but may have led to more luxuriant vegetation cover on continents, such as during the Valanginian, to the benefit of herbivores. During Early Cretaceous EECs, reef systems and carbonate platforms in general were particularly vulnerable. They were the first to disappear and the last to recover, often only after several million years. (C) 2011 Elsevier Ltd. All rights reserved.

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.

Development of bacteria-based bioassays for arsenic detection in natural waters.

Arsenic contamination of natural waters is a worldwide concern, as the drinking water supplies for large populations can have high concentrations of arsenic. Traditional techniques to detect arsenic in natural water samples can be costly and time-consuming; therefore, robust and inexpensive methods to detect arsenic in water are highly desirable. Additionally, methods for detecting arsenic in the field have been greatly sought after. This article focuses on the use of bacteria-based assays as an emerging method that is both robust and inexpensive for the detection of arsenic in groundwater both in the field and in the laboratory. The arsenic detection elements in bacteria-based bioassays are biosensor-reporter strains; genetically modified strains of, e.g., Escherichia coli, Bacillus subtilis, Staphylococcus aureus, and Rhodopseudomonas palustris. In response to the presence of arsenic, such bacteria produce a reporter protein, the amount or activity of which is measured in the bioassay. Some of these bacterial biosensor-reporters have been successfully utilized for comparative in-field analyses through the use of simple solution-based assays, but future methods may concentrate on miniaturization using fiberoptics or microfluidics platforms. Additionally, there are other potential emerging bioassays for the detection of arsenic in natural waters including nematodes and clams.

Biogeochemical evolution of a marine shore tailings deposit during bioremediation

Summary : Mining activities produce enormous amounts of waste material known as tailings which are composed of fine to medium size particles. These tailings often contain sulfides, which oxidation can lead to acid and metal contamination of water; therefore they need to be remediated. In this work a tailings bioremediation approach was investigated by an interdisciplinary study including geochemistry, mineralogy and microbiology. The aim of the work was to study the effect of the implementation of wetland above oxidizing tailings on the hydrogeology and the biogeochemical element cycles, and to assess the system evolution over time. To reach these goals, biogeochemical processes occurring in a marine shore tailings deposit were investigated. The studied tailings deposit is located at the Bahìa de Ite, Pacific Ocean, southern Peru, where between 1940 and 1996 the tailings were discharged from the two porphyry copper mines Cuajone and Toquepala. After the end of deposition, a remediation approach was initiated in 1997 with a wetland implementation above the oxidizing tailings. Around 90% of the tailings deposits (total 16 km2) were thus remediated, except the central delta area and some areas close to the shoreline. The multi-stable isotope study showed that the tailings were saturated with fresh water in spite of the marine setting, due to the high hydraulic gradient resulting from the wetland implementation. Submarine groundwater discharge (SGD) was the major source of SO4 2-, C1-, Na+, Fe2+, and Mn2+ input into the tailings at the original shelf-seawater interface. The geochemical study (aquatic geochemistry and X-Ray diffraction (XRD) and sequential extractions from the solid fraction) showed that iron and sulfur oxidation were the main processes in the non-remediated tailings, which showed a top a low-pH oxidation zone with strong accumulation of efflorescent salts at the surface due to capillary upward transport of heavy metals (Fe, Cu, Zn, Mn, Cd, Co, and Ni) in the arid climate. The study showed also that the implementation of the wetland resulted in very low concentrations of heavy metals in solution (mainly under the detection limit) due to the near neutral pH and more reducing conditions (100-150 mV). The heavy metals, which were taken from solution, precipitated as hydroxides and sulfides or were bound to organic matter. The bacterial community composition analysis by Terminal Restriction Fragment Length Polymorphism (T-RFLP) and cloning and sequencing of 16S rRNA genes combined with a detailed statistical analysis revealed a high correlation between the bacterial distribution and the geochemical variables. Acidophilic autotrophic oxidizing bacteria were dominating the oxidizing tailings, whereas neutrophilic and heterotrophic reducing bacteria were driving the biogeochemical processes in the remediated tailings below the wetland. At the subsurface of the remediated tailings, an iron cycling was highlighted with oxidation and reduction processes due to micro-aerophilic niches provided by the plant rhizosphere in this overall reducing environment. The in situ bioremediation experiment showed that the main parameter to take into account for the effectiveness was the water table and chemistry which controls the system. The constructed remediation cells were more efficient and rapid in metal removal when saturation conditions were available. This study showed that the bioremediation by wetland implementation could be an effective and rapid treatment for some sulfidic mine tailings deposits. However, the water saturation of the tailings has to be managed on a long-term basis in order to guarantee stability. Résumé : L'activité minière produit d'énormes quantités de déchets géologiques connus sous le nom de « tailings » composées de particules de taille fine à moyenne. Ces déchets contiennent souvent des sulfures dont l'oxydation conduit à la formation d'effluents acides contaminés en métaux, d'où la nécessité d'effectuer une remédiation des sites de stockage concernés. Le but de ce travail est dans un premier temps d'étudier l'effet de la bio-remédiation d'un dépôt de tailings oxydés sur l'hydrogéologie du système et les cycles biogéochimiques des éléments et en second lieu, d'évaluer l'évolution du processus de remédiation dans le temps. Le site étudié dans ce travail est situé dans la Bahía de Ite, au sud du Pérou, au bord de l'Océan Pacifique. Les déchets miniers en question sont déposés dans un environnement marin. De 1940 à 1996, les déchets de deux mines de porphyre cuprifère - Cuajone et Toquepala - ont été acheminés sur le site via la rivière Locumba. En 1997, une première remédiation a été initiée avec la construction d'une zone humide sur les tailings. Depuis, environ 90% de la surface du dépôt (16 km2) a été traité, les parties restantes étant la zone centrale du delta du Locumba et certaines zones proches de la plage. Malgré la proximité de l'océan, les études isotopiques menées dans le cadre de ce travail ont montré que les tailings étaient saturés en eau douce. Cette saturation est due à la pression hydraulique résultant de la mise en place des zones humides. Un écoulement d'eau souterrain sous-marin a été à détecté à l'interface entre les résidus et l'ancien fond marin. En raison de la géologie locale, il constitue une source d'entrée de SO4 2-, Cl-, Na+, FeZ+, et Mn2+ dans le système. L'analyse de la géochimie aquatique, la Diffraction aux Rayons X (XRD) et l'extraction séquentielle ont montré que l'oxydation du fer et .des sulfures est le principal processus se produisant dans les déchets non remédiés. Ceci a entraîné le développement d'une zone d'oxydation à pH bas induisant une forte accumulation des sels efflorescents, conséquence de la migration capillaire des métaux lourds (Fe, Cu, Zn, Mn, Cd, Co et Ni) de la solution vers la surface dans ce climat aride. Cette étude a montré également que la construction de la zone humide a eu comme résultats une précipitation des métaux dans des phases minérales en raison du pH neutre et des conditions réductrices (100-150mV). Les métaux lourds ont précipité sous la forme d'hydroxydes et de sulfures ou sont adsorbés à la matière organique. L'analyse de la composition de la communauté bactérienne à l'aide la technique T-RFLP (Terminal Restriction Fragment Length Polymorphism) et par le clonage/séquençage des gènes de l'ARNr 16S a été combinée à une statistique détaillée. Cette dernière a révélé une forte corrélation entre la distribution de bactéries spécifiques et la géochimie : Les bactéries autotrophes acidophiles dominent dans les déchets oxydés non remédiés, tandis que des bactéries hétérotrophes neutrophiles ont mené les processus microbiens dans les déchets remédiés sous la zone humide. Sous la surface de la zone humide, nos analyses ont également mis en évidence un cycle du fer par des processus d'oxydoréduction rendus possibles par la présence de niches micro-aérées par la rhizosphère dans cet environnement réducteur. L'expérience de bio-remédiation in situ a montré que les paramètres clés qui contrôlent l'efficacité du traitement sont le niveau de la nappe aquifère et la chimie de l'eau. Les cellules de remédiation se sont montrées plus efficaces et plus rapides lorsque le système a pu être saturé en eau. Finalement, cette étude a montré que la bio-remédiation de déchets miniers par la construction de zones humides est un moyen de traitement efficace, rapide et peu coûteux. Cependant, la saturation en eau du système doit être gérée sur le long terme afin de garantir la stabilité de l'ensemble du système.

Focused time-lapse inversion of radio and audio magnetotelluric data

Geoelectrical techniques are widely used to monitor groundwater processes, while surprisingly few studies have considered audio (AMT) and radio (RMT) magnetotellurics for such purposes. In this numerical investigation, we analyze to what extent inversion results based on AMT and RMT monitoring data can be improved by (1) time-lapse difference inversion; (2) incorporation of statistical information about the expected model update (i.e., the model regularization is based on a geostatistical model); (3) using alternative model norms to quantify temporal changes (i.e., approximations of l(1) and Cauchy norms using iteratively reweighted least-squares), (4) constraining model updates to predefined ranges (i.e., using Lagrange Multipliers to only allow either increases or decreases of electrical resistivity with respect to background conditions). To do so, we consider a simple illustrative model and a more realistic test case related to seawater intrusion. The results are encouraging and show significant improvements when using time-lapse difference inversion with non l(2) model norms. Artifacts that may arise when imposing compactness of regions with temporal changes can be suppressed through inequality constraints to yield models without oscillations outside the true region of temporal changes. Based on these results, we recommend approximate l(1)-norm solutions as they can resolve both sharp and smooth interfaces within the same model. (C) 2012 Elsevier B.V. All rights reserved.

Effect of irrigation water on the incidence of Salmonella spp. on lettuces produced by urban agriculture and sold on the markets in Dakar, Senegal

The aim of our survey was to assess the effect of irrigation water of the microbiological quality on the production chain of lettuce in the Dakar area. Microbiological analysis showed that 35% of irrigation water was contaminated by Salmonella spp. between the two water-types used for irrigation (groundwater and wastewater), no significant difference (p>0.05) in their degree of contamination was found. The incidence of different types of irrigation water on the contamination rate of lettuces from the farm (Pikine and Patte d'Oie) was not different either (p>0.05). However, the contamination rate of lettuce from markets of Dalifort and Grand-Yoff that were supplied by the area of Patte d'Oie was greater than those of Sham and Zinc supplied by Pikine (p<0.05). Comparison of serotypes of Salmonella isolated from irrigation water and lettuce showed that irrigation water may affect the microbiological quality of lettuce. Manures, frequently used as organic amendment in cultivating lettuce are another potential source of contamination. These results showed that lettuce may constitute effective vectors for the transmission of pathogens to consumers. Extensive treatment of the used wastewater and/or composting of manure could considerably reduce these risks.

Integration of high-resolution geophysical data for the characterization of alluvial aquifers : new approaches and implications for hydrological predictions

Abstract Accurate characterization of the spatial distribution of hydrological properties in heterogeneous aquifers at a range of scales is a key prerequisite for reliable modeling of subsurface contaminant transport, and is essential for designing effective and cost-efficient groundwater management and remediation strategies. To this end, high-resolution geophysical methods have shown significant potential to bridge a critical gap in subsurface resolution and coverage between traditional hydrological measurement techniques such as borehole log/core analyses and tracer or pumping tests. An important and still largely unresolved issue, however, is how to best quantitatively integrate geophysical data into a characterization study in order to estimate the spatial distribution of one or more pertinent hydrological parameters, thus improving hydrological predictions. Recognizing the importance of this issue, the aim of the research presented in this thesis was to first develop a strategy for the assimilation of several types of hydrogeophysical data having varying degrees of resolution, subsurface coverage, and sensitivity to the hydrologic parameter of interest. In this regard a novel simulated annealing (SA)-based conditional simulation approach was developed and then tested in its ability to generate realizations of porosity given crosshole ground-penetrating radar (GPR) and neutron porosity log data. This was done successfully for both synthetic and field data sets. A subsequent issue that needed to be addressed involved assessing the potential benefits and implications of the resulting porosity realizations in terms of groundwater flow and contaminant transport. This was investigated synthetically assuming first that the relationship between porosity and hydraulic conductivity was well-defined. Then, the relationship was itself investigated in the context of a calibration procedure using hypothetical tracer test data. Essentially, the relationship best predicting the observed tracer test measurements was determined given the geophysically derived porosity structure. Both of these investigations showed that the SA-based approach, in general, allows much more reliable hydrological predictions than other more elementary techniques considered. Further, the developed calibration procedure was seen to be very effective, even at the scale of tomographic resolution, for predictions of transport. This also held true at locations within the aquifer where only geophysical data were available. This is significant because the acquisition of hydrological tracer test measurements is clearly more complicated and expensive than the acquisition of geophysical measurements. Although the above methodologies were tested using porosity logs and GPR data, the findings are expected to remain valid for a large number of pertinent combinations of geophysical and borehole log data of comparable resolution and sensitivity to the hydrological target parameter. Moreover, the obtained results allow us to have confidence for future developments in integration methodologies for geophysical and hydrological data to improve the 3-D estimation of hydrological properties.

Stable isotope compositions of mammoth teeth from Niederweningen, Switzerland: Implications for the Late Pleistocene climate, environment, and diet

Oxygen and carbon isotope compositions of well-preserved mammoth teeth from the Middle Wurmian (40-70 ka) peat layer of Niederweningen, the most important mammoth site in Switzerland, were analysed to reconstruct Late Pleistocene palaeoclimatic and palaeoenvironmental conditions. Drinking water (delta(18)O values of approximately -12.3 +/- 0.9 parts per thousand were calculated front oxygen isotope compositions of mammoth tooth enamel apatite using a species-specific calibration for modern elephants. These delta(18)O(H2O) values reflect the mean oxygen isotope composition of the palaeo-precipitation and are similar to those directly measured for fate Pleistocene groundwater from aquifers in northern Switzerland and southern Germany. Using a present-day delta(18)O(H2)o-precipitation-air temperature relation for Switzerland, a mean annual air temperature (MAT) of around 4.3 +/- 2.1 degrees C can be calculated for the Middle Wurmian at this site. This MAT is in good agreement with palaeotemperature estimates on the basis of Middle Wurmian groundwater recharge temperatures and beetle assemblages. Hence, the climatic conditions in this region were around 4 degrees C cooler during the Middle Wurmian interstadial phase, around 45-50ka BP, than they are today. During this period the mammoths from Niederweningen lived in an open tundra-like, C(3) plant-dominated environment as indicated by enamel (delta(13)C values of -11.5 +/- 0.3 parts per thousand and pollen and macroplant fossils found in the embedding peat. The low variability of enamel delta(13)C and delta(18)O values from different mammoth teeth reflects similar environmental conditions and supports a relatively small time frame for the fossil assemblage. (C) 2006 Elsevier Ltd and INQUA. All rights reserved.

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.

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.

Transport of a biocontrol Pseudomonas fluorescens through 2.5-m deep outdoor lysimeters and survival in the effluent water

Application of wild-type or genetically-modified bacteria to the soil environment entails the risk of dissemination of these organisms to the groundwater. To measure vertical transport of bacteria under natural climatic conditions, Pseudomonas fluorescens strain CHA0 was released together with bromide as a mobile tracer at the surface of large outdoor lysimeters. Two experiments, one starting in autumn 1993 and the other in spring 1994 were performed. Shortly after a heavy rainfall in late spring 1994, the released bacteria were detected for the first time in effluent water from the 2.5-m-deep lysimeters in both experiments, i.e. 210 d and 21 d, respectively, after inoculation. Only a 10−9 to 10−8 fraction of the inoculum was recovered as culturable cells in the effluent water, but a larger fraction of the CHA0 cells was in a non-culturable state as detected with immunofluorescence microscopy. As much as 50% of the mobile tracer percolated through the lysimeters, indicating that, compared with bromide, bacterial cells were retained in soil. In the second part of this study, persistence of CHA0 in groundwater microcosms consisting of lysimeter effluent water was studied for 380 d. Survival of the inoculant as culturable cells was better under anaerobic than under aerobic conditions. However, a large fraction of the cells became non-culturable in both cases. When the experiment was performed with filter-sterilized effluent water, the total count of introduced bacteria did not decline with time. In conclusion, the biocontrol strain was transported in low numbers to a potential groundwater level under natural climatic conditions, but could persist for an extended period in groundwater microcosms.

Local and global error models to improve uncertainty quantification

In groundwater applications, Monte Carlo methods are employed to model the uncertainty on geological parameters. However, their brute-force application becomes computationally prohibitive for highly detailed geological descriptions, complex physical processes, and a large number of realizations. The Distance Kernel Method (DKM) overcomes this issue by clustering the realizations in a multidimensional space based on the flow responses obtained by means of an approximate (computationally cheaper) model; then, the uncertainty is estimated from the exact responses that are computed only for one representative realization per cluster (the medoid). Usually, DKM is employed to decrease the size of the sample of realizations that are considered to estimate the uncertainty. We propose to use the information from the approximate responses for uncertainty quantification. The subset of exact solutions provided by DKM is then employed to construct an error model and correct the potential bias of the approximate model. Two error models are devised that both employ the difference between approximate and exact medoid solutions, but differ in the way medoid errors are interpolated to correct the whole set of realizations. The Local Error Model rests upon the clustering defined by DKM and can be seen as a natural way to account for intra-cluster variability; the Global Error Model employs a linear interpolation of all medoid errors regardless of the cluster to which the single realization belongs. These error models are evaluated for an idealized pollution problem in which the uncertainty of the breakthrough curve needs to be estimated. For this numerical test case, we demonstrate that the error models improve the uncertainty quantification provided by the DKM algorithm and are effective in correcting the bias of the estimate computed solely from the MsFV results. The framework presented here is not specific to the methods considered and can be applied to other combinations of approximate models and techniques to select a subset of realizations