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 upper continental crust, an aquifer and its fluid: hydaulic and chemical data from 4 km depth in fractured crystalline basement rocks at the KTB test site

Detailed information on the hydrogeologic and hydraulic properties of the deeper parts of the upper continental crust is scarce. The pilot hole of the deep research drillhole (KTB) in crystalline basement of central Germany provided access to the crust for an exceptional pumping experiment of 1‐year duration. The hydraulic properties of fractured crystalline rocks at 4 km depth were derived from the well test and a total of 23100 m³ of saline fluid was pumped from the crustal reservoir. The experiment shows that the water‐saturated fracture pore space of the brittle upper crust is highly connected, hence, the continental upper crust is an aquifer. The pressure–time data from the well tests showed three distinct flow periods: the first period relates to wellbore storage and skin effects, the second flow period shows the typical characteristics of the homogeneous isotropic basement rock aquifer and the third flow period relates to the influence of a distant hydraulic border, probably an effect of the Franconian lineament, a steep dipping major thrust fault known from surface geology. The data analysis provided a transmissivity of the pumped aquifer T = 6.1 × 10^?6 m² sec^?1, the corresponding hydraulic conductivity (permeability) is K = 4.07 × 10^?8 m sec^?1 and the computed storage coefficient (storativity) of the aquifer of about S = 5 × 10^?6. This unexpected high permeability of the continental upper crust is well within the conditions of possible advective flow. The average flow porosity of the fractured basement aquifer is 0.6–0.7% and this range can be taken as a representative and characteristic values for the continental upper crust in general. The chemical composition of the pumped fluid was nearly constant during the 1‐year test. The total of dissolved solids amounts to 62 g l^?1 and comprise mainly a mixture of CaCl₂ and NaCl; all other dissolved components amount to about 2 g l^?1. The cation proportions of the fluid (X_Ca approximately 0.6) reflects the mineralogical composition of the reservoir rock and the high salinity results from desiccation (H₂O‐loss) due to the formation of abundant hydrate minerals during water–rock interaction. The constant fluid composition suggests that the fluid has been pumped from a rather homogeneous reservoir lithology dominated by metagabbros and amphibolites containing abundant Ca‐rich plagioclase.     

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.     

Mineralogy and fluid inclusion gas chemistry of production well mineral scale deposits at the Dixie Valley geothermal field,USA

At the Dixie Valley geothermal field, Nevada, USA, fluid boiling triggered the precipitation of carbonate scale minerals in concentric bands around tubing inserted into production well 28–33. When the tubing was removed, this mineral scale was sampled at 44 depth intervals between the wellhead and 1227 m depth. These samples provide a unique opportunity to evaluate the effects of fluid boiling on the scale mineralogy and geochemistry of the vapor and liquid phase. In this study, the mineralogy of the scale deposits and the composition of the fluid inclusion gases trapped in the mineral scales were analyzed. The scale consists mainly of calcite from 670–1112 m depth and aragonite from 1125 to 1227 m depth, with traces of quartz and Mg‐smectite. Mineral textures, including hopper growth, twinning, and fibrous growth in the aragonite and banded deposits of fine grained calcite crystals, are the result of progressive boiling. The fluid inclusion noncondensable gas was dominated by CO₂. However, significant variations in He relative to N₂ and Ar provide evidence that the geothermal reservoir consists of mixed source deeply circulating reservoir water and shallow, air saturated meteoric water. Gas analyses for many inclusions also showed higher CH₄ and H₂ relative to CO₂ than measured in gas sampled from this well, other production wells, and fumaroles. These inclusions are interpreted to have trapped CH₄‐ and H₂‐enriched gas resulting from early stages of boiling.     

Mathematical modeling of phase separation of seawater near an igneous dike

We provide a simplified treatment of phase separation of seawater near an igneous dike to obtain rough estimates of the thickness and duration of the two‐phase zone, the volume fractions of vapor and brine formed, and their distribution in the subsurface. Under the assumption that heat transfer occurs mainly by thermal conduction we show that, for a 2‐m wide dike, the maximum width of the two phase zone is approximately 20 cm and that a zone of halite is initially deposited near the dike wall. The two‐phase zone is mainly filled with vapor. For a value of thermal diffusivity of a = 10^?6 m² sec^?1, the two‐phase zone begins to disappear at the base of the system after 13 days, and disappears completely by 16 days. For a lower value of thermal diffusivity, the width of the two‐phase region does not change appreciably but its duration increases as a^?1. The width of the two‐phase zone determined by this simplified model agrees reasonably well with transient numerical solutions for the analogous two‐phase flow in a pure water system; however the duration of two‐phase flow is matched better using a smaller value of a. We compare the seafloor values of vapor salinity and temperature given by the model with vapor salinity data from the ‘A’ vent at 9–10°N on the East Pacific Rise (EPR) and argue that either non‐equilibrium thermodynamic behavior or near‐surface mixing of brine with vapor in the two‐phase region may explain the discrepancies between model predictions and data.     

Salt tectonics and shallow subseafloor fluid convection: models of coupled fluid‐heat‐salt transport

Thermohaline convection associated with salt domes has the potential to drive significant fluid flow and mass and heat transport in continental margins, but previous studies of fluid flow associated with salt structures have focused on continental settings or deep flow systems of importance to petroleum exploration. Motivated by recent geophysical and geochemical observations that suggest a convective pattern to near‐seafloor pore fluid flow in the northern Gulf of Mexico (GoMex), we devise numerical models that fully couple thermal and chemical processes to quantify the effects of salt geometry and seafloor relief on fluid flow beneath the seafloor. Steady‐state models that ignore halite dissolution demonstrate that seafloor relief plays an important role in the evolution of shallow geothermal convection cells and that salt at depth can contribute a thermal component to this convection. The inclusion of faults causes significant, but highly localized, increases in flow rates at seafloor discharge zones. Transient models that include halite dissolution show the evolution of flow during brine formation from early salt‐driven convection to later geothermal convection, characteristics of which are controlled by the interplay of seafloor relief and salt geometry. Predicted flow rates are on the order of a few millimeters per year or less for homogeneous sediments with a permeability of 10^?15 m², comparable to compaction‐driven flow rates. Sediment permeabilities likely fall below 10^?15 m² at depth in the GoMex basin, but such thermohaline convection can drive pervasive mass transport across the seafloor, affecting sediment diagenesis in shallow sediments. In more permeable settings, such flow could affect methane hydrate stability, seafloor chemosynthetic communities, and the longevity of fluid seeps.     

Retrograde P–T evolution and high temperature–low pressure fluid circulation in relation to late Hercynian intrusions: a mineralogical and fluid inclusion study of the Charroux-Civray plutonic complex (north-western Massif Central, France)

The calc‐alkaline plutonic complex from Charroux‐Civray (north‐western part of the French Massif Central) displays multiphase hydrothermal alteration. Plutonic rocks, as well as early retrograde Ca–Al silicate assemblages, which have crystallized during cooling and uplifting of the plutonic series, are affected by multiphase chlorite–phengite–illite–carbonate alteration linked to intense pervasive fluid circulation through microfractures. The petrographic study of alteration sequences and their associated fluid inclusions in microfissures of the plutonic rocks, as well as in mineral fillings of the veins, yields a reconstruction of the P–T–X evolution of the Hercynian basement after the crystallization of the main calc‐alkaline plutonic bodies. This reconstruction covers the uplift of the basement to its exposure and the subsequent burial by Mesozoic sediments. Cooling of the calc‐alkaline plutonic series started at solidus temperatures (～650°C), at a pressure of about 4 kbar (1 bar = 10⁵ N m^?2), as indicated by magmatic epidote stability, hornblende barometry and fluid inclusion data. Cooling continued under slightly decreasing pressure during uplift down to 2–3 kbar at 200–280°C (prehnite–pumpellyite paragenesis). Then, a hot geothermal circulation of CO₂‐bearing fluids was induced within the calc‐alkaline rocks leading to the formation of greisen‐like mineralizations. During this stage, temperatures around 400–450°C were still high for the inferred depths (～2 kbar). They imply abnormal heat flows and thermal gradients of 60–80°C km^?1. The hypothesis of the existence of one large or a succession of smaller peraluminous plutons at depth, supported by geophysical data, suggests that localized heat flows were linked to concealed leucogranite intrusions. As uplift continued, greisen mineralization was subsequently affected by the chlorite–phengite–dolomite assemblage, correlated with aqueous and nitrogen‐bearing fluid circulations in the temperature range of 400–450°C. In a later stage, a continuous temperature decrease at constant pressure (～0.5 kbar) led to the alteration of the dolomite–illite–chlorite type in the 130–250°C temperature range.     

Br/Cl signature of hydrothermal fluids: liquid–vapour fractionation of bromine revisited

Br/Cl ratios of hydrothermal fluids are widely used as geochemical tracers in marine hydrothermal systems to prove fluid phase separation processes. However, previous results of the liquid–vapour fractionation of bromine are ambiguous. Here we report new experimental results of the liquid–vapour fractionation of bromine in the system H₂O–NaCl–NaBr at 380–450°C and 22.9–41.7 MPa. Our data indicate that bromine is generally more enriched than chlorine in the liquid phase. Calculated exchange coefficients K_D(Br‐Cl)^{liquid‐vapour} for the reaction Br^vapour + Cl^liquid = Br^liquid + Cl^vapour are between 0.94 ± 0.08 and 1.66 ± 0.14 within the investigated P–T range. They correlate positively with D_Cl^{liquid‐vapour} and suggest increasing bromine–chlorine fractionation with increasing opening of the liquid–vapour solvus, i.e. increasing distance to the critical curve in the H₂O–NaCl system. An empirical fit of the form K_D(Br‐Cl)^{liquid‐vapour} = a*ln[b*(D_Cl^{liquid‐vapour}?1) + e^1/a] yields a = 0.349 and b = 1.697. Based on this empirical fit and the well‐constrained phase relations in the H₂O–NaCl system we calculated the effect of fluid phase separation on the Br/Cl signature of a hydrothermal fluid with initial seawater composition for closed and open adiabatic ascents along the 4.5 and 4.8 J g^?1 K^?1 isentropes. The calculations indicate that fluid phase separation can significantly alter the Br/Cl ratio in hydrothermal fluids. The predicted Br/Cl evolutions are in accord with the Br/Cl signatures in low‐salinity vent fluids from the 9 to 10°N East Pacific Rise.     

Fracture-fill calcite as a record of microbial methanogenesis and fluid migration: a case study from the Devonian Antrim Shale, Michigan Basin

The Devonian Antrim Shale is an organic‐rich, naturally fractured black shale in the Michigan Basin that serves as both a source and reservoir for natural gas. A well‐developed network of major, through‐going vertical fractures controls reservoir‐scale permeability in the Antrim Shale. Many fractures are open, but some are partially sealed by calcite cements that retain isotopic evidence of widespread microbial methanogenesis. Fracture filling calcite displays an unusually broad spectrum of δ¹³C values (+34 to ?41‰ PDB), suggesting that both aerobic and anaerobic bacterial processes were active in the reservoir. Calcites with high δ¹³C values (>+15‰) record cementation of fractures from dissolved inorganic carbon (DIC) generated during bacterial methanogenesis. Calcites with low δ¹³C values (13C values between ?10 and ?30‰ can be attributed to variable organic matter oxidation pathways, methane oxidation, and carbonate rock buffering. Identification of ¹³C‐rich calcite provides unambiguous evidence of biogenic methane generation and may be used to identify gas deposits in other sedimentary basins. It is likely that repeated glacial advances and retreats exposed the Antrim Shale at the basin margin, enhanced meteoric recharge into the shallow part of the fractured reservoir, and initiated multiple episodes of bacterial methanogenesis and methanotrophic activity that were recorded in fracture‐fill cements. The δ¹⁸O values in both formation waters and calcite cements increase with depth in the basin (?12 to ?4‰ SMOW, and +21 to +27‰ PDB, respectively). Most fracture‐fill cements from outcrop samples have δ¹³C values between ?41 and ?15‰ PDB. In contrast, most cement in cores have δ¹³C values between +15 and +34‰ PDB. Radiocarbon and ²³⁰Th dating of fracture‐fill calcite indicates that the calcite formed between 33 and 390 ka, well within the Pleistocene Epoch.     

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.     

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.     

The link between fluids and rank variation in the South Wales Coalfield: evidence from fluid inclusions and stable isotopes

Fluid inclusion and stable isotope data from quartz and carbonate minerals in fracture fillings and ‘ironstone’ nodules from the South Wales Coalfield have been used to characterise the fluids generated during basin evolution and associated coalification. Carbonates grew first, probably at relatively shallow depths and low temperatures (<100°C). The carbonates exhibit a trend of increasing C‐isotopic values across the coalfield, ranging from δ¹³C = ?12‰ VPDB in the SE of the coalfield to 0‰ VPDB in the NW, possibly as a result of increasing methanogenesis in the deeper (NW) parts of the coalfield. Quartz formed at a later stage of basin formation, probably at temperatures between 150 and 200°C. Fluid inclusions in these minerals suggest that burial and coalification of the sediments were associated with mixed aqueous–petroleum fluids. Furthermore, the density of these petroleum fluids decreases towards the NW of the coalfield, where the rank of the associated coal increases to anthracite grade. The study confirms that the composition and temperature of these fluids closely correlate with the variations in coal rank, indicating a possible causal link. The data also give general support to models that propose regional fluid flow in the basin. and are consistent with the erosion of approximately 2 km of section which is not preserved today. A geothermal gradient (at maximum burial) of 45°C km^?1 is proposed, and thus no exceptionally anomalous thermal regime is required to explain coal rank variation.     

Structural controls on hydrothermal flow in a segmented rift system, Taupo Volcanic Zone, New Zealand

Vigorous hydrothermal convection transfers 10 times the average continental heat flow through the central Taupo Volcanic Zone (TVZ), a region of active extension (approximately 8 mm year^?1) and productive rhyolitic volcanism. Over 20 high‐temperature (>250°C) geothermal fields occur within Quaternary pyroclastic basins, with convective circulation to depths of 7–8 km presumably extending through basement rocks. Parallel‐striking normal faults, fractures and dikes dissect the convective regime, interacting with fluids to either enhance or restrict flow according to the relative permeability of structure and host rock. In the basement, high bulk permeability is maintained by focussed flow through faults and associated fractures well oriented for reactivation in the prevailing stress field. In contrast, distributed flow through fault‐bounded compartments prevails within Quaternary basins, masking any signal of deeper structural control. Exceptions occur where more competent rocks are exposed at the surface. As in narrow magmatic rifts elsewhere, the extensional fabric is partitioned into discrete rift segments linked along strike by accommodation zones. Eighty per cent of TVZ geothermal fields correlate spatially with rift architecture, with 60% located in accommodation zones. We suggest that segmented rift fabrics generate bulk permeability anisotropy that is to some extent predictable, with rift segments characterized by enhanced axial flow, and accommodation zones characterized by locally enhanced vertical permeability that is tectonically maintained. This provides a plausible explanation for the common occurrence of geothermal fields within accommodation zones and their notable absence within densely faulted rift segments. Maintenance of structural permeability in zones of active hydrothermal precipitation necessarily requires repeated brittle failure. Geothermal plumes therefore exploit tectonically maintained permeability within accommodation zones, with rift segments functioning mostly as drawdown regions. The influence of rift architecture on flow paths has important implications for geothermal extraction and epithermal mineral exploration within the TVZ and other structurally segmented hydrothermal systems, both active and extinct.     

Molecular transport of methane, ethane and nitrogen and the influence of diffusion on the chemical and isotopic composition of natural gas accumulations

Laboratory experiments have been performed to determine diffusion coefficients of natural gas components (methane, ethane and nitrogen) and isotope fractionation effects under simulated in situ pressure (up to 45 MPa effective stress) and temperature conditions (50–200°C) in water‐saturated pelitic and coarse‐grained rocks. Effective diffusion coefficients of molecular nitrogen (0.39 × 10^?11 to 21.6 × 10^?11 m² sec^?1 at 90°C) are higher than those for methane (0.18 × 10^?11 to 18.2 × 10^?11 m² sec^?1 at 90°C). Diffusive flux rates expressed in mass units are generally higher for N₂ than for CH₄. Both methane and (to a lesser extent) nitrogen diffusion coefficients decrease with increasing total organic carbon (TOC) content of the rock samples because of sorption processes on the organic matter. This effect decreases with increasing temperature. Effective diffusion coefficients increase upon a temperature increase from 50 to 200°C by a factor of four. Effective diffusion coefficients and steady‐state diffusive flux decrease with effective stress. Stationary diffusive fluxes drop by 50–70% for methane and 45–62% for nitrogen while effective diffusion coefficients are reduced by 38% (CH₄) and 32–48% (N₂), respectively. Isotope fractionation coefficients of diffusive transport are higher for methane (?1.56 and ?2.77‰) than for ethane (?0.84 and ?1.62‰). Application of the experimental results to geological systems show that diffusive transport has only a low transport efficiency. Significant depletion of natural gas reservoirs by molecular diffusion is only expected in cases of very poor caprock qualities (in terms of thickness and/or porosity) and over extended periods of geological time. Under these circumstances, the chemical and isotopic composition of a gas reservoir will change and maturity estimates based on these parameters may be deceptive. To account for these potential effects, nomograms have been developed to estimate diffusive losses and apply maturity corrections.     

Geochemical fingerprinting of hydraulic fracturing fluids from Qusaiba Hot Shale and formation water from Paleozoic petroleum systems,Saudi Arabia

Provenance studies of produced water are essential to trace flow dynamics and reservoir compartmentalization in petroleum systems and to quantify fluid recovery rates from unconventional fracturing. Produced water from a hydraulically fractured well in the Qusaiba Hot Shale in the Northern Exploration Area, Saudi Arabia, was daily monitored and analyzed for water chemistry, and environmental (δ²H, δ¹³C, δ¹⁸O_H2O, δ¹⁸O_SO4, δ³⁴S_SO4, δ³⁷Cl, ⁸⁷Sr/⁸⁶Sr) and cosmogenic isotopes (³H, ¹⁴C, ³⁶Cl), to differentiate from reference fluids of supply water, fracturing fluids, and formation water from adjacent Paleozoic units. Initially, recovered water is composed of fracturing fluids and subsequently replaced by a homogeneous cut of pristine formation water. Formation water is composed of dominant meteoric water (approximately 84 vol%) and minor fossil evaporated seawater. The young ¹⁴C‐apparent age between 6000 and 6700 years BP and depleted δ¹⁸O/δ²H values for the meteoric component confirm the infiltration of surface water into the Qusaiba Hot Shale interval or adjacent units during the Early Holocene Pluvial Period under cooler and wetter climatic conditions than present, which suggest the presence of a very recent, dynamic hydraulic flow system. ³⁶Cl/Cl ratios between 102 × 10^?15 and 31 × 10^?15 are ambiguous and can be attributed to atmospheric recharge close to the coast, mixing of ³⁶Cl‐enriched Quaternary meteoric recharge with ³⁶Cl‐depleted fossil seawater, and/or hypogene production by U‐Th‐enriched host rock. Produced waters from Qusaiba Hot Shale are within the compositional range of Na‐Cl‐type formation water from Paleozoic reservoir units in northern Saudi Arabia with salinities from 30 000 to 130 000 mg l^?1. As a novel technological approach for exploration wells in Northern Saudi Arabia, multi‐isotopic methods were successfully implemented to quantify flowback volumes from hydraulic fracturing, and to fingerprint pristine formation water or pore water in Paleozoic systems on their provenance, residence time, migration pathways, and secondary alteration processes.     

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     

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.     

Simulation of the impact of faults on CO2 injection into sandstone reservoirs

Numerical simulations of multiphase CO₂ behavior within faulted sandstone reservoirs examine the impact of fractures and faults on CO₂ migration in potential subsurface injection systems. In southeastern Utah, some natural CO₂ reservoirs are breached and CO₂‐charged water flows to the surface along permeable damage zones adjacent to faults; in other sites, faulted sandstones form barriers to flow and large CO₂‐filled reservoirs result. These end‐members serve as the guides for our modeling, both at sites where nature offers ‘successful’ storage and at sites where leakage has occurred. We consider two end‐member fault types: low‐permeability faults dominated by deformation‐band networks and high‐permeability faults dominated by fracture networks in damage zones adjacent to clay‐rich gouge. Equivalent permeability (k) values for the fault zones can range from <10^?14 m² for deformation‐band‐dominated faults to >10^?12 m² for fracture‐dominated faults regardless of the permeability of unfaulted sandstone. Water–CO₂ fluid‐flow simulations model the injection of CO₂ into high‐k sandstone (5 × 10^?13 m²) with low‐k (5 × 10^?17 m²) or high‐k (5 × 10^?12 m²) fault zones that correspond to deformation‐band‐ or fracture‐dominated faults, respectively. After 500 days, CO₂ rises to produce an inverted cone of free and dissolved CO₂ that spreads laterally away from the injection well. Free CO₂ fills no more than 41% of the pore space behind the advancing CO₂ front, where dissolved CO₂ is at or near geochemical saturation. The low‐k fault zone exerts the greatest impact on the shape of the advancing CO₂ front and restricts the bulk of the dissolved and free CO₂ to the region upstream of the fault barrier. In the high‐k aquifer, the high‐k fault zone exerts a small influence on the shape of the advancing CO₂ front. We also model stacked reservoir seal pairs, and the fracture‐dominated fault acts as a vertical bypass, allowing upward movement of CO₂ into overlying strata. High‐permeability fault zones are important pathways for CO₂ to bypass unfaulted sandstone, which leads to reduce sequestration efficiency. Aquifer compartmentalization by low‐permeability fault barriers leads to improved storativity because the barriers restrict lateral CO₂ migration and maximize the volume and pressure of CO₂ that might be emplaced in each fault‐bound compartment. As much as a 3.5‐MPa pressure increase may develop in the injected reservoir in this model domain, which under certain conditions may lead to pressures close to the fracture pressure of the top seal.     

Diffuse hydrothermal methane output and evidence of methanotrophic activity within the soils at Sousaki (Greece)

Methane soil flux measurements have been made in 38 sites at the geothermal system of Sousaki (Greece) with the closed chamber method. Fluxes range from ?47.6 to 29 150 mg m^?2 day^?1, and the diffuse CH₄ output of the system has been estimated at 19 t a^?1. Contemporaneous CO₂ flux measurements showed a moderate positive correlation between CO₂ and CH₄ fluxes. Comparison of the CO₂/CH₄ soil flux ratios with the CO₂/CH₄ ratio of the gases of the main gas manifestations provided evidence for methanotrophic activity within the soil. Laboratory CH₄ consumption experiments confirmed the presence of methanotrophic microorganisms in soil samples collected at Sousaki. Consumption was generally in the range from ?4.9 to ?38.9 pmolCH₄ h^?1 g^?1 but could sometimes reach extremely high values (?33 000 pmolCH₄ h^?1 g^?1). These results are consistent with recent studies on other geothermal systems that revealed the existence of thermoacidophilic bacteria exerting methanotrophic activity in hot, acid soils, thereby reducing methane emissions to the atmosphere.     

The stress regime of the Western Canadian Sedimentary Basin

All available information relevant to in situ stress orientations and magnitudes in the Western Canadian Sedimentary Basin (WCSB) were examined to provide a better understanding of how regional stress fields may affect geothermal development. The smallest principal stress is horizontal over most of the Western Canadian Sedimentary Basin, and it varies in magnitude across the region. Horizontal stress trajectories show that S_Hmax axes are generally aligned SW–NE. A total of 1643 measurements of microfracture and minifracture closure pressures, leak‐off pressures and fracture breakdown pressures have been harnessed to map S_Hmin gradients across the basin at depths of 156–500, 500–1000, 1000–4185 and 2000–4185 m. Vertical stress magnitudes, calculated in 91 wells, showed that at constant depth, S_V increases towards the Canadian Rocky Mountains. Resultant regional stress maps show consistent trends in orientation of stress axes. As a result, predictions can be made that propagation axes of subsurface hydraulic fractures will be dominantly SW–NE, except over the Peace River Arch area, where they will trend more towards SSW–NNE. Engineered geothermal systems in the WCSB can be optimised by drilling horizontal wells parallel to S_Hmin.