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.     

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.     

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.     

Fracture mineralization and fluid flow evolution: an example from Ordovician–Devonian carbonates,southwestern Ontario,Canada

Petrography, geochemistry (stable and radiogenic isotopes), and fluid inclusion microthermometry of matrix dolomite, fracture‐filling calcite, and saddle dolomite in Ordovician to Devonian carbonates from southwestern Ontario, Canada, provide useful insights into fluid flow evolution during diagenesis. The calculated δ¹⁸O_fluid, ΣREE, and REE_SN patterns of matrix and saddle dolomite suggest diverse fluids were involved in dolomitization and/or recrystallization of dolomite. The ⁸⁷Sr/⁸⁶Sr ratios of dolomite of each succession vary from values in the range of coeval seawater to values more radiogenic than corresponding seawater, which indicate diagenetic fluids were influenced by significant water/rock interaction. High salinities (22.4–26.3 wt. % NaCl + CaCl₂) of Silurian and Ordovician dolomite–hosted fluid inclusions indicate involvement of saline waters from dissolution of Silurian evaporites. High fluid inclusion homogenization temperatures (>100°C) in all samples from Devonian to Ordovician show temperatures higher than maximum burial (60–90°C) of their host strata and suggest involvement of hydrothermal fluids in precipitation and/or recrystallization of dolomite. A thermal anomaly over the mid‐continent rift during Devonian to Mississippian time likely was the source of excess heat in the basin. Thermal buoyancy resulting from this anomaly was the driving force for migration of hydrothermal fluids through regional aquifers from the center of the Michigan Basin toward its margin. The decreasing trend of homogenization temperatures from the basin center toward its margin further supports the interpreted migration of hydrothermal fluids from the basin center toward its margin. Hydrocarbon‐bearing fluid inclusions in late‐stage Devonian to Ordovician calcite cements with high homogenization temperatures (>80°C) and their ¹³C‐depleted values (approaching ?32‰ PDB) indicate the close relationship between hydrothermal fluids and hydrocarbon migration.     

Modern and early-middle Holocene shells of the freshwater mollusc Unio, from Çatalhöyük in the Konya Basin, Turkey: preliminary palaeoclimatic implications from molluscan isotope data

Carbon and oxygen isotope ratios in the shells of the freshwater Unio mollusc yield information on the isotopic composition of the water in which the shell was formed, which in turn relates to climatic conditions prevailing during the bivalves’ life span. Here we analysed shells from one modern Unio, from a modern lake shore in Anatolia, and 4 subfossil Unio shells from Çatalhöyük (dated between 7200 BC and 5000 BC, Neolithic and Chalcolithic periods). Sequential carbon and oxygen isotope analysis along the surface of the shells provides information on seasonal or shorter-term variability of lake waters during the lifetime of the organisms. δ¹⁸O values of the modern shell are consistent with it being collected from a dammed portion of the Seyhan River that does not undergo intensive summer evaporation. This is in contrast to many of the surface water bodies in Anatolia which suffer extensive evaporation in the arid summers. δ¹⁸O values of the subfossil shells from the Çatalhöyük middens indicate that the bivalves came from lakes which evaporated extensively during the summer months but were replenished by high winter rainfall, suggesting that the subfossil shells lived in relatively small lakes or ponds. Stable isotope analysis along the growth of freshwater bivalves is one of the few methods for investigating seasonal water fluctuations in the past. These results may alter current interpretations on the environment of Çatalhöyük during the Neolithic, although more work is needed to confirm these initial findings.     

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.     

Comparing closed system, flow-through and fluid infiltration geochemical modelling: examples from K-alteration in the Ernest Henry Fe-oxide–Cu–Au system

Potassic alteration of rocks adjacent to, and within the Ernest Henry Fe‐oxide–Cu–Au deposit is used here as a test case to investigate fluid–rock interactions using various equilibrium dynamic geochemical modelling approaches available in the HCh code. Reaction of a simple K–Fe–(Na,Ca) brine (constrained by published fluid inclusion analysis) with an albite‐bearing felsic volcanic rock, resulted in predicted assemblages defined by (i) K‐feldspar–muscovite–magnetite, (ii) biotite–K‐feldspar–magnetite, (iii) biotite–quartz–albite and (iv) albite–biotite–actinolite–pyroxene with increasing rock buffering (decreasing log w/r). Models for isothermal–isobaric conditions (450°C and 2500 bars) were compared with models run over a T–P gradient (450 to 200°C and 2500 to 500 bars). Three principal equilibrium dynamic simulation methods have been used: (i) static closed system, where individual steps are independent of all others, (ii) flow‐through and flush, where a part of the result is passed as input further along the flow line, and (iii) fluid infiltration models that simulate fluid moving through a rock column. Each type is best suited to a specific geological fluid–rock scenario, with increasing complexity, computation requirements and approximation to different parts of the natural system. Static closed system models can be used to quickly ascertain the broad alteration assemblages related to changes in the water/rock ratio, while flow‐through models are better suited to simulating outflow of reacted fluid into fresh rock. The fluid infiltration model can be used to simulate spatially controlled fluid metasomatism of rock, and we show that, given assumptions of porosity relationships and spatial dimensions, this model is a first‐order approximation to full reactive transport, without requiring significant computational time. This work presents an overview of the current state of equilibrium dynamic modelling technology using the HCh code with a view to applying these techniques to predictive modelling in exploration for mineral deposits. Application to the Ernest Henry Fe‐oxide–Cu–Au deposit demonstrates that isothermal fluid–rock reaction can account for some of the alteration zonation around the deposit.     

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.     

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.     

Moulds,Graminhos and Ribbands: a pilot study of the construction of saveiros in Valença and the Baía de Todos os Santos area,Brazil

The survival of late medieval Mediterranean techniques to conceive and build ships and boats in Brazil was noted by John Patrick Sarsfield in the 1980s, but his study of the Valença shipwrights was interrupted by his untimely death in 1990. This paper summarizes Sarsfield's account of these shipbuilding techniques, examines that published by Lev Smarcevski (1996), and provides some preliminary results of the pilot stage of a project to further research traditional shipbuilding in Valença and the Baía de Todos os Santos region.     

Complex raw materials and the supply system: Mineralogical and geochemical study of the jade artefacts of the Longshan Culture (2400–2000 bce) from Sujiacun site in coastal Shandong,China

X‐ray fluorescence (XRF), X‐ray diffraction (XRD) and Raman spectroscopy were employed to analyse the chemical compositions and phase structures of nine pieces of jade artefacts unearthed from Sujiacun, a Longshan Culture (2400–2000 bce ) site in coastal Shandong, eastern China. The results of the analyses indicated these samples were primarily made from multiple raw mineral materials, including antigorite, actinolite, clinochlore, turquoise and muscovite. No nephrite was recovered from the Sujiacun site. This is strikingly different from other top‐ranking settlements of Longshan Culture, such as Dantu, Liangchengzhen and Xizhufeng, in which nephrite was the main material of jade artefacts. This may reflect the distinctions in the supply system of jade raw materials among different rankings of Longshan period sites. In addition, laser ablation inductively coupled plasma mass spectrometry (LA‐ICP‐MS) was used to analyse the composition of the rare earth elements of serpentine jade. These preliminary results were compared with published data on the composition of serpentine jade and it was found that the Sujiacun serpentine jade artefacts were likely sourced from the Taishan jade deposit.     

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.