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.

     

A coupling of mixed and discontinuous Galerkin finite-element methods for poroelasticity

In this paper, we formulate a finite-element procedure for approximating the coupled fluid and mechanics in Biot’s consolidation model of poroelasticity. We approximate the flow variables by a mixed finite-element space and the displacement by a family of discontinuous Galerkin methods. Theoretical convergence error estimates are derived and, in particular, are shown to be independent of the constrained specific storage coefficient, c _o . This suggests that our proposed algorithm is a potentially effective way to combat locking, or the nonphysical pressure oscillations, which sometimes arise in numerical algorithms for poroelasticity.     

Hydrogeologic assessment of escalating groundwater exploitation in the Indus Basin, Pakistan

Groundwater development has contributed significantly to food security and reduction in poverty in Pakistan. Due to rapid population growth there has been a dramatic increase in the intensity of groundwater exploitation leading to declining water tables and deteriorating groundwater quality. In such prevailing conditions, the hydrogeological appraisal of escalating groundwater exploitation has become of paramount importance. Keeping this in view, a surface water–groundwater quantity and quality model was developed to assess future groundwater trends in the Rechna Doab (RD), a sub-catchment of the Indus River Basin. Scenario analysis shows that if dry conditions persist, there will be an overall decline in groundwater levels of around 10 m for the whole of RD during the next 25 years. The lower parts of RD with limited surface water supplies will undergo the highest decline in groundwater levels (10 to 20 m), which will make groundwater pumping very expensive for farmers. There is a high risk of groundwater salinization due to vertical upconing and lateral movement of highly saline groundwater into the fresh shallow aquifers in the upper parts of RD. If groundwater pumping is allowed to increase at the current rate, there will be an overall decline in groundwater salinity for the lower and middle parts of RD because of enhanced river leakage.     

Nutrient Mass Flux between Florida Bay and the Florida Keys National Marine Sanctuary

There is a net discharge of water and nutrients through Long Key Channel from Florida Bay to the Florida Keys National Marine Sanctuary (FKNMS). There has been speculation that this water and its constituents may be contributing to the loss of coral cover on the Florida Keys Reef tract over the past few decades, as well as speculation that changes in freshwater flow in the upstream Everglades ecosystem associated with the Comprehensive Everglades Restoration Plan may exacerbate this phenomenon. The results of this study indicate that although there is a net export of approximately 3,850 (±404) ton N year^?1 and 63 (±7) ton P year^?1, the concentrations of these nutrients flowing out of Florida Bay are the same as those flowing in. This implies that no significant nutrient enrichment is occurring in the waters of the FKNMS in the vicinity of Long Key Channel. Because of the effect of restricted southwestward water flow through Florida Bay by shallow banks and small islands, the volume of relatively high-nutrient water from central and eastern portions of the bay exiting through the channel is small compared to the average tidal exchange. Nutrient loading of relatively enriched bay waters is mediated by tidal exchange and mixing with more ambient concentrations of the western Florida Bay and Hawk Channel. System-wide budgets indicate that the contribution of Florida Bay waters to the inorganic nitrogen pool of the Keys coral reef is small relative to offshore inputs.     

A combined mathematical geophysical model for prediction of pyrite oxidation and pollutant leaching associated with a coal washing waste dump

The waste produced by coal washing process produces many environmental problems. In this study, the pollution problems associated with the waste produced by Alborz Sharghi Coal Washing Plant was investigated by mathematical modeling. The study area is located at 11 km. to Razmjah coal region and 45 km. to Tehran-Mashhad road in the north part od Iran. To achieve the goal, a few samples were taken from different depths at three points on the waste dump in order to investigate pyrite oxidation and pollution generation. The samples were then analysed, using an AA-670 Shimadzu atomic absorption to determine the fraction of pyrite remained within the waste particles. A numerical finite volume model using Phoenics package has been developed to simulate pyrite oxidation and pollution generation from the Alborz Sharghi coal washing waste dump. The pyrite oxidation reaction is described by the shrinking-core model. Gaseous diffusion is the main mechanism for the transport of oxygen through the waste. The results of numerical modelling were compared with the field observations and close agreement was achieved. A simple mathematical model incorporating advection and hydrodynamic dispersion processes was also presented in order to verify the results of geophysical time-laps method showing transportation of the pollutants through the downstream of the waste dump. Both mathematical model and geophysical time-laps method are agreed in the identification of pollutant transport emanated from the waste dump. The results of such investigations can be used for designing an effective environmental management program.     

Automated Kerogen Classification in Microscope Images of Dispersed Kerogen Preparation

We develop the classification part of a system that analyses transmitted light microscope images of dispersed kerogen preparation. The system automatically extracts kerogen pieces from the image and labels each piece as either inertinite or vitrinite. The image pre-processing analysis consists of background removal, identification of kerogen material, object segmentation, object extraction (individual images of pieces of kerogen) and feature calculation for each object. An expert palynologist was asked to label the objects into categories inertinite and vitrinite, which provided the ground truth for the classification experiment. Ten state-of-the-art classifiers and classifier ensembles were compared: Naïve Bayes, decision tree, nearest neighbour, the logistic classifier, multilayered perceptron (MLP), support vector machines (SVM), AdaBoost, Bagging, LogitBoost and Random Forest. The logistic classifier was singled out as the most accurate classifier, with an accuracy greater than 90. Using a 10 times 10-fold cross-validation provided within the Weka software, we found that the logistic classifier was significantly better than five classifiers (p<0.05) and indistinguishable from the other four classifiers. The initial set of 32 features was subsequently reduced to 6 features without compromising the classification accuracy. A further evaluation of the system alerted us to the possible sensitivity of the classification to the ground truth that might vary from one human expert to another. The analysis also revealed that the logistic classifier made most of the correct classifications with a high certainty.     

Influence of initial plastic anisotropy of overconsolidated clays on ground behaviour during tunneling

In order to evaluate the influence of the initial plastic anisotropy in the excavation of a tunnel, a “bubble” bounding surface model for structured soils, formulated by Kavvadas and Belokas (Proc. 10th IACMAG Conf., 2001), was implemented in the explicit finite difference code FLAC. Two different initial stress (K ₀) conditions were considered. The size and shape of the initial bounding surfaces were specified to be consistent with the initial stress field. The distorted and rotated shape of the bounding surface, supported by experimental results, defines the anisotropy of shear strength, which is shown to have a significant influence on the displacements. There is also considerable sensitivity of the soil model to the initial stress field.     

Microlites and "nanolites" in rhyolitic glass: microstructural and chemical characterization

 We used transmission electron microscopy (TEM) and scanning electron microscopy (SEM) to study magmatic crystals in the Ben Lomond rhyolite lava dome, Taupo Volcanic Center, New Zealand. Using TEM and SEM to investigate the size distributions of these crystals, we identified three size populations: microphenocrysts (>1.2 μm wide), microlites (>0.6 μm wide), and smaller crystals (<0.6 μm wide) which we term "nanolites". The predominant mineral phases of the microlites and nanolites are augites, pigeonites, and hypersthenes. The compositions and microstructures within these pyroxenes indicate disequilibrium crystallization at approximately 850–900  °C and undercoolings as high as 300  °C from equilibrium crystallization temperatures. Complex microstructures resulting from subsolidus reactions in augite and pigeonite are consistent with moderate cooling rates within the upper obsidian layer of the Ben Lomond rhyolite dome. This study demonstrates the existence of sub-micron magmatic crystals in a rhyolite and illustrates the potential of TEM to provide unique information about the crystallization and cooling histories of glassy volcanic rocks. Received: May 8, 1995 / Accepted: November 27, 1995     

LREE distribution in perovskite,apatite and titanite from South West Ugandan xenoliths and kamafugite lavas

Summary In-situ microprobe LREE analyses of perovskite and titanite (La, Ce, Nd), and apatite (La, Ce), from SW Ugandan clinopyroxenite xenoliths and kamafugite lavas indicate that LREE distribution in these minerals is determined by a number of factors related to their different parageneses: In particular LREE content is affected by whether the LREE-bearing minerals have crystallised from metasomatic carbonate or from silicate (i.e. metasomatic or magmatic) melts in the mantle. In this situation LREE partition favours carbonate over silicate melts. Distribution of LREE in perovskite and apatite crystallised from magmatic mantle melts or mantle-derived lavas is chiefly determined by preference of LREE for perovskite > apatite > titanite. LREE zoning in perovskite is influenced by changes in melt structure: increasing melt polymerisation enhancing mineral_LREE/melt_LREE partition into perovskite rims in magmatic xenoliths; decreasing melt polymerisation depleting LREE in lava perovskite rims. This zoning is reinforced by perovskite competition with apatite for LREE: perovskite (cores/rims) co-crystallising with apatite is reduced in LREE. There are 37 instances of perovskitewith Ce below detection while La and Nd levels are normal. These occur in both xenoliths and lavas; in grain zones or whole grains. Likewise Ce alone of the LREE is below detection in six out of ten titanite analyses. These observations are interpreted as evidence for increased f_{O 2}, Ce⁴ + being excluded from these mineral structures. Recognition of these various processes can elucidate the interpretation of bulk rock and bulk mineral LREE signatures in kamafugite volcanism.

---

LREE Verteilung in Perovskit, Apatit und Titanit aus Xenolithen und kamafugitischen Laven Südwest-Ugandas

Zusammenfassung In-situ LREE Analysen von Perovskit und Titanit (La, Ce, Nd) und Apatit (La, Ce) aus Klinopyroxenit-Xenolithen und kamafugitischen Laven Südwest-Ugandas zeigen, daß die LREE Verteilung in diesen Mineralen durch eine Vielzahl von Faktoren, die mit Unterschieden in den Paragenesen zusammenhängen, bestimmt wird: Der LREE-Gehalt wird im besonderen davon bestimmt, ob die LREE-führenden Minerale aus metasomatischen Karbonat- oder aus (metasomatischen oder magmatischen) Silikatschmelzen im Mantel auskristallisierten. Dabei erfolgt die LREE Fraktionierung zu Gunsten der Karbonatschmelzen. Die LREE-Verteilung von Perovskit und Apatit, die aus magmatischen Mantelschmelzen oder -laven kristallisierten, wird vorrangig durch den bevorzugten Einbau der LREE in Perovskit > Apatit > Titanit kontrolliert. Der LREE Zonarbau von Perovskit wird durch die Änderungen der Schmelzstruktur beinflußt: Verstärkte Schmelzpolymerisation führt zu verstärkter Mineral_LFEE/Schmelze_LREE Fraktionierung in den Perovskiträndern magmatischer Xenolithe, eine Abnahme der Schmelzpolymerisation hingegen resultiert in einer Abreicherung der LREE in den Perovskiträndern. Diese Art der Zonierung wird durch den Wettbewerb von Perovskit mit Apatit um die LREE verstärkt. Perovskit (Kerne/Ränder), der mit Apatit gemeinsam auskristallisierte, ist ärmer an LREE. 37 Fälle, in denenCe nicht nachweisbar war, La und Nd aber in normaler Konzentration auftreten, wurden sowohl in den Xenolithen als auch in den Laven gefunden; und zwar entweder in Kornbereichen oder in ganzen Körnern. Vergleichsweise liegt Ce nur in sechs von zehn Titanitproben unterhalb der Nachweisgrenze. Diese Beobachtungen werden als Hinweise auf erhöhte SauerstoffFugazitäten, bei denen Ce^4– aus der Mineralstruktur ausgeschlossen wird, angesehen.Ein Verständnis dieser verschiedenen Prozesse kann zur besseren Interpretation von LREE Gesamtgesteins- und Gesamtmineral-Signaturen in Kamafugiten beitragen.

---

---

With 3 Figures