Application of environmental groundwater tracers at the Sulphur Bank Mercury Mine, California, USA

Boron, chloride, sulfate, δD, δ¹⁸O, and ³H concentrations in surface water and groundwater samples from the Sulphur Bank Mercury Mine (SBMM), California, USA were used to examine geochemical processes and provide constraints on evaporation and groundwater flow. SBMM is an abandoned sulfur and mercury mine with an underlying hydrothermal system, adjacent to Clear Lake, California. Results for non-³H tracers (i.e., boron, chloride, sulfate, δD, and δ¹⁸O) identify contributions from six water types at SBMM. Processes including evaporation, mixing, hydrothermal water input and possible isotopic exchange with hydrothermal gases are also discerned. Tritium data indicate that hydrothermal waters and other deep groundwaters are likely pre-bomb (before ~1952) in age while most other waters were recharged after ~1990. A boron-based steady-state reservoir model of the Herman Impoundment pit lake indicates that 71–79% of its input is from meteoric water with the remainder from hydrothermal contributions. Results for groundwater samples from six shallow wells over a 6–month period for δD and δ¹⁸O suggests that water from Herman Impoundment is diluted another 3% to more than 40% by infiltrating meteoric water, as it leaves the site. Results for this investigation show that environmental tracers are an effective tool to understand the SBMM hydrogeologic regime.

---

Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.

---

Résumé Les concentrations en bore, chlorure, sulfate, δD, δ¹⁸O, et ³H d’échantillons d’eaux de surface et souterraine prélevés dans le banc de soufre de la mine de mercure (SBMM en anglais) en Californie, USA, ont été utilisées pour étudier les processus géochimiques et pour fournir des contraintes à l’évaporation et à l’écoulement des eaux souterraines. La SBMM est une mine de soufre et de mercure abandonnée, adjacente au lac Clear en Californie et sous laquelle se trouve un système hydrothermal.Les résultats des traceurs autres que le tritium (bore, chlorure, sulfate, δD, et δ¹⁸O) ont permis d’identifier des contributions de six types d’eaux à SBMM. Des processus tels que l’évaporation, le mélange, l’entrée d’eau hydrothermale et de possibles échanges isotopiques avec des gaz hydrothermaux ont également été identifiés. Les données de tritium montrent que les eaux hydrothermales et d’autres eaux profondes sont probablement d’age antérieure à la bombe (avant ~1952), alors que la plupart des autres eaux sont issues de la recharge après ~1990. Un modèle de réservoir représentant le lac situé dans la partie Herman Impoundment de l’ancienne mine, en régime permanent et basé sur le bore, montre que 71–79% de l’eau provient des précipitations, le reste provenant de contributions hydrothermales. Les résultats de δD et δ¹⁸O pour des échantillons d’eau souterraine de six puits peu profonds sur une période de 6 mois suggèrent que l’eau de Hermann Impoundment est encore diluée entre 3% jusqu’à plus de 40% lorsqu’elle quitte le site du fait de l’infiltration d’eau météorique. Les résultats de cette étude montrent que les traceurs environnementaux constituent un outil efficace pour comprendre le régime hydrogéologique de la SBMM.

---

Resumen Se han utilizado datos de concentraciones de boro, cloruros, sulfatos, δD, δ¹⁸O, y ³H en muestras de aguas superficiales y subterráneas procedentes de la Mina de Mercurio Sulphur Bank (SBMM), California, USA, para estudiar los procesos geoquímicos y caracterizar la evaporación y el flujo del agua subterránea. SBMM es una mina abandonada de azufre y mercurio con un sistema hidrotermal subyacente, cercano al Lago Clear, California.Los resultados de los trazadores que no son ³H (por ejemplo, boro, cloruros, sulfatos δD, y δ¹⁸O) identifican las contribuciones de seis tipos de agua en la SBMM. Se han identificado diversos procesos, como evaporación, mezcla, entrada de agua hidrotermal y la posibilidad de intercambios isotópicos con gases hidrotermales. Los datos de tritio indican que las aguas hidrotermales y otras aguas subterráneas profundas son probablemente pre-bombas (previas a ~1952) en edad mientras que la mayoría de las otras aguas se han recargado después de ~1990. Un modelo estacionario basado en el boro de la fosa del lago Herman Impoundment indica que el 71–79% de su entrada procede de agua meteórica con restos de contribución hidrotermal. Los resultados para las muestras de aguas subterráneas procedentes de seis pozos superficiales en un periodo de 6 meses para δD y δ¹⁸O sugieren que el agua procedente de Herman Impoundment se diluye en otro 3% hasta más del 40% debido a la infiltración de agua meteórica., que sale del sitio. Los resultados de esta investigación muestran que los trazadores ambientales son una herramienta efectiva para entender el régimen hidrogeológico de la SBMM.

     

Evaluation of distributed recharge in an upland semi-arid karst system: the West Bank Mountain Aquifer,Middle East

Assessment of recharge in a structurally complex upland karst limestone aquifer situated in a semi-arid environment is difficult. Resort to surrogate indicators such as measurement of spring outflow and borehole discharge, is a common alternative, and attempts to apply conventional soil moisture deficit analysis may not adequately account for the intermittent spate conditions that arise in such environments. A modelling approach has been made using the West Bank Mountain Aquifer system in the Middle East as a trial. The model uses object oriented software which allows various objects to be switched on and off. Each of the main recharge processes identified in the West Bank is incorporated. The model allows either conventional soil moisture deficit analysis calculations or wetting threshold calculations to be made as appropriate, and accommodates both direct recharge and secondary recharge. Daily time steps enable recharge and runoff routing to be calculated for each node. Model runs have enabled a series of simulations for each of the three aquifer basins in the West Bank and for the whole of the West Bank. These provide recharge estimates comparable to those prepared by earlier workers by conventional means. The model is adaptable and has been successfully used in other environments.     

Glacio-lacustrine stratigraphy,aquifer characterization and contaminant transport: a case study in South Lake Tahoe,California, USA

A hydrogeologic model that has been used by many researchers and consultants to describe an area of South Lake Tahoe, California, USA impacted by MTBE contamination describes a relatively homogeneous unconfined aquifer comprised of poorly sorted glacial outwash deposits, within which water-supply wells are able to exert significant alteration in natural groundwater flow. A re-examination of the area’s hydrogeology is presented, which supports a layered heterogeneous aquifer system constructed of alternating fine and coarser-grained glacio-lacustrine depositional units. This re-evaluation was accomplished through a review of lithologic logs across an area of approximately 1 km², combined with observations of significant hydraulic head differences and knowledge of the depositional environments controlled by Pleistocene Lake Tahoe high stands. Many of the fine-grained units observed at depths from 6 to 15 m, although relatively thin, are generally continuous and serve as significant barriers to groundwater flow. The vertical migration of contamination across these fine-grained units to deeper groundwaters was facilitated by cross-screened monitoring wells installed as part of site investigation activities. This conclusion highlights the importance of geologic characterization and proper monitoring well construction at contaminated site investigations.

---

Electronic supplementary material   The online version of this article (doi:) contains supplementary material, which is available to authorized users.

Accurate local upscaling with variable compact multipoint transmissibility calculations

We propose a new single-phase local upscaling method that uses spatially varying multipoint transmissibility calculations. The method is demonstrated on two-dimensional Cartesian and adaptive Cartesian grids. For each cell face in the coarse upscaled grid, we create a local fine grid region surrounding the face on which we solve two generic local flow problems. The multipoint stencils used to calculate the fluxes across coarse grid cell faces involve the six neighboring pressure values. They are required to honor the two generic flow problems. The remaining degrees of freedom are used to maximize compactness and to ensure that the flux approximation is as close as possible to being two-point. The resulting multipoint flux approximations are spatially varying (a subset of the six neighbors is adaptively chosen) and reduce to two-point expressions in cases without full-tensor anisotropy. Numerical tests show that the method significantly improves upscaling accuracy as compared to commonly used local methods and also compares favorably with a local–global upscaling method.     

Two-dimensional modelling of contaminant transport through saturated porous media using the radial point interpolation method (RPIM)

A new numerical tool is presented which models the two-dimensional contaminant transport through saturated porous media using a meshfree method called the radial point interpolation method (RPIM) with polynomial reproduction. In RPIM, an approximate solution is constructed entirely in terms of a set of nodes and no characterisation of the interrelationship of the nodes is needed. An advection-dispersion equation with sorption is considered to illustrate the applicability of the RPIM. The Galerkin weak form of the governing equation is formulated using two-dimensional meshfree shape functions constructed using thin plate spline radial basis functions. A computer program is developed for the implementation of the RPIM procedure. Three numerical examples are presented and the results are compared with those obtained from the analytical solution and finite element method. The experimental results are also used to validate the approach. The proposed RPIM has generated results with no oscillations and they are insensitive to Peclet constraints.     

Application of FLATModel, a 2D finite volume code, to debris flows in the northeastern part of the Iberian Peninsula

FLATModel is a two-dimensional shallow-water approximation code with corrections and modifications that create a simulation tool adapted to debris-flows behaviour. FLATModel uses the finite volume method with the numerical implementation of the Godunov scheme and includes correction terms regarding the effect of flow over high slopes and curvature. Additionally, the stop-and-go phenomenon, the basal entrainment and a correction regarding the front inclination of the final deposit are incorporated into FLATModel. In addition, different flow resistance laws were integrated in the numerical code including Bingham, Herschel–Bulkley and Voellmy fluid model. Firstly, our numerical model was validated using analytical solutions of a dam-break scenario and published data on a laboratory experiment. Secondly, three real events, which occurred in the northeastern part of the Iberian Peninsula, were back-calculated. Although field observations of the three events are not very detailed, the back-analyses revealed interesting patterns on the flow dynamics, and the numerical results generally showed good agreement with field data. Comparing the different flow resistance laws, the Voellmy fluid model presents the best behaviour regarding both the flow behaviour and the deposit characteristics. Preliminary simulation runs incorporating the effect of basal entrainment offered satisfactory results, although the final volume is rather sensitive on the selected friction angle of channel-bed material. The outcomes regarding the correction of the calculated front inclination of the final deposit showed that this implementation strongly improves the simulation results and better represents steep fronts of final deposits.     

Debris-flow event in the Frangerello Stream-Susa Valley (Italy)—calibration of numerical models for the back analysis of the 16 October, 2000 rainstorm

Three debris-flow simulation model software have been applied to the back analysis of a typical alpine debris flow that caused significant deposition on an urbanized alluvial fan. Parameters used in the models were at first retrieved from the literature and then adjusted to fit field evidence. In the case where different codes adopted the same parameters, the same input values were used, and comparable outputs were obtained. Results of the constitutive laws used (Bingham rheology, Voellmy fluid rheology and a quadratic rheology formulation which adds collisional and turbulent stresses to the Bingham law) indicate that no single rheological model appears to be valid for all debris flows. The three applied models appear to be capable of reasonable reproduction of debris-flow events, although with different levels of detail. The study shows how different software can be used to predict the debris-flow motion for various purposes from a first screening, to predict the runout distance and deposition of the solid material and to the different behaviour of the mixtures of flows with variation of maximum solid concentration.     

Coalification in Carboniferous sediments from the Lötschberg base tunnel

Vitrinite reflectance (Rr), proximate analysis and carbon isotope composition (δ¹³C) have been used to characterise coal samples from two zones of Late Carboniferous sediments (Gastern and Ferden) in the Aar massif where they are penetrated by the Lötschberg base tunnel (constructed between 1999 and 2005). Samples are characterised by variable ash yields (21.7–93.9%; dry basis); those with ash yields of less than ~50% and with volatile matter content (V;dry ash-free basis) within the limits 2 < V% ≤ 8 are anthracite. Values of Rr range from 3.89% to 5.17% and indicate coalification to the rank of anthracite and meta-anthracite in both Gastern and Ferden Carboniferous zones. Samples of anthracite and shale from the Gastern Carboniferous exhibit a relatively small range in δ¹³C values (–24.52‰ to –23.38‰; mean: –23.86‰) and are lighter than anthracite samples from the Ferden Carboniferous (mean: –22.20‰). The degree of coalification in the Gastern and Ferden Carboniferous zones primarily depends on the maximum rock temperature (T) attained as a result of burial heating. Vitrinite reflectance based estimates of T range from ~290° –360 °C. For a proposed palaeogeothemal gradient of 25 ° C/km at the time of maximum coalification the required overburden is attributable to relatively thin autochtonous Mesozoic/Cenozoic sedimentary cover of the Aar massif and Gastern granite and deep tectonic burial beneath advancing Helvetic, Ultrahelvetic and Prealpine (Penninic) nappes in Early Oligocene to Miocene.     

Zircon M257 ‐ a Homogeneous Natural Reference Material for the Ion Microprobe U‐Pb Analysis of Zircon

We introduce and propose zircon M257 as a future reference material for the determination of zircon U‐Pb ages by means of secondary ion mass spectrometry. This light brownish, flawless, cut gemstone specimen from Sri Lanka weighed 5.14 g (25.7 carats). Zircon M257 has TIMS‐determined, mean isotopic ratios (2s uncertainties) of 0.09100 ± 0.00003 for ²⁰⁶pb/²³⁸U and 0.7392 ± 0.0003 for ²⁰⁷pb/²³⁵U. Its ²⁰⁶pb/²³⁸U age is 561.3 ± 0.3 Ma (unweighted mean, uncertainty quoted at the 95% confidence level); the U‐Pb system is concordant within uncertainty of decay constants. Zircon M257 contains ～ 840 μg g^?1 U (Th/U ～ 0.27). The material exhibits remarkably low heterogeneity, with a virtual absence of any internal textures even in cathodoluminescence images. The uniform, moderate degree of radiation damage (estimated from the expansion of unit‐cell parameters, broadening of Raman spectral parameters and density) corresponds well, within the “Sri Lankan trends”, with actinide concentrations, U‐Pb age, and the calculated alpha fluence of 1.66 × 10¹⁸ g^?1. This, and a (U+Th)/He age of 419 ± 9 Ma (2s), enables us to exclude any unusual thermal history or heat treatment, which could potentially have affected the retention of radiogenic Pb. The oxygen isotope ratio of this zircon is 13.9%o VSMOW suggesting a metamorphic genesis in a marble or calc‐silicate skarn.     

The Glarus thrust: excursion guide and report of a field trip of the Swiss Tectonic Studies Group (Swiss Geological Society, 14.–16. 09. 2006)

This excursion guide results form a field trip to the Glarus nappe complex organized by the Swiss Tectonic Studies Group in 2006. The aim of the excursion was to discuss old and recent concepts related to the evolution of the Glarus thrust. The major aspects were (i) the interplay between deformation, fluid flow and geochemical alteration, (ii) episodic versus continuous deformation and fluid flow, and (iii) the link between large-scale structures, microstructures, and geochemical aspects. Despite 150 years of research in the Glarus nappe complex and the new results discussed during the excursion, there exist controversies that still are unsolved.