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.     

Fluidized sandstone intrusions as an indicator of Paleostress orientation, Santa Cruz, California

The late Miocene sandstone intrusions of northern Santa Cruz County, California, are the largest subaerial exposures of clastic intrusions on earth. The intrusions are sourced from a sandstone, underlying mudstone, accumulated in an outer shelf to upper slope environment. Dikes are the most frequent intrusion type, reach the greatest thickness and tend to strike north‐east and dip steeply. One giant dike is more than 150 m wide. Sills are least frequent, locally > 8 m thick and have no clear preferred geographical distribution. Clustered intrusions are commonly < 10 cm thick and mostly composed of dikes of various attitudes. The majority of the intrusions probably were injected shallowly as some extrude onto the seafloor. The local seafloor extrusion also indicates injection during the deposition of the Santa Cruz Mudstone (7–9 Ma). The intrusions are concentrated at the basin margin. Fluid pressure at the centre of the basin and perhaps hydrocarbons were communicated to the basin margin through the then sand, causing fluid overpressures that contributed to the fluidization and intrusion into the overlying mudstone. Primarily north‐east‐striking, steeply dipping dikes and secondarily, shallowly dipping sills are most significant in terms of regional connectivity of intrusions and physical dilation of the formation. The orientation of the dikes and sills indicates a regional stress field with a horizontal NE–SW maximum and a NW–SE minimum compressive direction. The simultaneous development of dikes and sills suggests similar magnitudes of the minimum and intermediate principal stresses. Preferential weakness along bedding contributed to the development of sills. Palaeomagnetic data indicate no significant block rotation around a vertical axis. The maximum principal stress direction indicated by the intrusions is about 55° to the San Gregorio Fault and about 70° to the San Andreas Fault during the late Miocene. This stress field is similar to the modern stress field and suggests moderate fault weakness.     

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.     

Bleached mudstone,iron concretions,and calcite veins: a natural analogue for the effects of reducing CO2‐bearing fluids on migration and mineralization of iron,sealing properties,and composition of mudstone cap rocks

This study investigates the Wangfu Depression of the Songliao Basin, China, as a natural analogue site for Fe migration (bleaching) and mineralization (formation of iron concretions) caused by reducing CO₂‐bearing fluids that leak along fractures after carbon capture, utilization, and storage. We also examined the origin of fracture‐filling calcite veins, the properties of self‐sealing fluids, the influence of fluids on the compositions of mudstone and established a bleaching model for the study area. Our results show that iron concretions are the oxidative products of precursor minerals (pyrite and siderite) during uplift and are linked to H₂S and CO₂ present in early stage fluids. The precipitation of calcite veins is the result of CO₂ degassing and is related to CO₂, CH₄, and minor heavy hydrocarbons in the main bleaching fluids. In our model, fluids preferentially enter high‐permeability fracture systems and result in the bleaching of surrounding rocks and precipitation of calcite veins. The infilling of calcite veins significantly decreases the permeability of fractures and forces the fluids to slowly enter and bleach the mudstone rocks. The Fe²⁺ released during bleaching migrates to elsewhere with the solutions or is reprecipitated in the calcite veins and iron concretions. The formation of calcite veins reduces the fracture space and effectively prevents fluid flow. The fluids have an insignificant effect on minerals within the mudstone. In terms of the chemistry of the mudstone, only the contents of Fe₂O₃, U, and Mo change significantly, with the content of U increasing in the mudstone and the contents of Fe₂O₃ and Mo decreasing during bleaching.     

Evolution of deformation and fault‐related fluid flow within an ancient methane seep system (Eocene,Varna, Bulgaria)

Faults are often important in fuelling methane seep systems; however, little is known on how different components in fault zones control subsurface fluid circulation paths and how they evolve through time. This study provides insight into fault‐related fluid flow systems that operated in the shallow subsurface of an ancient methane seep system. The Pobiti Kamani area (NE Bulgaria) encloses a well‐exposed, fault‐related seep system in unconsolidated Lower Eocene sandy deposits of the Dikilitash Formation. The Beloslav quarry and Beloslav N faults displace the Dikilitash Formation and are typified by broad, up to 80 m wide, preferentially lithified hanging wall damage zones, crosscut by deformation bands and deformation band zones, smaller slip planes and fault‐related joints. The formation of a shallow plumbing system and chimney‐like concretions in the Dikilitash Formation was followed by at least two phases of fault‐related methane fluid migration. Widespread fluid circulation through the Dikilitash sands caused massive cementation of the entire damage zones in the fault hanging walls. During this phase, paths of ascending methane fluids were locally obstructed by decimetre‐thick, continuous deformation band zones that developed in the partly lithified sands upon the onset of deformation. Once the entire damage zone was pervasively cemented, deformation proceeded through the formation of slip planes and joints. This created a new network of more localized conduits in close vicinity to the main fault plane and around through‐going slip planes. ¹³C‐depleted crustiform calcite cements in several joints record the last phase of focused methane fluid ascent. Their formation predated Neogene uplift and later meteoric water infiltration along the joint network. This illustrates how fault‐related fluid pathways evolved, over time, from ‘plumes’ in unconsolidated sediments above damage zones, leading to chimney fields, over widespread fluid paths, deflected by early deformation structures, to localized paths along fracture networks near the main fault.     

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.     

Seismic evidence for fluid migration pathways from an overpressured system in the South China Sea

Overpressured systems and intense, anomalously hot fluid expulsion in the Yinggehai Basin of the South China Sea offer an opportunity to understand the history of fluid flow and the process of hydrocarbon accumulation in overpressured environments. Fluid migration pathways from overpressured compartments in the basin are largely controlled by the distribution of faults and fractures. Episodic opening of these faults are related to the dynamics of an overpressured system and tectonic movements during basin evolution. At the crests of diapiric structures, fluid expulsion is seismically imaged as chimney‐ or plume‐like features, low to middle seismic amplitudes, and intermittently chaotic and blank reflecting seismic facies. These fluid pathways are controlled by vertical faults, which commonly penetrate overpressured and overlying normally pressured zones. Fluid expulsion is also observed near the main faults, such as the No. 1 Fault at the north‐eastern margin of the basin. Investigation by sidescan sonar on onshore and offshore Hainan Island indicates that there are more than 100 gas seepages adjacent to the No. 1 Fault. Migration pathways in the diapiric structures are controlled by three types of fault and fracture. Penetrative faults formed by dextral strike‐slip movement of the Red River faults commonly occur in the centre of the diapirs, and may have been a triggering factor for the diapirism, and controlled their distribution. Hydrofractures occur in certain mud‐rich layers and may have been generated by hydraulic fracturing. Radial normal faults occur at the top of diapirs and were formed by the intrusive process. These fluid migration pathways played an important role in regional hydrocarbon accumulation.     

Indications of formation water flow and compartmentalization on the flank of a salt structure derived from salinity and seismic data

Vertical and lateral variations in lithology, salinity, temperature, and pressure determined from wireline LAS logs, produced water samples, and seismic data on the south flank of a salt structure on the continental shelf, offshore Louisiana indicate three hydrogeologic zones in the study area: a shallow region from 0 to 1.1 km depth with hydrostatically pressured, shale‐dominated Pleistocene age sediments containing pore waters with sea water (35 g l^?1) or slightly above sea water salinity; a middle region from 1.1 to 3.2 km depth with near hydrostatically pressured, sand‐dominated Pliocene age sediments that contain pore waters that range from seawater salinity to up to 5 times sea water salinity (180 g l^?1); and a deep section below 3.2 km depth with geopressured, shale‐dominated Miocene age sediments containing pore waters that range from sea water salinity to 125 g l^?1. Salt dissolution has generated dense, saline waters that appear to be migrating down dip preferentially through the thick Pliocene sandy section. Sand layers that come in contact with salt contain pore waters with high salinity. Isolated sands have near sea water salinity. Salinity information in conjunction with seismic data is used to infer fluid compartmentalization. Both vertical and lateral lithologic barriers to fluid flow at tens to hundreds of meters scale are observed. Fluid compartmentalization is also evident across a supradomal normal fault. Offset of salinity contours are consistent with the throw of the fault, which suggests that saline fluids migrated before fault formation.     

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.     

Polyphasic hydrothermal and meteoric fluid regimes during the growth of a segmented fault involving crystalline and carbonate rocks (Barcelona Plain,NE Spain)

A polyphasic tectonic‐fluid system of a fault that involves crystalline and carbonate rocks (Hospital fault, Barcelona Plain) has been inferred from regional to thin section scale observations combined with geochemical analyses. Cathodoluminescence, microprobe analyses and stable isotopy in fracture‐related cements record the circulation of successive alternations of hydrothermal and low‐temperature meteoric fluids linked with three main regional tectonic events. The first event corresponds to the Mesozoic extension, which had two rifting stages, and it is characterized by the independent tectonic activity of two fault segments, namely southern and northern Hospital fault segments. During the Late Permian‐Middle Jurassic rifting, these segments controlled the thickness and distribution of the Triassic sediments. Also, dolomitization was produced in an early stage by Triassic seawater at shallow conditions. During increasing burial, formation of fractures and their dolomite‐related cements took place. Fault activity during the Middle Jurassic–Late Cretaceous rifting was localized in the southern segment, and it was characterized by hydrothermal brines, with temperatures over 180°C, which ascended through this fault segment precipitating quartz, chlorite, and calcite. The second event corresponds to the Paleogene compression (Chattian), which produced exhumation, folding and erosion, favouring the percolation of low‐temperature meteoric fluids which produced the calcitization of the dolostones and of the dolomite cements. The third event is linked with the Neogene extension, where three stages have been identified. During the syn‐rift stage, the southern segment of the Hospital fault grew by tip propagation. In the relay zone, hydrothermal brines with temperature around 140°C upflowed. During the late postrift, the Hospital fault acted as a unique segment and deformation occurred at shallow conditions and under a low‐temperature meteoric regime. Finally, and possibly during the Messinian compression, NW‐SE strike‐slip faults offset the Hospital fault to its current configuration.     

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.     

Hydrothermal fluid flow models of stratiform ore genesis in the McArthur Basin, Northern Territory, Australia

This paper explores the role of basin‐scale fluid migration in stratiform Pb–Zn ore formation in the southern McArthur Basin, Australia. Mathematical models are presented for coupled brine migration and heat transport in the basin. The models account for: (i) topographically driven flow (forced convection) during periods when parts of the McArthur Basin were subaerial and elevated above the central Batten Fault Zone; (ii) density‐driven flow (free convection) during periods when the basin was mostly submarine; and (iii) transient flows associated with fault rupture during periods of transpression. These hydrologic models help to compare and contrast a variety of hypotheses concerning deep fluid migration and the origin of base metal ores in the McArthur Basin. The numerical results exhibit a strong structural control on fluid flow caused by the north‐trending fault systems that characterize the Batten Fault Zone. As a result, fluids descend to depths of a few kilometers along the western side, migrate laterally to the east through the clastic and volcanic aquifers of the upper Tawallah and lowest McArthur Groups, and then ascend along the eastern side of the fault zone. This recharge–discharge pattern dominates all of the hydrogeologic models. The basin‐wide flow pattern suggests that Na–Ca–Cl brines acquired base metals in the deepest levels of the basin stratigraphy as the fluids migrated eastwards through the aquifer system. Upward flow was relatively rapid along the Emu Fault Zone, so much so that fluid temperatures likely approached 130°C in the muddy sediments near the sea floor due to upward flow and venting at the HYC (‘Here’s Your Chance'). Transient pulses of flow characterized periods of transpressional stress and subsequent faulting may have punctuated the basin history. Large‐scale free convection, however, characterized notably long periods of diagenesis and ore mineralization during the Proterozoic in the McArthur Basin.     

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.     

Changes in chemical composition caused by water–rock interactions across a strike‐slip fault zone: case study of the Atera Fault,Central Japan

In some cases, water–rock interactions in fault zones can affect radionuclide migration. Here, we analyzed the chemical compositions of well‐exposed fault rocks from the strike‐slip Atera Fault, Central Japan, in order to understand the variability and behavior of major and selected trace elements. The fault zone has a 1.2‐m‐wide, smectite‐rich fault core and paired damage zones that developed within welded tuff on one side of the core and within granite on the other side. The 30‐cm‐wide, kaolinite‐rich fault gouge is developed in granite cataclasite, and it shows indications of the latest fault activity, while the 1.2‐m‐wide fault core appears to be older. Hydrogen and oxygen isotope ratios in the clay‐rich fault gouges, and carbon and oxygen isotope ratios in carbonates indicate that the two major clay‐rich zones formed in bedrock near the surface, consistent with observed deformation structures. Based on chemical analyses, we identified (1) a slight depletion in SiO₂, Na₂O, K₂O, and light rare earth elements at the edges of the 1.2‐m‐wide fault core, (2) a clear depletion in SiO₂, Na₂O, K₂O, and all rare earth elements except Eu in the 30‐cm‐wide fault gouge, and (3) an increase in CaO, MnO, and heavy rare earth elements across the entire 1.2‐m‐wide fault core. Findings (1) and (2) reflect water–rock interactions in the 1.2‐m‐wide fault core and in the 30‐cm‐wide fault gouge that resulted in the formation of smectite and kaolinite. Finding (3) reflects carbonate precipitation caused by the addition of basalt fragments from a nearby site to the 1.2‐m‐wide fault core during faulting, and subsequent sorption reactions of heavy rare earth elements via processes such as complexation with the carbonates.     

Simulation of the impact of faults on CO2 injection into sandstone reservoirs

Numerical simulations of multiphase CO₂ behavior within faulted sandstone reservoirs examine the impact of fractures and faults on CO₂ migration in potential subsurface injection systems. In southeastern Utah, some natural CO₂ reservoirs are breached and CO₂‐charged water flows to the surface along permeable damage zones adjacent to faults; in other sites, faulted sandstones form barriers to flow and large CO₂‐filled reservoirs result. These end‐members serve as the guides for our modeling, both at sites where nature offers ‘successful’ storage and at sites where leakage has occurred. We consider two end‐member fault types: low‐permeability faults dominated by deformation‐band networks and high‐permeability faults dominated by fracture networks in damage zones adjacent to clay‐rich gouge. Equivalent permeability (k) values for the fault zones can range from <10^?14 m² for deformation‐band‐dominated faults to >10^?12 m² for fracture‐dominated faults regardless of the permeability of unfaulted sandstone. Water–CO₂ fluid‐flow simulations model the injection of CO₂ into high‐k sandstone (5 × 10^?13 m²) with low‐k (5 × 10^?17 m²) or high‐k (5 × 10^?12 m²) fault zones that correspond to deformation‐band‐ or fracture‐dominated faults, respectively. After 500 days, CO₂ rises to produce an inverted cone of free and dissolved CO₂ that spreads laterally away from the injection well. Free CO₂ fills no more than 41% of the pore space behind the advancing CO₂ front, where dissolved CO₂ is at or near geochemical saturation. The low‐k fault zone exerts the greatest impact on the shape of the advancing CO₂ front and restricts the bulk of the dissolved and free CO₂ to the region upstream of the fault barrier. In the high‐k aquifer, the high‐k fault zone exerts a small influence on the shape of the advancing CO₂ front. We also model stacked reservoir seal pairs, and the fracture‐dominated fault acts as a vertical bypass, allowing upward movement of CO₂ into overlying strata. High‐permeability fault zones are important pathways for CO₂ to bypass unfaulted sandstone, which leads to reduce sequestration efficiency. Aquifer compartmentalization by low‐permeability fault barriers leads to improved storativity because the barriers restrict lateral CO₂ migration and maximize the volume and pressure of CO₂ that might be emplaced in each fault‐bound compartment. As much as a 3.5‐MPa pressure increase may develop in the injected reservoir in this model domain, which under certain conditions may lead to pressures close to the fracture pressure of the top seal.     

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.     

Physical properties of carbonate fault rocks,fucino basin (Central Italy): implications for fault seal in platform carbonates

We documented the porosity, permeability, pore geometry, pore type, textural anisotropy, and capillary pressure of carbonate rock samples collected along basin‐bounding normal faults in central Italy. The study samples consist of one Mesozoic platform carbonate host rock with low porosity and permeability, four fractured host rocks of the damage zones, and four fault rocks of the fault cores. The four fractured samples have high secondary porosity, due to elongated, connected, soft pores that provide fluid pathways in the damage zone. We modeled this zone as an elastic cracked medium, and used the Budiansky–O'Connell correlation to compute its permeability from the measured elastic moduli. This correlation can be applied only to fractured rocks with large secondary porosity and high‐aspect ratio pores. The four fault rock samples are made up of survivor clasts embedded in fine carbonate matrices and cements with sub‐spherical, stiff pores. The low porosity and permeability of these rocks, and their high values of capillary pressure, are consistent with the fault core sealing as much as 77 and 140 m of gas and oil columns, respectively. We modeled the fault core as a granular medium, and used the Kozeny–Carmen correlation, assigning the value of 5 to the Kozeny constant, to compute its permeability from the measured porosities and pore radii. The permeability structure of the normal faults is composed of two main units with unique hydraulic characteristics: a granular fault core that acts as a seal to cross‐fault fluid flow, and an elastic cracked damage zone that surrounds the core and forms a conduit for fluid flow. Transient pathways for along‐fault fluid flow may form in the fault core during seismic faulting due to the formation of opening‐mode fractures within the cemented fault rocks.     

Migration of Mn-rich fluids through normal faults and fine-grained terrigenous sediments during early development of the Neogene Vallès-Penedès half-graben (NE Spain)

The Miocene siliciclastic sediments infilling the Vallès‐Penedès half‐graben are affected by two sets of structures developed during the extensional tectonics that created the basin. The first set, represented by extension fractures infilled with mud and sands, is attributed to seismically induced liquefaction. The second set, represented by normal faults, corresponds to a high‐permeability horsetail extensional fracture mesh developed near the surface in the hanging walls of normal faults. The incremental character of the vein‐fills indicates episodic changes in the tectonic stress state and fault zone permeability. Two episodes of fluid migration are recorded. The first episode occurred prior to consolidation and lithification when shallow burial conditions allowed oxidizing meteoric waters to flow horizontally through the more porous and permeable sandy layers. Development of clastic dikes allowed local upward flow and dewatering of the sandy beds. Liquefaction and expulsion of fluids were probably driven by seismic shaking. During the first episode of fluid migration there was no cementation of the sandstone or within the fractures, probably because little fluid was mobilized by the predominantly compaction‐driven flow regime. The second episode of fluid migration occurred synchronously with normal fault development, during which time the faults acted as fluid conduits. Fluids enriched in manganese, probably leached from local manganese oxyhydroxides soon after sedimentation, moved laterally and produced cementation in the sandstone layers, eventually arriving at the more porous and permeable fault pathways that connected compartments of different porosities and permeabilities. Carbonate probably precipitated in fractures saturated with meteoric water near the ground surface at a transitional redox potential. Once the faults became occluded by calcite cement, shortly after fault development, they became barriers to both vertical and horizontal fluid flow.     

Basin‐scale fluid movement patterns revealed by veins: Wessex Basin,UK

Calcite veins at outcrop in the Mesozoic, oil‐bearing Wessex Basin, UK, have been studied using field characterization, petrography, fluid inclusions and stable isotopes to help address the extent, timing and spatial and stratigraphic variability of basin‐scale fluid flow. The absence of quartz shows that veins formed at low temperature without an influence of hydrothermal fluids. Carbon isotopes suggest that the majority of vein calcite was derived locally from the host rock but up to one quarter of the carbon in the vein calcite came from CO₂ from petroleum source rocks. Veins become progressively enriched in source‐rock‐derived CO₂ from the outer margin towards the middle, indicating a growing influence of external CO₂. The carbon isotope data suggest large‐scale migration of substantial amounts of CO₂ around the whole basin. Fluid inclusion salinity data and interpreted water‐δ¹⁸O data show that meteoric water penetrated deep into the western part of the basin after interacting with halite‐rich evaporites in the Triassic section before entering fractured Lower and Middle Jurassic rocks. This large‐scale meteoric invasion of the basin probably happened during early Cenozoic uplift. A similar approach was used to reveal that, in the eastern part of the basin close to the area that underwent most uplift, uppermost Jurassic and Cretaceous rocks underwent vein formation in the presence of marine connate water suggesting a closed system. Stratigraphically underlying Upper Jurassic mudstone and Lower Cretaceous sandstone, in the most uplifted part of the basin, contain veins that resulted from intermediate behaviour with input from saline meteoric water and marine connate waters. Thus, while source‐rock‐derived CO₂ seems to have permeated the entire section, water movement has been more restricted. Oil‐filled inclusions in vein calcite have been found within dominant E‐W trending normal faults, suggesting that these may have facilitated oil migration.     

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.