Quantifying carbon fixation in trace minerals from processed kimberlite: A comparative study of quantitative methods using X-ray powder diffraction data with applications to the Diavik Diamond Mine, Northwest Territories, Canada
Introduction
Emission of anthropogenic greenhouse gases (e.g., CO2, CH4, and N2O) has been implicated as a cause of current warming of the Earth’s climate. Carbon dioxide is by far the most significant of these greenhouse gases, representing 77% of total anthropogenic emissions in 2004 (IPCC, 2007). By 2005, the global atmospheric concentration of CO2 had increased to 379 ppm from a pre-1750 (i.e., pre-industrial) value of 280 ± 20 ppm. Approximately 2/3 of this increase is attributed to combustion of fossil fuels and ⅓ to changes in land-use patterns during the past 260 a (IPCC, 2007). It has been suggested that strategies for decarbonizing energy sources, increasing energy efficiency, and trapping and storing CO2 must be developed and implemented in order to stabilize concentrations of atmospheric CO2 and curtail the most damaging effects of anthropogenic climate change (e.g., Hoffert et al., 2002, Lackner, 2003, Pacala and Socolow, 2004, Broecker, 2007, IPCC, 2007).
Approximately 90% of C on Earth is fixed within carbonate minerals (Sundquist, 1985, Sundquist, 1993) and it is expected that these minerals will be the ultimate sink for most anthropogenic CO2 on a timescale of 1 Ma (Kump et al., 2000). Storage of CO2 in carbonate minerals is recognized as a safe and effective method for the sequestration of anthropogenic C (Seifritz, 1990, Lackner et al., 1995, Lackner, 2003). Precipitation of carbonate minerals in situ by dissolution of silicate mine residues is one potential implementation of this process. The development of secondary carbonate minerals has been documented in tailings at several mine sites in Canada: at the Kidd Creek Cu–Zn mine near Timmins, Ontario (Al et al., 2000), the Lower Williams Lake U mine near Elliot Lake, Ontario (Paktunc and Davé, 2002), and in chrysotile mine tailings at Thetford, Québec (Huot et al., 2003), Clinton Creek, Yukon Territory (Wilson et al., 2004), and Cassiar, British Columbia (Wilson et al., 2005). The C bound within secondary carbonate minerals in the tailings at some of these mines may not have had an atmospheric source. However, Wilson et al. (2009) demonstrate that tailings from the historical chrysotile mines at Clinton Creek, Yukon Territory and Cassiar, British Columbia are trapping and storing CO2 from the atmosphere. Accelerating the uptake of CO2 into tailings from active mines could reduce or offset the net greenhouse gas emissions of many mining operations.
Bulk geochemical methods for CO2 abundance cannot distinguish among various carbonate minerals nor can they discern the difference between atmospheric, bedrock, biological, and industrial sources of C within minerals. However, the sources of bound C can be distinguished using radiocarbon and stable isotopes of C and O (Wilson et al., 2009). Also, automated point-counting techniques (e.g., mineral liberation analysis) cannot be used to quantify fine-grained minerals or hydrous minerals that are easily vaporized by an electron beam. As an alternative to point counting, the amount of CO2 trapped within fine-grained, hydrous carbonate minerals can be estimated from weight-percent abundances determined with quantitative phase analysis using X-ray powder diffraction (XRPD) data.
At the Diavik Diamond Mine, Northwest Territories, Canada, efflorescent films of Ca-, Na- and Mg-carbonate minerals form in the tailings from the fine and coarse Processed Kimberlite Containment facilities (PKC). These minerals precipitate at the surface of kimberlite waste that is beached along the edge of a central pond used for storing process water. Based on the authors’ observations, the most common secondary carbonate mineral, and the one best preserved at depth, is nesquehonite (MgCO3·3H2O). Also, geochemical modelling by Rollo and Jamieson (2006) suggests that C mineralization may be occurring in waste kimberlite at the nearby EKATI Diamond Mine.
Processed kimberlite from Diavik contains a variety of minerals that are characterized by one or more of the following: (1) extensive solid solution, (2) structural disorder, and (3) severe preferred orientation. Furthermore, processed kimberlite at Diavik generally contains abundant serpentine and forsterite with minor to trace amounts of many other phases, resulting in complicated XRPD patterns that consist of many overlapped peak profiles. The combination of these factors presents a challenge for quantifying C mineralization with the Rietveld method (Rietveld, 1969). This is chiefly because the Rietveld method requires that the crystal structures and chemistry of the phases being analyzed be known and also gives less reliable results for minerals present at low abundances (e.g., Raudsepp et al., 1999).
Alternatively, Chung’s (1974) method of normalized reference intensity ratios (RIR) and similar methods have been used successfully to quantify trace abundances of minerals in multi-phase mixtures (e.g., Bish and Chipera, 1991, Omotoso et al., 2006). Calibration curves produced according to the internal standard method (Alexander and Klug, 1948) have also been used successfully to measure trace abundances of minerals (e.g., Sanchez and Gunter, 2006). These three methods are typically used in isolation on very different systems of minerals and, as such, it is difficult to recommend the use of one method over another for quantification of low abundances of minerals. Here the ability of Chung’s RIR method, the internal standard method, and the Rietveld method to measure minor to trace amounts of nesquehonite are assessed using weighed mixtures of pure mineral standards, prepared to simulate processed kimberlite. Based on the results of the tests, a procedure for accurate quantification of CO2 trapping within trace minerals in kimberlite mine tailings is outlined. This procedure is subsequently applied to samples of natural processed kimberlite from the fine and coarse PKC at Diavik. Results of quantitative phase analysis are used to estimate trapping of CO2 within nesquehonite at Diavik (from stoichiometry of this mineral) and to determine the contribution of nesquehonite to neutralization potential of the tailings according to the method of Jambor et al. (2007).
Section snippets
Locality and sampling strategy
The Diavik Diamond Mine is located on East Island, in Lac de Gras, approximately 300 km NE of Yellowknife, Northwest Territories, Canada (Fig. 1). There are four mineable kimberlite pipes on the Diavik property; two of which, A154 North and A154 South, are currently being mined from a single open pit. A third pipe, A418, is located to the south of the A154 pipes and began production in 2008. The fourth pipe, A21, is not currently being developed. The kimberlites at Lac de Gras, including those
Sample preparation and data collection
Thirteen mixtures of pure mineral samples were prepared to simulate processed kimberlite (mine tailings) (Table 1). Minerals used in these mixtures represent those commonly found in the tailings at Diavik: lizardite, high-Mg forsterite, diopside, almandine–pyrope series garnet, phlogopite, calcite, quartz, oligoclase, and nesquehonite. The abundance of nesquehonite in the synthetic tailings was varied from 0.10 to 5.00 wt.% in order to assess the applicability of the Rietveld, 1969, Chung, 1974
Rietveld refinement results
Rietveld refinements were initially done on data collected from serrated specimens of synthetic kimberlite tailings. The results for Series 1 and 2 mixtures gave significant overestimates for the abundance of lizardite while underestimating the abundances of forsterite and the minor phases. This effect bore a resemblance to the pattern of misestimates that would have resulted from a weighing error for the corundum spike (e.g., Gualtieri, 2000, Wilson et al., 2006). More precisely, it appeared
Implications for neutralization potential and carbon dioxide sequestration
Jambor et al. (2007) have demonstrated that the neutralization potential (NP) of a geological sample can be computed from Rietveld refinement results. In their study, they calculated the NP values of geological materials using appropriate values for NP of the individual minerals in each sample. These values were drawn from a library of NP data for pure mineral phases that had been determined using the Sobek method (Sobek et al., 1978). Most of the NP results computed by Jambor et al. (2007)
Conclusions
Trace amounts of secondary carbonate minerals in kimberlite mine tailings can be quantified accurately using either the Rietveld method or the method of normalized reference intensity ratios for abundances ⩾0.5%. A calibration curve, constructed using the internal standard method, provides more accurate results for abundances <0.5%.
Surface roughness effects, when combined with structureless pattern fitting and the use of an internal standard, can produce a pattern of systematic misestimates in
Acknowledgements
We would like to thank John Jambor for his tireless dedication to mineralogy and for the creativity and social responsibility he showed in his work. A more personal thanks to John for being a good friend to the Department of Earth and Ocean Sciences at UBC, for giving so generously of his knowledge and his time, and for sharing his enthusiasm with us.
We acknowledge the support of the Natural Sciences and Engineering Research Council of Canada (NSERC) and Diavik Diamond Mines Inc. (DDMI) through
References (89)
- et al.
Coupled thermal, hydraulic and geochemical evolution of pyritic tailings in unsaturated column experiments
Geochim. Cosmochim. Acta
(2007) - et al.
Carbonate-mineral/water interactions in sulfide-rich mine tailings
Geochim. Cosmochim. Acta
(2000) - et al.
Microstructure characterization of high energy ball-milled nanocrystalline V2O5 by Rietveld analysis
Mater. Chem. Phys.
(2003) - et al.
Carbon dioxide disposal in carbonate minerals
Energy
(1995) - et al.
C13 and O18 compositions in some fresh-water carbonates associated with ultramafic rocks and serpentinites: western United States
Geochim. Cosmochim. Acta
(1971) - et al.
The hydromagnesite playas of Atlin, British Columbia, Canada: a biogeochemical model for CO2 sequestration
Chem. Geol.
(2009) - et al.
Interaction of diamond mine waste and surface water in the Canadian Arctic
Appl. Geochem.
(2006) - et al.
Structure refinements of solid solutions in the system CaO–MgO–2SiO2–Al2O3
J. Korean Ceram. Soc.
(1986) - et al.
Basic aspects of X-ray absorption in quantitative diffraction analysis of powder mixtures
Anal. Chem.
(1948) - et al.
Variation of displacement parameters in structure refinements of low albite
Am. Mineral.
(1990)
Single-crystal X-ray structure study of synthetic pyrope–almandine garnets at 100 and 293 K
Am. Mineral.
Environmental geochemistry of kimberlite materials: Diavik Diamonds Project, Lac de Gras, Northwest Territories, Canada
Explor. Mining Geol.
Transmission electron microscopy of subsolidus oxidation and weathering of olivine
Contrib. Mineral. Petrol.
Detection of trace amounts of erionite using X-ray powder diffraction: erionite in tuffs of Yucca Mountain, Nevada, and Central Turkey
Clays Clay Miner.
Quantitative phase analysis using the Rietveld method
J. Appl. Crystallogr.
Quantitative mineralogical analysis using the Rietveld full-pattern fitting method
Am. Mineral.
The effect of grain or particle size on X-ray reflections from mixed powders and alloys, considered in relation to the quantitative determination of crystalline substances by X-ray methods
Philos. Mag.
CO2 arithmetic
Science
The electron distribution in corundum. A study of the utility of merging single-crystal and powder diffraction data
Acta Crystallogr., Sect. A: Found. Crystallogr.
A fundamental parameters approach to X-ray line-profile fitting
J. Appl. Crystallogr.
Multireflection RIR and intensity normalizations for quantitative analyses: applications to feldspars and zeolites
Powder Diffr.
Quantitative interpretation of X-ray diffraction patterns of mixtures. II. Adiabatic principle of X-ray diffraction analysis of mixtures
J. Appl. Crystallogr.
Computer simulation of structures and cohesive properties of micas
Am. Mineral.
Shear deformation experiments at 11 GPa–1400 °C in the multianvil apparatus
Eur. J. Mineral.
Accuracy in Rietveld quantitative phase analysis of Portland cements
J. Appl. Crystallogr.
Rietveld quantitative amorphous content analysis
J. Appl. Crystallogr.
Measurement of particle absorption by X-ray fluorescence
Acta Cryst.
Correction of intensities for preferred orientation in powder diffractometry: application of the March model
J. Appl. Crystallogr.
The crystal structure of nesquehonite, MgCO3·3H2O, from Lavrion, Greece
Mineral. Petrol.
Crystal structures of the low-temperature quartz-type phases of silica and germanium dioxide at elevated pressure
Z. Kristallogr.
The formation of weathering products on the LEW 85320 ordinary chondrite: evidence from carbon and oxygen stable isotope compositions and implications for carbonates in SNC meteorites
Meteoritics
Accuracy of XRPD QPA using the combined Rietveld – RIR method
J. Appl. Crystallogr.
Quantitative phase analysis from neutron powder diffraction data using the Rietveld method
J. Appl. Crystallogr.
Advanced technology paths to global climate stability: energy for a greenhouse planet
Science
Nesquehonite and dypingite from Shiraki, Mie Prefecture, Japan
Chigaku Kenkyu
Climate Change 2007: The Physical Basis
Static tests of neutralization potentials of silicate and aluminosilicate minerals
Environ. Geol.
Measured and computed neutralization potentials from static tests of diverse rock types
Environ. Geol.
Application of symmetrized harmonics expansion to correction of the preferred orientation effect
J. Appl. Crystallogr.
Cited by (71)
Experimental investigation of multiple industrial wastes for carbon dioxide removal strategies
2023, International Journal of Greenhouse Gas ControlA review of mineral carbonation by alkaline solidwaste
2022, International Journal of Greenhouse Gas ControlMineral carbonation of ultramafic tailings: A review of reaction mechanisms and kinetics, industry case studies, and modelling
2022, Cleaner Engineering and Technology