# Rapid migration of CO2-rich micro-fluids in calcite matrices

Pierpaolo Zuddas, Stefano Salvi, Olivier Lopez, Giovanni Degiudici, Paolo Censi

► To cite this version:

Pierpaolo Zuddas, Stefano Salvi, Olivier Lopez, Giovanni Degiudici, Paolo Censi. Rapid migration of CO2-rich micro-fluids in calcite matrices. Scientific Reports, 2018, 8, pp.14080. 10.1038/s41598-018-32461-8. hal-01890351

HAL Id: hal-01890351

https://hal.sorbonne-universite.fr/hal-01890351

Submitted on 8 Oct 2018

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
SCIENTIFIC REPORTS

<figure><img src="image_2.png" /></figure>

OPEN

# Rapid migration of CO₂-rich micro-fluids in calcite matrices

Pierpaolo Zuddas<fnref n="1" />, Stefano Salvi², Olivier Lopez³, Giovanni DeGiudici⁴ & Paolo Censi⁵

Received: 11 April 2018

Accepted: 4 September 2018

Published online: 20 September 2018

The transport of supercritical fluids is a determining factor for several geological processes and fundamental in predicting natural resource accumulation and distribution. Calcite, ubiquitous in most geological environments, may contain supercritical CO₂ trapped under the form of fluid inclusions that may move through grain boundaries affecting the rock physical properties. However, despite macroscopic evidence for this process, until recent it was not possible to characterize this process at the nano-scale due to the difficulty of such observations. In this study, we report nanometer-scale observations on calcite crystal surfaces and demonstrate that stress with absence of visible deformation produces fluid leakage from fluid inclusions. Atomic Force Microscopy scanning experiments on freshly cleaved calcite crystals containing visible fluid inclusions revealed the spontaneous formation of nanometer-scale hillocks on flat crystal terraces in only a few minutes, without evidence of surface dissolution. The fact the hillocks formed on flat surface in a short time was unexpected and suggests deposition of material from the inner crystal to the surface through small-scale fluid migration. We estimated the rate of this fluid mobility is by several orders of magnitude higher than the diffusion rate through vacancies estimated in calcite crystals showing that CO₂-rich fluids through micro-pore and nano-pore spaces is in reality much higher than previously assumed using current predictive models.

Fluids play a key role in enhancing mineral reactions and rock deformation, influencing a wide range of geological processes such as fluid-rock interaction, solution transfer, electrical conductivity, seismic wave velocity and attenuation<fnref n="1" />. Observations of minerals in a rock at high magnification has revealed that tiny aliquots of these paleofluids can be retained in microscopic cavities that form fluids inclusions wetting grain surfaces and forming a thin, continuous film that may short-circuit transport. In non-fractured rocks, however, fluids can easily move as they are located at the grain boundary spaces². One significant weakness in our fluid transport knowledge concerns the mechanisms controlling the fluid distribution among grains of solid rock with which they are in chemical and mechanical equilibrium. The fluid flow between grains and mineral pores has been documented at elevated temperatures and pressures³ while it is assumed to be negligible at 25°C and 1 atmosphere⁴. Some earlier evidence to the contrary is, however, constituted by the self-diffusion of calcium isotopes documented in some soil carbonates⁵. In the case of CO₂-bearing aqueous fluids, the extreme conditions required for studying displacement processes at the micrometric and nanometric scale (i.e. the scale of pores and grain boundaries), are such that only a limited number of studies have been published on the subject⁶,⁷. These small-scale phenomena, interpreted as interfacial processes, have been often neglected in continuum-scale modelling such a viscous fingering⁶ and wetting alteration by CO₂-rich fluids⁷,⁸. The rheological properties of thin complex fluids is of paramount importance for gas flow in rock matrices⁹ and geophysical models of non active tectonic areas¹⁰,¹¹. Earlier work¹²⁻¹⁴ suggests that fluids confined at the micrometric (or nanometric) scale may move fast and have a key role in fluid-solid reactivity. Fluid movement and solid permeability at grain boundaries have often been explained in terms of interfacial energy⁷, while when a small amount of fluid is trapped by surface tension, the capillarity pressure is sheared. At this scale, fluids must be present in the form of absorbed molecular films, non-equilibrium or wetting film, as well as isolated fluid inclusions maintained by steric force originated from hydrated layer at the mineral surface¹³,¹⁵.

In this contribution, we report nanometer-scale observations made by Atomic Force microscopy (AFM) on natural calcite crystals containing micrometric fluid inclusions and demonstrate stress-induced production of fluid leakage from what most likely are nanometric fluid inclusions, in a very short time under standard temperature and pressure conditions. Our results yield the first independent determination of nano-scale fluid movement

<footnote n="1">¹Sorbonne Université, CNRS, ISteP, 4 place Jussieu, 75005, Paris, France. ²GET, Université de Toulouse, CNRS, IRD, CNES, Toulouse, France. ³Statoil, Trondheim, Norway. ⁴Department of Chemistry and Earth Sciences, University of Cagliari, Cagliari, Italy. ⁵DIsTEM, University of Palermo, Palermo, Italy. Correspondence and requests for materials should be addressed to P.Z. (email: pierpaolo.zuddas@sorbonne-universite.fr)</footnote>
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Optical microscope photomicrographs (plane polarized transmitted light) of a 120-µm-thick doubly-polished wafer (120-µm in thickness) prepared from a calcite Iceland spar. Images (a,b) show the same area at different focal depths (a nearest the surface and b: a few microns below) to illustrate a healed fracture defined by a plane of fluid inclusions (FI). The sketch in the inset illustrates the FI trapping mechanism. Note somewhat larger fluid inclusions off the plane. Reflection of internal features is due to the double refraction of calcite. Image (c) depicts a micrometric fluid inclusion in a crystal, showing a hairline fracture departing from it due to (partial) thermal deprecitation of the inclusion.</figcaption></figure>

inside a calcite matrix. The potential significance of the results to fluid transport in rock is also briefly discussed. This study points the way for future similar research on other minerals, helpful for modelling fluid transport of matter, fault lubrication during earthquakes¹⁶ and geological sequestration of CO₂¹⁷.

## Results

**Heterogeneity of natural calcite.** We monitored the state of a freshly cleaved calcite surface using optical-quality Iceland spar from Chihuahua, Mexico (99.5% CaCO₃). Optical microscope observations of the calcite crystal denoted the presence of several fracture planes of few millimeters in length, as well as fluid inclusions organized along similar planes or occurring isolated or in small groups (Fig. 1a,b). The size of visible fluid inclusions ranges from <1 micron to relatively large (~20 micron), with shapes varying from a rhombohedrum to irregular morphology. Based on the number of phases present at room temperature and on microthermometric data, we identified two types of fluid inclusion: Type-1 fluid inclusions contain liquid H₂O and a relatively high-density homogeneous CO₂ phase (in a few fluid inclusions, CO₂ consisted of distinct liquid and vapor phases); Type-2 fluid inclusions, in addition to aqueous and carbonic phases, contain halite and, in a few cases, other unidentified solids. Both fluid inclusion types occur together in planes and groups. On cooling, the homogeneous CO₂ separated into liquid and vapor phases that, upon heating, homogenized to the liquid phases at temperatures ranging from -5° to 27°C. The corresponding bulk density ranging from 0.85 to 1.15 g/cm³, which would produce internal pressures as high as 50 MPa at room temperature. The nature of dissolved salts in the aqueous solutions of both inclusion types, estimated by microthermometry (see method section), consists mostly of Na⁺ with salinity in the range of 5–26 wt. % NaCl equivalent. Some halite-bearing inclusions have higher salinity, the exact degree of which could not be determined because most inclusions deprecitated before total homogenization upon heating, likely due to high CO₂ pressure. High Resolution Transmission Electron Microscopy observations reveal the presence of rupture of the periodicity of calcite producing several grain junctions delimiting variable size of mono-crystal domains ranging from 10 to 500 nm in diameter.

**AFM observation and evidences for fast fluid movement at standard conditions.** Examination of over 100 freshly cleaved calcite surfaces by AFM showed the presence of atomically flat terraces separated by cleavage steps of various heights. During real-time monitoring of samples, after only two minutes of scanning,
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Evolution of the state of the (1014) calcite surface during 18 min of AFM scanning at a scan frequency of 6 Hz and an applied force of 7 nN. The scan surface area is 5 *5 µm². After 2 min of scanning, randomly distributed hillocks appear over the investigated surface area. The average height is 9 ± 2 nm while size and geometry change. After 13 minutes of scanning, the growth stops.</figcaption></figure>

we observed nanometer-scale hillocks forming spontaneously on the calcite surface. These mounds did not concentrate on surface imperfection sites such as steps or kinks, which are recognized sites of excess energy, but appeared to be randomly distributed. Also, although hillocks grew with a preferential horizontal direction, they did not present any particular morphology or crystallographic pattern (Fig. 2). Their height averaged around 0.9 ± 0.2 nm, whereas their horizontal size and geometry varied slightly (0.01 to 0.5 µm). After 13 minutes of scanning, growth of the earliest-formed hillocks stopped, however, within this time range new hillock nuclei had
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Evolution of the hillock volume as a function of time at different scan frequencies (in Hz). At lower frequencies, hillocks grow rapidly during the first 10 min. This fast growth rate decreases at higher scan frequencies. Lower AFM scan frequencies correspond to a higher time of contact between tip and surface.</figcaption></figure>

appeared at other sites on the observed mineral surface. In this study 80% of the observation reflect the feature reported in Fig. 2. The fact that hillocks formed spontaneously on flat terraces in only few minutes, without evidence of surface dissolution, and under standard atmospheric conditions, was unexpected.

In an earlier study¹⁸ hillocks appeared only after 2–5 hours of AFM scanning on a mineral surface that presents evidence of dissolution pits. Changes of the cleaved calcite surface have also been observed in several X-ray reflectivity studies of air-calcite interaction¹⁹,²⁰ and were explained by the absorption of water molecules coming from the surrounding air moisture²¹. However, our AFM experiments put into question the fact that the origin of the water molecules is exclusively due to the dampness in the air. The extreme rapidity of the phenomenon that we observed, and the absence of dissolution pits, eliminates the possibility that hillocks formed by dissolution/re-precipitation of small parts of the calcite surface. In addition, hillock growth was observed exclusively in the surface areas that were investigated, suggesting a possible role of internal characteristics of our calcite sample. We therefore set out to test the hypothesis that fluid inclusions of sub-micron size and located a few microns below the calcite surface are the source of at least some of this matter. In our samples, fluid inclusions consist of an aqueous solution enriched in $Cl^-$ and $Na^+$, coexisting with liquid $CO_2$, confined at an internal pressure estimated at around 50 MPa. The strength of engagement of the AFM tip under our experimental conditions is by several orders of magnitude lower than the calcite Shear Modulus (modulus of rigidity measuring the stiffness of calcite macroscopic crystals, which is about $10^{10}$ Pa)²². The surface free energy of ‘dried’ large calcite crystals is equal to 0.59 J/mol, whereas that of ‘hydrated’ crystals is 0.33 J/mol²³ and may decrease by one order of magnitude for smaller nanometric sized particules²⁴. Assuming that the fluid inclusions are in equilibrium with the mineral walls, the stress produced by the tip could be responsible for the disruption of such conditions. Thus, sub-micron fluid inclusions can move through the calcite crystal (via dislocations) to the mineral surface and cause the formation of the observed hillocks. Migration of several-micron-sized fluid inclusion has been postulated by several authors in minerals such as quartz, via dislocation or by dissolution-precipitation of the crystal walls²⁵,²⁶.

**Estimating the fast micro-fluid movement.** We carried out a number of additional experiments to investigate the interaction between the AFM tip and the mineral surface. Since scan frequency of an AFM tip over a surface is inversely proportional to the tip-surface residence time, by decreasing the scan frequency we could increase the contact time (for a same scanned surface), thus enhancing the mechanical strain undergone by the mineral. Resulting data shows that a 50% decrease in scan frequency produced a 600% increase in hillock volume (Fig. 3). Moreover, hillock growth stopped more rapidly with higher scan frequencies. We also tested the influence of temperature by varying this parameter in the AFM cell from 5 to $50 \pm 4$ °C, and found that hillock volumes increased in high temperature runs, within similar observation time spans. We evaluated the volume variation, by estimating the dimensions of discrete feature in the AFM images²⁷ and found that a 30 °C temperature increase produced a 4-fold increase in hillock volume (Fig. 4). At a temperature of 5 °C, however, hillocks did not form despite several hours of scanning. These results indicate that mechanical strain controls hillock formation and that this process is strongly enhanced at higher temperatures.

The observed hillock formation can be explained by the fact that $CO_2$-rich fluid inclusions located just under the surface are near their rupture limit because of the elevated internal pressure, thus even a very small mechanical strain will cause them to decrepitate and lose their fluid. This phenomenon is potentially enhanced by heat, as any temperature increase would cause the internal pressure to rise dramatically due to conservation of internal density²⁸. Fluid inclusion decrepitation in response to heat-induced overpressure is a well-known phenomenon in fluid-inclusion studies; hairline fracture propagating from the decrepitated micrometric fluid inclusions can sometimes be observed (cf. Fig. 1c). Alternatively, it is possible that fluid movement to the surface took place along preferential pathways such as micro or nano-fractures already present²⁹,³⁰. However, this would not be consistent with the fact that hillocks were randomly distributed over the calcite surface. Assuming that (i) the fluid inclusions involved in the fluid movement are located at least 0.01 to 0.1 µm below the calcite surface, (ii) the pathway is vertical, and (iii) precipitation at the surface is instantaneous, our real-time measurements allow us to
<figure><img src="image_6.png" /><figcaption>**Figure 4.** AFM images of a freshly cleaved calcite surface at different temperatures (a, 5°C, b, 25°C and c, 55°C) obtained with a scan frequency of 9 Hz. (d) Represents the hillock volume formed after one hour of scanning.</figcaption></figure>

estimate a rate of fluid transport at ~5 × 10⁻¹⁰±⁰.⁵ m·s⁻¹. This value does not represent a diffusion coefficient in the sense of the Ficks Law but is simply an estimation of the rate of inner fluid-mobility in our calcite matrix³⁰. The transport rate that we calculated is about 4–5 orders of magnitude higher than that estimated by low-temperature extrapolation of solid-state diffusion vacancy experiments³¹. These results highlight the extreme mobility of fluids at the micro-nano fluidic scale in calcite crystals at standard pressure and temperature conditions.

AFM observations alone cannot provide information on the composition of these surface features. However, the fact that the hillocks grow on the surface may correspond to an eventual accumulation by epitaxial growth³². The CO₂ coming from the fluid inclusions could mix with the thin water layer in equilibrium with ambient air at the surface³³, inducing precipitation of calcium carbonate hillocks. Alternatively, these hillocks could consist of halite nanocrystals resulting from evaporation of the NaCl-rich aqueous fluid coming from fluid inclusions. However, halite precipitates formed from heat-induced deprecitation of fluid inclusions generally display varied morphologies, ranging from cubic crystals to toothpaste-like tubes²⁸, features that were not recognized in this study. We also carried out additional experiments under dry atmosphere conditions (using toluene) and found much smaller hillocks, suggesting the importance of ambient water in the formation of the hillocks and thus validating the former hypothesis.

**Interfacial transport at the micro-fluidic scale.** Our experimental observations show that stress applied to a calcite crystal can produce fluid migration, at least within the more external parts of the crystal. The question that arises, then, is the limiting conditions for fluid mobility and when it could be a dominant transport mechanism. In a stable fluid inclusion, the liquid resistance to driving forces increases when the channel size decreases¹³. However, in case of micrometric or nanometric channel sizes, transport coefficients should take interfacial properties into account³⁴,³⁵. For simple fluid/solid interface with pressure flow represented by a slit of thickness h and slip length b, the mean velocity increases by a factor 1 + 6b/h as compared to the no-slip surfaces, suggesting that the size of the channel h drives the flow. However, our observations suggest that the interfacial structure at
the solid-liquid interface could be an alternative way to generate flow. The electrical interaction between tip and surface atoms generates a stress over the mineral surface. Applying a stress E, the fluid acquires a plug flow comparable to a velocity profile, with a velocity $v_{EO}$ proportional to E in agreement with the Smoluchowski equation³⁴:

$$v_{E_0} = -\varepsilon\zeta E/\eta \quad (1)$$

where $\varepsilon$ is the dielectric permittivity of the solvent; $\eta$ its viscosity; $\zeta$ is the Zeta potential of the surface, which is traditionally assumed to match the electrostatic potential at the position where the velocity profile vanishes. This interface transport, which originates within the Debye layer at the interface, quantifies the width of the interfacial regions at the charged surface. The Debye layer, of size, $\lambda_{D}$, results from the competition between ion attraction at the charged surface and the entropic effects. The value of $\lambda_{D}$, is generally function of the salt concentration in the fluids ($\lambda_{D}$ is about 3–4 nm for a concentration of salt of 0.01 M). The interfacial Debye layer is electrically charged, while the remaining bulk solvent is electrically neutral. Therefore, when applying the stress, the electric driving force is located within the Debye layer only. The driving force per unit surface acting on the liquid thus becomes $Fe = \sigma E$, with $\sigma$ the charge of the Debye layer, which is exactly opposite to the surface charge $\sigma = \varepsilon V_0/\lambda_{D}$ for weak surface potential $V_0$. This leads to $Fe \sim -\varepsilon(V_0/\lambda_{D})E$. Now, the fluid velocity results from the balance between this electrical driving force in the Debye layer and the viscous stress at the surface. For a no-slip boundary condition, the latter is $F_{\eta} \sim \eta v_{EO}/\lambda_{D}$, since the velocity varies on a scale given by the Debye length. If slippage is exhibited at the surface, the velocity varies on the size $\lambda D + b$ and the viscous stress becomes accordingly $F_{\eta} \sim \eta v_{EO}/(\lambda_{D} + b)$, lower than the no-slip boundary conditions. Gathering these results, we get:

$$v \frac{v_{EO}}{\lambda_{D} + b} \approx -\varepsilon \frac{V_0}{\lambda_{D}} E \quad (2)$$

This leads to the Smoluchowski formula and provides an expression for the zeta potential in the form:

$$\zeta = V_0 (1 + b/\lambda_D) \quad (3)$$

Equation 3 shows that for a non-slipping surface, $V_0$ can be identified as the potential at the plane of shear $V_0 = V_{(Zs)}$, with Zs being of the order of one liquid layer. As expected, a strong amplification of charge transport is therefore demonstrated on a slipping surface. Generalizing this behavior to all interfacial transport phenomena, the flow fluid velocity in confined media results proportional to the applied gradient (stress, electrical potential) and depends on the size L, of the interface where the fluids interacts specifically with the mineral surface. Since our observed phenomenon is proportional to the applied stress that produces a surface energy gradient, we propose that slippage amplifies the transport according to the ratio b/L, which in turn controls the interfacial transport of fluids at the micro and nanometer scale. While structural defects at the external surface of minerals have been identified as limiting steps in transport and reactivity in non-confined fluids³², fast movement of CO₂-rich fluids are related to fluid-fluid and fluid-calcite interfaces.

## Discussion

Our unexpected findings have potentially far-reaching implications for developing new fluid transport prediction models. Traditional approaches to fluid transport assume the presence of inter-grain pathways in polycrystalline materials such as micro and nano-pore tubes³⁶,³⁷. Fluid transport is usually observed in poly-phase material but individual mineral grains can also exhibit clear indications of multi-path migration even when visible evidence of such paths is lacking. Quantifying this process at the micro-nanoscale for a mineral at standard external conditions, we show that fluid mobility is in reality much higher than previously assumed using current predictive models⁸,³⁸⁻⁴⁰. Our findings also suggest that mechanical stress plays a more significant role in fluid transport than previously assumed. Classical reasoning assumes local partition equilibrium at the fluid-solid interface, or, equivalently, a continuity of the chemical potential. Based on our observations and on the approach outlined above, it would seem more consistent to introduce a discontinuous flux across the interface. The magnitude of this discontinuity characterizes the kinetics of pressure solution rather than intergranular diffusion. The rheological properties of small-size fluids and their flow in carbonate rocks are of paramount importance in many geological situations. In the case of lubrication of faults, calcite is an important, and in some cases dominant mineral in seismically active regions worldwide, leading an extra anisotropy in P- or S-wave velocities. Fast fluid movement at the microscale and nanoscale have significant consequences in interpreting anisotropy in fluid-bearing rock systems during carbon crust cycling⁴¹⁻⁴³ and geological CO₂ sequestration¹⁷,⁴⁴. High mobility of fluids in crystals may also have important implications for fluid-mediated metal transport. The origin of native gold grains observed at the surface of sulfide grains in orogenic-type deposits⁴⁵ could be explained by movement of micron-size fluid inclusions commonly observed in sulfides rather than by dissolution/precipitation processes involving external fluids⁴⁶. Relatively low-strain (lower than Shear Modulus) could easily trigger micro-fluids containing dissolved gold (and other metals) to migrate out of the crystals, precipitating native metal like gold on the surface with no deformation of the parent sulfide minerals⁴⁷, thus explaining field-scale petrographic observations.

## Methods

**Fluid inclusion analysis.** The fluid inclusions were identified by preparing wafers of 150 µm thickness from the calcite crystals, polished on both sides to optimise transparency. The fluid composition was estimated non-destructively, by measuring the temperatures of phase change inside the fluid inclusions, using a Linkam THMGS 600 heating-freezing stage and following the procedures in Roedder (1984) and Shepherd et al. (1985). The stage was calibrated using synthetic material to accuracies of ±0.1 °C in the temperature range covered in this study. Salinity (expressed as wt.% eq. NaCl), bulk composition, and density data were calculated by reducing raw
thermometric data with the BULK and ISOC software packages of Bakker)⁴⁸. The equations of state used was that of Bowers & Helgeson (1983) for the H₂O–CO₂–NaCl system.

**Experimental Details.** Atomic Force Microscopy (AFM) experiments were carried out under controlled atmospheric conditions with 70% residual humidity, using a Digital Instruments Nanoscope III device. In order to eliminate the possibility of residual dust dissolution in the water film on the calcite surface (cf. below), freshly cleaved grains were cleaned with a N₂ flux just before experiments. Therefore, in order to avoid any interaction between tip and surface, the normal force applied to the samples was set to lower than 10 nN. We carried out all standard tests for imaging calcite samples: rotation, scaling, reproducibility on other sites and samples.

**References**

1. Ague, J. J. Fluid Flow in the Deep Crust. In *Treatise on Geochemistry* (eds Holland, H. D. and Turekian, K. K.), 2nd Edition, vol 4, Elsevier, Oxford, 203–247 (2014)

2. Desbois, G., Urai, J. L. & de Bresser, J. H. P. Fluid distribution in grain boundaries of natural fine-grained rock salt deformed at low differential stress (Qom Kuh salt fountain, central Iran): Implications for rheology and transport properties. *J. of Struct. Geol.* **43**, 128–143 (2012).

3. Farver, J. R. Oxygen self-diffusion in calcite: Dependence on temperature and water fugacity. *Earth Planet. Sci. Lett.* **121**, 575–587 (1994).

4. Wark, D. A. & Watson, E. B. Interdiffusion of H₂O and CO₂ in metamorphic fluids at similar to 490 to 690 degrees C and 1GPa. *Geochim. Cosmochim. Acta* **68**, 2693–2698 (2004).

5. Lahav, N. & Bolt, G. H. Self-diffusion of 45Ca into certain carbonates. *Soil Sci.* **97**, 293–299 (1964).

6. Zhang, C. Y., Oostrom, M., Wietsma, T. W., Grate, J. W. & Warner, M. G. Influence of viscous and capillary forces on immiscible fluid displacement: Pore-scale experimental study in a water-wet micromodel demonstrating viscous and capillary fingering. *Energy Fuels* **25**(8), 3493–3505 (1964).

7. Kim, Y., Wan, J., Kneafsey, T. & Tokunaga, T. K. Dewetting of silica surfaces upon reactions with supercritical CO₂ and brines: pore scale studies in micromodels? *Environ. Sci. Technol.* **46**(7), 4228–4235 (2012).

8. Dehoff, K. J., Oostrom, M., Zhang, C. & Grate, J. W. Evaluation of two-phase relative permeability and capillarity pressure relations for unstable displacements in a pore network. *Vadose Zone J.* **11**, 4 (2012).

9. Zhang, P., Liming, H., Meegoda, J. N. & Gao, S. Micro/Nano-pore network analysis of gas flow in shale matrix. *Scientific Reports* **5**, 13501 (2015).

10. Bailey, R. C. Trapping of aqueous fluids in the deep crust. *Geophys. Research Lett.* **17**(8), 1129–1132 (1990).

11. Bailey, R. C. Fluid trapping in mid-crustal reservoirs by H₂O-CO₂ mixtures. *Nature* **371**, 238–240 (1994).

12. Stevenson, D. J. On the role of surface tension in the migration of melt and fluids. *Geophys. Res. Lett.* **13**(11), 1149–1152 (1986).

13. de Meer, S., Nakashima, S. & Spiers, C. J. Structure and diffusive properties of fluid-filled grain boundaries: an *in-situ* study using infrared (micro) spectroscopy. *Earth and Planetary Science Letters* **232**(3-4), 403–414 (2005).

14. Yardley, B. W. D. & Bodnar, R. J. Fluids in the continental crust. *Geochemical Perspectives* **3**(1), 1–127 (2014).

15. Smith, S. A. F., Nielsen, S. & Di Toro, G. Strain localization and the onset of dynamic weakening in calcite fault gouge. *Earth and Planetary Science Letters* **413**, 25–36 (2015).

16. Di Toro, G. et al. Fault lubrification during earthquakes. *Nature* **471**, 494–498 (2011).

17. Matter, J. M. et al. Rapid carbon mineralization for permanent disposal of anthropogenic carbon dioxide emissions. *Science* **352**(6291), 1321–1324 (2016).

18. Stipp, S. L. S., Konnerup-Madsen, J., Franzreb, K., Kulik, A. & Mathieu, H. J. Spontaneous movement of ions through calcite at standard temperature and pressure. *Nature* **396**, 356–359 (1998).

19. Stipp, S., Gutmannsbauer, W. & Lehmann, T. The dynamic nature of calcite surfaces in air. *American Mineralogist* **81**(1-2), 1–8 (2015).

20. Fenter, P. et al. Surface speciation of calcite observed *in situ* by high-resolution X-ray reflectivity. *Geochim. Cosmochim. Acta* **64**, 1221–1228 (2000).

21. Geissbuhler, P. et al. Three-dimensional structure of the calcite-water interface by surface X-ray scattering. *Surf. Science* **573**, 191–203 (2004).

22. Jackson, I. & Paterson, M. S. Shear modulus and internal friction of calcite rocks at seismic frequencies: pressure, frequency and grain size dependence. *Physic of Earth and Plan. Int.* **45**, 349–367 (1987).

23. De Leeuw, N. H. & Parker, S. C. Surface Structure and Morphology of Calcium Carbonate Polymorphs Calcite, Aragonite, and Vaterite: An Atomistic Approach. *J. Phys. Chem. B* **102**, 2914–2922 (1998).

24. Forbes, T. Z., Radha, A. V. & Navrotsky, A. The energetics of nanophase calcite. *Geochim. Cosmochim. Acta* **75**, 7893–7905 (2011).

25. Bakker, R. J. & Jansen, B. H. Preferential water leakage from fluid inclusions by means of mobile dislocations. *Nature* **345**, 58–60 (1990).

26. Heggie, M. A molecular water pump in quartz dislocations. *Nature* **355**, 337–339 (1992).

27. Bickmore, B. R., Rufe, E., Barrett, S. & Hochella, M. F. Measuring discrete feature dimensions in AFM images with image SXM. *Geological Materials Research* **1**(5), 8 (1999).

28. Roedder, E. Fluid inclusions. *Reviews in Mineralogy (Mineral. Society of America)* **12**, 646 (1984).

29. Nakamura, M., Yurimoto, H. & Watson, B. Grain growth control of isotope exchange between rocks and fluids. *Geology* **33**(10), 829–832 (2005).

30. Watson, E. B. & Baxter, E. F. Diffusion in solid-Earth systems. *Earth Planet. Sci. Letter* **253**, 307–327 (2007).

31. Poirier, J. P. Creep of crystals. *Cambridge Earth Sciences Serie*, 274 (1985).

32. Putnis, A. Why mineral interface matter. *Science* **343**(6178), 1441–1442 (2014).

33. Sturchio, N. C. et al. Lead adsorption at the calcite-water interface: Synchrotron X-ray standing wave and X-ray reflectivity studies. *Geochim. Cosmochim. Acta* **61**, 251–26 (1997).

34. Bocquet, L. & Barrat, J. L. Flow boundary conditions from nano- to micro-scales. *Soft Matter* **3**, 685–693 (2007).

35. Bocquet, L. & Charlaix, E. Nanofluidics, from bulk to interfaces. *Chem. Soc. Rev* **39**, 1073–1095 (2010).

36. Siffert, B. & Fimbel, P. Parameters affecting the sign and the magnitude of the electrokinetic potential of calcite. *Coll. Surf.* **11**, 377–389 (1984).

37. Baxter, E. F. & DePaolo, D. J. Field Measurement of Slow Metamorphic Reaction Rates at Temperatures at 500–600°C. *Science* **288**, 1411–1414 (2000).

38. Blunt, M. J., Jackson, M. D., Piri, M. & Valavnte, P. H. Detailed physics, predictive capability and macroscopic consequences for pore-network models of multiphase flow. *Adv. Water Resour.* **25**(8-12), 1069–1089 (2002).

39. Fernandez, J. F., Rangel, R. & Rivero, J. Crossover length from invasion percolation to diffusion-limited aggregation in porous-media. *Phys. Rev. Lett.* **67**(21), 2958–2961 (1991).

40. Skelton, A. Flux rates for water and carbon during greenschist facies metamorphism. *Geology* **39**, 43–46 (2011).

41. Evans, K. Metamorphic carbon fluxes: how much and how fast? *Geology* **39**, 95–96 (2011).
42. Frezzotti, M. L., Selverstone, J., Sharp, Z. D. & Compagnoni, R. Carbonate dissolution during subduction revealed by diamont-bearing rocks from Alps. *Nature Geosciences* **4**, 703–706 (2011).

43. Malusà, M. G. et al. Active carbon sequestration in the Alpine mantle wedge and implication for long-term climate trends. *Scientific Report* **8**, 4740–4748 (2018).

44. Vialle, S. et al. Percolation of CO₂-rich fluids in a limestone sample: Evolution of hydraulic, electrical, chemical, and structural properties. *J. Geophys. Res. Solid Earth* **119**, 2828–2847 (2014).

45. Salvi, S. et al. The Pampe gold deposit (Ghana): Constraints on sulfide evolution during gold mineralization. *Ore Geology Reviews*, v. **78**, 673–686 (2016).

46. Fougerouse, D. et al. The golden ark: Arsenopyrite crystal plasticity and the retention of gold through high strain and metamorphism. *Terra Nova* **28**, 181–187 (2016).

47. Velásquez, G. et al. Formation and Deformation of Pyrite and Implications for Gold Mineralization in the El Callao District, Venezuela. *Economic Geology* **109**, 457–486 (2014).

48. Bakker, R. J. Package FLUIDS 1. Computer programs for analysis of fluid inclusion data and for modelling bulk fluid properties. *Chemical Geology* **194**, 3–23 (2003).

**Acknowledgements**

Funding was provided by E.U. project Greenhouse-gas Removal and Sorbonne Université grants. We thank J.P. Poirier (French Academy of Science) for discussions.

**Author Contributions**

P.Z., O.L. and G.D.G. design the work and performed the experiments. S.S. and P.C. performed the fluid chemical analysis. P.Z., S.S., O.L., G.D.G. and P.C. analyzed the data and wrote the manuscript.

**Additional Information**

**Competing Interests:** The authors declare no competing interests.

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_7.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2018