Permeability of continental crust influenced by internal and external forcing

The permeability of continental crust is so highly variable that it is often considered to defy systematic characterization. However, despite this variability, some order has been gleaned from globally compiled data. What accounts for the apparent coherence of mean permeability in the continental crust (and permeability–depth relations) on a very large scale? Here we argue that large‐scale crustal permeability adjusts to accommodate rates of internal and external forcing. In the deeper crust, internal forcing – fluxes induced by metamorphism, magmatism, and mantle degassing – is dominant, whereas in the shallow crust, external forcing – the vigor of the hydrologic cycle – is a primary control. Crustal petrologists have long recognized the likelihood of a causal relation between fluid flux and permeability in the deep, ductile crust, where fluid pressures are typically near‐lithostatic. It is less obvious that such a relation should pertain in the relatively cool, brittle upper crust, where near‐hydrostatic fluid pressures are the norm. We use first‐order calculations and numerical modeling to explore the hypothesis that upper‐crustal permeability is influenced by the magnitude of external fluid sources, much as lower‐crustal permeability is influenced by the magnitude of internal fluid sources. We compare model‐generated permeability structures with various observations of crustal permeability.     

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 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.     

Exploring the potential linkages between oil‐field brine reinjection,crystalline basement permeability,and triggered seismicity for the Dagger Draw Oil field,southeastern New Mexico,USA, using hydrologic modeling

We used hydrologic models to explore the potential linkages between oil‐field brine reinjection and increases in earthquake frequency (up to M_d 3.26) in southeastern New Mexico and to assess different injection management scenarios aimed at reducing the risk of triggered seismicity. Our analysis focuses on saline water reinjection into the basal Ellenburger Group beneath the Dagger Draw Oil field, Permian Basin. Increased seismic frequency (>M_d 2) began in 2001, 5 years after peak injection, at an average depth of 11 km within the basement 15 km to the west of the reinjection wells. We considered several scenarios including assigning an effective or bulk permeability value to the crystalline basement, including a conductive fault zone surrounded by tighter crystalline basement rocks, and allowing permeability to decay with depth. We initially adopted a 7 m (0.07 MPa) head increase as the threshold for triggered seismicity. Only two scenarios produced excess heads of 7m five years after peak injection. In the first, a hydraulic diffusivity of 0.1 m² s^?1 was assigned to the crystalline basement. In the second, a hydraulic diffusivity of 0.3 m² s^?1 was assigned to a conductive fault zone. If we had considered a wider range of threshold excess heads to be between 1 and 60 m, then the range of acceptable hydraulic diffusivities would have increased (between 0.1–0.01 m² s^?1 and 1–0.1 m² s^?1 for the bulk and fault zone scenarios, respectively). A permeability–depth decay model would have also satisfied the 5‐year time lag criterion. We also tested several injection management scenarios including redistributing injection volumes between various wells and lowering the total volume of injected fluids. Scenarios that reduced computed excess heads by over 50% within the crystalline basement resulted from reducing the total volume of reinjected fluids by a factor of 2 or more.     

Spring temperatures in the Sagehen Basin, Sierra Nevada, CA: implications for heat flow and groundwater circulation

Heat flow in the Sierra Nevada, CA, is low despite its young geologic age. We investigate the possibility that advective heat transport by groundwater flow leads to an underestimate of heat flow in the Sierras based purely on borehole measurements. Using temperature and discharge measurements at springs in Sagehen Basin, we find that groundwater removes the equivalent of approximately 20–40 mW m⁻² of geothermal heat from the basin. This is comparable with other heat flow measurements in the region and indicates that, in this basin, at least, groundwater does transport a significant amount of geothermal heat within the basin. Additionally, we use estimates of the mean residence time of water discharged at the springs along with hourly temperature records in springs to provide constraints on groundwater flow depths within the basin. An analytical model based on these constraints indicates that the heat removed by groundwater may represent 20% to >90% of the total heat flow in the basin. Without better constraints on the regional hydrogeology and the depth of circulation, we cannot determine whether the heat discharged at the springs represents a change in the mode of heat transfer, i.e. from conduction to advection at shallow depths (<100 m) or whether this is a component of heat transfer that should be added to measured conductive values. If the latter is true, and Sagehen Basin is representative of the Sierras, basal heat flow in the Sierra Nevada may be higher than previously thought.     

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.     

Modeling fluid flow in a heterogeneous, fault-controlled hydrothermal system

Previous studies have shown that most hydrothermal systems discharging at the land surface are associated with faulting, and that the location, temperature and rate of discharge of these systems are controlled by the geometry and style of the controlling fault(s). Unfortunately, the transport of heat and fluid in fault-controlled hydrothermal systems is difficult to model realistically; although heterogeneity and anisotropy are assumed to place important controls on flow in faults, few data or observations are available to constrain the distribution of hydraulic properties within active faults. Here, analytical and numerical models are combined with geostatistical models of spatially varying hydraulic properties to model the flow of heat and fluid in the Borax Lake fault of south-east Oregon, USA. A geometric mean permeability within the fault of 7 × 10⁻¹⁴ m² with 2× vertical/horizontal anisotropy in correlation length scale is shown to give the closest match to field observations. Furthermore, the simulations demonstrate that continuity of flow paths is an important factor in reproducing the observed behavior. In addition to providing some insight into possible spatial distributions of hydraulic properties at the Borax Lake site, the study highlights one potential avenue for integrating field observations with simulation results in order to gain greater understanding of fluid flow in faults and fault-controlled hydrothermal and petroleum reservoirs.     

AN ANGLO-SAXON DISC BROOCH FROM BLETCHINGDON, OXFORDSHIRE, WITH STYLE II ANIMAL ART

Summary.   The brooch was found in 2003 near Bletchingdon, Oxfordshire. ¹ Its gilt Style II animal ornament and garnet setting indicate that it is a high status dress-fastener, datable to the period AD 575–610. Analysis of its form and ornament suggests that it was a local experiment in the application of a new animal style.     

Large‐scale chemical stratification of fluids in the crust: hydraulic and chemical data from the geothermal research site Urach,Germany

The Urach 3 research borehole in SW Germany has been drilled through a sedimentary cover sequence and reached gneisses of the Variscan crystalline basement at 1604 m below surface. An additional 2840 m has been drilled through fractured basement rocks. The borehole has been used for hydraulic tests in the context of a ‘hot dry rock’ (HDR) project. The sedimentary cover ranges from the Carboniferous to the Middle Jurassic (Dogger) in age and comprises mostly clastic sediments in the Paleozoic and limestone and shale in the Mesozoic. Water composition data from 10 different depths include samples from all major lithological units. The total dissolved solids (TDS) increases from the surface to about 650 m where it reaches 4.1 g l^?1 in Triassic limestone. In lower Triassic sandstones, TDS increases very sharply to 28.5 g l^?1 and the water is saturated with pure CO₂ gas. With increasing depth, TDS does not change much in the clastic sediments of the Permian and Carboniferous. The crystalline basement is marked by a very sharp increase in TDS to 55.5 g l^?1 at about 1770 m depth. TDS increases within the basement to more than 78.5 g l^?1 at about 3500 m depth. The data suggest that there is limited vertical chemical communication over long periods of time. The CO₂ gas cap in the lower Triassic sandstones requires a gastight cover. The chemical stratification of the fluids relates to the permeability structure of the crust at the Urach site and fits well with hydraulic and thermal data from the site.     

Source of ore fluids at the Kalahari Goldridge deposit, Kraaipan greenstone belt, South Africa: evidence from Sr, C and O isotope signatures in carbonates

The Kalahari Goldridge deposit is located in the Archaean Kraaipan greenstone belt in the north-west province of South Africa. Gold mineralization in this deposit is hosted within banded iron formation which is flanked by a mafic schist in the footwall and clastic metasedimentary units in the hanging wall. Data from carbonate minerals from mineralized veins and bulk rock from the A and D zone ore bodies have helped to define the ultimate origin of the ore-forming fluids and their migration history. Carbon isotope ratios of carbonates from both the A and D zone ore bodies have tight clustering from −7.6 to −5.3‰ that indicates a unique origin for the ore-forming fluids associated with the mineralization at Kalahari Goldridge. The δ¹⁸O values of the carbonates have been influenced by temperature gradients and variable degrees of fluid–rock interaction promoting oxygen isotope exchange between ore fluid and host rocks. Minimum ⁸⁷Sr/⁸⁶Sr ratio values of 0.70354 in mineralized veins are most consistent with ore-forming fluids being relatively pristine with a mantle origin. Strontium and the corresponding ore-forming fluids were most likely derived from mantle-derived magmatic rocks probably represented by the meta-basaltic rocks that underlie the ferruginous package in the Kraaipan greenstone belt. Strontium isotopic composition of vein carbonates show considerable variation in ⁸⁷Sr/⁸⁶Sr ratios ranging from 0.70354 to 0.73914. This is consistent with an ore fluid composition that has been modified by the addition of radiogenic Sr possibly during passage of fluid through siliciclastic country rock concomitant with the observed hydrothermal alteration.     

Upper crustal fluid flow in the outboard region of the Southern Alps,New Zealand

The currently active fluid regime within the outboard region of the Southern Alps, New Zealand was investigated using a combination of field observations, carbon‐ and oxygen‐stable isotopes from fault‐hosted calcites and interpretation of magnetotelluric (MT) data. Active faulting in the region is dominated by NE striking and N striking, oppositely dipping thrust fault pairs. Stable isotopic analyses of calcites hosted within these fault zones range from 10 to 25‰δ¹⁸O and from ?33 to 0‰δ¹³C. These values reflect mixing of three parent fluids: meteoric water, carbon‐exchanged groundwater and minor deeper rock‐exchanged fluids, at temperatures of 10–90°C in the upper 3–4 km of the crust. A broad, ‘U‐shaped’ zone of high electrical conductivity (maximum depth c. 28 km) underlies the central Southern Alps. In the ductile region of the crust, the high‐conductivity zone is subhorizontal. Near‐vertical zones of high conductivity extend upward to the surface on both sides of the conductive zone. On the outboard side of the orogen, the conductive zone reaches the surface coincident with the trace of the active Forest Creek Faults. Near‐surface flow is shown to dominate the outboard region. Topographically driven meteoric water interacts, on a kilometre scale, with either carbon‐exchanged groundwater or directly with organic material within Pliocene gravels, resulting in a distinctive low ¹³C signal within fault‐hosted calcites of the outboard region. The high‐strain zone in the lower crust focuses the migration of deeply sourced fluids upward to the base of the brittle–ductile transition. Interconnected fluid is imaged as a narrow vertical zone of high conductivity in the upper crust, implying continuous permeability and possibly buoyancy‐driven flow of deeply sourced fluids to higher levels of the crust where they are detected by the isotopic analysis of the fault‐hosted calcites.     

Structural controls on hydrothermal flow in a segmented rift system, Taupo Volcanic Zone, New Zealand

Vigorous hydrothermal convection transfers 10 times the average continental heat flow through the central Taupo Volcanic Zone (TVZ), a region of active extension (approximately 8 mm year^?1) and productive rhyolitic volcanism. Over 20 high‐temperature (>250°C) geothermal fields occur within Quaternary pyroclastic basins, with convective circulation to depths of 7–8 km presumably extending through basement rocks. Parallel‐striking normal faults, fractures and dikes dissect the convective regime, interacting with fluids to either enhance or restrict flow according to the relative permeability of structure and host rock. In the basement, high bulk permeability is maintained by focussed flow through faults and associated fractures well oriented for reactivation in the prevailing stress field. In contrast, distributed flow through fault‐bounded compartments prevails within Quaternary basins, masking any signal of deeper structural control. Exceptions occur where more competent rocks are exposed at the surface. As in narrow magmatic rifts elsewhere, the extensional fabric is partitioned into discrete rift segments linked along strike by accommodation zones. Eighty per cent of TVZ geothermal fields correlate spatially with rift architecture, with 60% located in accommodation zones. We suggest that segmented rift fabrics generate bulk permeability anisotropy that is to some extent predictable, with rift segments characterized by enhanced axial flow, and accommodation zones characterized by locally enhanced vertical permeability that is tectonically maintained. This provides a plausible explanation for the common occurrence of geothermal fields within accommodation zones and their notable absence within densely faulted rift segments. Maintenance of structural permeability in zones of active hydrothermal precipitation necessarily requires repeated brittle failure. Geothermal plumes therefore exploit tectonically maintained permeability within accommodation zones, with rift segments functioning mostly as drawdown regions. The influence of rift architecture on flow paths has important implications for geothermal extraction and epithermal mineral exploration within the TVZ and other structurally segmented hydrothermal systems, both active and extinct.     

Spatial variations in the salinity of pore waters in northern deep water Gulf of Mexico sediments: implications for pathways and mechanisms of solute transport

Spatial variations in the salinity of pore waters in sedimentary basins can provide important insight into basin-scale hydrogeologic processes. Although there have been numerous studies of brine seeps in the deep water Gulf of Mexico, much less is known about porewater salinities in the vast areas between seeps. A study has been made of spatial variation in pore water salinities in sediments in an approximately 500 km by 200 km area of the northern deep water (water depth >500 m) Gulf of Mexico sedimentary basin (GOM) to provide insight into pathways and mechanisms of solute transport in this portion of the basin. A second objective was to document salinities in the upper 500 m of the sedimentary section, the approximate depth to which methane hydrates, a potential future energy resource, may be stable. Elevated salinities would reduce the P – T stability range of hydrates. Salinities were calculated from borehole logs using a dual electrical conductivity model. Even though much of the northern GOM is underlain by allochthonous salt most of the undisturbed shallow sedimentary section has not been permeated by hypersaline waters. These waters are limited to areas near brine seeps. Hypersaline waters having salinities in excess of 100 g l⁻¹ become more common at subseafloor depths of 2 km and greater. A field study at Green Canyon 65 and published numerical simulations of fluid flow above tabular salt bodies suggest that brines produced by salt dissolution migrate laterally and pond above salt and/or within minibasins and that the dominant mechanism of vertical solute transport is a combination of compaction-driven advection and diffusion, not large-scale thermohaline overturn. Superimposed on this diffuse upward flux of dissolved salt is the more focused and localized expulsion of saline fluids up faults.     

Heterogeneity effects on possible salinity-driven free convection in low-permeability strata

Although studies of free convection commonly focus on highly permeable strata, but numerical analyses indicate that density-driven free convection may also occur in heterogeneous low-permeability strata. Traditional Rayleigh number criteria are overly conservative in predicting thermohaline convection in these systems; so, numerical models are used to make inferences on the process. Simulations with stochastic realizations of permeability fields show that dense plumes can take preferential pathways to sink through generally low-permeability strata; patch analysis using percolation theory shows that the threshold permeability for the onset of free convection can be as low as 10⁻¹⁶ m² even with a mean permeability of 10⁻¹⁸ m². Threshold permeability for the percolation pathways decreases with increasing concentration gradient, vertical correlation length and the mean and variance of the permeability. The connectedness of relatively high-permeability zones is important in initiating and controlling plume fingers of free convection in both single-layer and sand-shale sequence models. Permeable units above and below are conducive to free convection through intervening low-permeability strata if buoyancy gradients exist. This heterogeneity is on scales that are difficult to sample by drilling and too localized to be simulated in regional models but may be significant in solute transport in these systems.     

Chemical modification of pyroclastic rock by hot water: an experimental investigation of mass transport at the fluid–solid interface

Hydrothermal water–(pyroclastic) rock interactions were examined using flow-through experiments to deduce the effect of mass transport phenomena on the reaction process. A series of experiments were conducted over the temperature range 75–250°C, with a constant temperature for each experiment, and at saturated vapour pressure, to estimate the apparent rate constants as a function of temperature.
Based on the chemistry of analysed solutions, the water–rock interaction in the experiments was controlled by diffusion from the reaction surface and by the existence of a surface layer at the rock–fluid interface, which regulated the chemical reaction rate. The reaction progress depended to a high degree on flow velocity and temperature conditions, with element abundances in the fluid significantly affected by these factors. Mass transport coefficients for diffusion from the rock surface to the bulk solution have been estimated. Ca is selectively depleted under lower temperature conditions ( T  < 150°C), whereas Na is greatly depleted under higher temperature conditions ( T  > 150°C), and K reaction rates are increased when flow velocity increases. Using these conditions, specific alkali and alkali earth cations were selectively leached from mineral surfaces. The 'surface layer' comprised a 0.5–1.8 mm boundary film on the solution side (the thickness of this layer has no dependence on chemical character) and a reaction layer. The reaction layer was composed of a Si, Al-rich cation-leached layer, whose thickness was dependent on temperature, flow velocity and reaction length. The reaction layer varied in thickness from about 10⁻⁴ to 10⁻⁷ mm under high temperature/low fluid velocity and low temperature/high fluid velocity conditions, respectively.     

Rapid Na, Cu exchange between synthetic fluid inclusions and external aqueous solutions: evidence from LA-ICP-MS analysis

Three sets of equilibration experiments (Set 1 to Set 3) were performed in cold-seal pressure vessels to investigate the compositional modification of quartz-hosted fluid inclusions after entrapment. Each set of experiment consisted of two stages. In a pre-run, fluid inclusions containing 5–10 wt% NaCl and selected trace elements were synthesized at 700°C/140–200 MPa. These samples were then loaded into new capsules together with Cu-bearing solutions and some mineral buffers, and re-equilibrated at 600–800°C and 70–140 MPa for 6–8 days. LA-ICP-MS analysis of individual fluid inclusions reveals that in re-equilibration experiments in which the outer fluid was composed of KCl (±FeCl₂) up to 83% of the original Na content of pre-existing fluid inclusions were lost, and up to 5660 ppm Cu were gained. Other elements with larger ionic radii (i.e. K, Fe, Ba, Sr) were not exchanged, demonstrating that the inclusions remained physically intact and that Na and Cu were transported through quartz by diffusion. The observed Na loss from pre-existing fluid inclusions correlates positively with Cu gain, with about 1 Cu atom being gained per 10 Na atoms lost. Thus, Na and Cu (plus probably H) were exchanged by interdiffusion. Remarkably, this processes resulted in up to 10 times higher Cu concentrations in re-equilibrated inclusions than were present in the outer fluid, i.e. Cu diffused 'uphill'. Large variations of Cu concentrations relative to the concentration of other elements are common also in natural fluid inclusion assemblages. However, no evidence for a correlation between Cu content and Na content was found so far, suggesting that Cu diffusion in natural samples may be dominated by processes other than Na–Cu interdiffusion.     

Modelling the Lost City hydrothermal field: influence of topography and permeability structure

The Lost City hydrothermal field (LCHF) is hosted in serpentinite at the crest of the Atlantis Massif, an oceanic core complex close to the mid‐Atlantic Ridge. It is remarkable for its longevity and for venting low‐temperature (40–91°C) alkaline fluids rich in hydrogen and methane. IODP Hole U1309D, 5 km north of the LCHF, penetrated 1415 m of gabbroic rocks and contains a near‐conductive thermal gradient close to 100°C km^?1. This is remarkable so close to an active hydrothermal field. We present hydrothermal modelling using a topographic profile through the vent field and IODP site U1309. Long‐lived circulation with vent temperatures similar to the LCHF can be sustained at moderate permeabilities of 10^?14 to 10^?15 m² with a basal heatflow of 0.22 W m^?2. Seafloor topography is an important control, with vents tending to form and remain in higher topography. Models with a uniform permeability throughout the Massif cannot simultaneously maintain circulation at the LCHF and the near‐conductive gradient in the borehole, where permeabilities <10^?16 m² are required. A steeply dipping permeability discontinuity between the LCHF and the drill hole is required to stabilize venting at the summit of the massif by creating a lateral conductive boundary layer. The discontinuity needs to be close to the vent site, supporting previous inferences that high permeability is most likely produced by faulting related to the transform fault. Rapid increases in modelled fluid temperatures with depth beneath the vent agree with previous estimates of reaction temperature based on geochemical modelling.     

Spatial coupling between spilitization and carbonation of basaltic sills in SW Scottish Highlands: evidence of a mineralogical control of metamorphic fluid flow

In a geochemical and petrological analysis of overprinting episodes of fluid–rock interaction in a well‐studied metabasaltic sill in the SW Scottish Highlands, we show that syn‐deformational access of metamorphic fluids and consequent fluid–rock interaction is at least in part controlled by preexisting mineralogical variations. Lithological and structural channelling of metamorphic fluids along the axis of the Ardrishaig Anticline, SW Scottish Highlands, caused carbonation of metabasaltic sills hosted by metasedimentary rocks of the Argyll Group in the Dalradian Supergroup. Analysis of chemical and mineralogical variability across a metabasaltic sill at Port Cill Maluaig shows that carbonation at greenschist to epidote–amphibolites facies conditions caused by infiltration of H₂O‐CO₂ fluids was controlled by mineralogical variations, which were present before carbonation occurred. This variability probably reflects chemical and mineralogical changes imparted on the sill during premetamorphic spilitization. Calculation of precarbonation mineral modes reveals heterogeneous spatial distributions of epidote, amphibole, chlorite and epidote. This reflects both premetamorphic spilitization and prograde greenschist facies metamorphism prior to fluid flow. Spilitization caused albitization of primary plagioclase and spatially heterogeneous growth of epidote ± calcic amphibole ± chlorite ± quartz ± calcite. Greenschist facies metamorphism caused breakdown of primary pyroxene and continued, but spatially more homogeneous, growth of amphibole + chlorite ± quartz. These processes formed diffuse epidote‐rich patches or semi‐continuous layers. These might represent precursors of epidote segregations, which are better developed elsewhere in the SW Scottish Highlands. Chemical and field analyses of epidote reveal the evidence of local volume fluctuations associated with these concentrations of epidote. Transient permeability enhancement associated with these changes may have permitted higher fluid fluxes and therefore more extensive carbonation. This deflected metamorphic fluid such that its flow direction became more layer parallel, limiting propagation of the reaction front into the sill interior.     

The influence of poorly interconnected fault zone flow paths on spring geochemistry

Thermal springs commonly occur along faults because of the enhanced vertical permeability afforded by fracture zones. Field and laboratory studies of fault zone materials document substantial heterogeneities in fracture permeabilities. Modeling and field studies of springs suggest that spatial variations in permeability strongly influence spring locations, discharge rates and temperatures. The impact of heterogeneous permeability on spring geochemistry, however, is poorly documented. We present stable isotope and water chemistry data from a series of closely spaced thermal springs associated with the Hayward Fault, California. We suggest that substantial spatial variations observed in δ¹⁸O and chloride values reflect subsurface fluid transport through a poorly connected fracture network in which mixing of subsurface waters remains limited. Our measurements provide insight into the effect of fracture zone heterogeneities on spring geochemistry, offer an additional tool to intuit the nature of tectonically induced changes in fault zone plumbing, and highlight the need to consider local variations when characterizing fracture zone fluid geochemistry from spring systems with multiple discharge sites.