Geochemical evidence for mixing of magmatic fluids with seawater, Nisyros hydrothermal system, Greece

The chemical and isotopic compositions (deltaD(H2O), delta(18)O(H2O), delta(18)O(CO2), delta(13)C(CO2), delta(34)S, and He/N-2 and He/Ar ratios) of fumarolic gases from Nisyros, Greece, indicate that both arc-type magmatic water and local seawater feed the hydrothermal system. Isotopic composition of the deep fluid is estimated to be +4.9+/-0.5parts per thousand for delta(18)O and -11+/-5parts per thousand for deltaD corresponding to a magmatic water fraction of 0.7. Interpretation of the stable water isotopes was based on liquid-vapor separation conditions obtained through gas geothermometry. The H-2-Ar, H-2-N-2, and H-2-H2O geothermometers suggest reservoir temperatures of 345+/-15 degreesC, in agreement with temperatures measured in deep geothermal wells, whereas a vapor/liquid separation temperature of 260+/-30 degreesC is indicated by gas equilibria in the H2O-H-2-CO2-CO-CH4 system. The largest magmatic inputs seem to occur below the Stephanos-Polybotes Micros crater, whereas the marginal fumarolic areas of Phlegeton-Polybotes Megalos craters receive a smaller contribution of magmatic gases.

Detection and characterization of hydraulically active fractures in a carbonate aquifer: Results from self-potential, temperature and fluid conductivity logging along a 260-m-deep borehole in the Combioula hydrothermal system in the southwestern Swiss Alps

A geophysical and geochemical study has been conducted in a fractured carbonate aquifer located at Combioula in the southwestern Swiss Alps with the objective to detect and characterize hydraulically active fractures along a 260-m-deep borehole. Hydrochemical analyses, borehole diameter, temperature and fluid electrical conductivity logging data were integrated in order to relate electrokinetic self-potential signals to groundwater flow inside the fracture network. The results show a generally good, albeit locally variable correlation of variations of the self-potential signals with variations in temperature, fluid electrical conductivity and borehole diameter. Together with the hydrochemical evidence, which was found to be critical for the interpretation of the self-potential data, these measurements not only made it possible to detect the hydraulically active fractures but also to characterize them as zones of fluid gain or fluid loss. The results complement the available information from the corresponding litholog and illustrate the potential of electrokinetic self-potential signals in conjunction with temperature, fluid electrical conductivity and hydrochemical analyses for the characterization of fractured aquifers, and thus may offer a perspective for an effective quantitative characterization of this increasingly important class of aquifers and geothermal reservoirs.

Optimization of Large-Scale Hydrothermal System Operation

This paper presents the development of a mathematical model to optimize the management and operation of the Brazilian hydrothermal system. The system consists of a large set of individual hydropower plants and a set of aggregated thermal plants. The energy generated in the system is interconnected by a transmission network so it can be transmitted to centers of consumption throughout the country. The optimization model offered is capable of handling different types of constraints, such as interbasin water transfers, water supply for various purposes, and environmental requirements. Its overall objective is to produce energy to meet the country's demand at a minimum cost. Called HIDROTERM, the model integrates a database with basic hydrological and technical information to run the optimization model, and provides an interface to manage the input and output data. The optimization model uses the General Algebraic Modeling System (GAMS) package and can invoke different linear as well as nonlinear programming solvers. The optimization model was applied to the Brazilian hydrothermal system, one of the largest in the world. The system is divided into four subsystems with 127 active hydropower plants. Preliminary results under different scenarios of inflow, demand, and installed capacity demonstrate the efficiency and utility of the model. From this and other case studies in Brazil, the results indicate that the methodology developed is suitable to different applications, such as planning operation, capacity expansion, and operational rule studies, and trade-off analysis among multiple water users. DOI: 10.1061/(ASCE)WR.1943-5452.0000149. (C) 2012 American Society of Civil Engineers.

Modeling hydrothermal system: deriving observables and hydrothermal instability in volcanic and non-volcanic setting

Hydrothermal fluids are a fundamental resource for understanding and monitoring volcanic and non-volcanic systems. This thesis is focused on the study of hydrothermal system through numerical modeling with the geothermal simulator TOUGH2. Several simulations are presented, and geophysical and geochemical observables, arising from fluids circulation, are analyzed in detail throughout the thesis. In a volcanic setting, fluids feeding fumaroles and hot spring may play a key role in the hazard evaluation. The evolution of the fluids circulation is caused by a strong interaction between magmatic and hydrothermal systems. A simultaneous analysis of different geophysical and geochemical observables is a sound approach for interpreting monitored data and to infer a consistent conceptual model. Analyzed observables are ground displacement, gravity changes, electrical conductivity, amount, composition and temperature of the emitted gases at surface, and extent of degassing area. Results highlight the different temporal response of the considered observables, as well as the different radial pattern of variation. However, magnitude, temporal response and radial pattern of these signals depend not only on the evolution of fluid circulation, but a main role is played by the considered rock properties. Numerical simulations highlight differences that arise from the assumption of different permeabilities, for both homogeneous and heterogeneous systems. Rock properties affect hydrothermal fluid circulation, controlling both the range of variation and the temporal evolution of the observable signals. Low temperature fumaroles and low discharge rate may be affected by atmospheric conditions. Detailed parametric simulations were performed, aimed to understand the effects of system properties, such as permeability and gas reservoir overpressure, on diffuse degassing when air temperature and barometric pressure changes are applied to the ground surface. Hydrothermal circulation, however, is not only a characteristic of volcanic system. Hot fluids may be involved in several mankind problems, such as studies on geothermal engineering, nuclear waste propagation in porous medium, and Geological Carbon Sequestration (GCS). The current concept for large-scale GCS is the direct injection of supercritical carbon dioxide into deep geological formations which typically contain brine. Upward displacement of such brine from deep reservoirs driven by pressure increases resulting from carbon dioxide injection may occur through abandoned wells, permeable faults or permeable channels. Brine intrusion into aquifers may degrade groundwater resources. Numerical results show that pressure rise drives dense water up to the conduits, and does not necessarily result in continuous flow. Rather, overpressure leads to new hydrostatic equilibrium if fluids are initially density stratified. If warm and salty fluid does not cool passing through the conduit, an oscillatory solution is then possible. Parameter studies delineate steady-state (static) and oscillatory solutions.

Bajo de la Alumbrera copper-gold deposit: Stable isotope evidence for a porphyry-related hydrothermal system dominated by magmatic aqueous fluids

Alteration zones at the gold-rich Bajo de la Alumbrera porphyry copper deposit in northwestern Argentina are centered on several porphyritic intrusions. They are zoned from a central copper-iron sulfide and gold-mineralized potassic (biotite-K-feldspar +/- quartz) core outward to propylitic (chlorite-illite-epidote-calcite) assemblages. A mineralized intermediate argillic alteration assemblage (chlorite-illite +/- pyrite) has overprinted the potassic alteration zone across the top and sides of the deposit and is itself zoned outward into phyllic (quartzinuscovite-illite +/- pyrite) alteration. This study contributes new data to previously reported delta(18)O and delta D compositions of fluids responsible for the alteration at Bajo de la Alumbrera, and the data are used to infer likely ore-forming processes. Measured and calculated delta(18)O and delta D values of fluids (+8.3 to +10.2 and -33 to -81 parts per thousand, respectively) confirm a primary magmatic origin for the earliest potassic alteration phase. Lower temperature potassic alteration formed from magmatic fluids with lower delta D values (down to -123 parts per thousand). These depleted compositions are distinct from meteoric water and consistent with degassing and volatile exsolution of magmatic fluids derived from an underlying magma. Variability in the calculated composition of fluid associated with potassic alteration is explained in terms of phase separation (or boiling). if copper-iron sulfide deposition occurred during cooling (as proposed elsewhere), this cooling was largely a result of phase separation. Magmatic water was directly involved in the formation of overprinting intermediate argillic alteration assemblages at Bajo de la Alumbrera. Calculated delta(18)O and delta D values of fluids associated with this alteration range from +4.8 to +8.1 and -31 to -71 per mil, respectively Compositions determined for fluids associated with phyllic alteration (-0.8 to +10.2 and -31 to -119 parts per thousand) overlap with the values determined for the intermediate argillic alteration. We infer that phyllic alteration assemblages developed during two stages; the first was a high-temperature (400 degrees-300 degrees C) stage with D-depleted water (delta D = -66 to -119 parts per thousand). This compositional range may have resulted from magma degassing and/or the injection of new magmatic water into a compositionally evolved hydrothermal system. The isotopic variations also can be explained by increased fluid-rock interaction. The second stage of phyllic alteration occurred at a lower temperature (similar to 200 degrees C), and variations in the modeled isotopic compositions imply mixing of magmatic and meteoric waters. Ore deposition that occurred late in the evolution of the hydrothermal system was probably associated with further cooling of the magmatic fluid, in part caused by fluid-rock interaction and phase separation. Changing pH and/or oxygen fuoracity may have caused additional ore deposition. The ingress of meteoric water appears to postdate the bulk of mineralization and occurred as the system at Bajo de la Alumbrera waned.

A model of the hydrothermal system at Casa Diablo in Long Valley, California, based on resistivity profiles and soil mercury analyses

A description and model of the near-surface hydrothermal system at Casa Diablo, with its implications for the larger-scale hydrothermal system of Long Valley, California, is presented. The data include resistivity profiles with penetrations to three different depth ranges, and analyses of inorganic mercury concentrations in 144 soil samples taken over a 1.3 by 1.7 km area. Analyses of the data together with the mapping of active surface hydrothermal features (fumaroles, mudpots, etc.), has revealed that the relationship between the hydrothermal system, surface hydrothermal activity, and mercury anomalies is strongly controlled by faults and topography. There are, however, more subtle factors responsible for the location of many active and anomalous zones such as fractures, zones of high permeability, and interactions between hydrothermal and cooler groundwater. In addition, the near-surface location of the upwelling from the deep hydrothermal reservoir, which supplies the geothermal power plants at Casa Diablo and the numerous hot pools in the caldera with hydrothermal water, has been detected. The data indicate that after upwelling the hydrothermal water flows eastward at shallow depth for at least 2 km and probably continues another 10 km to the east, all the way to Lake Crowley.

Finite element analysis of heat transfer and mineralization in layered hydrothermal systems with upward throughflow

We use the finite element method to solve the coupled problem between convective pore-fluid flow, heat transfer and mineralization in layered hydrothermal systems with upward throughflow. In particular, we present the improved rock alteration index (IRAI) concept for predicting the most probable precipitation and dissolution regions of gold (Au) minerals in the systems. To validate the numerical method used in the computation, analytical solutions to a benchmark problem have been derived. After the numerical method is validated, it is used to investigate the pattern of pore-fluid Aom, the distribution of temperature and the mineralization pattern of gold minerals in a layered hydrothermal system with upward throughflow. The related numerical results have demonstrated that the present concept of IRAI is useful and applicable for predicting the most probable precipitation and dissolution regions of gold (Au) minerals in hydrothermal systems. (C) 2000 Elsevier Science S.A. All rights reserved.

Finite element modelling of rock alteration and metamorphic process in hydrothermal systems

We use the finite element method to simulate the rock alteration and metamorphic process in hydrothermal systems. In particular, we consider the fluid-rock interaction problems in pore-fluid saturated porous rocks. Since the fluid rock interaction takes place at the contact interface between the pore-fluid and solid minerals, it is governed by the chemical reaction which usually takes place very slowly at this contact interface, from the geochemical point of view. Due to the relative slowness of the rate of the chemical reaction to the velocity of the pore-fluid flow in the hydrothermal system to be considered, there exists a retardation zone, in which the conventional static theory in geochemistry does not hold true. Since this issue is often overlooked by some purely numerical modellers, it is emphasized in this paper. The related results from a typical rock alteration and metamorphic problem in a hydrothermal system have shown not only the detailed rock alteration and metamorphic process, but also the size of the retardation zone in the hydrothermal system. Copyright (C) 2001 John Wiley & Sons, Ltd.

Oligo-Miocene extensional tectonics and fluid flow across the Northern Snake Range detachment system, Nevada

The Northern Snake Range (Nevada) represents a spectacular example of a metamorphic core complex and exposes a complete section from the mylonitic footwall into the hanging wall of a fossil detachment system. Paired geochronological and stable isotopic data of mylonitic quartzite within the detachment footwall reveal that ductile deformation and infiltration of meteoric fluids occurred between 27 and 23 Ma. Ar-40/Ar-39 ages display complex recrystallization-cooling relationships but decrease systematically from 26.9 +/- 0.2 Ma at the top to 21.3 +/- 0.2 Ma at the bottom of footwall mylonite. Hydrogen isotope (delta D) values in white mica are very low (-150 to -145 %) within the top 80-90 m of detachment footwall, in contrast to values obtained from the deeper part of the section where values range from -77 to -64 %, suggesting that time-integrated interaction between rock and meteoric fluid was restricted to the uppermost part of the mylonitic footwall. Pervasive mica-water hydrogen isotope exchange is difficult to reconcile with models of Ar-40 loss during mylonitization solely by volume diffusion. Rather, we interpret the Ar-40/Ar-39 ages of white mica with low-delta D values to date syn-mylonitic hydrogen and argon isotope exchange, and we conclude that the hydrothermal system of the Northern Snake Range was active during late Oligocene (27-23 Ma) and has been exhumed by the combined effects of ductile strain, extensional detachment faulting, and erosion.

Bitumens in the late Variscan hydrothermal vein-type uranium deposit of Pribram, Czech Republic: Sources, radiation-induced alteration, and relation to mineralization

The late Variscan (275-278 Ma) Pribram uranium deposit is one of the largest known accumulations of uraniferous bitumens in hydrothermal veins. The deposit extends along the northwestern boundary of the Central Bohemian pluton (345-335 Ma) with low-grade metamorphosed Late Proterozoic and unmetamorphosed Cambrian rocks. From a net uranium production of 41,742 metric tons (t), more than 6,000 t were extracted from bitumen-uraninite ores during 43 years of exploration and mining. Three morphological varieties of solid bitumen are recognized: globular, asphaltlike, and cokelike. While the globular bitumen is uranium free, the other two types are uraniferous. The amount of bitumen in ore veins gradually decreases toward the contact with the plutonic body and increases with depth. Two types of bitumen microtextures are recognized using high-resolution transmission electron microscopy: amorphous and microporous, the former being less common in uraniferous samples. A lower Raman peak area ratio (1,360/1,575 cm(-1)) in mineralized bitumens (0.9) compared with uranium-free samples (2.0) indicates a lower degree of microtextural organization in the latter The H/C and O/C atomic ratios in uranium-free bitumens (0.9-1.1 and 0.09, respectively) are higher than those in mineralized samples (H/C = 0.3-0.8, O/C = 0.03-0.09). The chloroform extractable matter yield is Very low in uranium-free bitumens (0.30-0.35% of the total organic carbon,TOC) and decreases with uranium content increase. The extracted solid uraniferous bitumen infrared spectra show depletion in aliphatic CH2 and CH3 groups compared to uranium-free samples. The concentration of oxygen-bearing functional groups relative to aromatic bonds in the IR spectra of uranium-free and mineralized bitumen, however, do not differ significantly. C-13 NMR confirmed than the aromaticity of a uraniferous sample is higher (F-ar = 0.61) than in the uranium-free bitumen (F-ar = 0.51). Pyrolysates from uraniferous and nonuraniferous bitumens do not differ significantly, being predominantly cresol, alkylphenols, alkylbenzenes, and alkylnaphthalenes. The liquid pyrolysate yield decreases significantly with increasing uranium content. The delta(13)C Values of bulk uranium-free bitumens and low-grade uraniferous, asphaltlike bitumens range from -43.6 to 52.3 per mil. High-grade, cokelike, uraniferous bitumens are more C-13 depleted (54.5 to -58.4 parts per thousand). In contrast to the very light isotopic ratios of the high-grade uraniferous cokelike bitumen bulk carbon, the individual n-alkanes and isoprenoids (pristane and phytane) extracted from the same sample are significantly C-13 enriched. The isotopic composition of the C13-24 n-alkanes extracted from the high-grade uraniferous sample (delta(13)C = -28.0 to 32.6 parts per thousand) are heavier compared with the same compounds in a uranium-free sample (delta(13)C = 31.9 to 33.8 parts per thousand). It is proposed that the bitumen source was the isotopically light (delta(13)C = 35.8 to 30.2 parts per thousand) organic matter of the Upper Proterozoic host rocks that were pyrolyzed during intrusion of the Central Bohemian pluton. The C-13- depleted pyrolysates were mobilized from the innermost part of the contact-metamorphic aureole, accumulated in structural traps in less thermally influenced parts of the sedimentary complex and were later extracted by hydrothermal fluids. Bitumens at the Pribram deposit are younger than the main part of the uranium mineralization and were formed through water-washing and radiation-induced polymerization of both the gaseous and liquid pyrolysates. Direct evidence for pyrolysate reduction of uranium in the hydrothermal system is difficult to obtain as the chemical composition of the original organic fluid phase was modified during water-washing and radiolytic alteration. However, indirect evidence-e.g., higher O/C atomic ratios in uranium-free bitumens (0.1) relative to the Upper Proterozoic source rocks (0.02-0.05), isotopically very light carbon in associated whewellite (delta(13)C = 31.7 to -28.4 parts per thousand), and the striking absence of bitumens in the pre-uranium, hematite stage of the mineralization-indicates that oxidation of organic fluids may have contributed to lowering of aO(2) and uraninite precipitation.

Carbon and oxygen isotope study of hydrothermal carbonates in the zinc-lead deposits of the San Vicente district, central Peru: A quantitative modeling on mixing processes and CO2 degassing

Carbon and oxygen isotope studies of the host and gangue carbonates of Mississippi Valley-type zinc-lead deposits in the San Vicente District hosted in the Upper Triassic to Lower Jurassic dolostones of the Pucara basin (central Peru) were used to constrain models of the ore formation. A mixing model between an incoming hot saline slightly acidic radiogenic (Pb, Sr) fluid and the native formation water explains the overall isotopic variation (delta(13)C = - 11.5 to + 2.5 parts per thousand relative to PDB and delta(18)O = + 18.0 to + 24.3 parts per thousand relative to SMOW) of the carbonate generations. The dolomites formed during the main ore stage show a narrower range (delta(13)C = - 0.1 to + 1.7 parts per thousand and delta(18)O = + 18.7 to + 23.4 parts per thousand) which is explained by exchange between the mineralizing fluids and the host carbonates combined with changes in temperature and pressure. This model of fluid-rock interaction explains the pervasive alteration of the host dolomite I and precipitation of sphalerite I. The open-space filling hydrothermal white sparry dolomite and the coexisting sphalerite II formed by prolonged fluid-host dolomite interaction and limited CO2 degassing. Late void-filling dolomite III (or calcite) and the associated sphalerite III formed as the consequence of CO2 degassing and concomitant pH increase of a slightly acidic ore fluid. Widespread brecciation is associated to CO2 outgassing. Consequently, pressure variability plays a major role in the ore precipitation during the late hydrothermal events in San Vicente. The presence of native sulfur associated with extremely carbon-light calcites replacing evaporitic sulfates (e.g., delta(13)C = - 11.5 parts per thousand), altered native organic matter and heavier hydrothermal bitumen (from - 27.0 to - 23.0 parts per thousand delta(13)C) points to thermochemical reduction of sulfate and/or thiosulfate. The delta(13)C- and delta(18)O-values of the altered host dolostone and hydrothermal carbonates, and the carbon isotope composition of the associated organic matter show a strong regional homogeneity. These results coupled with the strong mineralogical and petrographic similarities of the different MVT occurrences perhaps reflects the fact that the mineralizing processes were similar in the whole San Vicente belt, suggesting the existence of a common regional mineralizing hydrothermal system with interconnected plumbing.

Characterization of magmatic sulfur in the Aegean island arc by means of the δ34S values of fumarolic H2S, elemental S, and hydrothermal gypsum from Nisyros and Milos islands

A delta(34)S value of +6.3 +/- 1.5% was estimated for the rhyodacitic degassing magma present underneath the hydrothermal system of Nisyros, based on the S isotope ratios of H2S in fumarolic vapors. This value was estimated by modeling the irreversible water-rock mass transfers occurring during the generation of the hydrothermal liquid which separates these fumarolic vapors. The S isotope ratio of the rhyodacitic degassing magma of Nisyros is consistent with fractional crystallization of a parent basaltic magma with an initial delta(34)S value of +4% (+/-at least 1.5%). This positive value could be explained by mantle contamination due to by either transference of fluids derived from subducted materials or involvement of altered oceanic crust, whereas contribution of biogenic sulfides from sediments seems to be negligible or nil. This conclusion agrees with the lack of N-2 and CO2 from thermal decomposition of organic matter contained in subducted sediments, which is a characteristic of the whole Aegean arc system. Since hydrothermal S at Milos and Santorini has isotope ratios similar to those determined at Nisyros, it seems likely that common controlling processes are active throughout the Aegean island arc. (C) 2002 Elsevier, Science B.V. All rights reserved.

Crystal growth and photoluminescence of europium-doped strontium titanate prepared by a microwave hydrothermal method

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)