The direct crystallization of dolomite from an aqueous solution at temperatures between 60–220 °C was followed in situ through time-resolved synchrotron-based energy-dispersive X-ray diffraction combined with offline high-resolution imaging, X-ray diffraction, and infrared spectroscopy. Crystalline CaMg(CO3)2 phases form through a three-stage process. In the first stage, a nanoparticulate magnesium-deficient, amorphous calcium carbonate (Mg-ACC) with a nominal formula of Ca0.606Mg0.394CO3·1.37H2O forms. After a temperature-dependent induction time, during stage 2 the Mg-ACC partially dehydrates and orders prior to its rapid (<5 min) crystallization to non-stoichiometric proto-dolomite. This occurs via the dissolution of Mg-ACC, followed by the secondary nucleation of proto-dolomite from solution. The proto-dolomite crystallization proceeds via spherulitic growth that follows a growth front nucleation mechanism with a de-nuovo and continuous formation of nanocrystalline proto-dolomite subunits that form spherical aggregates. In stage three of the reaction, the proto-dolomite transforms to highly crystalline and stoichiometric dolomite on a much longer timescale (hours to days), via an Ostwald-ripening mechanism. Such a three-stage crystallization can explain microbially induced proto-dolomites observed in modern hypersaline settings and may also be the route by which the Cryogenian cap dolomite deposits of the Neoproterozoic formed.
Dolomite, CaMg(CO3)2, is one of the most common carbonate minerals in the geologic record (Machel 2004) and the most oversaturated carbonate mineral in the modern ocean (Tucker and Wright 1990). However, dolomite rarely forms in modern environmental systems (MacKenzie and Andersson 2013), which has led to considerable interest in understanding its formation pathway. There are two main mechanisms by which dolomite can form in the natural environment: (1) “directly from solution” or primary dolomite, and (2) dolomites that formed during diagenesis or burial, most often through secondary replacement (dolomitization) of previously precipitated calcite or aragonite (Machel 2004; Mckenzie and Vasconcelos 2009). The largest dolomite deposits in the geologic record (Machel 2004) are in the Phanerozoic and are diagenetic in origin. In modern environments dolomite currently forms through the “direct” route, but only in a restricted number of hypersaline environments [i.e., evaporitic lakes or sabkhas; e.g., Jones (1965), Lindtke et al. (2011), Meister et al. (2011)]. In such settings, microbial metabolic processes (predominantly sulfate reduction or methanogenesis) induce high alkalinity. Together with the moderate to high Mg:Ca concentration ratios (1:1 to 100:1), this creates high supersaturation levels that commonly lead to the formation of “dumbbell shaped” aggregates of dolomite (Warthmann et al. 2000; Sánchez-Román et al. 2011). A similar “direct” formation mechanism, which also produces a radial growth morphology, has been suggested for dolomite cements in Cryogenian (~850 to ~650 Mya) reef carbonates (Hood et al. 2011; Hood and Wallace 2012). The presence of directly precipitated dolomite cements during this period, has led to the oceans being classified as “aragonite-dolomite-seas.” In contrast the Phanerozoic oceans are classified as “calcite-aragonite seas” (Sandberg 1983), with the abundant dolomite found in Phanerozoic sediments having formed via various diagenetic and burial processes, but not through direct precipitation from seawater. The switch from dominantly “direct” dolomite formation to indirect “dolomitization” as the formation pathway may be related to the changes in ocean chemistry over time. The abundance of dolomite in the Phanerozoic, yet its scarcity in modern systems is not easily explained because we lack a fundamental understanding of how low-temperature dolomites form in marine settings and whether such formation reactions lead to primary or early secondary diagenetic dolomites. This has lead to the paradox termed the “Dolomite Problem,” which is a long-standing debate with important mineralogical and palaeo-environmental implications, as it links the formation of dolomite with the long-term evolution of Earth’s ocean chemistry. However, there are key knowledge gaps in this field, which limit a thorough interpretation of dolomite formation in the geological record. One of these is a lack of a mechanistic understanding of the “direct” dolomite formation pathway, and the factors (e.g., temperature, supersaturation, etc.) that may favor or inhibit this process.
This is particularly relevant because, as yet, inorganic dolomite has not been synthesized at ambient conditions directly from solution (Land 1998). Most experimental inorganic dolomite formation studies have quantified the “dolomitization” of calcite (Gaines 1977; Katz and Matthews 1977; Miura and Kawabe 2000; Nordeng and Sibley 1994; Kaczmarek and Sibley 2007), other Ca-Mg carbonates (Katz and Matthews 1977; Zempolich and Baker 1993; Sibley et al. 1994), or followed proto-dolomite crystallization (Malone et al. 1996) or dolomite growth (Gaines 1977; Arvidson and Mackenzie 1999; Higgins and Hu 2005; Hu et al. 2005, 2006) on various carbonate seed crystals. These experiments all mimic the diagenetic dolomite formation processes. Only a few experimental studies have attempted to quantify the fundamental reactions controlling the formation of inorganic dolomite “directly” from solution (Ohde and Kitano 1978; Kelleher and Redfern 2002; Schmidt et al. 2005). These inorganic experiments showed that a crystalline carbonate [but typically not stoichiometric CaMg(CO3)2] formed from an amorphous Ca-Mg carbonate precursor. This precursor precipitated from highly supersaturated solutions, and transformed to dolomite via proto-dolomite at moderate to high temperatures (e.g., 40–200 °C; Kelleher and Redfern 2002; Schmidt et al. 2005). Proto-dolomite has no ordering of Mg and Ca within its structure and often forms micrometer-sized spheres made up of aggregates of small particles (<200 nm; e.g., Malone et al. 1996). Such proto-dolomites are usually distinguished from crystalline dolomite by the lack of superstructure peaks in the powder diffraction pattern (Kelleher and Redfern 2002). Malone et al. (1996) showed that at 50–200 °C and over long time periods proto-dolomite transforms to crystalline dolomite via a dissolution and reprecipitation reaction. This transformation lead to dolomite crystals with well developed crystallographic faces, but significant crystallographic ordering of the Ca and Mg was only achieved at 200 °C. Several microbial studies (e.g., Kenward et al. 2009; Sánchez-Román et al. 2011) have also shown the formation of proto-dolomite and dolomite at ambient conditions, but it is unclear whether an amorphous carbonate precursor formed prior to crystallization. These inorganic and microbial studies have conceptually described the multi-step crystallization of proto-dolomite and dolomite, however, the mechanisms by which these reactions take place have not yet been quantified.
Comparing the direct dolomite formation pathway with that of calcite also from an amorphous precursor (amorphous calcium carbonate, ACC) reveals many similarities. Magnesium-free ACC crystallizes rapidly to vaterite via spherulitic growth (Bots et al. 2012), followed by the slower transformation of vaterite to calcite via dissolution and re-precipitation (Rodriguez-Blanco et al. 2011). Spherulitic crystal growth is common in various systems [e.g., polymers supercooled from a molten state, Keith and Padden (1963); crystallization of viscous magmas, Lofgren (1971); see also review of Shtukenberg et al. (2012) and references therein]. Spherulitic growth has been shown to proceed via a growth front nucleation mechanism (Gránásy et al. 2005; Shtukenberg et al. 2012), which requires a high crystallization driving force (i.e., high supersaturation) that leads to the formation of spherical nanocrystal aggregates. For the carbonate system it has been demonstrated that the amorphous ACC precursor transforms to vaterite via spherulitic growth (Beck and Andreassen 2010; Andreassen et al. 2010; Bots et al. 2012), and that this reaction occurs at high superaturations and leads to spheroidal vaterite aggregates. Interestingly, such spheroidal morphologies are similar to those observed for proto-dolomite/dolomite formed in microbial systems (i.e., spheroidal or dumbbell-shaped dolomite aggregates) suggesting that microbially induced Ca-Mg-carbonates may have also formed via spherulitic growth (Warthmann et al. 2000; Sánchez-Román et al. 2011). Furthermore, the radial growth morphology of the “direct” dolomite cements in the Neoproterozoic cap carbonates may also be a product of spherulitic growth indicating that they also may have formed from an amorphous precursor (Hood et al. 2011; Hood and Wallace 2012).
Thus, we hypothesize that an amorphous precursor is a necessary prerequisite for the crystallization of proto-dolomite/dolomite and that the spherulitic crystal growth mechanism observed in the pure CaCO3 system also controls the “direct” formation of proto-dolomite and dolomite in the Mg-Ca-CO3 system. To address this we characterized the crystallization of an amorphous Ca-Mg carbonate precursor that was precipitated from a supersaturated, aqueous Ca-Mg carbonate solution with time via both conventional and in situ, time-resolved synchrotron-based methods. We then used these data to derive kinetic and mechanistic information for the nucleation, growth, and/or crystallization of proto-dolomite and dolomite. We did this to address three primary questions: (1) what is the mechanism of proto-dolomite formation from the amorphous precursor; (2) what is the mechanism of proto-dolomite transformation to crystalline dolomite; and (3) what are the key factors (e.g., saturation states, temperature), which control these mineral formation and crystallization reactions. The answers to these questions will provide a mechanistic basis for the crystallization pathway for the direct proto-dolomite and dolomite formation in both modern and ancient systems, e.g., microbially mediated proto-dolomite formation in sabka environments and Cryogenian reef dolomites of the Neoproterozoic.
Experiments were carried out by mixing 1, 1.5, or 2.5 M Na2CO3, CaCl2, and MgCl2 solutions (ratio of CO32−:Ca2+:Mg2+ = 2:1:1) at room temperature in a stirred, polytetrafluoroethylene-lined hydrothermal reaction cell. Crystallization was then induced by heating the sealed cells to temperatures between 25 and 220 °C for between 30 min and 10 days. Mixing the starting solutions induced the instantaneous formation of a white, gel-like precipitate. The transformation of this white gel was followed online with synchrotron-based diffraction and offline with conventional techniques. The changes in the reaction cell contents were monitored in situ and in a time resolved manner via changes in energy-dispersive X-ray diffraction (ED-XRD) patterns collected on the white-beam station 16.4 at the Synchrotron Radiation Source (SRS), Daresbury Laboratory, U.K. At temperatures between 60 and 220 °C, ED-XRD scans of the crystallizing gel were recorded at 1 min intervals to follow the formation and changes in Bragg peak intensities throughout each reaction for up to 2 h. The reactions were considered complete when the intensity of all Bragg peaks stopped increasing. The ED-XRD patterns were fitted using XFIT (Cheary and Coelho 1992) and the peak areas and the overall background intensity as a function of time calculated and normalized to values between 0 and 1 to express the degree of reactions (α). The increase in peak area with time from each experiment was fitted to various crystallization kinetic models (Johnson and Mehl 1939; Bischoff and Fyfe 1968; Ogino et al. 1990) to derive the kinetic parameters. More detailed descriptions of the online, in situ experimental protocols and data analyses procedures can be found in Cahill et al. (2000), Shaw et al. (2005), and Davidson et al. (2008).
In equivalent offline experiments (with no ED-XRD data collection) the white gel-like precipitate was reacted at 25–220 °C for up to 6 h (1 M initial concentrations) and at 25 to 60 °C for up to 10 days (1, 1.5, and 2.5 M initial concentrations). These offline experiments were quenched and separated from the aqueous solutions by vacuum-filtration (0.2 μm polycarbonate Cyclopore filters) and the solids (initial gel, intermediates, and end-products) were dried with isopropanol following the procedure described in Rodriguez-Blanco et al. (2008). These solids were analyzed with powder X-ray diffraction using a Bruker D8 powder X-ray diffractometer fitted with a GE (111) monochromator (PXRD, CuKα1, λ = 1.541 Å; 2θ = 5–70°; 0.001° step size; 0.5 s dwell time/step). In all cases, a silicon internal standard (Aldrich; ICDD PDF 27-1402) was used for 2θ calibration. From the PXRD patterns, the molar fraction of MgCO3 (χMgCO3) in each crystalline end product was determined using the expression from Arvidson and Mackenzie (1999)
where MMgCO3 and MCaCO3 correspond to the moles of MgCO3 and CaCO3 per formula unit of dolomite, and d(104) corresponds to the d-spacing (Å) of the (104) reflection of dolomite. To determine the magnesium content in the initial white gel-like precipitate, nine replicates of the gel were dried and digested in HCl and the calcium and magnesium concentration analyzed by ion chromatography following Bots et al. (2011). To discriminate between the proto-dolomite and dolomite, high-resolution PXRD scans (0.001° step size; 2 s dwell time/step) of the offline reaction end products were acquired over the 2θ range between 20 to 45°. This way the low-intensity (~1%) dolomite superstructure reflections, (101) and (015) were distinguishable (Goldschmidt et al. 1961; Arvidson and Mackenzie 1999). Finally, from these PXRD patterns, crystallite sizes were calculated via the Scherrer equation (Scherrer 1918) using the software code TOPAS (Coelho 2006).
The vibrational characteristics of the solids were also analyzed using an A2-Technology MicroLab mid-infrared spectrometer, with a diamond internal reflection sampling system. Spectra between 650 to 4000 cm−1 were accumulated by co-adding 512 scans with a 4 cm−1 resolution and spectral manipulation including baseline adjustment, normalization, and band component analysis was carried out using the Thermo-Nicolet OMNIC ESP 5.1 software package.
High-resolution photomicrographs of the solids were acquired using a field emission gun scanning electron microscope (FEG-SEM, LEO 1530 Gemini at 3 keV). From the FEG-SEM images, the particle size distributions in the white gel like precursor phase and the crystalline end products were determined by measuring the diameters of ~100 particles each in representative images. It should be noted that due to the small size of some of the particles (<50 nm) the measurement error from FEG-SEM image evaluations can reach 30%. The water content of the precursor phase was determined by thermogravimetric analyses (Mettler Toledo TGA/DSC1 Thermal Analyzer, 10°/min ramp between 30 and 800 °C) of 2 aliquots of the isopropanol-quenched starting material.
Last, the saturation indexes (SI) of the initial gel-like precipitate and the crystalline proto-dolomite and dolomite were calculated using the geochemical code PHREEQC (Parkhurst 1995) where the SI is defined as
with a corresponding to the activity of the Ca2+, Mg2+, and CO32− ions at various times of the reaction and Ksp corresponding to the solubility product of either dolomite or proto-dolomite (Robie et al. 1978; Carpenter 1980). To calculate the aqueous Ca and Mg concentrations after the formation of the initial gel-like precipitate, we assumed that the solubility of the initial Mg-rich gel was similar to pure amorphous calcium carbonate precursors (ACC; solubility from Brečević and Nielsen 1989). This was based on our own solid analysis and comparing with other carbonate studies in the literature (Davis et al. 2000). This way the aqueous post-gel precipitation solution composition (e.g., Ca2+ and Mg2+) could be calculated using the initial aqueous ion concentrations (0.250 M CaCl2, 0.250 M MgCl2, and 0.500 M Na2CO3), the ACC solubility and the analyzed composition and water content of the initial gel-like precipitate. Full details of the modeling methodology and the assumptions made in this process are available in the Supplementary Information1.
Offline PXRD analysis of the white gel-like phase that formed immediately upon mixing of the starting supersaturated solutions revealed broad diffraction peaks centered at approximately 30°, 45°, and 63° 2θ (Fig. 1, bottom pattern, marked with *) indicating that the initial phase was poorly ordered. The position and width of these peaks are equivalent to those observed for amorphous calcium carbonate (ACC, Rodriguez-Blanco et al. 2008). The magnesium content in the precursor formed by mixing the 1 M solutions (ionic strength of 2.0) was measured to be χMgCO3 = 0.394 (±0.005; n = 9). We will refer hereafter to this precursor phase as Mg-ACC, following the notation of Günther et al. (2005). TGA analyses of our Mg-ACC revealed a water content of 1.37 ± 0.04 mol per mole of Mg-ACC (Table 1). Combining the χMgCO3 and water contents allowed us to calculate a Ca0.606Mg0.394CO3·1.37H2O formula for our Mg-ACC. Based on PHREEQC modeling (Supplementary Information1), the Mg-ACC in our experiments precipitated from solutions with high-ionic strengths (1, 1.5, and 2.5 M) that were highly supersaturated with respect to all Ca-Mg-carbonate phases. For example, starting with 1 M solutions, the SI of ACC at 25 °C was 4.02, while for proto-dolomite and dolomite the SI was 8.62 and 9.16, respectively. Once the Mg-ACC precipitated at 25 °C the SI for proto-dolomite and dolomite dropped to 2.46 and 3.00, respectively (see text and Supplementary Table 11).
The offline experiments at 60 to 200 °C, with initial reactant concentrations of 1 M, showed that after 6 h the Mg-ACC crystallized (upper pattern in Fig. 1) to either proto-dolomite (60 and 100 °C) or dolomite (≥140 °C). Dolomite was differentiated from proto-dolomite through the presence of the (101) and (015) superstructure peaks in the PXRD pattern (pDol and Dol in insets in Fig. 1). Analysis of the high-temperature reaction products (>140 °C) after 30 and 60 min revealed only proto-dolomite peaks, with the superstructure peaks absent. Between 1 and 6 h the growth of the superstructure peaks indicated the transformation from proto-dolomite to dolomite. The initial Mg-ACC (χMgCO3,Mg-ACC = 0.39) transformed to a non-stoichiometric, slightly Mg-deficient proto-dolomite (χMgCO3,pDol = 0.44–0.45), which in turn transformed to a highly crystalline and stoichiometric dolomite (χMgCO3,Dol = ~0.50; Table 1). In contrast at 60 °C, even after 10 days, proto-dolomite did not crystallize to dolomite. Furthermore, at lower temperatures (40 and 25 °C) and the same initial 1 M concentration, proto-dolomite did not form, instead the reaction products were monohydrocalcite (CaCO3·H2O) and minor (~5%) hydromagnesite [Mg5(CO3)4(OH)2·4H2O] (Supplementary Fig. 11). When higher initial concentrations were used (i.e., 1.5 and 2.5 M), proto-dolomite formed at both 50 and 40 °C, yet at 25 °C again only monohydrocalcite and minor hydromagnesite (~5%) formed and both remained stable for long time periods (days).
Data from an in situ and time-resolved ED-XRD experiment of the crystallization of Mg-ACC at 180 °C is shown in Figure 2. At the start of the reaction a large background hump with no distinct diffraction peaks dominated. This background hump is a consequence of X-ray scattering from both the aqueous solution and the poorly ordered Mg-ACC. In all experiments, after a temperature-dependent induction period (Table 1) Bragg peaks at ~31° and 33° 2θ, identified as the (104) and (006) proto-dolomite/dolomite peaks began to grow (Fig. 2). With increasing temperature, the induction time (t0) for the initiation of the crystallization reaction decreased from 48 min at 60 °C to 10 min at 220 °C (Fig. 3; Table 1). Once peaks started to grow, their intensity increased very rapidly (~5 min) and they reached their maximum intensity simultaneously, after which no further changes in intensities were observed. A distinction between dolomite and proto-dolomite was not possible using the ED-XRD data because the low resolution of the energy-dispersive detection system and large background signal from the solution did not allow the low-intensity superstructure peaks to be distinguished. The intensity of the background hump (Fig. 2) in the time resolved ED-XRD data decreased throughout the reaction. This decrease started prior to the formation of the Bragg peaks, and continued until the intensity of the peaks stopped increasing. For example, in the 180 °C experiment (Fig. 4), the first appearance of the (104) peak was observed after 10 min, when the intensity of the background hump had already decreased to α = 0.5. At 15 min the Bragg peaks reached their maximum intensity while the background intensity reached a minimum.
The temperature-dependent crystallization of Mg-ACC described above was also confirmed by the changes in the FTIR spectra (Fig. 5). The spectra of Mg-ACC, proto-dolomite and dolomite all showed the characteristic bands for the carbonate ion (symmetric and asymmetric νx peaks for CO3, Fig. 5) typical for most calcium and magnesium carbonates (Farmer 1974). However, similar to the spectrum for pure ACC (Günther et al. 2005; Rodriguez-Blanco et al. 2011), in the spectrum for the precursor Mg-ACC (bottom pattern in Fig. 5), the ν4 symmetric in-plane bending vibration at 725 cm−1 was absent. The presence of a sharp ν4 vibration is normally an unambiguous signature for crystalline calcium or magnesium carbonates (upper pattern in Fig. 5; White 1974). The second clear distinguishing feature between the precursor Mg-ACC and the two crystalline phases was the presence of bands related to OH vibrations. These are the broad features centered around ~3300 cm−1, corresponding to O-H stretching and the sharper band at 1650 cm−1 corresponding to O-H bending vibrations, both are signatures for structural water. In the precursor Mg-ACC both these OH vibrations were prominent (bottom pattern in Fig. 5 and top spectra in inset), clearly showing that this phase was highly hydrated, also confirming the TGA data (~1.4 mol of water per mole of Mg-ACC). Upon crystallization to proto-dolomite and dolomite (upper pattern in main Fig. 5), the intensity of the broad OH-stretch and sharper OH-bend vibrations became weaker (inset in Fig. 5, normalized for comparison to the common 2525 cm−1 2ν2+ν4 symmetric C-O vibration). The prominent OH band in the Mg-ACC spectra decreased as proto-dolomite partially dehydrated (marked with * in inset Fig. 5), while dolomite was fully anhydrous (lower 3 patterns in inset Fig. 5).
The FEG-SEM images of the Mg-ACC revealed the presence of ~45 nm Ø spherical nanoparticles (Fig. 6b; Table 1) that formed loose aggregates. These transformed to proto-dolomite (Figs. 6c and 6d) that were characterized by small, nanocrystals that coalesced into larger spheroidal aggregates or co-joined micrometer-sized hemispheres (Fig. 6c). The sizes of the nanocrystalline subunits in these larger proto-dolomite spheroids were smaller than the precursor Mg-ACC, and averaged ~35 nm (Fig. 6d; Table 1). Finally, the higher temperature dolomite formed larger spheres (~5–50 μm, Fig. 6d), but these were made up of larger (>50 and <600 nm) and highly crystalline subunits (Fig. 6f) characterized by idiomorphic crystal faces. The size of the individual subunits in the proto-dolomite and dolomite evaluated from the FEG-SEM images are similar to the Scherrer crystallite sizes calculated from the PXRD data. These gradually increased from 23 nm for proto-dolomite formed at 60 °C to 139 nm for dolomite formed at 220 °C (Table 1).
The results described above show that starting from supersaturated aqueous solutions can lead to the formation of dolomite via a three-stage reaction. The first stage is the fast precipitation of a poorly ordered, hydrated, and nanoparticulate precursor, Mg-ACC. In stage two, this precursor breaks down to form proto-dolomite, which ultimately transforms to dolomite in stage three. Each of these three stages and a detailed mechanistic evaluation of the reaction sequence are described below.
Nature and stability of the Mg-ACC precursor
The formation of the Mg-ACC was instantaneous (within <1 s) in our experiments and thus its nucleation and growth were not evaluated in this study, we only characterized it and followed its crystallization. The PXRD and FTIR data for the Mg-ACC precursor phase are identical to those for ACC (Rodriguez-Blanco et al. 2008, 2011). The morphology (spherical) and size (44 ± 13 nm) of our starting Mg-ACC were also similar to that reported for ACC (e.g., Radha et al. 2010, sample ACC-5, Ø 50–100 nm; Bots et al. 2012, ACC, Ø ~ 40 nm). These suggest that our Mg-ACC is structurally similar to pure ACC. However, the ~1.4 mol of water in our Mg-ACC are higher than the average 1 mol of water for pure ACC in the literature (Koga et al. 1998; Koga and Yamane 2008; Michel et al. 2008; Radha et al. 2012). Our data clearly links the high-Mg content in our starting Mg-ACC with the higher water content of our Mg-rich precursor when compared to pure ACC. This is also is consistent with the recent study by Radha et al. (2012), who evaluated to crystallization energetics of various Mg containing ACC phases and inferred that a hydrated amorphous phase with a close to 0.5 Mg content is the least metastable and may be precursors for dolomite formation. This is linked with the fact that Mg is well known to have a significantly higher dehydration energy compared to calcium (Di Tommaso and de Leeuw 2010), and thus it is not surprising that it can retain more water when incorporated into the amorphous carbonate phase (Wang et al. 2009; Schmidt et al. 2005).
The decrease in background intensity in the time-resolved ED-XRD patterns to ~0.5 (Fig. 4) prior to the appearance of any diffraction peak of proto-dolomite suggests that the Mg-ACC started transforming before the first proto-dolomite crystals grew. Studies of both inorganic and biogenic ACC have shown that an initial hydrated and highly disordered ACC (i.e., ACC with low-structural order) forms from solution and that this transforms to a less disordered and anhydrous ACC prior to crystallization (Gower 2008; Bots et al. 2012). Radha et al. (2010) showed that the thermodynamic stability (enthalpy) of less-disordered and anhydrous ACC is greater than that of more-disordered and hydrated ACC, therefore the dehydration and ordering of the ACC is likely to be thermodynamically driven (lowering of enthalpy in the system). We have shown that the initially hydrated and highly disordered Mg-ACC forms from solution and we suggest that this Mg-ACC gradually orders and may partially dehydrate prior to its crystallization, and that this likely occurs via a similar mechanism to that observed for pure ACC (Bots et al. 2012). Less-disordered Mg-ACC has not been described previously but the dramatic decrease in background intensity in our ED-XRD patterns prior to its crystallization of proto-dolomite (Fig. 2) can be explained by the gradual formation of a more-ordered Mg-ACC. This is akin to the behavior observed for pure ACC (Radha et al. 2010; Bots et al. 2012) where prior to crystallization an increase in ordering of the amorphous phase was also likely accompanied by a lowering in water content and a final transformation to a fully anhydrous amorphous phase.
When compared to pure ACC, the current study showed that high concentrations of magnesium in the Mg-ACC dramatically delayed (300%) the onset of the crystallization. This is comparable to other studies where lower magnesium contents (10%) either delayed the crystallization onset (by up to 100%; Reddy and Nancollas 1976; Loste et al. 2003; Raz et al. 2003) or changed the crystallization pathway (Rodriguez-Blanco et al. 2012). For example, pure ACC, at ambient temperatures transforms within minutes to vaterite and subsequently to calcite (Rodriguez-Blanco et al. 2011), while 10% Mg in the ACC leads to a direct transformation to of the amorphous precursor to calcite with no vaterite intermediate but with a longer an induction time of ~15 min (Rodriguez-Blanco et al. 2012).
Comparing these low-Mg-ACC transformation induction times to data in the current study shows that the high-Mg content in our Mg-ACC (χMgCO3 = 0.39) significantly increases the induction time of crystallization, allowing the poorly order Mg-ACC phase to remain stable for days at room temperature, and more than 30 min at 60 °C (Fig. 3). This is consistent with studies that showed a direct proportionality between Mg content and crystallization induction time in the calcium carbonate system (e.g., calcite formation; Reddy and Nancollas 1976; Tracy et al. 1998a, 1998b; Loste et al. 2003; Raz et al. 2003). This stabilizing effect of Mg in ACC can be attributed to the presence of hydrated Mg within the framework of the Mg-ACC nanostructure, which retards proto-dolomite nucleation due to its slow dehydration kinetics (Raiteri and Gale 2010; Jiang et al. 2010). We previously described a similar effect also in a system where the Mg content in the initial aqueous solution was lower [Rodriguez-Blanco et al. (2012); 10% of the aqueous Ca2+ was replaced with Mg2+] compared to the current study (50% of the aqueous Ca2+ was replaced with Mg2+ in this study). In Rodriguez-Blanco et al. (2012) we demonstrated that the presence of Mg delayed the formation of the first crystalline phase compared to the pure system (Rodriguez-Blanco et al. 2011) and also that direct calcite formation from the amorphous precursor was favored.
Mechanism of formation for proto-dolomite
To derive the mechanism of crystallization for the precursor Mg-ACC our data offers various lines of evidence which include: the overall morphology, the crystallite size, the supersaturation, as well as the reaction rates all point to a spherulitic crystal growth mechanism for proto-dolomite. These will be now discussed in detail below.
First, the morphology of the proto-dolomite particles formed (i.e., spheroidal aggregates made up of nanocrystalline subunits) in all our experiments is indicative of a spherulitic growth mechanism (Andreassen 2005; Gránásy et al. 2005). Spherulitic growth of materials has been shown to occur in many systems, via a nucleation-controlled growth process where the continuous nucleation of new particles on the surface of existing particles occurs via non-crystallographic branching (i.e., there is no structural relationship between the nucleating particle and that of the particle on which it has formed; Shtukenberg et al. 2012). This “secondary” process is usually termed “growth front nucleation” (Gránásy et al. 2005), and results in the formation of spherulites a few micrometers in size consisting of aggregated and fairly monodispersed nanoparticulate crystallites. Gránásy et al. (2005) argues that small changes in the nucleation process can lead to spherulites with various morphologies. These can be spheres (classified as category 1) or dumbbells (category 2). The growth morphology of the proto-dolomite in our study (co-joined micrometer-sized hemispheres consisting of aggregated nanoparticles) is consistent with “category 2” spherulites. However, morphological evidence on its own is not satisfactory to provide definitive confirmation for a spherulitic growth mechanism. Our second line of evidence supporting this mechanism is the high supersaturation with respect to proto-dolomite throughout our reaction. Once the Mg-ACC had precipitated after mixing of the initial highly supersaturated solutions, the SI for proto-dolomite was between 2.46 and 3.00 (at 25 and 220 °C, respectively, Supplementary Table 11). Such a high supersaturation fosters a high crystallization driving force, which is another prerequisite for the continuous growth front nucleation process that controls spherulitic growth. Furthermore, the dissolution of minor amounts of Mg-ACC prior to the crystallization of proto-dolomite would increase the SI to higher levels, especially at higher temperatures (e.g., 220 °C) due to the reverse solubility of carbonates. As the reaction proceeds via a continual dissolution of the precursor Mg-ACC and the growth of proto-dolomite, the Mg-ACC is likely to be present until virtually all the proto-dolomite has formed, maintaining the high level of supersaturation throughout the spherulitic crystal growth reaction leading to proto-dolomite. Previous studies that followed the crystallization of calcium carbonates have shown that a saturation index >2–3 is required for spherulitic growth to occur. For example, in the case of the pure ACC to vaterite transformation that occurs also via spherulitic growth (Beck and Andreassen 2010; Andreassen et al. 2010; Bots et al. 2012) the SI for vaterite was usually equal or larger than 3.
The third line of evidence is the crystallization rate. Due to the nucleation-controlled growth mechanism, and the high level of supersaturation throughout the growth process, the formation rate during spherultic growth is usually extremely rapid, relative to conventional surface controlled growth. The time resolved ED-XRD data indicates that the entire crystallization reaction is very rapid. Regardless of the temperature (60–220 °C) the complete crystallization occurred in ~5 min, which is consistent with the growth front nucleation process occurring during spherulitic growth. Finally, the nucleation controlled growth mechanism during spherulitic growth leads to the formation of nanoparticles, as very little, or no, surface growth occurs following nucleation. The average crystallite/particle sizes determined for the proto-dolomite in our experiments of between 23 and 34 (<100 °C, Table 1) are consistent with a nucleation dominated particle formation reaction. These particle sizes are comparable with the crystallite size for vaterite also formed via spherulitic growth (~10 nm, Bots et al. 2012).
A kinetic model, which described the formation of three-dimensional spherulites via growth front nucleation, was derived by Gránásy et al. (2005). They indicate that the degree of reaction (α) for the crystalizing phase can be defined by the following equation
where k (s−1) is the rate constant, t (s) is time, t0 (s) is the induction time. This kinetic model produced the best fit to the data for the growth of proto-dolomite obtained from the time-resolved ED-XRD data (Fig. 3; Table 1), and was a significantly better fit than any kinetic models for surface-controlled growth or simple chemical kinetics. This fit to the kinetic model for spherultic growth provides additional evidence for this growth mechanism for proto-dolomite.
In summary, the overall morphology, the solution chemistry, crystallite sizes, as well as the derived kinetics for proto-dolomite formation all point to a spherulitic growth mechanism and no other crystallization mechanisms fits all the observations. This conclusion is consistent with previous studies of pure ACC crystallization in solution, which leads to the spherulitic growth of vaterite. However, this is the first study to show that proto-dolomite forms from Mg-ACC via the same spherulitic crystallization mechanism.
Proto-dolomite to dolomite transformation
The offline crystallization data indicated that following the spherulitic growth of proto-dolomite, a gradual transformation to dolomite occurred at temperatures >100 °C (Figs. 1 and 5). During this stage of the reaction no increase in PXRD peak intensities was observed (Figs. 2 and 4), suggesting no additional material formed. Malone et al. (1996) proposed that the transformation of proto-dolomite to ordered dolomite may involve multiple episodes of recrystallization and Ostwald-ripening. Kelleher and Redfern (2002) proposed Ostwald-ripening of a hydrous Ca-Mg carbonate phase (likely proto-dolomite) as the kinetic answer to the problem of dehydrating Mg to form ordered dolomite. This mechanism is consistent with the observation in this current study. The increase in crystallite and particle sizes (Table 1) and the development of well-defined crystal faces as proto-dolomite transforms to ordered dolomite, occurs via an Ostwald ripening, dissolution and reprecipitation mechanism. Ostwald ripening reactions occur when the supernatant solution is near saturation with respect to the crystallizing phase. This is in contrast to the very high level of saturation required for proto-dolomite spherulitic growth. Once all Mg-ACC crystallized to proto-dolomite via spherulitic growth, and the Ostwald ripening to ordered dolomite was initiated, the SI of the solution with respect to dolomite must have decreased dramatically due to its equilibration with the next most soluble phase, which is proto-dolomite. The small driving force required for the Ostwald ripening reaction is a consequence of the small difference in solubility between proto-dolomite and ordered dolomite. This second crystallization stage was slower than the spherulitic proto-dolomite growth (hours to weeks), but ultimately led to the formation of anhydrous and highly crystalline, euhedral dolomites with a stoichiometric Ca:Mg composition (Table 1). Malone et al. (1996) also crystallized proto-dolomite to ordered dolomite and showed that below 100 °C, even after 336 days only 50% of his initial proto-dolomite had recrystallized, while above 100 °C, the recrystallization rates were higher but only reached 100% at 200 °C after 40 days. In his recrystallization experiments, Malone et al. (1996) regardless of temperature and time, only obtained ordered dolomite that was not non-stoichiometric. The temperature-dependent proto-dolomite to ordered dolomite transformation was confirmed in the current study, yet we obtained highly ordered and stoichiometric dolomites and obtained a full crystallization at much faster rates (i.e., >100 °C within 6 days). This difference is most likely a consequence of the variations in supersaturation levels during this stage of the reaction, which were higher in the current study compared to those in the Malone et al. (1996) study, who re-equilibrated his synthesized proto-dolomite with a new solution of seawater ionic strength. Nevertheless, both studies show a proto-dolomite to dolomite transition, which is driven by a surface dissolution and reprecipitation mechanism as also suggested by Kelleher and Redfern (2002).
How does the formation of proto-dolomite/dolomite from Mg-ACC compare to the crystallization of pure ACC?
The fast precipitation of Mg-ACC followed by spherulitic growth of proto-dolomite and ripening to dolomite is similar to the mechanism described for the pure calcium carbonate system. In the absence of magnesium, ACC transforms to vaterite via spherulitic growth and this subsequently ripens and re-crystallizes to calcite (Bots et al. 2012; Rodriguez-Blanco et al. 2011). Our initial hypothesis that an amorphous precursor is a prerequisite for the formation of proto-dolomite via spherulitic growth is correct. We further assert that the first crystalline phase to form from all compositions of ACC (with or without Mg) always occurs via a spherulitic growth mechanism. The prime reason for this universal process is the large difference in solubility between the ACC (or Mg-ACC of variable compositions) and all crystalline Mg-Ca-carbonate polymorphs. This spherulitic growth pathway has so far been documented for the crystallization of pure ACC (Andreassen 2005; Bots et al. 2012; Rodriguez-Blanco et al. 2011; Sand et al. 2012), for a 30% Mg-ACC precursor to monohydrocalcite (Rodriguez-Blanco et al. 2014) as well as in this study where the initial solutions contained 50% Mg. A comprehensive evaluation of the formation of other Ca/MgCO3 phases from amorphous precursors needs to be undertaken to prove this hypothesis.
Do our experimental results mimic the observed “direct” microbial mediated proto-dolomite/dolomite formation in modern hypersaline settings?
In such environments biogenic proto-dolomite forms at ambient to moderately high temperatures (e.g., alkaline lakes or shabkhas T between 27–50 °C), high alkalinities (up to 500 mM), and high-Mg:Ca concentration ratios (>1:1 to 100:1) (Folk and Land 1975; Last and De Deckker 1990; Meister et al. 2011 and references therein). These conditions lead to supersaturated conditions suitable for the precipitation of amorphous carbonates. Comparing our experimental results with bacterially mediated dolomite formation studies (e.g., Kenward et al. 2009; Sánchez-Román et al. 2011; Meister et al. 2011), reveals many similarities including the resulting particle morphologies, the supersaturation/alkalinity and temperature. In our experiments, we showed that the reaction from the amorphous precursor to a first crystalline phase occurs in all cases via fast, secondary nucleation dominated crystallization. Such a reaction results in nanocrystalline particles that develop a spherulitic morphology, a shape that is akin to the “dumb-bells and cauliflower structures” described by Mckenzie and Vasconcelos (2009) for bacterially mediated proto-dolomites. Thus we suggest that such a spherulitic growth mechanism is the controlling prerequisite for the formation of the intricately shaped biomineralized proto-dolomite and dolomite observed in modern hypersaline settings. In our inorganic experiments the initial amorphous Mg-ACC precipitates from solutions that were highly supersaturated with respect to proto-dolomite (SI > 8), while the saturation indexes for proto-dolomite and dolomite dropped to close to 3 once Mg-ACC had formed. Similarly, in bacterially mediated reactions, SI values of 3.36 were reported (Kenward et al. 2009) which are comparable to those reported here and by others for spherulitic growth carbonate systems (Andreassen 2005; Bots et al. 2012). Furthermore, our Mg-ACC only crystallized to proto-dolomite at temperatures >40 °C, yet in microbial systems proto-dolomite and dolomite have been shown to form at temperatures as low as 25 °C. This indicates that besides the high supersaturation and high alkalinity, another so far not understood factor has to affect the dolomite crystallization process. This factor may be the presence of organics. Zhang et al. (2012), based on their experiments with carboxymethyl cellulose and agar, suggested that extracellular polysaccharides produced by microorganisms could be key in the formation of proto-dolomite at ambient temperatures. Yet, based on the above described chemical as well as the morphological similarities (spheroidal nanocrystalline aggregates), our data indicates that in microbially mediated systems the initial formation of proto-dolomite must proceed via the same spherulitic transformation pathway of the amorphous precursor leading to proto-dolomite and then to dolomite via re-crystallization. We suggest that in bacterial systems dolomite forms due to the alkalinity produced by the microbial processes, and that the degree of order and also the Mg content in the resulting mineral phase increases from proto-dolomite to dolomite during ripening, both processes we have quantified in this current study for the inorganic system.
Finally, what can we learn from our experimental results with regard to the inferred “direct” dolomite precipitation processes in ancient natural systems?
Hood et al. (2011) and Hood and Wallace (2012) suggest that the environmental conditions that dominated and promoted the “direct precipitation” of dolomite from solution in the Neoproterozoic seas were characterized by high alkalinities, elevated Mg:Ca ratios (>7) and high supersaturations. These conditions are similar to those in our experiments indicating that the dolomite cements most likely formed via the precipitation of an Mg-rich ACC, and its transformation to proto-dolomite via spherulitic growth as we have documented in this study. The Cryogenian dolomite cements were characterized by radially shaped growth morphologies, similar to the spherical morphology observed in this study and in the microbial studies. We suggest that the radial morphology is a remnant of a cross section through spherulites formed during the spherulitic growth of proto-dolomite. This leads us to infer that if the assertion of a “direct” dolomite cement formation scenario is correct, the spherulitic growth mechanism could have been the controlling pathway to produce the Cryogenian cap dolomites described by Hood et al. (2011), regardless of whether their formation occurred via a microbially mediated processes as suggested by Font et al. (2010).
Overall, our combined data sets showed that the occurrence, structure, composition, hydration, and morphology of proto-dolomite and dolomite can be explained by a three-step process: (1) the precipitation of a highly hydrated, magnesium deficient, and amorphous carbonate precursor, Mg-ACC; (2) the dehydration of this Mg-ACC precursor and concomitant growth of the proto-dolomite via a spherulitic mechanism; and (3) the re-crystallization of proto-dolomite via Ostwald ripening to form dolomite. Our results also demonstrate that in biomineralization reactions, growth morphologies can be a consequence of spherulitic growth that is in turn promoted by high supersaturation. Last, our results have important palaeo-environmental implications as they shed light on the past ocean chemical conditions necessary for the formation of ancient dolomites in the geological record and help explain the formation of proto-dolomite/dolomite in modern settings.
This research was supported by the Marie Curie EU-FP6 MIN-GRO Research and Training Network under contract MRTN-CT-2006-035488 and CCLRC beamtime award (Grant No. 50115 to Liane G. Benning). The authors would like to thank Dave Taylor from the Synchrotron Radiation Source, Daresbury Laboratory, U.K., for help during beamtime.