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.     

Thermal convection in faulted extensional sedimentary basins: theoretical results from finite-element modeling

Thermal convection has the potential to be a significant and widespread mechanism of fluid flow, mass transport, and heat transport in rift and other extensional basins. Based on numerical simulation results, large‐scale convection can occur on the scale of the basin thickness, depending on the Rayleigh number for the basin. Our analysis indicates that for syn‐rift and early post‐rift settings with a basin thickness of 5 km, thermal convection can occur for basal heat flows ranging from 80 to 150 mW m^?2, when the vertical hydraulic conductivity is on the order of 1.5 m year^?1 and lower. The convection cells have characteristic wavelengths and flow patterns depending on the thermal and hydraulic boundary conditions. Steeply dipping extensional faults can provide pathways for vertical fluid flow across large thicknesses of basin sediments and can modify the dynamics of thermal convection. The presence of faults perturbs the thermal convective flow pattern and can constrain the size and locations of convection cells. Depending on the spacing of the faults and the hydraulic properties of the faults and basin sediments, the convection cells can be spatially organized to align with adjacent faults. A fault‐bounded cell occurs when one convection cell is constrained to occupy a fault block so that the up‐flow zone converges into one fault zone and the down‐flow zone is centred on the adjacent fault. A fault‐bounded cell pair occurs when two convection cells occupy a fault block with the up‐flow zone located between the faults and the down‐flow zones centred on the adjacent faults or with the reverse pattern of flow. Fault‐bounded cells and cell pairs can be referred to collectively as fault‐bounded convective flow. The flow paths in fault‐bounded convective flow can be lengthened significantly with respect to those of convection cells unperturbed by the presence of faults. The cell pattern and sense of circulation depend on the fault spacing, sediment and fault permeabilities, lithologic heterogeneity, and the basal heat flow. The presence of fault zones also extends the range of conditions for which thermal convection can occur to basin settings with Rayleigh numbers below the critical value for large‐scale convection to occur in a basin without faults. The widespread potential for the occurrence of thermal convection suggests that it may play a role in controlling geological processes in rift basins including the acquisition and deposition of metals by basin fluids, the distribution of diagenetic processes, the temperature field and heat flow, petroleum generation and migration, and the geochemical evolution of basin fluids. Fault‐bounded cells and cell pairs can focus mass and heat transport from longer flow paths into fault zones, and their discharge zones are a particularly favourable setting for the formation of sediment‐hosted ore deposits near the sea floor.     

The application of structured‐light illumination microscopy to hydrocarbon‐bearing fluid inclusions

Structured‐light illumination (SLI)‐based microscopy offers geologists a new perspective for screening of hydrocarbon‐bearing (HCFI) and small aqueous fluid inclusion (AFI) assemblages. This optical‐sectioning technique provides rapid, confocal‐like imaging, using relatively simple and inexpensive instrumentation. The 3D fluorescent images of HCFI planes and large single HCFIs permit the visualization of the relationships between HCFI assemblages, examination of HCFI morphology, and volume estimates of the fluorescent components within HCFIs. By the use of normal white light illumination, SLI image capture, and varying acquisition time it is also possible to image AFI because of the random movements of vapour bubbles within the inclusions. This allows the near‐simultaneous visualization of hydrocarbon and AFI which is of significant importance for the study of sedimentary basins and petroleum reservoirs. SLI is a unique and accessible 3D petrographic tool, with practical advantages over conventional epifluorescence and confocal laser scanning microscopy.     

The significance of pockmarks to understanding fluid flow processes and geohazards

Underwater gas and liquid escape from the seafloor has long been treated as a mere curiosity. It was only after the advent of the side‐scan sonar and the subsequent discovery of pockmarks that the scale of fluid escape and the moon‐like terrain on parts of the ocean floor became generally known. Today, pockmarks ranging in size from the ‘unit pockmark’ (1–10 m wide, < 0.6 m deep) to the normal pockmark (10–700 m wide, up to 45 m deep) are known to occur in most seas, oceans, lakes and in many diverse geological settings. In addition to indicating areas of the seabed that are ‘hydraulically active’, pockmarks are known to occur on continental slopes with gas hydrates and in association with slides and slumps. However, possibly their potentially greatest significance is as an indicator of deep fluid pressure build‐up prior to earthquakes. Whereas only a few locations containing active (bubbling) pockmarks are known, those that become active a few days prior to major earthquakes may be important precursors that have been overlooked. Pockmark fields and individual pockmarks need to be instrumented with temperature and pressure sensors, and monitoring should continue over years. The scale of such research calls for a multinational project in several pockmark fields in various geological settings.     

Fluid flow and stability of the US continental slope offshore New Jersey from the Pleistocene to the present

We predict that portions of the New Jersey continental slope were unstable approximately 0.5 million years ago. This instability was caused by rapid sediment loading during a Pleistocene sea‐level lowstand and by flow focusing in underlying, permeable Miocene strata. The simulated instability is consistent with soft‐sediment deformation and small slumps in Pleistocene strata of the Hudson Apron. Stability of the New Jersey margin has increased since 0.3 Ma because sedimentation rate has decreased. Today, the modelled factor of safety (FS) for the upper slope is approximately 1.5 whereas in the lower slope it exceeds 3. We predict that sedimentation rate is a dominant factor on slope stability. When rapid and asymmetric loading of a highly permeable sedimentary layer occurs, the location of instability can shift seaward to regions where sedimentation rates are low. Stability calculations use pressures and effective stresses predicted by a coupled sedimentation‐fluid flow model. This hydrodynamic analysis demonstrates how the interplay of sedimentation and fluid migration affects the distribution, timing, and size of sedimentary failures.     

Pore pressure evolution and fluid flow during visco‐elastic single‐layer buckle folding

Pore pressure and fluid flow during the deformational history of geologic structures are directly influenced by tectonic deformation events. In this contribution, 2D plane strain finite element analysis is used to study the influence of different permeability distributions on the pore pressure field and associated flow regimes during the evolution of visco‐elastic single‐layer buckle folds. The buckling‐induced fluid flow regimes indicate that flow directions and, to a lesser degree, their magnitudes vary significantly throughout the deformation and as a function of the stratigraphic permeability distribution. The modelling results suggest that the volumetric strain and the permeability distribution significantly affect the resulting flow regime at different stages of fold development. For homogeneous permeability models (k > 10^?21 m²), low strain results in a mostly pervasive fluid flow regime and is in agreement with previous studies. For larger strain conditions, fluid focusing occurs in the buckling layer towards the top of the fold hinge. For low permeabilities (<10^?21 m²), local focused flow regimes inside the buckling layer emerge throughout the deformation history. For models featuring a low‐permeability layer embedded in a high‐permeability matrix or sandwiched between high‐permeability layers, focused flow regimes inside the folded layer result throughout the deformation history, but with significant differences in the flow vectors of the surrounding layers. Fluid flow vectors induced by the fold can result in different, even reversed, directions depending on the amount of strain. In summary, fluid flow regimes during single‐layer buckling can change from pervasive to focused and fluid flow vectors can be opposite at different strain levels, that is the flow vectors change significantly through time. Thus, a complete understanding of fluid flow regimes associated with single‐layer buckle folds requires consideration of the complete deformation history of the fold.     

Long‐term chemical evolution and modification of continental basement brines – a field study from the Schwarzwald,SW Germany

Highly saline, deep‐seated basement brines are of major importance for ore‐forming processes, but their genesis is controversial. Based on studies of fluid inclusions from hydrothermal veins of various ages, we reconstruct the temporal evolution of continental basement fluids from the Variscan Schwarzwald (Germany). During the Carboniferous (vein type i), quartz–tourmaline veins precipitated from low‐salinity (<4.5wt% NaCl + CaCl₂), high‐temperature (≤390°C) H₂O‐NaCl‐(CO₂‐CH₄) fluids with Cl/Br mass ratios = 50–146. In the Permian (vein type ii), cooling of H₂O‐NaCl‐(KCl‐CaCl₂) metamorphic fluids (T ≤ 310°C, 2–4.5wt% NaCl + CaCl₂, Cl/Br mass ratios = 90) leads to the precipitation of quartz‐Sb‐Au veins. Around the Triassic–Jurassic boundary (vein type iii), quartz–haematite veins formed from two distinct fluids: a low‐salinity fluid (similar to (ii)) and a high‐salinity fluid (T = 100–320°C, >20wt% NaCl + CaCl₂, Cl/Br mass ratios = 60–110). Both fluids types were present during vein formation but did not mix with each other (because of hydrogeological reasons). Jurassic–Cretaceous veins (vein type iv) record fluid mixing between an older bittern brine (Cl/Br mass ratios ~80) and a younger halite dissolution brine (Cl/Br mass ratios >1000) of similar salinity, resulting in a mixed H₂O‐NaCl‐CaCl₂ brine (50–140°C, 23–26wt% NaCl + CaCl₂, Cl/Br mass ratios = 80–520). During post‐Cretaceous times (vein type v), the opening of the Upper Rhine Graben and the concomitant juxtaposition of various aquifers, which enabled mixing of high‐ and low‐salinity fluids and resulted in vein formation (multicomponent fluid H₂O‐NaCl‐CaCl₂‐(SO₄‐HCO₃), 70–190°C, 5–25wt% NaCl‐CaCl₂ and Cl/Br mass ratios = 2–140). The first occurrence of highly saline brines is recorded in veins that formed shortly after deposition of halite in the Muschelkalk Ocean above the basement, suggesting an external source of the brine's salinity. Hence, today's brines in the European basement probably developed from inherited evaporitic bittern brines. These were afterwards extensively modified by fluid–rock interaction on their migration paths through the crystalline basement and later by mixing with younger meteoric fluids and halite dissolution brines.     

Freezing and melting behaviors of H2O‐NaCl‐CaCl2 solutions in fused silica capillaries and glass‐sandwiched films: implications for fluid inclusion studies

Fluid inclusions of the H₂O‐NaCl‐CaCl₂ system are notorious for their metastable behavior during cooling and heating processes, which can render microthermometric measurement impossible or difficult and interpretation of the results ambiguous. This study addresses these problems through detailed microscopic examination of synthetic solutions during cooling and warming runs, development of methods to enhance nucleation of hydrates, and comparison of microthermometric results with different degrees of metastability with values predicted for stable conditions. Synthetic H₂O‐NaCl‐CaCl₂ solutions with different NaCl/(NaCl + CaCl₂) ratios were prepared and loaded in fused silica capillaries and glass‐sandwiched films for microthermometric studies; pure solutions were used with the capillaries to simulate fluid inclusions, whereas alumina powder was added in the solutions to facilitate ice and hydrate crystallization in the sandwiched samples. The phase changes observed and the microthermometric data obtained in this study have led to the following conclusions that have important implications for fluid inclusion studies: (i) most H₂O‐NaCl‐CaCl₂ inclusions that appear to be completely frozen in the first cooling run to ?185°C actually contain large amounts of residual solution, as also reported in some previous studies; (ii) inability of H₂O‐NaCl‐CaCl₂ inclusions to freeze completely may be related to their composition (low NaCl/(NaCl + CaCl₂) ratios) and lack of solid particles; (iii) crystallization of hydrates, which is important for cryogenic Raman spectroscopic studies of fluid inclusion composition, can be greatly enhanced by finding an optimum combination of cooling and warming rates and temperatures; and (iv) even if an inclusion is not completely frozen, the melting temperatures of hydrohalite and ice are still valid for estimating the fluid composition.     

Culturally modified skeletal remains from the site of Huamelulpan,Oaxaca, Mexico

A series of crania from the site of Huamelulpan, Oaxaca, Mexico (400 BC to AD 800), were examined. Four showed notable cultural modifications. One exhibited a healed trephination, while the other three were perforated through the frontal. The cultural context and significance of these modifications is discussed, especially in relationship to the site of Monte Albán, where trephination was more common than anywhere else in Mesoamerica. The post-mortem cranial perforations appear to be connected with the practice of ancestor veneration. © 1997 John Wiley & Sons, Ltd.     

Infiltration of basinal fluids into high-grade basement, South Norway: sources and behaviour of waters and brines

Quartz veins hosted by the high‐grade crystalline rocks of the Modum complex, Southern Norway, formed when basinal fluids from an overlying Palaeozoic foreland basin infiltrated the basement at temperatures of c. 220°C (higher in the southernmost part of the area). This infiltration resulted in the formation of veins containing both two‐phase and halite‐bearing aqueous fluid inclusions, sometimes with bitumen and hydrocarbon inclusions. Microthermometric results demonstrate a very wide range of salinities of aqueous fluids preserved in these veins, ranging from c. 0 to 40 wt% NaCl equivalent. The range in homogenization temperatures is also very large (99–322°C for the entire dataset) and shows little or no correlation with salinity. A combination of aqueous fluid microthermometry, halogen geochemistry and oxygen isotope studies suggest that fluids from a range of separate aquifers were responsible for the quartz growth, but all have chemistries comparable to sedimentary formation waters. The bulk of the quartz grew from relatively low δ¹⁸O fluids derived directly from the basin or equilibrated in the upper part of the basement (T < 200°C). Nevertheless, some fluids acquired higher salinities due to deep wall‐rock hydration reactions leading to salt saturation at high temperatures (>300°C). The range in fluid inclusion homogenization temperatures and densities, combined with estimates of the ambient temperature of the basement rocks suggests that at different times veins acted as conduits for influx of both hotter and colder fluids, as well as experiencing fluctuations in fluid pressure. This is interpreted to reflect episodic flow linked to seismicity, with hotter dry basement rocks acting as a sink for cooler fluids from the overlying basin, while detailed flow paths reflected local effects of opening and closing of individual fractures as well as reaction with wall rocks. Thermal considerations suggest that the duration of some flow events was very short, possibly in the order of days. As a result of the complex pattern of fracturing and flow in the Modum basement, it was possible for shallow fluids to penetrate basement rocks at significantly higher temperatures, and this demonstrates the potential for hydrolytic weakening of continental crust by sedimentary fluids.