Hydrothermal,multiphase convection of H2O‐NaCl fluids from ambient to magmatic temperatures: a new numerical scheme and benchmarks for code comparison

Thermohaline convection of subsurface fluids strongly influences heat and mass fluxes within the Earth's crust. The most effective hydrothermal systems develop in the vicinity of magmatic activity and can be important for geothermal energy production and ore formation. As most parts of these systems are inaccessible to direct observations, numerical simulations are necessary to understand and characterize fluid flow. Here, we present a new numerical scheme for thermohaline convection based on the control volume finite element method (CVFEM), allowing for unstructured meshes, the representation of sharp thermal and solute fronts in advection‐dominated systems and phase separation of variably miscible, compressible fluids. The model is an implementation of the Complex Systems Modelling Platform CSMP++ and includes an accurate thermodynamic representation of strongly nonlinear fluid properties of salt water for magmatic‐hydrothermal conditions (up to 1000°C, 500 MPa and 100 wt% NaCl). The method ensures that all fluid properties are taken as calculated on the respective node using a fully upstream‐weighted approach, which greatly increases the stability of the numerical scheme. We compare results from our model with two well‐established codes, HYDROTHERM and TOUGH2, by conducting benchmarks of different complexity and find good to excellent agreement in the temporal and spatial evolution of the hydrothermal systems. In a simulation with high‐temperature, high‐salinity conditions currently outside of the range of both HYDROTHERM and TOUGH2, we show the significance of the formation of a solid halite phase, which introduces heterogeneity. Results suggest that salt added by magmatic degassing is not easily vented or accommodated within the crust and can result in dynamic, complex hydrologies.     

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     

Fluid flow in shear zones: insights from the geometry and evolution of ore bodies at Renco gold mine,Zimbabwe

The geometry of mineral deposits can give insights into fluid flow in shear zones. Lode gold ore bodies at Renco Mine, in the Limpopo Belt, Zimbabwe, occur as siliceous breccias and mylonites within amphibolite facies shear zones that dip either gently or steeply. The two sets of ore bodies formed synchronously from hydrothermal fluids. The ore bodies are oblate, but have well‐defined long axes. Larger ore bodies are more oblate. High‐grade gold ore shoots have long axes that plunge down dip; this direction is perpendicular to the long axes of the low‐grade ore bodies. The centres of the high‐grade ore bodies align within the low‐grade ore bodies along strike in both gently and steeply dipping groups. The range of sizes and shapes of the ore bodies are interpreted as a growth sequence. Geometrical models are proposed for the gently and steeply dipping ore bodies, in which individual ore bodies grow with long axes plunging down dip, and merge to form larger, more oblate ore bodies. The models show that when three or more ore bodies coalesce, the long axis of the merged ore body is perpendicular to the component ore bodies, and that ore bodies in the deposit may have a range of shapes due to both growth of individual ore bodies, and their coalescence. The long axes of the high‐grade ore bodies are parallel to the shear directions of both the gently and steeply dipping dip slip shear zones, which were the directions of greatest permeability and fluid flow. The larger, lower grade bodies, which may have formed by coalescence, are elongate perpendicular to these directions.     

Downward hydrocarbon migration predicted from numerical modeling of fluid overpressure in the Paleozoic Anticosti Basin,eastern Canada

The Anticosti Basin is a large Paleozoic basin in eastern Canada where potential source and reservoir rocks have been identified but no economic hydrocarbon reservoirs have been found. Potential source rocks of the Upper Ordovician Macasty Formation overlie carbonates of the Middle Ordovician Mingan Formation, which are underlain by dolostones of the Lower Ordovician Romaine Formation. These carbonates have been subjected to dissolution and dolomitization and are potential hydrocarbon reservoirs. Numerical simulations of fluid‐overpressure development related to sediment compaction and hydrocarbon generation were carried out to investigate whether hydrocarbons generated in the Macasty Formation could migrate downward into the underlying Mingan and Romaine formations. The modeling results indicate that, in the central part of the basin, maximum fluid overpressures developed above the Macasty Formation due to rapid sedimentation. This overpressured core dissipated gradually with time, but the overpressure pattern (i.e. maximum overpressure above source rock) was maintained during the generation of oil and gas. The downward impelling force associated with fluid‐overpressure gradients in the central part of the basin was stronger than the buoyancy force for oil, whereas the buoyancy force for gas and for oil generated in the later stage of the basin is stronger than the overpressure‐related force. Based on these results, it is proposed that oil generated from the Macasty Formation in the central part of the basin first moved downward into the Mingan and Romaine formations, and then migrated laterally up‐dip toward the basin margin, whereas gas throughout the basin and oil generated in the northern part of the basin generally moved upward. Consequently, gas reservoirs are predicted to occur in the upper part of the basin, whereas oil reservoirs are more likely to be found in the strata below the source rocks. Geofluids (2010) 10 , 334–350     

Weakening processes in thrust faults: insights from the Monte Perdido thrust fault (southern Pyrenees,Spain)

The Monte Perdido thrust fault (southern Pyrenees) consists of a 6‐m‐thick interval of intensely deformed clay‐bearing rocks. The fault zone is affected by a pervasive pressure solution seam and numerous shear surfaces. Calcite extensional‐shear veins are present along the shear surfaces. The angular relationships between the two structures indicate that shear surfaces developed at a high angle (70°) to the local principal maximum stress axis σ₁. Two main stages of deformation are present. The first stage corresponds to the development of calcite shear veins by a combination of shear surface reactivation and extensional mode I rupture. The second stage of deformation corresponds to chlorite precipitation along the previously reactivated shear surfaces. The pore fluid factor λ computed for the two deformation episodes indicates high fluid pressures during the Monte Perdido thrust activity. During the first stage of deformation, the reactivation of the shear surface was facilitated by a suprahydrostatic fluid pressure with a pore fluid factor λ equal to 0.89. For the second stage, the fluid pressure remained still high (with a λ value ranging between 0.77 and 0.84) even with the presence of weak chlorite along the shear surfaces. Furthermore, evidence of hydrostatic fluid pressure during calcite cement precipitation supports that incremental shear surface reactivations are correlated with cyclic fluid pressure fluctuations consistent with a fault‐valve model.     

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.     

Factors affecting fluid flow in strike–slip fault systems: coupled deformation and fluid flow modelling with application to the western Mount Isa Inlier, Australia

Deformation and focused fluid flow within a mineralized system are critical in the genesis of hydrothermal ore deposits. Dilation and integrated fluid flux due to coupled deformation and fluid flow in simple strike–slip fault geometries were examined using finite difference analysis in three dimensions. A series of generic fault bend and fault jog geometries consistent with those seen in the western Mount Isa Inlier were modelled in order to understand how fault geometry parameters influence the dilation and integrated fluid flux. Fault dip, fault width, bend/jog angle, and length were varied, and a cross-cutting fault and contrasting rock types were included. The results demonstrate that low fault dips, the presence of contrasts in rock type, and wide faults produce highest dilation and integrated fluid flux values. Increasing fault bend lengths and angles increases dilation and integrated fluid flux, but increasing fault jog length or angle has the opposite effect. There is minimal difference between the outputs from the releasing and restraining fault bend and jog geometries. Model characteristics producing greater fluid flows and/or gradients can be used in a predictive capacity in order to focus exploration on regions with more favorable fault geometries, provided that the mineralized rocks had Mohr–Coulomb rheologies similar to the ones used in the models.     

Nari (calcrete) outcrop contribution to ancient agricultural terraces in the Southern Shephelah,Israel: insights from digital terrain analysis and a geoarchaeological field survey

A field survey revealed that Byzantine and Early Arab (ca. 5th to 8th century C.E.) agricultural systems in the semi-arid region of the Shephelah (central Israel) were similar to runoff agricultural systems in the arid region of the Negev (southern Israel). This similarity led to the hypothesis that systems in the Shephelah also function as runoff farms. This hypothesis is not trivial since runoff values in semi-arid regions are generally low due to intensive but short rainfall events, and due to the presence of sink patches that absorb runoff on slope surface. The aim of the current research is to examine whether runoff potential in a representative agricultural system in the Shephelah is sufficient for sustaining runoff farming. A geoarchaeological field survey and digital terrain analysis show that large Nari (calcrete) outcrops on the footslopes generate high runoff values that improve water potential. Hydrological simulations and calculations show that 230 mm of direct rainfall generates a water potential equivalent to 300 mm of direct rainfall. In view of these results, it is reasonable to conclude that the presence of Nari enabled runoff agricultural farming in the Shephelah region, even in drought years.