- Research article
- Open Access
Mineral precipitation-induced porosity reduction and its effect on transport parameters in diffusion-controlled porous media
© Chagneau et al. 2015
- Received: 7 January 2015
- Accepted: 5 August 2015
- Published: 3 September 2015
In geochemically perturbed systems where porewater and mineral assemblages are unequilibrated the processes of mineral precipitation and dissolution may change important transport properties such as porosity and pore diffusion coefficients. These reactions might alter the sealing capabilities of the rock by complete pore-scale precipitation (cementation) of the system or by opening new migration pathways through mineral dissolution. In actual 1D continuum reactive transport codes the coupling of transport and porosity is generally accomplished through the empirical Archie’s law. There is very little reported data on systems with changing porosity under well controlled conditions to constrain model input parameters. In this study celestite (SrSO4) was precipitated in the pore space of a compacted sand column under diffusion controlled conditions and the effect on the fluid migration properties was investigated by means of three complementary experimental approaches: (1) tritiated water (HTO) tracer through diffusion, (2) computed micro-tomography (µ-CT) imaging and (3) post-mortem analysis of the precipitate (selective dissolution, SEM/EDX).
The through-diffusion experiments reached steady state after 15 days, at which point celestite precipitation ceased and the non-reactive HTO flux became constant. The pore space in the precipitation zone remained fully connected using a 6 µm µ-CT spatial resolution with 25 % porosity reduction in the approx. 0.35 mm thick dense precipitation zone. The porosity and transport parameters prior to pore-scale precipitation were in good agreement with a porosity of 0.42 ± 0.09 (HTO) and 0.40 ± 0.03 (µ-CT), as was the mass of SrSO4 precipitate estimated by µ-CT at 25 ± 5 mg and selective dissolution 21.7 ± 0.4 mg, respectively. However, using this data as input parameters the 1D single continuum reactive transport model was not able to accurately reproduce both the celestite precipitation front and the remaining connected porosity. The model assumed there was a direct linkage of porosity to the effective diffusivity using only one cementation value over the whole porosity range of the system investigated.
The 1D single continuous model either underestimated the remaining connected porosity in the precipitation zone, or overestimated the amount of precipitate. These findings support the need to implement a modified, extended Archie’s law to the reactive transport model and show that pore-scale precipitation transforms a system (following Archie’s simple power law with only micropores present) towards a system similar to clays with micro- and nanoporosity.
- Mineral precipitation
- Reactive transport
- Porous media
- Through diffusion
- Archie’s law
Understanding how changing porosity as a results of chemical reactions affects permeability and tortuosity is of paramount importance for describing various hydrogeological processes, such as fluid circulation in acid mine drainage remediation  or hydrothermal systems [2, 3], management of aquifer recharge , efficiency of permeable reactive barriers  and the long-term geochemical evolution of host rock considered for high-level nuclear waste repositories [6, 7]. In some cases a porosity and flow reduction (clogging) effect is desired and intended to inhibit fluid and solute migration (e.g. mitigation strategies for CO2 leakage along wellbores  or for microbial-enhanced-hydrocarbon-recovery ). Describing clogging phenomena in a reactive transport code is challenging , particularly over the long periods (1,000 to 1,000,000 years) that are of interest in most safety assessment cases. In addition to validating numerical schemes, experiments validating the laws implemented for coupling the geochemistry (e.g. dissolution-precipitation reactions) to the transport parameters (induced porosity changes) and overall fluid flow processes (e.g. Carman-Kozeny relationship or Archie’s law) are required.
In this study, we attempted to apply Archie’s law in a 1D single continuum reactive transport model for a simple sand system, reducing the pore space by pore scale precipitation. Numerical model inputs such as porosity, effective diffusion coefficient and the cementation factor applied in the Archie’s law have been derived from the experimental data obtained prior to the porosity reduction experiments. Afterwards, a systematic study of diffusion-controlled precipitation of a simple salt, celestite (SrSO4), in the poorly compacted sand was performed. Double-reservoir diffusion experiments are widely used in the literature to characterize the transport properties of porous materials [13, 18, 19] or in so called two-way U-tube experiments to control the precipitation of various mineral phases [20, 21] in porous silica gel. Our experiments presented here are a combination based on those two approaches mentioned above. The second method (3D imaging coupled to pore-morphology modeling) became of great interest in the past two decades in many domains, including CO2 sequestration [22, 23] and soil sciences , as a powerful and non-destructive analysis tool for characterizing a porous material and predicting how it will evolve under various petro-physical, geochemical, and hydrogeological conditions [25–28].
Our experimental approach consisted of three complementary techniques believed to characterize the reactive system in the most complete way: (1) conservative tracer experiments before and during the formation of a precipitation front, (2) computed micro-tomography (µ-CT) coupled to pore-morphology modeling before and after the formation of a precipitation front, and (3) post analysis of the precipitate (dissolution, SEM-EDX). The experimental data were used as input parameters for a 1D reactive transport model based on the CRUNCHFLOW code . The calculated flux, mineral mass formed and porosity reduction were then compared to experimental results obtained after celestite precipitation.
Initial transport parameters
Summary of the average initial transport parameters of compacted sea sand obtained from three analytical techniques: MIP, tracing experiments and µCT
×10−10 m2 s−1
0.42 ± 0.02
0.42 ± 0.09
4.48 ± 0.17
2.11 ± 0.42
0.40 ± 0.03
2.21 ± 0.43
2.72 ± 0.39
Evolution of transport parameters during porosity reduction (clogging)
The precipitation front that was visible after 7 days to the naked eye formed a sharp disk located near the middle of the diffusion columns, slightly closer to the Sr reservoir. They appeared as white disks less than 1 mm thick, perpendicular to the direction of diffusion. The exact position was 21 ± 1 and 18 ± 2.5 mm from the filter supports of the Sr-reservoir for Prec1 and Prec3, respectively. The front in diffusion cell Prec1 was imaged with µ-CT at two different resolutions (5.89 and 15.2 µm), and the 3D volume was implemented in a pore-morphology model to estimate its influence on the initial transport parameters (porosity, diffusivity, and tortuosity). The Laplace field calculations consist of applying a non-reactive particle concentration gradient over the porous material and running a random-walk simulation, tracing the particles over a selected running time.
We saw the effect of porosity reduction on the tracer diffusivity in diffusion cells Prec2 and Prec3 after about 7 days of reaction (Fig. 2). The HTO flux while the pore space was being reduced by celestite precipitation significantly decreased compared with the breakthrough curves obtained in the materials’ initial state. The HTO flux increased during the first 6 days in a fashion similar to the behavior in the unreacted media. Afterwards, it decreased quickly until it reached steady state at a lower level at approximately day 15. The steady state flux had a value of 8 ± 1 Bq m−2 s−1, about 4.6 times lower than the fluxes obtained for the control (no precipitation) experiments (37 ± 2 Bq m−2 s−1). The fact that there was measurable flux even after 66 days of celestite precipitation indicates that at least part of the pore space in the reaction zone remained open and connected.
Characterizing the precipitate
µCT imaging of a precipitation front
Porosity and precipitate mass quantification by low resolution µCT
We attempted to estimate the mass of precipitate in the pore space and the remaining porosity using a simple pixel-counting treatment. We estimated a volume of celestite of 5.1 ± 0.8 mm3 and a corresponding mass of 20.2 ± 3.2 mg. This mass is in good agreement with the dissolved volume of celestite in the 5 mm thick sample containing the precipitation front detected in the post mortem analysis of the same diffusion cell, which gave 21.7 ± 0.4 mg.
The initial porosity of the sand column was estimated with a similar pixel-counting approach for the pore space in the regions far from the precipitation zone. An initial porosity was estimated at 0.44 ± 0.05, again in reasonable agreement with all other experimental results (0.40–0.42). From the volume of precipitated celestite and the initial porosity, we were able to estimate the minimum porosity reached in the 350 µm thick dense precipitation zone. Considering that approximately 75 % of the precipitate was included in this dense precipitation zone, we estimated that the porosity in this area was ≥0.28 ± 0.03. Unfortunately, it was not possible to estimate the connectivity of this residual porosity with the low spatial resolution of the micro-tomography.
Porosity and precipitate mass quantification by high resolution µCT
The mass of celestite precipitate and the porosity profile were determined with simple voxel-counting on the segmented volume of the high-spatial resolution dataset. This method gave a mass of precipitate of 30.2 mg, a value higher than the estimate determined by dissolving the precipitate during post-mortem analysis of the same diffusion cell (21.7 ± 0.4 mg). This difference in values might be due to the over estimate of the gray values threshold that we selected for celestite prior to the beam hardening correction. A similar voxel counting directly on the measured images without filtering gave a value of 25 ± 5 mg. The discrepancy observed cannot be attributed to incomplete precipitate dissolution during the post mortem analysis, since SEM-EDX analysis showed no residual celestite after acid treatment (Fig. 3). The discrepancy might be due to precipitate microporosity, since a close look at the untreated CT images reveals that the precipitate has variable grey values, corresponding to possible variable densities (see Fig. 4). On the SEM images of the Prec2 and Prec3 precipitates, the aggregates appeared to show micro- or even nanoporosity in between the individual celestite crystals (Fig. 3).
The small increase in porosity before and after the precipitation front observed on the porosity profile suggests variability in the sand compaction. This variability was also observed in a reference diffusion cell (in which no celestite precipitation occurred) measured with CT (Fig. 5). The variation in this reference cell was quasi-periodic, with an amplitude of ±0.03 for an average porosity of 0.35.
Transport parameter quantification by high resolution µCT
The average porosity for the undisturbed porous material in the three Cartesian directions was calculated for ROI-top and ROI-bottom. The pore-morphology model interpreted neighboring voxels as belonging to the same connected group if they were in contact by either their face, edge or vertex. Isolated groups of voxels connected with the pore value were seen as isolated pores. The total porosity was 0.41 ± 0.01, with no isolated pores detected (Fig. 6). This value is in good agreement with the HTO diffusion experiments (0.42 ± 0.09), MIP (0.42 ± 0.02), the µ-CT calculations on the reference diffusion cells (0.40 ± 0.03) and the voxel counting approach on the low (0.44 ± 0.05) and high (0.39 ± 0.03) resolution data sets.
The average effective diffusivity in the three Cartesian directions was calculated with a bulk (Laplace) diffusion approach. This is a continuum mechanics approach using the Laplace equation in pores applying the Neumann boundary conditions at the pore-surface interface. De was found to be (5.03 ± 0.02) × 10−10 m2 s−1 for the ROI-top and ROI-bottom. This value is in good agreement with the HTO tracing experiments (4.48 ± 0.17) × 10−10 m2 s−1).
The average tortuosity factor of the undisturbed material in the three Cartesian directions was calculated for ROI-top and ROI-bottom using the ratio of the pore-morphology model derived porosity and effective diffusivity. It was 1.83 ± 0.02, which is again in good agreement with the geometric factor calculated based on HTO diffusion (1.92 ± 0.44).
In the precipitation zone, the porosity calculated was only slightly lower, as it was an average of the entire ROI-middle, which included portions of the unreacted material. However, the tortuosity increased and the diffusivity decreased significantly. The model indicated that no isolated pores occurred, suggesting that the porosity remaining in the precipitation zone is fully connected.
Prior to reactive transport modeling, the CRUNCHFLOW model was calibrated on the pure HTO-diffusion experiments in which no precipitation occurred. The model was run for each independent experiment, with a mesh size of 0.1 mm. The calculated diffusion curves fit the experimental data well (Fig. 2, solid line) and allowed us to determine the effective diffusion coefficient of HTO, the porosity of the porous media and the cementation factor m in Archie’s Law.
The celestite precipitation rate was then introduced to the model, and an extensive parameter sensitivity study was carried out in which the kinetic regime, mesh size and celestite surface area were varied (see Additional file 1). It was shown that the results of the model depended on mesh size when larger than 0.1 mm, as previously reported [32, 33]. Therefore, the mesh size was limited to 0.1 mm for all further calculations (Additional file 1: Figure S1 and Figure S2). The modeled diffusion and flux curves, mass of precipitate, and porosity profiles were only slightly dependent on the celestite’s kinetic regime (5.13 × 10−8 mol m−2 s−1 for kinetics or 1.0 × 10−3 mol m−2 s−1 for quasi-equilibrium). Furthermore, published values for celestite surface area vary between 0.44 and 72.7 m2 g−1 (e.g. [34, 35]). In addition, two commercially available powders were measured with N2-BET and were found to have surface areas of 0.525 and 0.998 m2 g−1, respectively. The influence on model predictions was tested with the smallest (0.44 m2 g−1), the average (40 m2 g−1) and the largest (72.7 m2 g−1) surface area values. Results are reported in Fig. S3 of the Additional file 1. The model is only slightly dependent on the celestite surface area. The curves for 40 and 72.7 m2 g−1 are very similar. The curve for the 0.44 m2 g−1 case is quite different, as the peak is broader: the precipitation takes about 5 additional days to influence the flux drop. However, the minimum flux value is reached very shortly after the other cases, and the flux value reached is also very similar, as is the specific surface area (see Additional file 1). The main calibration parameters left were the minimum porosity and cementation factor. The minimum porosity is a value that can be set in the model input, such that further precipitation would produce no further porosity reduction.
Summary of the initial experimental conditions for the Sr2+ sorption experiments on purified sea sand
Contact time (h)
S/L ratio (g/L)
0.05 to 0.09
XRF analysis of the raw and cleaned sand
Major elements (wt%)
Minor elements (ppm)
Conditions of through-diffusion experiments in cells Prec1, Prec2 and Prec3 (clogging by celestite precipitation from 0.5 M SrCl2 and Na2SO4 solutions in reservoirs) and celestite masses determined
Background (mol L−1)
HTO tracer (kBq L−1)
Duration of run (days)
Acid selective dissolution (mg)
µCT mass estimation (mg)
21.7 ± 0.4
31.6 ± 0.4
1.74 × 10−4
35.6 ± 0.4
Correspondingly, the final flux predicted by the model of 4.17 Bq m−2 s−1 is about two times lower than the experimentally determined steady-state flux after precipitation of 8.3 ± 0.9 Bq m−2 s−1 observed.
The use of a 1D single continuum model is controversial in reactive and highly heterogeneous media , even in simple systems such as the sea sand presented here. Molins et al.  have shown an overestimation of the dissolution in continuum-scale model simulations by a factor of ~1.9 for well-defined calcite dissolution experiments under advective transport conditions. Similarly, Tartakovski, et al.  have demonstrated that conventional low-resolution advective-dispersion equations produced erroneous results. By an adaptive high-resolution model based on homogeneous and heterogeneous reaction terms in Darcy-scale advection dispersion equations with grid sizes in the mixing zone smaller than the size of the sand grains, these authors were able to qualitatively reproduce the essential features of their experiment. These terms involve transport and mixing indices that account for highly non-uniform pore-scale concentration distributions and highly localized reactions and can serve as an alternative to computationally expensive high-resolution simulations.
Our reactive transport modeling of porosity reduction has shown that the exponential factor m in the Archie’s law formula, which regulates the transport response based on porosity reduction, has the strongest impact on the temporal and spatial evolution of the effective diffusion transport parameters affected by porosity reduction. Furthermore, the exponential factor m deviates strongly in different papers published on the same system as detailed below.
For example Shao et al.  used a cementation factor m of 2.0, while Steefel and Lichtner  used m = 1.0 to interpret the same data for the Maqarin natural analogue site. This suggests the need for extensions to the classical Archie’s law formulation concerning diffusive transport for systems evolving from high to low porosity as demonstrated here by mineral precipitation. First observations of the deviation from the Archie’s law were made by Ioannidis et al.  in porous media with porosity lower than about 0.2, where these authors observed a strong deviation of the formation factor on porosity. By fixing the microporosity between 0.38 and 0.23, theoretical work by Tyagi et al.  were able to show that the nanopores have an increasing effect on the overall transport parameters of the porous media.
We have presented an approach of inter-diffusion of Sr and sulfate forcing supersaturation to describe mineral precipitation-induced porosity reduction in porous media under purely diffusion-controlled conditions. This approach is based on: (1) classical radiotracer inter-diffusion experiments involving Sr and sulfate driving supersaturation and precipitation of celestite; (2) µ-CT coupled to 3D pore-morphology modeling; and (3) post analysis consisting of selective dissolution and SEM-EDX data to quantitatively determine the amount of precipitated celestite. The combination of techniques gave consistent and complementary results.
Celestite porosity reduction of sea sand was incomplete (the initial porosity was decreased by only 25 %), as evidenced by the non-zero steady-state HTO flux at the end of the experiment and the fact that micro-CT mapping indicates an open and connected porosity of 0.30. The celestite precipitate was primarily localized in a ~0.35 mm wide disk-shaped precipitation front surrounded by small areas of disperse precipitation in approximately the middle of the diffusion cell. The heterogeneity created by the precipitation front has a very strong influence on the transport parameters, despite the fact that the remaining porosity is fully connected. This phenomenon was evidenced by the strong effect of the porosity reduction on HTO flux and by the Laplace field calculation performed on the µ-CT 3D volumes.
The 1D single continuum reactive transport model successfully reproduced the experimental HTO through-diffusion data obtained prior to the sea sand porosity reduction with an average porosity of 0.40. However, the same model failed to reproduce the HTO through-diffusion during celestite precipitation, indicating that a single continuum model is not capable of describing the system. When the minimum porosity of 0.30 determined with µ-CT was used, the modeled tracer resembles more a homogeneously higher compacted system, which clearly contradicts the experimental observations. The best fit for the HTO through-diffusion flux curve was obtained using a minimum porosity of 0.01. However, the amount of precipitate obtained from this calculation was about 4.5 times greater than the experimental result. This can partly be explained because the reactive transport model was based on the assumption that all of the porosity (0.40–0.44) was available for precipitation. However, the simulated width of the zone affected by the precipitation was comparable with the experimental observations of diffusion cell Prec1.
From these observations we see that the system was well characterized, but that a simple 1D single continuum model cannot reproduce the pore reduction (clogging) experiments satisfactorily. A more complex dual continuum 2D/3D pore-scale reactive transport model including a modified diffusivity/porosity relationship and possibly complex precipitation kinetics (e.g. heterogeneous nucleation and inhibition processes) may provide a better description. Indeed, the SEM observations seemed to indicate an intricate mix between heterogeneous nucleation and surface growth processes, as was shown by the wide range of shapes and sizes of precipitated celestite. The question of why no more precipitation seems to occur in the column after a very short time (less than 10 days) should also be addressed in terms of inhibition processes.
This more complex model is needed to reproduce the experimental data presented here, and to predict reliably transport behavior in more complex, natural systems on long time scales, unreachable to the experimentalist.
We chose the porous material and diffusion cell’s dimensions based on several feasibility criteria; (1) the material grain size and pore diameters had to be large enough to be well resolved by µ-CT, but not so large as to allow for advective fluid flow, (2) the experiments had to run over a reasonable time period (1–2 months) (this limited the diffusion cell length) and (3) the cell geometry, material and wall thickness had to be suitable for the µ-CT equipment. The materials we chose and the diffusion cells are described in this section.
Porous material characterization
The porous material used was the commercially available Merck purified sea sand with an average particle size of 100–300 µm (product no. 1.07711.5000). The sand was already purified by the manufacturer with acidic treatment and calcination. In another treatment step the sea sand was put in MilliQ water with a volume to mass ratio of 2:1. A small aliquot of 100 µL HNO3 was added to the suspension (pH ~3–4), which was then vigorously agitated for a few seconds. The turbid supernatant was poured out of the vessel a few seconds after agitation, giving enough time for the biggest particles to settle leaving the fines in suspension. This procedure was repeated 7-10 times until the supernatant appeared clear and a stable and a neutral pH was reached. The sand was then dried in an oven at 105 °C for at least 24 h (this sample is later referred as cleaned sand; Table 3). The N2-BET surface area of the powder was equal to 0.55 ± 0.02 m2 g−1. The mineralogical composition of the sand powder was characterized via XRD and XRF. The porosity and petro-physical properties of the compacted sand were determined with mercury intrusion porosimetry (MIP, before the formation of a precipitation front) and µ-CT (before and after the formation of a precipitation front).
X-ray diffraction (XRD) analysis on the bulk powder was performed on a Bruker D8 Advance equipped with a Cu K-alpha radiation source and a Sol-X energy dispersive detector. The analysis confirmed that the sand consists almost solely of silica (>99 % SiO2). The impurities were concentrated in order to be identified. For this purpose, a sample of the sea sand was suspended in MilliQ water and vigorously agitated for a few seconds. The largest particles settled within the first seconds following the agitation, but the supernatant remained turbid for about 3 h afterwards, due to the presence of numerous fine particles. A drop of the supernatant was sampled a few seconds after agitation and left to dry at room temperature on a XRD sample holder. The analysis of this sample revealed the presence of at least two minor mineral phases from the clay and mica groups.
The XRD analysis was confirmed by X-ray Fluorescence (XRF) measurements on the bulk powder (Table 3). The XRF were performed with a MagiXPRO, Fa. Philips, with a Rh anode at 3.6 kW. The trace elements were measured in undiluted powder pellets, and the major elements were measured at 3.2 kW on melt pellets diluted 14 times. The powder mainly consisted of silicon (97.8 wt% SiO2). Aluminum was identified as major impurity with 0.9 wt% Al2O3. Assuming the presence of Al solely in the form of kaolinite, this would correspond to 2.25 wt% kaolinite in the sea sand. The other major elements found were K, Ti, Fe and Ca, in decreasing order of abundance. As depicted from Table 3 the removal of fines from the raw material had no effect on the major element distribution.
Ba, Ce, Cr, Nd, Rb, Sr and Zr were present in trace amounts (10–100 ppm) and did not change during the purification step. Based on the Sr XRF analysis and the total mass of ~6.1 ± 0.3 g sea sand inserted in each column, the Sr amount was estimated to be 0.12 mg which would result, under the assumption of all Sr being bonded in celestite, in a total mass of 0.25 mg SrSO4 present in the total column. Such a small amount would not significantly influence the mass of precipitate obtained from SrSO4 dissolution.
Petrophysical properties and porosity of the compacted material
For the filling, the sand material was thoroughly cleaned with MilliQ water in order to remove impurities present as fine particles. To achieve this, the powder was suspended in twice its volume of MilliQ water and vigorously agitated. The turbid supernatant was discarded, and this procedure repeated until no visible turbidity was left in the supernatant. The sand grains were then dried in an oven at 105 °C for at least 24 h. After drying, the cleaned sand was equilibrated with a solution of either 1 mmol L−1 or 0.5 mol L−1 NaCl background electrolyte for a minimum of 3 days. The wet sand was scooped into the half-opened columns of the diffusion cells, and manually compacted by repeatedly tapping the cell onto the preparation table. Once the columns were filled and the supernatant water decanted, the cells were closed and the reservoirs were filled with the same background electrolyte as was used to saturate the pore space. An approximate porosity could be calculated for the material from the mass of sand compacted in the diffusion cells (about 6.1 ± 0.3 g) and the volume of the sand column (50 mm long by 10 mm diameter). Considering that quartz was the only constituent with a density of 2.65 g cm−3, the porosity of the compacted sand for these cells was calculated at 0.42 ± 0.03.
The porosity of the manually compacted powder was more accurately measured with mercury intrusion porosimetry (MIP) and computed microtomography (µ-CT). MIP was performed at the Federal Institute of Material Research and Testing (BAM, Berlin, Germany) with an Autopore III apparatus (Micromeritics, USA), following a modified approach after Minagawa and coworkers . Two samples were prepared in glass tubes 10 mm long with 7 mm inner diameter, manually compacted in the same manner as the diffusion cells, and closed by two glass frits of 200 µm average pore size (ROBU glass filter, porosity 0, Hattert, Germany).
CT measurements were performed using a CT-alpha 160 X-ray device, (Procon, Germany100 keV X-ray energy at 125 µA tube current, with a 1 mm Al-filter on a diamond-coated target using ten 2 s exposure times on a Hamamatsu flat panel detector) directly on the same column used later for this experiment. The 800 projections obtained for each of the 8 samples (1 with a precipitation front, and 7 without) were reconstructed using the Octopus code (Inside Matters, Ghent, Belgium) in stacks of 2,048 images of 2,048 × 2,048 pixels (2,0483 voxels), oriented perpendicular to the diffusion direction. The reconstructed 3D volumes were segmented with image treatment software supplied by Avizo (FEI Visualization Sciences Group, Bordeaux, France). The segmentation process identifies and labels each component of an imaged system. For example, in the sea sand-celestite system, the result of the segmentation was the creation of a 3D volume with only three attenuation values (or colors): one for the pore space, one for the sand grains, and one for the celestite. For six of these samples, the segmentation was performed both with and without non-local mean filtering (noise reduction) . The segmented volumes were implemented in the commercial pore-morphology model code GeoDICT (Math2Market, Kaiserslautern, Germany), which can estimate porosity, tortuosity, effective diffusivity and other petrophysical properties.
A separate imaging procedure was applied to one of the data sets (Prec1) obtained with µ-CT on the celestite precipitation front. The bright grey values for celestite, by contrast with other components of the porous material, induced artifacts during the beam hardening correction of the high resolution data set (5.89 µm voxel size). Therefore, it was pre-treated by a Matlab script to segment the dense material prior to the beam hardening correction. After the Matlab correction, a non-local mean filter was applied to the volume in 3D, and the pores and sand grains were then segmented as usual with the Avizo code.
Retardation of Sr by adsorption on sand
To characterize the effects of Sr2+ ion retardation by sorption during the celestite precipitation experiments were maintained under diffusion-controlled conditions. Sr sorption on the purified sea sand was measured in batch experiments. Aliquots of the sea-sand powder were suspended in 20 mL of a 1 or 100 mmol L−1 NaCl background electrolyte solution, and spiked with SrCl2 ranging from 1 µ L−1 to 0.1 mmol L−1. We varied the solid to liquid ratio (between 10 and 150), the contact times, and the pH of the solution. All experimental conditions are summarized in Table 2. The suspensions were 0.22 µm membrane-filtered, and the equilibrated solutions were analyzed with ICP-OES.
Initial transport parameters via HTO through-diffusion experiments
The initial porosity and diffusivity of the sand was determined by through-diffusion of a conservative tracer (tritiated water, HTO) in 10 replicates (diffusion cells named Trac1 to Trac10).
Neither aqueous speciation or precipitation/co-precipitation processes, and/or sorption processes are expected to affect the mobility of tritium in soil/water systems . The experiments were performed in homogeneous background electrolyte concentration, so that the HTO was the only species with a concentration gradient. Background stock solutions of NaCl were prepared at 1 mmol L−1 or 0.5 mol L−1, respectively. One reservoir was filled with the spiked solution (approx. 4.20 MBq L−1) while the reservoir at the opposing end was filled with the non-active solution of the same ionic strength. In the case of simple diffusion (no pore space reduction), the spiked reservoir is referred to as the high concentration (HC) reservoir, and the unspiked reservoir is referred to as the low concentration (LC) reservoir. Both were sampled twice a week in order to measure the flux of tracer from the high to the low concentration reservoir. After each sampling, the solutions were renewed to maintain the boundary conditions in the HC and LC reservoirs constant. The concentration of tracer in the LC reservoir did not exceed 1 % of the amount in the HC reservoir, and the difference in concentration in the HC reservoir did not exceed 5 % between two renewals, as recommended by van Loon and coworkers .
HTO through-diffusion during celestite precipitation
Celestite precipitation was induced by counter-diffusion experiments [20, 21], where the constituting cation and anion diffused from the diffusion cell reservoirs at either end of the column. The Sr-reservoir contained the cation at a concentration of 0.5 M SrCl2, and the opposite SO4-reservoir contained the anion in equal concentration of 0.5 M Na2SO4. The celestite supersaturation increased in the center of the column between the two diffusion fronts, until the salt began to precipitate once supersaturation was achieved.
Celestite (SrSO4) was selected as the clogging mineral phase for these experiments because of (1) its density is significantly higher than that of sand, making it possible to distinguish it from the quartz in the sand, (2) its rapid precipitation kinetics, and (3) well established knowledge concerning thermodynamic and mineralogic/crystallographic data available. The initial pH of the stock solutions in the reservoirs was 5.5, and it remained stable throughout the duration of experiments. As for simple diffusion, the Sr-reservoir was spiked with tritiated water as tracer and the reservoirs were sampled and renewed twice a week.
Three experiments were run in three different diffusion cells (Prec1, Prec2, and Prec3). The compacted sand in Prec1 was initially saturated with 1 mM NaCl, while Prec2 and Prec3 had 0.5 M NaCl. In Prec3, sodium silicate (AppliChem, Darmstadt, CAS-No.: 1344-09-8) was added to the reservoirs, in a concentration that corresponds to equilibrium with quartz (1.74 × 10−4 M). This was to ensure homogeneous conditions in the entire diffusion cell including the reservoirs. In Prec1 no radioactive tracer was added, as the experiment was designed to be measured with µ-CT outside the hot laboratory. The initial conditions are summarized in Table 4. All diffusion experiments were conducted with the columns oriented horizontally to eliminate any possibility of gravity-driven flow (Fig. 9).
The activity of the tracer in both reservoirs was measured by liquid scintillation counting (LSC) on a Tri-Carb 3110 TR device (Perkin Elmer, USA). Sample aliquots of 0.1–1 mL were mixed with 10 mL of LSC cocktail (Ultima Gold XR, Perkin Elmer, USA). The alpha- and beta-emission spectra were measured for 30–60 min over the 0–2,000 keV energy range.
Post mortem analysis
The mass of celestite precipitated in the pore space of diffusion cells Prec1, Prec2 and Prec3 was quantified by dissolving the secondary phase and measuring the Sr and Si concentration released in solution. The reservoirs were removed and the compacted sand was pushed out of the column with a piston and sectioned every 5 mm. Each column provided 10 samples of about 0.6 g of material. The solid samples were rinsed with MilliQ water in a vacuum filter. Large volumes of cleaning solution were used compared with pore solution (about 10 mL of MilliQ water for about 0.2 mL of pore solution). The contact time between the samples and the cleaning solution was minimized (a few seconds) to avoid any considerable dissolution of celestite (dissolution rate of celestite at 32 °C, pH 5.7 was determined to be 10−6.36 mol m−2 s−1; ). The samples were then dried in an oven at 105 °C for at least 24 h.
The dissolution protocol used was adapted from a study on celestite dissolution kinetics . The samples were suspended in 100 mL of 0.5 M HCl (pH ~2). The dissolution was performed at 70 °C for 6 h in open 250 mL Duran® glass beaker using an agitator at 500 rpm. The remaining sand suspension was filtered at 0.22 µm, and two aliquots of each extract solution analyzed by ICP-OES (PerkinElmer Optima 4300 DV) for their Sr and Si content. The detection limits were 0.1 and 0.01 mg L−1 for Si and Sr, respectively. The remaining solid (sand matrix) after extraction was rinsed with MilliQ water and dried.
The efficiency of the method and the forms, morphology and purity of the precipitated phase were checked with SEM-EDX (FEI Quanta 650 FEG) on three solid subsamples from diffusion cells Prec1, Prec2 and Prec3 containing a celestite precipitation front before and after reaction.
The experimental data were modeled with CRUNCHFLOW , using the gimrt mode and the thermoddem.dbs database , which applies the Debye-Hückel equation for activity corrections. For the multicomponent diffusion calculations, the diffusion coefficients in water (D w ) at 25 °C for the main primary dissolved species were taken into account. For species with no influence on the precipitation processes, the Dw was set equal to that of the tracer (Dw(HTO) = 2.24 × 10−9 m2 s−1). For the other species, Dw was set at Cl− (Dw(Cl) = 2.03 × 10−9 m2 s−1), Na+ (Dw(Na) = 1.33 × 10−9 m2 s−1), SO4 2− (Dw(SO4) = 1.07 × 10−9 m2 s−1), Sr2+ (Dw(Sr) = 0.794 × 10−9 m2 s−1), and H4SiO4 (Dw(Si) = 0.79 × 10−9 m2 s−1) as implemented in the phreeqc.dat database file. Fundamental problems arise concerning multi-component diffusion because electrolyte solutions must be electrically neutral and as a result of D w differences, the flux of one ion must be linked to the flux of other ions to maintain this electro-neutrality. The complete theoretical development and general description of the so called irreversible processes of multicomponent transport is presented by Miller  where ion fluxes are related to electrochemical potential gradients for all ions in a system. Onsager phenomenological coefficients are defined by relating the flux of each ion to the gradients of all of the thermodynamic forces in the system . The electrical cross-coupling between ions can markedly modify the diffusion process [55, 56]. This phenomenon, however, may be negligible in specific cases such as diffusion of tracers. As the required Onsager phenomenological coefficients are available for only a limited number of systems, assumptions have been made for practical simplifications. Such assumptions include use of infinite dilute ionic mobilities, negligence of off-diagonal Onsager coefficients and use of ideal chemical potential derivatives .
The mineral phases considered in the model were α-quartz as the porous material, and celestite as the precipitating phase. The celestite precipitation rate at 25 °C (experimental conditions 25° ± 3 °C) was set to 5.13 × 10−8 mol m2 s−1, and its activation energy to 8.126 kcal mol−1 . A specific surface area of 40 m2 g−1 was used to calculate the precipitation rate for celestite. This value corresponds to an average of the values found in the literature [59–62].
A Dirichlet boundary condition (constant Sr and SO4 concentration in the reservoirs held throughout the experimental duration) was set. The modeling was carried out with 601 cells 0.1 mm thick in the x direction. Cells 1–50 and 551–600 correspond to the filter sets, and the last cell (601) was dedicated to the virtually forced precipitation of the HTO tracer. The experimentally determined porosity and cementation factor of simple HTO through-diffusion experiments of the sea sand material were directly implemented in the model. The temperature was set at 25 °C, and the pH was 5.5, in agreement with the experimental conditions.
TS, FC, AC and BM conceived and designed the diffusion experiments. AC and SH performed the diffusion experiments. AC, SH, FC and TS analyzed the data. FE and AC performed the µCT measurements and post processed the data. AC, FE and FC contributed to the 1D reactive transport and pore scale modeling work. AC, TS, FC, FE and MK co-wrote the paper. All authors read and approved the final manuscript.
The research conducted at INE leading to the experimental results was performed within the KIT-BRGM bilateral partnership. The results presented in this article were collected during the GL-Transfer Project granted by ANDRA within the framework of the ANDRA/BRGM scientific partnership. We are very grateful for the co-funding received from the European Union’s European Atomic Energy Community’s (Euratom) 7th Framework Programme FP7/2007-2011 under Grant Agreement No. 249624 (CATCLAY project). The authors thank Christophe Tournassat for earlier comments on the manuscript. For data requests please contact directly Thorsten Schäfer (email@example.com).
Compliance with ethical guidelines
Competing interests The authors declare that they have no competing interests.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided 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 Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Ayora C, Caraballo MA, Macias F, Rotting TS, Carrera J, Nieto JM (2013) Acid mine drainage in the Iberian Pyrite Belt: 2. Lessons learned from recent passive remediation experiences. Environ Sci Pollut Res 20:7837–7853View ArticleGoogle Scholar
- Guarracino L, Rotting T, Carrera J (2014) A fractal model to describe the evolution of multiphase flow properties during mineral dissolution. Adv Water Resour 67:78–86View ArticleGoogle Scholar
- Steefel CI, Lasaga AC (1994) A coupled model for transport of multiple chemical-species and kinetic precipitation dissolution reactions with application to reactive flow in single-phase hydrothermal systems. Am J Sci 294:529–592View ArticleGoogle Scholar
- Masciopinto C (2013) Management of aquifer recharge in Lebanon by removing seawater intrusion from coastal aquifers. J Environ Manage 130:306–312View ArticleGoogle Scholar
- Indraratna B, Pathirage PU, Rowe RK, Banasiak L (2014) Coupled hydro-geochemical modelling of a permeable reactive barrier for treating acidic groundwater. Comput Geotech 55:429–439View ArticleGoogle Scholar
- Xie ML, Kolditz O, Moog HC (2011) A geochemical transport model for thermo-hydro-chemical (THC) coupled processes with saline water. Water Resour Res 47:14Google Scholar
- Liu SH, Jacques D, Govaerts J, Wang L (2014) Conceptual model analysis of interaction at a concrete-Boom Clay interface. Phys Chem Earth 70–71:150–159View ArticleGoogle Scholar
- Wasch LJ, Wollenweber J, Tambach TJ (2013) Intentional salt clogging: a novel concept for long-term CO2 sealing. Greenh Gases 3:491–502View ArticleGoogle Scholar
- Surasani VK, Li L, Ajo-Franklin JB, Hubbard C, Hubbard SS, Wu YX (2013) Bioclogging and Permeability Alteration by L-mesenteroides in a Sandstone Reservoir: a Reactive Transport Modeling Study. Energy Fuels 27:6538–6551View ArticleGoogle Scholar
- Xie M, Mayer KU, Claret F, Alt-Epping P, Jacques D, Steefel C et al (2015) Implementation and evaluation of permeability-porosity and tortuosity-porosity relationships linked to mineral dissolution-precipitation. Comput Geosci 19:655–671View ArticleGoogle Scholar
- Archie GE (1942) Electrical resistivity log as an aid in determining some reservoir characteristics. T Am I Min Met Eng 146:54–62Google Scholar
- Navarre-Sitchler A, Steefel CI, Yang L, Tomutsa L, Brantley SL (2009) Evolution of porosity and diffusivity associated with chemical weathering of a basalt clast. J Geophys Res 114:F02016Google Scholar
- Boving TB, Grathwohl P (2001) Tracer diffusion coefficients in sedimentary rocks: correlation to porosity and hydraulic conductivity. J Contam Hydrol 53:85–100View ArticleGoogle Scholar
- Appelo CAJ, Van Loon LR, Wersin P (2010) Multicomponent diffusion of a suite of tracers (HTO, Cl, Br, I, Na, Sr, Cs) in asingle sample of Opalinus Clay. Geochim Cosmochim Ac 74:1201–1219View ArticleGoogle Scholar
- NEA (2009) Natural tracer profiles across argillaceous formations: the CLAYTRAC project, Report NEA No. 6253, NEA-OECD. Nuclear Energy Agency—Organisation for Economic Co-Operation and Development, ParisGoogle Scholar
- Mazurek M, Alt-Epping P, Bath A, Gimmi T, Niklaus Waber H, Buschaert S et al (2011) Natural tracer profiles across argillaceous formations. Appl Geochem 26:1035–1064View ArticleGoogle Scholar
- Gao Z, Hu Q, Liang H (2013) Gas diffusivity in porous media: determination by mercury intrusion porosimetry and correlation to porosity and permeability. J Porous Media 16:607–617View ArticleGoogle Scholar
- Van Loon LR, Glaus MA, Müller W (2007) Anion exclusion effects in compacted bentonites: towards a better understanding of anion diffusion. Appl Geochem 22:2536–2552View ArticleGoogle Scholar
- Vanbavel CHM (1952) Gaseous diffusion and porosity in porous media. Soil Sci 73:91–104View ArticleGoogle Scholar
- Prieto M, FernandezGonzalez A, Putnis A, FernandezDiaz L (1997) Nucleation, growth, and zoning phenomena in crystallizing (Ba, Sr)CO3, Ba(SO4, CrO4), (Ba, Sr)SO4, and (Cd, Ca)CO3 solid solutions from aqueous solutions. Geochim Cosmochim Ac 61:3383–3397View ArticleGoogle Scholar
- Prieto M, Fernandez-Gonzalez A, Martin-Diaz R (2002) Sorption of chromate ions diffusing through barite-hydrogel composites: implications for the fate and transport of chromium in the environment. Geochim Cosmochim Ac 66:783–795View ArticleGoogle Scholar
- Sell K, Enzmann F, Kersten M, Spangenberg E (2013) Microtomographic quantification of hydraulic clay mineral displacement effects during a CO2 Sequestration experiment with saline aquifer sandstone. Environ Sci Technol 47:198–204View ArticleGoogle Scholar
- Molins S, Trebotich D, Yang L, Ajo-Franklin JB, Ligocki TJ, Shen CP, Steefel CI (2014) Pore-Scale Controls on Calcite Dissolution Rates from Flow-through Laboratory and Numerical Experiments. Environ Sci Technol 48:7453–7460View ArticleGoogle Scholar
- Khan F, Enzmann F, Kersten M, Wiegmann A, Steiner K (2012) 3D simulation of the permeability tensor in a soil aggregate on basis of nanotomographic imaging and LBE solver. J Soil Sediment 12:86–96View ArticleGoogle Scholar
- Andra H, Combaret N, Dvorkin J, Glatt E, Han J, Kabel M et al (2013) Digital rock physics benchmarks-Part I: imaging and segmentation. Comput Geosci-Uk 50:25–32View ArticleGoogle Scholar
- Andra H, Combaret N, Dvorkin J, Glatt E, Han J, Kabel M et al (2013) Digital rock physics benchmarks-part II: computing effective properties. Comput Geosci-Uk 50:33–43View ArticleGoogle Scholar
- Saenger EH, Enzmann F, Keehm Y, Steeb H (2011) Digital rock physics: effect of fluid viscosity on effective elastic properties. J Appl Geophys 74:236–241View ArticleGoogle Scholar
- Hilpert M, Miller CT (2001) Pore-morphology-based simulation of drainage in totally wetting porous media. Adv Water Resour 24:243–255View ArticleGoogle Scholar
- Steefel CI, Lasaga AC (1994) A coupled model for transport of multiple chemical species and kinetic precipitation/dissolution reactions with application to reactive flow in single phase hydrothermal systems. Am J Sci 294:529–592View ArticleGoogle Scholar
- EPA US (1999) Understanding Variation in Partition Coefficient, Kd, Values Volume II: Review of Geochemistry and Available Kd Values for Cadmium, Cesium, Chromium, Lead, Plutonium, Radon, Strontium, Thorium, Tritium (3H), and Uranium. EPA 402-R-99-004B, EPA, WashingtonGoogle Scholar
- Chagneau A, Tournassat C, Steefel CI, Bourg IC, Gaboreau S, Esteve I et al (2015) Complete Restriction of 36Cl– Diffusion by Celestite Precipitation in Densely Compacted Illite. Environ Sci Technol Lett 2:139–143View ArticleGoogle Scholar
- Marty NCM, Tournassat C, Burnol A, Giffaut E, Gaucher EC (2009) Influence of reaction kinetics and mesh refinement on the numerical modelling of concrete/clay interactions. J Hydrol 364:58–72View ArticleGoogle Scholar
- Hayek M, Kosakowski G, Churakov S (2011) Exact analytical solutions for a diffusion problem coupled with a precipitation-dissolution reaction and feedback of porosity change. Water Resour Res 47:W07545View ArticleGoogle Scholar
- Suhnel O, Handlirova M (1984) Precipitation of strontium sulfate. Cryst Res Technol 19:477–490View ArticleGoogle Scholar
- Liu S, Yu J, Cheng B, Zhao L, Zhang Q (2005) Facile fabrication of monodispersed mesoporous celestine particles with peanut-shaped morphology. J Cryst Growth 279:461–465View ArticleGoogle Scholar
- Chen CC, Coleman ML, Katz LE (2006) Bridging the gap between macroscopic and spectroscopic studies of metal ion sorption at the oxide/water interface: Sr(II), Co(II), and Pb(II) Sorption to Quartz. Environ Sci Technol 40:142–148View ArticleGoogle Scholar
- Parkman RH, Charnock JM, Livens FR, Vaughan DJ (1998) A study of the interaction of strontium ions in aqueous solution with the surfaces of calcite and kaolinite. Geochim Cosmochim Acta 62:1481–1492View ArticleGoogle Scholar
- Vojtech O, Süssmilch J, Selecka H (1974) Sorption of radioactive strontium on binary-systems of some alkaline-earth sulfates. J Radioanal Chem 21:341–348View ArticleGoogle Scholar
- Tartakovsky AM, Redden G, Lichtner PC, Scheibe TD, Meakin P (2008) Mixing-induced precipitation: experimental study and multiscale numerical analysis. Water Resour Res 44:W06S04View ArticleGoogle Scholar
- Shao H, Kosakowski G, Berner U, Kulik DA, Mäder U, Kolditz O (2013) Reactive transport modeling of the clogging process at Maqarin natural analogue site. Physi Chem Earth Parts A/B/C 64:21–31View ArticleGoogle Scholar
- Steefel CI, Lichtner PC (1998) Multicomponent reactive transport in discrete fractures: II: Infiltration of hyperalkaline groundwater at Maqarin, Jordan, a natural analogue site. J Hydrol 209:200–224View ArticleGoogle Scholar
- Ioannidis MA, Kwiecien MJ, Chatzis I (1997) Electrical conductivity and percolation aspects of statistically homogeneous porous media. Transp Porous Media 29:61–83View ArticleGoogle Scholar
- Tyagi M, Gimmi T, Churakov SV (2013) Multi-scale micro-structure generation strategy for up-scaling transport in clays. Adv Water Resour 59:181–195View ArticleGoogle Scholar
- Saloman EB, Hubbell JH, Scofield JH (1988) X-ray attenuation cross-sections for energies 100 Ev To 100 Kev and elements Z = 1 TO Z = 92. Atom Data Nucl Data Tables 38:1–196View ArticleGoogle Scholar
- Minagawa H, Nishikawa Y, Ikeda I, Miyazaki K, Takahara N, Sakamoto Y, Komai T, Narita H (2008) Characterization of sand sediment by pore size distribution and permeability using proton nuclear magnetic resonance measurement. J Geophys Res 113:B07210Google Scholar
- Buades A, Coll B, Morel JM (2005) A review of image denoising algorithms, with a new one. Multiscale Model. Simul 4:490–530View ArticleGoogle Scholar
- Van Loon LR, Soler JM, Jakob A, Bradbury MH (2003) Effect of confining pressure on the diffusion of HTO, Cl-36(-) and I-125(-) in a layered argillaceous rock (opalinus clay): diffusion perpendicular to the fabric. Appl Geochem 18:1653–1662View ArticleGoogle Scholar
- Van Loon LR, Soler JM, Bradbury MH (2003) Diffusion of HTO, 36Cl-, and 125I- in Opalinus Clay samples from Mont Terri: Effect of confining pressure. J Contam Hydrol 61:73–83View ArticleGoogle Scholar
- Dove PM, Czank CA (1995) Crystal chemical controls on the dissolution kinetics of the isostructural sulfates: Celestite, anglesite, and barite. Geochim Cosmochim Acta 59:1907–1915View ArticleGoogle Scholar
- Aydoğan S, Erdemoğlu M, Aras A, Uçar G, Özkan A (2006) Dissolution kinetics of celestite (SrSO4) in HCl solution with BaCl2. Hydrometallurgy 84:239–246View ArticleGoogle Scholar
- Steefel CI, Appelo CAJ, Arora B, Jacques D, Kalbacher T, Kolditz O et al (2014) Reactive transport codes for subsurface environmental simulation. Comput Geosci. doi:10.1007/s10596-014-9443-x Google Scholar
- Blanc P, Lassin A, Piantone P (2007) BRGM (Orléans, France). http://thermoddem.brgm.fr. Accessed 17 Dec 2014
- Miller DG (1959) Ternary isothermal diffusion and the validity of the onsager reciprocity relations. J Phys Chem 63:570–578View ArticleGoogle Scholar
- Onsager L (1945) Theories and problems of liquid diffusion. Ann NY Acad Sci 46:241–265View ArticleGoogle Scholar
- Ben-Yaakov S (1972) Diffusion of sea water ions—I. Diffusion of sea water into a dilute solution. Geochim Cosmochim Acta 36:1395–1406View ArticleGoogle Scholar
- Lasaga AC (1979) The treatment of multi-component diffusion and ion pairs in diagenetic fluxes. Am J Sci 279:324–346View ArticleGoogle Scholar
- Felmy AR, Weare JH (1991) Calculation of multicomponent ionic diffusion from zero to high concentration: I. The system Na-K-Ca-Mg-Cl-SO4-H2O at 25 °C. Geochim Cosmochim Acta 55:113–131View ArticleGoogle Scholar
- Blanc P, Lassin A, Piantone P, Azaroual M, Jacquemet N, Fabbri A et al (2012) Thermoddem: a geochemical database focused on low temperature water/rock interactions and waste materials. Appl Geochem 27:2107–2116View ArticleGoogle Scholar
- Giffaut E, Grivé M, Blanc P, Vieillard P, Colàs E, Gailhanou H et al (2014) Andra thermodynamic database for performance assessment: ThermoChimie. Appl Geochem 49:225–236View ArticleGoogle Scholar
- Hamdona SK, Hamza SM, Mangood AH (2010) The influence of polyphosphonates on the precipitation of strontium sulfate (celestite) from aqueous solutions. Desalination Water Treat 24:55–60View ArticleGoogle Scholar
- Liu S, Jacques D, Govaerts J, Wang L (2014) Conceptual model analysis of interaction at a concrete-Boom Clay interface. Phys Chem Earth 70–71:150–159View ArticleGoogle Scholar
- Marty NCM, Claret F, Lassin A, Tremosa J, Blanc P, Madé B et al (2015) A database of dissolution and precipitation rates for clay-rocks minerals. Appl Geochem 55:108–118View ArticleGoogle Scholar