Abstract

Minerals in veins and shear zones often show oxygen and hydrogen isotope ratios that are interpreted as recording interaction with meteoric water, at depths up to about 10 km. Downward fluid flow to these depths can only occur in the unlikely case of fluid pressures that are significantly lower than lithostatic overburden pressures. We therefore propose that fluid movement was upward instead of downward. In our model, the pore space within sediments and exhumed rocks below an unconformity is filled with meteoric and possibly seawater fluids. Burial of these rocks traps the fluids that can retain their meteoric isotopic composition as long as temperatures remain below about 300-350°C. Extension or rapid exhumation, such as that experienced by metamorphic core complexes, which results in decompression or fluid heating can release these old “meteoric” fluids, of which we find the isotopic fingerprint in veins and shear zone minerals.

1. Introduction

Fluids of different origins are characterised by their typical chemical and isotopic composition [1]. Meteoric fluids, for example, are characterised by very low salinity and (relatively) low δ18O and δD values. Fluid composition and isotope ratios can change by a variety of processes, such as fluid-rock interaction, mixing with other fluids, and gravitational fractionation [27]. Unfortunately, it is usually difficult to directly measure the isotopic composition of palaeofluids. Instead, minerals that precipitated from or interacted with these fluids can be used to infer their compositions, as these fluids may have left their “fingerprint” or “signature” in the isotopic composition of the minerals (e.g., [7]). If the mineral data indicate that the fluid had the composition of a meteoric fluid, the mineral can be said to carry a “meteoric signature.” The same applies to seawater, metamorphic, or magmatic signatures, among others. Similarly, a fluid has a “meteoric signature” if its composition is similar to that of meteoric water or of meteoric water that has been altered by, for example, fluid-rock interaction or mixing with other fluids (e.g., brines, metamorphic, and magmatic). The term “signature” is used here to emphasise that whether or not a palaeofluid is meteoric, or has a meteoric component or origin, is an interpretation of compositional data, normally including isotopes, and not a direct observation in case of data derived from the mineral record.

Meteoric signatures, mostly δ18O values, have been claimed to be visible in veins and shear zones that formed in deep, even ductile parts of the crust [828]. For want of other explanations for these meteoric signatures, they are usually explained as resulting from interaction with meteoric water. This raises the question of how these surface waters reached the ductile parts of the crust with °C and depths over 10 km [11, 22, 23]. Meteoric signatures are in particular found in mylonitic detachment zones of metamorphic core complexes, which are formed during rapid exhumation [1316, 20, 21, 29]. The meteoric signature in these studies is constituted by (i) lowered δ18O values in the detachment, relative to the adjacent less-deformed rocks, (ii) low (mostly <-100‰) δD values in micas, especially inside the detachment, and (3) low δD values in fluid inclusions in quartz. To explain these observations, most current authors favour downward flow of meteoric fluid along extensional faults or detachments in association with exhumation and shearing of the rocks in which the meteoric signature is found (e.g., [9, 20, 21, 3032]), even to depths of up to 20 km [17, 22]. The alternative model of Clark et al. [12] envisages that the rocks achieved their meteoric signature at near-surface conditions and were subsequently buried. The meteoric signature observed in the detachment would then be inherited and not the result of deep penetration of meteoric fluids. However, the observed isotopic shifts that authors such as Mulch et al. [21] and Gébelin et al. [29] use to infer a meteoric fluid influx into metamorphic core complex detachments are more prominent in the top of the detachments, where shearing is most intense, compared to footwall rocks further down. This would suggest that the rocks did not carry an inherited meteoric signature before the detachment activity but that the isotopic shifts occurred during shearing.

A problem with the downward flow of fluid is that deep circulation of fluids is only possible when the fluid pressure is near hydrostatic [9, 12]. Only then can fluids penetrate deeply due to topography differences or thermally induced convection, as envisioned by the aforementioned authors. Circulation of fluids from the surface to a depth of up to even 20 km and temperatures of around 400°C has indeed been modelled numerically (e.g., [17, 22, 23, 33]). To achieve such circulation and maintain it over hundreds of thousands of years or more, these authors used constant and high permeabilities of 10-15 to 10-14 m2 and typically high porosities of up to 10% [33] for the fault zones. Note that these permeabilities are high and expected for upper crustal rocks (<5 km depth; [34]) with high intergranular porosities (ca. 10%) or alternatively (micro-)fracturing or damage zones [3537]. Furthermore, a hydrostatic fluid pressure at 10 to 20 km depth implies that the difference between fluid and lithostatic pressure is about 1.7 to 3.4 kbar, unless there is a seal higher up in the system. As already pointed out by Bons et al. [38], it is highly unlikely that such a fluid underpressure can be maintained under ductile conditions at elevated temperatures.

The intense shear fabrics (mylonitic fabric, recrystallised quartz, and mica fish) of the mylonitic rocks of detachments in which evidence for meteoric fluids was found (e.g., the Columbia River Fault; [23]; or the South Tibet Detachment; [15]) are evidence of the ability of these rocks to flow ductilely. The difference in fluid and solid pressure would drive ductile flow of the rock around pores and fractures to equalise the pressure. Here, we present a simple numerical model to demonstrate that less than a million years are required to equalise 50% the fluid and solid pressure in a rock’s pore at a temperature above about 250 to 400°C, assuming a spherical pore surrounded by ductile quartz rock that flows in to equalise the pressure difference (see Methods and Results and Discussion for details). Note that (i) at 50% pressure equilibration, fluid pressure is significantly above hydrostatic pressure and hydraulic head gradients will drive fluid upward and (ii) this temperature range corresponds to the brittle-to-ductile transition (e.g., [39] and references therein).

Because of the arguments given above, it is highly unlikely, if not impossible, to drive meteoric fluids down to ductile regimes where the temperature is in the order of 300 to 400°C for any prolonged period of time. This implies that the faults and shear zones that now exhibit the meteoric water signature were not themselves involved in the downward transport of the meteoric water and that the downward transport did not occur at the same time of the recording of the meteoric signature in the rocks. Raimondo et al. [27] suggested that the fault rocks acquired their meteoric signature near the surface (where meteoric fluid flow and alterations are not controversial) and were then buried and later reactivated as midcrustal shear zones. As an alternative, we argue that the meteoric fluids were originally near-surface pore waters that were buried passively [40]. By using an oxygen isotope diffusion model, we show how such old meteoric fluids were able to keep their isotopic signatures for long time (tens to hundreds of million years) if their temperature was below about 350°C. We also discuss mechanisms for the upward release of such trapped fluids during exhumation and present fluid budgets for an idealised metamorphic core complex. The model we present here, in which fluid infiltration and release are separated in time, is a hydrodynamically feasible alternative to the classical paradigm in which downward flow and upward flow are simultaneous.

2. Methods

2.1. Model to Calculate the Rate of Fluid and Rock Pressure Equilibration

To constrain the time it would take to equalise the fluid and solid pressure, we here consider the idealised case of a perfect fluid-filled sphere surrounded by solid, ductile rock that flows in to equalise the pressure difference. A sphere is the most robust structure, and any other geometry of a pore would close faster than a sphere. Assuming, for simplicity, that the rock is an incompressible material with a rheology that can be represented by the quartz flow laws of Paterson and Luan [41] or Gleason and Tullis [42], we can employ the equations by Nye [43], derived again by Dabrowski et al. [44], to calculate the rate of pressure equilibration between fluid and rock. Or, in other words, we can estimate what the shrinking rate of a fluid-filled inclusion is.

Nye [43] derived the equations for the shrinking of an underpressured spherical inclusion for the case that is controlled by the rheology of the embedding matrix. The matrix is assumed to be an incompressible material with a power-law rheology that relates the strain intensity () to the stress ():

Here, is the stress exponent and is a rate constant. Note that flow laws for rock-forming minerals are usually expressed in a different way with instead a rate constant . The conversion will be discussed below. Strain intensity and stress are derived from the strain rate () and the deviatoric stress tensor ():

For an isolated spherical inclusion with radius , the radial stress () is now given by where is the stress at the surface of the inclusion, at distance , assuming zero radial stress at infinity. thus equals the pressure difference, , between solid and fluid in the inclusion. The rate of contraction is

Flow laws for rock-forming minerals, such as quartz, are typically derived from uniaxial shortening experiments and reported as

The subscript refers to the direction of uniaxial shortening, is the preexponential rate factor, is the activation energy, is the universal gas constant, and is the absolute temperature. To obtain the rate factor that is used above, we use and to obtain

To calculate the shrinkage rate of an inclusion, one needs to take into account the change in pressure within the inclusion. We use a compressibility of  GPa, with

Here, is the reference volume at which . The change in pressure difference between the fluid and the matrix can now be related to the matrix rheology:

The integration constant equals , with the initial pressure difference at . From this, we finally derive the time, , needed to achieve a pressure difference :

3. Results and Discussion

3.1. Fluid and Rock Pressure Equilibration at the Brittle-to-Ductile Transition

We present models based on the method presented above to estimate the rate of equilibration of a fluid-filled spherical pore within a ductile quartz rock to show that the time to equilibrate the pressure difference between a fluid at hydrostatic pressure and the solid matrix at lithostatic pressure by 50% depends strongly on the temperature and, hence, on depth (Figure 1). Using the three flow laws proposed by Paterson and Luan [41] or Gleason and Tullis [42] (see Table 1), we see that this time quickly reduces from many millions of years to less than a million year in the temperature range somewhere between about 250 and 400°C, which coincides with the brittle-ductile transition condition. Note that at 50% pressure equilibration, fluid pressure is significantly above hydrostatic pressure and hydraulic head gradients will drive fluid upward. This clearly demonstrates that it is nearly impossible to drive meteoric fluids down to crustal zones where ductile deformation dominates (and where temperatures are on the order of 250 to 400°C) for prolonged periods of time. Therefore, downward and upward flow of meteoric fluids could not take place simultaneously along extensional faults, and an alternative model in which meteoric fluid infiltration and release are separated in time is required.

3.2. Origin and Burial of Meteoric Fluid

The German KTB drill hole shows that meteoric fluid can penetrate deep below the erosion surface. Fluid pressure is hydrostatic down to 9 km where the temperature is about 275°C [45]. The close to hydrostatic pressure at this depth and temperature is consistent with the results shown in Figure 1, where time to pressure equilibration starts to reduce significantly just above this temperature, except when using the “softest” flow law of Paterson and Luan [41]. Porosities are about 1%, with local excursions up to ≥5% in fault zones, at least down to 4 km depth. Below 4 km, the porosity decreases somewhat [45], possibly down to about 0.1% [46]. With a geothermal gradient of 28°C/km, the decline in porosity in the KTB drill hole occurs at about 110°C. The depth of fluid infiltration is probably strongly temperature-dependent (see above). Infiltration may therefore be double the depth in old cratonic basement rocks with a low geothermal gradient. Below an unconformity, one can thus expect several kilometres of thick zone permeated by a significant amount of meteoric water [3, 40]. Renewed sedimentation will bury the infiltrated zone, together with its pore fluid. Burial may be enhanced by an orogenic event, with the emplacement of nappes on top of the unconformity. Pore fluid will become trapped during burial as pore space becomes disconnected and close. At the end of burial, a package of rock with a more or less pressure-equilibrated meteoric fluid is brought down to deeper crustal levels.

3.3. Preserving the Meteoric Isotope Ratios in a Pore Fluid

A crucial question to address now is how long the meteoric isotope signature of the pore fluid can be preserved. Since the postulated deep infiltration of meteoric fluids is mostly based on δ18O values, we only consider these here.

We first consider oxygen and quartz, which is one of the most abundant crustal minerals. For a rock volume at rest, we assume that the isotopic equilibration is controlled by diffusion of the isotope in the minerals that surround the pores. A pore of size must exchange isotopes with the surrounding mineral in a zone with a width in the same order of . A pore is thus approximately isotopically equilibrated in a time () that is needed to diffuse the mean distance , given by , where is the diffusion coefficient. Using the preexponential constant ( m2/s) and activation energy (243 kJ/mol) of Farver and Yund [47], we see that it takes many millions of years to diffuse >1 μm when the temperature is below about 300°C (Figure 2).

This simple approximation of the isotopic equilibrium time for fluid in rocks shows that fluids can retain their oxygen-isotopic signature for long periods of time (millions of years) below the brittle-ductile transition at about 300-350°C, which is also the temperature range where dynamic recrystallisation sets in. This, however, only applies when the rocks remained at rest. Deformation, such as cataclasis, would enhance isotopic exchange and equilibration at lower temperatures (as envisaged by [27]). This process only affects a small portion of the total rock volume, allowing most of the fluid to keep their original signature at low temperatures. Compaction of a rock by pressure solution [48, 49] involves dissolution of solid and its reprecipitation. This would lead to a quick equilibration of isotopes in the pore fluid and the solid. Whereas this is likely in the case of initially high-porosity sediments, metamorphic or igneous basement rocks do not compact, as their initial porosity is already low. It should further be noted that the equilibrium rates used here are several orders of magnitude slower than those reported in Person et al. [23]. The difference is caused by these authors using a reaction-controlled equilibrium rate, which does not apply if the rock is at rest and the exchange of isotopes is controlled by the solid-state diffusion towards the fluid-mineral interface. Once rocks reach the ductile regime, dynamic recrystallisation, including grain boundary migration, can greatly enhance fluid-rock interaction.

The same rationale can be applied to hydrogen isotopic ratios. A difference with oxygen, however, is that hydrogen is much less abundant in rocks than oxygen. About 8% of all hydrogen atoms reside in the pore fluid in the case of a pure biotite or muscovite schist with 1% porosity. In most cases, therefore, the hydrogen reservoir in the pore fluid is of the same order of magnitude as that of the rock, and the hydrogen isotopic ratio of a stagnant pore fluid will never be completely overprinted by that of the rock.

3.4. Exhumation and Release of Fluids

A characteristic of shear zones and faults with a meteoric signature is that they formed during (rapid) exhumation. The detachment faults of core complexes are prime examples (e.g., [23]). Exhumation can lead to the release of pore fluids, due to the difference in compressibility of fluid and rock [50]. Upon exhumation, the lithostatic pressure in a rock volume is reduced instantaneously. Fluid pressure does not immediately record this decompression, as the fluid-filled pores adjust their volume more slowly. A decrease in effective lithostatic pressure can also be induced by an extensional stress field, as the effective pressure is the average of the overburden and horizontal stresses. Fluid originally in pressure equilibrium with their host will thus become overpressured during exhumation, especially during extension (Figure 3).

The effect of overpressure on fluid-filled pores is not simply the inverse of that of underpressure. Whereas closing an underpressured pore during burial is a stable process, the opening of an overpressured pore to equilibrate solid and fluid pressure is unstable, because rocks cannot sustain significant fluid overpressure [51]. Overpressure is likely to result in (intergranular) cracks that rapidly create a permeability through which the overpressured fluid can flow upward. This leads to interaction with other overpressured pores and the possibility of a runaway effect [52, 53]. This runaway effect would be enhanced if the newly created permeability allows the fluid pressure gradient to decrease towards hydrostatic. This reduces the fluid pressure upstream but increases it downstream (at a shallower depth). This effect leads to the upward propagation of hydrofractures (Figure 4) [54].

Summarising, extension and exhumation lead to relative overpressure in any fluid that is trapped in porosity and has reached approximate pressure equilibrium with the lithostatic pressure when the crust was at rest. Overpressure leads to the creation of permeability, and part of the fluid will migrate upward, driven by the high hydraulic gradient.

3.5. Fluid Budgets

To gain an impression of the amount of meteoric fluid that could pass through a fault or shear zone, we consider a schematic core complex with one major detachment fault. The detachment is situated on the unconformity that separates fluid-impregnated basement rocks (the “infiltration zone”) from a cover series (Figure 5(a)). Movement along the detachment and updoming produce a core complex with the exhumed basement reaching the surface (Figure 5(b)). The offset along the fault is visualised by two split markers just below and above the detachment, separated by the amount of offset on the detachment. The figure shows that a section as long as the offset on the fault of the basement passes each point along the fault. The section passing below the upper black marker is marked in darker grey.

The basement, being stretched and exhumed, releases trapped fluid, which will move upward towards the detachment. The hanging wall does not experience any increase in permeability, as it is moving horizontally or downward, which decreases its permeability. Flow is thus focused into the detachment fault. This way, the detachment experiences a high flux of fluids that came from deeper and hotter levels in the crust. Most fluid-rock exchange is expected at the detachment, as is observed in, for example, the Columbia River Detachment [21].

All the fluid released by the grey rectangle passes below the upper black marker in Figure 5. If the offset along the detachment is 15 km and the thickness of the infiltration zone is 5 km, with a 1% porosity, the total available fluid budget is 75·106 m3 per 1 m width of the fault zone. Not all of this fluid needs to be liberated by the decompression resulting from extension, doming, and movement along the detachment. A conservative estimate would be about 4%, which is equivalent to the fraction of fluid release needed to equilibrate fluid and rock pressure if the fluid is overpressured by 100 MPa (equivalent to about 4 km exhumation, see [50]). If this 4% of the available fluid is released and focused in a 100 m wide shear zone, the integrated fluid flux is 3·104 m3/m2. At 5% porosity in the shear zone, and a duration of 200,000 years for the event [23], this is equivalent to a flow rate of 3 m/y. This qualitative example shows that significant fluid fluxes can be generated that are in the order of those postulated for the Columbia River Detachment [23], but with a radically different source of the meteoric fluids.

The integrated fluid flux is highest at the top of the shear zone, as deeper levels of the shear zone arrived at the reference point later. It is expected that any alteration of the shear zone increases towards the top of the shear zone and then decreases abruptly into the hanging wall, which is exactly as observed in the Columbia River Detachment [23].

4. Conclusions

Meteoric fluid signatures are often found in midcrustal shear zone minerals. This does not necessarily mean that these fluids infiltrated at the time that these shear zones were active, as is often assumed. Meteoric signatures can be preserved for long periods of time, as long as the temperature does not exceed about 300-350°C for oxygen in a quartz-dominated rock. This means that meteoric water can infiltrate down from an erosional surface to be buried and subsequently remobilised many million years later. The crucial point is that the meteoric signature found in shear zones is then from fluid flowing upward, not downward.

Data Availability

All the models and data in this article are based on the equations published here and parameters described in the main text and Table 1. The models can be reproduced with standard spreadsheet software.

Conflicts of Interest

The authors declare that there is no conflict of interest regarding the publication of this paper.

Acknowledgments

This study was developed within the framework of DGMK (German Society for Petroleum and Coal Science and Technology) project 718, funded by the companies ExxonMobil Production Deutschland GmbH, GDF SUEZ E&P Deutschland GmbH, RWE Dea AG, and Wintershall Holding GmbH. EGR acknowledges funding by the AGAUR (Agència de Gestió d’Ajuts Universitaris i de Recerca) of the Generalitat de Catalunya (“Beatriu de Pinós” fellowship 2017SGR-824) and the Spanish Ministry of Science, Innovation and Universities (“Ramón y Cajal” fellowship RYC2018-026335-I and research project PGC2018-093903-B-C22).