# The Dependency of Diffusion Coefficients and Geometric Factor on the Size of the Diffusing Molecule: Observations for Different Clay-Based Materials

Authors: Elke Jacops; Marc Aertsens; Norbert Maes; Christophe Bruggeman; Rudy Swennen; Bernd Krooss; Alexandra Amann-Hildenbrand; Ralf Littke
Journal: Geofluids
 (2017)

---

## Abstract

In order to investigate in more detail the relation between the size of diffusing molecules and their diffusion coefficients (and geometric factors), diffusion experiments with gases of different size and tritiated water (HTO) have been performed on different clayey samples (Boom Clay, Eigenbilzen Sands, Opalinus Clay, Callovo-Oxfordian Clay, and bentonite with different dry densities). We observed that, for unreactive gases in clayey materials, the effective diffusion coefficient varies with the size of the diffusing molecule and this variation can be described by an exponential or a power law function. The variation of the geometric factor can also be described by an exponential function. The observed experimental relations can be used to estimate diffusion coefficients; by measuring experimentally in clay the effective diffusion coefficient of two unreactive dissolved gases with a different size, the diffusion coefficients of other dissolved gases (with a size in between the two measured gases) can be estimated by using the fitted exponential relationship.


---

## Body

## 1. Introduction

Clay-based materials are considered by many countries in their concepts for the safe disposal of high- and intermediate-level radioactive waste, either as the material of choice in the engineered barrier system or because of the choice of argillaceous formations to host the repository. Examples of European countries where argillaceous formations are being explored as potential host formations are Belgium, Switzerland, and France [1–3]. The use of clay-based materials as engineered barriers is studied in, for example, Switzerland [4], France [5], and Sweden [6]. The clays under consideration all have a high sorption capacity for many radionuclides [7–9], a low hydraulic conductivity [10, 11], and interesting self-sealing properties [12, 13].In Belgium, no formal decision has been taken yet on a host formation, but for R&D purposes, the Belgian Radioactive Waste Management Organization (ONDRAF/NIRAS) considers Boom Clay (BC) as a potential host formation for a geological disposal facility. In France, the National Agency for Radioactive Waste Management (Andra) selected the Callovo-Oxfordian Clay Formation in the east of France as a potential host formation [14] and in Switzerland the National Cooperative for the Disposal of Radioactive Waste (NAGRA) proposed the Jurassic Opalinus Clay (OPA) as a host. As engineered barrier, mainly bentonite is studied. A frequently used type of bentonite is MX80, which is studied as backfill material by Andra [5, 15], NAGRA [4], and the Swedish Nuclear Fuel and Waste Management Co. (SKB) [16, 17].In the context of nuclear waste disposal, the transport of dissolved gases in compacted clays is an area, which receives a high amount of interest. First, the production of gas within a geological repository is unavoidable. Mainly anaerobic corrosion of metals will lead to the production of hydrogen. If the rate of gas generation is larger than the diffusive flux into the clay, a free gas phase will form, which might have a negative effect on the performance of the barriers. In order to compute a comprehensive and reliable balance between gas generation versus gas dissipation, correct estimates for gas diffusion coefficients of dissolved gases are essential. Moreover, the produced hydrogen may be converted to other gases like CH4 due to, for example, microbial activity. Thus, also the diffusion coefficient of methane needs to be established.Secondly, naturally occurring noble gases such as He and Ar can act as natural tracers whose profiles can be used to constrain transport properties on the scale of the formation [18–20]. As diffusion is considered to be the dominant transport mechanism, it should be possible to represent these natural tracer profiles by diffusion models, but the availability of reliable diffusion coefficients for He and Ar is limited.Measuring reliable diffusion coefficients of gases is not evident [21] and only a limited amount of data is available in the literature [22]. Recently, diffusion coefficients of dissolved gases (He and Ar) have been measured in Boom Clay, Callovo-Oxfordian Clay, and Opalinus Clay and results have been reported by Jacops et al. [22] and Jacops et al. [23].In case of hydrogen, Jacops et al. [24] have shown that diffusion experiments with hydrogen often suffer from microbial activity: here methanogenic microbes convert H2 into CH4, making the accurate determination of diffusion coefficients impossible [24]. Similar observations have been reported by Vinsot et al. [25].When performing scoping calculations on the diffusive mobility and possible build-up of dissolved gases in a geological repository or in a geological formation, reliable gas diffusion coefficients obtained from laboratory experiments are often not available. Hence, these diffusion coefficients are often estimated from measured values of other species like HTO. In case of hydrogen, the diffusion coefficient for helium is often used as a surrogate. In this approach, it is implicitly assumed that the geometric factor of the formation (a factor which describes the effect of the porous network on diffusion and for which the value obtained for HTO is used) is equal for all other gases and species considered (e.g., [26]).Different approaches for estimating the geometric factor exist: for example, (i) from diffusion experiments (mostly with HTO [26, 27]), (ii) by models (e.g., [28–31]), or (iii) from diffusion simulations on reconstructed clay structures [32, 33]. Archie’s law [28] is a well-known empirical relation relating diffusivity to porosity (but neglecting constrictivity). Variations on Archie’s law are proposed, for example, by Weissberg [34] and Boudreau [29]. Saripalli et al. [30] described a method to calculate tortuosity and constrictivity from the specific surface determined from, for example, N2 adsorption measurements. Chou et al. [31] discussed models to calculate tortuosity for variously saturated soil samples, based on their water content. Both Robinet et al. [32] and Keller et al. [33] calculated geometric factors from diffusion simulations on reconstructed mesostructures of, respectively, the Callovo-Oxfordian and the Opalinus Clay. Their developed approach allows the determination of the geometric factor as a function of mineralogy. All these methods calculate the geometric factor as an intrinsic parameter of a specific material, without taking into account the possible effect of the size of the diffusing molecule.In clay, a relation between the size of the diffusing molecule (expressed by the kinetic diameter) and its diffusion coefficient has been observed [22, 23]. The effect of the size of diffusing molecules on its diffusion coefficient in Callovo-Oxfordian Clay has also been discussed by Dagnelie et al. [35]. For anions, a relation was shown between the anion size (expressed as M1/3 with M the molecular mass) and the ratio of aqueous diffusion coefficients D0/D0HTO. For cations, a similar relation exists but size is expressed as 1/hydrated radius.For mortars with a different sand content, the measured geometric factors for Li+, Cl−, and HTO are relatively similar, leading to the conclusion that the formation factor can be used to determine the order of magnitude of the effective diffusion coefficient of other diffusion species [36]. Analyses of these data indicate a variability of the geometric factor with the size of the diffusing molecule, but this is less pronounced due to the similar size of Li+, Cl−, and HTO. The dependence of the geometric factor on the size of the diffusion molecule is demonstrated in diffusion experiments with He, O2, and Xe in mortars (unpublished data).The main objective of this paper is to investigate in more detail the relation between the size of the diffusing gas molecules and their diffusion coefficients and hence geometric factors in different clayey materials. We also investigate whether this relation can be used to estimate diffusion coefficients of gases, based on their size.This objective is achieved by performing diffusion experiments with gases of different size, on different clayey samples (Boom Clay, Eigenbilzen Sands, Opalinus Clay, Callovo-Oxfordian Clay, and bentonite (Volclay KWK with different dry densities)).


## 2. Materials and Methods

### 2.1. Clay Samples

The Boom Clay is a marine sediment that was deposited in the early Oligocene (Rupelian), 29 to 32 million years ago in the North Sea Basin, at water depths between 50 and 100 m [37]. Shortly after deposition, the accumulated sediment became reducing which is reflected in the common occurrence of framboidal pyrite. The Boom Clay comprises clay minerals (up to 60%, dominated by illite, mixed layered illite-smectite, kaolinite, and traces of chlorite), as well as quartz, K-feldspar, Na-plagioclase, pyrite, and carbonates [38]. It consists of different lithological subunits; more specifically a rhythmic alteration of silty and more clay-rich layers has been observed, as well as the presence of organic- and carbonate-rich layers. Based on these lithological variations, the Boom Clay has been subdivided in four members: the Boeretang Member, the Putte Member, the Terhagen Member, and the Belsele-Waas Member [37]. The sandy unit lying above the Boom Clay is named the Eigenbilzen Sands. The latter consists of dark green, glauconite-rich, clayey, fine-grained to medium-grained sands, with bioturbations [39]. The amount of fine sand increases significantly but the alternation of silty and clayey intervals as observed in the Boom Clay remains [37]. Samples from this formation are in this work described as “clayey sand.”One of the topics under investigation is the effect of variations in the clay and silt/sand content on the diffusion parameters [23]. Therefore, samples have been selected based on their location in the lithostratigraphic column. The clayey samples originate from the Putte or Terhagen Member, whereas the clayey sand samples originate from the Eigenbilzen Sands. Whether the samples are a representative selection or not is investigated by mineralogical analyses, by grain size analyses, and by measuring the hydraulic conductivity. More information on the origin of the used samples and their orientation with respect to bedding plane can be found in Table 1.Table 1
Overview of used samples with the SCK ID, core number, type, and orientation with respect to the bedding.
Material
SCK ID
Core #
Type
Orientation
Boom Clay
K2
ON-Mol-1 84b
Clay
⊥
Boom Clay
K4
ON-Mol-1 127b
Clay
//
Eigenbilzen Sands
K14
ON-Mol-1 36a
Clayey sand
//
Eigenbilzen Sands
K15
ON-Mol-1 37b
Clayey sand
//
Eigenbilzen Sands
K16
ON-Mol-1 35b
Clayey sand
⊥
Eigenbilzen Sands
K17
ON-Mol-1 39b
Clayey sand
⊥
Volclay KWK
Bentonite 1.4
Dry density 1.4 g/cm3
Volclay KWK
Bentonite 1.6
Dry density 1.6 g/cm3
Callovo-Oxfordian Clay
COX
EST 49109
Claystone
⊥
Opalinus Clay
OPA
Schlattingen 860.32 m depth
Claystone
⊥All samples from the Boom Clay and Eigenbilzen Sands used in this study are taken from the ON-Mol1 borehole which was drilled in 1997, in the town of Mol, in the northeast of Belgium (Lambert coordinatesX (m) 200191.278 andY (m) 211651.761). The samples have been stored under anoxic conditions (vacuum packed in aluminium-coated PE-foil) at 4°C. Prior to their use, they were visually inspected for signs of oxidation which is the main quality threat and the main indicator for perturbations and deteriorations, which might have occurred during storage.The Callovo-Oxfordian Clay (COX) is an indurated, Middle Jurassic mudstone from marine origin, which was deposited 150–160 million years ago. It has been intensively studied by Andra on a 250 km2 area. In this area, Andra has performed several drilling campaigns and constructed an Underground Research Laboratory (URL) in Bure, at a depth of 490 m in the COX layer. The thickness of the COX formation is about 135 m at the URL. The COX consists mainly of clay minerals (illite, interstratified smectite/illite, and others), calcite, dolomite/ankerite, quartz, feldspars, and minor amounts of accessory phases, such as pyrite [40]. The used sample (EST 49109; depth, 478.52 m) was taken from the OHZ6560 borehole, which has been drilled from a technical gallery at the URL. The sample axis is oriented perpendicular to the bedding plane.The Opalinus Clay is a fine-grained sedimentary rock, which was deposited 172 million years ago. The latter consists mainly of clay minerals (illite, illite/smectite, kaolinite, and others), calcite, dolomite/ankerite, quartz, feldspars, pyrite, and organic matter. Research on the Opalinus Clay is mainly performed in the Mont Terri Underground Research Lab, located in the canton Jura. The studied sample was, however, taken from the Schlattingen borehole (Northeast of Switzerland, canton Thurgau), at a depth of 860.32 m below surface. As the burial history for the Mont Terri site and Schlattingen is different, with deeper burial in Schlattingen [41], the physical properties of the Opalinus Clay are different at both locations. The sample axis is oriented perpendicular to the bedding plane.The Volclay KWK (a MX 80 type bentonite) is a fine-grained sodium bentonite with montmorillonite as the main component [42]. This type of bentonite was selected because it was characterized within the NF-PRO project [43] and a batch was readily available at SCK•CEN. Besides, a similar type of bentonite was used as a sealing material within the PRACLAY heater experiment which is performed in the HADES URL [44]. Two different bentonite samples were used, each compacted at different dry densities (1.4 and 1.6 g/cm3).


### 2.2. Sample Preparation

For the experimental setup of the design of the diffusion cell, see Jacops et al. [21]. The Boom Clay sample (diameter 80 mm, length 30 mm) was loaded into the diffusion cell with a hydraulic press. Confinement (constant volume) was achieved by sealing the diffusion cell (welding). Due to the high plasticity and swelling capacity [12] perfect sealing between clay and cell was achieved. Porous stainless steel filter plates (diameter 80 mm; thickness 2 mm; porosity 40%; GKN Sinter Metals) were used at both sides. This provided a good contact between the gas-saturated water and the clay.As Callovo-Oxfordian Clay and Opalinus Clay contain lower proportions of swelling clay minerals, their swelling capacity/plasticity is therefore much lower compared to the Boom Clay. Hence, the samples had to be embedded in a resin (Sikadur 52 Injection Normal) in order to seal the interface clay sample-cell. The embedded samples (80 mm diameter and 25 mm height for Opalinus Clay and 70 mm diameter and 30 mm height for the Callovo-Oxfordian Clay) were placed between an upper and lower flange, and subsequently the flanges were welded to the diffusion cell (confinement with constant volume). As the flanges were provided with a circulation loop, good contact between the water containing the dissolved gas and the clay sample was ensured. More information on the preparation of these samples can be found in Jacops et al. [23] and Jacops et al. [22].The Volclay KWK was compacted to two different dry densities, that is, 1.4 and 1.6 g/cm3. Prior to the compaction, the water content of the batch was measured by drying at 105°C for at least 24 hours. The correct amount of bentonite was compacted into a cylinder of 20 by 20 mm. Next, the compacted samples were saturated with a synthetic pore water of 0.05 M NaClO4.


### 2.3. Sample Characterization

For mineralogical analyses samples were grinded first with mortar and pestle and afterwards (after adding ethanol) in a McCrone Micronizing Mill. After drying, the residue was crushed again with mortar and pestle and finally it was packed in the sample holder [45]. Mineralogy was determined by X-ray diffraction (XRD) analyses as described by [46]. The mineralogy of the Opalinus Clay sample was determined by the University of Bern [47].The specific surface area was measured by nitrogen adsorption at 77 K in a TriStar 3020 (Micromeritics), using the Brunauer-Emmett-Teller theory. Prior to the measurement, ca. 3 g sample was degassed at 110°C. Next, the pressure of N2 was gradually increased and the amount of adsorbed gas was measured as a function of the relative pressure, which plots the adsorption isotherm.Once the diffusion experiments ended, the porosity was determined by drying the samples at 105°C until the mass was stable. Based on the mass loss during drying and the density, the porosity could be calculated [48].


### 2.4. Hydraulic Conductivity

Hydraulic conductivity was measured prior to the diffusion experiment. The used technique is based on Darcy's law and applies a constant pressure difference [49]. The sample was considered to be saturated when water outflow and inflow were equal in at least five successive measurements.


### 2.5. Diffusion of Dissolved Gases

The methodology of through-diffusion experiments has been described in detail in Jacops et al. [21]. Only the main aspects are repeated here. The basis is the double through-diffusion technique [50]. Both sides of the test core were connected to a water vessel filled with circa 500 mL oxygen-free synthetic pore water and pressurised with a gas (circa 500 mL gas at 1 MPa). The composition of the pore water used for Boom Clay, and Eigenbilzen Sands, Callovo-Oxfordian Clay, and Opalinus Clay can be found in Jacops et al. [22] and for the bentonite samples, degassed demineralised water with 0.05 M NaClO4 was used. Both vessels were filled with a different gas, but the total pressure was the same. In this way, no advective flux could occur and the clay sample remained fully water saturated. According to Henry's law, equilibrium was obtained between the free gas in the gas phase and the dissolved gas in the water. The water at both sides was then circulated over the filters, which were in contact with the clay core, allowing the dissolved gases to diffuse through the clay core, towards the reservoir on the opposing side. The change in gas composition in both reservoirs was measured as a function of time by gas chromatography.Sampling of the gas phase was performed on a regular basis (generally once per week) until 10 data points were obtained in the regime of approximately constant outlet flux of the diffusion process. The gas composition was analysed with a CP4900 micro GC (equipped with a Molsieve 5A and a Pora Plot U column and TCD detectors, Agilent, USA) or a CG4 Compact GC (equipped with a RT-Qbond column, a Molsieve 5A column and TCD detectors, Interscience, The Netherlands). Both GCs were operated with the EZChrom CDS software.The experiment was performed in a temperature-controlled room (21 ± 2°C).The diffusion experiments with several gases on the same sample were performed consecutively: the content of the inlet and outlet vessel was replaced each time while the sample remained in the same position.


### 2.6. Diffusion of HTO

In order to allow comparison of our results with available results in literature, diffusion experiments with HTO were performed with a setup similar to the one used to measure through-diffusion of dissolved gases. A more detailed description can be found in Jacops et al. [23]. The water vessels were filled with 250 mL synthetic pore water and 250 mL Ar at 1 MPa. A HTO spiking solution of approximately 50 kBq was added in the high concentration compartment. Both the decrease in activity in the high concentration compartment and the increase in the low concentration compartment were measured by taking at each side two samples of 1 mL. The water in both compartments was not replaced during the experiment. Sampling occurred regularly (once or twice per week) until 10 data points were obtained in the regime of approximately constant outlet flux. Samples were analysed for their radioactivity with a Packard TRI-CARB 2100TR liquid scintillation counter.


## 2.1. Clay Samples

The Boom Clay is a marine sediment that was deposited in the early Oligocene (Rupelian), 29 to 32 million years ago in the North Sea Basin, at water depths between 50 and 100 m [37]. Shortly after deposition, the accumulated sediment became reducing which is reflected in the common occurrence of framboidal pyrite. The Boom Clay comprises clay minerals (up to 60%, dominated by illite, mixed layered illite-smectite, kaolinite, and traces of chlorite), as well as quartz, K-feldspar, Na-plagioclase, pyrite, and carbonates [38]. It consists of different lithological subunits; more specifically a rhythmic alteration of silty and more clay-rich layers has been observed, as well as the presence of organic- and carbonate-rich layers. Based on these lithological variations, the Boom Clay has been subdivided in four members: the Boeretang Member, the Putte Member, the Terhagen Member, and the Belsele-Waas Member [37]. The sandy unit lying above the Boom Clay is named the Eigenbilzen Sands. The latter consists of dark green, glauconite-rich, clayey, fine-grained to medium-grained sands, with bioturbations [39]. The amount of fine sand increases significantly but the alternation of silty and clayey intervals as observed in the Boom Clay remains [37]. Samples from this formation are in this work described as “clayey sand.”One of the topics under investigation is the effect of variations in the clay and silt/sand content on the diffusion parameters [23]. Therefore, samples have been selected based on their location in the lithostratigraphic column. The clayey samples originate from the Putte or Terhagen Member, whereas the clayey sand samples originate from the Eigenbilzen Sands. Whether the samples are a representative selection or not is investigated by mineralogical analyses, by grain size analyses, and by measuring the hydraulic conductivity. More information on the origin of the used samples and their orientation with respect to bedding plane can be found in Table 1.Table 1
Overview of used samples with the SCK ID, core number, type, and orientation with respect to the bedding.
Material
SCK ID
Core #
Type
Orientation
Boom Clay
K2
ON-Mol-1 84b
Clay
⊥
Boom Clay
K4
ON-Mol-1 127b
Clay
//
Eigenbilzen Sands
K14
ON-Mol-1 36a
Clayey sand
//
Eigenbilzen Sands
K15
ON-Mol-1 37b
Clayey sand
//
Eigenbilzen Sands
K16
ON-Mol-1 35b
Clayey sand
⊥
Eigenbilzen Sands
K17
ON-Mol-1 39b
Clayey sand
⊥
Volclay KWK
Bentonite 1.4
Dry density 1.4 g/cm3
Volclay KWK
Bentonite 1.6
Dry density 1.6 g/cm3
Callovo-Oxfordian Clay
COX
EST 49109
Claystone
⊥
Opalinus Clay
OPA
Schlattingen 860.32 m depth
Claystone
⊥All samples from the Boom Clay and Eigenbilzen Sands used in this study are taken from the ON-Mol1 borehole which was drilled in 1997, in the town of Mol, in the northeast of Belgium (Lambert coordinatesX (m) 200191.278 andY (m) 211651.761). The samples have been stored under anoxic conditions (vacuum packed in aluminium-coated PE-foil) at 4°C. Prior to their use, they were visually inspected for signs of oxidation which is the main quality threat and the main indicator for perturbations and deteriorations, which might have occurred during storage.The Callovo-Oxfordian Clay (COX) is an indurated, Middle Jurassic mudstone from marine origin, which was deposited 150–160 million years ago. It has been intensively studied by Andra on a 250 km2 area. In this area, Andra has performed several drilling campaigns and constructed an Underground Research Laboratory (URL) in Bure, at a depth of 490 m in the COX layer. The thickness of the COX formation is about 135 m at the URL. The COX consists mainly of clay minerals (illite, interstratified smectite/illite, and others), calcite, dolomite/ankerite, quartz, feldspars, and minor amounts of accessory phases, such as pyrite [40]. The used sample (EST 49109; depth, 478.52 m) was taken from the OHZ6560 borehole, which has been drilled from a technical gallery at the URL. The sample axis is oriented perpendicular to the bedding plane.The Opalinus Clay is a fine-grained sedimentary rock, which was deposited 172 million years ago. The latter consists mainly of clay minerals (illite, illite/smectite, kaolinite, and others), calcite, dolomite/ankerite, quartz, feldspars, pyrite, and organic matter. Research on the Opalinus Clay is mainly performed in the Mont Terri Underground Research Lab, located in the canton Jura. The studied sample was, however, taken from the Schlattingen borehole (Northeast of Switzerland, canton Thurgau), at a depth of 860.32 m below surface. As the burial history for the Mont Terri site and Schlattingen is different, with deeper burial in Schlattingen [41], the physical properties of the Opalinus Clay are different at both locations. The sample axis is oriented perpendicular to the bedding plane.The Volclay KWK (a MX 80 type bentonite) is a fine-grained sodium bentonite with montmorillonite as the main component [42]. This type of bentonite was selected because it was characterized within the NF-PRO project [43] and a batch was readily available at SCK•CEN. Besides, a similar type of bentonite was used as a sealing material within the PRACLAY heater experiment which is performed in the HADES URL [44]. Two different bentonite samples were used, each compacted at different dry densities (1.4 and 1.6 g/cm3).


## 2.2. Sample Preparation

For the experimental setup of the design of the diffusion cell, see Jacops et al. [21]. The Boom Clay sample (diameter 80 mm, length 30 mm) was loaded into the diffusion cell with a hydraulic press. Confinement (constant volume) was achieved by sealing the diffusion cell (welding). Due to the high plasticity and swelling capacity [12] perfect sealing between clay and cell was achieved. Porous stainless steel filter plates (diameter 80 mm; thickness 2 mm; porosity 40%; GKN Sinter Metals) were used at both sides. This provided a good contact between the gas-saturated water and the clay.As Callovo-Oxfordian Clay and Opalinus Clay contain lower proportions of swelling clay minerals, their swelling capacity/plasticity is therefore much lower compared to the Boom Clay. Hence, the samples had to be embedded in a resin (Sikadur 52 Injection Normal) in order to seal the interface clay sample-cell. The embedded samples (80 mm diameter and 25 mm height for Opalinus Clay and 70 mm diameter and 30 mm height for the Callovo-Oxfordian Clay) were placed between an upper and lower flange, and subsequently the flanges were welded to the diffusion cell (confinement with constant volume). As the flanges were provided with a circulation loop, good contact between the water containing the dissolved gas and the clay sample was ensured. More information on the preparation of these samples can be found in Jacops et al. [23] and Jacops et al. [22].The Volclay KWK was compacted to two different dry densities, that is, 1.4 and 1.6 g/cm3. Prior to the compaction, the water content of the batch was measured by drying at 105°C for at least 24 hours. The correct amount of bentonite was compacted into a cylinder of 20 by 20 mm. Next, the compacted samples were saturated with a synthetic pore water of 0.05 M NaClO4.


## 2.3. Sample Characterization

For mineralogical analyses samples were grinded first with mortar and pestle and afterwards (after adding ethanol) in a McCrone Micronizing Mill. After drying, the residue was crushed again with mortar and pestle and finally it was packed in the sample holder [45]. Mineralogy was determined by X-ray diffraction (XRD) analyses as described by [46]. The mineralogy of the Opalinus Clay sample was determined by the University of Bern [47].The specific surface area was measured by nitrogen adsorption at 77 K in a TriStar 3020 (Micromeritics), using the Brunauer-Emmett-Teller theory. Prior to the measurement, ca. 3 g sample was degassed at 110°C. Next, the pressure of N2 was gradually increased and the amount of adsorbed gas was measured as a function of the relative pressure, which plots the adsorption isotherm.Once the diffusion experiments ended, the porosity was determined by drying the samples at 105°C until the mass was stable. Based on the mass loss during drying and the density, the porosity could be calculated [48].


## 2.4. Hydraulic Conductivity

Hydraulic conductivity was measured prior to the diffusion experiment. The used technique is based on Darcy's law and applies a constant pressure difference [49]. The sample was considered to be saturated when water outflow and inflow were equal in at least five successive measurements.


## 2.5. Diffusion of Dissolved Gases

The methodology of through-diffusion experiments has been described in detail in Jacops et al. [21]. Only the main aspects are repeated here. The basis is the double through-diffusion technique [50]. Both sides of the test core were connected to a water vessel filled with circa 500 mL oxygen-free synthetic pore water and pressurised with a gas (circa 500 mL gas at 1 MPa). The composition of the pore water used for Boom Clay, and Eigenbilzen Sands, Callovo-Oxfordian Clay, and Opalinus Clay can be found in Jacops et al. [22] and for the bentonite samples, degassed demineralised water with 0.05 M NaClO4 was used. Both vessels were filled with a different gas, but the total pressure was the same. In this way, no advective flux could occur and the clay sample remained fully water saturated. According to Henry's law, equilibrium was obtained between the free gas in the gas phase and the dissolved gas in the water. The water at both sides was then circulated over the filters, which were in contact with the clay core, allowing the dissolved gases to diffuse through the clay core, towards the reservoir on the opposing side. The change in gas composition in both reservoirs was measured as a function of time by gas chromatography.Sampling of the gas phase was performed on a regular basis (generally once per week) until 10 data points were obtained in the regime of approximately constant outlet flux of the diffusion process. The gas composition was analysed with a CP4900 micro GC (equipped with a Molsieve 5A and a Pora Plot U column and TCD detectors, Agilent, USA) or a CG4 Compact GC (equipped with a RT-Qbond column, a Molsieve 5A column and TCD detectors, Interscience, The Netherlands). Both GCs were operated with the EZChrom CDS software.The experiment was performed in a temperature-controlled room (21 ± 2°C).The diffusion experiments with several gases on the same sample were performed consecutively: the content of the inlet and outlet vessel was replaced each time while the sample remained in the same position.


## 2.6. Diffusion of HTO

In order to allow comparison of our results with available results in literature, diffusion experiments with HTO were performed with a setup similar to the one used to measure through-diffusion of dissolved gases. A more detailed description can be found in Jacops et al. [23]. The water vessels were filled with 250 mL synthetic pore water and 250 mL Ar at 1 MPa. A HTO spiking solution of approximately 50 kBq was added in the high concentration compartment. Both the decrease in activity in the high concentration compartment and the increase in the low concentration compartment were measured by taking at each side two samples of 1 mL. The water in both compartments was not replaced during the experiment. Sampling occurred regularly (once or twice per week) until 10 data points were obtained in the regime of approximately constant outlet flux. Samples were analysed for their radioactivity with a Packard TRI-CARB 2100TR liquid scintillation counter.


## 3. Theory and Modelling

### 3.1. Diffusion Theory

The theoretical aspects of diffusion are discussed in detail in Jacops et al. [23] and only some aspects are discussed below.In general, two transport parameters can be obtained from diffusion experiments: the apparent diffusion coefficientDapp (m2/s) and the capacity factor ηR, being the product of the accessible porosity η (dimensionless) and the retardation factorR (dimensionless). The capacity factor is the ratio of the tracer concentration Cb in the bulk sample and the corresponding concentration C in the pore fluid (solution): ηR=Cb/C.From these two basic parameters, one can calculate the effective diffusion coefficientDeff (m2/s) (1)Deff=ηRDappand, assuming the accessible porosity equals the total porosity ηtot, the pore diffusion coefficient Dp (m2/s)(2)Dp=RDapp.For unretarded species without ion exclusion (e.g., HTO, noble gases), the retardation factor is generally taken equal to one (R=1), leading to Dapp=Dp and a capacity factor equal to the total porosity ηtot.The geometric factorG relates the diffusion coefficient in a porous medium to the corresponding diffusion coefficient D0 (m2/s) in water and is defined by(3)G=D0Dp=D0RDapp=ηD0Deff.Although the geometric factor is a black box factor, it is often split into two dimensionless factors: tortuosity τ and constrictivity δ,(4)G=τδ.In this expression, used in the present paper and Grathwohl [51], tortuosity is defined as the square of the ratio of the effective path length (along the path) of a diffusing component to the shortest distance (end-to-end distance) of that path. Constrictivity takes into account the reduction of the effective diffusion coefficient due to a drag by the pore wall and depends on the ratio of the solute diameter to the pore diameter.Grathwohl [51] lists some empirical expressions for constrictivity, for example,(5)δ=exp⁡-4.6σRpore,with σ the size (diameter) (m) of the diffusing molecule and Rpore the pore width (diameter) (m). Substituting (5) into (4) leads to(6)G=τexp⁡4.6σRpore.Evidently, reducing the pore size distribution of clay to a single pore size Rpore is an enormous simplification as clays have in reality a wide pore size distribution. In expression (6), like in reality, the influence of the size of the diffusing molecule σ depends on the ratio of molecular size to pore width. Although empirical (and certainly not valid over the entire σ/Rpore range), expression (6) also takes into account that sufficiently large molecules (σ→∞) can no longer diffuse through the clay anymore (G→∞ corresponding to a zero effective/apparent/pore diffusion coefficient). This has been observed, for example, for the transport of colloids (natural organic matter) in Boom Clay [52], where small sized diffusing molecules can move through the clay while larger ones cannot. This is related to a percolation transition [53]. If the size of the diffusing molecule rises, fewer pores are available for the transport (diffusion) of that molecule. In a first stage, this will lead to longer travel paths (“tortuosity”) since only pores with a pore size larger than the size of the diffusing molecule are available. When still increasing the molecular size, there is no longer a connected path (of pores available for transport) between the larger pores and the diffusion coefficient becomes zero. In the remainder of this paper, expression (6) is used as a simple empirical expression, trying to capture the evolution of theG-factor as a function of size of the diffusing molecule (but clearly wrong when σ/Rpore≫1), which does not imply the approval of expression (4) and the corresponding definition of tortuosity and constrictivity.


### 3.2. Single Pore Size Hydraulic Conductivity Model

Due to the single pore size approximation in theG-factor, also a single pore size model is presented for advective flow. The flow Qcore (m3/s) out of a cylindrical clay core with radius Rcore (m) and length Lcore (m) is provided by Darcy's law,(7)Qcore=πRcore2ΔpLcoreKwith Δp being the pressure difference (in meter water column) over the core.Poiseuille's law gives the flowQpore (m3/s) through a cylindrical pore with radius Rpore (m)(8)Qpore=1μπRpore48ΔpLcorewith μ being the viscosity (kg/(m s)) (a value μ=0.7978 mPa s=8.135×10-8 m s was used (Edge Engineers)).IfCpore is the number of pores per unit surface (m−2), the flow Qcore also equals(9)Qcore=πRcore2CporeQporeand the core porosity ηtot is given by(10)ηtot=πRpore2Cpore.Combining expressions (7)–(10) leads to the next expression for the hydraulic conductivity(11)K=18μRpore2ηtot=π8μRpore4Cporeshowing that, for two clays with the same porosity, the hydraulic conductivity is the highest for the clay with the largest pore size. It is also clear that, in case of a pore size distribution, a major contribution to the hydraulic conductivity comes from the largest pores.


### 3.3. Size of the Diffusing Molecules

The size of a diffusing molecule can be characterized in several ways. As it would be not appropriate to use values which have been determined from diffusion measurements (hence leading to circular reasoning), some extra clarification on the selection of the kinetic diameters is given in this paragraph. The most widely used measure for the size of a (small) gas molecule is the kinetic diameter(σ) [54]. In this approach, the interaction potential between two (gas) molecules is assumed a Lennard-Jones VLJ(r) potential (with r distance between the centre of two molecules) [55].The potentialVLJ(r) is zero at r=σ and the distance σ (called the kinetic diameter) is considered representative for the size of the gas molecule. For the kinetic diameters used in this work, the parameters of the Lennard-Jones potential are derived from viscosity measurements at different temperatures.However, the selected kinetic diameters are representative for dilute gas mixtures and are thus only an approximation for gases dissolved in water. More information on this topic can be found in [23].


### 3.4. Size Dependency of the Diffusion Coefficient of Gases in Water and Clay

Empirically, the gas diffusion coefficientsD0 in water were fitted as a function of the kinetic diameter σ (expressed in 10-10m) by an exponential equation (12)D0=D0,expexp⁡-n0,expσleading to D0,exp=5.5±2.3×10-8 m2/s and n0,exp=(8.8±1.2)×109 m−1 [23]. Another possibility is a power law fit (13)D0=D0,powσ-n0,powleading to D0,pow=1.1±0.46×10-7 m2/s and n0,pow=3.0±0.4 [23]. Because expressions (12) and (13) do not take into account details like molecular geometry or polarity of a molecule, they only describe a global evolution without predicting correctly every single point (like, e.g., Archie's law).Similar to the diffusion coefficientD0 in free water, for gases (and HTO) trend lines for the dependence of the effective diffusion coefficient Deff on the kinetic diameter σ are fitted by assuming similar expressions as (12) and (13) [23]: (14)Deff=Deff,expexp⁡-nexpσ,(15)Deff=Deff,powσ-npow.Substituting (12) and (14) into (1) provides for unretarded tracers without ion exclusion (the retardation factor R=1) an expression for the geometric factor:(16)G=ηD0,expDeff,expexp⁡nexp-n0,expσ.This expression reduces to the Gratwohl expression (6) by putting(17)τ=ηD0,expDeff,exp(18)Rpore=4.6nexp-n0,exp,where a positive pore size Rpore requires that nexp>n0,exp.A similar combination as (16) can be made with expressions (13) and (15) leading to (and assuming npow>=n0,pow)(19)G=ηD0,powDeff,powσnpow-n0,pow.An expression similar to (17) can be defined for the power law (19), but no characteristic distance similar to Rpore exists. Also note that with the power law, G(σ=0)=0 (assuming npow>=n0,pow) instead of G(σ=0)=τ for the exponential relation (16).


### 3.5. Modelling

A detailed description on the modelling approach can be found in Jacops et al. [23].The diffusion experiments are modelled by fitting the solutions of the diffusion equation with the appropriate boundary and initial conditions. For both the gas diffusion and the HTO diffusion experiments, the diffusion equation is solved by COMSOL coupled with MATLAB for optimization.A through-diffusion experiment allows fitting both the apparent diffusion coefficientDapp and the capacity factor ηR, which are in general both unknown. Besides, also the initial concentration, which is known but prone to measurement errors, is fitted and this is described as the 3-parameter fit. Diffusion coefficients for HTO, reported in this paper have been obtained with this 3-parameter fit. An alternative fit strategy is also possible: assumingR = 1, only Dapp is fitted and ηR is a constant equal to ηtot (measured on either a sample or an average value from literature). In some gas through-diffusion experiments, the (quasi)stationary state is reached very fast, leading to a nearly total correlation between the capacity factor, the apparent diffusion coefficient, and the initial concentration. In this case, the one parameter fit is the only option and only the effective diffusion coefficient can be determined. Therefore, all diffusion coefficients for dissolved gases are obtained from a one-parameter fit. The one- and three-parameter fits can lead to small differences in the optimal values [23].The water content of a core was measured at the end of the diffusion experiments, and if this information was not available, a reference value obtained from literature was used (see also Table3).As discussed above, confining filters are used in the diffusion experiments. As discussed by Glaus et al. [56], Birgersson and Karnland [57], and Aertsens et al. [48] the filters might have an influence on the results as they contribute to the diffusive resistance of the experiment. Therefore, two fitting strategies are possible: one assuming the same transport parameter values for clay core and filter (case “Filter = Clay” – “F = C”) and one with fixed (=the measured but not very accurate values) filter parameter values (case “Fixed Filter,” “FF”). However, the only element for which the values of the transport parameters in the filters are measured is HTO [48]. For Boom Clay oriented perpendicular to the bedding plane, the values of the transport parameters are roughly similar to those of the filter so it can be a good approximation to assume the same transport parameters for the filter as for the clay core (case F = C). However, also Boom Clay cores oriented parallel to the bedding plane and other types of cores are investigated and, in these cases, the diffusion coefficient of the filter differs considerably from the HTO diffusion coefficient in the samples. Therefore, the most correct approach is to use in the model the known HTO filter diffusion coefficient (case fixed filter (FF)).Diffusion coefficients of dissolved gases in the filters are not available and therefore only the value of HTO could be used as an approximation. As this would introduce a large uncertainty into the model and given the large length of the samples (35 mm) compared to the length of the filters (2 × 2 mm), the effect of the filter on diffusion is not taken into account in the experiments with dissolved gases.


## 3.1. Diffusion Theory

The theoretical aspects of diffusion are discussed in detail in Jacops et al. [23] and only some aspects are discussed below.In general, two transport parameters can be obtained from diffusion experiments: the apparent diffusion coefficientDapp (m2/s) and the capacity factor ηR, being the product of the accessible porosity η (dimensionless) and the retardation factorR (dimensionless). The capacity factor is the ratio of the tracer concentration Cb in the bulk sample and the corresponding concentration C in the pore fluid (solution): ηR=Cb/C.From these two basic parameters, one can calculate the effective diffusion coefficientDeff (m2/s) (1)Deff=ηRDappand, assuming the accessible porosity equals the total porosity ηtot, the pore diffusion coefficient Dp (m2/s)(2)Dp=RDapp.For unretarded species without ion exclusion (e.g., HTO, noble gases), the retardation factor is generally taken equal to one (R=1), leading to Dapp=Dp and a capacity factor equal to the total porosity ηtot.The geometric factorG relates the diffusion coefficient in a porous medium to the corresponding diffusion coefficient D0 (m2/s) in water and is defined by(3)G=D0Dp=D0RDapp=ηD0Deff.Although the geometric factor is a black box factor, it is often split into two dimensionless factors: tortuosity τ and constrictivity δ,(4)G=τδ.In this expression, used in the present paper and Grathwohl [51], tortuosity is defined as the square of the ratio of the effective path length (along the path) of a diffusing component to the shortest distance (end-to-end distance) of that path. Constrictivity takes into account the reduction of the effective diffusion coefficient due to a drag by the pore wall and depends on the ratio of the solute diameter to the pore diameter.Grathwohl [51] lists some empirical expressions for constrictivity, for example,(5)δ=exp⁡-4.6σRpore,with σ the size (diameter) (m) of the diffusing molecule and Rpore the pore width (diameter) (m). Substituting (5) into (4) leads to(6)G=τexp⁡4.6σRpore.Evidently, reducing the pore size distribution of clay to a single pore size Rpore is an enormous simplification as clays have in reality a wide pore size distribution. In expression (6), like in reality, the influence of the size of the diffusing molecule σ depends on the ratio of molecular size to pore width. Although empirical (and certainly not valid over the entire σ/Rpore range), expression (6) also takes into account that sufficiently large molecules (σ→∞) can no longer diffuse through the clay anymore (G→∞ corresponding to a zero effective/apparent/pore diffusion coefficient). This has been observed, for example, for the transport of colloids (natural organic matter) in Boom Clay [52], where small sized diffusing molecules can move through the clay while larger ones cannot. This is related to a percolation transition [53]. If the size of the diffusing molecule rises, fewer pores are available for the transport (diffusion) of that molecule. In a first stage, this will lead to longer travel paths (“tortuosity”) since only pores with a pore size larger than the size of the diffusing molecule are available. When still increasing the molecular size, there is no longer a connected path (of pores available for transport) between the larger pores and the diffusion coefficient becomes zero. In the remainder of this paper, expression (6) is used as a simple empirical expression, trying to capture the evolution of theG-factor as a function of size of the diffusing molecule (but clearly wrong when σ/Rpore≫1), which does not imply the approval of expression (4) and the corresponding definition of tortuosity and constrictivity.


## 3.2. Single Pore Size Hydraulic Conductivity Model

Due to the single pore size approximation in theG-factor, also a single pore size model is presented for advective flow. The flow Qcore (m3/s) out of a cylindrical clay core with radius Rcore (m) and length Lcore (m) is provided by Darcy's law,(7)Qcore=πRcore2ΔpLcoreKwith Δp being the pressure difference (in meter water column) over the core.Poiseuille's law gives the flowQpore (m3/s) through a cylindrical pore with radius Rpore (m)(8)Qpore=1μπRpore48ΔpLcorewith μ being the viscosity (kg/(m s)) (a value μ=0.7978 mPa s=8.135×10-8 m s was used (Edge Engineers)).IfCpore is the number of pores per unit surface (m−2), the flow Qcore also equals(9)Qcore=πRcore2CporeQporeand the core porosity ηtot is given by(10)ηtot=πRpore2Cpore.Combining expressions (7)–(10) leads to the next expression for the hydraulic conductivity(11)K=18μRpore2ηtot=π8μRpore4Cporeshowing that, for two clays with the same porosity, the hydraulic conductivity is the highest for the clay with the largest pore size. It is also clear that, in case of a pore size distribution, a major contribution to the hydraulic conductivity comes from the largest pores.


## 3.3. Size of the Diffusing Molecules

The size of a diffusing molecule can be characterized in several ways. As it would be not appropriate to use values which have been determined from diffusion measurements (hence leading to circular reasoning), some extra clarification on the selection of the kinetic diameters is given in this paragraph. The most widely used measure for the size of a (small) gas molecule is the kinetic diameter(σ) [54]. In this approach, the interaction potential between two (gas) molecules is assumed a Lennard-Jones VLJ(r) potential (with r distance between the centre of two molecules) [55].The potentialVLJ(r) is zero at r=σ and the distance σ (called the kinetic diameter) is considered representative for the size of the gas molecule. For the kinetic diameters used in this work, the parameters of the Lennard-Jones potential are derived from viscosity measurements at different temperatures.However, the selected kinetic diameters are representative for dilute gas mixtures and are thus only an approximation for gases dissolved in water. More information on this topic can be found in [23].


## 3.4. Size Dependency of the Diffusion Coefficient of Gases in Water and Clay

Empirically, the gas diffusion coefficientsD0 in water were fitted as a function of the kinetic diameter σ (expressed in 10-10m) by an exponential equation (12)D0=D0,expexp⁡-n0,expσleading to D0,exp=5.5±2.3×10-8 m2/s and n0,exp=(8.8±1.2)×109 m−1 [23]. Another possibility is a power law fit (13)D0=D0,powσ-n0,powleading to D0,pow=1.1±0.46×10-7 m2/s and n0,pow=3.0±0.4 [23]. Because expressions (12) and (13) do not take into account details like molecular geometry or polarity of a molecule, they only describe a global evolution without predicting correctly every single point (like, e.g., Archie's law).Similar to the diffusion coefficientD0 in free water, for gases (and HTO) trend lines for the dependence of the effective diffusion coefficient Deff on the kinetic diameter σ are fitted by assuming similar expressions as (12) and (13) [23]: (14)Deff=Deff,expexp⁡-nexpσ,(15)Deff=Deff,powσ-npow.Substituting (12) and (14) into (1) provides for unretarded tracers without ion exclusion (the retardation factor R=1) an expression for the geometric factor:(16)G=ηD0,expDeff,expexp⁡nexp-n0,expσ.This expression reduces to the Gratwohl expression (6) by putting(17)τ=ηD0,expDeff,exp(18)Rpore=4.6nexp-n0,exp,where a positive pore size Rpore requires that nexp>n0,exp.A similar combination as (16) can be made with expressions (13) and (15) leading to (and assuming npow>=n0,pow)(19)G=ηD0,powDeff,powσnpow-n0,pow.An expression similar to (17) can be defined for the power law (19), but no characteristic distance similar to Rpore exists. Also note that with the power law, G(σ=0)=0 (assuming npow>=n0,pow) instead of G(σ=0)=τ for the exponential relation (16).


## 3.5. Modelling

A detailed description on the modelling approach can be found in Jacops et al. [23].The diffusion experiments are modelled by fitting the solutions of the diffusion equation with the appropriate boundary and initial conditions. For both the gas diffusion and the HTO diffusion experiments, the diffusion equation is solved by COMSOL coupled with MATLAB for optimization.A through-diffusion experiment allows fitting both the apparent diffusion coefficientDapp and the capacity factor ηR, which are in general both unknown. Besides, also the initial concentration, which is known but prone to measurement errors, is fitted and this is described as the 3-parameter fit. Diffusion coefficients for HTO, reported in this paper have been obtained with this 3-parameter fit. An alternative fit strategy is also possible: assumingR = 1, only Dapp is fitted and ηR is a constant equal to ηtot (measured on either a sample or an average value from literature). In some gas through-diffusion experiments, the (quasi)stationary state is reached very fast, leading to a nearly total correlation between the capacity factor, the apparent diffusion coefficient, and the initial concentration. In this case, the one parameter fit is the only option and only the effective diffusion coefficient can be determined. Therefore, all diffusion coefficients for dissolved gases are obtained from a one-parameter fit. The one- and three-parameter fits can lead to small differences in the optimal values [23].The water content of a core was measured at the end of the diffusion experiments, and if this information was not available, a reference value obtained from literature was used (see also Table3).As discussed above, confining filters are used in the diffusion experiments. As discussed by Glaus et al. [56], Birgersson and Karnland [57], and Aertsens et al. [48] the filters might have an influence on the results as they contribute to the diffusive resistance of the experiment. Therefore, two fitting strategies are possible: one assuming the same transport parameter values for clay core and filter (case “Filter = Clay” – “F = C”) and one with fixed (=the measured but not very accurate values) filter parameter values (case “Fixed Filter,” “FF”). However, the only element for which the values of the transport parameters in the filters are measured is HTO [48]. For Boom Clay oriented perpendicular to the bedding plane, the values of the transport parameters are roughly similar to those of the filter so it can be a good approximation to assume the same transport parameters for the filter as for the clay core (case F = C). However, also Boom Clay cores oriented parallel to the bedding plane and other types of cores are investigated and, in these cases, the diffusion coefficient of the filter differs considerably from the HTO diffusion coefficient in the samples. Therefore, the most correct approach is to use in the model the known HTO filter diffusion coefficient (case fixed filter (FF)).Diffusion coefficients of dissolved gases in the filters are not available and therefore only the value of HTO could be used as an approximation. As this would introduce a large uncertainty into the model and given the large length of the samples (35 mm) compared to the length of the filters (2 × 2 mm), the effect of the filter on diffusion is not taken into account in the experiments with dissolved gases.


## 4. Results

### 4.1. Sample Characterization

Both the Boom Clay and the Eigenbilzen Sands are mainly composed of quartz and 2 : 1 clay minerals. As shown in Table2, the main differences between the samples of the Boom Clay and the Eigenbilzen Sands are related to the main components quartz and 2 : 1 clay minerals. Given the high quartz content and the lower content of 2 : 1 clay minerals, the samples of the Eigenbilzen Sands can be considered “sandy” and are categorized as sandy clay.Table 2
Quantitative mineralogical composition of the different samples (in mass%) and the specific surface area (m2/g). Minerals which were below the detection limit are indicated with BD.
Boom Clay
Eigenbilzen Sands
Bentonite
COX
OPA
Core 84b K2
Core 127b K4
Core 36a K14
Core 37b K15
Core 35b K16
Core 39b K17
Quartz
31
28
59
58
60
54
3
24
28
K-feldspar
8
5
8
9
10
8
BD
5
3
Plagioclase
3
1
6
6
5
5
3
3
BD
Calcite
0.2
2
BD
BD
BD
BD
BD
21
8
Ankerite/dolomite
BD
BD
BD
BD
BD
BD
BD
4
BD
Pyrite
2
2
0.6
0.5
1
0.5
BD
1
0.6
Gypsum
∗
0.6
0.5
1
1
0
0
BD
BD
BD
Anatase
0.6
0.7
BD
BD
BD
BD
BD
BD
BD
Kaolinite
8
9
2
4
2
3
BD
0
BD
2 : 1 Al clay
34
41
20
18
21
27
93
39
59
†
Muscovite
9
8
BD
BD
BD
BD
BD
BD
BD
Chlorite
2
2
3
3.5
0
2.5
BD
3
BD
Opal A
2
1
BD
BD
BD
BD
BD
BD
BD
Specific surface area
38
45
14
12
8
20
+
28
21
Gypsum
∗ is a secondary phase; it is an oxidation product of pyrite after exposure of the sample to air; this† value contains all clay minerals (counted together); not+ measured.Table 3
HTO Diffusion coefficients, fitted capacity factor, and total porosity assuming fixed filter transport parameter values (Fixed Filter, FF).
Material
Core
D
a
p
p
,
F
F
η
R
a
p
p
F
F
η
t
o
t
D
e
f
f
,
F
F
(m2/s)
(-)
(m2/s)
Boom Clay
84b (K2)
⊥
3.4
E
-
10
0.56
0.40
1.9
E
-
10
Boom Clay
127b (K4)
//
8.7
E
-
10
0.32
0.38
2.8
E
-
10
Eigenbilzen Sand
36a (K14)
//
1.2
E
-
09
0.27
0.39
3.2
E
-
10
Eigenbilzen Sand
37b (K15)
//
7.6
E
-
10
0.58
0.40
4.4
E
-
10
Eigenbilzen Sand
35b (K16)
⊥
7.6
E
-
10
0.44
0.41
3.3
E
-
10
Eigenbilzen Sand
39b (K17)
⊥
5.5
E
-
10
0.55
0.40
3.0
E
-
10
Bentonite 1.4
2.5
E
-
10
0.62
0.47
∗
1.6
E
-
10
Bentonite 1.6
1.4
E
-
10
0.58
0.4
∗
8.1
E
-
11
COX
⊥
NM
NM
0.16
NM
OPA
OPA 1
⊥
7.3
E
-
11
0.17
0.096
1.2
E
-
11
Porosity
∗ value from literature [60].The composition of the COX sample corresponds well to the reference values presented in [58]. The composition of the Opalinus Clay samples corresponds well to results obtained for samples located nearby our sample [47].


### 4.2. Hydraulic Conductivity

The hydraulic conductivity(K) values vary over several orders of magnitude (Table 4). The values of the Eigenbilzen Sands (the clayey sand samples ON-Mol1 core 35b, 36a, 37b, and 39b) display the highest values, which are in line with the reported average value of 3.1 × 10−10 m/s for the Eigenbilzen Sands from Aertsens et al. [59], except sample K16 which shows a considerable higher K value and for which we have no clear explanation. The hydraulic conductivity of the “clayey” Boom Clay samples 84b and 127b (discussed in [23]) is roughly two orders of magnitude lower. Still an order of magnitude lower is the hydraulic conductivity of the COX and OPA samples [22] and both bentonite samples. Evidently, the 1.6 density bentonite has a lower hydraulic conductivity than the 1.4 bentonite.Table 4
Overview of measured hydraulic conductivity and diffusion coefficientsDeff and the uncertainty on the fit (95% confidence interval).
Kinetic diameterσ(1E-10m) D0 (1E-9 m2/s)
He
HTO
Ne
H2
Ar
CH4
Xe
C2H6
2.58
2.75
2.79
2.97
3.42
3.82
4.06
4.42
7.28
2.20
4.03
5.11
2.44
1.84
1.47
1.38
Core
Code
K
η
t
o
t
D
e
f
f
D
eff
D
e
f
f
D
e
f
f
D
e
f
f
D
e
f
f
D
e
f
f
D
e
f
f
(m/s)
(-)
(
1
E
-
11
m
2
/
s
)
(
1
E
-
11
m
2
/
s
)
(
1
E
-
11
m
2
/
s
)
(
1
E
-
11
m
2
/
s
)
(
1
E
-
11
m
2
/
s
)
(
1
E
-
11
m
2
/
s
)
(
1
E
-
11
m
2
/
s
)
(
1
E
-
11
m
2
/
s
)
BC 84b
K2
⊥
3.3
E
-
12
0.4
46.8
±
1.7
18.7
±
0.6
17.5
±
0.3
NM
6.9
±
0.2
9.7
±
0.3
6.1
±
0.2
4.6
±
0.2
BC 127b
K4
//
3.4
E
-
12
0.38
74.7
±
2.0
27.8
±
0.9
22.9
±
1.0
51.2
±
1.1
14.5
±
0.2
15.5
±
0.5
6.6
±
0.9
5.9
±
0.1
BC 36a
K14
//
4.4
E
-
10
0.37
67.9
±
0.7
32.1
±
3.0
33.6
±
1.1
NM
NM
23.1
±
0.6
NM
14.3
±
0.1
BC 37b
K15
//
1.1
E
-
10
0.40
82.3
±
2.1
43.9
±
2.3
40.0
±
2.1
NM
NM
27.8
±
0.5
NM
15.2
±
0.2
BC 35b
K16
⊥
1.2
E
-
08
0.41
74.5
±
2.1
33.1
±
1.8
35.4
±
1.3
NM
NM
37.2
±
0.5
NM
13.7
±
0.4
BC 39b
K17
⊥
4.2
E
-
10
0.40
58.2
±
1.1
29.9
±
2.7
30.6
±
0.6
NM
NM
24.4
±
0.5
NM
8.9
±
0.3
Bent 1.4
Bent 1.4
1.5
E
-
13
0.47
17.7
±
1.1
15.6
±
1.2
10.3
±
0.4
NM
NM
2.8
±
0.1
NM
1.1
±
0.02
Bent 1.6
Bent 1.6
5.8
E
-
14
0.40
27.0
±
0.9
8.07
±
0.5
8.7
±
0.7
NM
NM
0.9
±
0.1
NM
0.3
±
0.04
COX 1
COX 1
⊥
3.4
E
-
13
0.16
∗
8.1
±
0.2
NM
2.1
±
0.1
NM
0.7
±
0.02
N
M
NM
0.2
±
0.01
OPA 1
OPA 1
⊥
1.2
E
-
13
0.096
6.8
±
0.3
1.2
±
0.04
0.6
±
0.0
NM
0.4
±
0.1
F
F
F
The
∗ measured porosity for the COX sample is 0.16. However, for compatibility reasons, the reference value of 0.18 is taken as this value has also been used in Jacops et al. [22]; NM indicates that a gas was “not measured,” while F indicates a “failed” experiment.Assuming that a soil consists of spherical grains, the pores largely consist of the space between the grains. In a very rough approximation, the pore size is related to the grain size and the sorting of the sedimented particles. Because the grain size of sand is larger than that of clay grains, according to expression (11) the hydraulic conductivity of the Boom Clay samples is expected to be smaller than the hydraulic conductivity of the samples of the Eigenbilzen Sands, which is confirmed by the measurements. Due to the lower porosity of COX and OPA, a lower hydraulic conductivity can be expected than for Boom Clay. Despite their relatively high porosity, the low hydraulic conductivity of the bentonites can be understood by their fine grain size and because of their composition (mainly montmorillonite). The interlayers in montmorillonite, with a thickness of a thin water layer, can be considered as “pores,” reducing the typical pore size.


### 4.3. Diffusion Experiments

#### 4.3.1. HTO Diffusion Results

As described above, the HTO experiments are fitted by using fixed filter (FF) values. The results are shown in Table3. Please note that the value for HTO is now obtained with fixed filter values, while the OPA value in Jacops et al. [22] was obtained considering filter = clay.For the capacity factors, there is for some samples a considerable difference between the fitted capacity factors and the total porosityηtot (for this parameter, sometimes a reference value is used, and sometimes it is measured by water content), in particular for OPA. However, the error on the fits is small (Table 4); hence the fits are considered to be good.


#### 4.3.2. Diffusion Coefficients and Geometric Factors of Dissolved Gases

All measured effective diffusion coefficients and the corresponding geometrical factors for dissolved gases and HTO are reported in Tables4 and 5. The COX and OPA (He and Ar) results have already been reported in [22], as well as the Boom Clay samples ON-Mol1 84b and 127b in [23]. Please note that the fits for COX were performed with a porosity of 18%, in order to be consistent with Jacops et al. [22].Table 5
Geometrical factors obtained from the aqueous and effective diffusion coefficients from Table4. Note that similar to the previous table, the reference value of 0.18 is used for the COX1 total porosity.
Core
Code
η
t
o
t (-)
He
HTO
Ne
H2
Ar
CH4
Xe
C2H6
G (-)
BC 84b
K2
⊥
0.4
6.2
4.7
9.2
14.1
7.6
9.6
11.9
BC 127b
K4
//
0.38
3.7
3.0
6.7
3.8
6.4
4.5
8.4
8.9
BC 36a
K14
//
0.37
4.0
2.5
4.4
2.9
3.6
BC 37b
K15
//
0.4
3.5
2.0
4.0
2.6
3.6
BC 35b
K16
⊥
0.41
4.0
2.7
4.7
2.0
4.1
BC 39b
K17
⊥
0.4
5.0
2.9
5.3
3.0
6.2
Bent 1.4
Bent 1.4
0.47
19.3
6.6
18.4
30.7
60.3
Bent 1.6
Bent 1.6
0.4
10.8
10.9
18.5
79.9
161.7
COX 1
COX 1
⊥
0.18
16.1
34.6
60.2
131.1
OPA 1
OPA 1
⊥
0.096
10.2
17.5
62.1
54.0Opalinus Clay has the lowest porosity of all the samples used in this paper. It is evidently the first candidate showing possible percolation transition problems. These occurred for all gases larger than Argon: CH4, Xe, and C2H6:(i)
For CH4 the measurements were scattered and all below 100 ppm. As 100 ppm is the lower limit for reliable CH4 peak detection and measurement, all measurements contain a significant portion of noise and are therefore considered nonsignificant.(ii)
For Xe, no breakthrough was measured which could also be related to the detection limit of 100 ppm.(iii)
For C2H6, where the detection limit is only 5 ppm, a clear breakthrough curve was detected, but it could not be fitted by a simple diffusion model (see Figure 1): a quasi-stationary state is reached very soon, but after about 150–200 days, the outlet concentration starts to increase much faster. Because the system is probably very close to the percolation threshold, a possible explanation might be diffusion over another percolating, longer pathway.Figure 1
Evolution of the outlet concentration versus time in the C2H6 diffusion experiment in Opalinus Clay. A clear breakthrough after a couple of days, leading to a quasi-stationary state is followed by an sudden unexplained increase after more than 150 days.When comparing the measured diffusion coefficients in the Boom Clay samples and the samples from the Eigenbilzen Sands, one can observe that despite a comparable porosity, the diffusion coefficients in the Eigenbilzen Sands are slightly higher than in Boom Clay. This agrees with previous results from pulse injection experiments [59], showing on average higher HTO dispersion coefficients in the Eigenbilzen Sands than in Boom Clay. Evidently, the present results confirm that diffusion coefficients are anisotropic and higher for cores parallel to the bedding compared to cores perpendicular to the bedding. Despite its higher porosity compared to the Boom Clay and the Eigenbilzen Sands, the diffusion coefficients in bentonite 1.4 are smaller (as was also the case for the hydraulic conductivity). For helium in bentonite, we observed that diffusion coefficient in bentonite 1.4 is considerably lower than for bentonite 1.6. In fact, one would expect the opposite and we cannot explain this observation. The other three gases show the expected behaviour: lower diffusion coefficients in bentonite 1.6 compared to bentonite 1.4. The lowest diffusion coefficients were observed for the clays with the lowest porosity: COX and OPA.Also, note that, for the samples from the Eigenbilzen Sands and for the Boom Clay samples, the diffusion coefficients for HTO and Ne (two molecules with a nearly equal kinetic diameter) are about the same. Due the very different values for the HTO and Ne aqueous diffusion coefficients (D0), the corresponding geometric factors differ considerably.


#### 4.3.3. The Effective Diffusion Coefficients and the Geometric Factor as a Function of the Size of the Diffusing Molecule

Table6 shows the fitted values according to expression (14) from the measured effective gas diffusion coefficients and the corresponding graphs are shown in Figure 2. Although HTO is included in Table 4, it is not included in the fits (shown in Table 6) as we only describe the data for dissolved gases.Table 6
Parameter values obtained by fitting the effective diffusion coefficient versus molecular size by the exponential expression (14). The numbers in the top line are taken from Jacops et al. [23]. Rpore is calculated from expression (18) and Rpore,K is calculated from expression (11). From the fitted Rpore value, the hydraulic conductivity is calculated (Kcalc) and compared to the experimental value Kexp.
Water
D
0
,
e
x
p
n
0
,
e
x
p
(m2/s)
(m−1)
5.5
E
-
08
±
2.3
E
-
08
8.8
E
+
09
±
1.2
E
+
09
Code
η
t
o
t
D
e
f
f
,
e
x
p
n
e
x
p
τ
R
p
o
r
e
R
p
o
r
e
,
K
K
c
a
l
c
K
e
x
p
K
e
x
p
/
K
c
a
l
c
(-)
(m2/s)
(m−1)
( )
(m)
(m)
(m/s)
(m/s)
K2
⊥
0.4
2.0
E
-
09
±
2.1
E
-
09
8.7E + 09 ± 3.0E + 09
11.1
−1.53E − 07
2.32
E
-
09
1.4
E
-
08
3.3
E
-
12
0.0002
K4
//
0.38
5.7
E
-
09
±
5.4
E
-
09
1.1
E
+
10
±
2.7
E
+
09
3.7
2.51
E
-
09
2.41
E
-
09
3.7
E
-
12
3.4
E
-
12
0.92
K14
//
0.37
2.4
E
-
09
±
1.6
E
-
09
6.4E + 09 ± 1.9E + 09
8.4
−1.91E − 09
2.78
E
-
08
2.1
E
-
12
4.4
E
-
10
212
K15
//
0.4
3.4
E
-
09
±
2.4
E
-
09
6.9E + 09 ± 2.0E + 09
6.4
−2.51E − 09
1.34
E
-
08
3.9
E
-
12
1.1
E
-
10
28
K16
0.41
2.8
E
-
09
±
2.7
E
-
09
6.6E + 09 ± 2.7E + 09
8.0
−2.08E − 09
1.38
E
-
07
2.7
E
-
12
1.2
E
-
08
4397
K17
⊥
0.4
3.7
E
-
09
±
2.9
E
-
09
8.2E + 09 ± 2.2E + 09
5.9
−8.21E − 09
2.61
E
-
08
4.1
E
-
11
4.2
E
-
10
10.1
Bent 1.4
0.47
6.4
E
-
09
±
2.0
E
-
09
1.4
E
+
10
±
8.7
E
+
08
4.0
8.17
E
-
10
4.56
E
-
10
4.8
E
-
13
1.5
E
-
13
0.31
Bent 1.6
0.4
4.5
E
-
08
±
4.4
E
-
08
2.2
E
+
10
±
2.8
E
+
09
0.5
3.53
E
-
10
3.07
E
-
10
7.7
E
-
14
5.8
E
-
14
0.76
COX 1
⊥
0.18
2.3
E
-
09
±
3.0
E
-
09
1.6
E
+
10
±
3.8
E
+
09
4.3
6.03
E
-
10
1.11
E
-
09
1.0
E
-
13
3.4
E
-
13
3.4
OPA 1
⊥
0.096
9.6
E
-
11
±
5.7
E
-
10
9.3
E
+
09±1.9E+10
55.0
9.20
E
-
09
9.02
E
-
10
1.2
E
-
11
1.2
E
-
13
0.010Figure 2
Effective diffusion coefficients for dissolved gases in Boom Clay (a) and other clayey materials (b), fitted with expression (14). The diffusion coefficientD is either the diffusion coefficientD0 in pure water or the effective diffusion coefficient Deff in a porous medium.
(a)
(b)Most of the fits of Table6 are fair but show a relatively large error. The OPA fit is not reliable: for both fit parameters, the error is larger than 100%.For all Boom Clay and Eigenbilzen Sands samples, the exponential factornexp is on average (7.9 ± 2) × 10−10 which is similar to the exponential factor n0,exp of the fit of D08.8±1.2×10-9. For the other clayey samples, the factor nexp is always substantially larger than n0,exp.As discussed earlier, for gases, the geometric factorG can be calculated from diffusion experiments. Similar to the evolution of Deff with size; also the evolution of G with size can be described with an exponential or a power law function and this is shown in Figure 3.Figure 3
Geometric factors for dissolved gases in Boom Clay (a) and other clayey materials (b).
(a)
(b)From Figure3 and Table 5, one can observe that for the Boom Clay sample 127b (K4) the geometric factor tends to increase with the kinetic diameter, but the difference in G between the smallest (He) and largest (C2H6) molecule is smaller than a factor 3.For the Eigenbilzen Sands samples, the geometric factor is rather a constant value. This can also be deduced from the exponential factors of the fit ofDeff (Table 6): if n0,exp and nexp are similar (nexp-n0,exp≈0), according to expression (16), the variation of G with molecular size is small, in agreement with the experimental results (see Table 5).For the other clayey samples, the geometric factor clearly increases quite spectacularly with the kinetic diameter: theG ratio between smallest and largest molecule ranges from a factor 3 (bentonite 1.4) up to a factor 16 (Bentonite 1.6). This is also reflected in the exponential factors nexp which differ much more from the n0,exp factor compared to the Boom Clay and Eigenbilzen Sands samples (see Table 6).From the fits in Table6, the tortuosity τ (expression (17)) and Rpore (expression (18)) can be calculated. The calculated tortuosity factors τ are, apart from the unreliable OPA value and the small Bent 1.6 value, quite similar. But the calculated Rpore values are sometimes negative which is physically impossible (Boom Clay sample K2 and the four samples of the Eigenbilzen Sands). For these samples, the values of n0,exp and nexp are within each other’s error interval which means that in some cases n0,exp>nexp, which leads according to expression (18) to negative Rpore values. For these samples, the difference in Deff between the largest and smallest molecule is about one order of magnitude. So a slight change in one of the Deff values can strongly influence nexp. This observation points to a limitation of this approach: it can only be used when nexp differs significantly from n0,exp. From the measured K values, Rpore,K can be calculated according to expression (11) and these results are also shown in Table 6. When comparing both calculated Rpore values, it is clear that when a reliable exponential fit is obtained, and when nexp differs significantly from n0,exp (so no negative Rpore values) both calculated Rpore values are relatively similar.According to expression (6) which relates the geometric factor to the size of the diffusing molecule and the pore size, the geometrical factor always increases as a function of the size of the diffusing molecule. For a sufficiently small molecular size σ (with respect to the pore size Rpore) the exponential factor exp⁡(σ/Rpore) can be approximated by one, so that the geometrical factor is approximately constant. This means that, according to expression (6), the geometric factor for a very small diffusing molecule (kinetic diameter σ→0) is roughly similar in all these clays and the pore size Rpore is not relevant. This applies most to the smallest molecule, helium, for which indeed the variability of G over the different samples is small. From expression (6), one can also conclude that if, for the same sample, all diffusing molecules are sufficiently small compared to the pore size, the geometric factor is also approximately constant. This is the case for the samples of the Eigenbilzen Sands: as Rpore,K (calculated from expression (11)) is around 2 × 10−8m, the pores are much larger than the diffusing molecules.When the diffusing molecules are not small compared to the pore size, the pore size becomes more important. When, in this case, the size of the diffusing molecule increases, the geometric factors in the different clays diverge, and the geometric factor for different gases in a single sample is not constant. For both bentonites and COX, both fittedRpore values are small (see Table 6) (similar to the size of the diffusing molecules) and the geometric factor clearly rises with increasing molecular size (see Figure 3). For Boom Clay, sample K4 which is expected to have larger pores compared to COX and the bentonites, the geometric factor also increases with the size of the diffusing molecule, although less steeply.Despite the enormously simplifying assumption to replace a whole pore distribution by just one pore size in both the hydraulic conductivity model and the effective diffusion coefficient versus diffusing molecule size model, both models lead to similar values for the typical pore size. Using the fittedRpore values to predict the hydraulic conductivity with expression (11) gives quite a good agreement between the calculated hydraulic conductivity Kcalc and the experimental value Kexp for both bentonites, COX and Boom Clay sample K4 (see Table 6). It is straightforward that, for negative Rpore values (sample K2 and the Eigenbilzen Sands), this prediction cannot be good, and it has already been mentioned that the OPA fits are not reliable.The effective diffusion coefficients can also be fitted as a function of size by using the power law expression (19), leading to the optimal values summarized in Table 7. Apart from the fact that for the power law no typical pore size like Rpore exists, the conclusions are similar as for the exponential fit: a bad fit for OPA, a steep increase of the geometric factor as a function of molecular size for both bentonites and COX, a slower increase for both Boom Clay samples, and a hard-to-explain (apart from measurement errors) decrease (npow-n0,pow<0) for the Eigenbilzen Sands (indicated in italic).Table 7
Parameter values obtained by fitting the effective diffusion coefficient versus size of the diffusing molecule by the power law expression (15). The numbers in the top line are taken from [23]. “Tortuosity” τ is defined similarly as expression (17).
Water
D
0
,
p
o
w
n
0
,
p
o
w
(m2/s)
( )
1.1
E
-
07
±
4.6
E
-
08
3.0
±
0.4
Code
η
t
o
t
D
e
f
f
,
p
o
w
n
p
o
w
τ
n
p
o
w
-
n
0
,
p
o
w
(-)
(m2/s)
(m−1)
( )
( )
K2
⊥
0.4
5.0
E
-
09
±
5.8
E
-
09
3.2
±
0.9
8.4
0.2
K4
//
0.38
1.3
E
-
08
±
1.4
E
-
08
3.7
±
0.9
3.1
0.7
K14
//
0.37
4.0
E
-
09
±
3.1
E
-
09
2.2 ± 0.6
10.0
−0.8
K15
//
0.4
5.9
E
-
09
±
4.8
E
-
09
2.4 ± 0.7
6.6
−0.6
K16
⊥
0.41
4.8
E
-
09
±
5.5
E
-
09
2.3 ± 0.9
8.3
−0.7
K17
⊥
0.4
7.0
E
-
09
±
6.8
E
-
09
2.9 ± 0.8
5.5
−0.2
Bent 1.4
0.47
1.8
E
-
08
±
7.2
E
-
09
5.0
±
0.3
2.8
1.9
Bent 1.6
0.4
2.3
E
-
07
±
1.9
E
-
07
7.5
±
0.7
0.2
4.5
COX 1
⊥
0.18
1.0
E
-
08
±
1.3
E
-
08
5.8
±
1.1
1.9
2.8
OPA 1
⊥
0.096
1.1
E
-
10
±
7.1
E
-
10
2.7
±
5.9
94.2
−0.4


### 4.4. How to Estimate Gas Diffusion Coefficients as a Function of the Size of the Diffusing Molecule?

Based on the results presented in this work, one can state that the approach of using one geometric factor for estimatingDeff of species with a different molecular size is not always correct and the validity depends on the type of sample and the size of the diffusing molecule. For instance, for the Eigenbilzen Sand samples, the geometric factor from HTO could be used to calculate Deff for, for example, C2H6. On the contrary, for COX or bentonite, using the geometric factor of HTO to calculate Deff for C2H6 would lead to a serious overestimation. In order to illustrate this, we use for COX a geometric factor of 17 for HTO, as suggested by Jacops et al. [22]. If this value would be used to predict Deff for C2H6, a value of 1.5 × 10−11 m2/s would be obtained while the measured value was 0.2 × 10−11 m2/s so this leads to an overestimation of factor 7.A more correct way to estimate diffusion coefficients would be by using the exponential or power law relation betweenDeff and the molecular size (expressions (14) and (15)). Each expression has two parameters, meaning that for two gases the effective diffusion coefficient needs to be determined experimentally. Evidently, it is best to do this for a large molecule (e.g., C2H6) and a small one. He is the smallest of all molecules considered here, but its diffusion coefficients are always higher than expected (located above the fitted curve in Figure 2). This leads to the observation that the measured Deff is on average 1.7 times larger than the calculated Deff (see Table 8). For the second smallest molecule (neon), the calculated average Deff is similar to the measured Deff; hence Ne is recommended for being used in the fits. An example of such a prediction with the exponential expression (14) is shown in Table 9. Evidently, because Ne and C2H6 are used to determine the fit parameters, for both of these gases, the predicted value corresponds exactly to the experimental value. For the other gases, the prediction can be improved by multiplying the predicted diffusion coefficient with the averages of the ratio of the experimental to fitted value for all measured clays (see Table 8). This manipulation allows correcting for the specific bias of this ratio, which is different for each gas. After this correction, for sample K4 (one of the experiments for which most measurements were carried out), apart from He, the predicted and measured effective diffusion coefficients differ up to less than 30% (see Table 9).Table 8
Ratios between the experimentally determinedDeff (values see Table 4) and the fitted Deff (using the fits given in Table 6), as well as the averages for each element. Due to the deviating value of He for sample OPA1, the corresponding Deff is not included in the average.
Core
Code
He
Ne
H2
Ar
CH4
Xe
C2H6
Ratio Exp/Fit (-)
BC 84b
K2
2.3
1.0
0.7
1.4
1.1
1.1
BC 127b
K4
2.0
0.8
2.1
1.0
1.6
0.9
1.1
BC 36a
K14
1.4
0.8
1.1
1.0
BC 37b
K15
1.4
0.8
1.2
1.0
BC 35b
K16
1.4
0.8
1.6
0.9
BC 39b
K17
1.3
0.8
1.5
0.9
Bent 1.4
Bent 1.4
1.1
0.9
1.1
1.0
Bent 1.6
Bent 1.6
1.7
0.8
0.8
1.2
COX 1
COX 1
2.4
0.9
0.9
1.2
OPA 1
OPA 1
7.8
0.9
1.1
Average
1.7
0.8
2.1
0.9
1.3
1.0
1.0
Standard deviation
0.4
0.1
0.2
0.3
0.2
0.1
For
∗ He, OPA1 not included.Table 9
Predicting the gas effective diffusion coefficients of sample K4 after measuring the Ne and C2H6 (both in italic) effective diffusion coefficients. In a first phase, the prediction is based only on expression (14). Subsequently, this prediction is improved by multiplication with a “correction factor” (see Table 8), which takes into account the gas specific deviation between the experimental and the predicted value.
He
Ne
H2
Ar
CH4
Xe
C
2
H
6
Kinetic diameterσ
(1E-10 m)
2.58
2.79
2.97
3.42
3.82
4.06
4.42
D
e
f
f (experiment)
(1E-11 m2/s)
74.7
22.9
51.2
14.5
15.5
6.6
5.9
D
e
f
f (prediction)
(
1
E
-
11 m2/s)
27.2
22.9
19.7
13.6
9.7
8.0
5.9
RatioDeff (exp)/Deff (pred)
( )
2.7
1.0
2.6
1.1
1.6
0.8
1.0
Correction factor
( )
1.7
0.8
2.1
0.9
1.3
1.0
1.0
D
e
f
f (pred) ∗ Corr factor
(1E-11 m2/s)
45.8
41.2
12.2
12.4
7.7
RatioDeff (exp)/(Deff (pred) ∗ Corr factor)
( )
1.6
1.2
1.2
1.2
0.9By using this approach, the diffusion coefficient of, for instance, H2 (which is difficult to measure due to microbial activity [24]), can be estimated based on its molecular size. But note that the diffusion coefficient for H2 was (due to experimental difficulties) measured for only one sample; hence the correction factor might be only an estimate.


## 4.1. Sample Characterization

Both the Boom Clay and the Eigenbilzen Sands are mainly composed of quartz and 2 : 1 clay minerals. As shown in Table2, the main differences between the samples of the Boom Clay and the Eigenbilzen Sands are related to the main components quartz and 2 : 1 clay minerals. Given the high quartz content and the lower content of 2 : 1 clay minerals, the samples of the Eigenbilzen Sands can be considered “sandy” and are categorized as sandy clay.Table 2
Quantitative mineralogical composition of the different samples (in mass%) and the specific surface area (m2/g). Minerals which were below the detection limit are indicated with BD.
Boom Clay
Eigenbilzen Sands
Bentonite
COX
OPA
Core 84b K2
Core 127b K4
Core 36a K14
Core 37b K15
Core 35b K16
Core 39b K17
Quartz
31
28
59
58
60
54
3
24
28
K-feldspar
8
5
8
9
10
8
BD
5
3
Plagioclase
3
1
6
6
5
5
3
3
BD
Calcite
0.2
2
BD
BD
BD
BD
BD
21
8
Ankerite/dolomite
BD
BD
BD
BD
BD
BD
BD
4
BD
Pyrite
2
2
0.6
0.5
1
0.5
BD
1
0.6
Gypsum
∗
0.6
0.5
1
1
0
0
BD
BD
BD
Anatase
0.6
0.7
BD
BD
BD
BD
BD
BD
BD
Kaolinite
8
9
2
4
2
3
BD
0
BD
2 : 1 Al clay
34
41
20
18
21
27
93
39
59
†
Muscovite
9
8
BD
BD
BD
BD
BD
BD
BD
Chlorite
2
2
3
3.5
0
2.5
BD
3
BD
Opal A
2
1
BD
BD
BD
BD
BD
BD
BD
Specific surface area
38
45
14
12
8
20
+
28
21
Gypsum
∗ is a secondary phase; it is an oxidation product of pyrite after exposure of the sample to air; this† value contains all clay minerals (counted together); not+ measured.Table 3
HTO Diffusion coefficients, fitted capacity factor, and total porosity assuming fixed filter transport parameter values (Fixed Filter, FF).
Material
Core
D
a
p
p
,
F
F
η
R
a
p
p
F
F
η
t
o
t
D
e
f
f
,
F
F
(m2/s)
(-)
(m2/s)
Boom Clay
84b (K2)
⊥
3.4
E
-
10
0.56
0.40
1.9
E
-
10
Boom Clay
127b (K4)
//
8.7
E
-
10
0.32
0.38
2.8
E
-
10
Eigenbilzen Sand
36a (K14)
//
1.2
E
-
09
0.27
0.39
3.2
E
-
10
Eigenbilzen Sand
37b (K15)
//
7.6
E
-
10
0.58
0.40
4.4
E
-
10
Eigenbilzen Sand
35b (K16)
⊥
7.6
E
-
10
0.44
0.41
3.3
E
-
10
Eigenbilzen Sand
39b (K17)
⊥
5.5
E
-
10
0.55
0.40
3.0
E
-
10
Bentonite 1.4
2.5
E
-
10
0.62
0.47
∗
1.6
E
-
10
Bentonite 1.6
1.4
E
-
10
0.58
0.4
∗
8.1
E
-
11
COX
⊥
NM
NM
0.16
NM
OPA
OPA 1
⊥
7.3
E
-
11
0.17
0.096
1.2
E
-
11
Porosity
∗ value from literature [60].The composition of the COX sample corresponds well to the reference values presented in [58]. The composition of the Opalinus Clay samples corresponds well to results obtained for samples located nearby our sample [47].


## 4.2. Hydraulic Conductivity

The hydraulic conductivity(K) values vary over several orders of magnitude (Table 4). The values of the Eigenbilzen Sands (the clayey sand samples ON-Mol1 core 35b, 36a, 37b, and 39b) display the highest values, which are in line with the reported average value of 3.1 × 10−10 m/s for the Eigenbilzen Sands from Aertsens et al. [59], except sample K16 which shows a considerable higher K value and for which we have no clear explanation. The hydraulic conductivity of the “clayey” Boom Clay samples 84b and 127b (discussed in [23]) is roughly two orders of magnitude lower. Still an order of magnitude lower is the hydraulic conductivity of the COX and OPA samples [22] and both bentonite samples. Evidently, the 1.6 density bentonite has a lower hydraulic conductivity than the 1.4 bentonite.Table 4
Overview of measured hydraulic conductivity and diffusion coefficientsDeff and the uncertainty on the fit (95% confidence interval).
Kinetic diameterσ(1E-10m) D0 (1E-9 m2/s)
He
HTO
Ne
H2
Ar
CH4
Xe
C2H6
2.58
2.75
2.79
2.97
3.42
3.82
4.06
4.42
7.28
2.20
4.03
5.11
2.44
1.84
1.47
1.38
Core
Code
K
η
t
o
t
D
e
f
f
D
eff
D
e
f
f
D
e
f
f
D
e
f
f
D
e
f
f
D
e
f
f
D
e
f
f
(m/s)
(-)
(
1
E
-
11
m
2
/
s
)
(
1
E
-
11
m
2
/
s
)
(
1
E
-
11
m
2
/
s
)
(
1
E
-
11
m
2
/
s
)
(
1
E
-
11
m
2
/
s
)
(
1
E
-
11
m
2
/
s
)
(
1
E
-
11
m
2
/
s
)
(
1
E
-
11
m
2
/
s
)
BC 84b
K2
⊥
3.3
E
-
12
0.4
46.8
±
1.7
18.7
±
0.6
17.5
±
0.3
NM
6.9
±
0.2
9.7
±
0.3
6.1
±
0.2
4.6
±
0.2
BC 127b
K4
//
3.4
E
-
12
0.38
74.7
±
2.0
27.8
±
0.9
22.9
±
1.0
51.2
±
1.1
14.5
±
0.2
15.5
±
0.5
6.6
±
0.9
5.9
±
0.1
BC 36a
K14
//
4.4
E
-
10
0.37
67.9
±
0.7
32.1
±
3.0
33.6
±
1.1
NM
NM
23.1
±
0.6
NM
14.3
±
0.1
BC 37b
K15
//
1.1
E
-
10
0.40
82.3
±
2.1
43.9
±
2.3
40.0
±
2.1
NM
NM
27.8
±
0.5
NM
15.2
±
0.2
BC 35b
K16
⊥
1.2
E
-
08
0.41
74.5
±
2.1
33.1
±
1.8
35.4
±
1.3
NM
NM
37.2
±
0.5
NM
13.7
±
0.4
BC 39b
K17
⊥
4.2
E
-
10
0.40
58.2
±
1.1
29.9
±
2.7
30.6
±
0.6
NM
NM
24.4
±
0.5
NM
8.9
±
0.3
Bent 1.4
Bent 1.4
1.5
E
-
13
0.47
17.7
±
1.1
15.6
±
1.2
10.3
±
0.4
NM
NM
2.8
±
0.1
NM
1.1
±
0.02
Bent 1.6
Bent 1.6
5.8
E
-
14
0.40
27.0
±
0.9
8.07
±
0.5
8.7
±
0.7
NM
NM
0.9
±
0.1
NM
0.3
±
0.04
COX 1
COX 1
⊥
3.4
E
-
13
0.16
∗
8.1
±
0.2
NM
2.1
±
0.1
NM
0.7
±
0.02
N
M
NM
0.2
±
0.01
OPA 1
OPA 1
⊥
1.2
E
-
13
0.096
6.8
±
0.3
1.2
±
0.04
0.6
±
0.0
NM
0.4
±
0.1
F
F
F
The
∗ measured porosity for the COX sample is 0.16. However, for compatibility reasons, the reference value of 0.18 is taken as this value has also been used in Jacops et al. [22]; NM indicates that a gas was “not measured,” while F indicates a “failed” experiment.Assuming that a soil consists of spherical grains, the pores largely consist of the space between the grains. In a very rough approximation, the pore size is related to the grain size and the sorting of the sedimented particles. Because the grain size of sand is larger than that of clay grains, according to expression (11) the hydraulic conductivity of the Boom Clay samples is expected to be smaller than the hydraulic conductivity of the samples of the Eigenbilzen Sands, which is confirmed by the measurements. Due to the lower porosity of COX and OPA, a lower hydraulic conductivity can be expected than for Boom Clay. Despite their relatively high porosity, the low hydraulic conductivity of the bentonites can be understood by their fine grain size and because of their composition (mainly montmorillonite). The interlayers in montmorillonite, with a thickness of a thin water layer, can be considered as “pores,” reducing the typical pore size.


## 4.3. Diffusion Experiments

### 4.3.1. HTO Diffusion Results

As described above, the HTO experiments are fitted by using fixed filter (FF) values. The results are shown in Table3. Please note that the value for HTO is now obtained with fixed filter values, while the OPA value in Jacops et al. [22] was obtained considering filter = clay.For the capacity factors, there is for some samples a considerable difference between the fitted capacity factors and the total porosityηtot (for this parameter, sometimes a reference value is used, and sometimes it is measured by water content), in particular for OPA. However, the error on the fits is small (Table 4); hence the fits are considered to be good.


### 4.3.2. Diffusion Coefficients and Geometric Factors of Dissolved Gases

All measured effective diffusion coefficients and the corresponding geometrical factors for dissolved gases and HTO are reported in Tables4 and 5. The COX and OPA (He and Ar) results have already been reported in [22], as well as the Boom Clay samples ON-Mol1 84b and 127b in [23]. Please note that the fits for COX were performed with a porosity of 18%, in order to be consistent with Jacops et al. [22].Table 5
Geometrical factors obtained from the aqueous and effective diffusion coefficients from Table4. Note that similar to the previous table, the reference value of 0.18 is used for the COX1 total porosity.
Core
Code
η
t
o
t (-)
He
HTO
Ne
H2
Ar
CH4
Xe
C2H6
G (-)
BC 84b
K2
⊥
0.4
6.2
4.7
9.2
14.1
7.6
9.6
11.9
BC 127b
K4
//
0.38
3.7
3.0
6.7
3.8
6.4
4.5
8.4
8.9
BC 36a
K14
//
0.37
4.0
2.5
4.4
2.9
3.6
BC 37b
K15
//
0.4
3.5
2.0
4.0
2.6
3.6
BC 35b
K16
⊥
0.41
4.0
2.7
4.7
2.0
4.1
BC 39b
K17
⊥
0.4
5.0
2.9
5.3
3.0
6.2
Bent 1.4
Bent 1.4
0.47
19.3
6.6
18.4
30.7
60.3
Bent 1.6
Bent 1.6
0.4
10.8
10.9
18.5
79.9
161.7
COX 1
COX 1
⊥
0.18
16.1
34.6
60.2
131.1
OPA 1
OPA 1
⊥
0.096
10.2
17.5
62.1
54.0Opalinus Clay has the lowest porosity of all the samples used in this paper. It is evidently the first candidate showing possible percolation transition problems. These occurred for all gases larger than Argon: CH4, Xe, and C2H6:(i)
For CH4 the measurements were scattered and all below 100 ppm. As 100 ppm is the lower limit for reliable CH4 peak detection and measurement, all measurements contain a significant portion of noise and are therefore considered nonsignificant.(ii)
For Xe, no breakthrough was measured which could also be related to the detection limit of 100 ppm.(iii)
For C2H6, where the detection limit is only 5 ppm, a clear breakthrough curve was detected, but it could not be fitted by a simple diffusion model (see Figure 1): a quasi-stationary state is reached very soon, but after about 150–200 days, the outlet concentration starts to increase much faster. Because the system is probably very close to the percolation threshold, a possible explanation might be diffusion over another percolating, longer pathway.Figure 1
Evolution of the outlet concentration versus time in the C2H6 diffusion experiment in Opalinus Clay. A clear breakthrough after a couple of days, leading to a quasi-stationary state is followed by an sudden unexplained increase after more than 150 days.When comparing the measured diffusion coefficients in the Boom Clay samples and the samples from the Eigenbilzen Sands, one can observe that despite a comparable porosity, the diffusion coefficients in the Eigenbilzen Sands are slightly higher than in Boom Clay. This agrees with previous results from pulse injection experiments [59], showing on average higher HTO dispersion coefficients in the Eigenbilzen Sands than in Boom Clay. Evidently, the present results confirm that diffusion coefficients are anisotropic and higher for cores parallel to the bedding compared to cores perpendicular to the bedding. Despite its higher porosity compared to the Boom Clay and the Eigenbilzen Sands, the diffusion coefficients in bentonite 1.4 are smaller (as was also the case for the hydraulic conductivity). For helium in bentonite, we observed that diffusion coefficient in bentonite 1.4 is considerably lower than for bentonite 1.6. In fact, one would expect the opposite and we cannot explain this observation. The other three gases show the expected behaviour: lower diffusion coefficients in bentonite 1.6 compared to bentonite 1.4. The lowest diffusion coefficients were observed for the clays with the lowest porosity: COX and OPA.Also, note that, for the samples from the Eigenbilzen Sands and for the Boom Clay samples, the diffusion coefficients for HTO and Ne (two molecules with a nearly equal kinetic diameter) are about the same. Due the very different values for the HTO and Ne aqueous diffusion coefficients (D0), the corresponding geometric factors differ considerably.


### 4.3.3. The Effective Diffusion Coefficients and the Geometric Factor as a Function of the Size of the Diffusing Molecule

Table6 shows the fitted values according to expression (14) from the measured effective gas diffusion coefficients and the corresponding graphs are shown in Figure 2. Although HTO is included in Table 4, it is not included in the fits (shown in Table 6) as we only describe the data for dissolved gases.Table 6
Parameter values obtained by fitting the effective diffusion coefficient versus molecular size by the exponential expression (14). The numbers in the top line are taken from Jacops et al. [23]. Rpore is calculated from expression (18) and Rpore,K is calculated from expression (11). From the fitted Rpore value, the hydraulic conductivity is calculated (Kcalc) and compared to the experimental value Kexp.
Water
D
0
,
e
x
p
n
0
,
e
x
p
(m2/s)
(m−1)
5.5
E
-
08
±
2.3
E
-
08
8.8
E
+
09
±
1.2
E
+
09
Code
η
t
o
t
D
e
f
f
,
e
x
p
n
e
x
p
τ
R
p
o
r
e
R
p
o
r
e
,
K
K
c
a
l
c
K
e
x
p
K
e
x
p
/
K
c
a
l
c
(-)
(m2/s)
(m−1)
( )
(m)
(m)
(m/s)
(m/s)
K2
⊥
0.4
2.0
E
-
09
±
2.1
E
-
09
8.7E + 09 ± 3.0E + 09
11.1
−1.53E − 07
2.32
E
-
09
1.4
E
-
08
3.3
E
-
12
0.0002
K4
//
0.38
5.7
E
-
09
±
5.4
E
-
09
1.1
E
+
10
±
2.7
E
+
09
3.7
2.51
E
-
09
2.41
E
-
09
3.7
E
-
12
3.4
E
-
12
0.92
K14
//
0.37
2.4
E
-
09
±
1.6
E
-
09
6.4E + 09 ± 1.9E + 09
8.4
−1.91E − 09
2.78
E
-
08
2.1
E
-
12
4.4
E
-
10
212
K15
//
0.4
3.4
E
-
09
±
2.4
E
-
09
6.9E + 09 ± 2.0E + 09
6.4
−2.51E − 09
1.34
E
-
08
3.9
E
-
12
1.1
E
-
10
28
K16
0.41
2.8
E
-
09
±
2.7
E
-
09
6.6E + 09 ± 2.7E + 09
8.0
−2.08E − 09
1.38
E
-
07
2.7
E
-
12
1.2
E
-
08
4397
K17
⊥
0.4
3.7
E
-
09
±
2.9
E
-
09
8.2E + 09 ± 2.2E + 09
5.9
−8.21E − 09
2.61
E
-
08
4.1
E
-
11
4.2
E
-
10
10.1
Bent 1.4
0.47
6.4
E
-
09
±
2.0
E
-
09
1.4
E
+
10
±
8.7
E
+
08
4.0
8.17
E
-
10
4.56
E
-
10
4.8
E
-
13
1.5
E
-
13
0.31
Bent 1.6
0.4
4.5
E
-
08
±
4.4
E
-
08
2.2
E
+
10
±
2.8
E
+
09
0.5
3.53
E
-
10
3.07
E
-
10
7.7
E
-
14
5.8
E
-
14
0.76
COX 1
⊥
0.18
2.3
E
-
09
±
3.0
E
-
09
1.6
E
+
10
±
3.8
E
+
09
4.3
6.03
E
-
10
1.11
E
-
09
1.0
E
-
13
3.4
E
-
13
3.4
OPA 1
⊥
0.096
9.6
E
-
11
±
5.7
E
-
10
9.3
E
+
09±1.9E+10
55.0
9.20
E
-
09
9.02
E
-
10
1.2
E
-
11
1.2
E
-
13
0.010Figure 2
Effective diffusion coefficients for dissolved gases in Boom Clay (a) and other clayey materials (b), fitted with expression (14). The diffusion coefficientD is either the diffusion coefficientD0 in pure water or the effective diffusion coefficient Deff in a porous medium.
(a)
(b)Most of the fits of Table6 are fair but show a relatively large error. The OPA fit is not reliable: for both fit parameters, the error is larger than 100%.For all Boom Clay and Eigenbilzen Sands samples, the exponential factornexp is on average (7.9 ± 2) × 10−10 which is similar to the exponential factor n0,exp of the fit of D08.8±1.2×10-9. For the other clayey samples, the factor nexp is always substantially larger than n0,exp.As discussed earlier, for gases, the geometric factorG can be calculated from diffusion experiments. Similar to the evolution of Deff with size; also the evolution of G with size can be described with an exponential or a power law function and this is shown in Figure 3.Figure 3
Geometric factors for dissolved gases in Boom Clay (a) and other clayey materials (b).
(a)
(b)From Figure3 and Table 5, one can observe that for the Boom Clay sample 127b (K4) the geometric factor tends to increase with the kinetic diameter, but the difference in G between the smallest (He) and largest (C2H6) molecule is smaller than a factor 3.For the Eigenbilzen Sands samples, the geometric factor is rather a constant value. This can also be deduced from the exponential factors of the fit ofDeff (Table 6): if n0,exp and nexp are similar (nexp-n0,exp≈0), according to expression (16), the variation of G with molecular size is small, in agreement with the experimental results (see Table 5).For the other clayey samples, the geometric factor clearly increases quite spectacularly with the kinetic diameter: theG ratio between smallest and largest molecule ranges from a factor 3 (bentonite 1.4) up to a factor 16 (Bentonite 1.6). This is also reflected in the exponential factors nexp which differ much more from the n0,exp factor compared to the Boom Clay and Eigenbilzen Sands samples (see Table 6).From the fits in Table6, the tortuosity τ (expression (17)) and Rpore (expression (18)) can be calculated. The calculated tortuosity factors τ are, apart from the unreliable OPA value and the small Bent 1.6 value, quite similar. But the calculated Rpore values are sometimes negative which is physically impossible (Boom Clay sample K2 and the four samples of the Eigenbilzen Sands). For these samples, the values of n0,exp and nexp are within each other’s error interval which means that in some cases n0,exp>nexp, which leads according to expression (18) to negative Rpore values. For these samples, the difference in Deff between the largest and smallest molecule is about one order of magnitude. So a slight change in one of the Deff values can strongly influence nexp. This observation points to a limitation of this approach: it can only be used when nexp differs significantly from n0,exp. From the measured K values, Rpore,K can be calculated according to expression (11) and these results are also shown in Table 6. When comparing both calculated Rpore values, it is clear that when a reliable exponential fit is obtained, and when nexp differs significantly from n0,exp (so no negative Rpore values) both calculated Rpore values are relatively similar.According to expression (6) which relates the geometric factor to the size of the diffusing molecule and the pore size, the geometrical factor always increases as a function of the size of the diffusing molecule. For a sufficiently small molecular size σ (with respect to the pore size Rpore) the exponential factor exp⁡(σ/Rpore) can be approximated by one, so that the geometrical factor is approximately constant. This means that, according to expression (6), the geometric factor for a very small diffusing molecule (kinetic diameter σ→0) is roughly similar in all these clays and the pore size Rpore is not relevant. This applies most to the smallest molecule, helium, for which indeed the variability of G over the different samples is small. From expression (6), one can also conclude that if, for the same sample, all diffusing molecules are sufficiently small compared to the pore size, the geometric factor is also approximately constant. This is the case for the samples of the Eigenbilzen Sands: as Rpore,K (calculated from expression (11)) is around 2 × 10−8m, the pores are much larger than the diffusing molecules.When the diffusing molecules are not small compared to the pore size, the pore size becomes more important. When, in this case, the size of the diffusing molecule increases, the geometric factors in the different clays diverge, and the geometric factor for different gases in a single sample is not constant. For both bentonites and COX, both fittedRpore values are small (see Table 6) (similar to the size of the diffusing molecules) and the geometric factor clearly rises with increasing molecular size (see Figure 3). For Boom Clay, sample K4 which is expected to have larger pores compared to COX and the bentonites, the geometric factor also increases with the size of the diffusing molecule, although less steeply.Despite the enormously simplifying assumption to replace a whole pore distribution by just one pore size in both the hydraulic conductivity model and the effective diffusion coefficient versus diffusing molecule size model, both models lead to similar values for the typical pore size. Using the fittedRpore values to predict the hydraulic conductivity with expression (11) gives quite a good agreement between the calculated hydraulic conductivity Kcalc and the experimental value Kexp for both bentonites, COX and Boom Clay sample K4 (see Table 6). It is straightforward that, for negative Rpore values (sample K2 and the Eigenbilzen Sands), this prediction cannot be good, and it has already been mentioned that the OPA fits are not reliable.The effective diffusion coefficients can also be fitted as a function of size by using the power law expression (19), leading to the optimal values summarized in Table 7. Apart from the fact that for the power law no typical pore size like Rpore exists, the conclusions are similar as for the exponential fit: a bad fit for OPA, a steep increase of the geometric factor as a function of molecular size for both bentonites and COX, a slower increase for both Boom Clay samples, and a hard-to-explain (apart from measurement errors) decrease (npow-n0,pow<0) for the Eigenbilzen Sands (indicated in italic).Table 7
Parameter values obtained by fitting the effective diffusion coefficient versus size of the diffusing molecule by the power law expression (15). The numbers in the top line are taken from [23]. “Tortuosity” τ is defined similarly as expression (17).
Water
D
0
,
p
o
w
n
0
,
p
o
w
(m2/s)
( )
1.1
E
-
07
±
4.6
E
-
08
3.0
±
0.4
Code
η
t
o
t
D
e
f
f
,
p
o
w
n
p
o
w
τ
n
p
o
w
-
n
0
,
p
o
w
(-)
(m2/s)
(m−1)
( )
( )
K2
⊥
0.4
5.0
E
-
09
±
5.8
E
-
09
3.2
±
0.9
8.4
0.2
K4
//
0.38
1.3
E
-
08
±
1.4
E
-
08
3.7
±
0.9
3.1
0.7
K14
//
0.37
4.0
E
-
09
±
3.1
E
-
09
2.2 ± 0.6
10.0
−0.8
K15
//
0.4
5.9
E
-
09
±
4.8
E
-
09
2.4 ± 0.7
6.6
−0.6
K16
⊥
0.41
4.8
E
-
09
±
5.5
E
-
09
2.3 ± 0.9
8.3
−0.7
K17
⊥
0.4
7.0
E
-
09
±
6.8
E
-
09
2.9 ± 0.8
5.5
−0.2
Bent 1.4
0.47
1.8
E
-
08
±
7.2
E
-
09
5.0
±
0.3
2.8
1.9
Bent 1.6
0.4
2.3
E
-
07
±
1.9
E
-
07
7.5
±
0.7
0.2
4.5
COX 1
⊥
0.18
1.0
E
-
08
±
1.3
E
-
08
5.8
±
1.1
1.9
2.8
OPA 1
⊥
0.096
1.1
E
-
10
±
7.1
E
-
10
2.7
±
5.9
94.2
−0.4


## 4.3.1. HTO Diffusion Results

As described above, the HTO experiments are fitted by using fixed filter (FF) values. The results are shown in Table3. Please note that the value for HTO is now obtained with fixed filter values, while the OPA value in Jacops et al. [22] was obtained considering filter = clay.For the capacity factors, there is for some samples a considerable difference between the fitted capacity factors and the total porosityηtot (for this parameter, sometimes a reference value is used, and sometimes it is measured by water content), in particular for OPA. However, the error on the fits is small (Table 4); hence the fits are considered to be good.


## 4.3.2. Diffusion Coefficients and Geometric Factors of Dissolved Gases

All measured effective diffusion coefficients and the corresponding geometrical factors for dissolved gases and HTO are reported in Tables4 and 5. The COX and OPA (He and Ar) results have already been reported in [22], as well as the Boom Clay samples ON-Mol1 84b and 127b in [23]. Please note that the fits for COX were performed with a porosity of 18%, in order to be consistent with Jacops et al. [22].Table 5
Geometrical factors obtained from the aqueous and effective diffusion coefficients from Table4. Note that similar to the previous table, the reference value of 0.18 is used for the COX1 total porosity.
Core
Code
η
t
o
t (-)
He
HTO
Ne
H2
Ar
CH4
Xe
C2H6
G (-)
BC 84b
K2
⊥
0.4
6.2
4.7
9.2
14.1
7.6
9.6
11.9
BC 127b
K4
//
0.38
3.7
3.0
6.7
3.8
6.4
4.5
8.4
8.9
BC 36a
K14
//
0.37
4.0
2.5
4.4
2.9
3.6
BC 37b
K15
//
0.4
3.5
2.0
4.0
2.6
3.6
BC 35b
K16
⊥
0.41
4.0
2.7
4.7
2.0
4.1
BC 39b
K17
⊥
0.4
5.0
2.9
5.3
3.0
6.2
Bent 1.4
Bent 1.4
0.47
19.3
6.6
18.4
30.7
60.3
Bent 1.6
Bent 1.6
0.4
10.8
10.9
18.5
79.9
161.7
COX 1
COX 1
⊥
0.18
16.1
34.6
60.2
131.1
OPA 1
OPA 1
⊥
0.096
10.2
17.5
62.1
54.0Opalinus Clay has the lowest porosity of all the samples used in this paper. It is evidently the first candidate showing possible percolation transition problems. These occurred for all gases larger than Argon: CH4, Xe, and C2H6:(i)
For CH4 the measurements were scattered and all below 100 ppm. As 100 ppm is the lower limit for reliable CH4 peak detection and measurement, all measurements contain a significant portion of noise and are therefore considered nonsignificant.(ii)
For Xe, no breakthrough was measured which could also be related to the detection limit of 100 ppm.(iii)
For C2H6, where the detection limit is only 5 ppm, a clear breakthrough curve was detected, but it could not be fitted by a simple diffusion model (see Figure 1): a quasi-stationary state is reached very soon, but after about 150–200 days, the outlet concentration starts to increase much faster. Because the system is probably very close to the percolation threshold, a possible explanation might be diffusion over another percolating, longer pathway.Figure 1
Evolution of the outlet concentration versus time in the C2H6 diffusion experiment in Opalinus Clay. A clear breakthrough after a couple of days, leading to a quasi-stationary state is followed by an sudden unexplained increase after more than 150 days.When comparing the measured diffusion coefficients in the Boom Clay samples and the samples from the Eigenbilzen Sands, one can observe that despite a comparable porosity, the diffusion coefficients in the Eigenbilzen Sands are slightly higher than in Boom Clay. This agrees with previous results from pulse injection experiments [59], showing on average higher HTO dispersion coefficients in the Eigenbilzen Sands than in Boom Clay. Evidently, the present results confirm that diffusion coefficients are anisotropic and higher for cores parallel to the bedding compared to cores perpendicular to the bedding. Despite its higher porosity compared to the Boom Clay and the Eigenbilzen Sands, the diffusion coefficients in bentonite 1.4 are smaller (as was also the case for the hydraulic conductivity). For helium in bentonite, we observed that diffusion coefficient in bentonite 1.4 is considerably lower than for bentonite 1.6. In fact, one would expect the opposite and we cannot explain this observation. The other three gases show the expected behaviour: lower diffusion coefficients in bentonite 1.6 compared to bentonite 1.4. The lowest diffusion coefficients were observed for the clays with the lowest porosity: COX and OPA.Also, note that, for the samples from the Eigenbilzen Sands and for the Boom Clay samples, the diffusion coefficients for HTO and Ne (two molecules with a nearly equal kinetic diameter) are about the same. Due the very different values for the HTO and Ne aqueous diffusion coefficients (D0), the corresponding geometric factors differ considerably.


## 4.3.3. The Effective Diffusion Coefficients and the Geometric Factor as a Function of the Size of the Diffusing Molecule

Table6 shows the fitted values according to expression (14) from the measured effective gas diffusion coefficients and the corresponding graphs are shown in Figure 2. Although HTO is included in Table 4, it is not included in the fits (shown in Table 6) as we only describe the data for dissolved gases.Table 6
Parameter values obtained by fitting the effective diffusion coefficient versus molecular size by the exponential expression (14). The numbers in the top line are taken from Jacops et al. [23]. Rpore is calculated from expression (18) and Rpore,K is calculated from expression (11). From the fitted Rpore value, the hydraulic conductivity is calculated (Kcalc) and compared to the experimental value Kexp.
Water
D
0
,
e
x
p
n
0
,
e
x
p
(m2/s)
(m−1)
5.5
E
-
08
±
2.3
E
-
08
8.8
E
+
09
±
1.2
E
+
09
Code
η
t
o
t
D
e
f
f
,
e
x
p
n
e
x
p
τ
R
p
o
r
e
R
p
o
r
e
,
K
K
c
a
l
c
K
e
x
p
K
e
x
p
/
K
c
a
l
c
(-)
(m2/s)
(m−1)
( )
(m)
(m)
(m/s)
(m/s)
K2
⊥
0.4
2.0
E
-
09
±
2.1
E
-
09
8.7E + 09 ± 3.0E + 09
11.1
−1.53E − 07
2.32
E
-
09
1.4
E
-
08
3.3
E
-
12
0.0002
K4
//
0.38
5.7
E
-
09
±
5.4
E
-
09
1.1
E
+
10
±
2.7
E
+
09
3.7
2.51
E
-
09
2.41
E
-
09
3.7
E
-
12
3.4
E
-
12
0.92
K14
//
0.37
2.4
E
-
09
±
1.6
E
-
09
6.4E + 09 ± 1.9E + 09
8.4
−1.91E − 09
2.78
E
-
08
2.1
E
-
12
4.4
E
-
10
212
K15
//
0.4
3.4
E
-
09
±
2.4
E
-
09
6.9E + 09 ± 2.0E + 09
6.4
−2.51E − 09
1.34
E
-
08
3.9
E
-
12
1.1
E
-
10
28
K16
0.41
2.8
E
-
09
±
2.7
E
-
09
6.6E + 09 ± 2.7E + 09
8.0
−2.08E − 09
1.38
E
-
07
2.7
E
-
12
1.2
E
-
08
4397
K17
⊥
0.4
3.7
E
-
09
±
2.9
E
-
09
8.2E + 09 ± 2.2E + 09
5.9
−8.21E − 09
2.61
E
-
08
4.1
E
-
11
4.2
E
-
10
10.1
Bent 1.4
0.47
6.4
E
-
09
±
2.0
E
-
09
1.4
E
+
10
±
8.7
E
+
08
4.0
8.17
E
-
10
4.56
E
-
10
4.8
E
-
13
1.5
E
-
13
0.31
Bent 1.6
0.4
4.5
E
-
08
±
4.4
E
-
08
2.2
E
+
10
±
2.8
E
+
09
0.5
3.53
E
-
10
3.07
E
-
10
7.7
E
-
14
5.8
E
-
14
0.76
COX 1
⊥
0.18
2.3
E
-
09
±
3.0
E
-
09
1.6
E
+
10
±
3.8
E
+
09
4.3
6.03
E
-
10
1.11
E
-
09
1.0
E
-
13
3.4
E
-
13
3.4
OPA 1
⊥
0.096
9.6
E
-
11
±
5.7
E
-
10
9.3
E
+
09±1.9E+10
55.0
9.20
E
-
09
9.02
E
-
10
1.2
E
-
11
1.2
E
-
13
0.010Figure 2
Effective diffusion coefficients for dissolved gases in Boom Clay (a) and other clayey materials (b), fitted with expression (14). The diffusion coefficientD is either the diffusion coefficientD0 in pure water or the effective diffusion coefficient Deff in a porous medium.
(a)
(b)Most of the fits of Table6 are fair but show a relatively large error. The OPA fit is not reliable: for both fit parameters, the error is larger than 100%.For all Boom Clay and Eigenbilzen Sands samples, the exponential factornexp is on average (7.9 ± 2) × 10−10 which is similar to the exponential factor n0,exp of the fit of D08.8±1.2×10-9. For the other clayey samples, the factor nexp is always substantially larger than n0,exp.As discussed earlier, for gases, the geometric factorG can be calculated from diffusion experiments. Similar to the evolution of Deff with size; also the evolution of G with size can be described with an exponential or a power law function and this is shown in Figure 3.Figure 3
Geometric factors for dissolved gases in Boom Clay (a) and other clayey materials (b).
(a)
(b)From Figure3 and Table 5, one can observe that for the Boom Clay sample 127b (K4) the geometric factor tends to increase with the kinetic diameter, but the difference in G between the smallest (He) and largest (C2H6) molecule is smaller than a factor 3.For the Eigenbilzen Sands samples, the geometric factor is rather a constant value. This can also be deduced from the exponential factors of the fit ofDeff (Table 6): if n0,exp and nexp are similar (nexp-n0,exp≈0), according to expression (16), the variation of G with molecular size is small, in agreement with the experimental results (see Table 5).For the other clayey samples, the geometric factor clearly increases quite spectacularly with the kinetic diameter: theG ratio between smallest and largest molecule ranges from a factor 3 (bentonite 1.4) up to a factor 16 (Bentonite 1.6). This is also reflected in the exponential factors nexp which differ much more from the n0,exp factor compared to the Boom Clay and Eigenbilzen Sands samples (see Table 6).From the fits in Table6, the tortuosity τ (expression (17)) and Rpore (expression (18)) can be calculated. The calculated tortuosity factors τ are, apart from the unreliable OPA value and the small Bent 1.6 value, quite similar. But the calculated Rpore values are sometimes negative which is physically impossible (Boom Clay sample K2 and the four samples of the Eigenbilzen Sands). For these samples, the values of n0,exp and nexp are within each other’s error interval which means that in some cases n0,exp>nexp, which leads according to expression (18) to negative Rpore values. For these samples, the difference in Deff between the largest and smallest molecule is about one order of magnitude. So a slight change in one of the Deff values can strongly influence nexp. This observation points to a limitation of this approach: it can only be used when nexp differs significantly from n0,exp. From the measured K values, Rpore,K can be calculated according to expression (11) and these results are also shown in Table 6. When comparing both calculated Rpore values, it is clear that when a reliable exponential fit is obtained, and when nexp differs significantly from n0,exp (so no negative Rpore values) both calculated Rpore values are relatively similar.According to expression (6) which relates the geometric factor to the size of the diffusing molecule and the pore size, the geometrical factor always increases as a function of the size of the diffusing molecule. For a sufficiently small molecular size σ (with respect to the pore size Rpore) the exponential factor exp⁡(σ/Rpore) can be approximated by one, so that the geometrical factor is approximately constant. This means that, according to expression (6), the geometric factor for a very small diffusing molecule (kinetic diameter σ→0) is roughly similar in all these clays and the pore size Rpore is not relevant. This applies most to the smallest molecule, helium, for which indeed the variability of G over the different samples is small. From expression (6), one can also conclude that if, for the same sample, all diffusing molecules are sufficiently small compared to the pore size, the geometric factor is also approximately constant. This is the case for the samples of the Eigenbilzen Sands: as Rpore,K (calculated from expression (11)) is around 2 × 10−8m, the pores are much larger than the diffusing molecules.When the diffusing molecules are not small compared to the pore size, the pore size becomes more important. When, in this case, the size of the diffusing molecule increases, the geometric factors in the different clays diverge, and the geometric factor for different gases in a single sample is not constant. For both bentonites and COX, both fittedRpore values are small (see Table 6) (similar to the size of the diffusing molecules) and the geometric factor clearly rises with increasing molecular size (see Figure 3). For Boom Clay, sample K4 which is expected to have larger pores compared to COX and the bentonites, the geometric factor also increases with the size of the diffusing molecule, although less steeply.Despite the enormously simplifying assumption to replace a whole pore distribution by just one pore size in both the hydraulic conductivity model and the effective diffusion coefficient versus diffusing molecule size model, both models lead to similar values for the typical pore size. Using the fittedRpore values to predict the hydraulic conductivity with expression (11) gives quite a good agreement between the calculated hydraulic conductivity Kcalc and the experimental value Kexp for both bentonites, COX and Boom Clay sample K4 (see Table 6). It is straightforward that, for negative Rpore values (sample K2 and the Eigenbilzen Sands), this prediction cannot be good, and it has already been mentioned that the OPA fits are not reliable.The effective diffusion coefficients can also be fitted as a function of size by using the power law expression (19), leading to the optimal values summarized in Table 7. Apart from the fact that for the power law no typical pore size like Rpore exists, the conclusions are similar as for the exponential fit: a bad fit for OPA, a steep increase of the geometric factor as a function of molecular size for both bentonites and COX, a slower increase for both Boom Clay samples, and a hard-to-explain (apart from measurement errors) decrease (npow-n0,pow<0) for the Eigenbilzen Sands (indicated in italic).Table 7
Parameter values obtained by fitting the effective diffusion coefficient versus size of the diffusing molecule by the power law expression (15). The numbers in the top line are taken from [23]. “Tortuosity” τ is defined similarly as expression (17).
Water
D
0
,
p
o
w
n
0
,
p
o
w
(m2/s)
( )
1.1
E
-
07
±
4.6
E
-
08
3.0
±
0.4
Code
η
t
o
t
D
e
f
f
,
p
o
w
n
p
o
w
τ
n
p
o
w
-
n
0
,
p
o
w
(-)
(m2/s)
(m−1)
( )
( )
K2
⊥
0.4
5.0
E
-
09
±
5.8
E
-
09
3.2
±
0.9
8.4
0.2
K4
//
0.38
1.3
E
-
08
±
1.4
E
-
08
3.7
±
0.9
3.1
0.7
K14
//
0.37
4.0
E
-
09
±
3.1
E
-
09
2.2 ± 0.6
10.0
−0.8
K15
//
0.4
5.9
E
-
09
±
4.8
E
-
09
2.4 ± 0.7
6.6
−0.6
K16
⊥
0.41
4.8
E
-
09
±
5.5
E
-
09
2.3 ± 0.9
8.3
−0.7
K17
⊥
0.4
7.0
E
-
09
±
6.8
E
-
09
2.9 ± 0.8
5.5
−0.2
Bent 1.4
0.47
1.8
E
-
08
±
7.2
E
-
09
5.0
±
0.3
2.8
1.9
Bent 1.6
0.4
2.3
E
-
07
±
1.9
E
-
07
7.5
±
0.7
0.2
4.5
COX 1
⊥
0.18
1.0
E
-
08
±
1.3
E
-
08
5.8
±
1.1
1.9
2.8
OPA 1
⊥
0.096
1.1
E
-
10
±
7.1
E
-
10
2.7
±
5.9
94.2
−0.4


## 4.4. How to Estimate Gas Diffusion Coefficients as a Function of the Size of the Diffusing Molecule?

Based on the results presented in this work, one can state that the approach of using one geometric factor for estimatingDeff of species with a different molecular size is not always correct and the validity depends on the type of sample and the size of the diffusing molecule. For instance, for the Eigenbilzen Sand samples, the geometric factor from HTO could be used to calculate Deff for, for example, C2H6. On the contrary, for COX or bentonite, using the geometric factor of HTO to calculate Deff for C2H6 would lead to a serious overestimation. In order to illustrate this, we use for COX a geometric factor of 17 for HTO, as suggested by Jacops et al. [22]. If this value would be used to predict Deff for C2H6, a value of 1.5 × 10−11 m2/s would be obtained while the measured value was 0.2 × 10−11 m2/s so this leads to an overestimation of factor 7.A more correct way to estimate diffusion coefficients would be by using the exponential or power law relation betweenDeff and the molecular size (expressions (14) and (15)). Each expression has two parameters, meaning that for two gases the effective diffusion coefficient needs to be determined experimentally. Evidently, it is best to do this for a large molecule (e.g., C2H6) and a small one. He is the smallest of all molecules considered here, but its diffusion coefficients are always higher than expected (located above the fitted curve in Figure 2). This leads to the observation that the measured Deff is on average 1.7 times larger than the calculated Deff (see Table 8). For the second smallest molecule (neon), the calculated average Deff is similar to the measured Deff; hence Ne is recommended for being used in the fits. An example of such a prediction with the exponential expression (14) is shown in Table 9. Evidently, because Ne and C2H6 are used to determine the fit parameters, for both of these gases, the predicted value corresponds exactly to the experimental value. For the other gases, the prediction can be improved by multiplying the predicted diffusion coefficient with the averages of the ratio of the experimental to fitted value for all measured clays (see Table 8). This manipulation allows correcting for the specific bias of this ratio, which is different for each gas. After this correction, for sample K4 (one of the experiments for which most measurements were carried out), apart from He, the predicted and measured effective diffusion coefficients differ up to less than 30% (see Table 9).Table 8
Ratios between the experimentally determinedDeff (values see Table 4) and the fitted Deff (using the fits given in Table 6), as well as the averages for each element. Due to the deviating value of He for sample OPA1, the corresponding Deff is not included in the average.
Core
Code
He
Ne
H2
Ar
CH4
Xe
C2H6
Ratio Exp/Fit (-)
BC 84b
K2
2.3
1.0
0.7
1.4
1.1
1.1
BC 127b
K4
2.0
0.8
2.1
1.0
1.6
0.9
1.1
BC 36a
K14
1.4
0.8
1.1
1.0
BC 37b
K15
1.4
0.8
1.2
1.0
BC 35b
K16
1.4
0.8
1.6
0.9
BC 39b
K17
1.3
0.8
1.5
0.9
Bent 1.4
Bent 1.4
1.1
0.9
1.1
1.0
Bent 1.6
Bent 1.6
1.7
0.8
0.8
1.2
COX 1
COX 1
2.4
0.9
0.9
1.2
OPA 1
OPA 1
7.8
0.9
1.1
Average
1.7
0.8
2.1
0.9
1.3
1.0
1.0
Standard deviation
0.4
0.1
0.2
0.3
0.2
0.1
For
∗ He, OPA1 not included.Table 9
Predicting the gas effective diffusion coefficients of sample K4 after measuring the Ne and C2H6 (both in italic) effective diffusion coefficients. In a first phase, the prediction is based only on expression (14). Subsequently, this prediction is improved by multiplication with a “correction factor” (see Table 8), which takes into account the gas specific deviation between the experimental and the predicted value.
He
Ne
H2
Ar
CH4
Xe
C
2
H
6
Kinetic diameterσ
(1E-10 m)
2.58
2.79
2.97
3.42
3.82
4.06
4.42
D
e
f
f (experiment)
(1E-11 m2/s)
74.7
22.9
51.2
14.5
15.5
6.6
5.9
D
e
f
f (prediction)
(
1
E
-
11 m2/s)
27.2
22.9
19.7
13.6
9.7
8.0
5.9
RatioDeff (exp)/Deff (pred)
( )
2.7
1.0
2.6
1.1
1.6
0.8
1.0
Correction factor
( )
1.7
0.8
2.1
0.9
1.3
1.0
1.0
D
e
f
f (pred) ∗ Corr factor
(1E-11 m2/s)
45.8
41.2
12.2
12.4
7.7
RatioDeff (exp)/(Deff (pred) ∗ Corr factor)
( )
1.6
1.2
1.2
1.2
0.9By using this approach, the diffusion coefficient of, for instance, H2 (which is difficult to measure due to microbial activity [24]), can be estimated based on its molecular size. But note that the diffusion coefficient for H2 was (due to experimental difficulties) measured for only one sample; hence the correction factor might be only an estimate.


## 5. Conclusion

Summarizing, for unreactive gases in clayey materials, the effective diffusion coefficient varies with the size of the diffusing molecule and this variation can be described by an exponential or a power law function. Besides, the geometric factor can be described by a decreasing exponential function of the ratio between the size of the diffusing molecule and a characteristic pore sizeRpore. This means that, in case of a large characteristic pore size (what corresponds to a higher hydraulic conductivity) and a small diffusing molecule, the geometric factor (being the ratio of the aqueous diffusion coefficient and the apparent diffusion coefficient) is approximately constant or slowly increasing as a function of the size of the diffusing molecule. For an increasing size of the diffusing molecule, the increase of the geometric factor is expected to become more steep, proportional to the increase of the path length caused by the nonavailability of pore sizes smaller than the molecular size. Such a steep increase is observed for small diffusing molecules in clays with a small characteristic pore size (corresponding to a low hydraulic conductivity). If the size of the diffusing molecule still increases, there is no longer a path of sufficiently large pores connecting both sides of the clay, and the diffusion coefficient becomes zero. Experimental problems due to such a percolation threshold were observed in Opalinus Clay for gas molecules larger than argon. For samples of the Eigenbilzen Sands, we observed that an increase in sand content results in an increased hydraulic conductivity, an increase in diffusivity, and less variability ofG with the size of the diffusing molecule.By measuring experimentally the effective (or apparent) diffusion coefficient of two unreactive gases in clay, the remaining diffusion coefficients can be estimated by an exponential (or alternatively a power law) relationship.



---

*Source: 8652560-2017-12-27.xml*