# Rapid CO2 mineralisation into calcite at the CarbFix storage site quantified using calcium isotopes

Philip a E Pogge von Strandmann, Kevin W Burton, Sandra O Snæbjörnsdóttir, Bergur Sigfússon, Edda S Aradóttir, Ingvi Gunnarsson, Helgi A Alfredsson, Kiflom G Mesfin, Eric H Oelkers, Sigurður R Gislason

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

## ► To cite this version:

Philip a E Pogge von Strandmann, Kevin W Burton, Sandra O Snæbjörnsdóttir, Bergur Sigfússon, Edda S Aradóttir, et al.. Rapid CO2 mineralisation into calcite at the CarbFix storage site quantified using calcium isotopes. Nature Communications, 2019, 10 (1), pp.1983. 10.1038/s41467-019-10003-8.
hal-04904381

HAL Id: hal-04904381

https://ut3-toulouseinp.hal.science/hal-04904381v1

Submitted on 21 Jan 2025

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

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

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

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

# Rapid CO₂ mineralisation into calcite at the CarbFix storage site quantified using calcium isotopes

Philip A.E. Pogge von Strandmann<fnref n="1" />, Kevin W. Burton<fnref n="2" />, Sandra O. Snæbjörnsdóttir³, Bergur Sigfússon⁴, Edda S. Aradóttir⁴, Ingvi Gunnarsson⁴, Helgi A. Alfredsson³, Kiflom G. Mesfin³,⁵, Eric H. Oelkers<fnref n="1" />,³,⁶ & Sigurður R. Gislason³

The engineered removal of atmospheric CO₂ is now considered a key component of mitigating climate warming below 1.5 °C. Mineral carbonation is a potential negative emissions technique that, in the case of Iceland's CarbFix experiment, precipitates dissolved CO₂ as carbonate minerals in basaltic groundwater settings. Here we use calcium (Ca) isotopes in both pre- and post-CO₂ injection waters to quantify the amount of carbonate precipitated, and hence CO₂ stored. Ca isotope ratios rapidly increase with the pH and calcite saturation state, indicating calcite precipitation. Calculations suggest that up to 93% of dissolved Ca is removed into calcite during certain phases of injection. In total, our results suggest that 165 ± 8.3 t CO₂ were precipitated into calcite, an overall carbon storage efficiency of 72 ± 5%. The success of this approach opens the potential for quantification of similar mineral carbonation efforts where drawdown rates cannot be estimated by other means.

<footnote n="1">¹LOGIC, Institute of Earth and Planetary Sciences, University College London and Birkbeck, University of London, Gower Street, London WC1E 6BT, UK.</footnote>

<footnote n="2">²Department of Earth Sciences, University of Durham, Durham DH1 3LE, UK. ³Institute of Earth Sciences, University of Iceland, Reykjavik 600169-2039, Iceland. ⁴Reykjavik Energy, Bæjarhálsi 1, Reykjavik 551298-3029, Iceland. ⁵HS Orka, Svartsengi, 240 Grindavik, Iceland. ⁶GET, CNRS UMR-5563, 14, Avenue Édouard Belin, Toulouse 31400, France. Correspondence and requests for materials should be addressed to P.A.E.P.v.S. (email: p.strandmann@ucl.ac.uk)</footnote>
The engineered removal of CO₂ from the atmosphere (negative emissions) is now thought to be crucial for keeping global warming under the limit of 1.5 °C¹. A plethora of different potential geo-engineering techniques has been proposed, with different costs, potential and environmental consequences of implementation². Indeed, the overall scale of necessary carbon capture and storage (CCS) suggests that a combination of different techniques would be most practical. An effective and, critically, efficient approach is the enhancement of natural carbon storage processes, in particular those associated with the long- and short-term carbon cycle, including carbonate formation, alkalinity drawdown and organic carbon enhancement³,⁴.

Of these, one particular method is mineral carbonation, which is the engineered enhancement of carbonate precipitation and hence CO₂ storage. The CarbFix programme in Iceland has now progressed from a pilot (original CarbFix)⁵ to the industrial scale (CarbFix2)⁶. In simple terms, CO₂ dissolved in water is injected into natural basaltic aquifers near the Hellisheiði geothermal power plant. Basalt, which at this site contains ~10 wt% CaO⁷ plus other carbonate-forming cations such as Mg or Fe, is highly reactive. The concept is that the precipitation of carbonates at depth occurs from carbon sourced from the CO₂ injection, and from calcium (or Mg or Fe) sourced from dissolution of the surrounding basalts. In the pilot experiment, carbonate precipitated within 2 years of injection at 20–50 °C and a few months at 60 to 260 °C⁵,⁶,⁸. This technique has the advantage that carbon is stored as solid carbonate minerals, a significantly more robust substrate than supercritical CO₂ stored via other CCS techniques, and that this mineralised carbon is relatively stable over long timescales⁵,⁹.

Calculation of the carbon stored during CO₂ injection at the pilot CarbFix site was made through the use of conservative (non-reactive) and reactive tracers, which were added to the injected CO₂-charged water (including radioactive ¹⁴C)⁵,⁸. Carbon storage within carbonate minerals was calculated by determining the difference between the amount of injected dissolved inorganic carbon (DIC), ¹⁴C and conservative tracers, with that sampled from a downstream monitoring well⁵,⁸. The stability of different minerals (different carbonates, as well as silicate secondary minerals) through the subsurface sequestration site has also been modelled, including prediction of groundwater conditions (temperature, pH, concentration) under which carbonate will form⁷,⁸,¹⁰. However, monitoring the extent and flux of carbonate and other secondary mineral precipitation using dissolved elemental concentrations alone is difficult because supply rates are not readily quantified, and, moreover, there are no diagnostic elemental signals that can uniquely indicate which metal-bearing phases have precipitated. Using concentration of solids from drill cores would also be extremely difficult, due to the relatively small amount of solid precipitated over millions of m³ and hundreds of Mt of basaltic rock. This highlights the difficulty of determining carbonation efficiency in field studies compared to laboratory experiments (e.g. ref. ¹¹).

Over recent years the stable isotopes of a number of elements have been used with considerable success as tracers of weathering processes in groundwaters and rivers. Calcium isotopes (here reported as δ⁴⁴/⁴⁰Ca relative to the standard NIST (National Institute of Standards and Technology) 915a) have an obvious potential as a tracer of carbonate precipitation processes. There is no known isotope fractionation associated with the dissolution of primary minerals, but there is significant fractionation associated with the precipitation of secondary Ca-rich phases¹². These are dominated by carbonates, but also include minerals such as Ca-smectites and -zeolites. In all cases precipitation leaves waters enriched in the heavy isotopes¹²⁻¹⁶. Therefore, calcium isotopes

have the potential to serve as a diagnostic tracer of carbonate precipitation rates (and hence CO₂ sequestration rates) at the CarbFix site.

Here we present Ca isotope data for the pre- and post-carbon injection groundwaters collected from the original CarbFix pilot site, to determine whether carbonate formation (and hence CO₂ sequestration) rates can be directly quantified. Ca isotope ratios increase with pH and calcite saturation states, indicating rapid calcite precipitation. The most rapid calcite mineralisation and CO₂ drawdown occurred within 1–2 months of the cessation of carbon injection. Overall, the Ca isotopes show that 165 ± 8.3 (2σ) t CO₂ of the 220 t injected were precipitated into calcite, confirming calcite as the primary carbon storage phase, and giving a carbon storage efficiency of 72 ± 5%.

## Results

**The CarbFix site.** The location of the original CarbFix pilot-study injection site is in SW Iceland, ~30 km east of Reykjavik, and 3 km SW of the Hellisheiði geothermal power plant, operated by Reykjavik Energy. The injection site has eight monitoring wells ranging from 50 to 2000 m depth (Supplementary Figs. 1–3)—six of these are located downstream from the HN-02 injection well. The shallow wells source water from 200 to 300 m depth, above a low-permeable hyaloclastite formation, while the deeper wells are encased to 400 m depth and sample the target formation (Supplementary Figs. 1–3). The injection site and its groundwaters have been characterised⁷, and a number of wells were sampled multiple times (up to three) over 3 years (2008–2010) prior to the CO₂ injection.

During the original CarbFix injection, water was pumped from well HN-01, and co-injected with CO₂ into well HN-02. This occurred in two phases: first, starting in late January 2012, 175 t of pure CO₂ was injected together with the water¹⁷. This phase ran continuously until 9 March 2012; the second phase was from mid-June to early August 2012, in which 73 t of a gas mixture containing 75 mol% CO₂, 24 mol% H₂S and 1 mol% H₂ was injected. The plume of injected material was monitored at the closest monitoring well HN-04, 500 m from HN-02 at depth. Changes in DIC, pH and tracer concentration were observed starting within 2 weeks of injection. Tracer injection suggests that ~95 ± 3% of the injected carbon was mineralised within 2 years⁵,⁸. The rapid conversion of the dissolved CO₂ to carbonate minerals is likely due to the novel method of dissolved CO₂ injection, the rapid dissolution rate of basalt, providing the necessary cations, such as Ca, Mg or Fe, the mixing of the injected water with alkaline groundwaters and the dissolution of pre-existing carbonates at the onset of CO₂ injection⁵. Analysis of solids recovered from the monitoring well and pumps shows the precipitation of calcite, but no aragonite precipitation⁵.

**Groundwater Ca isotope results.** We report Ca isotope ratios from 13 pre-injection aquifer samples characterising the injection site, including both shallow and deep groundwaters. Following this, we also report Ca isotope ratios from a post-injection time series of 19 samples from a single monitoring well. Calcium concentrations range from 1.2 to 21.7 µg/ml in the pre-injection samples⁷, and from 1.4 to 16.3 µg/ml in the post-injection samples⁸. All data are presented in Tables 1 and 2. Basalts have calcium isotope ratios of δ⁴⁴/⁴⁰Ca ~0.8‰¹⁴,¹⁵. Icelandic rivers have δ⁴⁴Ca values of 0.9–1.3‰, and previously measured groundwaters have values of 0.51–1.5‰¹⁴,¹⁵. In both cases, the CarbFix site groundwaters extend to heavier Ca isotope ratios. The shallow pre-injection groundwaters have relatively low δ⁴⁴Ca values, and also higher Ca/Na ratios than the deeper groundwaters (Fig. 1).
<table><caption>Table 1 Pre-injection conditions, elemental concentrations (from Alfredsson et al.⁷) and Ca isotope ratios</caption><thead><tr><th>Borehole</th><th>Sample</th><th>Sampling date</th><th>Depth</th><th>T (°C)</th><th>pH</th><th>Na (µg/ml)</th><th>Si (µg/ml)</th><th>Ca (µg/ml)</th><th>Sr (ng/ml)</th><th>Li (ng/ml)</th><th>Calcite SI</th><th>δ<sup>44</sup>/<sup>40</sup>Ca</th><th>2 s.e.</th></tr></thead><tbody><tr><td>HK-7b</td><td>08HAA08</td><td>08/07/2008</td><td>Shallow</td><td>12.4</td><td>7.65</td><td>9.17</td><td>12.0</td><td>11.3</td><td>20.1</td><td></td><td>-1.13</td><td></td><td>0.05</td></tr><tr><td>HN-4</td><td>08HAA01</td><td>01/07/2008</td><td>Deep</td><td>32.3</td><td>9.43</td><td>48.6</td><td>25.1</td><td>1.64</td><td>1.62</td><td>0.40</td><td>0.12</td><td>1.95</td><td>0.05</td></tr><tr><td>HK-34</td><td>08HAA03</td><td>04/07/2008</td><td>Deep</td><td>25</td><td>9.63</td><td>49.6</td><td>22.7</td><td>2.17</td><td>3.98</td><td></td><td></td><td></td><td>0.04</td></tr><tr><td>HK-31</td><td>08HAA05</td><td>04/07/2008</td><td>Deep</td><td>17.4</td><td>9.29</td><td>41.8</td><td>20.7</td><td>3.75</td><td>11.9</td><td>0.31</td><td>0.29</td><td>2.04</td><td>0.04</td></tr><tr><td>HK-26</td><td>08HAA06</td><td>04/07/2008</td><td>Deep</td><td>18.8</td><td>8.44</td><td>72.5</td><td>17.8</td><td>5.72</td><td>28.3</td><td></td><td></td><td></td><td></td></tr><tr><td>HK-26</td><td>09HAA17</td><td>29/05/2009</td><td>Deep</td><td>16.5</td><td>8.51</td><td>77.0</td><td>16.7</td><td>5.36</td><td>26.5</td><td>0.59</td><td>0.02</td><td>1.22</td><td>0.05</td></tr><tr><td>HK-12</td><td>09HAA18</td><td>29/05/2009</td><td>Shallow</td><td>5.3</td><td>8.33</td><td>8.82</td><td>8.63</td><td>8.39</td><td>14.8</td><td>0.16</td><td>-0.78</td><td>1.11</td><td>0.05</td></tr><tr><td>HK-31</td><td>09HAA19</td><td>29/05/2009</td><td>Deep</td><td>18.9</td><td>9.41</td><td>48.6</td><td>21.1</td><td>3.69</td><td>11.9</td><td>0.44</td><td>0.37</td><td>1.91</td><td>0.05</td></tr><tr><td>HK-25</td><td>09HAA20</td><td>29/05/2009</td><td>Shallow</td><td>7.4</td><td>8.09</td><td>8.00</td><td>9.0</td><td>8.02</td><td>13.9</td><td>0.28</td><td>-1.00</td><td>1.09</td><td>0.04</td></tr><tr><td>HK-34</td><td>09HAA21</td><td>29/05/2009</td><td>Deep</td><td>27.5</td><td>9.79</td><td>55.0</td><td>24.6</td><td>1.30</td><td>1.74</td><td>0.23</td><td>0.15</td><td>2.04</td><td>0.03</td></tr><tr><td>rpt</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>1.98</td><td>0.06</td></tr><tr><td>HK-7b</td><td>09HAA22</td><td>29/05/2009</td><td>Shallow</td><td>11.7</td><td>7.67</td><td>12.2</td><td>13.0</td><td>21.7</td><td>38.9</td><td>0.31</td><td>-0.50</td><td>1.14</td><td>0.03</td></tr><tr><td>HN-1</td><td>09HAA23</td><td>29/05/2009</td><td>Deep</td><td>24.7</td><td>9.26</td><td>43.5</td><td>14.1</td><td>5.37</td><td>18.1</td><td>0.16</td><td>0.51</td><td>2.07</td><td>0.04</td></tr><tr><td>HN-4</td><td>09HAA24</td><td>05/06/2009</td><td>Deep</td><td>34.5</td><td>9.56</td><td>55.6</td><td>25.2</td><td>1.32</td><td>1.27</td><td>0.27</td><td>0.10</td><td>1.89</td><td>0.04</td></tr><tr><td>HK-26</td><td>10HAA26</td><td>25/06/2010</td><td>Deep</td><td>17.2</td><td>8.65</td><td>74.5</td><td>16.6</td><td>5.25</td><td>26.1</td><td>0.32</td><td>0.05</td><td>1.62</td><td>0.04</td></tr><tr><td>HK-31</td><td>10HAA28</td><td>25/06/2010</td><td>Deep</td><td>16.9</td><td>9.55</td><td>48.6</td><td>23.1</td><td>3.61</td><td>12.0</td><td>0.41</td><td>0.42</td><td>1.91</td><td>0.06</td></tr><tr><td>HN-4</td><td>10HAA29</td><td>25/06/2010</td><td>Deep</td><td>30.4</td><td>9.69</td><td>54.7</td><td>30.6</td><td>1.29</td><td>1.14</td><td></td><td></td><td></td><td></td></tr><tr><td>HK-34</td><td>10HAA30</td><td>25/06/2010</td><td>Deep</td><td>28.0</td><td>9.88</td><td>51.0</td><td>27.1</td><td>1.21</td><td>1.67</td><td>0.24</td><td>0.03</td><td>1.75</td><td>0.03</td></tr></tbody></table>

Borehole names are the same as in Figs. S1-3. The shallow wells, HK-7b, HK-12, and HK-25, draw water from the dominate aquifers in the upper system, whereas the deeper wells, HN-2, HN-4, HK-34, HN-1, HK-31, and HK-26, draw water from the highest discharge aquifers below 400 m since these wells are cased down to that depth⁷

<table><caption>Table 2 Post-injection conditions, elemental concentrations (from Snæbjörnsdóttir et al.⁸) and Ca isotope ratios</caption><thead><tr><th>Sample</th><th>Sampling date</th><th>pH</th><th>Na (mg/ml)</th><th>Ca (µg/ml)</th><th>Sr (ng/ml)</th><th>Li (ng/ml)</th><th>Calcite SI</th><th>δ<sup>44</sup>/<sup>40</sup>Ca</th><th>2 s.e.</th></tr></thead><tbody><tr><td>12KGM08</td><td>09/02/2012</td><td>8.98</td><td>50.8</td><td>5.44</td><td>3.38</td><td>0.28</td><td>0.44</td><td>1.54</td><td>0.03</td></tr><tr><td>12KGM11</td><td>16/02/2012</td><td>7.94</td><td>52.0</td><td>7.81</td><td>4.67</td><td>0.29</td><td>-0.28</td><td>1.12</td><td>0.05</td></tr><tr><td>12KGM19</td><td>27/02/2012</td><td>7.18</td><td>54.3</td><td>9.56</td><td>5.62</td><td>0.31</td><td>-0.88</td><td>1.05</td><td>0.04</td></tr><tr><td>12KGM25</td><td>08/03/2012</td><td>6.79</td><td>53.2</td><td>14.4</td><td>8.58</td><td>0.35</td><td>-1.07</td><td>0.89</td><td>0.05</td></tr><tr><td>12KGM33</td><td>26/03/2012</td><td>6.71</td><td>54.6</td><td>16.3</td><td>10.0</td><td>0.36</td><td>-1.08</td><td>1.09</td><td>0.04</td></tr><tr><td>12KGM44</td><td>18/04/2012</td><td>7.70</td><td>54.6</td><td>12.0</td><td>7.62</td><td>0.35</td><td>-0.28</td><td>1.29</td><td>0.04</td></tr><tr><td>12KGM49</td><td>04/05/2012</td><td>9.00</td><td>55.8</td><td>6.50</td><td>4.23</td><td>0.32</td><td>0.57</td><td>2.13</td><td>0.06</td></tr><tr><td>12KGM60</td><td>30/05/2012</td><td>8.81</td><td>53.7</td><td>8.75</td><td>5.87</td><td>0.32</td><td>0.57</td><td>2.11</td><td>0.04</td></tr><tr><td>12SOS01</td><td>28/06/2012</td><td>7.36</td><td>55.0</td><td>11.9</td><td>8.14</td><td>0.36</td><td>-0.60</td><td>1.07</td><td>0.04</td></tr><tr><td>12SOS09</td><td>17/07/2012</td><td>8.28</td><td>55.8</td><td>10.9</td><td>7.86</td><td>0.36</td><td>0.25</td><td>1.91</td><td>0.05</td></tr><tr><td>rpt</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>1.94</td><td>0.04</td></tr><tr><td>12SOS15</td><td>31/07/2012</td><td>8.32</td><td>55.9</td><td>11.3</td><td>8.13</td><td>0.35</td><td>0.37</td><td>2.13</td><td>0.05</td></tr><tr><td>rpt</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>1.99</td><td>0.06</td></tr><tr><td>12SOS21</td><td>14/08/2012</td><td>7.25</td><td>57.0</td><td>14.3</td><td>9.96</td><td>0.43</td><td>-0.79</td><td>0.80</td><td>0.05</td></tr><tr><td>12SOS28</td><td>28/08/2012</td><td>7.50</td><td>57.9</td><td>13.1</td><td>9.57</td><td>0.44</td><td>0.02</td><td>1.64</td><td>0.06</td></tr><tr><td>12SOS34</td><td>24/09/2012</td><td>8.23</td><td>60.5</td><td>12.1</td><td>9.53</td><td>0.36</td><td>0.29</td><td>2.18</td><td>0.02</td></tr><tr><td>12SOS39</td><td>29/10/2012</td><td>8.26</td><td>61.6</td><td>11.0</td><td>9.45</td><td>0.46</td><td>0.22</td><td>2.18</td><td>0.04</td></tr><tr><td>13SOS01</td><td>07/01/2013</td><td>8.73</td><td>62.7</td><td>6.85</td><td>7.42</td><td>0.42</td><td>0.41</td><td>2.02</td><td>0.04</td></tr><tr><td>13SOS10</td><td>16/04/2013</td><td>8.76</td><td>62.4</td><td>7.19</td><td>7.38</td><td>0.44</td><td>0.47</td><td>2.10</td><td>0.03</td></tr><tr><td>13SOS17</td><td>10/06/2013</td><td>8.86</td><td>60.4</td><td>7.13</td><td>7.43</td><td>0.56</td><td>0.54</td><td>1.54</td><td>0.03</td></tr><tr><td>14SOS11</td><td>17/03/2014</td><td>9.08</td><td>63.7</td><td>3.49</td><td>4.12</td><td>0.56</td><td>0.31</td><td>1.39</td><td>0.05</td></tr></tbody></table>

These samples are a time series from the monitoring well HN-4 with a constant temperature of 35 °C

Pre-injection samples have δ<sup>44</sup>/<sup>40</sup>Ca values of 1.09 to 2.07‰ (relative to SRM-915a), and exhibit a positive relationship with pH (Fig. 2). The post-injection δ<sup>44</sup>/<sup>40</sup>Ca values exhibit similar overall isotope ratios of 0.8–2.1‰, but have considerably more scatter than the samples from the monitoring well before the CO₂ injection (1.89–1.95‰ in HN-04). The post-injection samples also show a positive co-variation with pH (Fig. 2).

Discussion

Prior to the CO₂ injection, the Ca isotopic compositions of the aquifers show a significant variation (1.09 to 2.07‰ relative to SRM-915a), trending from similar to somewhat higher values than observed in Icelandic rivers¹⁴,¹⁵ (Fig. 1). These studies

suggest that rivers' δ<sup>44</sup>Ca are largely controlled by mixing between the dissolution of basalt (δ<sup>44</sup>Ca ~ 0.8‰) and calcite (which is undersaturated and hence dissolving) plus hydrothermal and meltwater inputs. Additional isotope fractionation is caused by the formation of Ca-bearing secondary minerals, such as heulandite or stilbite.

Fractionation processes can be examined by comparing Ca isotope ratios to Ca/Na ratios. In basaltic environments, Na is the most mobile of the major cations¹⁸, consequently it most readily goes into solution compared to less mobile cations, which are preferentially taken up by secondary minerals. Only the shallow groundwaters have generally similar Ca/Na values to Icelandic rivers (Fig. 1). In contrast, the deeper aquifers tend to have lower pre-injection Ca/Na, suggesting that relatively more Ca is being
removed from these waters into secondary phases than in most rivers. In these groundwaters there is no correlation between δ<sup>44</sup>Ca and the calculated saturation index of zeolites such as heulandite<sup>8,10,19</sup>. Equally, there is no co-variation between δ<sup>44</sup>Ca and Mg isotope ratios, where the latter are thought to be controlled by the precipitation of smectites<sup>20</sup>. This lack of correlation also suggests that any precipitating secondary silicates are Mg-rich, but Ca-poor. However, an additional process that is highly likely to be controlling Ca isotopes in these waters is the precipitation of carbonate, especially in the higher pH aquifers. Calcite has been observed precipitating as a secondary mineral from these pre-injection waters<sup>7,10</sup>, as well as in other basaltic settings<sup>21</sup>, and indeed the precipitation of carbonate is the main goal of the CarbFix project.

There is a strong relationship between calcite saturation (calculated using PHREEQC<sup>8</sup>) and δ<sup>44</sup>Ca (Fig. 2b), largely controlled by a similar co-variation between δ<sup>44</sup>Ca and pH (Fig. 2a). Importantly, when calcite is undersaturated, δ<sup>44</sup>Ca is low (close to basaltic values, and similar to Icelandic rivers<sup>14,15</sup>), which likely represents secondary Ca-silicate precipitation. In contrast, when calcite is

supersaturated, and therefore likely precipitating, δ<sup>44</sup>Ca in the waters rapidly increases. Calcite is known to preferentially take up light Ca isotopes. The precise value for the fractionation factor is debated, but is around 0.6–0.8‰ (α = 0.999–0.9995) for inorganic calcite, with a degree of temperature-dependence<sup>13,22</sup>. Using the relationship from Gussone et al.<sup>7</sup> and the measured temperatures, the pre-injection samples are likely to have fractionation factors into calcite of ~0.9991–0.9995, due to temperatures ranging from 5 to 35 °C. Using a standard Rayleigh fractionation equation and assuming that all Ca initially originates from isotopically conservative basalt dissolution (this method is detailed in the Methods, and has previously been used to calculate similar parameters in speleothem systems<sup>16</sup>; equilibrium fractionation using this assumption cannot reproduce the higher δ<sup>44</sup>Ca values), it is calculated that around 28% of Ca is taken up into calcite from solution at pH < 8.5, but this reaches ~90% at pH 9.5–10. In other words, at high pH, most of the dissolved Ca is taken up by calcite precipitation, following these calculations.

The CarbFix injection has a clear effect on Ca isotopes: the pre-injection waters from well HN-04 (the primary monitoring well) had δ<sup>44</sup>Ca values of 1.89–1.95‰ between 2008 and 2009. In comparison, the same well exhibits δ<sup>44</sup>Ca values of 0.8–2.1‰ post-injection. The effect of the CO<sub>2</sub> addition in both injection phases was initially to decrease pH by over 2 units<sup>8</sup>. Following this, pH recovered (Fig. 3), because the dissolution of the host basalts and fluid mixing neutralised the pH of the injected fluids<sup>8</sup>, allowing the precipitation of calcite<sup>5,8</sup> (but no observed aragonite<sup>5</sup>). Hence, there are two periods of calcite undersaturation, matching the two phases of injection. For the rest of the post-injection period, calcite is supersaturated, and is demonstrably precipitating on and within the pumps and pipes in the monitoring well, and is the only recovered phase containing the injected radioactive carbon (<sup>14</sup>C)<sup>5,8</sup>.

Similar to the pre-injection samples, the post-injection δ<sup>44</sup>Ca values correlate with calcite saturation (Fig. 2). Due to an increased groundwater flow rate during injection (see Methods), water and elemental fluxes are significantly higher in post-compared to pre-injection waters. Assuming that the Ca isotope ratios in the post-injection samples are dominantly controlled by calcite precipitation, identical calculations can be performed as for the pre-injection aquifers. Due to the high reactivity and high Ca concentrations of basalt, the vast majority of the calcium (and other cations necessary for carbonate precipitation, such as Mg or Fe) are supplied from the dissolution of basalt. As such, the Ca dissolved in pre-injection groundwaters is insignificant for the purposes of enhanced carbonate precipitation and the following

<figure><img src="image_4.png" /><figcaption>**Fig. 1** Ca isotopes as a function of Ca/Na ratios. Elemental/Na ratios theoretically represent the uptake of the element in question by secondary minerals, relative to the highly mobile Na cation. The red diamonds with grey outlines represent the shallow groundwaters, while the pure red diamonds represent the injection well. Black crosses are previously measured Icelandic rivers and groundwaters<sup>14,15</sup>. Pre-injection samples are from a range of monitoring wells, while post-injection samples are a time series from a single well (HN-4)</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>**Fig. 2** Relationship between Ca isotopes and groundwater pH and calcite saturation indices. **a** Ca isotope ratios plotted against water pH; **b** and against the calcite saturation index. Both plots show that Ca isotope ratios increase as the conditions for calcite precipitation improve. Pre-injection waters are the red diamonds, and post-injection waters the blue squares. The horizontal dashed black line represents the Ca isotope ratio of basalt, while the vertical dashed black line represents saturation (SI = 0). The error bars represent 2 s.e. internal analytical uncertainty</figcaption></figure>
<figure><img src="image_6.png" /><figcaption>**Fig. 3** The trends of the calcite saturation state and Ca isotopes with time during the carbon injection phase from monitoring borehole HN-4. **a** Ca isotope ratios (blue line) and calcite saturation index (dotted black line). The blue shaded area represents the 2 s.e. analytical uncertainty on the isotope measurements. **b** shows the evolution of pH (dotted black line) and the calculated CO₂ precipitation rate based on Ca isotopes (red line). The red shaded area represents the 2σ propagated uncertainty on the individual precipitation rates. The grey shaded areas represent the carbon injection periods</figcaption></figure>

mass balance calculations⁵,⁹. In this case, because the samples all stem from the same borehole, the temperature was constant at 35 °C⁸. Assuming that only calcite precipitation fractionated Ca isotopes towards higher values in this setting, and taking a temperature-dependent isotope fractionation factor¹³,²² of 1000 $ln \alpha = -1.02 \pm 0.25 + (0.015 \pm 0.013) \times T$ (°C), this gives a fractionation factor of $\alpha = 0.9995$. The large pH variations due to the injection of CO₂ lead to 0 to 93% of Ca taken up into precipitated calcite. Hence, at periods of high pH, Ca storage behaviour is similar to that reported to occur to CO₂, in that almost 95% is sequestered into carbonates⁵,⁸.

The precipitation rates and total mass of calcite formed during and after the CarbFix CO₂ injection can be determined from the initial isotope ratio of the basalt, the isotope ratios and Ca concentrations of the measured fluid samples, the fluid-calcite fractionation factor and the water flow rate. Details of this calculation are provided in the Methods. A significant advantage of determining these rates from stable metal isotope measurements in such systems is that it potentially provides information on the precipitation rates of distinct minerals simply by measuring the isotope ratios of the groundwaters and the flow rate, and without the necessity of adding tracers to the injected carbon.

In the original CarbFix experiment, approximately 230 t CO₂ was injected in two stages. Matter et al.⁵ estimated a CO₂ sequestration efficiency of around 95 ± 3% based on ¹⁴C and 98 ± 4% based on DIC. Hence, ~220 t of CO₂ were precipitated during the initial 2 years, with modelling suggesting this dominantly was in calcite (with modelling also suggesting that trace Fe-carbonate might form at the lowest pH¹⁹). To corroborate this, only calcite was found to contain the ¹⁴C tracer⁵,⁸. The Ca isotopes suggest that approximately 165 ± 8.3 (2σ) t CO₂ were precipitated into calcite, a sequestration efficiency of 72 ± 5%, confirming calcite as the dominant carbon mineralisation phase (Fig. 3—see Methods for error propagation formulae). The greatest precipitation rates occur within 2 months following the cessation of injection. The

estimate based on Ca isotopes is thus lower than that based on injected carbon tracers and a conservative mass balance approach⁵. In effect, this suggests that more carbon is being removed from solution than is being mineralised by calcite. This may be due to the early formation of Fe-Mg carbonate minerals over calcite at relatively low pH¹⁹. Overall, these results show that Ca isotopes can be used to determine the amount of CO₂ stored in calcite in the subsurface without the use of additional tracers.

An interesting corollary to the Ca isotope method is its ability to time-resolve calcite precipitation rates (as opposed to carbon removal rates from the injected tracers), which shows both that precipitation rates are highest 1–2 months after cessation of CO₂ injection, and that $\delta^{44}\text{Ca}$ values decrease towards the end of the experiment, starting about a year after the final injection phase (Fig. 3). In both cases, the observations suggest that the calcite precipitation rates are not always directly coupled to pH or calculated saturation indices (Fig. 3), likely because precipitation is driven by both fluid saturation and reactive surface area, which can vary subtly along the flowpath and with time. For example, the slower rate of precipitation near the end of the original Carbfix experiment (down to ~2 mol/h CO₂ compared to a height of >30 mol/h) is despite pH and calcite saturation remaining relatively high. It seems likely that this reflects the slowing of the reaction as the injected carbon (i.e. the reactants) becomes exhausted, especially given the efficiency of carbon removal from solution⁸, and the system returning to its previous state (i.e. supersaturated for calcite, but with low precipitation due to low flow rates and fluxes⁷). Overall, therefore, a tracer of precipitation rates that can act in a time-resolved manner (such as Ca isotopes) can be very useful in understanding and quantifying subsurface carbon storage.

In summary, the original CarbFix industrial pilot experiment involved pumping dissolved CO₂ into basaltic groundwaters to precipitate carbonate, a carbon storage method known as mineral carbonation. Initial estimates, based on conservative tracers of carbon, suggested that 95 ± 3% of the added CO₂ (~220 t) was
removed through mineralisation. Using calcium isotopes, we estimate the amount of Ca lost from the groundwaters due to carbonate formation. The data here suggest that 165 ± 8.3 (2σ) t CO₂ were precipitated into calcite (confirming calcite as the primary carbon storage phase), a storage efficiency of 72 ± 5% (2σ). These two estimates are therefore just outside of uncertainty of each other, confirming modelling that suggests that carbon initially was mineralised into Mg-Fe carbonates at low pH. Overall, the study indicates that the efficient storage of carbon in the CarbFix experiment has been confirmed through independent methods. This method also allows a time-resolved determination of calcite precipitation rates that is distinct from pH or saturation calculations. Further, Ca isotopes appear to be a useful tool for quantifying mineral carbonation in settings where tracer injection cannot take place, when permissions (or finances) for using chemical or radioactive do not exist, or during natural weathering involving rapid carbonation.

## Methods

**Analytical methods.** Calcite SI (saturation indices) are taken from Alfredsson et al.⁷ for pre-injection samples, and from Snæbjörnsdóttir et al.⁸ for post-injection samples. In both cases they were calculated using the PHREEQC programme, using measured pH, temperature, DIC and ionic strengths of individual samples.

Around 10 µg of calcium was purified through a two-stage column procedure, the first column containing AG50 X12 resin, which removes most matrix elements, and the second containing Sr-spec resin, which removes any Sr. Splits collected before and after the Ca elution peak contained <0.5% of the Ca, indicating column yields of >99.5%. Analyses were performed on a Nu Instruments MC-ICP-MS at Oxford, relative to the standard SRM-915a. Sr isobaric interference was monitored at mass 43.5, and used to correct the ⁴²Ca, ⁴³Ca, and ⁴⁴Ca intensities. This column and mass spectrometric methodology has been detailed in a series of studies²³⁻²⁶. Seawater measured by this method yielded δ⁴⁴/⁴⁰Ca values of 1.92 ± 0.14‰ (2 s.d., n = 16), in keeping with other studies¹⁴,²⁷.

**Precipitation calculations.** Given an initial isotope ratio (basalt dissolution provides almost all of the dissolved Ca), the isotope ratio of the measured solutions, and a fractionation factor, the fraction of an element remaining in solution relative to that taken into the solid can be calculated using a standard Rayleigh fractionation equation (Eq. 1):

$$ \delta = \delta^i + 1000(\alpha - 1)\ln f, \quad (1) $$

where δ is the δ⁴⁴/⁴⁰Ca of the groundwater and δ is the δ⁴⁴/⁴⁰Ca of the starting compositions of the input Ca, in this case the Ca isotope ratio of basalt. α is the isotopic fractionation factor, and f is the fraction of Ca remaining in solution¹⁶. This equation is solved for f, which allows calculation of the amount of Ca precipitated, when using the Ca concentration of the solution. This approach also has the advantage that any dissolution of pre-existing calcite, which may occur immediately after injection when pH is low (Fig. 3), will result in very high values of f. Hence, the values of f = 0.84–1.00 determined for these periods means that very little apparent calcite precipitation occurred (Supplementary Fig. 4). In fact, if these values are set to >1 (implying Ca gain from dissolution), the impact on the final amount of sequestered CO₂ is <0.1 t.

In turn, f values can be converted to a calcite (and hence CO₂) precipitation rate (e.g. mol/h), by factoring in the water flow rate, according to the standard equation²⁸ (Eq. 2):

$$ [x] = \left[ \frac{(x)_{\text{solution}}}{f} - (x)_{\text{solution}} \right] D, \quad (2) $$

where, in this case, [x] is the Ca concentration in carbonate, (x)ₐ₉₉₉₇ is the concentration in the groundwater, and D is the x/Ca partition coefficient (here =1). The flux was obtained by multiplying the concentration (in mol/kg) by the flow rate (in l/h). Carbonate and CO₂ precipitation were calculated according to (using here a simplified equation based on wollastonite):

$$ \mathrm{CO}_{2}+\mathrm{H}_{2} \mathrm{O}+\mathrm{CaSiO}_{3}=\mathrm{CaCO}_{3}+\mathrm{H}_{2} \mathrm{O}+\mathrm{SiO}_{2} $$

This is slightly different from the normal silicate weathering equation because the pH is considerably lower in the injected medium, and hence the ocean-based carbonate cycle does not apply.

The natural groundwater flow rate was extremely slow (~25 m/yr), but, during the CarbFix experiment, water was pumped into and out of different boreholes to speed the water flow rate (input pump rate 7200 l/h; output pump rate 3500 l/h). This then allows us to calculate mineral precipitation rates. In principle, this could be a significant advantage of measuring stable metal isotopes in such systems, because it potentially provides information on the precipitation rates of different minerals simply by measuring the isotope ratios of the groundwaters and the flow rate. In all cases, the isotope ratio of the injected water was not taken into account

(well HN-1). This is because the pre-injection Ca fluxes are insignificant compared to post-injection ones, due to the increased flow rates. Ca/Na ratios also approximately an order of magnitude lower in the pre-injection waters (Fig. 1). Also, the Ca isotope ratio is significantly different from the post-injection values.

Uncertainty was propagated from a combination of analytical uncertainty (isotopic and concentration) and the uncertainty on the isotopic fractionation factor (also effectively analytical uncertainty). The propagated uncertainty was then calculated according to (Eq. 3):

$$ \sigma_P^2 = (t_2 - t_1)^2 \sigma_1^2 / 4 + (t_n - t_{n-1})^2 \sigma_n^2 / 4 + \sum_{i=2}^{n-1} (t_{i+1} - t_{i-1})^2 \sigma_i^2 / 4, \quad (3) $$

where t is the time at the ith data point. Of the reported uncertainty (Fig. 3), approximately 40% is from analyses (similar to other Ca isotope studies that also calculated values of f¹⁶), and the remainder largely from uncertainty on the fractionation factor.

## Data availability

The authors declare that the data supporting the findings of this study are available within the paper (and its supplementary information files).

Received: 12 November 2018 Accepted: 8 April 2019

Published online: 30 April 2019

## References

1. Masson-Delmotte, V. et al. An IPCC Special Report on the impacts of global warming of 1.5°C above pre-industrial levels and related global greenhouse gas emission pathways, in the context of strengthening the global response to the threat of climate change, sustainable development, and efforts to eradicate poverty. World Meteorological Organisation, Geneva (2018).

2. Shepherd, J.G. et al., Working Group on Geoengineering the Climate Geoengineering the climate: science, governance and uncertainty (RS Policy document, 10/29) London, GB. Royal Society 98pp (2009).

3. Hartmann, J. et al. Enhanced chemical weathering as a geoengineering strategy to reduce atmospheric carbon dioxide, supply nutrients, and mitigate ocean acidification. Rev. Geophys. 51, 113–150 (2013).

4. Renforth, P., Pogge von Strandmann, P. A. E. & Henderson, G. M. The dissolution of olivine added to soil: implications for enhanced weathering. Appl. Geochem. 61, 109–118 (2015).

5. Matter, J. M. et al. Rapid carbon mineralization for permanent disposal of anthropogenic carbon dioxide emissions. Science 352, 1312–1314 (2016).

6. Gunnarsson, I. et al. The rapid and cost-effective capture and subsurface mineral storage of carbon and sulfur at the CarbFix2 site. Int. J. Greenh. Gas Control 79, 117–126 (2018).

7. Alfredsson, H. A. et al. The geology and water chemistry of the Hellisheiði, SW-Iceland carbon storage site. Int. J. Greenh. Gas Control 12, 399–418 (2013).

8. Snæbjörnsdóttir, S. O. et al. The chemistry and saturation states of subsurface fluids during the in situ mineralisation of CO₂ and H₂S at the CarbFix site in SW-Iceland. Int. J. Greenh. Gas Control 58, 87–102 (2017).

9. Gislason, S. R. & Oelkers, E. H. Carbon storage in basalt. Science 344, 373–374 (2014).

10. Aradottir, E. S. P., Sonnenthal, E. L., Bjornsson, G. & Jonsson, H. Multidimensional reactive transport modeling of CO₂ mineral sequestration in basalts at the Hellisheiði geothermal field, Iceland. Int. J. Greenh. Gas Control 9, 24–40 (2012).

11. Wang, F., Dreisinger, D., Jarvis, M., Hitchins, T. & Dyson, D. Quantifying kinetics of mineralization of carbon dioxide by olivine under moderate conditions. Chem. Eng. J. 360, 452–463 (2019).

12. Fantle, M. S. & Tipper, E. T. Calcium isotopes in the global biogeochemical Ca cycle: implications for development of a Ca isotope proxy. Earth-Sci. Rev. 129, 148–177 (2014).

13. Marriott, C. S., Henderson, G. M., Belshaw, N. S. & Tudhope, A. W. Temperature dependence of δ⁷Li, δ⁴⁴Ca and Li/Ca during growth of calcium carbonate. Earth Planet. Sci. Lett. 222, 615–624 (2004).

14. Hindshaw, R. S., Bourdon, B., Pogge von Strandmann, P. A. E., Vigier, N. & Burton, K. W. The stable calcium isotopic composition of rivers draining basaltic catchments in Iceland. Earth Planet. Sci. Lett. 374, 173–184 (2013).

15. Jacobson, A. D., Andrews, M. G., Lehn, G. O. & Holmden, C. Silicate versus carbonate weathering in Iceland: new insights from Ca isotopes. Earth Planet. Sci. Lett. 416, 132–142 (2015).

16. Owen, R. A. et al. Calcium isotopes in caves as a proxy for aridity: modern calibration and application to the 8.2 kyr event. Earth Planet. Sci. Lett. 443, 129–138 (2016).
17. Sigfusson, B. et al. Solving the carbon-dioxide buoyancy challenge: the design and field testing of a dissolved CO₂ injection system. *Int. J. Greenh. Gas Control* **37**, 213–219 (2015).

18. Gislason, S. R., Arnorsson, S. & Armannsson, H. Chemical weathering of basalt in southwest Iceland: Effects of runoff, age of rocks and vegetative/ glacial cover. *Am. J. Sci.* **296**, 837–907 (1996).

19. Snæbjörnsdóttir, S. O., Gislason, S. R., Galeczka, I. M. & Oelkers, E. H. Reaction path modelling of in-situ mineralisation of CO₂ at the CarbFix site at Hellisheidi, SW-Iceland. *Geochim. Cosmochim. Acta* **220**, 348–366 (2018).

20. Oelkers, E. H. et al. Using stable Mg isotope signatures to assess the fate of magnesium during the in situ mineralisation of CO₂ and H₂S at the CarbFix site in SW-Iceland. *Geochim. Cosmochim. Acta* **245**, 542–555 (2019).

21. Galeczka, I. et al. The effect of the 2002 glacial flood on dissolved and suspended chemical fluxes in the Skaftá river, Iceland. *J. Volcanol. Geotherm. Res.* **301**, 253–276 (2015).

22. Gussone, N. et al. Calcium isotope fractionation in calcite and aragonite. *Geochim. Cosmochim. Acta* **69**, 4485–4494 (2005).

23. Blättler, C. L., Jenkyns, H. C., Reynard, L. M. & Henderson, G. M. Significant increases in global weathering during Oceanic Anoxic Events 1a and 2 indicated by calcium isotopes. *Earth Planet. Sci. Lett.* **309**, 77–88 (2011).

24. Reynard, L. M., Day, C. C. & Henderson, G. M. Large fractionation of calcium isotopes during cave-analogue calcium carbonate growth. *Geochim. Cosmochim. Acta* **75**, 3726–3740 (2010).

25. Halicz, L., Galy, A., Belshaw, N. S. & O'Nions, K. High-precision measurement of calcium isotopes in carbonates and related materials by multiple collector inductively coupled plasma mass spectrometry (MC-ICP-MS). *J. Anal. Spectrom.* **14**, 1835–1838 (1999).

26. Perez-Fernandez, A., Berninger, U.-N., Mavromatis, V., Pogge von Strandmann, P. A. E. & Oelkers, E. H. Ca and Mg isotope fractionation during the stoichiometric dissolution of dolomite at temperatures from 51 to 126 °C and 5 bars CO₂ pressure. *Chem. Geol.* **467**, 76–88 (2017).

27. Hippler, D. et al. Calcium isotopic composition of various reference materials and seawater. *Geostand. Newsl.* **27**, 13–19 (2003).

28. Pogge von Strandmann, P. A. E., Olsson, J., Luu, T.-H., Gislason, S. R. & Burton, K. W. Using Mg isotopes to estimate natural calcite compositions and precipitation rates during the 2010 Eyjafjallajökull eruption. *Front. Earth Sci.* **7**, https://doi.org/10.3389/feart.2019.00006 (2019).

## Acknowledgements

Analyses and PPvS were funded by NERC Advanced Fellowship NE/1020571/2 and ERC Consolidator grant CONTROLPASTCO2 682760. Chris Coath (Bristol University) is thanked for discussions and formulae on error propagation. The CarbFix pilot infrastructure, injection, sampling, sample distribution and interpretations were funded by the European Commission through the projects CarbFix (EC coordinated action 283148); Min-GRO (MC-RTN-35488); Nordic fund 11029-NORDICCS; the Icelandic GEORG

Geothermal Research fund (09-02-001); the U.S. Department of Energy under award number DE-FE0004847) and Reykjavik Energy. CarbFix has further received grants from the European Union's Horizon 2020 research and innovation program under grant agreements No. 764760 (CarbFix2) and 764810(S4CE). We are indebted to Matin Stute and Jennifer Hall at Columbia University and Juerg Matter at Columbia University and Southampton University; Einar Örn Thrastarson, Trausti Kristinsson, Vordis Eiriksdottir, Halldor Bergmann and Thorsteinn A. Thorgeirsson at Reykjavik Energy; Vigdis Harðardottir, Finnbogi Oskarsson, Kristjan Hrafn Sigurðsson and Steindor Nielsson at ISOR; and Thorsteinn Jonsson, Sveinbjörn Steinthorsson, Iwona Galezcka, Eydis S. Eiriksdottir, Deirdre Clark, Chris Grimm and Flora Brocza at the University of Iceland for helping the injection and sampling campaign.

## Author contributions

P.A.E.P.v.S. carried out the analyses, interpreted the data and wrote the manuscript, with contributions from K.W.B., S.O.S., B.S., E.S.A., I.G., H.A.A., K.G.M., E.H.O. and S.R.G., who are also all involved in the CarbFix project. S.O.S. and S.R.G. provided the samples.

## Additional information

Supplementary Information accompanies this paper at https://doi.org/10.1038/s41467-019-10003-8.

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

**Reprints and permission** information is available online at http://npg.nature.com/reprintsandpermissions/

**Journal Peer Review Information:** *Nature Communications* thanks Fei Wang and other anonymous reviewers for their contribution to the peer review of this work. Peer reviewer reports are available

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

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

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

© The Author(s) 2019