Low-temperature nucleation in a kinetic Ising model with soft stochastic dynamics

We study low-temperature nucleation in kinetic Ising models by analytical and simulational methods, confirming the general result for the average metastable lifetime, = A*exp(beta*Gamma) (beta = 1/kT) [E. Jordao Neves and R.H. Schonmann, Commun. Math. Phys. 137, 209 (1991)]. Contrary to common belief, we find that both A and Gamma depend significantly on the stochastic dynamic. In particular, for a ``soft'' dynamic, in which the effects of the interactions and the applied field factorize in the transition rates, Gamma does NOT simply equal the energy barrier against nucleation, as it does for the standard Glauber dynamic, which does not have this factorization property.Comment: 4 pages RevTex4, 2 figures. Phys. Rev. Lett., in pres

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))

Strontium as a tracer of weathering processes in a silicate catchment polluted by acid atmospheric inputs, Strengbach, France

This paper determines the weathering and atmospheric contributions of Ca in surface water from a small spruce forested silicate catchment (N–E France) receiving acid atmospheric inputs. The bedrock is a granite with K-feldspar and albite as dominant phases. The calcium content in plagioclase is low and the Ca/Na ratio in surface water is high, reflecting other sources of calcium from those expected from the weathering of major mineral phases. The biotite content is low. Only traces of apatite were detected while no calcite was found in spite of a major hydrothermal event having affected the granite. The strontium isotopic ratio 87Sr/86Sr and Sr content was used as a tracer of weathering and was determined in minerals and bulk bedrock, open field precipitation, throughfall, soil solution, spring and stream water. The Sr isotopic ratio of the reacting weathering end-member was predicted by simulating the alteration of the granite minerals by incorporating strontium into the water–rock interaction kinetic code KINDIS. In the early stages of water–rock interaction, K-feldspar and biotite strongly influence the isotopic composition of the weathering solution whereas, the Na-rich plagioclase appears to be the main long-term reactive weathering end-member. Approximately 50% of dissolved Sr in streamwater are atmospherically derived. The 87Sr/86Sr ratios of exchangeable Sr in the fine fraction at 1-m depth from a soil profile indicate that the amount of exchangeable Sr seems essentially controlled by atmospheric inputs. The exception is the deep saprolite where weathering processes could supply the Sr (and Ca). Na-Plagioclase weathering obviously control the chemistry and the isotopic composition of surface waters. The weathering of trace mineral plays a secondary role, the exception is for apatite when plagioclase is absent. Our hydrochemical, mineralogical and isotopic investigations show that a major part of the strong Ca losses detected in catchment hydrochemical budgets that result from the neutralization of acid precipitation has an atmospheric origin. Consequently, in the long term, in such areas, the availability of such an exchangeable base cation might be strongly limited and surface waters consequently acidified

Holistic approach to dissolution kinetics : linking direction-specific microscopic fluxes, local mass transport effects and global macroscopic rates from gypsum etch pit analysis

Dissolution processes at single crystal surfaces often involve the initial formation and expansion of localized, characteristic (faceted) etch-pits at defects, in an otherwise comparatively unreactive surface. Using natural gypsum single crystal as an example, a simple but powerful morphological analysis of these characteristic etch pit features is proposed that allows important questions concerning dissolution kinetics to be addressed. Significantly, quantitative mass transport associated with reactive microscale interfaces in quiescent solution (well known in the field of electrochemistry at ultramicroelectrodes) allows the relative importance of diffusion compared to surface kinetics to be assessed. Furthermore, because such mass transport rates are high, much faster surface kinetics can be determined than with existing dissolution methods. For the case of gypsum, surface processes are found to dominate the kinetics at early stages of the dissolution process (small etch pits) on the cleaved (010) surface. However, the contribution from mass transport becomes more important with time due to the increased area of the reactive zones and associated decrease in mass transport rate. Significantly, spatial heterogeneities in both surface kinetics and mass transport effects are identified, and the morphology of the characteristic etch features reveal direction-dependent dissolution kinetics that can be quantified. Effective dissolution velocities normal to the main basal (010) face are determined, along with velocities for the movement of [001] and [100] oriented steps. Inert electrolyte enhances dissolution velocities in all directions (salting in), but a striking new observation is that the effect is direction-dependent. Studies of common ion effects reveal that Ca2+ has a much greater impact in reducing dissolution rates compared to SO42−. With this approach, the new microscopic observations can be further analysed to obtain macroscopic dissolution rates, which are found to be wholly consistent with previous bulk measurements. The studies are thus important in bridging the gap between microscopic phenomena and macroscopic measurements

Mechanisms of goethite dissolution in the presence of desferrioxamine B and Suwannee River fulvic acid at pH 6.5

Siderophores are Fe3+ specific low MW chelating ligands secreted by microorganisms in response to Fe stress. Low MW organic acids such as oxalate have been shown to enhance siderophore mediated dissolution of Fe3+ oxides. However, the effect of fulvic acid presence on siderophore function remains unknown. We used batch dissolution experiments to investigate Fe release from goethite in the goethite-fulvic acid desferrioxamine B (goethite-SRFA-DFOB) ternary system. Experiments were conducted at pH 6.5 while varying reagent addition sequence. FTIR and UV-Vis spectroscopy were employed to characterise the Fe-DFOB, Fe-SRFA and DFOB–SRFA complexes. Iron released from goethite in the presence of SRFA alone was below detection limit. In the presence of both SRFA and DFOB, dissolved Fe increased with reaction time, presence of the DFOB-SRFA complex, and where SRFA was introduced prior to DFOB. FTIR data show that in the ternary system, Fe3+ is complexed primarily to oxygen of the DFOB hydroxamate group, whilst the carboxylate C=O of SRFA forms an electrostatic association with the terminal NH3+ of DFOB. We propose that SRFA sorbed to goethite lowers the net positive charge of the oxide surface, thus facilitating adsorption of cationic DFOB and subsequent Fe3+ chelation and release. Furthermore, the sorbed SRFA weakens Fe-O bonds at the goethite surface, increasing the population of kinetically labile Fe. This work demonstrates the positive, though indirect role of SRFA in increasing the bioavailability of Fe3+

Toward catchment hydro-biogeochemical theories

Headwater catchments are the fundamental units that connect the land to the ocean. Hydrological flow and biogeochemical processes are intricately coupled, yet their respective sciences have progressed without much integration. Reaction kinetic theories that prescribe rate dependence on environmental variables (e.g., temperature and water content) have advanced substantially, mostly in well-mixed reactors, columns, and warming experiments without considering the characteristics of hydrological flow at the catchment scale. These theories have shown significant divergence from observations in natural systems. On the other hand, hydrological theories, including transit time theory, have progressed substantially yet have not been incorporated into understanding reactions at the catchment scale. Here we advocate for the development of integrated hydro-biogeochemical theories across gradients of climate, vegetation, and geology conditions. The lack of such theories presents barriers for understanding mechanisms and forecasting the future of the Critical Zone under human- and climate-induced perturbations. Although integration has started and co-located measurements are well under way, tremendous challenges remain. In particular, even in this era of "big data," we are still limited by data and will need to (1) intensify measurements beyond river channels and characterize the vertical connectivity and broadly the shallow and deep subsurface; (2) expand to older water dating beyond the time scales reflected in stable water isotopes; (3) combine the use of reactive solutes, nonreactive tracers, and isotopes; and (4) augment measurements in environments that are undergoing rapid changes. To develop integrated theories, it is essential to (1) engage models at all stages to develop model-informed data collection strategies and to maximize data usage; (2) adopt a "simple but not simplistic," or fit-for-purpose approach to include essential processes in process-based models; (3) blend the use of process-based and data-driven models in the framework of "theory-guided data science." Within the framework of hypothesis testing, model-data fusion can advance integrated theories that mechanistically link catchments' internal structures and external drivers to their functioning. It can not only advance the field of hydro-biogeochemistry, but also enable hind- and fore-casting and serve the society at large. Broadly, future education will need to cultivate thinkers at the intersections of traditional disciplines with hollistic approaches for understanding interacting processes in complex earth systems.This article is categorized under:Science of Water > Method

Tracking the weathering of basalts on Mars using lithium isotope fractionation models

An edited version of this paper was published by AGU. Copyright (2015) American Geophysical UnionLithium (Li), the lightest of the alkali elements, has geochemical properties that include high aqueous solubility (Li is the most fluid mobile element) and high relative abundance in basalt-forming minerals (values ranking between 0.2 and 12 ppm). Li isotopes are particularly subject to fractionation because the two stable isotopes of lithium - 7Li and 6Li - have a large relative mass difference (∼15%) that results in significant fractionation between water and solid phases. The extent of Li isotope fractionation during aqueous alteration of basalt depends on the dissolution rate of primary minerals - the source of Li - and on the precipitation kinetics, leading to formation of secondary phases. Consequently, a detailed analysis of Li isotopic ratios in both solution and secondary mineral lattices could provide clues about past Martian weathering conditions, including weathering extent, temperature, pH, supersaturation, and evaporation rate of the initial solutions in contact with basalt rocks. In this paper, we discuss ways in which Martian aqueous processes could have lead to Li isotope fractionation. We show that Li isotopic data obtained by future exploration of Mars could be relevant to highlighting different processes of Li isotopic fractionation in the past, and therefore to understanding basalt weathering and environmental conditions early in the planet's historyData supporting our models and calculations are available as supporting information. The research leading to these results is a contribution from the Project ‘icyMARS’’, funded by the European Research Council, Starting Grant no 307496. This work was also partially supported by the European FEDER program and the Spanish Ministry of Science (MICINN) through the project CGL2011–30079. Comments by R. James and four anonymous reviewers helped us to clarify and strengthen our wor

Primary weathering rates, water transit times, and concentration-discharge relations:a theoretical analysis for the critical zone

The permeability architecture of the critical zone exerts a major influence on the hydrogeochemistry of the critical zone. Water flow path dynamics drive the spatiotemporal pattern of geochemical evolution and resulting streamflow concentration-discharge (C-Q) relation, but these flow paths are complex and difficult to map quantitatively. Here we couple a new integrated flow and particle tracking transport model with a general reversible Transition State Theory style dissolution rate law to explore theoretically how C-Q relations and concentration in the critical zone respond to decline in saturated hydraulic conductivity (K-s) with soil depth. We do this for a range of flow rates and mineral reaction kinetics. Our results show that for minerals with a high ratio of equilibrium concentration (C-eq) to intrinsic weathering rate (R-max), vertical heterogeneity in K-s enhances the gradient of weathering-derived solute concentration in the critical zone and strengthens the inverse stream C-Q relation. As C-eq/R-max decreases, the spatial distribution of concentration in the critical zone becomes more uniform for a wide range of flow rates, and stream C-Q relation approaches chemostatic behavior, regardless of the degree of vertical heterogeneity in K-s. These findings suggest that the transport-controlled mechanisms in the hillslope can lead to chemostatic C-Q relations in the stream while the hillslope surface reaction-controlled mechanisms are associated with an inverse stream C-Q relation. In addition, as C-eq/R-max decreases, the concentration in the critical zone and stream become less dependent on groundwater age (or transit time)