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.     

Fluid dynamics and fluid‐structural relationships in the Red Lake mine trend,Red Lake greenstone belt,Ontario, Canada

The Red Lake mine trend, a deformation zone in the Archean Red Lake greenstone belt that hosts the world‐class Campbell‐Red Lake gold deposit, is characterized by abundant foliation‐parallel iron‐carbonate ± quartz veins with banded colloform‐crustiform structures and cockade breccias overprinted by silicification and gold mineralization. There is an apparent incompatibility between the cavity‐fill structures of the veins and breccias (typically developed at shallow crustal depths) and the upper greenschist to lower amphibole facies metamorphic conditions recorded in the host rocks (indicating relatively deep environments). This, together with the development of veins along the foliation plane, represents an enigmatic problem that may be related to the interplay between fluid dynamics and stress field. We approach this problem through systematic study of fluid inclusion planes (FIPs) in the vein minerals, including the orientations of the FIPs and the pressure–temperature conditions inferred from fluid inclusion microthermometry. We find that fluid inclusions in the main stage vein minerals (pregold mineralization ankerite and quartz and syn‐ore quartz) are predominantly carbonic without a visible aqueous phase, whereas many inclusions in the postore stage contain an aqueous phase. Most FIPs are subvertical, and many are subparallel to the foliation. High fluid pressure coupled with the high wetting angles of the water‐poor, carbonic fluids may have been responsible for the abundance of brittle deformation features. The development of subvertical FIPs is interpreted to indicate episodic switching of the maximum principal compressive stress (σ₁) from subhorizontal (perpendicular to the foliation) to subvertical (parallel to the foliation) orientation. The subvertical σ₁ is favorable for the formation of foliation‐parallel veins, as fractures are preferentially opened along the foliation in such a stress regime, the origin of which may be linked to the fluid source.     

Fluid inclusion and stable isotope constraints on the genesis of the Jian copper deposit,Sanandaj–Sirjan metamorphic zone,Iran

The Jian copper deposit, located on the eastern edge of the Sanandaj–Sirjan metamorphic zone, southwest of Iran, is contained within the Surian Permo‐Triassic volcano‐sedimentary complex. Retrograde metamorphism resulted in three stages of mineralization (quartz ± sulfide veins) during exhumation of the Surian metamorphic complex (Middle Jurassic time; 159–167 Ma), and after the peak of the metamorphism (Middle to Late Triassic time; approximately 187 Ma). The early stage of mineralization (stage 1) is related to a homogeneous H₂O–CO₂ (XCO₂ > 0.1) fluid characterized by moderate salinity (<10 wt.% NaCl equivalent) at high temperature and pressure (>370°C, >3 kbar). Early quartz was followed by small amounts of disseminated fine‐grained pyrite and chalcopyrite. Most of the main‐ore‐stage (stage 2) minerals, including chalcopyrite, pyrite and minor sphalerite, pyrrhotite, and galena, precipitated from an aqueous‐carbonic fluid (8–18 wt.% NaCl equivalent) at temperatures ranging between 241 and 388°C during fluid unmixing process (CO₂ effervescence). Fluid unmixing in the primary carbonaceous fluid at pressures of 1.5–3 kbar produced a high XCO₂ (>0.05) and a low XCO₂ (<0.01) aqueous fluid in ore‐bearing quartz veins. Oxygen and hydrogen isotope compositions suggest mineralization by fluids derived from metamorphic dehydration (δ¹⁸O_fluid = +7.6 to +10.7‰ and δD = ?33.1 to ?38.5‰) during stage 2. The late stage (stage 3) is related to a distinct low salinity (1.5–8 wt.% NaCl equivalent) and temperatures of (120–230°C) aqueous fluid at pressures below 1.5 kbar and the deposition of post‐ore barren quartz veins. These fluids probably derived from meteoric waters, which circulated through the metamorphic pile at sufficiently high temperatures and acquire the characteristics of metamorphic fluids (δ¹⁸O_fluid = +4.7 to +5.1‰ and δD = ?52.3 to ?53.9‰) during waning stages of the postearly Cimmerian orogeny in Surian complex. The sulfide‐bearing quartz veins are interpreted as a small‐scale example of redistribution of mineral deposits by metamorphic fluids. This study suggests that mineralization at the Jian deposit is metamorphogenic in style, probably related to a deep‐seated mesothermal system.     

Geological setting, nature of ore fluids and sulphur isotope geochemistry of the Fu Ning Carlin-type gold deposits, Yunnan Province, China

Carlin‐type gold deposits in southern China are present in Palaeozoic to Mesozoic siliciclastic and carbonate rocks. The border region of Yunnan, Guizhou and Guangxi Provinces contains gold deposits on the south‐western margin of the Pre‐Cambrian South China Craton in south‐eastern Yunnan Province. The Fu Ning gold deposits host epigenetic, micron‐sized disseminated gold in: (i) Middle Devonian (D₁p) black carbonaceous mudstone at the Kuzhubao gold deposit and (ii) fault breccia zones at the contact between Triassic gabbro (β ) and the Devonian mudstone (D₁p) at the Bashishan gold deposit. The deposits are associated with zones of intense deformation with enhanced permeability and porosity that focused hydrothermal fluid flow, especially where low‐angle N‐S striking thrust faults are cut by NW striking strike‐slip and/or NE striking normal faults. Major sulphide ore minerals in the Fu Ning gold deposits are pyrite, arsenopyrite, arsenic‐rich pyrite, stibnite and minor iron‐poor sphalerite. Gangue minerals are quartz, sericite, calcite, ankerite and chlorite. Hypogene ore grades range from 1 to 7 g t^?1 Au and up to 18 g t^?1 Au at the Kuzhubao gold deposit and are generally less than 3 g t^?1 Au at the Bashishan gold deposit. Sub‐microscopic gold mineralization is associated with finely disseminated arsenic‐rich pyrite in the Stage III mineral assemblage. Two types of primary fluid inclusions have been recorded: Type I liquid–vapour inclusions with moderate‐to‐high liquid/vapour ratios, and Type II inclusions containing moderate liquid/vapour ratios with CO₂ as determined from laser Raman analysis. Temperature of homogenization (T_h) data collected from these primary fluid inclusions in gold‐ore Stage III quartz ranged from 180 to 275°C at the Kuzhubao gold deposit and 210 to 330°C at the Bashishan gold deposit. Salinity results indicate that there were possibly two fluids present during gold deposition, including: (i) an early fluid with 0.8–6.5 wt.% NaCl equivalent, similar to salinity in shear‐zone‐hosted gold deposits with metamorphic derived fluids; and (ii) a late fluid with 11.8–13.4 wt.% NaCl equivalent, indicating possible derivation from connate waters and/or brine sources. CO₂ and trace CH₄ were only detected by laser Raman spectrometry in gold‐ore‐stage primary fluid inclusions. Results of sulphur isotope studies showed that δ³⁴S values for pyrite and arsenopyrite associated with gold‐ore mineralization during Stage III at the Kuzhubao and Bashishan gold deposits are isotopically similar and moderately heavy with a range from +9 to +15 per mil, and also fall into the range of δ³⁴S values reported for Carlin‐type gold deposits. Sulphur isotopes suggest that the Fu Ning gold deposits were formed from connate waters and/or basinal brines. Fluid geochemistry data from the Fu Ning gold deposits suggest a Carlin‐type genetic model, involving fluid mixing between: (i) deep CO₂‐rich metamorphic fluids, (ii) moderately saline, reduced connate waters and/or basinal brines; and (iii) evolved meteoric waters.     

Effects of episodic fluid flow on hydrocarbon migration in the Newport‐Inglewood Fault Zone,Southern California

Fault permeability may vary through time due to tectonic deformations, transients in pore pressure and effective stress, and mineralization associated with water‐rock reactions. Time‐varying permeability will affect subsurface fluid migration rates and patterns of petroleum accumulation in densely faulted sedimentary basins such as those associated with the borderland basins of Southern California. This study explores the petroleum fluid dynamics of this migration. As a multiphase flow and petroleum migration case study on the role of faults, computational models for both episodic and continuous hydrocarbon migration are constructed to investigate large‐scale fluid flow and petroleum accumulation along a northern section of the Newport‐Inglewood fault zone in the Los Angeles basin, Southern California. The numerical code solves the governing equations for oil, water, and heat transport in heterogeneous and anisotropic geologic cross sections but neglects flow in the third dimension for practical applications. Our numerical results suggest that fault permeability and fluid pressure fluctuations are crucial factors for distributing hydrocarbon accumulations associated with fault zones, and they also play important roles in controlling the geologic timing for reservoir filling. Episodic flow appears to enhance hydrocarbon accumulation more strongly by enabling stepwise build‐up in oil saturation in adjacent sedimentary formations due to temporally high pore pressure and high permeability caused by periodic fault rupture. Under assumptions that fault permeability fluctuate within the range of 1–1000 millidarcys (10^?15–10^?12 m²) and fault pressures fluctuate within 10–80% of overpressure ratio, the estimated oil volume in the Inglewood oil field (approximately 450 million barrels oil equivalent) can be accumulated in about 24 000 years, assuming a seismically induced fluid flow event occurs every 2000 years. This episodic petroleum migration model could be more geologically important than a continuous‐flow model, when considering the observed patterns of hydrocarbons and seismically active tectonic setting of the Los Angeles basin.     

High‐temperature magmatic fluid exsolved from magma at the Duobuza porphyry copper–gold deposit,Northern Tibet

The Duobuza porphyry copper–gold deposit (proven Cu resources of 2.7 Mt, 0.94% Cu and 13 t gold, 0.21 g t^?1 Au) is located at the northern margin of the Bangong‐Nujiang suture zone separating the Qiangtang and Lhasa Terranes. The major ore‐bearing porphyry consists of granodiorite. The alteration zone extends from silicification and potassic alteration close to the porphyry stock to moderate argillic alteration and propylitization further out. Phyllic alteration is not well developed. Sericite‐quartz veins only occur locally. High‐temperature, high‐salinity fluid inclusions were observed in quartz phenocrysts and various quartz veins. These fluid inclusions are characterized by sylvite dissolution between 180 and 360°C and halite dissolution between 240 and 540°C, followed by homogenization through vapor disappearance between 620 and 960°C. Daughter minerals were identified by SEM as chalcopyrite, halite, sylvite, rutile, K–feldspar, and Fe–Mn‐chloride. They indicate that the fluid is rich in ore‐forming elements and of high oxidation state. The fluid belongs to a complex hydrothermal system containing H₂O – NaCl – KCl ± FeCl₂ ± CaCl₂ ± MnCl₂. With decreasing homogenization temperature, the fluid salinity tends to increase from 34 to 82 wt% NaCl equiv., possibly suggesting a pressure or Cl/H₂O increase in the original magma. No coexisting vapor‐rich fluid inclusions with similar homogenization temperatures were found, so the brines are interpreted to have formed by direct exsolution from magma rather than trough boiling off of a low‐salinity vapor. Estimated minimum pressure of 160 MPa imply approximately 7‐km depth. This indicates that the deposit represents an orthomagmatic end member of the porphyry copper deposit continuum. Two key factors are proposed for the fluid evolution responsible for the large size of the gold‐rich porphyry copper deposit of Duobuza: (i) ore‐forming fluids separated early from the magma, and (ii) the hydrothermal fluid system was of magmatic origin and highly oxidized.     

Three-dimensional geometry, ore distribution and time-integrated mass transfer through the quartz–tourmaline–gold vein network of the Sigma deposit (Abitibi belt, Canada)

We present a reconstruction of the three‐dimensional (3D) geometry and gold grade distribution of shear zone‐hosted, Au‐mineralized, quartz–tourmaline veins of the Sigma deposit (Abitibi belt). Host shears and veins form a network of anastomosing, steeply dipping structures associated with smaller subhorizontal extensional veins. Our reconstruction has been carried out using the exceptionally large geological database of the mine. From this database, we extracted the geometric position, thickness and gold grades of geometrically best‐defined steep veins contained in a representative subvolume of the deposit. These data allowed the 3D representation of 53 veins, which have been constructed by fitting surfaces through the geometrical data and by contouring thickness and gold grade. The geometry of the network is mainly characterized by: (i) a few large segmented veins, with sinuous and helicoidal shape, and typical vertical dimension of >100 m; (ii) a large number of smaller vertical veins, some of which splay off the steep veins with high dip angles; (iii) subhorizontal extension veins (joints) located at, or close to, the tips of steep veins. The absolute thickness of the vertically short veins is the same as that of the large veins, suggesting that they formed simultaneously, but only a few of them interconnect to form vertically continuous bodies. Patchy, vertically elongated zones of high dilation are present in the large veins, and are poorly correlated with Au‐rich zones. They presumably represent former high‐permeability zones of the network. The highest gold grades occur at the interconnections between the large veins and small splays or subhorizontal joints. This indicates the important role of vein interconnection for fluid flow and gold precipitation within the network. Combining the calculation of the volume of the network with the estimation of tourmaline abundance in the veins, we calculate that 2.1 × 10⁶ m³ of tourmaline and 3.2 × 10⁶ m³ of quartz precipitated during Au deposition.     

Relationship between deformation bands and petroleum migration in an exhumed reservoir rock,Los Angeles Basin,California, USA

An oil‐bearing sandstone unit within the Monterey Formation is exposed in the Los Angeles Basin along the Newport‐Inglewood fault zone in southern California. The unit preserves structures, some original fluids, and cements that record the local history of deformation, fluid flow, and cementation. The structures include two types of deformation bands, which are cut by later bitumen veins and sandstone dikes. The bands formed by dilation and by shear. Both types strike on average parallel to the Newport‐Inglewood fault zone (317°–332°) and show variable dip angles and directions. Generally the older deformation bands are shallow, and the younger bands are steep. The earlier set includes a type of deformation band not previously described in other field examples. These are thin, planar zones of oil 1–2 mm thick sandwiched between parallel, carbonate‐cemented, positively weathering ribs. All other deformation bands appear to be oil‐free. The undeformed sandstone matrix also contains some hydrocarbons. The oil‐cored bands formed largely in opening mode, similar to dilation bands. The oil‐cored bands differ from previously described dilation bands in the degree of carbonate cementation (up to 36% by volume) and in that some exhibit evidence for plane‐parallel shear during formation. Given the mostly oil‐free bands and oil‐rich matrix, deformation bands must have formed largely before the bulk of petroleum migration and acted as semi‐permeable baffles. Oil‐cored bands provide field evidence for early migration of oil into a potential reservoir rock. We infer a hydrofracture mechanism, probably from petroleum leaking out of a stratigraphically lower overpressured reservoir. The deformation bands described here provide a potential field example of a mechanism inferred for petroleum migration in modern systems such as in the Gulf of Mexico.     

Carbon dioxide controlled earthquake distribution pattern in the NW Bohemian swarm earthquake region,western Eger Rift,Czech Republic – gas migration in the crystalline basement

The ascent of magmatic carbon dioxide in the western Eger (Oh?e) Rift is interlinked with the fault systems of the Variscian basement. In the Cheb Basin, the minimum CO₂ flux is about 160 m³ h^?1, with a diminishing trend towards the north and ceasing in the main epicentral area of the Northwest Bohemian swarm earthquakes. The ascending CO₂ forms Ca‐Mg‐HCO₃ type waters by leaching of cations from the fault planes and creates clay minerals, such as kaolinite, as alteration products on affected fault planes. These mineral reactions result in fault weakness and in hydraulically interconnected fault network. This leads to a decrease in the friction coefficient of the Coulomb failure stress (CFS) and to fault creep as stress build‐up cannot occur in the weak segments. At the transition zone in the north of the Cheb Basin, between areas of weak, fluid conductive faults and areas of locked faults with frictional strength, fluid pressure can increase resulting in stress build‐up. This can trigger strike‐slip swarm earthquakes. Fault creep or movements in weak segments may support a stress build‐up in the transition area by transmitting fluid pressure pulses. Additionally to fluid‐driven triggering models, it is important to consider that fluids ascending along faults are CO₂‐supersaturated thus intensifying the effect of fluid flow. The enforced flow of CO₂‐supersaturated fluids in the transitional zone from high to low permeability segments through narrowings triggers gas exsolution and may generate pressure fluctuations. Phase separation starts according to the phase behaviour of CO₂‐H₂O systems in the seismically active depths of NW Bohemia and may explain the vertical distribution of the seismicity. Changes in the size of the fluid transport channels in the fault systems caused, or superimposed, by fault movements, can produce fluid pressure increases or pulses, which are the precondition for triggering fluid‐induced swarm earthquakes.     

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.     

Evolution of porosity, permeability and fluid pressure in dilatant faults post-failure: implications for fluid flow and mineralization

Mineralization of brittle fault zones is associated with sudden dilation, and the corresponding changes in porosity, permeability and fluid pressure, that occur during fault slip events. The resulting fluid pressure gradients cause fluid to flow into and along the fault until it is sealed. The volume of fluid that can pass through the deforming region depends on the degree of dilation, the porosity and permeability of the fault and wall rocks, and the rate of fault sealing. A numerical model representing a steep fault cutting through a horizontal seal is used to investigate patterns of fluid flow following a dilatant fault slip event. The model is initialized with porosity, permeability and fluid pressure representing the static mechanical state of the system immediately after such an event. Fault sealing is represented by a specified evolution of porosity, coupled to changes in permeability and fluid pressure, with the rate of porosity reduction being constrained by independent estimates of the rate of fault sealing by pressure solution. The general pattern of fluid flow predicted by the model is of initial flow into the fault from all directions, followed by upward flow driven by overpressure beneath the seal. The integrated fluid flux through the fault after a single failure event is insufficient to account for observed mineralization in faults; mineralization would require multiple fault slip events. Downward flow is predicted if the wall rocks below the seal are less permeable than those above. This phenomenon could at least partially explain the occurrence of uranium deposits in reactivated basement faults that cross an unconformity between relatively impermeable basement and overlying sedimentary rocks.     

Rates of retrograde metamorphism and their implications for crustal rheology

The rate of reaction of a natural hornblende garnet granulite with water under a range of mid- to lower crustal conditions has been investigated experimentally. In runs of between 7 and 84 days small but measurable amounts of water were consumed, and sheet silicates (300°C, 300 MPa and 400°C, 400 MPa) and/or secondary actinolite (400°C, 400 MPa and 500°C, 500 MPa) were observed to have grown. When normalized to the surface area of the starting materials, hydration rates were in the range of 2–5 × 10⁻⁸ g m⁻² sec⁻¹. These reaction rates imply that a film of water that infiltrated a planar crack with a half width of 100 μm would be completely consumed within c. 100 years. These results imply that where water infiltrates the crust along faults or underlying shear zones in response to a deformation, it will remain as a free phase for only a finite period of time, which in some cases will be less than the repeat time for major earthquakes in the fault system. Thus, the rheology of fault zones and shear zones is likely to be cyclical, with the zone becoming stronger with time as water is consumed, and then weakened by infiltration of water after each rupture.     

Fluid evolution along multistage composite fault systems at the southern margin of the Lower Palaeozoic Anglo-Brabant fold belt, Belgium

Mineralised vein systems have been investigated at nine localities at the southern margin of the Anglo‐Brabant fold belt in Belgium. During the late Silurian to early Middle Devonian Caledonian orogeny, shear zones formed, inferred to be associated with granitoid basement blocks in the subsurface. The circulation of a metamorphic fluid, possibly originating in the Cambrian core of the fold belt, along these shear zones resulted in the formation of mesozonal orogenic mineralisation at the southern margin of the Anglo‐Brabant fold belt. The fluid had a composition dominated by H₂O–CO₂–X–NaCl–KCl. The shear zones form part of a greater fault zone, the Nieuwpoort–Asquempont fault zone, which is characterised by normal faulting that started before the Givetian and by the reactivation of the shear zones. Two fluid generations are associated with this normal faulting. First, a low salinity H₂O–NaCl(–KCl) fluid migrated through the Palaeozoic rocks after the Silurian. Based on the isotopic composition, this fluid could be a late‐metamorphic Caledonian fluid or a younger fluid that originated from the Rhenohercynian basin and interacted with Lower Devonian rocks along its migration path. Second, a high salinity H₂O–NaCl–CaCl₂ fluid was identified in the fault systems. Similar fluids have been found in southern and eastern Belgium, where they produced Mississippi Valley‐type Zn–Pb deposits. These fluids are interpreted as evaporative brines that infiltrated the Lower Palaeozoic basement, from where they were expelled during extensional tectonism in the Mesozoic.     

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.     

A comprehensive review of hydrocarbons and genetic model of the sandstone‐hosted Dongsheng uranium deposit,Ordos Basin,China

The Dongsheng uranium deposit, the largest in situ leach uranium mine in the Ordos Basin, geometrically forms a roll‐front type deposit that is hosted in the Middle Jurassic Zhiluo Formation. The genesis of the mineralization, however, has long been a topic of great debate. Regional faults, epigenetic alterations in surface outcrops, natural oil seeps, and experimental findings support a reducing microenvironment during ore genesis. The bulk of the mineralization is coffinite. Based on thin‐section petrography, some of the coffinite is intimately intergrown with authigenic pyrite (ore‐stage pyrite) and is commonly juxtaposed with some late diagenetic sparry calcite (ore‐stage calcite) in primary pores, suggesting simultaneous precipitation. Measured homogenization temperatures of greater than 100°C from fluid inclusions indicate circulation of low‐temperature hydrothermal fluids in the ore zone. The carbon isotopic compositions of late calcite cement (δ¹³C_VPDB = ?31.0 to ?1.4‰) suggest that they were partly derived from sedimentary organic carbon, possibly from deep‐seated petroleum fluids emanating from nearby faults. Hydrogen and oxygen isotope data from kaolinite cement (δD = ?133 to ?116‰ and δ¹⁸O_SMOW = 12.6–13.8‰) indicate that the mineralizing fluids differed from magmatic and metamorphic fluids and were more depleted in D (²H) than modern regional meteoric waters. Such a strongly negative hydrogen isotopic signature suggests that there has been selective modification of δD by CH₄±H₂S±H₂ fluids. Ore‐stage pyrite lies within a very wide range of δ³⁴S (?39.2 to 26.9‰), suggesting that the pyrite has a complex origin and that bacterially mediated sulfate reduction cannot be precluded. Hydrocarbon migration and its role in uranium reduction and precipitation have here been unequivocally defined. Thus, a unifying model for uranium mineralization can be established: Early coupled bacterial uranium mineralization and hydrocarbon oxidation were followed by later recrystallization of ore phases in association with low‐temperature hydrothermal solutions under hydrocarbon‐induced reducing conditions.     

The application of failure mode diagrams for exploring the roles of fluid pressure and stress states in controlling styles of fracture-controlled permeability enhancement in faults and shear zones

Permeability enhancement associated with deformation processes in faults and shear zones plays a key role in facilitating fluid redistribution between fluid reservoirs in the crust. Especially in high fluid flux hydrothermal systems, fracture-controlled permeability can be relatively short-lived, unless it is repeatedly regenerated by ongoing deformation. Failure mode diagrams in pore fluid factor and differential stress space, here termed λ–σ failure mode diagrams, provide a powerful tool for analysing how fluid pressure and stress states drive failure, associated permeability enhancement and vein styles during deformation in faults and shear zones. During fault-valve behaviour in the seismogenic regime, relative rates of recovery of pore fluid factor, differential stress and fault cohesive strength between rupture events impact on styles of veining and associated, fracture-controlled permeability enhancement in faults and shear zones. Examples of vein-rich fault zones are used to illustrate how constraints can be placed, not just on fluid pressure and stress states at failure, but also on the fluid pressurization and loading paths associated with failure and transitory permeability enhancement in faults and shear zones. This provides insights about when, during the fault-valve cycle, various types of veins can form. The use of failure mode diagrams also provides insights about the relative roles of optimally oriented faults and misoriented faults as hydraulically conductive structures. The analysis highlights the dynamics of competition between fluid pressures and loading rates in driving failure and repeated permeability regeneration in fracture-controlled, hydrothermal systems.     

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.     

Numerical model of pore‐pressure diffusion associated with the initiation of the 2010–2011 Guy–Greenbrier,Arkansas earthquakes

We model pore‐pressure diffusion caused by pressurized waste‐fluid injection at two nearby wells and then compare the buildup of pressure with the observed initiation and migration of earthquakes during the early part of the 2010–2011 Guy–Greenbrier earthquake swarm. Pore‐pressure diffusion is calculated using MODFLOW 2005 that allows the actual injection histories (volume/day) at the two wells to diffuse through a fractured and faulted 3D aquifer system representing the eastern Arkoma basin. The aquifer system is calibrated using the observed water‐level recovery following well shut‐in at three wells. We estimate that the hydraulic conductivities of the Boone Formation and Arbuckle Group are 2.2 × 10^?2 and 2.03 × 10^?3 m day^?1, respectively, with a hydraulic conductivity of 1.92 × 10^?2 m day^?1 in the Hunton Group when considering 1.72 × 10^?3 m day^?1 in the Chattanooga Shale. Based on the simulated pressure field, injection near the relatively conductive Enders and Guy–Greenbrier faults (that hydraulically connect the Arbuckle Group with the underlying basement) permits pressure diffusion into the crystalline basement, but the effective radius of influence is limited in depth by the vertical anisotropy of the hydraulic diffusivity. Comparing spatial/temporal changes in the simulated pore‐pressure field to the observed seismicity suggests that minimum pore‐pressure changes of approximately 0.009 and 0.035 MPa are sufficient to initiate seismic activity within the basement and sedimentary sections of the Guy–Greenbrier fault, respectively. Further, the migration of a second front of seismicity appears to follow the approximately 0.012 MPa and 0.055 MPa pore‐pressure fronts within the basement and sedimentary sections, respectively.     

Stress‐dependent transport properties of fractured arkosic sandstone

We measure the fluid transport properties of microfractures and macrofractures in low‐porosity polyphase sandstone and investigate the controls of in situ stress state on fluid flow conduits in fractured rock. For this study, the permeability and porosity of the Punchbowl Formation sandstone, a hydrothermally altered arkosic sandstone, were measured and mapped in stress space under intact, microfractured, and macrofractured deformation states. In contrast to crystalline and other sedimentary rocks, the distributed intragranular and grain‐boundary microfracturing that precedes macroscopic fracture formation has little effect on the fluid transport properties. The permeability and porosity of microfractured and intact sandstone depend strongly on mean stress and are relatively insensitive to differential stress and proximity to the frictional sliding envelope. Porosity variations occur by elastic pore closure with intergranular sliding and pore collapse caused by microfracturing along weakly cemented grain contacts. The macroscopic fractured samples are best described as a two‐component system consisting (i) a tabular fracture with a 0.5‐mm‐thick gouge zone bounded by 1 mm thick zones of concentrated transgranular and intragranular microfractures and (ii) damaged sandstone. Using bulk porosity and permeability measurements and finite element methods models, we show that the tabular fracture is at least two orders of magnitude more permeable than the host rock at mean stresses up to 90 MPa. Further, we show that the tabular fracture zone dilates as the stress state approaches the friction envelope resulting in up to a three order of magnitude increase in fracture permeability. These results indicate that the enhanced and stress‐sensitive permeability in fault damage zones and sedimentary basins composed of arkosic sandstones will be controlled by the distribution of macroscopic fractures rather than microfractures.     

Numerical modeling of the origin of calcite mineralization in the Refugio‐Carneros fault,Santa Barbara Basin,California

Many faults in active and exhumed hydrocarbon‐generating basins are characterized by thick deposits of carbonate fault cement of limited vertical and horizontal extent. Based on fluid inclusion and stable isotope characteristics, these deposits have been attributed to upward flow of formation water and hydrocarbons. The present study sought to test this hypothesis by using numerical reactive transport modeling to investigate the origin of calcite cements in the Refugio‐Carneros fault located on the northern flank of the Santa Barbara Basin of southern California. Previous research has shown this calcite to have low δ¹³C values of about ?40 to ?30‰PDB, suggesting that methane‐rich fluids ascended the fault and contributed carbon for the mineralization. Fluid inclusion homogenization temperatures of 80–125°C in the calcite indicate that the fluids also transported significant quantities of heat. Fluid inclusion salinities ranging from fresh water to seawater values and the proximity of the Refugio‐Carneros fault to a zone of groundwater recharge in the Santa Ynez Mountains suggest that calcite precipitation in the fault may have been induced by the oxidation of methane‐rich basinal fluids by infiltrating meteoric fluids descending steeply dipping sedimentary layers on the northern basin flank. This oxidation could have occurred via at least two different mixing scenarios. In the first, overpressures in the central part of the basin may have driven methane‐rich formation waters derived from the Monterey Formation northward toward the basin flanks where they mixed with meteoric water descending from the Santa Ynez Mountains and diverted upward through the Refugio‐Carneros fault. In the second scenario, methane‐rich fluids sourced from deeper Paleogene sediments would have been driven upward by overpressures generated in the fault zones because of deformation, pressure solution, and flow, and released during fault rupture, ultimately mixing with meteoric water at shallow depth. The models in the present study were designed to test this second scenario, and show that in order for the observed fluid inclusion temperatures to be reached within 200 m of the surface, moderate overpressures and high permeabilities were required in the fault zone. Sudden release of overpressure may have been triggered by earthquakes and led to transient pulses of accelerated fluid flow and heat transport along faults, most likely on the order of tens to hundreds of years in duration. While the models also showed that methane‐rich fluids ascending the Refugio‐Carneros fault could be oxidized by meteoric water traversing the Vaqueros Sandstone to form calcite, they raised doubts about whether the length of time and the number of fault pulses needed for mineralization by the fault overpressuring mechanism were too high given existing geologic constraints.