Effect of solution saturation state and temperature on diopside dissolution

Steady-state dissolution rates of diopside are measured as a function of solution saturation state using a titanium flow-through reactor at pH 7.5 and temperature ranging from 125 to 175°C. Diopside dissolved stoichiometrically under all experimental conditions and rates were not dependent on sample history. At each temperature, rates continuously decreased by two orders of magnitude as equilibrium was approached and did not exhibit a dissolution plateau of constant rates at high degrees of undersaturation. The variation of diopside dissolution rates with solution saturation can be described equally well with a ion exchange model based on transition state theory or pit nucleation model based on crystal growth/dissolution theory from 125 to 175°C. At 175°C, both models over predict dissolution rates by two orders of magnitude indicating that a secondary phase precipitated in the experiments. The ion exchange model assumes the formation of a Si-rich, Mg-deficient precursor complex. Lack of dependence of rates on steady-state aqueous calcium concentration supports the formation of such a complex, which is formed by exchange of protons for magnesium ions at the surface. Fit to the experimental data yields [Formula: see text] where the Mg-H exchange coefficient, n = 1.39, the apparent activation energy, E(a )= 332 kJ mol(-1), and the apparent rate constant, k = 10(41.2 )mol diopside cm(-2 )s(-1). Fits to the data with the pit nucleation model suggest that diopside dissolution proceeds through retreat of steps developed by nucleation of pits created homogeneously at the mineral surface or at defect sites, where homogeneous nucleation occurs at lower degrees of saturation than defect-assisted nucleation. Rate expressions for each mechanism (i) were fit to [Formula: see text] where the step edge energy (α) for homogeneously nucleated pits were higher (275 to 65 mJ m(-2)) than the pits nucleated at defects (39 to 65 mJ m(-2)) and the activation energy associated with the temperature dependence of site density and the kinetic coefficient for homogeneously nucleated pits (E(b-homogeneous )= 2.59 × 10(-16 )mJ K(-1)) were lower than the pits nucleated at defects (E(b-defect assisted )= 8.44 × 10(-16 )mJ K(-1))

A 4D view on the evolution of metamorphic dehydration reactions

Metamorphic reactions influence the evolution of the Earth's crust in a range of tectonic settings. For example hydrous mineral dehydration in a subducting slab can produce fluid overpressures which may trigger seismicity. During reaction the mechanisms of chemical transport, including water expulsion, will dictate the rate of transformation and hence the evolution of physical properties such as fluid pressure. Despite the importance of such processes, direct observation of mineral changes due to chemical transport during metamorphism has been previously impossible both in nature and in experiment. Using time-resolved (4D) synchrotron X-ray microtomography we have imaged a complete metamorphic reaction and show how chemical transport evolves during reaction. We analyse the dehydration of gypsum to form bassanite and H2O which, like most dehydration reactions, produces a solid volume reduction leading to the formation of pore space. This porosity surrounds new bassanite grains producing fluid-filled moats, across which transport of dissolved ions to the growing grains occurs via diffusion. As moats grow in width, diffusion and hence reaction rate slow down. Our results demonstrate how, with new insights into the chemical transport mechanisms, we can move towards a more fundamental understanding of the hydraulic and chemical evolution of natural dehydrating systems

Corrosion of the International Simple Glass under acidic to hyperalkaline conditions

Assessment of glass dissolution kinetics, under disposal relevant temperature and pH environments, is required to credibly estimate radionuclide release rates from vitrified radioactive waste. Leaching of the International Simple Glass (ISG) under acidic to hyperalkaline conditions was examined. Forward rate measurements have been obtained using the dynamic leaching SPFT protocol and rate parameters for B, Na and Si in the basic regime; errors in rates predicted using these parameters at high pH and temperature are significant because the fitting uses logarithmic data. Longer term behaviour under hyperalkaline conditions, representative of some disposal environments, was investigated using the PCT and MCC-1 static leaching protocols with Ca(OH)2 solutions for up to 120 days (PCT) and 720 days (MCC-1). In hyperalkaline conditions dissolution was incongruent for all elements and the presence of alternating zirconia-rich and zirconia-poor alteration layers was observed on all leached monoliths, indicating the occurrence of a self-organisation phenomenon during leaching

Modeling acid-gas generation from boiling chloride brines

<p>Abstract</p> <p>Background</p> <p>This study investigates the generation of HCl and other acid gases from boiling calcium chloride dominated waters at atmospheric pressure, primarily using numerical modeling. The main focus of this investigation relates to the long-term geologic disposal of nuclear waste at Yucca Mountain, Nevada, where pore waters around waste-emplacement tunnels are expected to undergo boiling and evaporative concentration as a result of the heat released by spent nuclear fuel. Processes that are modeled include boiling of highly concentrated solutions, gas transport, and gas condensation accompanied by the dissociation of acid gases, causing low-pH condensate.</p> <p>Results</p> <p>Simple calculations are first carried out to evaluate condensate pH as a function of HCl gas fugacity and condensed water fraction for a vapor equilibrated with saturated calcium chloride brine at 50-150°C and 1 bar. The distillation of a calcium-chloride-dominated brine is then simulated with a reactive transport model using a brine composition representative of partially evaporated calcium-rich pore waters at Yucca Mountain. Results show a significant increase in boiling temperature from evaporative concentration, as well as low pH in condensates, particularly for dynamic systems where partial condensation takes place, which result in enrichment of HCl in condensates. These results are in qualitative agreement with experimental data from other studies.</p> <p>Conclusion</p> <p>The combination of reactive transport with multicomponent brine chemistry to study evaporation, boiling, and the potential for acid gas generation at the proposed Yucca Mountain repository is seen as an improvement relative to previously applied simpler batch evaporation models. This approach allows the evaluation of thermal, hydrological, and chemical (THC) processes in a coupled manner, and modeling of settings much more relevant to actual field conditions than the distillation experiment considered. The actual and modeled distillation experiments do not represent expected conditions in an emplacement drift, but nevertheless illustrate the potential for acid-gas generation at moderate temperatures (<150°C).</p

Geophysical monitoring and reactive transport modeling of ureolytically-driven calcium carbonate precipitation

Ureolytically-driven calcium carbonate precipitation is the basis for a promising in-situ remediation method for sequestration of divalent radionuclide and trace metal ions. It has also been proposed for use in geotechnical engineering for soil strengthening applications. Monitoring the occurrence, spatial distribution, and temporal evolution of calcium carbonate precipitation in the subsurface is critical for evaluating the performance of this technology and for developing the predictive models needed for engineering application. In this study, we conducted laboratory column experiments using natural sediment and groundwater to evaluate the utility of geophysical (complex resistivity and seismic) sensing methods, dynamic synchrotron x-ray computed tomography (micro-CT), and reactive transport modeling for tracking ureolytically-driven calcium carbonate precipitation processes under site relevant conditions. Reactive transport modeling with TOUGHREACT successfully simulated the changes of the major chemical components during urea hydrolysis. Even at the relatively low level of urea hydrolysis observed in the experiments, the simulations predicted an enhanced calcium carbonate precipitation rate that was 3-4 times greater than the baseline level. Reactive transport modeling results, geophysical monitoring data and micro-CT imaging correlated well with reaction processes validated by geochemical data. In particular, increases in ionic strength of the pore fluid during urea hydrolysis predicted by geochemical modeling were successfully captured by electrical conductivity measurements and confirmed by geochemical data. The low level of urea hydrolysis and calcium carbonate precipitation suggested by the model and geochemical data was corroborated by minor changes in seismic P-wave velocity measurements and micro-CT imaging; the latter provided direct evidence of sparsely distributed calcium carbonate precipitation. Ion exchange processes promoted through NH4+ production during urea hydrolysis were incorporated in the model and captured critical changes in the major metal species. The electrical phase increases were potentially due to ion exchange processes that modified charge structure at mineral/water interfaces. Our study revealed the potential of geophysical monitoring for geochemical changes during urea hydrolysis and the advantages of combining multiple approaches to understand complex biogeochemical processes in the subsurface

Decadal to monthly timescales of magma transfer and reservoir growth at a caldera volcano

International audienceCaldera-forming volcanic eruptions are low-frequency, highimpact events capable of discharging tens to thousands of cubic kilometres of magma explosively on timescales of hours to days, with devastating effects on local and global scales1. Because no such eruption has been monitored during its long build-up phase, the precursor phenomena are not well understood. Geophysical signals obtained during recent episodes of unrest at calderas such as Yellowstone, USA, and Campi Flegrei, Italy, are difficult to interpret, and the conditions necessary for large eruptions are poorly constrained2,3. Here we present a study of pre-eruptive magmatic processes and their timescales using chemically zoned crystals from the 'Minoan' caldera-formingeruption of Santorini volcano,Greece4, which occurred in the late 1600s BC. The results provide insights into how rapidly large silicic systems may pass from a quiescent state to one on the edge of eruption5,6. Despite the large volume of erupted magma4 (40-60 cubic kilometres), and the 18,000-year gestation period between the Minoan eruption and the previous major eruption, most crystals in the Minoan magma record processes that occurred less than about 100 years before the eruption. Recharge of the magma reservoir by large volumes of silicic magma (and some mafic magma) occurred during the century before eruption, and mixing between different silicicmagmabatches was still taking place during the final months. Final assembly of large silicic magma reservoirs may occur on timescales that are geologically very short by comparison with the preceding repose period, with major growth phases immediately before eruption. These observations have implications for the monitoring of long-dormant, but potentially active, caldera systems

Atmospheric oxygen regulation at low Proterozoic levels by incomplete oxidative weathering of sedimentary organic carbon

It is unclear why atmospheric oxygen remained trapped at low levels for more than 1.5 billion years following the Paleoproterozoic Great Oxidation Event. Here, we use models for erosion, weathering and biogeochemical cycling to show that this can be explained by the tectonic recycling of previously accumulated sedimentary organic carbon, combined with the oxygen sensitivity of oxidative weathering. Our results indicate a strong negative feedback regime when atmospheric oxygen concentration is of order pO2∼0.1 PAL (present atmospheric level), but that stability is lost at pO2<0.01 PAL. Within these limits, the carbonate carbon isotope (δ13C) record becomes insensitive to changes in organic carbon burial rate, due to counterbalancing changes in the weathering of isotopically light organic carbon. This can explain the lack of secular trend in the Precambrian δ13C record, and reopens the possibility that increased biological productivity and resultant organic carbon burial drove the Great Oxidation Event

Role of hydrodynamic factors in controlling the formation and location of unconformity-related uranium deposits: insights from reactive-flow modeling

The role of hydrodynamic factors in controlling the formation and location of unconformity-related uranium (URU) deposits in sedimentary basins during tectonically quiet periods is investigated. A number of reactive-flow modeling experiments at the deposit scale were carried out by assigning different dip angles and directions to a fault and various permeabilities to hydrostratigraphic units). The results show that the fault dip angle and direction, and permeability of the hydrostratigraphic units govern the convection pattern, temperature distribution, and uranium mineralization. Avertical fault results in uranium mineralization at the bottom of the fault within the basement, while a dipping fault leads to precipitation of uraninite below the unconformity either away from or along the plane of the fault, depending on the fault permeability. A more permeable fault causes uraninite precipitates along the fault plane,whereas a less permeable one gives rise to the precipitation of uraninite away from it. No economic ore mineralization can form when either very low or very high permeabilities are assigned to the sandstone or basement suggesting that these units seem to have an optimal window of permeability for the formation of uranium deposits. Physicochemical parameters also exert an additional control in both the location and grade of URU deposits. These results indicate that the difference in size and grade of different URU deposits may result from variation in fluid flow pattern and physicochemical conditions, caused by the change in structural features and hydraulic properties of the stratigraphic units involved