An experimental technique has been developed and applied to the problem of determining effective diffusion coefficients and partition coefficients of Sr in low permeability geological materials. This technique, the micro-reactor simulated channel method (MRSC), allows rapid determination of contaminant transport parameters with resulting values comparable to those determined by more traditional methods and also creates product surfaces that are amenable for direct chemical analysis. An attempt to further constrain mass flux was completed by detailed ion beam analysis of polished tuff surfaces (tuff is a polycrystalline polyminerallic aggregate dominated by silicate phases) that had been reacted with Sr solutions at concentrations of 10−5, 10−3 and 10−1 mol l−1. Ion beam analysis was carried out using beams of both protons (using particle induced X-ray emission and elastic backscattering spectrometry or EBS) and alpha-particles (using Rutherford backscattering spectrometry). The ion beam analyses showed that increased solution concentrations resulted in increased surface concentrations and that in the highest concentration experiment, Sr penetrated to at least 4 μm below the primary interface. The Sr surface concentrations determined by EBS were 0.06 (±0.05), 0.87 (±0.30) and 2.40 (±1.0) atomic weight % in the experiments with starting solution concentrations of 10−5, 10−3, and 10−1 mol l−1, respectively.


Retention of radionuclides by geological materials is an important but often poorly constrained process (Seitz et al., 1987). Measuring diffusive flux in fluid–rock systems is difficult due to the extremely slow rates of diffusion and generally low solute concentrations. However, diffusive flux is expected to be an important mass transport process in radioactive waste repositories and therefore estimates of the diffusive flux of contaminant radionuclides are critical. Thus two obstacles exist in making such measurements: (1) conventional experiments need to be run over long periods of time under stable conditions; and (2) the analysis of adsorbed or precipitated contaminants must be able to determine low concentrations ideally with spatial resolution.

Interaction between the rock matrix and dissolved radionuclides can be characterized using the effective diffusion coefficient (De) and the distribution coefficient (Kd). Previous laboratory scale experiments have been completed along with a theoretical analysis of the results in order to quantify the values of these coefficients (Skagius and Neretnieks, 1988). Despite their utility, certain limitations exist in the case of standard lab-scale diffusivity determinations. Firstly, the conventional technique is relatively time-consuming and months are needed to complete the measurement, particularly for strongly sorbing nuclides (Skagius and Neretnieks, 1982; Tsukamoto and Ohe, 1991). Long elapse times cause difficulties in maintaining constant experimental conditions and thereby create the risk of large uncertainties. In addition, different sample geometries are typically used depending on the parameter to be determined, e.g. crushed samples may be used for determining Kd whereas intact samples may be used for determining De. The dependence of both Kd (Rancon, 1986) and De (Tsukamoto and Ohe, 1991) on size has been observed and therefore an additional correction due to differing sample geometry is required to maintain consistency between parameters.

Therefore a fast method has been developed, hereafter referred to as the micro-reactor simulated-channel method (MRSC), to determine both De and Kd simultaneously using non-crushed rock samples (Okuyama et al., 2008). The concept of this method is similar to that of the microchemical reactor, which has a very thin fluid flow channel (10–200 μm) and has the advantage of increasing the rate of surface reactions because of the high surface area to liquid volume ratio. This feature enables us to make very fast determinations even in a small reactor volume.

In the experiments, an aqueous fluid containing dissolved Sr has been pumped through a simulated rock fracture. Strontium is an expected radioactive contaminant in many nuclear waste repository scenarios, in particular the 90Sr isotope is the radioisotope of concern and has a half-life of 28.90 years. For convenience in these experiments we have used non-radioactive Sr. A solution bearing Sr as a proxy for radioactive contaminants has been forced through the simulated channel (or fracture) created between two polycrystalline wafers. Solute concentration was monitored as a function of time; Sr concentrations in the effluent were decreased relative to the inlet solution, and hence a quantity of Sr remained attached to the solids within the reactor. Concentrations remaining behind are much higher than would be predicted from simple surface adsorption onto the exposed wafers, and hence it is postulated that some Sr has diffused into the boundaries between the grains that compose the polycrystalline wafers.

Therefore in order to determine the Sr inventory in the near surface of the solid interface, we developed an analytical regimen that included ion beam techniques. In particular, this approach was developed for the following two reasons. (1) It allowed us to take advantage of the extremely low background radiation levels at the SrKα energy in particle induced X-ray emission (PIXE) analysis in order to accurately quantify Sr taken up by mineral surfaces in the micro-reactor experiments (Wogelius et al., 1992, 1997). (2) We were able to use the depth profiling capability of elastic (non-Rutherford) backscattering spectrometry (EBS) and Rutherford backscattering spectrometry (RBS) to non-destructively constrain diffusion profiles below the surface. The distinction between RBS and EBS is solely in the applicable scattering cross-sections, which for RBS can be derived analytically from the Coulomb potential (Gurbich, 2010).

Our detailed approach involved preliminary analysis by environmental scanning electron microscope (ESEM) to determine grain morphology and check that levels of Sr were indeed enhanced after reaction. Samples were then to be analysed via PIXE with simultaneous EBS measurements. The EBS measurement was then repeated using alpha particles (RBS) in an attempt to give a fully constrained concentration depth profile with improved elemental resolution. Results were intended to give constraints on the grain boundary diffusion coefficients and total mass transfer into the solid matrix under our experimental conditions.

Experimental materials and methods

Micro-reactor with simulated channel

The new apparatus consists of an injection syringe pump, a reaction unit, an auto sampler, and a drainage tank (Fig. 1a). Rock samples were a natural tuff known as Itaya clinoptilolite (Iijima, 1980). All specimens were cut parallelepipeds, with dimensions of 30 × 20 × 3 mm, and the surfaces smoothed with abrasive paper (#200; Sankyo Rikagaku). The side edge surfaces were then coated with an epoxy resin to avoid evaporation of the test solution during experiments. The micro fluid channel was constructed using a fluoroplastic spacer. The reactor was connected with inlet and outlet ports (diameter 0.5 mm) for the injection or drainage of strontium-bearing solutions or distilled water. A test without the rock sample in place was carried out to verify that the surfaces of the fluoroplastic spacers and base plates were non-reactive.

Intact rock samples were cut and sized, then water-saturated in distilled water under a vacuum for 24 h prior to the experiments. A series of experiments was carried out at a constant flow rate of 8.3 × 10−11 m3 s−1. The accuracy of the measurements was better than 5%. The experimental apparatus was assembled in a glove box and all the experiments were performed at 25–2°C. Diffusive flux is perpendicular (normal) to the fluid rock interface whereas advective flux is along the axis of the channel, normal to diffusion. Figure 1b indicates how samples were prepared for ion beam analysis. The inlet side was marked with a black dot and the sample was cut longitudinally. Analyses presented here were of the top surface, that is, the surface that had been directly in contact with the aqueous solution.

Ion beam analysis analytical details

The PIXE and RBS analyses were completed at the Ion Beam Centre, University of Surrey. A 2 MV tandem accelerator (Simon et al., 2004) was used to produce the incident proton beam. The PIXE and EBS(H) used a beam of 2.5 MeV protons operating at 0.2 nA current. The resulting spectra were measured simultaneously for each area analysed. Detection of X-rays was by energy-dispersive spectroscopy (EDS) via a Si(Li) detector at either 30 or 60 mm distance with a 130 μm Be filter placed between sample and detector to suppress low energy background. The sample-normal to detector angle was 45° and the beam spot was approximately 3 × 3 μm. The PIXE precision was ~1% of the reported value except for Al, where the precision was ~10% (unless otherwise noted). In the case of the backscattered proton signal, the detector surface area was 150 mm2 and the detector distance was 50 mm, giving a solid angle of 60 mSr. The scattering angle was 150°. Subsequent RBS(α) measurements used a 1.554 MeV beam using two detectors with scattering angles of 148.2° (‘IBM’ geometry) and 172.8° (‘Cornell’ geometry) and solid angles of 6.5 and 1.25 mSr, respectively. The beam current was ~10 nA and the nominal beam size (normal incidence) was 1 mm. The RBS errors in both cases are ~5% for Si, O, and C; ~20% for Ca, Fe, Al; ~40% for Sr, Ba, K, and Na (unless otherwise noted).

Samples were attached to aluminium stubs by the epoxy enclosures and analysed under vacuum (1 × 10−5 atm). The total charge (Q) and detector dead time were recorded and used in the calculation algorithms. The PIXE/EBS analysis was done using Dan32 (Grime, 1996) and used the GUPIX code as a PIXE calculator (Blaauw et al., 2002). The RBS analysis was done using the DataFurnace code (Jeynes et al., 2003). This practice of analysing multiple IBA datasets self-consistently has been recently reviewed in Jeynes et al. (2012).

Environmental scanning electron microscopy and X-ray diffraction

The energy dispersive spectrometry (EDS) analysis was performed with an ESEM using 20 kV accelerating voltage, 32° take off angle and 20% deadtime, with 0.2 Torr pressure in the column. The XRD was completed in a standard powder diffractometer using CuKα radiation and a scan speed of 1° per minute with an angular resolution of 0.02°.

Starting material

The dominant mineral phases in the material were clinoptilolite (44–48 wt.%) and mordenite (10–43 wt.%) with quartz, opal and plagioclase also reported (Takasaka et al.,1989). Figure 2 shows an XRD pattern of the starting material along with a comparison to two standard reference powder patterns. An excellent match exists between the clinoptilolite reference pattern and the sample used in these experiments. The stoichiometry of the clinoptilolite reference patterns is Na4.24K0.28Ca1.8Mg0.1 6Al8.16Si27.84O72·24.88 H2O (79-1461). The reference pattern stoichiometries can be compared to the ideal stoichiometry as reported in Deer et al. (1992): (Na,K)6[Al6Si30O72]·24H2O. The resultant RBS spectra also allowed us to calculate stoichiometries for comparison to these reference materials as presented below. In addition the XRD data indicated that a significant amount of quartz was present in the starting material (accounting for the peak at 55° 2θ, d = 1.6695 Å).

Reactant fluids

Three experiments were completed at different dissolved Sr concentrations. We designate these simply as: (1) sample A, reacted with 10−5 mol l−1 Sr solution; (2) sample B, reacted with 10−3 mol l−1 Sr solution; and (3) sample C, reacted with 10−1 mol l−1 Sr solution.

Two different samples were prepared for sample B: one for the MRSC experiment alone (sample B1) and the other both for the MRSC experiment and later beam analyses (sample B2). A stock solution of Sr was prepared by dissolving reagent grade SrCl2 into deionized water. The solution pH was then adjusted to 6.0 by addition of either HCl or NaOH. Deuterated water (30 wt.%) as a non-sorbing solute tracer was included for sample B1. The deuterium concentrations in the outlet were determined by visible light absorption spectroscopy as a function of time (Okuyama et al., 2008).

Numerical analysis

In the following calculations, both advection and diffusion in the micro-channel were considered but only diffusion was considered for the rock matrix because its hydraulic conductivity is very small compared to that of the micro fluid channel. The micro channel was modelled as a hollow cylindrical tube in order to apply a radially symmetric computational approach, and the rock sample was regarded as a hollow cylinder mounted on this tube. The equation of the system was expressed as the one-dimensional diffusion-advection equation for the micro fluid channel coupled with the matrix diffusion equation in two-dimensional cylindrical coordinate system. The governing equations were thus formulated as follows (Sudicky and Frind, 1982).

For the micro flow channel (0xL):  
For the rock matrix (Rar≤Rb, 0≤x≤L):  
where C(x,t) is the Sr concentration in the micro fluid channel (mol m−3); Cm(x,r,t) is the Sr concentration in the pore water of the rock matrix (mol m−3); r is the radial distance from the surface of the micro channel (m); x is the longitudinal distance from the injection point (m); t is the elapsed time (s); UX is the flow velocity in the micro channel (m s−1); DL is the longitudinal dispersion coefficient (m2 s−1); De is the matrix diffusion coefficient (m2 s−1); ε is the porosity of the rock (dimensionless); ρd is the density of the rock (kg m−3); b is the depth of the micro channel (m); Kd is the sorption coefficient onto the rock matrix (m3 kg−1); Kc is the surface sorption coefficient (m3 m−2); L is the length of the rock sample and the micro channel (m); H is the height of the rock sample (m); W is the width of the rock sample (m); Ra is the inner radius of the hollow cylindrical rock sample (m); and Rb is the outer radius of the hollow cylindrical rock sample (m).
The initial concentration conditions were set equal to zero and a constant concentration boundary condition at the inlet was assumed. The conditions at both right and left ends of the sample were assumed to be a reflection boundary and then equations 1 and 2 are solved by using the finite control volume method (Patanker, 1980). A Langmuir sorption isotherm for Sr was assumed as has been used for other clinoptilolites (Li et al., 2002), thus the sorption coefficient may be derived by equation 3:  
where q is the sorption amount (mol kg−1); C is the liquid phase concentration (mol m−3) q is the maximum sorption amount (mol kg−1); and KL is the Langmuir sorption coefficient (m3 kg−1).

The maximum sorption amount (q) was set equal to the half of the cation exchange capacity (CEC) of a high grade sample (Takasaka et al., 1989) because of the impurities found in the XRD analysis. The Langmuir sorption coefficient (KL) was determined from the results of the batch sorption experiments, which gave a Kd value of 58.7 (m3 kg−1). This figure is similar to that of another Japanese tuff (Akiba and Hashimoto, 1990). From this value, KL was simply calculated from equation 3 as 78.5 m3 mol−1.

Results and discussion

Breakthrough curves

The breakthrough curves both of D2O and Sr are shown in Fig. 3. Remarkably different breakthrough patterns were observed because of the different sorption characteristics of Sr2+ and D2O on the sample. The effective diffusion coefficient in the rock matrix (De) of D2O given by least squares fitting was 1.2 × 10−10 m2 s−1.

Papelis and Um (2003) reported the De of non-sorbing species Br in a tuffaceous core sample with similar porosity of 28.8% to be 1.75 × 10−10 m2 s−1. This value is quite similar to that of the present analysis for D2O. However, it is difficult to obtain a reliable De value for Sr from the breakthrough curve because concentrations are low (near the detection limit). Instead, the diffusion coefficient of Sr2+ was calculated using equation 4 (van Brakel and Heertjes, 1974) which considers the geometric effect of the pores in the matrix. The calculated breakthrough curve of D2O is depicted in Fig. 3 and the value of F is determined as 0.064 by equation 4:  
where F is the formation factor; Dν is the diffusion coefficient of Sr2+ in free water; De,D is the effective diffusion coefficient of D2O; and Dν,D is the self diffusion coefficient of D2O.

The properties used for the calculation of the breakthrough curves are summarized in Table 1. The matrix diffusion coefficient for Sr in this tuff is approximately a factor of two smaller than that for D2O.

Analytical results

The ESEM maps and EDS analysis were completed on sample C showing that most grains were under 5 μm in length (in their longest dimension). Along with the textural information, elemental mapping shows most regions have a constant Al/Si ratio of approximately ¼, consistent with clinoptilolite [different from quartz (0) or feldspar (⅓ to 1)], thus corroborating the XRD analysis. Compositional heterogeneity in the surface is quite high, however, and Fig. 4 shows a silicon map with small bright spots. These bright spots give an EDS pattern consistent with a pure SiO2 phase which has been identified (and corroborated) via XRD analysis of the unreacted samples as quartz. The Sr was difficult to detect via ESEM-EDS.

Particle induced X-ray emission

The PIXE analyses have an approximately 103 decrease in background relative to electron microprobe analysis as well as improved spatial resolution. The EBS with a proton beam (H-EBS) can be completed simultaneously with PIXE and allows depth resolution for flat surfaces plus direct analysis for light atoms such as carbon and oxygen. This means that mineral stoichiometry, metal concentration and element depth profiles can be determined at the same time and from the same point, line or area in a sample. This is of particular utility in studying trace-element profiles because it means that the depth profile of a specific component can be measured along with depth profiles of major elements (Grime et al., 1993). In this way heterogeneities in the matrix both in the plane and normal to the surface plane may be accounted for and quantified. For comparison to the H-EBS data, α-RBS using heavier alpha particles can typically give improved elemental resolution and hence improve depth profiling for specific components.

Samples A (10−5 mol l−1), B (10−3 mol l−1), and C (10−1 mol l−1) were analysed. Sample C showed clear surface enrichment in Sr. Figure 5a shows a typical PIXE spectrum from this sample. For comparison, the Fig. 5b, show typical PIXE spectra taken from samples B and A, respectively. The major-element chemistry of all three analyses is essentially identical, as shown by the relative intensities of the K, Ca and Fe peaks. The Sr concentrations are completely different for each sample, however, as shown by the characteristic peak intensities at approximately 14 keV. Quantitative analysis of Sr PIXE data obtained from area analyses (15 × 15 μm) for each sample are presented on Table 2, and indicate that sample C has a Sr loading of 1.61 at. wt.%, whereas sample A has levels that are below the detection limit. Although the absolute errors increase from A to C, the relative errors decrease as the concentration of Sr rises well above the detection limits for both the PIXE and EBS measurements.

The PIXE concentration maps were also used to check for Sr heterogeneity. Most of the recorded characteristic Sr intensity in sample A can be ascribed to elevated Sr concentrations in a few small plagioclase feldspar grains included within the starting material. A map and point analyses of a feldspar grain and of the surrounding nearby clinoptilolite are shown in Fig. 6. These analyses clearly show that small quantities of Sr are indeed contained by the plagioclase (which also has higher Ca levels than the clinoptilolite) within the starting material. This is not a volumetrically important phase in the tuff and the Sr concentrations are low (110 ppm), but this nonetheless needs to be accounted for when considering the Sr inventory of these samples. We also note that the Sr signal is homogeneous and relatively high over the entire surface area of sample C, and hence we assume that the Sr intensity is overwhelmingly due to surface adsorption in this case and is unaffected by Sr in the starting material. Sample B shows some degree of compositional heterogeneity for Sr.

The EBS analyses resolved the presence of Sr on the surfaces of all three samples. A calculated fit to the sample C spectrum is presented in Fig. 7a and gives a surface concentration of 2.4 wt.% Sr, with the profile showing that Sr has penetrated to a depth of at least 4 μm below the sample surface. Using the bulk composition of the clinoptilolite reference material as a starting point, the α-RBS calculation not only gave results for Sr, but also showed surface contamination with carbon and refined the stoichiometry of the clinoptilolite to be Na3.4K2.2Ca0.1Fe0.4Al3.7Si25.8O72.5·24H2O. This best fit stoichiometry is close to the reference pattern and gives us confidence that the H-EBS fit is robust. Figure 7b,c also presents the best fit calculations of the H-EBS data for samples B and A, respectively. Surface Sr was detected in both cases by this technique, although in the case of sample A these concentrations were barely above the detection limit (Table 2). Perhaps most importantly, fits to the EBS and RBS data confirm that Sr has penetrated to at least 4 μm below the surface of the tuff in samples B and C.

Summary and conclusions

Table 3 presents a comparison of the concentrations of major and trace elements determined for the surface of the reacted sample C by the various analytical techniques employed in this study plus a reference composition for clinoptilolite taken from the literature. We note that all of the techniques employed were standardless. Given the heterogeneities in the starting material some variability is expected, however the accuracy and precision for all of the major elements is extremely good, except the ESEM-EDS analysis for Sr. The ion beam techniques agree within error for the surface loading of Sr in this experiment, and this, along with the De and Kd values determined above, therefore allows us to constrain surface loading and mass transfer for samples A, B and C.

A new technique, the micro-reactor simulated fracture method, was developed by adopting the concept of a micro chemical reactor for simultaneous determination of the diffusion coefficient and sorption coefficient of rock samples. The advantages of the new method were verified by comparing with conventional techniques. Strontium was chosen for the analyses and the breakthrough curves were obtained for a parallelepiped tuff sample. Good agreement was found between values obtained by the different methods within measurement errors. Furthermore, ion beam techniques were used to produce depth-resolved chemical analyses with low detection limits of reacted surfaces in order to fully constrain the inventory of Sr in the near surface of the reacted rock. It is clear that increased solution concentrations resulted in increased adsorbate density. The PIXE and EBS/RBS easily quantified Sr at the surfaces, however grains such as those imaged in Fig. 6 greatly complicated the analysis and created a heterogeneous background. Future work should be aware of such complications. In particular, for Sr diffusion measurements with natural materials investigators need to be aware of the high concentrations of Sr associated with some feldspar and carbonate minerals. The RBS depth profiling verified that Sr penetrated deeply below the surface of the channel (approximately 4 μm), however high sample roughness precluded accurately resolved depth profiling. Granite or another similar rock type that can be chemo-mechanically polished, with large crystals and well defined grain boundaries, will be used in future work.


Both RAW and TO gratefully acknowledge EPSRC grant EP/D032210/1 to the University of Surrey Ion Beam Centre and the assistance of Drs Geoff Grime and Karen Kirkby; PMM was supported via a University of Manchester Research Associate position granted with the assistance of Prof. R.A.D. Pattrick; BZ is supported by the EPSRC/Diamond consortium. Financial support was also given by the Central Research Institute of the Electrical Power Industry of Japan.

P2Freely available online through the publisher-supported open access option.