The role of background electrolytes on the kinetics and mechanism of calcite dissolution

The influence of background electrolytes on the mechanism and kinetics of calcite dissolution was investigated using in situ Atomic Force Microscopy (AFM). Experiments were carried out far from equilibrium by passing alkali halide salt (NaCl, NaF, NaI, KCl and LiCl) solutions over calcite cleavage surfaces. This AFM study shows that all the electrolytes tested enhance the calcite dissolution rate. The effect and its magnitude is determined by the nature and concentration of the electrolyte solution. Changes in morphology of dissolution etch pits and dissolution rates are interpreted in terms of modification in water structure dynamics (i.e. in the activation energy barrier of breaking water-water interactions), as well as solute and surface hydration induced by the presence of different ions in solution. At low ionic strength, stabilization of water hydration shells of calcium ions by non-paired electrolytes leads to a reduction in the calcite dissolution rate compared to pure water. At high ionic strength, salts with a common anion yield similar dissolution rates, increasing in the order Cl⁻ < I⁻ < F⁻ for salts with a common cation due to an increasing mobility of water around the calcium ion. Changes in etch pit morphology observed in the presence of F⁻ and Li⁺ are explained by stabilization of etch pit edges bonded by like-charged ions and ion incorporation, respectively. As previously reported and confirmed here for the case of F⁻, highly hydrated ions increased the etch pit nucleation density on calcite surfaces compared to pure water. This may be related to a reduction in the energy barrier for etch pit nucleation due to disruption of the surface hydration layer.     

Variation in calcite dissolution rates:: A fundamental problem?

A comparison of published calcite dissolution rates measured far from equilibrium at a pH of ∼ 6 and above shows well over an order of magnitude in variation. Recently published AFM step velocities extend this range further still. In an effort to understand the source of this variation, and to provide additional constraint from a new analytical approach, we have measured dissolution rates by vertical scanning interferometry. In areas of the calcite cleavage surface dominated by etch pits, our measured dissolution rate is 10^−10.95 mol/cm²/s (PCO₂ 10^−3.41 atm, pH 8.82), 5 to ∼100 times slower than published rates derived from bulk powder experiments, although similar to rates derived from AFM step velocities. On cleavage surfaces free of local etch pit development, dissolution is limited by a slow, “global” rate (10^−11.68 mol/cm²/s). Although these differences confirm the importance of etch pit (defect) distribution as a controlling mechanism in calcite dissolution, they also suggest that “bulk” calcite dissolution rates observed in powder experiments may derive substantial enhancement from grain boundaries having high step and kink density. We also observed significant rate inhibition by introduction of dissolved manganese. At 2.0 μM Mn, the rate diminished to 10^−12.4 mol/cm²/s, and the well formed rhombic etch pits that characterized dissolution in pure solution were absent. These results are in good agreement with the pattern of manganese inhibition in published AFM step velocities, assuming a step density on smooth terraces of ∼9 μm⁻¹.     

Solubility of nanocrystalline scorodite and amorphous ferric arsenate: Implications for stabilization of arsenic in mine wastes

Solubility experiments were performed on nanocrystalline scorodite and amorphous ferric arsenate. Nanocrystalline scorodite occurs as stubby prismatic crystals measuring about 50 nm and having a specific surface area of 39.88 ± 0.07 m²/g whereas ferric arsenate is amorphous and occurs as aggregated clusters measuring about 50–100 nm with a specific surface area of 17.95 ± 0.19 m²/g. Similar to its crystalline counterpart, nanocrystalline scorodite has a solubility of about 0.25 mg/L at around pH 3–4 but has increased solubilities at low and high pH (i.e. <2 and >6). Nanocrystalline scorodite dissolves incongruently at about pH > 2.5 whereas ferric arsenate dissolution is incongruent at all the pH ranges tested (pH 2–5). It appears that the solubility of scorodite is not influenced by particle size. The dissolution rate of nanocrystalline scorodite is 2.64 × 10⁻¹⁰ mol m⁻² s⁻¹ at pH 1 and 3.25 × 10⁻¹¹ mol m⁻² s⁻¹ at pH 2. These rates are 3–4 orders of magnitude slower than the oxidative dissolution of pyrite and 5 orders of magnitude slower than that of arsenopyrite. Ferric arsenate dissolution rates range from 6.14 × 10⁻⁹ mol m⁻² s⁻¹ at pH 2 to 1.66 × 10⁻⁹ mol m⁻² s⁻¹ at pH 5. Among the common As minerals, scorodite has the lowest solubility and dissolution rate. Whereas ferric arsenate is not a suitable compound for As control in mine effluents, nanocrystalline scorodite that can be easily precipitated at ambient pressure and temperature conditions would be satisfactory in meeting the regulatory guidelines at pH 3–4.     

Magnesium inhibition of calcite dissolution kinetics

We present evidence of inhibition of calcite dissolution by dissolved magnesium through direct observations of the (104) surface using atomic force microscopy (AFM) and vertical scanning interferometry (VSI). Far from equilibrium, the pattern of magnesium inhibition is dependent on solution composition and specific to surface step geometry. In CO₂-free solutions (pH 8.8), dissolved magnesium brings about little inhibition even at concentrations of 0.8 × 10⁻³ molal. At the same pH, magnesium concentrations of less than 0.05 × 10⁻³ molal in carbonate-buffered solutions generate significant inhibition, although no changes in surface and etch pit morphology are observed. As concentrations exceed magnesite saturation, the dissolution rate shows little additional decrease; however, selective pinning of step edges results in unique etch pit profiles, seen in both AFM and VSI datasets. Despite the decreases in step velocity, magnesium addition in carbonated solutions also appears to activate the surface by increasing the nucleation rate of new defects. These relationships suggest that the modest depression of the bulk rate measured by VSI reflects a balance between competing reaction mechanisms that simultaneously depress the rate through selective inhibition of step movement, but also enhance reactivity on terraces by lowering the energy barrier to new etch pit formation.     

Biotite weathering and nutrient uptake by ectomycorrhizal fungus, Suillus tomentosus, in liquid-culture experiments

Ectomycorrhiza-forming fungi (EMF) alter the nutrient-acquisition capabilities of vascular plants, and may play an important role in mineral weathering and the partitioning of products of weathering in soils under nutrient-limited conditions. In this study, we isolated the weathering function of Suillus tomentosus in liquid-cultures with biotite micas incubated at room temperature. We hypothesized that the fungus would accelerate weathering by hyphal attachment to biotite surfaces and transmission of nutrient cations via direct exchange into the fungal biomass. We combined a mass-balance approach with scanning electron microscopy (SEM) and atomic force microscopy (AFM) to estimate weathering rates and study dissolution features on biotite surfaces. Weathering of biotite flakes was about 2-3 orders of magnitude faster in shaken liquid-cultures with fungus compared to shaken controls without fungus, but with added inorganic acids. Adding fungus in nonshaken cultures caused a higher dissolution rate than in inorganic pH controls without fungus, but it was not significantly faster than organic pH controls without fungus. The K⁺, Mg²⁺ and Fe²⁺ from biotite were preferentially partitioned into fungal biomass in the shaken cultures, while in the nonshaken cultures, K⁺ and Mg²⁺ was lost from biomass and Fe²⁺ bioaccumulated much less. Fungal hyphae attached to biotite surfaces, but no significant surface changes were detected by SEM. When cultures were shaken, the AFM images of basal planes appeared to be rougher and had abundant dissolution channels, but such channel development was minor in nonshaken conditions. Even under shaken conditions the channels only accounted for only 1/100 of the total dissolution rate of 2.7 × 10⁻¹⁰ mol of biotite m⁻² s⁻¹. The results suggest that fungal weathering predominantly occurred not by attachment and direct transfer of nutrients via hyphae, but because of the acidification of the bulk liquid by organic acids, fungal respiration (CO₂), and complexation of cations which accelerated dissolution of biotite. Results further suggest that both carbohydrate source (abundant here) and a host with which nutrients are exchanged (missing here) may be required for EMF to exert an important weathering effect in soils. Unsaturated conditions and physical dispersal of nutrient-rich minerals in soils may also confer a benefit for hyphal growth and attachment, and promote the attachment-mediated weathering which has been observed elsewhere on soil mineral surfaces.     

Calcite dissolution: Effects of trace cations naturally present in Iceland spar calcites

In situ dissolution experiments on a set of pure, optical quality Iceland spar calcite samples from four different localities showed etch pit step retreat rates to be inversely proportional to total inherent trace cation composition. Atomic absorption spectroscopy (AAS) revealed Fe²⁺, Mg²⁺, Mn²⁺ and Sr²⁺ in amounts varying from a few to hundreds of ppm. We used a very simple experimental set-up, with an Atomic Force Microscope (AFM) fluid cell and a droplet of MilliQ water. As the calcite dissolved and approached equilibrium with the solution, trace cations were released, which were then present for interaction with the dissolving surface. We monitored continuous free-drift dissolution, in situ, on fresh cleavage surfaces for up to 40 min. Dissolution produced one-layer-deep, rhombic etch pits that continually expanded as we collected images. The rhombohedral symmetry of calcite defines two obtuse and two acute edges on the cleavage surface of etch pits and these, as expected from previous work, had different dissolution rates. Despite identical experimental conditions for all samples, we observed lower step retreat rates for both obtuse and acute edges on calcite characterised by relatively high trace cation composition. Increased cation concentration, particularly Mn, was also correlated with rounding of obtuse-obtuse corners, resulting in obtuse step retreat rates similar to those for acute sides. Physcial limitations of the AFM technique were taken into account when measuring step rate retreat and results were collected only from single-layer etch pits, which represent crystalline calcite with minimal defects. Dissolution rates presented here are thus lower than previous reports for studies of deep etch pits and where the physical limitations of imaging may not have been considered. In addition to molecular-level proof that divalent cations inherent at ppm levels in the calcite affect the dissolution process, these results show that pure, optical quality Iceland spar calcite should not be considered pure in the chemical sense. The results imply that dissolution rates determined for ideal systems with pure, synthetic or natural, materials may be considered as the boundary condition for dissolution in real systems in nature, where cations are always present both in the solution and in the initial solid.     

Does reactive surface area depend on grain size? Results from pH 3, 25 °C far-from-equilibrium flow-through dissolution experiments on anorthite and biotite

Laboratory determined mineral weathering rates need to be normalised to allow their extrapolation to natural systems. The principle normalisation terms used in the literature are mass, and geometric- and BET specific surface area (SSA). The purpose of this study was to determine how dissolution rates normalised to these terms vary with grain size. Different size fractions of anorthite and biotite ranging from 180-150 to 20-10 μm were dissolved in pH 3, HCl at 25 °C in flow through reactors under far from equilibrium conditions. Steady state dissolution rates after 5376 h (anorthite) and 4992 h (biotite) were calculated from Si concentrations and were normalised to initial- and final- mass and geometric-, geometric edge- (biotite), and BET SSA. For anorthite, rates normalised to initial- and final-BET SSA ranged from 0.33 to 2.77 × 10⁻¹⁰ mol_feldspar m⁻² s⁻¹, rates normalised to initial- and final-geometric SSA ranged from 5.74 to 8.88 × 10⁻¹⁰ mol_feldspar m⁻² s⁻¹ and rates normalised to initial- and final-mass ranged from 0.11 to 1.65 mol_feldspar g⁻¹ s⁻¹. For biotite, rates normalised to initial- and final-BET SSA ranged from 1.02 to 2.03 × 10⁻¹² mol_biotite m⁻² s⁻¹, rates normalised to initial- and final-geometric SSA ranged from 3.26 to 16.21 × 10⁻¹² mol_biotite m⁻² s⁻¹, rates normalised to initial- and final-geometric edge SSA ranged from 59.46 to 111.32 × 10⁻¹² mol_biotite m⁻² s⁻¹ and rates normalised to initial- and final-mass ranged from 0.81 to 6.93 × 10⁻¹² mol_biotite g⁻¹ s⁻¹. For all normalising terms rates varied significantly (p ? 0.05) with grain size. The normalising terms which gave least variation in dissolution rate between grain sizes for anorthite were initial BET SSA and initial- and final-geometric SSA. This is consistent with: (1) dissolution being dominated by the slower dissolving but area dominant non-etched surfaces of the grains and, (2) the walls of etch pits and other dissolution features being relatively unreactive. These steady state normalised dissolution rates are likely to be constant with time. Normalisation to final BET SSA did not give constant ratios across grain size due to a non-uniform distribution of dissolution features. After dissolution coarser grains had a greater density of dissolution features with BET-measurable but unreactive wall surface area than the finer grains. The normalising term which gave the least variation in dissolution rates between grain sizes for biotite was initial BET SSA. Initial- and final-geometric edge SSA and final BET SSA gave the next least varied rates. The basal surfaces dissolved sufficiently rapidly to influence bulk dissolution rate and prevent geometric edge SSA normalised dissolution rates showing the least variation. Simple modelling indicated that biotite grain edges dissolved 71-132 times faster than basal surfaces. In this experiment, initial BET SSA best integrated the different areas and reactivities of the edge and basal surfaces of biotite. Steady state dissolution rates are likely to vary with time as dissolution alters the ratio of edge to basal surface area. Therefore they would be more properly termed pseudo-steady state rates, only appearing constant because the time period over which they were measured (1512 h) was less than the time period over which they would change significantly.     

Interferometric study of the dolomite dissolution: a new conceptual model for mineral dissolution

The dissolution rate and mechanism of three different cleavage faces of a dolomite crystal from Navarra (near Pamplona), Spain, were studied in detail by vertical scanning interferometry techniques. A total of 37 different regions (each about 124 × 156 μm in size) on the three sample surfaces were monitored as a function of time during dissolution at 25°C and pH 3. Dissolution produced shallow etch pits with widths reaching 20 μm during 8 h of dissolution. Depth development as a function of time was remarkably similar for all etch pits on a given dolomite surface.On the basis of etch pit distribution and volume as a function of time, the calculated dissolution rate increases from near zero to 4 × 10⁻¹¹ mol cm⁻² s⁻¹ over 5 h. The time variation is different for each of the three cleavage surfaces studied. In addition, the absolute dissolution rates of different parts of the dolomite crystal surface can be computed by using a reference surface. The different surfaces yield an “average” rate of 1.08 × 10⁻¹¹ mol cm⁻² s⁻¹ with a standard deviation of 0.3 × 10⁻¹¹ mol cm⁻² s⁻¹ based on about 60 analyses. The mean absolute rate of the dolomite surface is about 10 times slower than the rate calculated from etch pit dissolution alone. On the other hand, earlier batch rate data that used BET surface areas yield rates that are at least 30 to 60 times faster than our directly measured mean dissolution rate for the same pH and temperature.A conceptual model for mineral dissolution has been inferred from the surface topography obtained by the interferometry investigations. In this model, mineral dissolution is not dominated by etch pit formation itself but rather by extensive dissolution stepwaves that originate at the outskirts of the etch pits. These stepwaves control the overall dissolution as well as the dependence on temperature and saturation state.     

Magnesium-isotope fractionation during low-Mg calcite precipitation in a limestone cave - Field study and experiments

The chemical and isotopic composition of speleothem calcite and particularly that of stalagmites and flowstones is increasingly exploited as an archive of past environmental change in continental settings. Despite intensive research, including modelling and novel approaches, speleothem data remain difficult to interpret. A possible way foreword is to apply a multi-proxy approach including non-conventional isotope systems. For the first time, we here present a complete analytical dataset of magnesium isotopes (δ²⁶Mg) from a monitored cave in NW Germany (Bunker Cave). The data set includes δ²⁶Mg values of loess-derived soil above the cave (−1.0 ± 0.5‰), soil water (−1.2 ± 0.5‰), the carbonate hostrock (−3.8 ± 0.5‰), dripwater in the cave (−1.8 ± 0.2‰), speleothem low-Mg calcite (stalactites, stalagmites; −4.3 ± 0.6‰), cave loam (−0.6 ± 0.1‰) and runoff water (−1.8 ± 0.1‰) in the cave, respectively. Magnesium-isotope fractionation processes during weathering and interaction between soil cover, hostrock and solute-bearing soil water are non-trivial and depend on a number of variables including solution residence times, dissolution rates, adsorption effects and potential neo-formation of solids in the regolith and the carbonate aquifer. Apparent Mg-isotope fractionation between dripwater and speleothem low-Mg calcite is about 1000lnα_Mg-cc-Mg(aq) = −2.4‰. A similar Mg-isotope fractionation (1000lnα_Mg-cc-Mg(aq) ≈ −2.1‰) is obtained by abiogenic precipitation experiments carried out at aqueous Mg/Ca ratios and temperatures close to cave conditions. Accordingly, ²⁶Mg discrimination during low-Mg calcite formation in caves is highly related to inorganic fractionation effects, which may comprise dehydration of Mg²⁺ prior to incorporation into calcite, surface entrapment of light isotopes and reaction kinetics. Relevance of kinetics is supported by a significant negative correlation of Mg-isotope fractionation with the precipitation rate for inorganic precipitation experiments.     

Trace/minor element:calcium ratios in cultured benthic foraminifera. Part II: Ontogenetic variation

Ontogenetic (developmental stage) measurements of Mg/Ca and Sr/Ca were made on the benthic foraminifer Bulimina aculeata, which were cultured under controlled physicochemical conditions of temperature, pH, alkalinity, salinity, and trace- and minor-element concentrations. We utilized two methods of ontogenetic sampling—whole specimens progressively increasing in length and laser microdissection of a single specimen with subsequent analysis of dissected portions. A novel high-resolution laser-microdissection (HRLM) method allowed for precise (10 μm) cuts of the foraminiferal tests (shells) along the geometrically complex sutures distinguishing individual chambers. This new microdissection method limited sample loss and cross-contamination between foraminiferal chambers. Little or no variation in D_Sr was observed at different foraminiferal developmental stages. Conversely, D_Mg was enriched during a mid-developmental stage of whole-specimen samples (150-225 μm D_Mg = 1.6 × 10⁻³) compared to earlier and later stages (<150 μm, >225 μm D_Mg = 8.3 × 10⁻⁴). Further analysis of HRLM ontogenetic samples showed a larger, age-dependent D_Mg signature variation. This increase in shell Mg/Ca may contribute substantially to the measured inter-individual variability in Mg/Ca temperature prediction for cultured B. aculeata. Due to relatively large Mg/Ca inter- and intra-individual variability, measuring similar-size foraminiferal samples may improve the precision of paleotemperature prediction. Additionally, partial dissolution of the highest ontogenetically Mg-enriched calcite (D_Mg = 1.3 × 10⁻²-1.6 × 10⁻²) may occur in undersaturated bottom-water environments or during reductive cleaning procedures. Thus, the calcite phases remaining after partial dissolution by either natural or laboratory cleaning processes may not accurately represent the calcification environment.     

In situ atomic force microscopy measurements of biotite basal plane reactivity in the presence of oxalic acid

We have used a direct imaging technique, in situ atomic force microscopy (AFM), to observe the dissolution of the basal biotite surface by oxalic acid over a range of temperatures close to ambient conditions, using a specially designed AFM liquid cell and non-invasive intermittent contact mode of operation. From the 3-dimensional nanometre-resolution data sets, we observe a process characterised by the slow formation of shallow etch pits in the (0 0 1) surface and fast growth of etch pits from the resulting steps, which represent proxies for the {h k 0} surface. Measurements of dissolution rates as a function of temperature allow a determination of an apparent activation energy (E_a,app) for the process, via mass-loss calculations from image analysis. We obtain a value of E_a,app = 49 ± 2 kJ mol⁻¹, which is consistent with separate calculations based on planar area etch pit growth, and measurements of etch pit perimeters, indicating that this value of E_a,app is representative of {h k 0} surface dissolution. The measurement of etch pit perimeters also enables an estimation of apparent activation energy as a function of step density indicating substantially higher apparent activation energy, up to E_a,app = 140 kJ mol⁻¹, on extrapolation towards a pristine surface with no defects. We suggest that this higher value of E_a,app represents the slow formation of etch pits into the (0 0 1) surface.     

Role of molecular oxygen in the dissolution of siderite and rhodochrosite

The dissolution of siderite (FeCO₃) and rhodochrosite (MnCO₃) under oxic and anoxic conditions is investigated at 298 K. The anoxic dissolution rate of siderite is 10^−8.65 mol m⁻² s⁻¹ for 5.5 < pH < 12 and increases as [H⁺]^0.75 for pH < 5.5. The pH dependence is consistent with parallel proton-promoted and water hydrolysis dissolution pathways. Atomic force microscopy (AFM) reveals a change in pit morphology from rhombohedral pits for pH > 4 to pits elongated at one vertex for pH < 4. Under oxic conditions the dissolution rate decreases to below the detection limit of 10⁻¹⁰ mol m⁻² s⁻¹ for 6.0 < pH < 10.3, and hillock precipitation preferential to steps is observed in concurrent AFM micrographs. X-ray photoelectron spectroscopy (XPS) and thermodynamic analysis identify the precipitate as ferrihydrite. At pH > 10.3, the oxic dissolution rate is as high as 10^−7.5 mol m⁻² s⁻¹, which is greater than under the corresponding anoxic conditions. A fast electron transfer reaction between solution O₂ or [Fe³⁺(OH)₄]⁻ species and surficial >Fe^II hydroxyl groups is hypothesized to explain the dissolution kinetics. AFM micrographs do not show precipitation under these conditions. Anoxic dissolution of rhodochrosite is physically observed as rhombohedral pit expansion for 3.7 < pH < 10.3 and is chemically explained by parallel proton- and water-promoted pathways. The dissolution rate law is 10^−4.93[H⁺] + 10^−8.45 mol m⁻² s⁻¹. For 5.8 < pH < 7.7 under oxic conditions, the AFM micrographs show a tabular precipitate growing by preferential expansion along the a-axis, though the macroscopic dissolution rate is apparently unaffected. For pH > 7.7 under oxic conditions, the dissolution rate decreases from 10^−8.45 to 10^−9.0 mol m⁻² s⁻¹. Flattened hillock precipitates grow across the entire surface without apparent morphological influence by the underlying rhodochrosite surface. XPS spectra and thermodynamic calculations implicate the precipitate as bixbyite for 5.8 < pH < 7.7 and MnOOH (possibly feitnkechtite) for pH >7.7.     

Microbial dissolution of calcite at T = 28 °C and ambient pCO2

This study used batch reactors to quantify the mechanisms and rates of calcite dissolution in the presence and absence of a single heterotrophic bacterial species (Burkholderia fungorum). Experiments were conducted at T = 28°C and ambient pCO₂ over time periods spanning either 21 or 35 days. Bacteria were supplied with minimal growth media containing either glucose or lactate as a C source, NH₄⁺ as an N source, and H₂PO₄⁻ as a P source. Combining stoichiometric equations for microbial growth with an equilibrium mass-balance model of the H₂O-CO₂-CaCO₃ system demonstrates that B. fungorum affected calcite dissolution by modifying pH and alkalinity during utilization of ionic N and C species. Uptake of NH₄⁺ decreased pH and alkalinity, whereas utilization of lactate, a negatively charged organic anion, increased pH and alkalinity. Calcite in biotic glucose-bearing reactors dissolved by simultaneous reaction with H₂CO₃ generated by dissolution of atmospheric CO₂ (H₂CO₃ + CaCO₃ → Ca²⁺ + 2HCO₃⁻) and H⁺ released during NH₄⁺ uptake (H⁺ + CaCO₃ → Ca²⁺ + HCO₃⁻). Reaction with H₂CO₃ and H⁺ supplied ∼45% and 55% of the total Ca²⁺ and ∼60% and 40% of the total HCO₃⁻, respectively. The net rate of microbial calcite dissolution in the presence of glucose and NH₄⁺ was ∼2-fold higher than that observed for abiotic control experiments where calcite dissolved only by reaction with H₂CO₃. In lactate bearing reactors, most H⁺ generated by NH₄⁺ uptake reacted with HCO₃⁻ produced by lactate oxidation to yield CO₂ and H₂O. Hence, calcite in biotic lactate-bearing reactors dissolved by reaction with H₂CO₃ at a net rate equivalent to that calculated for abiotic control experiments. This study suggests that conventional carbonate equilibria models can satisfactorily predict the bulk fluid chemistry resulting from microbe-calcite interactions, provided that the ionic forms and extent of utilization of N and C sources can be constrained. Because the solubility and dissolution rate of calcite inversely correlate with pH, heterotrophic microbial growth in the presence of nonionic organic matter and NH₄⁺ appears to have the greatest potential for enhancing calcite weathering relative to abiotic conditions.     

Connections between surface complexation and geometric models of mineral dissolution investigated for rhodochrosite

Mineral dissolution rates have been rationalized in the literature by surface complexation models (SCM) and morphological and geometric models (GM), and reconciliation of these conceptually different yet separately highly successful models is an important goal. In the current work, morphological alterations of the surface are observed in real time at the microscopic level by atomic force microscopy (AFM) while dissolution rates are simultaneously measured at the macroscopic level by utilizing the AFM fluid cell as a classic flow-through reactor. Rhodochrosite dissolution is studied from pH = 2 to 11 at 298 K, and quantitative agreement is found between the dissolution rates determined from microscopic and macroscopic observations. Application of a SCM model for the interpretation of the kinetic data indicates that the surface concentration of >CO₃H regulates dissolution for pH < 7 while the surface concentration of >MnOH₂⁺ regulates dissolution for pH > 7. A GM model explains well the microscopic observations, from which it is apparent that dissolution occurs at steps associated with anisotropic pit expansion. On the basis of the observations, we combine the SCM and GM models to propose a step-site surface complexation model (SSCM), in which the dissolution rates are quantitatively related to the surface chemical speciation of steps. The governing SSCM equation is as follows: R = χ_1/2(k_co + k_ca)[>CO₃H] + χ_1/2(k_mo + k_ma)[>MnOH₂⁺ ], where R is the dissolution rate (mol m⁻² s⁻¹), 2χ_1/2 is the fraction of surface sites located at steps, [>CO₃H] and [>MnOH₂⁺ ] are surface concentrations (mol m⁻²), and k_co, k_ca, k_mo, and k_ma are the respective dissolution rate coefficients (s⁻¹) for the >CO₃H and the >MnOH₂⁺ surface species on obtuse and acute steps. We find k_co = 2.7 s⁻¹, k_ca = 2.1 × 10⁻¹ s⁻¹, k_mo = 4.1 × 10⁻² s⁻¹, k_ma = 3.7 × 10⁻² s⁻¹, and χ_1/2 = 0.015 ± 0.005. The rate coefficients quantify the net result of complex surface step processes, including double-kink initiation and single-kink propagation. We propose that the SSCM model may have general applicability for dissolution far from equilibrium of flat mineral surfaces of ionic crystals, at least those that dissolve by step retreat.     

The role of sulfate groups in controlling CaCO3 polymorphism

The nucleation and growth of CaCO₃ phases from aqueous solutions with SO₄²⁻:CO₃²⁻ ratios from 0 to 1.62 and a pH of ∼10.9 were studied experimentally in batch reactors at 25 °C. The mineralogy, morphology and composition of the precipitates were characterized by X-ray diffraction, Fourier transform infrared spectroscopy, scanning electron microscopy and microanalyses. The solids recovered after short reaction times (5 min to 1 h) consisted of a mixture of calcite and vaterite, with a S content that linearly correlates with the SO₄²⁻:CO₃²⁻ ratio in the aqueous solution. The solvent-mediated transformation of vaterite to calcite subsequently occurred. After 24 h of equilibration, calcite was the only phase present in the precipitate for aqueous solutions with SO₄²⁻:CO₃²⁻ ? 1. For SO₄²⁻:CO₃²⁻ > 1, vaterite persisted as a major phase for a longer time (>250 h for SO₄²⁻:CO₃²⁻ = 1.62). To study the role of sulfate in stabilizing vaterite, we performed a molecular simulation of the substitution of sulfate for carbonate groups into the crystal structure of vaterite, aragonite and calcite. The results obtained show that the incorporation of small amounts (<3 mole%) of sulfate is energetically favorable in the vaterite structure, unfavorable in calcite and very unfavorable in aragonite. The computer modeling provided thermodynamic information, which, combined with kinetic arguments, allowed us to put forward a plausible explanation for the observed crystallization behavior.     

The major-ion composition of silurian seawater

One-hundred fluid inclusions in Silurian marine halite were analyzed in order to determine the major-ion composition of Silurian seawater. The samples analyzed were from three formations in the Late Silurian Michigan Basin, the A-1, A-2, and B Evaporites of the Salina Group, and one formation in the Early Silurian Canning Basin (Australia), the Mallowa Salt of the Carribuddy Group. The results indicate that the major-ion composition of Silurian seawater was not the same as present-day seawater. The Silurian ocean had lower concentrations of Mg²⁺, Na⁺, and SO₄²⁻, and much higher concentrations of Ca²⁺ relative to the ocean’s present-day composition. Furthermore, Silurian seawater had Ca²⁺ in excess of SO₄²⁻. Evaporation of Silurian seawater of the composition determined in this study produces KCl-type potash minerals that lack the MgSO₄-type late stage salts formed during the evaporation of present-day seawater. The relatively low Na⁺ concentrations in Silurian seawater support the hypothesis that oscillations in the major-ion composition of the oceans are primarily controlled by changes in the flux of mid-ocean ridge brine and riverine inputs and not global or basin-scale, seawater-driven dolomitization. The Mg²⁺/Ca²⁺ ratio of Silurian seawater was ∼1.4, and the K⁺/Ca²⁺ ratio was ∼0.3, both of which differ from the present-day counterparts of 5 and 1, respectively. Seawaters with Mg²⁺/Ca²⁺ <2 facilitate the precipitation of low-magnesian calcite (mol % Mg < 4) marine ooids and submarine carbonate cements whereas seawaters with Mg²⁺/Ca²⁺ >2 (e.g., modern seawater) facilitate the precipitation of aragonite and high-magnesian calcite. Therefore, the early Paleozoic calcite seas were likely due to the low Mg²⁺/Ca²⁺ ratio of seawater, not the pCO₂ of the Silurian atmosphere.     

Thickness and structure of the water film deposited from vapour on calcite surfaces

Synchrotron X-ray reflectivity (SXR) was used to measure the thickness of the water film that adsorbs on a {} cleavage surface of calcite (CaCO₃) in a sample chamber where relative humidity could be controlled within the range from <4% to 90%. Gases used to carry water vapour were initially either 100% N₂ or 100% CO₂. The product water film was remarkably constant in thickness at 15.5 Å (±1 Å) and independent of humidity. When N₂ was used as the carrier gas, this film displayed a gap in its electron density at between 0.6 and 2 Å distance from the calcite surface, depending on humidity. This implies that a change in the arrangement of water molecules occurs in direct proximity to the surface. This electron density discontinuity was measurably further from the calcite surface, at 3.4 Å, when CO₂ was used as the carrier gas. Except for this thin low density region proximate to the calcite surface, the density of the adsorbed water layer was 0.9 g cm⁻³, therefore suggesting a significant degree of ordering. Atomic force microscopy (AFM) images were completed in conjunction with the SXR measurements on similarly prepared calcite cleavage surfaces. AFM showed that terraces may be atomically flat over 1 μm or more. SXR corroborated this observation, with results showing that carefully cleaved surfaces had a starting root mean square (r.m.s.) roughness of ∼1.2 Å. Diffuse scatter measurements constrained the correlation lengths of these surfaces to be at least 18,000 Å. For comparison with the cleaved samples, a surface was also prepared by chemo-mechanical Syton polishing. This surface gave an r.m.s roughness by SXR that was an order of magnitude higher, equal to 12.1 Å. In this case, diffuse scatter resolved a correlation length of 950 Å, and revealed a fractal dimension that was higher than for the cleaved surface. On Syton polished samples, the water film determined by SXR was about twice as thick as for freshly cleaved surfaces, with a density of 1.0 g cm⁻³, equal to that of bulk water. However, surface roughness was too large to allow resolution of any gap in the electron density within the water layer proximate to the solid surface. Our AFM observations also confirm previous reports of calcite surface recrystallization. The electron density of the solid surface determined by SXR is indistinguishable from that of calcite, indicating that any material recrystallized within the adsorbed water film is compositionally indistinguishable from the calcite substrate.     

Defect distribution and dissolution morphologies on low-index surfaces of α-quartz

The dissolution of prismatic and rhombohedral quartz surfaces by KOH/H₂O solutions was investigated by atomic force microscopy. Rates of dissolution of different classes of surface features (e.g., steps, voids, and dislocation etch pits) were measured. The prismatic surface etched almost two orders of magnitude faster than the rhombohedral surface, mostly due to the difference in the number and the rate of dissolution of extended defects, such as dislocations. Because of the presence of imperfect twin boundaries, defect densities on the prismatic surface were estimated at 50-100 μm⁻², whereas the rhombohedral surface possessed only ∼0.5-1.0 μm⁻², mostly in the form of crystal voids. Crystal voids etched almost one order of magnitude faster on the prismatic surface than on the rhombohedral surface due to differences in the number and the density of steps formed by voids on the different surfaces. In the absence of extended defects, both surfaces underwent step-wise dissolution at similar rates. Average rates of step retreat were comparable on both surfaces (∼3-5 nm/h on the prismatic surface and ∼5-10 nm/h on the rhombohedral surface). Prolonged dissolution left the prismatic surface reshaped to a hill-and-valley morphology, whereas the rhombohedral surface dissolved to form coalescing arrays of oval-shaped etch pits.     

Holographic interferometry study of the dissolution and diffusion of gypsum in water

We have performed holographic interferometry measurements of the dissolution of the (0 1 0) plane of a cleaved gypsum single crystal in pure water. These experiments have provided the value of the dissolution rate constant k of gypsum in water and the value of the interdiffusion coefficient D of its aqueous species in water. D is 1.0 × 10⁻⁹ m² s⁻¹, a value close to the theoretical value generally used in dissolution studies. k is 4 × 10⁻⁵ mol m⁻² s⁻¹. It directly characterizes the microscopic transfer rate at the solid-liquid interface, and is not an averaged value deduced from quantities measured far from the surface as in macroscopic dissolution experiments. It is found to be two times lower than the value obtained from macroscopic experiments.