Modelling the effects of salt solutions on the hydration of calcium ions

Classical molecular dynamics simulations of several aqueous alkali halide salt solutions have been used to determine the effect of electrolytes on the structure of water and the hydration properties of calcium ions. Compared with the simulations of Ca2+ ions in pure liquid water, the frequency of water exchange in the first hydration shell of calcium, which is a fundamental process in controlling the reactivity of calcium(II) aqua-ions, is drastically reduced in the presence of other electrolytes in solution. The strong stabilization of the hydration shell of Ca2+ occurs not only when the halide anions are directly coordinated to calcium, but also when the alkali and halide ions are placed at or outside the second coordination shell of Ca2+, suggesting that the reactivity of the first solvation shell of the calcium ion can be influenced by the specific affinity of other ions in solution for the water molecules coordinated to Ca2+. Analysis of the hydrogen-bonded structure of water in the vicinity of the calcium ion shows that the average number of hydrogen bonds per water molecules, which is 1.8 in pure liquid water, decreases as the concentration of alkali–halide salts in solution increases, and that the temporal fluctuations of hydrogen bonds are significantly larger than those obtained for Ca2+ in pure liquid water. This effect has been explained in terms of the dynamics of reorganization of the O–H X (X = F, Cl and Br) hydrogen bond. This work shows the importance of solution composition in determining the hydrogen-bonding network and ligand-exchange dynamics around metal ions, both in solution and at the mineral–water interfaces, which in turn has implications for interactions occurring at the mineral–water interface, ultimately controlling the mobilization of ions in the environment as well as in industrial processes

Relative rates of fluid advection, elemental diffusion and replacement govern reaction front patterns

Replacement reactions during fluid infiltration into porous media, rocks and buildings are known to have important implications for reservoir development, ore formation as well as weathering. Natural observations and experiments have shown that in such systems the shape of reaction fronts can vary significantly ranging from smooth, rough to highly irregular. It remains unclear what process-related knowledge can be derived from these reaction front patterns. In this contribution we show a numerical approach to test the effect of relative rates of advection, diffusion, and reaction on the development of reaction fronts patterns in granular aggregates with permeable grain boundaries. The numerical model takes (i) fluid infiltration along permeable grain boundaries, (ii) reactions and (iii) elemental diffusion into account. We monitor the change in element concentration within the fluid, while reactions occur at a pre-defined rate as a function of the local fluid concentration. In non-dimensional phase space using Péclet and Damköhler numbers, results show that there are no rough fronts without advection (Péclet10−3). As advection becomes more dominant and reaction slower, roughness develops across several grains with a full microstructure mimicking replacement in the most extreme cases. The reaction front patterns show an increase in roughness with increasing Péclet number from Péclet 10 to 100 but then a decrease in roughness towards higher Péclet numbers controlled by the Damköhler number. Our results indicate that reaction rates are crucial for pattern formation and that the shape of reaction fronts is only partly due to the underlying transport mechanism

Reaction mechanism for the replacement of calcite by dolomite and siderite: Implications for geochemistry, microstructure and porosity evolution during hydrothermal mineralisation

Carbonate reactions are common in mineral deposits due to CO2-rich mineralising fluids. This study presents the first in-depth, integrated analysis of microstructure and microchemistry of fluid-mediated carbonate reaction textures at hydrothermal conditions. In doing so, we describe the mechanisms by which carbonate phases replace one another, and the implications for the evolution of geochemistry, rock microstructures and porosity. The sample from the 1.95 Moz Junction gold deposit, Western Australia, contains calcite derived from carbonation of a metamorphic amphibole—plagioclase assemblage that has further altered to siderite and dolomite. The calcite is porous and contains iron-rich calcite blebs interpreted to have resulted from fluid-mediated replacement of compositionally heterogeneous amphiboles. The siderite is polycrystalline but nucleates topotactically on the calcite. As a result, the boundaries between adjacent grains are low-angle boundaries (<10°), which are geometrically similar to those formed by crystal–plastic deformation and recovery. Growth zoning within individual siderite grains shows that the low-angle boundaries are growth features and not due to deformation. Low-angle boundaries develop due to the propagation of defects at grain faces and zone boundaries and by impingement of grains that nucleated with small misorientations relative to each other during grain growth.The cores of siderite grains are aligned with the twin planes in the parent calcite crystal showing that the reactant Fe entered the crystal along the twin boundaries. Dolomite grains, many of which appear to in-fill space generated by the siderite replacement, also show alignment of cores along the calcite twin planes, suggesting that they did not grow into space but replaced the calcite. Where dolomite is seen directly replacing calcite, it nucleates on the Fe-rich calcite due to the increased compatibility of the Fe-bearing calcite lattice relative to the pure calcite. Both reactions are interpreted as fluid-mediated replacement reactions which use the crystallography and elemental chemistry of the calcite. Experiments of fluid-mediated replacement reactions show that they proceed much faster than diffusion-based reactions. This is important when considering the rates of reactions relative to fluid flow in mineralising systems

Reaction induced fracturing during replacement processes

Replacement processes are common transformation mechanisms in minerals and rocks at a variety of conditions and scales. The underlying mechanisms are, in general, poorly understood, but both mechanical and chemical processes are thought to be important. Replacement of leucite (KAlSi2O6) by analcime (NaAlSi2 O6·H2O) is common in silica-poor igneous rocks. A 10% increase in volume is associated with the replacement process, and this generates stresses that eventually cause fracturing of the reacting leucite. Experimentally reacted leucite samples display characteristic fracturing patterns that include both spalling of concentric 'onion-skin'-like layers near the reacting interface and the formation of cross-cutting, often hierarchically arranged, sets of fractures that divide the remaining leucite into progressively smaller domains. These structures may explain the 'patchy' alteration patterns observed in natural leucite samples and similar, so-called, mesh-textures associated with the serpentinization of olivine grains during hydration of mafic or ultramafic rocks. They are also strikingly similar to larger scale patterns formed during spheroidal weathering processes. A simple discrete element model illustrates the mechanics that control the formation of such systems, and shows how these replacement processes may be accelerated due to the generation of new reactive surface area by hierarchical fracturing processes. © Springer-Verlag 2008

Zircon coronas around Fe-Ti oxides: A physical reference frame for metamorphic and metasomatic reactions

Ilmenite in coronitic gabbros from the Bamble and Kongsberg sectors, southern Norway, is surrounded by zircons ranging in diameters from a fraction of a micrometer to 10 µm across. The zircons are inert during subsequent metamorphism (amphibolite- to pumpellyite-prehnite facies) and metasomatism (scapolitization and albitization) and can be found as trails in silicates (phlogopite, talc, chlorite, amphibole, albite, and tourmaline) in the altered rocks. The trails link up to form polygons outlining the former oxide grain boundary. This 3-dimensional framework of zircons is used to (a) recognize metasomatic origin of rocks, (b) quantify the mobility of elements during mineral replacement, (c) establish the growth direction of reaction fronts and to identify the reaction mechanism as dissolution-reprecipitation. Zircon coronas on Fe-Ti oxides have been described from a number of terrains and appear to be common in mafic rocks (gabbros and granulites) providing a tool for a better understanding of metasomatic and metamorphic reactions. © Springer-Verlag 2008

Alteration of crystalline zircon solid solutions: a case study on zircon from an alkaline pegmatite from Zomba-Malosa, Malawi

A natural, altered zircon crystal from an alkaline pegmatite from the Zomba-Malosa Complex of the Chilwa Alkaline Province in Malawi has been studied by a wide range of analytical techniques to understand the alteration process. The investigated zircon shows two texturally and chemically different domains. Whereas the central parts of the grain (zircon I) appear homogeneous in backscattered electron images and are characterised by high concentrations of trace elements, particularly Th, U, and Y, the outer regions (zircon II) contain significantly less trace elements, numerous pores, and inclusions of thorite, ytttrialite, and fergusonite. Zircon II contains very low or undetectable concentrations of non-formula elements such as Ca, Al, and Fe, which are commonly observed in high concentrations in altered radiation-damaged zircon. U-Pb dating of both zircon domains by LA-ICPMS and SHRIMP yielded statistically indistinguishable U-Pb weighted average ages of 119.3 2.1 (2r) and 118 1.2 (2r) Ma, respectively, demonstrating that the zircon had not accumulated a significant amount of self-irradiation damage at the time of the alteration event.Electron microprobe dating of thorite inclusions in zircon II yielded a Th-U-total Pb model age of 122 5 (2r) Ma, supporting the age relationship between both zircon domains. The hydrothermal solution responsible for the alteration of the investigated zircon was alkaline and rich in CO3 2-, as suggested by the occurrence of REE carbonates and CO2-bearing fluid inclusions. The alteration of the crystalline, trace element-rich zircon is explained by an interface-coupled dissolution-reprecipitation mechanism. During such a process, the congruent dissolution of the trace element-rich parent zircon I was spatially and temporally coupled to the precipitation of the trace element-poor zircon II at an inward moving dissolution-precipitation front. The driving force for such a process was merely the difference between the solubility of the trace element-rich and -poor zircon in the hydrothermal solution. The replacement process and the occurrence of mineral inclusions and porosity in the product zircon II is explained by the thermodynamics of solid solution-aqueous solution systems