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.     

Computing geochemical mass transfer and water/rock ratios in submarine hydrothermal systems: implications for estimating the vigour of convection

A method of calculating chemical water/rock ratios is presented that enables the estimation of fluid velocities in open, flow‐through hydrologic systems. The approach is based on relating the gain/loss of a chemical species per kilogram of solid phase to the loss/gain of that species in the fluid phase, integrated across a specified length of the flowpath. After examining the underlying approximations of the approach using a one‐dimensional model of seawater moving through a basalt under nonisothermal conditions, the method is applied to representative zones within a two‐dimensional hydrothermal convective system. The method requires that regions within the flow system can be identified in which the direction of flow is steady for an extended period of time. Estimates of fluid velocity are spatial and temporal averages for the length of the flowpath used in the calculation. The location within the flow system and the nature of the alteration reactions determine which species can provide reliable values of the chemical water/rock ratio and useful estimates of fluid velocities. Over the length of the flowpath considered, the calculation of water/rock ratios works best when a species is controlled by a single reaction. Accurate estimates are favoured if the concentration profile of a species along the flowpath increases or decreases monotonically. If the length of the flowpath extends over more than one reaction zone, then erroneous estimates of the water/rock ratio and fluid velocity are more likely. Model calculations suggest that the quartz/silica system should provide reliable estimates for fluid velocity under a wide range of temperature and flow conditions, in particular in those regions of a system at or near quartz equilibrium, so that the aqueous silica concentration is buffered by quartz and correlated with the temperature distribution.     

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     

Hydrothermal convection in and around mineralized fault zones: insights from two- and three-dimensional numerical modeling applied to the Ashanti belt, Ghana

Among hydrogeological processes, free convection in faults has been cited as a possible cause of gold mineralization along major fault zones. Here, we investigate the effects of free convection to determine whether it can cause giant orogenic gold deposits and their regular spatial distribution along major fault/shear zones. The approach comprises: (i) coupled two- and three-dimensional numerical heat- and fluid-flow simulations of simplified geological models; and (ii) calculation of the rock alteration index (RAI) to delineate regions where precipitation/dissolution can occur. Then, comparing the deduced alteration patterns with temperature distribution, potential areas of gold mineralization, defined by T  > 200°C and RAI < 0, are predicted. The models are based on the orogenic Paleoproterozoic ore deposits of the Ashanti belt in western Africa. These deposits occur in the most permeable parts of the fault zone, where the lateral permeability contrast is the highest. For a simple geometry, with a fault zone adjacent to a sedimentary basin half as permeable, we note a transition from three-dimensional circulation within the fault to a two-dimensional convective pattern in the basin far from the fault. Moreover, whereas two-dimensional undulated isotherms dominate in the basin, three-dimensional corrugated isotherms result from the preferred convective pattern within the fault, thus enhancing a periodic distribution of thermal highs and lows. In our most elaborate three-dimensional model with an imposed lateral permeability gradient, the RAI distribution indicates that fluid circulation in fault zones gives rise to a spatial periodicity of alteration patterns consistent with field data.     

Heat transport at boiling,near‐critical conditions

Two‐phase flow and near‐critical phenomena are likely to enhance energy transport in high‐temperature hydrothermal systems. We present a series of two‐dimensional simulations of two‐phase flow of pure water at near‐critical conditions. The results show that at near‐critical conditions, two‐phase convection can be more efficient in transporting energy than single‐phase convection. The highest heat fluxes are attained when two‐phase heat‐pipes form near the bottom boundary, recharging the root of the upflow zone and thereby enabling the formation of broad upflow regions. When the system becomes more vapor‐dominated, it loses this ability, upflow zones become narrower and the energy efficiency drops to more moderate values.     

A vector of high‐temperature paleo‐fluid flow deduced from mass transfer across permeability barriers (quartz veins)

Quartz veins acted as impermeable barriers to regional fluid flow and not as fluid‐flow conduits in Mesoproterozoic rocks of the Mt Painter Block, South Australia. Systematically distributed asymmetric alteration selvedges consisting of a muscovite‐rich zone paired with a biotite‐rich zone are centered on quartz veins in quartz–muscovite–biotite schist. Geometric analysis of the orientation and facing of 126 veins at Nooldoonooldoona Waterhole reveals a single direction along which a maximum of all veins have a muscovite‐rich side, irrespective of their specific individual orientation. This direction represents a Mesoproterozoic fluid‐flow vector and the veins represent permeability barriers to the flow. The pale muscovite‐rich zones formed on the downstream side of the vein and the dark biotite‐rich zones mark the upstream side. The alteration couplets formed from mica schist at constant Zr, Ga, Sc, and involved increases in Si, Na, Al and decreases in K, Fe, Mg for pale alteration zones, and inverse alteration within dark zones. The asymmetry of the alteration couplets is best explained by the pressure dependence of mineral–fluid equilibria. These equilibria, in combination with a Darcian flow model for coupled advection and diffusion, and with permeability barriers imposed by the quartz veins, simulate the pattern of both fluid flow and differential, asymmetric metasomatism. The determined vector of fluid flow lies along the regional foliation and is consistent with the known distribution of regional alteration products. The presence of asymmetric alteration zones in rock containing abundant pre‐alteration veins suggests that vein‐rich material may have generally retarded regional fluid flow.     

Massive dolomitization of a Messinian reef in the Great Bahama Bank: a numerical modelling evaluation of Kohout geothermal convection

The hypothesis that Kohout thermal convection may have induced the massive dolomitization of the 60 m thick lowest more reefal unit in well Unda [top of Great Bahama Bank (GBB)] is evaluated through numerical modelling. A two‐dimensional (2‐D) section, including lithological and petrophysical data, together with datings for the sediments of the GBB, was used in the basin model TEMISPACK to reconstruct the history of the whole platform, with a focus on the reef unit. Simulations showed that during high sea‐level periods, Kohout convection is a valid mechanism in the settings of the GBB, although the convection cell remains flat in most cases because of high permeability anisotropy. This mechanism induces rapid fluid flow in the superficial as well as in the deeper parts of the platform, with velocities of at least two orders of magnitude higher than with compaction alone. Lithology appears as a strong control of fluid circulations at the margin scale through the permeability anisotropy, for which a critical value lies between values of 10 and 100. The reefal unit in Unda is part of a larger area determined by the lithologic distribution, in which flow velocities are significantly higher than in the rest of the platform. These velocities are high enough to bring the magnesium necessary to precipitate the observed amounts of dolomite, within durations in agreement with the available time of post‐reef deposition high sea level(s). However, neither fluid flow pattern nor flow velocities are able to explain the preferential massive dolomitization of the lower reef unit and the complete absence of dolomite in the upper one.     

Modeling fluid flow in a heterogeneous, fault-controlled hydrothermal system

Previous studies have shown that most hydrothermal systems discharging at the land surface are associated with faulting, and that the location, temperature and rate of discharge of these systems are controlled by the geometry and style of the controlling fault(s). Unfortunately, the transport of heat and fluid in fault-controlled hydrothermal systems is difficult to model realistically; although heterogeneity and anisotropy are assumed to place important controls on flow in faults, few data or observations are available to constrain the distribution of hydraulic properties within active faults. Here, analytical and numerical models are combined with geostatistical models of spatially varying hydraulic properties to model the flow of heat and fluid in the Borax Lake fault of south-east Oregon, USA. A geometric mean permeability within the fault of 7 × 10⁻¹⁴ m² with 2× vertical/horizontal anisotropy in correlation length scale is shown to give the closest match to field observations. Furthermore, the simulations demonstrate that continuity of flow paths is an important factor in reproducing the observed behavior. In addition to providing some insight into possible spatial distributions of hydraulic properties at the Borax Lake site, the study highlights one potential avenue for integrating field observations with simulation results in order to gain greater understanding of fluid flow in faults and fault-controlled hydrothermal and petroleum reservoirs.     

Chemical modification of pyroclastic rock by hot water: an experimental investigation of mass transport at the fluid–solid interface

Hydrothermal water–(pyroclastic) rock interactions were examined using flow-through experiments to deduce the effect of mass transport phenomena on the reaction process. A series of experiments were conducted over the temperature range 75–250°C, with a constant temperature for each experiment, and at saturated vapour pressure, to estimate the apparent rate constants as a function of temperature.
Based on the chemistry of analysed solutions, the water–rock interaction in the experiments was controlled by diffusion from the reaction surface and by the existence of a surface layer at the rock–fluid interface, which regulated the chemical reaction rate. The reaction progress depended to a high degree on flow velocity and temperature conditions, with element abundances in the fluid significantly affected by these factors. Mass transport coefficients for diffusion from the rock surface to the bulk solution have been estimated. Ca is selectively depleted under lower temperature conditions ( T  < 150°C), whereas Na is greatly depleted under higher temperature conditions ( T  > 150°C), and K reaction rates are increased when flow velocity increases. Using these conditions, specific alkali and alkali earth cations were selectively leached from mineral surfaces. The 'surface layer' comprised a 0.5–1.8 mm boundary film on the solution side (the thickness of this layer has no dependence on chemical character) and a reaction layer. The reaction layer was composed of a Si, Al-rich cation-leached layer, whose thickness was dependent on temperature, flow velocity and reaction length. The reaction layer varied in thickness from about 10⁻⁴ to 10⁻⁷ mm under high temperature/low fluid velocity and low temperature/high fluid velocity conditions, respectively.     

The potential for abiotic organic synthesis and biosynthesis at seafloor hydrothermal systems

Calculations are presented of the extent to which chemical disequilibria are generated when submarine hydrothermal fluids mix with sea water. These calculations involve quantifying the chemical affinity for individual reactions by comparing equilibrium compositions with the compositions of mixtures in which oxidation–reduction reactions are inhibited. The oxidation–reduction reactions that depart from equilibrium in these systems provide energy for chemotrophic microbial metabolism. Methanogenesis is an example of this phenomenon, in which the combination of carbon dioxide, hydrogen and methane induced by fluid mixing is far from equilibrium, which can be approached if more methane is generated. Similarly, the production of other organic compounds is also favorable under the same conditions that permit methanogenesis. Alkanes, alkenes, alcohols, aldehydes, carboxylic acids and amino acids are among the compounds that, if formed, would lower the energetic state of the chemical composition generated in mixed fluids. The resulting positive values of chemical affinity correspond to the thermodynamic drive required for abiotic organic synthesis. It is also possible that energy release accompanies biosynthesis by chemotrophic organisms. In this way, hydrothermal ecosystems differ radically from familiar ecosystems at Earth's surface. If captured, the energy released may be sufficient to drive biosynthesis of carbohydrates, purines, pyrimidines and other compounds which require energy inputs.     

Hydrothermal control on organic matter alteration and illite precipitation, Mt Isa Basin, Australia

Abundant illite precipitation in Proterozoic rocks from Northern Lawn Hill Platform, Mt Isa Basin, Australia, occurred in organic matter‐rich black shales rather than in sandstones, siltstones and organic matter‐poor shales. Sandstones and siltstones acted as impermeable rocks, as early diagenetic quartz and carbonate minerals reduced the porosity–permeability. Scanning and transmission electron microscopy (SEM and TEM) studies indicate a relation between creation of microporosity–permeability and organic matter alteration, suitable for subsequent mineral precipitation. K–Ar data indicate that organic matter alteration and the subsequent illite precipitation within the organic matter occurred during the regional hydrothermal event at 1172 ± 50 (2σ) Ma. Hot circulating fluids are considered to be responsible for organic matter alteration, migration and removal of volatile hydrocarbon, and consequently porosity–permeability creation. Those rocks lacking sufficient porosity–permeability, such as sandstones, siltstones and organic matter poor shales, may not have been affected by fluid movement. In hydrothermal systems, shales and mudstones may not be impermeable as usually assumed because of hydrocarbons being rapidly removed by fluid, even with relatively low total organic carbon.     

Factors controlling free thermal convection in faults in sedimentary basins: implications for the formation of zinc–lead mineral deposits

Stratiform sediment‐hosted Zn–Pb–Ag mineral deposits constitute about 40% of the Earth's zinc resources ( Allen 2001 ), and in most cases their genesis involves the discharge of basinal brines near or on the seafloor through syndepositional faults ( Sangster 2002 ). From the point of view of base metal exploration, it is therefore essential to identify all possible faults that formerly carried the upwelling ore‐forming solutions during mineralising events. This paper presents a numerical investigation of the relative importance of various physical parameters in controlling fluid discharge, recharge and heat transport in faults. A two‐dimensional, free convection of pure water, hydrogeological model is developed for the McArthur basin in northern Australia based on the surface geology, known stratigraphic and structural relationships and regional geophysical interpretations. Numerical experiments and sensitivity analyses reveal that faults with strong initial heat input, due to depth of penetration or magmatic activity, are the most likely candidates to carry discharge fluids to the sites of metal precipitation. Deeper, wider and more permeable faults are more likely to behave as the fluid discharge pathways, whereas shallow, narrow or less permeable faults act as marine water recharge pathways. Compared with these fault‐related factors, aquifer physical properties are less important in determining fluid flow patterns and the geothermal regime. These results are an important step in understanding hydrothermal fluid flow in sedimentary basins in order to develop effective exploration criteria for the location of stratiform Zn–Pb–Ag deposits.     

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.     

3D numerical modeling of hydrothermal processes during the lifetime of a deep geothermal reservoir

Understanding hydrothermal processes during production is critical to optimal geothermal reservoir management and sustainable utilization. This study addresses the hydrothermal (HT) processes in a geothermal research doublet consisting of the injection well E GrSk3/90 and production well Gt GrSk4/05 at the deep geothermal reservoir of Groß Schönebeck (north of Berlin, Germany) during geothermal power production. The reservoir is located between ?4050 to ?4250 m depth in the Lower Permian of the Northeast German Basin. Operational activities such as hydraulic stimulation, production (T = 150°C; Q = ?75 m³ h^?1; C = 265 g l^?1) and injection (T = 70°C; Q = 75 m³ h^?1; C = 265 g l^?1) change the HT conditions of the geothermal reservoir. The most significant changes affect temperature, mass concentration and pore pressure. These changes influence fluid density and viscosity as well as rock properties such as porosity, permeability, thermal conductivity and heat capacity. In addition, the geometry and hydraulic properties of hydraulically induced fractures vary during the lifetime of the reservoir. A three‐dimensional reservoir model was developed based on a structural geological model to simulate and understand the complex interaction of such processes. This model includes a full HT coupling of various petrophysical parameters. Specifically, temperature‐dependent thermal conductivity and heat capacity as well as the pressure‐, temperature‐ and mass concentration‐dependent fluid density and viscosity are considered. These parameters were determined by laboratory and field experiments. The effective pressure dependence of matrix permeability is less than 2.3% at our reservoir conditions and therefore can be neglected. The results of a three‐dimensional thermohaline finite‐element simulation of the life cycle performance of this geothermal well doublet indicate the beginning of thermal breakthrough after 3.6 years of utilization. This result is crucial for optimizing reservoir management. Geofluids (2010) 10 , 406–421     

Granite‐related overpressure and volatile release in the mid crust: fluidized breccias from the Cloncurry District,Australia

The source and transport regions of fluidized (transported) breccias outcrop in the Cloncurry Fe‐oxide–Cu–Au district. Discordant dykes and pipes with rounded clasts of metasedimentary calc–silicate rocks and minor felsic and mafic intrusions extend several kilometres upwards and outwards from the contact aureole of the 1530 Ma Williams Batholith into overlying schists and amphibolites. We used analytical equations for particle transport to estimate clast velocities (≥20 m sec^?1), approaching volcanic ejecta rates. An abrupt release of overpressured magmatic‐hydrothermal fluid is suggested by the localization of the base of the breccias in intensely veined contact aureoles (at around 10 km, constrained by mineral equilibria), incorporation of juvenile magmatic clasts, the scale and discordancy of the bodies, and the wide range of pressure variation (up to 150 MPa) inferred from CO₂ fluid inclusion densities and related decrepitation textures. The abundance of clasts derived from depth, rather than from the adjacent wallrocks, suggests that the pressure in the pipes was sufficient to restrict the inwards spalling of fragments from breccia walls; that is, the breccias were explosive rather than implosive, and some may have vented to the surface. At these depths, such extreme behaviour may have been achieved by release of dissolved fluids from crystallizing magma, in combination with a strongly fractured and fluid‐laden carapace, sitting under a strong, low permeability barrier. The relationship of these breccias to the Ernest Henry iron‐oxide–Cu–Au deposit suggests they may have been sources of fluids or mechanical energy for ore genesis, or alternately provided permeable pathways for later ore fluids.