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.

     

In situ estimation of the effective chemical diffusion coefficient of a rock matrix in a fractured aquifer

An in situ method of estimating the effective diffusion coefficient for a chemical constituent that diffuses into the primary porosity of a rock is developed by abruptly changing the concentration of the dissolved constituent in a borehole in contact with the rock matrix and monitoring the time-varying concentration. The experiment was conducted in a borehole completed in mudstone on the campus of the University of the Free State in Bloemfontein, South Africa. Numerous tracer tests were conducted at this site, which left a residual concentration of sodium chloride in boreholes that diffused into the rock matrix over a period of years. Fresh water was introduced into a borehole in contact with the mudstone, and the time-varying increase of chloride was observed by monitoring the electrical conductivity (EC) at various depths in the borehole. Estimates of the effective diffusion coefficient were obtained by interpreting measurements of EC over 34 d. The effective diffusion coefficient at a depth of 36 m was approximately 7.8×10^?6 m²/d, but was sensitive to the assumed matrix porosity. The formation factor and mass flux for the mudstone were also estimated from the experiment.     

Estimation of hydraulic permeability considering the micro morphology of rocks of the borehole YAXCOPOIL-1 (Impact crater Chicxulub,Mexico)

Internal surface, formation factor, Nuclear Magnetic Resonance (NMR)-T2 relaxation times and pore radius distributions were measured on representative core samples for the estimation of hydraulic permeability. Permeability is estimated using various versions of the classic Kozeny–Carman-equation (K–C) and a further development of K–C, the fractal PaRiS-model, taking into account the internal surface. In addition to grain and pore size distribution, directly connected to permeability, internal surface reflects the internal structure (“micro morphology”). Lithologies could be grouped with respect to differences in internal surface. Most melt rich impact breccia lithologies exhibit large internal surfaces, while Tertiary post-impact sediments and Cretaceous lithologies in displaced megablocks display smaller internal surfaces. Investigations with scanning electron microscopy confirm the correlation between internal surface and micro morphology. In addition to different versions of K–C, estimations by means of NMR, pore radius distributions and some gas permeability measurements serve for cross-checking and calibration. In general, the different estimations from the independent methods and the measurements are in satisfactory accordance. For Tertiary limestones and Suevites bulk with very high porosities (up to 35%) permeabilites between 10⁻¹⁴ and 10⁻¹⁶ m² are found, whereas in lower Suevite, Cretaceous anhydrites and dolomites, bulk permeabilites are between 10⁻¹⁵ and 10⁻²³ m².     

Hydrogeochemistry and isotope analyses used to determine groundwater recharge and flow in low-gradient catchments of the province of Buenos Aires,Argentina

The region known as Pampa Plain, in Argentina, is a vast area characterized by slopes of less than 0.05%. The surface sediment is silty sand, mainly Aeolian, referred to as Pampean Loess, which is a phreatic aquifer unit of utmost importance for the water supply of the region. On account of the slight gradient, hydrogeological analyses using only hydraulic measurements are difficult to perform, often leading to confusing results. Thus, the study presented relies on hydrochemical modeling and isotopic determinations as well. The study area comprises three catchments in the inter-mountainous area corresponding to El Moro, Tamangueyú and Seco creeks, and covering 2,570 km² in the province of Buenos Aires. Measurements of piezometric levels and samples of groundwater, surface water and rainwater were carried out between January and March 2005. Major ion results were analyzed by means of hydrochemical graphs and hydrogeochemical modelling using NETPATH. The resulting data show that it is possible to identify local changes in recharge, flow direction and stream/groundwater relationship by using hydrochemical and isotopic information, which may become a useful and more precise tool for the study of particularly flat landscapes.     

Estimates of heat flow and heat production and a thermal model of the São Francisco craton

An updated analysis of geothermal data from the highland area of eastern Brazil has been carried out and the characteristics of regional variations in geothermal gradients and heat flow examined. The database employed includes results of geothermal measurements at 45 localities. The results indicate that the Salvador craton and the adjacent metamorphic fold belts northeastern parts of the study area are characterized by geothermal gradients in the range of 6–17°C/km. The estimated heat flow values fall in the range of 28–53 mW/m², with low values in the cratonic area relative to the fold belts. On the other hand, the São Francisco craton and the intracratonic São Francisco sedimentary basin in the southwestern parts are characterized by relatively higher gradient values, in the range of 14–42°C/km, with the corresponding heat flow values falling in the range of 36–89 mW/m². Maps of regional variations indicate that high heat flow anomaly in the São Francisco craton is limited to areas of sedimentary cover, to the west of the Espinhaço mountain belt. Crustal thermal models have been developed to examine the implications of the observed intracratonic variations in heat flow. The thermal models take into consideration variation of thermal conductivity with temperature as well as change of radiogenic heat generation with depth. Vertical distributions of seismic velocities were used in obtaining estimates of radiogenic heat production in crustal layers. Crustal temperatures are calculated based on a procedure that makes simultaneous use of the Kirchoff and Generalized Integral Transforms, providing thereby analytical solutions in 2D and 3D geometry. The results point to temperature variations of up to 300°C at the Moho depth, between the northern Salvador and southern São Francisco cratons. There are indications that differences in rheological properties, related to thermal field, are responsible for the contrasting styles of deformation patterns in the adjacent metamorphic fold belts.     

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.     

High-resolution seismic analysis of the coastal Mecklenburg Bay (North German Basin): the pre-Alpine evolution

The pre-Alpine structural and geological evolution in the northern part of the North German Basin have been revealed on the basis of a very dense reflection seismic profile grid. The study area is situated in the coastal Mecklenburg Bay (Germany), part of the southwestern Baltic Sea. From the central part of the North German Basin to the northern basin margin in the Grimmen High area a series of high-resolution maps show the evolution from the base Zechstein to the Lower Jurassic. We present a map of basement faults affecting the pre-Zechstein. The pre-Alpine structural evolution of the region has been determined from digital mapping of post-Permian key horizons traced on the processed seismic time sections. The geological evolution of the North German Basin can be separated into four distinct periods in the Rerik study area. During Late Permian and Early Triassic evaporites and clastics were deposited. Salt movement was initiated after the deposition of the Middle Triassic Muschelkalk. Salt pillows, which were previously unmapped in the study area, are responsible for the creation of smaller subsidence centers and angular unconformities in the Late Triassic Keuper, especially in the vicinity of the fault-bounded Grimmen High. In this area, partly Lower Jurassic sediments overlie the Keuper unconformably. The change from extension to compression in the regional stress field remobilized the salt, leading to a major unconformity marked at the base of the Late Cretaceous.     

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.     

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.