The hydromagnesite playas of Atlin, British Columbia, Canada: A biogeochemical model for CO2 sequestration

Anthropogenic greenhouse gas emissions may be offset by sequestering carbon dioxide (CO₂) through the carbonation of magnesium silicate minerals to form magnesium carbonate minerals. The hydromagnesite [Mg₅(CO₃)₄(OH)₂·4H₂O] playas of Atlin, British Columbia, Canada provide a natural model to examine mineral carbonation on a watershed scale. At near surface conditions, CO₂ is biogeochemically sequestered by microorganisms that are involved in weathering of bedrock and precipitation of carbonate minerals. The purpose of this study was to characterize the weathering regime in a groundwater recharge zone and the depositional environments in the playas in the context of a biogeochemical model for CO₂ sequestration with emphasis on microbial processes that accelerate mineral carbonation.Regions with ultramafic bedrock, such as Atlin, represent the best potential sources of feedstocks for mineral carbonation. Elemental compositions of a soil profile show significant depletion of MgO and enrichment of SiO₂ in comparison to underlying ultramafic parent material. Polished serpentinite cubes were placed in the organic horizon of a coniferous forest soil in a groundwater recharge zone for three years. Upon retrieval, the cube surfaces, as seen using scanning electron microscopy, had been colonized by bacteria that were associated with surface pitting. Degradation of organic matter in the soil produced chelating agents and acids that contributed to the chemical weathering of the serpentinite and would be expected to have a similar effect on the magnesium-rich bedrock at Atlin. Stable carbon isotopes of groundwater from a well, situated near a wetland in the southeastern playa, indicate that  12% of the dissolved inorganic carbon has a modern origin from soil CO₂.The mineralogy and isotope geochemistry of the hydromagnesite playas suggest that there are three distinct depositional environments: (1) the wetland, characterized by biologically-aided precipitation of carbonate minerals from waters concentrated by evaporation, (2) isolated wetland sections that lead to the formation of consolidated aragonite sediments, and (3) the emerged grassland environment where evaporation produces mounds of hydromagnesite. Examination of sediments within the southeastern playa–wetland suggests that cyanobacteria, sulphate reducing bacteria, and diatoms aid in producing favourable geochemical conditions for precipitation of carbonate minerals.The Atlin site, as a biogeochemical model, has implications for creating carbon sinks that utilize passive microbial, geochemical and physical processes that aid in mineral carbonation of magnesium silicates. These processes could be exploited for the purposes of CO₂ sequestration by creating conditions similar to those of the Atlin site in environments, artificial or natural, where the precipitation of magnesium carbonates would be suitable. Given the vast quantities of Mg-rich bedrock that exist throughout the world, this study has significant implications for reducing atmospheric CO₂ concentrations and combating global climate change.     

A calculated petrogenetic grid for ultramafic rocks in the system CaO-FeO-MgO-Al2O3-SiO2-CO2-H2O at low pressures

Calculated phase equilibria among the minerals amphibole, chlorite, clinopyroxene, orthopyroxene, olivine, dolomite, magnesite, serpentine, brucite, calcite, quartz and fluid are presented for the system CaO–FeO–MgO–Al₂O₃–SiO₂–CO₂–H₂O (CaF-MASCH), with chlorite and H₂O–CO₂ fluid in excess and for a temperature range of 440°C–600°C and low pressures. The minerals chosen in CaFMASCH represent the great majority of phases encountered in metamorphosed ultramafic rocks. The changes in mineral compositions in terms of FeMg_-1 and (Mg, Fe)SiAl_-1Al_-1 are related to variations in the intensive parameters. For example, equilibria at high in the presence of chlorite involve minerals which are relatively aluminous compared with those at low . The calculated invariant, univariant and divariant equilibria are compared with naturally-occurring greenschist and amphibolite facies ultramafic mineral assemblages. The correspondence of sequences of mineral assemblages and the compositions of the minerals in the assemblages is very good.     

Secondary mineral phases of metallic lead in soils of shooting ranges from Örebro County,Sweden

In some countries Pb-containing shotgun pellets have become one of the sources of soil contamination in shooting range areas. Pb pellets from eight shooting ranges in central Sweden were mineralogically analysed and the results show that when the Pb pellets come into contact with soil, about 10% of them are decomposed and transformed into secondary lead minerals as encrustations. The encrustation consists of two concentric rims: a 50 to 150 m wide outer rim of hydrocerussite (Pb₃(CO₃)₂(OH)₂) and a 10 to 30 m wide inner-rim of massicot (PbO). Anglesite (PbSO₄) occurs locally in the inner rim. The growing relationship between lead mineral phases suggests that replacement took place. The podzols of the shooting ranges studied are favorable for the formation of lead carbonate. Lead carbonate provides effective controls on the retention of lead in the upper soil layers.     

TiO2 activity in metamorphosed pelitic and basic rocks: principles and applications to metamorphism in southeastern Canadian Cordillera

The activity of TiO₂ can be precisely defined as a function of pressure, temperature and activities of other components for common mineral assemblages in metapelites (ilmenite-quartz-garnet-plagioclase-Al₂SiO₅) and in metabasites (plagioclase-sphene-ilmenite-quartzgarnet). These mineral assemblages can be modelled by the equilibria: 1) 3ilmenite+Al₂SiO₅+2quartz=almandine+3TiO₂ 2) anorthite + 2sphene = grossular + 2TiO₂ + quartz 3) 3anorthite+3quartz+6ilmenite = grossular+ 6TiO₂+2almandine. These mineral assemblages can be used at (rutile saturation) and a given T to get maximum pressure limits of some metapelites and metabasites. When electron microprobe analyses of mineral grains adjacent to Ti-bearing phases are made, these data give maximum pressure estimates in reasonable agreement with other geobarometers. The activity of TiO₂ in many metapelites is very near rutile saturation, but for metabasites the activity of TiO₂ in some sillimanite zone rocks is as low as 0.6. The solubility of TiO₂ in biotite, hornblende and garnet is a complex function of T, P, the activities of components in coexisting minerals and crystal chemical constraints in these minerals. At a given P and T the solubility of TiO₂ in biotite and hornblende does not appear to be strongly dependent upon for sphene and ilmenite versus rutile-bearing assemblages.     

Isotopic fractionations associated with phosphoric acid digestion of carbonate minerals: Insights from first-principles theoretical modeling and clumped isotope measurements

Phosphoric acid digestion has been used for oxygen- and carbon-isotope analysis of carbonate minerals since 1950, and was recently established as a method for carbonate ‘clumped isotope’ analysis. The CO₂ recovered from this reaction has an oxygen isotope composition substantially different from reactant carbonate, by an amount that varies with temperature of reaction and carbonate chemistry. Here, we present a theoretical model of the kinetic isotope effects associated with phosphoric acid digestion of carbonates, based on structural arguments that the key step in the reaction is disproportionation of H₂CO₃ reaction intermediary. We test that model against previous experimental constraints on the magnitudes and temperature dependences of these oxygen isotope fractionations, and against new experimental determinations of the fractionation of ¹³C-¹⁸O-containing isotopologues (‘clumped’ isotopic species). Our model predicts that the isotope fractionations associated with phosphoric acid digestion of carbonates at 25 °C are 10.72‰, 0.220‰, 0.137‰, 0.593‰ for, respectively, ¹⁸O/¹⁶O ratios (1000 lnα^∗) and three indices that measure proportions of multiply-substituted isotopologues . We also predict that oxygen isotope fractionations follow the mass dependence exponent, λ of 0.5281 (where ). These predictions compare favorably to independent experimental constraints for phosphoric acid digestion of calcite, including our new data for fractionations of ¹³C-¹⁸O bonds (the measured change in Δ₄₇ = 0.23‰) during phosphoric acid digestion of calcite at 25 °C.We have also attempted to evaluate the effect of carbonate cation compositions on phosphoric acid digestion fractionations using cluster models in which disproportionating H₂CO₃ interacts with adjacent cations. These models underestimate the magnitude of isotope fractionations and so must be regarded as unsucsessful, but do reproduce the general trend of variations and temperature dependences of oxygen isotope acid digestion fractionations among different carbonate minerals. We suggest these results present a useful starting point for future, more sophisticated models of the reacting carbonate/acid interface. Examinations of these theoretical predictions and available experimental data suggest cation radius is the most important factor governing the variations of isotope fractionation among different carbonate minerals. We predict a negative correlation between acid digestion fractionation of oxygen isotopes and of ¹³C-¹⁸O doubly-substituted isotopologues, and use this relationship to estimate the acid digestion fractionation of for different carbonate minerals. Combined with previous theoretical evaluations of ¹³C-¹⁸O clumping effects in carbonate minerals, this enables us to predict the temperature calibration relationship for different carbonate clumped isotope thermometers (witherite, calcite, aragonite, dolomite and magnesite), and to compare these predictions with available experimental determinations. The success of our models in capturing several of the features of isotope fractionation during acid digestion supports our hypothesis that phosphoric acid digestion of carbonate minerals involves disproportionation of transition state structures containing H₂CO₃.     

Molecular models of a hydrated calcite mineral surface

Hydrated mineral surfaces play an important role in many processes in biological, geological, and industrial applications. An energy force field was developed for molecular mechanics and molecular dynamics simulations of hydrated carbonate minerals and was applied to investigate the behavior of water on the calcite surface. The force field is a significant development for large-scale molecular simulations of these systems, and provides good agreement with experimental and previous modeling results. Simulations indicate that water molecules are significantly ordered near the calcite surface. The predominant surface configuration (75-80%) results from coordination of a water molecule with a single calcium cation-carbonate anion pair, while the less common situation involves water coordination with two ion pairs. Surface restructuring and variation in coordination in the water layers results in distinct distances for water oxygens above the calcite surface—a two-component first monolayer (2.3 and 3.0 Å) and a secondary monolayer (5.0 Å). The different coordinations also alter lateral displacement, hydrogen bonding, and surface-normal orientation of the water molecules. The ordering of water molecules and the formation of a unique hydrogen bonding network at the calcite surface influence the physical properties of the interfacial water. Surface exchange of water molecules is observed by molecular dynamics simulation to occur at a rate of one exchange per 10 ps. Diffusion coefficients derived from mean square displacement analysis of atomic trajectories indicate a dependence of water transport based on the distance of the water molecules from the calcite surface.     

On the neutralization of acid rock drainage by carbonate and silicate minerals

The net result of acid-generating and-neutralizing reactions within mining wastes is termed acid rock drainage (ARD). The oxidation of sulfide minerals is the major contributor to acid generation. Dissolution and alteration of various minerals can contribute to the neutralization of acid. Definitions of alkalinity, acidity, and buffer capacity are reviewed, and a detailed discussion of the dissolution and neutralizing capacity of carbonate and silicate minerals related to equilibium conditions, dissolution mechanism, and kinetics is provided. Factors that determine neutralization rate by carbonate and silicate minerals include: pH, P_{CO 2}, equilibrium conditions, temperature, mineral composition and structure, redox conditions, and the presence of foreign ions. Similar factors affect sulfide oxidation. Comparison of rates shows sulfides react fastest, followed by carbonates and silicates. The differences in the reaction mechanisms and kinetics of neutralization have important implications in the prediction, control, and regulation of ARD. Current static and kinetic prediction methods upon which mine permitting, ARD control, and mine closure plans are based do not consider sample mineralogy or the kinetics of the acid-generating and-neutralizing reactions. Erroneous test interpretations and predictions can result. The importance of considering mineralogy for site-specific interpretation is highlighted. Uncertainty in prediction leads to difficulties for the mine operator in developing satisfactory and cost-effective control and remediation measures. Thus, the application of regulations and guidelines for waste management planning need to beflexible.     

Structure and dynamics of the water films confined between edges of pyrophyllite: A first principle study

Edge sites of clay minerals play a key role for pH dependent sorption of ions from solutions of electrolytes. Pyrophyllite, Al₂[Si₄O₁₀](OH)₂, is an important structural prototype for a variety of 2:1 dioctahedral phyllosilicates but in contrast to the other clays has no permanent structural charge. The structure of thin water films confined between most common edges of 1Tc pyrophyllite: (0 1 0), (1 1 0) and (1 0 0), was analyzed by means of ab initio molecular dynamic simulations. The system setup allowed for a full flexibility of the interfaces and a proton exchange between the edges of pyrophyllite and water molecules in solution. The structure of hydrated surfaces is compared with the recent predictions of static geometry optimizations for edge-vacuum interfaces. All surfaces studied reveal a strong hydrophilic character of edge similar to the hydrated silica surface and the facets of simple layered hydroxides. Spontaneous proton transfer between different surface sites were observed in molecular dynamics simulations of the (0 1 0) interface. The proton bound to the SiOH site was found to exchange with the AlOH group by the mechanism . The direction of the proton transfer agrees with the scale of relative proton affinities for surface sites obtained from the static calculations. Alternatively, the proton attached to the AlOH₂ site exchanges with the AlOH group. In both reactions, the protons are transferred through the chains of hydrogen bonds formed between water molecules in the solution and the surface sites. The observed mechanisms might be one of the basic schemes for the surface proton diffusion in compacted clays. Kinetics of the proton transfer at edge sites is limited by the rate of rearrangements of the water molecules near interface.     

High arsenic contents in groundwater of central Spain

The chemical and microbiological characteristics of groundwater from two provinces of central Spain were studied. In some zones of this area, the concentrations of As in groundwater exceed the guideline concentrations, set internationally between 10 g/l and 50 g/l, reaching levels over 100 g/l. A narrow correlation between the contents of arsenic and HCO₃^– was observed. These data suggest a possible mechanism of the As mobilization from aquifer sediments to groundwater: the bicarbonate ions could displace HAsO₄^2– adsorbed on aquifer oxyhydroxides. Sediments containing relatively high contents of adsorbed arsenic are deposited in surface water environments with low carbonate concentrations. Subsequently, the sediments become exposed to groundwater with highly dissolved carbonate content, and arsenic can be mobilized by displacement from mineral surfaces. In addition, the presence of Pseudomonas genera bacteria, which secrete siderophores (Fe chelating agents) could mobilize As adsorbed on Fe oxides through their dissolution. These combined microbiological and chemical processes might have increased the natural mobility of As.     

The redox chemistry of Pu(V)O2+ interaction with common mineral surfaces in dilute solutions and seawater

It has long been recognized that the Pu⁴⁺ ion can be readily adsorbed on solid surfaces, but it has been assumed that the generally more abundant Pu(V)O₂⁺ ion should have little affinity for surfaces. Our results indicate that Pu(V)O₂⁺ can be adsorbed from dilute solutions and seawater on goethite, aragonite, calcite, and δ-MnO₂. Adsorption on δ-MnO₂ is severely depressed in seawater, probably as a result of site competition with seawater cations.The sorption behavior of PuO₂⁺ is influenced by oxidation-reduction reactions occurring on the mineral surfaces. Adsorption on δ-MnO₂ results in oxidation of adsorbed Pu(IV) and Pu(V) to Pu(VI). However, adsorption on goethite results in a reaction in which Pu(IV) and Pu(VI) are formed on the mineral surface. The Pu(VI) is slowly reduced to Pu(IV), leaving Pu(IV) as the dominant surface Pu species. This reaction can be photochemically catalyzed. PuO₂⁺ adsorbed on carbonate minerals behaves similarly to Np(V)O₂⁺ and undergoes little change in oxidation state after adsorption.     

P CO2 in low-grade metamorphism; zeolite,carbonate, clay mineral,prehnite relations in the system CaO-Al2O3-SiO2-CO2-H2O

Equilibria for several reactions in the system CaO-Al₂O₃-SiO₂-CO₂-H₂O have been calculated from the reactions calcite+quartz=wollastonite+CO₂ (5) and calcite+Al₂SiO₅+quartz=anorthite+CO₂ (19) and other published experimental studies of equilibria in the systems Al₂O₃-SiO₂-H₂O and CaO-Al₂O₃-SiO₂-H₂O.The calculations indicate that the reactions laumontite+CO₂=calcite+kaolinite+2 quartz+2H₂O (1) and laumontite+calcite=prehnite+quartz+3H₂O+CO₂ (3) in the system CaO-Al₂O₃-SiO₂-CO₂-H₂O, are in equilibrium with an H₂O-CO₂ fluid phase having -0.0075 for P _fluid=P _total=2000 bars.These calculations limit the stability of zeolite assemblages to low p CO₂.Using the above reactions as model equilibria, several probelms of p CO₂ in low grade metamorphism are discussed. (a) the problem of producing zeolitic minerals from metasedimentary assemblages of carbonate, clay mineral, quartz. (b) the significance of calcite (or aragonite) associated with zeolite (or lawsonite) in low grade metamorphism and hydrothermal alteration. (c) the reaction of zeolites (or lawsonite) with calcite (or aragonite) to produce dense Ca-Al-hydrosilicates (eg. prehnite, zoisite, grossular).     

Timescales of interface-coupled dissolution-precipitation reactions on carbonates

In the Earth's upper crust, where aqueous fluids can circulate freely, most mineral transformations are controlled by the coupling between the dissolution of a mineral that releases chemical species into the fluid and precipitation of new minerals that contain some of the released species in their crystal structure, the coupled process being driven by a reduction of the total free-energy of the system. Such coupled dissolution-precipitation processes occur at the fluid-mineral interface where the chemical gradients are highest and heterogeneous nucleation can be promoted, therefore controlling the growth kinetics of the new minerals. Time-lapse nanoscale imaging using Atomic Force Microscopy (AFM) can monitor the whole coupled process under in situ conditions and allow identifying the time scales involved and the controlling parameters. We have performed a series of experiments on carbonate minerals (calcite, siderite, dolomite and magnesite) where dissolution of the carbonate and precipitation of a new mineral was imaged and followed through time. In the presence of various species in the reacting fluid (e. g. antimony, selenium, arsenic, phosphate), the calcium released during calcite dissolution binds with these species to form new minerals that sequester these hazardous species in the form of a stable solid phase. For siderite, the coupling involves the release of Fe²⁺ ions that subsequently become oxidized and then precipitate in the form of Fe^III oxyhydroxides. For dolomite and magnesite, dissolution in the presence of pure water (undersaturated with any possible phase) results in the immediate precipitation of hydrated Mg-carbonate phases. In all these systems, dissolution and precipitation are coupled and occur directly in a boundary layer at the carbonate surface. Scaling arguments demonstrate that the thickness of this boundary layer is controlled by the rate of carbonate dissolution, the equilibrium concentration of the precipitates and the kinetics of diffusion of species in a boundary layer. From these parameters a characteristic time scale and a characteristic length scale of the boundary layer can be derived. This boundary layer grows with time and never reaches a steady state thickness as long as dissolution of the carbonate is faster than precipitation of the new mineral. At ambient temperature, the surface reactions of these dissolving carbonates occur on time-scales of the order of seconds to minutes, indicating the rapid surface rearrangement of carbonates in the presence of aqueous fluids. As a consequence, many carbonate-fluid reactions in low temperature environments are controlled by local thermodynamic equilibria rather than by the global equilibrium in the whole system.     

Tonminerale in normalsedimenten,hydrothermal beeinflu\ten sedimenten und Erzschlämmen des Roten Meeres

From 14 deeps and other regions of the Red Sea totally 226 samples from 28 cores recovered during the VALDIVIA cruises (1971, 1972) were investigated according to their clay mineral content (<2μm resp. < 6.3 μm) after carbonate dissolution. Three facies groups are to distinguish:

1. normal sediments: dominance of chlorite, kaolinite, illite, small amounts of smectite and sepiolite. Two palygorskite types are present only in a few samples.
2. normal sediments with hydrothermal influence: clay mineral paragenesis similar like that of normal sediments; but increase of smectite and presence of goethite in each sample; partly small contents of talc.
3. heavy metal deposits: dominance of iron-bearing smectite, partly with amorphous components resp. pure ore mineral assemblages with authigenic silicates (talc, quartz, opal, chrysotile, sepiolite, palygorskite, chlorite).

Crystallinity of the clay and ore minerals is independent from sedimentary overburden. Sepiolite shows in small amounts a wide distribution; palygorskite₂ (d₁₁₀=11.3 Å) yields locally an increased concentration in the range of pteropod layers cemented by aragonite. The environment of ore deposits is characterized by iron-bearing smectite besides the ore minerals.     

Arsenate uptake by calcite: Macroscopic and spectroscopic characterization of adsorption and incorporation mechanisms

Batch uptake experiments and X-ray element mapping and spectroscopic techniques were used to investigate As(V) (arsenate) uptake mechanisms by calcite, including adsorption and coprecipitation. Batch sorption experiments in calcite-equilibrated suspensions (pH 8.3; PCO₂ = 10^−3.5 atm) reveal rapid initial sorption to calcite, with sorption rate gradually decreasing with time as available sorption sites decrease. An As(V)-calcite sorption isotherm determined after 24 h equilibration exhibits Langmuir-like behavior up to As concentrations of 300 μM. Maximum distribution coefficient values (K_d), derived from a best fit to a Langmuir model, are ∼190 L kg⁻¹.Calcite single crystals grown in the presence of As(V) show well-developed rhombohedral morphology with characteristic growth hillocks on surfaces at low As(V) concentrations (?5 μM), but habit modification is evident at As(V) concentrations ?30 μM in the form of macrostep development preferentially on the − vicinal surfaces of growth hillocks. Micro-X-ray fluorescence element mapping of surfaces shows preferential incorporation of As in the − vicinal faces relative to + vicinals. EXAFS fit results for both adsorption and coprecipitation samples confirm that As occurs in the 5+ oxidation state in tetrahedral coordination with oxygen, i.e., as arsenate. For adsorption samples, As(V) forms inner-sphere surface complexes via corner-sharing with Ca octahedra. As(V) coprecipitated with calcite substitutes in carbonate sites but with As off-centered, as indicated by two Ca shells, and with likely disruption of local structure. The results indicate that As(V) interacts strongly with the calcite surface, similar to often-cited analog phosphate, and uptake can occur via both adsorption and coprecipitation reactions. Therefore, calcite may be effective for partial removal of dissolved arsenate from aquatic and soil systems.     

Uptake of Cu2+ by the calcium carbonates vaterite and calcite as studied by continuous wave (cw) and pulse electron paramagnetic resonance

The investigation of the Cu²⁺ uptake by the calcium carbonate minerals vaterite and calcite with continuous wave and pulse electron paramagnetic resonance (EPR) yields information on a molecular scale about the relevant complexation reactions at the mineral–water interface. The structural assignment is based on changes in the coordination geometry of the copper complexes. Magnetic interactions of the unpaired Cu²⁺ electron with nuclei of ¹³C-labeled carbonate ligands and protons from water or hydroxyl ligands in the first and second coordination spheres of the cation are detected by pulse EPR techniques. Our results show that the Cu²⁺ ions are rapidly dehydrated upon adsorption on the mineral surface. The strong surface binding is due to monodentate coordination to three or four carbonate surface ions, comparable to chelate complexation in solution. The formation of square-planar or square-pyramidal copper complexes at exposed surface sites like kinks and steps yields a convincing explanation for the inhibition of calcium carbonate growth and dissolution. Upon recrystallization the Cu²⁺ ions are integrated into the calcite lattice where they exhibit a dynamic Jahn–Teller effect. The resulting local lattice distortions are expected to destabilize the Cu_xCa_(1−x)CO₃(s) solid solution. Our results support the concept of a dynamic calcium carbonate surface, covered by a thin, structured surface layer. The detailed structural information obtained for Cu²⁺ provides a better understanding of the interaction of other metal ions with calcium carbonate minerals.     

Can the long-term potential for carbonatization and safe long-term CO2 storage in sedimentary formations be predicted?

A sedimentary formation perturbated by supercritical CO₂ reacts by dissolving primary minerals and forming new secondary phases. In this process CO₂ may be trapped in stable carbonate minerals and may thereby be immobilized for long time spans. The potential for mineral trapping can be estimated by solving kinetic expressions for the reservoir minerals and possible secondary phases. This is, however, not trivial as kinetic data are uncertain or even lacking for the minerals of interest. Here, the rate equations most commonly used for CO₂ storage simulations have been solved, and the rate parameters varied, to obtain sensitivity on the total amount of CO₂ stored as mineral carbonate. As various expressions are in use to estimate growth rates of secondary carbonates, three formulations were compared, including one taking into account mineral nucleation preceding growth. The sensitivity studies were done on two systems, the Utsira Sand being representative for a cold quartz-rich sand (37 °C, 100 bar CO₂), and the Gulf Coast Sediment, being representative for a medium temperature quartz–plagioclase-rich system (75 °C, 300 bar CO₂).The simulations showed that the total predicted CO₂ mineral storage is especially sensitive to the choice of growth rate model and the reactive surface area. The largest sensitivity was found on α, fraction of total surface area available for reactions, with a reduction of one order of magnitude for all reacting phases leading to 3–4 times lower predicted CO₂ mineral storage. Because the reactive surface area is highly uncertain for natural systems, the range in predicted results may be even larger. The short-term predictions (<100–1000 a), such as the onset of carbonate growth, were highly sensitive to nucleation and growth rates. Moreover, the type of carbonate minerals formed was shown to be model dependent, with the simplest model predicting an unlikely carbonate assemblage at low temperature (i.e., formation of dolomite at 37 °C). Therefore, to use kinetic models to upscale short-term (<months) laboratory experiments in time, to identify the past reactions and physical conditions of natural CO₂ storage analogues, and finally to predict the potential for CO₂ trapping in existing and future storage projects, more knowledge has to be collected, especially on the reactive surface area of CO₂ storage reservoirs, and on the rate of secondary carbonate nucleation and growth.