Hydrogeology of the Krafla geothermal system,northeast Iceland

The Krafla geothermal system is located in Iceland's northeastern neovolcanic zone, within the Krafla central volcanic complex. Geothermal fluids are superheated steam closest to the magma heat source, two‐phase at higher depths, and sub‐boiling at the shallowest depths. Hydrogen isotope ratios of geothermal fluids range from ?87‰, equivalent to local meteoric water, to ?94‰. These fluids are enriched in ¹⁸O relative to the global meteoric line by +0.5–3.2‰. Calculated vapor fractions of the fluids are 0.0–0.5 wt% (~0–16% by volume) in the northwestern portion of the geothermal system and increase towards the southeast, up to 5.4 wt% (~57% by volume). Hydrothermal epidote sampled from 900 to 2500 m depth has δD values from ?127 to ?108‰, and δ¹⁸O from ?13.0 to ?9.6‰. Fluids in equilibrium with epidote have isotope compositions similar to those calculated for the vapor phase of two‐phase aquifer fluids. We interpret the large range in δD_EPIDOTE and δ¹⁸O_EPIDOTE across the system and within individual wells (up to 7‰ and 3.3‰, respectively) to result from variable mixing of shallow sub‐boiling groundwater with condensates of vapor rising from a deeper two‐phase reservoir. The data suggest that meteoric waters derived from a single source in the northwest are separated into the shallow sub‐boiling reservoir, and deeper two‐phase reservoir. Interaction between these reservoirs occurs by channelized vertical flow of vapor along fractures, and input of magmatic volatiles further alters fluid chemistry in some wells. Isotopic compositions of hydrothermal epidote reflect local equilibrium with fluids formed by mixtures of shallow water, deep vapor condensates, and magmatic volatiles, whose ionic strength is subsequently derived from dissolution of basalt host rock. This study illustrates the benefits of combining phase segregation effects in two‐phase systems during analysis of wellhead fluid data with stable isotope values of hydrous alteration minerals when evaluating the complex hydrogeology of volcano‐hosted geothermal systems.     

The Ixtacamaxtitlán kaolinite deposit and sinter (Puebla State, Mexico): a magmatic-hydrothermal system telescoped by a shallow paleoaquifer

The Ixtacamaxtitlán hydrothermal deposit is made up of a succession, from bottom to top, of: (1) a porphyritic subvolcanic body, crosscut by quartz veins, and a stockwork with subordinate sulfides (pyrite and chalcopyrite), showing propylitic alteration haloes overprinting a previous potassic alteration event (biotitization); (2) an overlying, kaolinized lithic‐rich rhyolitic tuff; and (3) a layered opal deposit with preserved sedimentary structures. This vertical zonation, coupled with the distribution of the alteration assemblages, lead us to the interpretation of the whole as a porphyry‐type deposit grading upwards to a barren, steam‐heated, acid‐leached, kaolinite blanket with a partially preserved silica sinter on top. Both the fluid inclusion study carried out on the veins and stockwork, and the stable isotopic analyses of the kaolinized bodies, suggest the presence of two major hydrothermal events. The older event is characterized by the occurrence of hot hypersaline fluids (up to 320°C and 36 wt% NaCl equivalents), likely of magmatic origin, closely associated with the emplacement of the underlying early Miocene porphyry‐type deposit. The later event is characterized by the presence of cooler and dilute fluids (up to 140°C and 4 wt% NaCl equivalents) and by advanced argillic alteration close to the paleosurface. The calculated isotopic composition of water in equilibrium with the kaolinitic sequence plots close to and underneath the meteoric water line, partially overlapping the Los Humeros present‐day geothermal fluids. This evidence coupled with the petrographic observations suggests that steam‐heated phreatic waters altered the lithic‐rich rhyolitic tuff. This would have occurred when acid vapors, exsolved from deeper hydrothermal fluids by boiling, reached the local paleowater table and condensed, after a sector collapse that changed the system from lithostatic to hydrostatic conditions.     

Hydrogeochemical dynamics affecting steam‐heated pools at El Chichón Crater (Chiapas – Mexico)

El Chichón is an active volcano located in the north‐western Chiapas, southern Mexico. The crater hosts a lake, a spring, named Soap Pool, emerging from the underlying volcanic aquifer and several mud pools/hot springs on the internal flanks of the crater which strongly interact with the current fumarolic system (steam‐heated pools). Some of these pools, the crater lake and a cold spring emerging from the 1982 pumice deposits, have been sampled and analysed. Water–volcanic gas interactions determine the heating (43–99°C) and acidification (pH 2–4) of the springs, mainly by H₂S oxidation. Significantly, in the study area, a significant NH₃ partial pressure has been also detected. Such a geochemically aggressive environment enhances alteration of the rock in situ and strongly increases the mineralization of the waters (and therefore their electrical conductivity). Two different mineralization systems were detected for the crater waters: the soap pool‐lake (Na⁺/Cl^? = 0.4, Na/Mg>10) and the crater mud pools (Na⁺/Cl^? > 10, Na/Mg < 4). A deep boiling, Na⁺‐K⁺‐Cl^?‐rich water reservoir generally influences the Soap Pool‐lake, while the mud pool is mainly dominated by water‐gas–rock interactions. In the latter case, conductivity of sampled water is directly proportional to the presence of reactive gases in solution. Therefore, chemical evolution proceeds through neutralization due to both rock alteration and bacterial oxidation of ammonium to nitrate. The chemical compositions show that El Chichón aqueous fluids, within the crater, interact with gases fed by a geothermal reservoir, without clear additions of deep magmatic fluids. This new geochemical dataset, together with previously published data, can be used as a base line with which to follow‐up the activity of this deadly volcano.     

Geothermal fluids circulation at Caldas do Moledo area, Northern Portugal: geochemical and isotopic signatures

A hydrogeological conceptual model of the Caldas do Moledo geothermal site is proposed that shows mixing between geothermal waters and local shallow groundwaters. Stable isotope values of Caldas do Moledo geothermal waters indicate recharge areas located at relatively high altitudes (850–1250 m a.s.l.). The NW–SE Vigo–Régua shear zone plays an important role in fluid recharge and circulation towards the NNE–SSW Régua–Verin fault system, forming a path for ascent of geothermal fluids. The apparent ¹⁴C age of geothermal fluids (15.66 ± 2.86 ka BP) was estimated in the total dissolved inorganic carbon (TDIC). Geothermometer calculations indicate that, assuming a conductive temperature gradient of 32°C per kilometer for northern Portugal, the maximum depth of circulation is roughly 1.8 ± 0.4 km. The K, Ca and SO₄ concentrations found in some Caldas do Moledo geothermal spring waters show mixing between deep geothermal and shallow groundwater systems. Local shallow groundwaters showing the highest SO₄ concentrations were found at low elevation areas, originating from fertilisers and pesticides applied to the Port wine vineyards in the Douro River valley. Geothermal waters from boreholes AC1 and AC2 do not show evidences of direct pollution from the spreading of such agrochemicals.     

Hydraulic and hydrochemical properties of deep sedimentary reservoirs of the Upper Rhine Graben,Europe

Hydraulic and hydrochemical data from several hundred wells mostly drilled by the oil and gas industry within the four deep carbonate and siliciclastic reservoirs of the Upper Rhine Graben area in France and Germany have been compiled, examined, validated and analysed with the aim to characterize fluids and reservoir properties. Due to enhanced temperatures in the subsurface of the Upper Rhine Graben, this study on hydraulic and hydrochemical properties has been motivated by an increasing interest in deep hydrogeothermal energy projects in the Rhine rift valley. The four examined geothermal reservoir formations are characterized by high hydraulic conductivity reflecting the active tectonic setting of the rift valley and its fractured and karstified reservoirs. The hydraulic conductivity decreases only marginally with depth in each of the reservoirs, because the Upper Rhine Graben is a young tectonically active structure. The generally high hydraulic conductivity of the reservoir rocks permits cross‐formation advective flow of thermal water. Water composition data reflect the origin and hydrochemical evolution of deep water. Shallow water to 500 m depth is, in general, weakly mineralized. The chemical signature of the water is controlled by fluid–rock geochemical interactions. With increasing depth, the total of dissolved solids (TDS) increases. In all reservoirs, the fluids evolve to a NaCl‐dominated brine. The high salinity of the reservoirs is partly derived from dissolution of halite in evaporitic Triassic and Cenozoic formations, and partly from the fluids residing in the crystalline basement. Water of all four reservoirs is saturated with respect to calcite and other minerals including quartz and barite.     

Real‐time quadrupole mass spectrometer analysis of gas in borehole fluid samples acquired using the U‐tube sampling methodology

Sampling of fluids in deep boreholes is challenging because of the necessity of minimizing external contamination and maintaining sample integrity during recovery. The U‐tube sampling methodology was developed to collect large volume, multiphase samples at in situ pressures. As a permanent or semi‐permanent installation, the U‐tube can be used for rapidly acquiring multiple samples or it may be installed for long‐term monitoring applications. The U‐tube was first deployed in Liberty County, TX to monitor crosswell CO₂ injection as part of the Frio CO₂ sequestration experiment. Analysis of gases (dissolved or separate phase) was performed in the field using a quadrupole mass spectrometer, which served as the basis for determining the arrival of the CO₂ plume. The presence of oxygen and argon in elevated concentrations, along with reduced methane concentration, indicates sample alteration caused by the introduction of surface fluids during borehole completion. Despite producing the well to eliminate non‐native fluids, measurements demonstrate that contamination persists until the immiscible CO₂ injection swept formation fluid into the observation wellbore.     

Origins of formation waters in the Llanos foreland basin of Colombia: geochemical variation and fluid flow history

Wells located in the Colombian Andean foreland often produce biodegraded hydrocarbons and relatively fresh water (total dissolved solids concentrations of 2000 mg l^?1 dominate across most of the basin). The ratio of water to hydrocarbon in the produced fluids is high and has often been interpreted by explorationists and simulation engineers as being due to massive meteoric water invasion. To challenge this hypothesis, previously published and existing data are compiled here, along with 68 new water samples collected from both the surface (rivers) and subsurface (production wells), which are analyzed for major elements, salinity, and stable isotope compositions. The data indicate a mixed origin for the water, which involves formation water, meteoric water, and a third source that we interpret as being related to diagenetic processes (such as clay dehydration). Analyses of the variation of the parameters from north to south and west to east allow us to define the mixtures as coming from waters of different origins and make it possible to determine the relative contribution of each source.     

Application of boron isotopes to the understanding of fluid–rock interactions in a hydrothermally stimulated oil reservoir in the Alberta Basin, Canada

Boron isotope ratios of reservoir minerals and fluids can be a useful geothermometer and monitor of fluid–rock interactions. In Cold Lake oil sands of northern Alberta, there is a large variation in δ¹¹B of the produced waters generated during steam injection and recovery of oil and water. The higher temperature waters (～ 200 °C) have isotopically light δ¹¹B values (+ 3‰ to + 14‰) and high B contents (～150 p.p.m.). It is inferred that the range of δ¹¹B values of the hydrothermal fluids results from reaction with the reservoir rock, and is a function of the temperature of the fluid–rock interaction. The distinct B geochemistry of the produced waters can be used to show that there is no detectable mixing of the oil recovery waters with the regional formation waters or shallow groundwater aquifers containing potable water. Examination of B isotope ratios of reservoir minerals, before and after steam injection, allows the evaluation of sources of B in the reservoir. The only significant phase containing B is pumice. It shows generally positive δ¹¹B values before steam injection and negative values after steam, with δ¹¹B as low as ? 28‰. Other possibly reactive phases include clay minerals and organic matter, but their abundance is not great enough to impact on the isotopic composition of the produced waters. This information makes it possible to evaluate the boron isotope fractionation equation derived from experimental data ( Williams LB (2000) Boron isotope geochemistry during burial diagenesis. PhD Dissertation. The University of Calgary, Alberta, Canada; Williams LB, Hervig RL, Holloway JR, Hutcheon I (2001a) Boron isotope geochemistry during diagenesis: Part 1. Experimental determination of fractionation during illitization of smectite. Geochimica et Cosmochimica Acta, in press). The results show that the fractionation curve predicts the difference between δ¹¹B of the pumice and hydrothermal fluids in the Cold Lake reservoir. This not only indicates that the reservoir fluids have approached boron isotope equilibrium with the reservoir rock, but also shows that B isotopes provide a useful geothermometer for hydrothermally stimulated oil reservoirs.     

Principal features of impact-generated hydrothermal circulation systems: mineralogical and geochemical evidence

Any hypervelocity impact generates a hydrothermal circulation system in resulting craters. Common characteristics of hydrothermal fluids mobilized within impact structures are considered, based on mineralogical and geochemical investigations, to date. There is similarity between the hydrothermal mineral associations in the majority of terrestrial craters; an assemblage of clay minerals–zeolites–calcite–pyrite is predominant. Combining mineralogical, geochemical, fluid inclusion, and stable isotope data, the distinctive characteristics of impact‐generated hydrothermal fluids can be distinguished as follows: (i) superficial, meteoric and ground water and, possibly, products of dehydration and degassing of minerals under shock are the sources of hot water solutions; (ii) shocked target rocks are sources of the mineral components of the solutions; (iii) flow of fluids occurs mainly in the liquid state; (iv) high rates of flow are likely (10^?4 to 10^?3 m s^?1); (v) fluids are predominantly aqueous and of low salinity; (vi) fluids are weakly alkaline to near‐neutral (pH 6–8) and are supersaturated in silica during the entire hydrothermal process because of the strong predominance of shock‐disordered aluminosilicates and fusion glasses in the host rocks; and (vii) variations in the properties of the circulating solutions, as well as the spatial distribution of secondary mineral assemblages are controlled by temperature gradients within the circulation cell and by a progressive cooling of the impact crater. Products of impact‐generated hydrothermal processes are similar to the hydrothermal mineralization in volcanic areas, as well as in modern geothermal systems, but impacts are always characterized by a retrograde sequence of alteration minerals.     

Indications of formation water flow and compartmentalization on the flank of a salt structure derived from salinity and seismic data

Vertical and lateral variations in lithology, salinity, temperature, and pressure determined from wireline LAS logs, produced water samples, and seismic data on the south flank of a salt structure on the continental shelf, offshore Louisiana indicate three hydrogeologic zones in the study area: a shallow region from 0 to 1.1 km depth with hydrostatically pressured, shale‐dominated Pleistocene age sediments containing pore waters with sea water (35 g l^?1) or slightly above sea water salinity; a middle region from 1.1 to 3.2 km depth with near hydrostatically pressured, sand‐dominated Pliocene age sediments that contain pore waters that range from seawater salinity to up to 5 times sea water salinity (180 g l^?1); and a deep section below 3.2 km depth with geopressured, shale‐dominated Miocene age sediments containing pore waters that range from sea water salinity to 125 g l^?1. Salt dissolution has generated dense, saline waters that appear to be migrating down dip preferentially through the thick Pliocene sandy section. Sand layers that come in contact with salt contain pore waters with high salinity. Isolated sands have near sea water salinity. Salinity information in conjunction with seismic data is used to infer fluid compartmentalization. Both vertical and lateral lithologic barriers to fluid flow at tens to hundreds of meters scale are observed. Fluid compartmentalization is also evident across a supradomal normal fault. Offset of salinity contours are consistent with the throw of the fault, which suggests that saline fluids migrated before fault formation.     

Temperature‐dependent Li isotope ratios in Appalachian Plateau and Gulf Coast Sedimentary Basin saline water

Lithium (Li) concentrations of produced water from unconventional (horizontally drilled and hydraulically fractured shale) and conventional gas wells in Devonian reservoirs in the Appalachian Plateau region of western Pennsylvania range from 0.6 to 17 mmol kg^?1, and Li isotope ratios, expressed as in δ⁷Li, range from +8.2 to +15‰. Li concentrations are as high as 40 mmol kg^?1 in produced waters from Plio‐Pleistocene through Jurassic‐aged reservoirs in the Gulf Coast Sedimentary Basin analyzed for this study, and δ⁷Li values range from about +4.2 to +16.6‰. Because of charge‐balance constraints and rock buffering, Li concentrations in saline waters from sedimentary basins throughout the world (including this study) are generally positively correlated with chloride (Cl), the dominant anion in these fluids. Li concentrations also vary with depth, although the extent of depth dependence differs among sedimentary basins. In general, Li concentrations are higher than expected from seawater or evaporation of seawater and therefore require water–mineral reactions that remove lithium from the minerals. Li isotope ratios in these produced waters vary inversely with temperature. However, calculations of temperature‐dependent fractionation of δ⁷Li between average shale δ⁷Li (?0.7‰) and water result in δ⁷Li_water that is more positive than that of most produced waters. This suggests that aqueous δ⁷Li may reflect transport of water from depth and/or reaction with rocks having δ⁷Li lighter than average shale.     

The sources and formation processes of brines from the Lunnan Ordovician paleokarst reservoir,Tarim Basin,northwest China

The most important petroleum exploration target in the Tarim Basin, northwest China, is the paleokarst reservoir. To understand the source and evolution of brine in this type of reservoir, a total of 37 formation‐water samples were collected from the Middle‐Lower Ordovician paleokarst reservoir in the Lunnan oilfield. The δD‐δ¹⁸O correlation and Cl/Br ratios reflect the mixture of two fluids: meteoric water and evaporated seawater. The different degree of mixture divided samples into two groups. Group 1 samples, from deep strata (5150–6667 m.b.s.l.) in the east of the field, with elevated δD (?53.5 to ?38.0‰), δ¹⁸O values (0.66–5.99‰), and lower Cl/Br ratios (336–478 for Cl/Br, except LN634‐1 and LN631‐1) were formed by evaporation of seawater plus a small contribution from meteoric water. Group 2 samples, from shallow strata (5038–6067 m.b.s.l.), in the west of the field, have contrasting features (?59.6 to ?48.5‰ for δD, ?0.47 to 2.17‰ for δ¹⁸O, and 501 to 871 for Cl/Br), which reflect a mixture of evaporated seawater with a high proportion of meteoric water. Both of the fluid types exchanged oxygen isotope with minerals. The investigation into cation composition reveals that, before entering into the current reservoir, waters suffered albitization of plagioclase; moreover, meteoric water dissolved evaporites and seawater experienced dolomitization. A mixing trend showed by strontium isotopes (0.709801–0.711628) gave further evidence for the mixture of two fluid types. Based on the correlation of geological history with our data, two infiltration models of meteoric waters can be constructed. According to the chemical and isotopic compositions of the waters, an east fluid regime (Group 1) and a west fluid regime (Group 2) have thus been defined. Better understanding of the subsurface fluid movement patterns may be helpful for the local exploration.     

Structure of grain boundaries in wet,synthetic polycrystalline,statically recrystallizing halite – evidence from cryo‐SEM observations

It is well known from nature and experiments that the presence of brine strongly affects the microstructural evolution and the mechanical and transport properties of halite. Existing interpretations of the grain boundary structure in deformed, wet, salt samples annealed statically at room temperature are based on indirect evidence from reflected light microscopy and conventional scanning electron microscopy. This paper presents direct observations of fluid‐filled grain boundaries using the cryogenic‐scanning electron microscope (cryo‐SEM) in which the grain boundary fluids were frozen before breaking the samples. The rapid cooling transforms the brine into two phases, i.e. ice and hydrohalite, which are easily recognized from characteristic segregation patterns. We studied samples of wet, synthetic, polycrystalline halite annealed under static conditions at room temperature. In coarse‐grained samples, fine‐scale segregation patterns were observed at the boundaries of the primary recrystallizing grains. These points indicate the existence of fluid films with a thickness in the range of 30 nm, but the finer scale structure of the fluid remains unknown. In fine‐grained samples, the distribution and reorganization of fluids with annealing time is recorded by the combination of contact healing and successive accumulation of fluids in triple junction tubes. The contact healing is attributed to the small initial grain size, such that the fluid film necks down by accumulating the fluids into previously existing triple junctions via neck growth. Electron backscatter diffraction measurements of both primary and secondary recrystallized grains indicate that they are euhedral, i.e. the grain growth morphology is controlled by the anisotropy of the grain boundary energy of the growing grain, which results in planar growth faces.     

A fluid inclusion record of magmatic/hydrothermal pulses in acid Salar Ignorado gypsum,northern Chile

Salar Ignorado is a shallow acid saline lake hosted by a small intervolcanic basin high in the Andes Mountains of northern Chile. Modern surface waters have 3.3–4.1 pH, 0.5–3% total dissolved solids (TDS) and are actively precipitating gypsum crystals. The gypsum crystals trap the acid saline water as fluid inclusions, providing a record of recent surface water characteristics. Salar Ignorado gypsum contains three distinct types of primary fluid inclusions, which result from growth of the gypsum from surface waters. Petrography and microthermometry were performed on 27 gypsum crystals from Salar Ignorado to gain an understanding of recent water chemistry of the salar. One 18.3‐cm‐long gypsum crystal, hosting primary fluid inclusions along 28 successive growth bands, was the focus for fluid inclusion studies and allowed a record of high‐resolution chemical trends. This crystal showed a change in parent fluids during growth, from low salinity, to high salinity, back to low salinity. At the bottom of the crystal, the lowest six fluid inclusion assemblages have salinities of 1.7–5.1 eq. wt. % NaCl. The next nine fluid inclusion assemblages have significantly higher salinity (18.6–27.4 eq. wt. % NaCl) inclusions. The twelve fluid inclusion assemblages near the top of the crystal have low salinity (0.9–8.3 eq. wt. % NaCl) like those at the bottom of the crystal. The high‐salinity fluid inclusions in the middle of this gypsum crystal are interpreted to have formed during a pulse of magmatic/hydrothermal fluids to the surface, perhaps during local active volcanism. Secondary evidence of a magmatic influence on surface waters includes hydrogen sulfide and high molecular weight solid hydrocarbons within some fluid inclusions. This study is among the first detailed fluid inclusion studies of gypsum and suggests that fluid inclusions in gypsum can be paleo‐hydrogeologic proxies.     

Pore‐scale analysis of trapped immiscible fluid structures and fluid interfacial areas in oil‐wet and water‐wet bead packs

The objective of this study is to obtain quantitative evidence of pore‐scale immiscible fluid distribution in oil‐wet and water‐wet porous media using X‐ray computed microtomography. Temporal and spatial saturation profiles, as well as surface and interfacial areas, are thoroughly analyzed through cycles of drainage and imbibition using samples with different wetting characteristics but similar pore structures. The population of individual immiscible fluid structures (‘blobs’) was also evaluated. The specific nonwetting phase surface areas of both porous media are found to be in close correlation with the specific solid surface area. On the other hand, the differing wetting strengths of the two porous media affect the curvature of the fluid–fluid interface and thus the specific meniscus interfacial area of the two porous media. Although the magnitude of the specific meniscus interfacial areas is different, they both trend toward a maximum at wetting phase saturations of 0.35–0.55. The differences in wetting characteristics are also apparent in the blob populations. The number of blobs in the oil‐wet porous media is three times greater than that of the water‐wet porous media at similar saturations; the increase in population is a result of the increase in the amount of smaller blobs inhabiting the smaller pore spaces. The surface areas of individual blobs as a function of the individual blob volumes are found to closely agree with the specific surface area of a sphere at blob volumes below the minimum individual grain volume and with the specific pore space surface area above this volume. These results show how wettability and saturation history influence the distribution of immiscible fluids within the pore space.     

The origins of salinity in metamorphic fluids

Evaluation of data on formation waters and metamorphic fluids sampled by drilling or preserved in fluid inclusions reveals little correlation between fluid salinity and metamorphic grade, but a strong link to original sedimentary setting. Sediments and metasediments deposited originally in shallow marine environments can contain fluids with a very wide range of salinities, but they are commonly near twice seawater salinity or higher. With increasing metamorphic grade, a very wide range of salinities may develop, with the highest levels tracking halite saturation. Oceanic and accretionary prism sequences yield low‐salinity fluids, close to seawater values, almost irrespective of metamorphic grade until extreme conditions are reached where removal of water may increase fluid salinity. The salinities of metamorphic fluids exert a fundamental control on both fluid phase equilibria and metal‐transporting capability, and appear, to a large degree, to reflect the original presence or absence of highly saline formation waters and/or evaporites in the initial sedimentary sequence.     

Boron isotope geochemistry of Na‐bicarbonate,Na‐chloride,and Ca‐chloride waters from the Northern Apennine Foredeep basin: other pieces of the sedimentary basin puzzle

The boron stable isotope ratio δ¹¹B of 12 water samples representative of three chemical facies (fresh Na‐bicarbonate, brackish Na‐chloride, saline, and brine Ca‐chloride) has been analyzed. Interpretation of the δ¹¹B data, along with the chemical compositions, reveals that Na‐carbonate waters from the Northern Apennine are of meteoric origin, with boron contributions from clay desorption and mixing with seawater‐derived fluids of Na‐chloride or Ca‐chloride compositions. The comparison of our new results with the literature data on other sedimentary basins of Mediterranean, and worldwide, confirms the contribution of Na‐bicarbonate waters to the genesis of mud volcano fluids. The Na‐chloride sample of Salvarola (SAL), which may represent the end‐member of the mud volcanoes, and the Ca‐chloride brine water from Salsomaggiore (SM) indicate boron release from clays compatible with the diagenetic process. The empirical equation: relating boron concentration and the stable isotope composition of the fluids observed in this study and the literature is proposed to trace the effect of diagenesis in sedimentary basins. A geothermometer associated to the diagenetic equation is also proposed: The application of this equation to obtain reservoir temperatures from δ¹¹B compositions of waters should be carefully evaluated against the results obtained from other chemical and isotopic geothermometers from other basins around the world.     

Deep geothermal groundwater flow in the Seferihisar–Balçova area,Turkey: results from transient numerical simulations of coupled fluid flow and heat transport processes

The Seferihisar–Balçova Geothermal system (SBG) is characterized by complex temperature and hydrochemical anomalies. Previous geophysical and hydrochemical investigations suggest that hydrothermal convection in the faulted areas of the SBG and recharge flow from the Horst may be responsible for the observed patterns. A numerical model of coupled fluid flow and heat transport processes has been built in order to study the possible fluid dynamics of deep geothermal groundwater flow in the SBG. The results support the hypothesis derived from interpreted data. The simulated scenarios provide a better understanding of the geophysical conditions under which the different fluid dynamics develop. When recharge processes are weak, the convective patterns in the faults can expand to surrounding reservoir units or below the seafloor. These fault‐induced drag forces can cause natural seawater intrusion. In the Melange of the Seferihisar Horst, the regional flow is modified by buoyant‐driven flow focused in the series of vertical faults. As a result, the main groundwater divide can shift. Sealing caprocks prevent fault‐induced cells from being overwhelmed by vigorous regional flow. In this case, over‐pressured, blind geothermal reservoirs form below the caprocks. Transient results showed that the front of rising hot waters in faults is unstable: the tip of the hydrothermal plumes can split and lead to periodical temperature oscillations. This phenomenon known as Taylor–Saffman fingering has been described in mid‐ocean ridge hydrothermal systems. Our findings suggest that this type of thermal pulsing can also develop in active, faulted geothermal systems. To some extent, the role of an impervious fault core on the flow patterns has been investigated. Although it is not possible to reproduce basin‐scale transport processes, this first attempt to model deep groundwater geothermal flow in the SBG qualitatively supported the interpreted data and described the different fluid dynamics of the basin. Geofluids (2010) 10 , 388–405     

Evaluation of recoverable energy potential from enhanced geothermal systems: a sensitivity analysis in a poro‐thermo‐elastic framework

This study presents application of an efficient approach to simulate fluid flow and heat transfer in naturally fractured geothermal reservoirs. Fluid flow is simulated by combining single continuum and discrete fracture approaches. The local thermal nonequilibrium approach is used to simulate heat transfer by conduction in the rock matrix and convection (including conduction) in the fluid. Fluid flow and heat transfer models are integrated within a coupled poro‐thermo‐elastic framework. The developed model is used to evaluate the long‐term response of a geothermal reservoir with specific boundary conditions and injection/production schedule. A comparative study and a sensitivity analysis are carried out to evaluate the capability of the integrated approach and understand the degree by which different reservoir parameters affect thermal depletion of Soultz geothermal reservoir, respectively. Also observed, there exists an optimum fracture permeability after which the reservoir stimulation does not change the recovery factor significantly. Estimation of fluid temperature by the assumption of local thermal nonequilibrium heat transfer between the fracture fluid and the rock matrix gives fluid temperature of about 3°C less than that of estimated by thermal equilibrium heat transfer at early stage of hot water production.     

Thermal waters of ‘tectonic origin’: the alkaline,Na‐HCO3 waters of the Rio Valdez geothermal area (Isla Grande de Tierra del Fuego,Argentina)

The geothermal area of Rio Valdez is located in the central portion of the Isla Grande de Tierra del Fuego (South Argentina), ten kilometers south of the southeastern sector of the Fagnano Lake. It consists of a series of thermal springs with low discharge rates (≤1 L/s) and temperatures in the range of 20–33°C distributed in an area of <1 km². The thermal springs are characterized by alkaline, Na‐HCO₃ waters with low salinity (0.53÷0.58 g/L), but relatively high fluoride contents (up to 19.4 mg/L). Their composition is the result of a slow circulation at depth, possibly through deep tectonic discontinuities connected with the Magallanes‐Fagnano Fault (MFF) system. According to geothermometric calculations, thermal waters reach temperatures in the range of 100–150°C and an almost complete chemical equilibrium with the alkali‐feldspars in the metavolcanic country rocks. The relatively high fluorine contents can be explained by the slow ascent and cooling of deep groundwaters followed by a progressive re‐equilibration with F‐bearing, hydrated Mg‐silicates, such as chlorite, which has been recognized as an abundant mineral in the metavolcanics of the Lemaire Formation and metapelites and metagraywackes of the Yahgán Formation. Finally, the isotopic composition of the investigated samples is consistent with the infiltration from local snow melting at altitudes in the range of 610–770 m asl. The comparison of our data with those collected in 1991 seems to suggest a possible progressive decline of the bulk thermal output in the near future. This possibility should be seriously considered before planning a potentially onerous exploitation of the resource. Presently, the only ways to exploit this geothermal resource by the population scattered in the area are the direct use of thermal waters and/or spa structures.