Effects of magnesium ions on near-equilibrium calcite dissolution: Step kinetics and morphology

Dissolution kinetics at the aqueous solution-calcite interface at 50 °C were investigated using in situ atomic force microscopy (AFM) to reveal the influence of magnesium concentration and solution saturation state on calcite dissolution kinetics and surface morphology. Under near-equilibrium conditions, dissolved Mg²⁺ displayed negligible inhibitory effects on calcite dissolution even at concentrations of . Upon the introduction of , the solution saturation state with respect to calcite, , acted as a “switch” for magnesium inhibition whereby no significant changes in step kinetics were observed at Ω_calcite<0.2, whereas a sudden inhibition from Mg²⁺ was activated at Ω_calcite?0.2. The presence of the Ω-switch in dissolution kinetics indicates the presence of critical undersaturation in accordance with thermodynamic principles. The etch pits formed in solutions with exhibited a unique distorted rhombic profile, different from those formed in Mg-free solutions and in de-ionized water. Such unique etch pit morphology may be associated with the anisotropy in net detachment rates of counter-propagating kink sites upon the addition of Mg²⁺.     

Sorption of phosphate onto calcite; results from batch experiments and surface complexation modeling

The adsorption of phosphate onto calcite was studied in a series of batch experiments. To avoid the precipitation of phosphate-containing minerals the experiments were conducted using a short reaction time (3 h) and low concentrations of phosphate (?50 μM). Sorption of phosphate on calcite was studied in 11 different calcite-equilibrated solutions that varied in pH, P_CO2, ionic strength and activity of Ca²⁺, and . Our results show strong sorption of phosphate onto calcite. The kinetics of phosphate sorption onto calcite are fast; adsorption is complete within 2-3 h while desorption is complete in less than 0.5 h. The reversibility of the sorption process indicates that phosphate is not incorporated into the calcite crystal lattice under our experimental conditions. Precipitation of phosphate-containing phases does not seem to take place in systems with ?50 μM total phosphate, in spite of a high degree of super-saturation with respect to hydroxyapatite (SI_HAP ? 7.83). The amount of phosphate adsorbed varied with the solution composition, in particular, adsorption increases as the activity decreases (at constant pH) and as pH increases (at constant activity). The primary effect of ionic strength on phosphate sorption onto calcite is its influence on the activity of the different aqueous phosphate species. The experimental results were modeled satisfactorily using the constant capacitance model with >CaPO₄Ca⁰ and either >CaHPO₄Ca⁺ or > as the adsorbed surface species. Generally the model captures the variation in phosphate adsorption onto calcite as a function of solution composition, though it was necessary to include two types of sorption sites (strong and weak) in the model to reproduce the convex shape of the sorption isotherms.     

Ammonium loss and nitrogen isotopic fractionation in biotite as a function of metamorphic grade in metapelites from western Maine, USA

Ammonium fixed in micas of metamorphic rocks is a sensitive indicator both of organic-inorganic interactions during diagenesis as well as of the devolatilization history and fluid/rock interaction during metamorphism. In this study, a collection of geochemically well-characterized biotite separates from a series of graphite-bearing Paleozoic greenschist- to upper amphibolite-facies metapelites, western Maine, USA, were analyzed for ammonium nitrogen () contents and isotopic composition (δ¹⁵N_NH4) using the HF-digestion distillation technique followed by the EA-IRMS technique. Biotite separates, sampled from 9 individual metamorphic zones, contain 3000 to 100 ppm with a wide range in δ¹⁵N from +1.6‰ to +9.1‰. Average contents in biotite show a distinct decrease from about 2750 ppm for the lowest metamorphic grade (∼500 °C) down to 218 ppm for the highest metamorphic grade (∼685 °C). Decreasing abundances in are inversely correlated in a linear fashion with increasing K⁺ in biotite as a function of metamorphic grade and are interpreted as a devolatilization effect. Despite expected increasing δ¹⁵N_NH4 values in biotite with nitrogen loss, a significant decrease from the Garnet Zones to the Staurolite Zones was found, followed by an increase to the Sillimanite Zones. This pattern for δ¹⁵N_NH4 values in biotite inversely correlates with Mg/(Mg + Fe) ratios in biotite and is discussed in the framework of isotopic fractionation due to different exchange processes between or , reflecting devolatilization history and redox conditions during metamorphism.     

First-principles estimates of equilibrium magnesium isotope fractionation in silicate, oxide, carbonate and hexaaquamagnesium(2+) crystals

Equilibrium mass-dependent magnesium isotope fractionation factors are estimated for a range of crystalline compounds including oxides, silicates, carbonates, and salts containing the solvation complex. Fractionation factors for the gas-phase species Mg and MgO are also presented. Fractionation factors are calculated with density functional perturbation theory (DFPT), using norm-conserving pseudopotentials. The results suggest that there will be substantial inter-mineral fractionation, particularly between tetrahedrally coordinated Mg²⁺ in spinel (MgAl₂O₄) and the more common octahedrally coordinated Mg²⁺-sites in silicate and carbonate minerals. Isotope fractionations calculated for Mg²⁺ in hexaaquamagnesium(2+) salts are in good agreement with previous fractionation models of based on large molecular clusters (Black et al., 2007), but show possibly more significant disagreement with a more recent study (Rustad et al., 2010). These models further suggest that solvated , in the form of , will have higher ²⁶Mg/²⁴Mg than coexisting magnesite and dolomite. Calculated fractionations are consistent with Mg-isotope fractionations observed in peridotite mineral separates and inorganic carbonate precipitates. Predicted large, temperature-sensitive spinel-silicate fractionations, in particular, may find use in determining equilibration temperatures of peridotites and other high-temperature rock types.     

An experimental study of Cobalt (II) complexation in Cl and H2S-bearing hydrothermal solutions

The speciation of cobalt (II) in Cl⁻ and H₂S-bearing solutions was investigated spectrophotometrically at temperatures of 200, 250, and 300 °C and a pressure of 100 bars, and by measuring the solubility of cobaltpentlandite at temperatures of 120-300 °C and variable pressures of H₂S. From the results of these experiments, it is evident that CoHS⁺ and predominate in the solutions except at 150 °C, for which the dominant chloride complex is CoCl₃⁻. The logarithms of the stability constant for CoHS⁺ show moderate variation with temperature, decreasing from 6.24 at 120 °C to 5.84 at 200 °C, and increasing to 6.52 at 300 °C. Formation constants for chloride species increase smoothly with temperature and at 300°C their logarithms reach 8.33 for , 6.44 for CoCl₃⁻, 4.94 to 5.36 for , and 2.42 for CoCl⁺. Calculations based on the composition of a model hydrothermal fluid (Ksp-Mu-Qz, KCl = 0.25 m, NaCl = 0.75 m, ΣS = 0.3 m) suggest that at temperatures ?200 °C, cobalt occurs dominantly as CoHS⁺, whereas at higher temperatures the dominant species is .     

Molybdenum, tungsten and manganese partitioning in the system pyrrhotite-Fe-S-O melt-rhyolite melt: Impact of sulfide segregation on arc magma evolution

Experiments were performed to determine the partitioning of molybdenum, tungsten and manganese among a rhyolitic melt (melt), pyrrhotite (po), and an immiscible Fe-S-O melt (Fe-S-O). Sulfide phases such as these may be isolated from a silicate melt along with other crystallizing phases during the evolution of arc magma, and partition coefficients are required to model the effect of this process on molybdenum and tungsten budgets.We developed an experimental design to take advantage of properties of the phases under study. Careful control of temperature allowed pyrrhotite and magnetite to be stable along with an Fe-S-O melt, and this phase assemblage allowed the composition of run-product pyrrhotite to be used to calculate both fS₂ and fO₂ for the experiments. At run temperature, (1042 ± 2 °C), a rhyolitic melt can be formed at low pressure, under nominally dry conditions, which removed the need for confining pressure as well as externally imposed fugacities. The silica-saturated melt allowed the charges to be contained in sealed evacuated silica tubes without danger of reaction, and with closed system behavior for molybdenum and tungsten.Experiments were run for durations up to 2000 min. Molybdenite (mb) and wolframite (wo) were added to the experiments as sources for molybdenum and tungsten, respectively. Manganese was added to the system as both a component of the starting rhyolitic pumice, and of Mn-bearing wolframite. Oxygen fugacity in these experiments was fixed at the Ni-NiO oxygen fugacity buffer. Sulfur fugacity was 10⁻¹ bar. Run products were analyzed by EPMA and LA-ICP-MS. Analysis of the run products yielded ( standard deviation of the mean): , , , and . The partition coefficients for manganese in this system are and .Simple Rayleigh fractionation modeling suggests that oxidized felsic melts produced through fractional crystallization may have lost as much as 14% of their initial molybdenum, but only 2% of their initial tungsten, through the removal of an Fe-S-O melt along with crystalline phases. Modeling consistent with conditions of oxygen and sulfur fugacity influenced by assimilation of sulfide (with low concentrations of molybdenum and tungsten) from, for example, sedimentary rock, results in evolved magmas significantly depleted in molybdenum, but only moderately depleted in tungsten. The molybdenum:tungsten ratio can vary by two orders of magnitude. These systematics may help to explain some of the variability in metal ratios of intrusion-related hydrothermal ore deposits.     

An integrated sulfur isotope model for Namibian shelf sediments

In this study the sulfur cycle in the organic-rich mud belt underlying the highly productive upwelling waters of the Namibian shelf is quantified using a 1D reaction-transport model. The model calculates vertical concentration and reaction rate profiles in the top 500 cm of sediment which are compared to a comprehensive dataset which includes carbon, sulfur, nitrogen and iron compounds as well as sulfate reduction (SR) rates and stable sulfur isotopes (³²S, ³⁴S). The sulfur dynamics in the well-mixed surface sediments are strongly influenced by the activity of the large sulfur bacteria Thiomargaritanamibiensis which oxidize sulfide (H₂S) to sulfate () using sea water nitrate () as the terminal electron acceptor. Microbial sulfide oxidation (SOx) is highly efficient, and the model predicts intense cycling between and H₂S driven by coupled SR and SOx at rates exceeding 6.0 mol S m⁻² y⁻¹. More than 96% of the SR is supported by SOx, and only 2-3% of the pool diffuses directly into the sediment from the sea water. A fraction of the produced by Thiomargarita is drawn down deeper into the sediment where it is used to oxidize methane anaerobically, thus preventing high methane concentrations close to the sediment surface. Only a small fraction of total H₂S production is trapped as sedimentary sulfide, mainly pyrite (FeS₂) and organic sulfur (S_org) (∼0.3 wt.%), with a sulfur burial efficiency which is amongst the lowest values reported for marine sediments (<1%). Yet, despite intense SR, FeS₂ and S_org show an isotope composition of ∼5 ‰ at 500 cm depth. These heavy values were simulated by assuming that a fraction of the solid phase sulfur exchanges isotopes with the dissolved sulfide pool. An enrichment in H₂S of ³⁴S towards the sediment-water interface suggests that Thiomargarita preferentially remove H₂³²S from the pore water. A fractionation of 20-30‰ was estimated for SOx (ε_SOx) with the model, along with a maximum fractionation for SR (ε_SR-max) of 100‰. These values are far higher than previous laboratory-based estimates for these processes. Mass balance calculations indicate negligible disproportionation of autochthonous elemental sulfur; an explanation routinely cited in the literature to account for the large fractionations in SR. Instead, the model indicates that repeated multi-stepped sulfide oxidation and intracellular disproportionation by Thiomargarita could, in principle, allow the measured isotope data to be simulated using much lower fractionations for ε_SOx (5‰) and ε_SR (78‰).     

Chemical controls on the solubility, speciation and mobility of lanthanum at near surface conditions: A geochemical modeling study

Recent experimental determinations of the solubility products of common rare earth minerals such as monazite and xenotime and stability constants for chloride, sulfate, carbonate and hydroxide complexes provide a basis to model quantitatively the solubility, and therefore the mobility, of rare earth elements (REE) at near surface conditions. Data on the mobility of REE and stabilities of REE complexes at near-neutral conditions are of importance to safe nuclear waste disposal, and environmental monitoring. The aim of this study is to understand REE speciation and solubility of a given REE in natural environments. In this study, a series of formation constants for La aqueous complexes are recommended by using the specific interaction theory (SIT) for extrapolation to infinite dilution. Then, a thermodynamic model has been employed for calculation of the solubility and speciation of La in soil solutions reacted with the La end-member of mineral monazite (LaPO₄), and other La-bearing solid phases including amorphous lanthanum hydroxide (La(OH)₃, am) and different La carbonates, as a function of various inorganic and organic ligand concentrations. Calculations were carried out at near-neutral pH (pH 5.5–8.5) and 25 °C at atmospheric CO₂ partial pressure. The model takes account of the species: La³⁺, LaCl²⁺, , , , , , , , , La(OH)²⁺, LaOx⁺, , LaAc²⁺ and (where Ox²⁻ = oxalate and Ac⁻ = acetate).The calculations indicate that the La species that dominate at pH 5.5–8.5 in the baseline model soil solution (BMSS) include La³⁺, LaOx⁺, , and in order of increasing importance as pH rises. The solubility of monazite in the BMSS remains less than 3 × 10⁻⁹ M, exhibiting a minimum of 2 × 10⁻¹² M at pH 7.5. The calculations quantitatively demonstrate that the concentrations of La controlled by the solubility of other La-bearing solid phases are many orders of magnitude higher than those controlled by monazite in the pH range from 5.5 to 8.5, suggesting that monazite is likely to be the solubility-controlling phase at this pH range. The calculations also suggest that significant mobility of La (and other REE) is unlikely because high water–rock ratios on the order of at least 10⁴ (mass ratio) are required to move 50% of the La from a soil. An increase in concentration of oxalate by one order of magnitude from that of the baseline model solution results in the dominance of LaOx⁺ at pH 5.5–7.5. Similarly, the increase in concentration of by one order of magnitude makes the dominant species at pH 5.5–7.5. Above pH 7.5, carbonate complexes are important. The increase in oxalate or concentrations by one order of magnitude can enhance the solubility of monazite by a factor of up to about 6 below neutral pH, in comparison with that in the baseline model soil solution. From pH 7.0 to 8.5, the solubility of monazite in the soil solutions with higher concentrations of oxalate or is similar, or almost identical, to that in the BMSS.     

Hydroxamate siderophore-promoted reactions between iron(II) and nitroaromatic groundwater contaminants

Recent studies show that ferrous iron (Fe^II), which is often abundant in anaerobic soil and groundwater, is capable of abiotically reducing many subsurface contaminants. However, studies also demonstrate that Fe^II redox reactivity in geochemical systems is heavily dependent upon metal speciation. This contribution examines the influence of hydroxamate ligands, including the trihydroxamate siderophore desferrioxamine B (DFOB), on Fe^II reactions with nitroaromatic groundwater contaminants (NACs). Experimental results demonstrate that ring-substituted NACs are reduced to the corresponding aniline products in aqueous solutions containing Fe^II complexes with DFOB and two monohydroxamate ligands (acetohydroxamic acid and salicylhydroxamic acid). Reaction rates are heavily dependent upon solution conditions and the identities of both the Fe^II-complexing hydroxamate ligand and the target NAC. Trends in the observed pseudo-first-order rate constants for reduction of 4-chloronitrobenzene (k_obs, s⁻¹) are quantitatively linked to the formation of Fe^II species with standard one-electron reduction potentials, (Fe^III/Fe^II), below −0.3 V. Linear free energy relationships correlate reaction rates with the (Fe^III/Fe^II) values of different electron-donating Fe^II complexes and with the apparent one-electron reduction potentials of different electron-accepting NACs, (ArNO₂). Experiments describing a redox auto-decomposition mechanism for Fe^II-DFOB complexes that occurs at neutral pH and has implications for the stability of hydroxamate siderophores in anaerobic environments are also presented. Results from this study indicate that hydroxamates and other Fe^III-stabilizing organic ligands can form highly redox-active Fe^II complexes that may contribute to the natural attenuation and remediation of subsurface contaminants.     

The influence of temperature and seawater composition on calcite crystal growth mechanisms and kinetics: Implications for Mg incorporation in calcite lattice

The composition of carbonate minerals formed in past and present oceans is assumed to be significantly controlled by temperature and seawater composition. To determine if and how temperature is kinetically responsible for the amount of Mg incorporated in calcite, we quantified the influence of temperature and specific dissolved components on the complex mechanism of calcite precipitation in seawater. A kinetic study was carried out in artificial seawater and NaCl-CaCl₂ solutions, each having a total ionic strength of 0.7 M. The constant addition technique was used to maintain [Ca²⁺] at 10.5 mmol kg⁻¹ while [] was varied to isolate the role of this variable on the precipitation rate of calcite.Our results show that the overall reaction of calcite precipitation in both seawater and NaCl-CaCl₂ solutions is dominated by the following reaction:

     

Modeling sulfur isotope fractionation and differential diffusion during sulfate reduction in sediments of the Cariaco Basin

Sulfur isotope composition (δ³⁴S) profiles in sediment pore waters often show an offset between sulfate and sulfide much greater in magnitude than S isotope fractionations observed in pure cultures. A number of workers have invoked an additional reaction, microbial disproportionation of sulfur intermediates, to explain the offset between experimental and natural systems. Here, we present an alternative explanation based on modeling of pore water sulfate and sulfide concentrations and stable isotope data from the Cariaco Basin (ODP Leg 165, Site 1002B). The use of unique diffusion coefficients for and , based on their unequal molecular masses, resulted in an increase in the computed fractionation by almost 10‰, when compared to the common assumption of equal diffusion coefficients for the two species. These small differences in diffusion coefficients yield calculated isotopic offsets between coeval sediment pore water sulfate and sulfide without disproportionation (up to 53.4‰) that exceed the largest fractionations observed in experimental cultures. Furthermore, the diffusion of sulfide within sediment pore waters leads to values that are even greater than those predicted by our model for sulfate reduction with unique diffusion coefficients. These diffusive effects on the sulfur isotope composition of pore water sulfate and sulfide can impact our interpretations of geologic records of sulfate and sulfide minerals, and should be considered in future studies.     

Experimental study of the incorporation of Li, Sc, Al and other trace elements into olivine

We performed a series of synthesis experiments at 1 atm pressure to investigate the substitution mechanisms of 1+ and 3+ ions into olivine. Forsterite crystals were grown from bulk compositions that contained the element of interest (e.g. Li) and different amounts of additional single trace elements. By working at constant (major element) liquid composition and temperature we eliminated all compositional effects other than those due to the trace elements. Mineral-melt pairs were then analysed to determine the compositional-dependence of the partition coefficient (D), which corresponds to , and where [element] refers to weight concentration of the element in the respective phase.We find that Li forms a stable coupled substitution with Sc and, at above ∼500 ppm Sc in the crystal, Li⁺ and Sc³⁺ ions form an ordered neutral complex ([LiSc]). This complex dissociates at lower trace element concentrations and a second, concentration-independent, mechanism begins to dominate. This second solution mechanism is most likely 2Li⁺ ⇔ Mg²⁺ where one of the Li atoms is in an interstitial position in the crystal lattice. Natural olivines show Li contents slightly greater than Sc (on an atomic basis), indicating that both substitution mechanisms are significant. Unlike Sc, Al does not appear to form a stable complex with Li in the olivine structure.Sodium is highly incompatible in olivine with of ∼0.00015-0.03. Olivine-liquid partitioning of Na⁺ is independent of Sc³⁺ or Al³⁺ concentration. This indicates that the coupled substitution of Na⁺ with any 3+ ions is unlikely. Instead, the relevant substitution mechanism appears to be 2Na⁺ ⇔ Mg²⁺. Although independent of 3+ ion concentration, is inversely correlated with the Li concentration of both melts and crystals, implying that Na competes (unsuccessfully) with Li to replace Mg in the olivine structure.Aluminium is highly incompatible in forsterite . Values of are similar for all phase pairs synthesised from starting materials containing between 10 and 100,000 ppm Al. This suggests that Al is principally incorporated in forsterite by replacing one Mg and one Si atom , where the Al atoms on octahedral (Mg) and tetrahedral (Si) sites are dissociated from one another.The incorporation of gallium into forsterite is influenced by the presence of Li. Where Li concentration in the crystal is much greater than that of Ga (on an atomic basis) we find an excellent correlation between and melt Li content. This relationship indicates that Ga³⁺ and Li⁺ replace 2Mg²⁺ on octahedral sites and that the Ga and Li atoms are, like Sc and Li, strongly associated in the crystal structure.The mechanism by which scandium is incorporated into forsterite is strongly governed by the presence Li. As discussed above, ordered complexes form readily in forsterite in Li-rich experiments. Under Li-absent but Sc-rich conditions (Sc in the crystal >∼500 ppm), is proportional to the concentration of Sc in the melt. This indicates that Sc incorporation is charge-balanced by the formation of magnesium vacancies , and that both species are associated . At lower Sc concentrations (<500 ppm in the crystal), the concentration-dependence of partitioning indicates that the complexes dissociate.Our results demonstrate that partitioning of 1+ and 3+ ions into olivine is complex and involves a range of point defects which yield strongly composition-dependent crystal-melt partition coefficients. Since physical and chemical properties of natural olivine, such as diffusion of ⁶Li and ⁷Li and H₂O solubility, depend on the concentrations of the defects identified in this study, our results provide an important insight into how determining substitution mechanisms can improve our understanding of large-scale mantle processes and properties.     

Complexation of Cu in Hydrothermal NaCl Brines: Ab initio molecular dynamics and energetics

Chloride complexation of Cu⁺ controls the solubility of copper(I) oxide and sulfide ore minerals in hydrothermal and diagenetic fluids. Solubility measurements and optical spectra of high temperature CuCl solutions have been interpreted as indicating the formation of CuCl, , and complexes. However, no other monovalent cation forms tri- and tetrachloro complexes. EXAFS spectra of high temperature Cu-Cl solutions, moreover, appear to show only CuCl and complexes at T > 100 °C. To reconcile these results, I investigated the nature and stability of Cu-Cl complexes using ab initio cluster calculations and ab initio (Car-Parrinello) molecular dynamics simulations for CuCl-NaCl-H₂O systems at 25 to 450 °C. Ab initio molecular dynamic simulations of 1 m CuCl in a 4 m Cl solution give a stable complex at 25 °C over 4 ps but show that the third Cl⁻ is weakly bound. When the temperature is increased along the liquid-vapour saturation curve to 125 °C, the complex dissociates into and Cl⁻; only forms at 325 °C and 1 kbar. Even in a 15.6 m Cl brine at 450 °C, only the complex forms over a 4 ps simulation run.Cluster calculations with a static dielectric continuum solvation field (COSMO) were used in an attempt directly estimate free energies of complex formation in aqueous solution. Consistent with the MD simulations, the complex is slightly stable at 25 °C but decreases in stability with decreasing dielectric constant (ε). The complex is predicted to be unstable at 25 °C and becomes increasingly unstable with decreasing dielectric constant. In hydrothermal fluids (ε < 30) both the and complexes are unstable to dissociation into and Cl⁻.The results obtained here are at odds with recent equations of state that predict and complexes are the predominant species in hydrothermal brines. In contrast, I predict that only complexes will be significant at T > 125 °C, even in NaCl-saturated brines. The high-temperature (T > 125 °C) optical spectra of CuCl solutions and solubility measurements of Cu minerals in Cl-brines need to be reinterpreted in terms of only the CuCl and complexes.     

Molecular simulation of the diffusion of uranyl carbonate species in aqueous solution

Potential-based molecular dynamics simulations of aqueous uranyl carbonate species (M_xUO₂(CO₃)_y^2+2x−2y with M = Mg, Ca, or Sr) were carried out to gain molecular-level insight into the hydration properties of these species. The simulation results were used to estimate the self-diffusion coefficients of these uranyl carbonate species, which often dominate uranyl speciation in groundwater systems. The diffusion coefficients obtained for the monoatomic alkaline-earth cations and polyatomic ions (uranyl, carbonate, and uranyl tri-carbonate) were compared with those calculated from the Stokes-Einstein (SE) equation and its variant formulation by Impey et al. (1983). Our results show that the equation of Impey et al. (1983), originally formulated for monovalent monoatomic ions, can be extended to divalent monoatomic ions, with some success in reproducing the absolute values and the overall trend determined from the molecular dynamics simulations, but not to polyatomic ions, for which the hydration shell is not spherically symmetrical. Despite the quantitative failure of both SE formulations, a plot of the diffusion coefficients of the uranyl carbonate complexes as a function of the inverse of the equivalent spherical radius showed that a general linear dependence is observed for these complexes as expected from the SE equation. The nature of the alkaline-earth cation in the uranyl carbonate complexes was not found to have a significant effect on the ion’s diffusion coefficient, which suggests that the use of a single diffusion coefficient for different alkaline-earth uranyl carbonate complexes in microscopic diffusion models is appropriate.The potential model reproduced well published quantum mechanical and experimental data of and of the individual constituent ions, and therefore is expected to offer reliable predictions of the structure of magnesium and strontium uranyl carbonate aqueous species, for which there is no structural data available to date. In addition, the interatomic distances reported for could help with the refinement of the interpretation of EXAFS data of these species, which is made difficult by the similar uranium-distant carbonate oxygen and uranium-calcium distances.An analysis of the dynamics of water exchange around the alkaline-earth cations revealed that the presence of the uranyl tri-carbonate molecule has a strong influence on the geometry of the cation’s first hydration shell, which, in turn, can considerably affect the water exchange kinetics depending on whether the imposed geometry matches that around the isolated alkaline-earth cation. This result shows that the alkaline-earth uranyl carbonate complexes have distinct water exchange dynamics, which may lead to different reactivities. Finally, significant changes in water residence time were also predicted when replacing carbonate for water ligands in the uranyl coordination shell.