Isotopic fractionation of noble gases by diffusion in liquid water: Molecular dynamics simulations and hydrologic applications

Despite their great importance in low-temperature geochemistry, diffusion coefficients of noble gas isotopes in liquid water (D) have been measured only for the major isotopes of helium, neon, krypton and xenon. Data on the diffusion coefficients of minor noble gas isotopes are essentially non-existent and so typically have been estimated by a kinetic-theory model in which D varies as the inverse square root of the isotopic mass (m): D ∝ m^−0.5. To examine the validity of the kinetic-theory model, we performed molecular dynamics (MD) simulations of the diffusion of noble gases in ambient liquid water. Our simulation results agree with available experimental data on the solvation structure and diffusion coefficients of the major noble gas isotopes and reveal for the first time that the isotopic mass-dependence of all noble gas self-diffusion coefficients has the power-law form D ∝ m^−β with 0 < β < 0.2. Thus our results call into serious question the widespread assumption that the ‘square-root’ model can be applied to estimate the kinetic fractionation of noble gas isotopes caused by diffusion in ambient liquid water. To illustrate the importance of this finding, we used the diffusion coefficients determined in our MD simulations to reconsider the geochemical modeling of ²⁰Ne/²²Ne and ³⁶Ar/⁴⁰Ar isotopic ratios in three representative hydrologic studies. Our new modeling results indicate that kinetic isotopic fractionation by diffusion may play a significant role in noble gas transport processes in groundwater.     

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.     

Diffusion of HTO,Cl and I in Upper Toarcian argillite samples from Tournemire: Effects of initial iodide concentration and ionic strength

Diffusion parameters for HTO, ³⁶Cl⁻_, and ¹²⁵I⁻ were determined on Upper Toarcian argillite samples from the Tournemire Underground Research Laboratory (Aveyron, France) using the through diffusion technique. The direction of diffusion was parallel to the bedding plane. The purpose of the present study was 3-fold; it was intended (i) to confirm the I⁻ interaction with Upper Toarcian argillite and to verify the effects of initial I⁻ concentration on this affinity, as previously observed by means of radial diffusion experiments, (ii) to highlight any discrepancy between Cl⁻ and I⁻ diffusivity, and (iii) to investigate the effect of an increase of the ionic strength of the solution on the anionic tracers’ diffusive behaviour. The results show that the effective diffusion coefficient (D_e) and diffusion accessible porosity (ε_a) values obtained with an ionic strength (I.S.) synthetic pore water of 0.01 eq L⁻¹ are: D_e = 2.35–2.50 × 10⁻¹¹ m² s⁻¹ and ε_a = 12.0–15.0% for HTO, and D_e = 14.5–15.5 × 10⁻¹³ m² s⁻¹ and ε_a = 2.5–2.9% for ³⁶Cl⁻. Because of anionic exclusion effects, anions diffuse slower and exhibit smaller diffusion accessible porosities than HTO, taken as a water tracer. The associated effective diffusion coefficient (D_e) and rock capacity factor (α) obtained for ¹²⁵I are: D_e = 7.00–8.60 × 10⁻¹³ m² s⁻¹ and α = 4.3–7.2%. Such values make it possible to calculate low ¹²⁵I distribution ratios (0.0057 < R_D < 0.0192 mL g⁻¹) which confirm the trend indicating that the ¹²⁵I rock capacity factor increases with the decrease of the initial I⁻ concentration. Additional through-diffusion experiments were carried out with a higher ionic strength synthetic pore water (I.S. = 0.11 eq L⁻¹). No evolution of HTO diffusion parameters was observed. The anionic tracers’ effective diffusion coefficient increased by a factor of two but no clear evolution of their accessible porosity was observed. Such a paradox could be related to the particularly small mean pore size of the Upper Toarcian argillite of Tournemire. The most significant finding of this study is the large discrepancy (factor of two) between the values of the effective diffusion coefficient for ¹²⁵I and ³⁶Cl. Whatever the ionic strength of the synthetic solution used, ¹²⁵I exhibited D_e values two times lower than those of ³⁶Cl. A detailed explanation for this difference cannot be given at present even if a hypothesis based on ion-pairing or on steric-exclusion cannot be excluded. This makes questionable the assumption usually made for quantifying ¹²⁵I sorption and postulating that ³⁶Cl and ¹²⁵I would diffuse in the same porosity. In other terms, at Tournemire, ¹²⁵I sorption could be more pronounced than previously indicated.     

Olivine/melt transition metal partitioning, melt composition, and melt structure—Melt polymerization and Q-speciation in alkaline earth silicate systems

The two most abundant network-modifying cations in magmatic liquids are Ca²⁺ and Mg²⁺. To evaluate the influence of melt structure on exchange of Ca²⁺ and Mg²⁺ with other geochemically important divalent cations (m-cations) between coexisting minerals and melts, high-temperature (1470-1650 °C), ambient-pressure (0.1 MPa) forsterite/melt partitioning experiments were carried out in the system Mg₂SiO₄-CaMgSi₂O₆-SiO₂ with ?1 wt% m-cations (Mn²⁺, Co²⁺, and Ni²⁺) substituting for Ca²⁺ and Mg²⁺. The bulk melt NBO/Si-range (NBO/Si: nonbridging oxygen per silicon) of melt in equilibrium with forsterite was between 1.89 and 2.74. In this NBO/Si-range, the NBO/Si(Ca) (fraction of nonbridging oxygens, NBO, that form bonds with Ca²⁺, Ca²⁺-NBO) is linearly related to NBO/Si, whereas fraction of Mg²⁺-NBO bonds is essentially independent of NBO/Si. For individual m-cations, rate of change of K_D(m−Mg) with NBO/Si(Ca) for the exchange equilibrium, m^melt + Mg^olivine ? m^olivine + Mg^melt, is linear. K_D(m−Mg) decreases as an exponential function of increasing ionic potential, Z/r² (Z: formal electrical charge, r: ionic radius—here calculated with oxygen in sixfold coordination around the divalent cations) of the m-cation. The enthalpy change of the exchange equilibrium, ΔH, decreases linearly with increasing Z/r² [ΔH = 261(9)-81(3)·Z/r² (Å⁻²)]. From existing information on (Ca,Mg)O-SiO₂ melt structure at ambient pressure, these relationships are understood by considering the exchange of divalent cations that form bonds with nonbridging oxygen in individual Qⁿ-species in the melts. The negative ∂K_D(m−Mg)/∂(Z/r²) and ∂(ΔH)/∂(Z/r²) is because increasing Z/r² is because the cations forming bonds with nonbridging oxygen in increasingly depolymerized Qⁿ-species where steric hindrance is decreasingly important. In other words, principles of ionic size/site mismatch commonly observed for trace and minor elements in crystals, also govern their solubility behavior in silicate melts.     

A threshold in soil formation at Earth’s arid-hyperarid transition

The soils of the Atacama Desert in northern Chile have long been known to contain large quantities of unusual salts, yet the processes that form these soils are not yet fully understood. We examined the morphology and geochemistry of soils on post-Miocene fans and stream terraces along a south-to-north (27° to 24° S) rainfall transect that spans the arid to hyperarid transition (21 to ∼2 mm rain y⁻¹). Landform ages are ? 2 My based on cosmogenic radionuclide concentrations in surface boulders, and Ar isotopes in interbedded volcanic ash deposits near the driest site indicate a maximum age of 2.1 My. A chemical mass balance analysis that explicitly accounts for atmospheric additions was used to quantify net changes in mass and volume as a function of rainfall. In the arid (21 mm rain y⁻¹) soil, total mass loss to weathering of silicate alluvium and dust (−1030 kg m⁻²) is offset by net addition of salts (+170 kg m⁻²). The most hyperarid soil has accumulated 830 kg m⁻² of atmospheric salts (including 260 kg sulfate m⁻² and 90 kg chloride m⁻²), resulting in unusually high volumetric expansion (120%) for a soil of this age. The composition of both airborne particles and atmospheric deposition in passive traps indicates that the geochemistry of the driest soil reflects accumulated atmospheric influxes coupled with limited in-soil chemical transformation and loss. Long-term rates of atmospheric solute addition were derived from the ion inventories in the driest soil, divided by the landform age, and compared to measured contemporary rates. With decreasing rainfall, the soil salt inventories increase, and the retained salts are both more soluble and present at shallower depths. All soils generally exhibit vertical variation in their chemistry, suggesting slow and stochastic downward water movement, and greater climate variability over the past 2 My than is reflected in recent (∼100 y) rainfall averages. The geochemistry of these soils shows that the transition from arid to hyperarid rainfall levels marks a fundamental geochemical threshold: in wetter soils, the rate and character of chemical weathering results in net mass loss and associated volumetric collapse after 10⁵ to 10⁶ years, while continuous accumulation of atmospheric solutes in hyperarid soils over similar timescales results in dramatic volumetric expansion. The specific geochemistry of hyperarid soils is a function of atmospheric sources, and is expected to vary accordingly at other hyperarid sites. This work identifies key processes in hyperarid soil formation that are likely to be independent of location, and suggests that analogous processes may occur on Mars.     

Isotopic fractionation of the major elements of molten basalt by chemical and thermal diffusion

Samples produced in piston cylinder experiments were used to document the thermal isotopic fractionation of all the major elements of basalt except for aluminum and the fractionation of iron isotopes by chemical diffusion between a natural basalt and rhyolite. The thermal isotopic fractionations are summarized in terms of a parameter Ω_i defined as the fractionation in per mil per 100 °C per atomic mass units difference between the isotopes. For molten basalt we report Ω_Ca = 1.6, Ω_Fe = 1.1, Ω_Si = 0.6, Ω_O = 1.5. In an earlier paper we reported Ω_Mg = 3.6. These fractionations represent a steady state balance between thermal diffusion and chemical diffusion with the mass dependence of the thermal diffusion coefficient being significantly larger than the mass dependence of the chemical diffusion coefficients for isotopes of the same element. The iron isotopic measurements of the basalt-rhyolite diffusion couple showed significant fractionation that are parameterized in terms of a parameter β_Fe = 0.03 when the ratio of the diffusion coefficients D₅₄ and D₅₆ of ⁵⁴Fe and ⁵⁶Fe is expressed in terms of the atomic mass as D₅₄/D₅₆ = (56/54)^β_Fe. This value of β_Fe is smaller than what we had measured earlier for lithium, magnesium and calcium (i.e., β_Li = 0.215, β_Ca = 0.05, β_Mg = 0.05) but still significant when one takes into account the high precision with which iron isotopic compositions can be measured (i.e., ±0.03‰) and that iron isotope fractionations at magmatic temperatures from other causes are extremely small. In a closing section we discuss technological and geological applications of isotopic fractionations driven by either or both chemical and thermal gradients.     

Monomethylmercury sources in a tropical artificial reservoir

The distribution and speciation of mercury (Hg) in the water column, the inputs (wet deposition and tributaries) and the outputs (atmospheric evasion and outlet) of an artificial partially anoxic tropical lake (Petit-Saut reservoir, French Guiana) were investigated on a seasonal basis in order to appraise the cycling and transformations of this metal. The total mercury (HgT) concentrations in the oxygenated epilimnetic waters averaged 5 ± 3 pmol L⁻¹ in the unfiltered samples (HgT_UNF) and 4 ± 2 pmol L⁻¹ in the dissolved (HgT_D) phase (<0.45 μm). On average, the monomethylmercury (MMHg) constituted 8%, 40% and 18% of the HgT in the dissolved phase, the particulate suspended matter and in the unfiltered samples, respectively. Covariant elevated concentrations of particulate MMHg and chlorophyll a in the epilimnion suggest that phytoplankton is an active component for the MMHg transfer in the lake. In the anoxic hypolimnion the HgT_UNF averages 13 ± 6 pmol L⁻¹ and the HgT_D 8 ± 4 pmol L⁻¹. The averages of MMHg_P and MMHg_D in hypolimnetic waters were two and three times the corresponding values of the epilimnion, 170 ± 90 pmol g⁻¹ and 0.9 ± 0.5 pmol L⁻¹, respectively. In the long dry and wet seasons, at the flooded forest and upstream dam sampling stations, the vertical profiles of MMHg_D concentrations accounted for two distinct maxima: one just below the oxycline and the other near the benthic interface. Direct wet atmospheric deposition accounted for 14 moles yr⁻¹ HgT_UNF, with 0.7 moles yr⁻¹ as MMHg_UNF, while circa 76 moles yr⁻¹ of HgT_UNF, with 4.7 moles yr⁻¹ as MMHg_UNF, coming from tributaries. Circa 78 moles (∼17% as MMHg) are annually exported through the dam, while 23 moles yr⁻¹ of Hg⁰ evolve in the atmosphere. A mass balance calculation suggests that the endogenic production of MMHg_UNF attained 8.1 moles yr⁻¹, corresponding to a methylation rate of 0.06% d⁻¹. As a result, the Petit-Saut reservoir is a large man-made reactor that has extensively altered mercury speciation in favor of methylated species.     

Dissolution kinetics as a function of the Gibbs free energy of reaction: An experimental study based on albite feldspar

Here we report on an experimental investigation of the relation between the dissolution rate of albite feldspar and the Gibbs free energy of reaction, ΔG_r. The experiments were carried out in a continuously stirred flow-through reactor at 150 °C and pH_(150 °C) 9.2. The dissolution rates R are based on steady-state Si and Al concentrations and sample mass loss. The overall relation between ΔG_r and R was determined over a free energy range of −150 < ΔG_r < −15.6 kJ mol⁻¹. The data define a continuous and highly non-linear, sigmoidal relation between R and ΔG_r that is characterized by three distinct free energy regions. The region furthest from equilibrium, delimited by −150 < ΔG_r < −70 kJ mol⁻¹, represents an extensive dissolution rate plateau with an average rate . In this free energy range the rates of dissolution are constant and independent of ΔG_r, as well as [Si] and [Al]. The free energy range delimited by −70 ? ΔG_r ? −25 kJ mol⁻¹, referred to as the ‘transition equilibrium’ region, is characterized by a sharp decrease in dissolution rates with increasing ΔG_r, indicating a very strong inverse dependence of the rates on free energy. Dissolution nearest equilibrium, defined by ΔG_r > −25 kJ mol⁻¹, represents the ‘near equilibrium’ region where the rates decrease as chemical equilibrium is approached, but with a much weaker dependence on ΔG_r. The lowest rate measured in this study, R = 6.2 × 10⁻¹¹ mol m⁻² s⁻¹ at ΔG_r = −16.3 kJ mol⁻¹, is more than two orders of magnitude slower than the plateau rate. The data have been fitted to a rate equation (adapted from Burch et al. [Burch, T. E., Nagy, K. L., Lasaga, A. C., 1993. Free energy dependence of albite dissolution kinetics at 80 °C and pH 8.8. Chem. Geol.105, 137-162]) that represents the sum of two parallel reactions

R=k₁[1-exp(-ng^m₁)]+k₂[1-exp(-g)]^m₂,     

Effects of dry density and exchangeable cations on the diffusion process of sodium ions in compacted montmorillonite

As a part of the safety assessment of the geological disposal of high-level radioactive waste, the effects of dry density and exchangeable cations on the diffusion process of Na⁺ ions in compacted bentonite were studied from the viewpoint of the activation energy for diffusion. The apparent self-diffusion coefficients of Na⁺ ions in compacted Na-montmorillonite and in a Na- and Ca-montmorillonite mixture were determined by one-dimensional, non-steady diffusion experiments at different temperatures and dry densities. A unique change in activation energy as a function of dry density was found for the Na⁺ ions in compacted Na-montmorillonite. The activation energy suddenly decreased from 18.1 to 14.1 kJ mol^− 1 as the dry density increased from 0.9 to 1.0 Mg m^− 3, whereas it increased to 24.7 kJ mol^− 1 as the dry density increased to 1.8 Mg m^− 3. Examination of the effect of exchangeable cations on the activation energies determined that the activation energies were almost constant, approximately 25 kJ mol^− 1, for the montmorillonite specimens at a dry density of 1.8 Mg m^− 3. However, three different activation energy values were obtained at a dry density of 1.0 Mg m^− 3. These findings cannot be explained by the conventional diffusion model (the pore water diffusion model), which suggests that the predominant diffusion process alternates among pore water diffusion, interlayer diffusion, and external surface diffusion.     

Hydration of rhyolitic glass during weathering as characterized by IR microspectroscopy

The mechanism and rate of hydration of rhyolitic glass during weathering were studied. Doubly polished thin sections of two rhyolites with different duration of weathering (Ohsawa lava: 26,000 yr, Awanomikoto lava: 52,000 yr) were prepared. Optical microscope observation showed that altered layers had developed along the glass surfaces. IR spectral line profile analysis was conducted on the glass sections from the surface to the interior for a length of 250 μm and the contents of molecular H₂O (H₂O_m), OH species (OH) and total water (H₂O_t) were determined. The diffusion profile of H₂O_m in Ohsawa lava extends beyond the layer observed by optical microscope. The content of H₂O_m in the hydrated region is much higher than that of OH species. Thus, the reaction from H₂O_m to OH appears to be little and H₂O_m is the dominant water species moving into the glass during weathering. Based on the concentration profiles, the diffusion coefficients of H₂O_m(D_H2Om) and H₂O_t(D_H2Ot) were determined to be 2.8 × 10⁻¹⁰ and 3.4 × 10⁻¹⁰ μm² s⁻¹ for Ohsawa lava, and 5.2 × 10⁻¹¹ and 4.1 × 10⁻¹¹ μm² s⁻¹ for Awanomikoto lava, respectively. The obtained D_H2Om during weathering are more than 2-3 orders of magnitude larger than the diffusion coefficient at ∼20 °C that is extrapolated from the diffusivity data for >400 °C. This might suggest that the mechanism of water transport is different at weathering conditions and >400 °C.     

Halogen diffusion in a basaltic melt

The diffusion of the halogens fluorine, chlorine and bromine was measured in a hawaiitic melt from Mt. Etna at 500 MPa and 1.0 GPa, 1250 to 1450 °C at anhydrous conditions; the diffusion of F and Cl in the melt was also studied with about 3 wt% of dissolved water. Experiments were performed using the diffusion-couple technique in a piston cylinder. Most experiments were performed with only one halogen diffusing between the halogen-enriched and halogen-poor halves of the diffusion couple, but a few experiments with a mixture of halogens (F, Cl and Br) were also performed in order to investigate the possibility of interactions between the halogens during diffusion. Fluorine and chlorine diffusivity show a very similar behavior, slightly diverging at low temperature. Bromine diffusion is a factor of about 2-5 lower than the other halogens in this study. Diffusion coefficients for fluorine range between 2.3 × 10⁻¹¹ and 1.4 × 10⁻¹⁰ m² s⁻¹, for chlorine between 1.1 × 10⁻¹¹ and 1.3 × 10⁻¹⁰ and for bromine between 9.4 × 10⁻¹² and 6.8 × 10⁻¹¹ m² s⁻¹. No pressure effect was detected at the conditions investigated. In experiments involving mixed halogens, the diffusivities appear to decrease slightly (by a factor of ∼3), and are more uniform among the three elements. However, activation energies for diffusion do not appear to differ between experiments with individual halogens or when they are all mixed together. The effect of water increases the diffusion coefficients of F and Cl by no more than a factor of 3 compared to the anhydrous melt (D_F = 4.0 × 10⁻¹¹ to 1.6 × 10⁻¹⁰ m² s⁻¹; D_Cl = 3.0 × 10⁻¹¹ to 1.9 × 10⁻¹⁰ m² s⁻¹). Comparing our results to the diffusion coefficients of other volatiles in nominally dry basaltic melts, halogen diffusivities are about one order of magnitude lower than H₂O, similar to CO₂, and a factor of ∼5 higher than S. The contrasting volatile diffusivities may affect the variable extent of volatile degassing upon melt depressurization and vesiculation, and can help our understanding of the compositions of rapidly grown magmatic bubbles.     

Metal-silicate partitioning of cesium: Implications for core formation

Here we present the first set of metal-silicate partitioning data for Cs, which we use to examine whether the primitive mantle depletion of Cs can be attributed to core segregation. Our experiments independently varied pressure from 5 to 15 GPa, temperature from 1900 to 2400 °C, metallic sulfur content from pure Fe to pure FeS, silicate melt polymerization, expressed as a ratio of non-bridging oxygens to tetrahedrally coordinated cations (nbo/t) from 1.26 to 3.1, and fO₂ from two to four log units below the iron-wüstite buffer. The most important controls on the partitioning behavior of alkalis were the metallic sulfur content, expressed as X_S, and the nbo/t of the silicate liquid. Normalization of X_S to 0.5 yielded the following expressions for D-values as a function of nbo/t: log D_Na = −2.0 + 0.44 × (nbo/t), log D_K = −2.4 + 0.67 × ( nbo/t), and log D_Cs = −3.2 + 1.17 × (nbo/t). Normalization of nbo/t to 2.7 resulted in the following equations for D-values as a function of S content: log D_Na = −4.1 + 6.4 × X_S, log D_K = −7.7 + 13.9 × X_S, and log D_Cs = −12.1 + 23.3 × X_S.There appears to be a negative pressure effect up to 15 GPa, but it should be noted that this trend was not present before normalization, and is based on only two measurements. There is a positive trend in cesium’s metal-silicate partition coefficient with increasing temperature. D_Cs exhibits the largest change and increased by a factor of three over 500 °C. The effect of oxygen fugacity has not been precisely determined but in general, lowering fO₂ by two log units resulted in a rise in all D-values of approximately an order of magnitude. In general, the sensitivity of partition coefficients to changing parameters increased with atomic number.The highest D-value for Cs observed in this study is 0.345, which was obtained at nbo/t of 2.7 and a metal phase of pure FeS. This metallic composition has far more S than has been suggested for any credible core-forming metal. We therefore conclude that the depletion of Cs in Earth’s mantle is either caused by radically different behavior of Cs at pressures higher than 15 GPa or is not related to core formation. Even so, we have shown that a planet with a sufficient S inventory may incorporate significant amounts of alkali elements into its core.     

Benthic nitrogen cycling traversing the Peruvian oxygen minimum zone

Benthic nitrogen (N) cycling was investigated at six stations along a transect traversing the Peruvian oxygen minimum zone (OMZ) at 11°S. An extensive dataset including porewater concentration profiles and in situ benthic fluxes of nitrate (NO₃⁻), nitrite (NO₂⁻) and ammonium (NH₄⁺) was used to constrain a 1-D reaction-transport model designed to simulate and interpret the measured data at each station. Simulated rates of nitrification, denitrification, anammox and dissimilatory nitrate reduction to ammonium (DNRA) by filamentous large sulfur bacteria (e.g. Beggiatoa and Thioploca) were highly variable throughout the OMZ yet clear trends were discernible. On the shelf and upper slope (80-260 m water depth) where extensive areas of bacterial mats were present, DNRA dominated total N turnover (?2.9 mmol N m⁻² d⁻¹) and accounted for ?65% of NO₃⁻ + NO₂⁻ uptake by the sediments from the bottom water. Nonetheless, these sediments did not represent a major sink for dissolved inorganic nitrogen (DIN = NO₃⁻ + NO₂⁻ + NH₄⁺) since DNRA reduces NO₃⁻ and, potentially NO₂⁻, to NH₄⁺. Consequently, the shelf and upper slope sediments were recycling sites for DIN due to relatively low rates of denitrification and high rates of ammonium release from DNRA and ammonification of organic matter. This finding contrasts with the current opinion that sediments underlying OMZs are a strong sink for DIN. Only at greater water depths (300-1000 m) did the sediments become a net sink for DIN. Here, denitrification was the major process (?2 mmol N m⁻² d⁻¹) and removed 55-73% of NO₃⁻ and NO₂⁻ taken up by the sediments, with DNRA and anammox accounting for the remaining fraction. Anammox was of minor importance on the shelf and upper slope yet contributed up to 62% to total N₂ production at the 1000 m station. The results indicate that the partitioning of oxidized N (NO₃⁻, NO₂⁻) into DNRA or denitrification is a key factor determining the role of marine sediments as DIN sinks or recycling sites. Consequently, high measured benthic uptake rates of oxidized N within OMZs do not necessarily indicate a loss of fixed N from the marine environment.     

Geotechnical characterisation of stratocone crater wall sequences, White Island Volcano, New Zealand

Geotechnical characterisation is undertaken for 3 broad units comprising the bulk of the stratigraphy identified on White Island Volcano, Bay of Plenty, New Zealand, an active island stratovolcano. Field and laboratory measurements were used to describe rock mass characteristics for jointed lava flow units, and ring shear tests were undertaken to derive residual strength parameters for joint infilling materials within the lavas. Rock Mass Rating (RMR) and Geological Strength Index (GSI) values were calculated and converted to Mohr-Coulomb strength parameters using the Hoek-Brown criterion. Backanalysis of known landslide scarps was used to derive strength parameters for brecciated rock masses and hydrothermally altered rock masses. Andesite lava flows have high intact strength (σ_ci = 184 ± 50 MN m^− 2; γ = 24.7 ± 0.3 kN m^− 3) and typically 3 wide, infilled joint sets, one parallel to flow direction and two steeply inclined, with spacings of 0.3-1.7 m. Joints are rough, with estimated friction angles for clean joints of ?_j = 42-47°. Joint infill materials are clayey silts derived from weathering of wall rocks and primary volcanic sources; they have low plastic (54%) and liquid (84%) limits and residual strength values of c_r = 0 kN m^− 2 and ?_r = 23.9 ± 3.1°. RMR values range from 70 to 73, giving calculated strength parameters of c′ = 1161-3391 kN m^− 2 and ?′ = 50.5-62.3°. Backanalysis suggests brecciated rock masses have c′ = 0 kN m^− 2 and ?′ = 35.4°, whereas GSI observations in the field suggest higher cohesion (c′ = 306-719 kN m^− 2) and a range of friction angles bracketing the backanalysed result (?′ = 30.6-41.7°). Hydrothermally altered rock masses have c′ = 369 kN m^− 2 and ?′ = 14.9°, indicating considerable loss of strength, especially frictional resistance, compared with the fresh lava units. Values measured at outcrop scale in this study are in keeping with other published values for similar volcanic edifices; backanalysed data suggest weaker rock mass properties than those determined at outcrop. This is interpreted as a scale issue, whereby rock mass characteristics of a large rock mass (crater wall scale) are weaker than those of small outcrops, due in part to the overestimation of friction angle from measurements on small exposures.     

Diffusive reequilibration of quartz-hosted silicate melt and fluid inclusions: Are all metal concentrations unmodified?

Experiments were conducted to determine the extent and mechanism by which the composition of quartz-hosted silicate melt inclusions (SMI) and aqueous fluid inclusions (FI) can undergo post-entrapment modification via diffusion. Quartz slabs containing assemblages of SMI and FI were reacted with synthetic HCl bearing and metalliferous aqueous fluids at T = 500-720 °C and P = 150-200 MPa. SMI from the single inclusion assemblages were analyzed by laser ablation inductively coupled plasma mass spectrometry (LA-ICPMS) and electron probe microanalysis (EPMA) before and after the experiments. Analyses revealed that rapid diffusion of the univalent cations Na⁺, Li⁺, Ag⁺, Cu⁺ and H⁺ occurred through the quartz from the surroundings, resulting in significant changes in the concentrations of these elements in the inclusions. Concentrations of other elements with an effective ionic radius larger than that of Ag⁺, or multiple valence states were not modified in the inclusions during the experiments. Our results warn inclusion‘‘ researchers that the interpretation of Na, Li, Cu and Ag concentrations from quartz-hosted SMI and FI should be treated critically.     

Characterization and cycling of atmospheric mercury along the central US Gulf Coast

Concentrations of atmospheric Hg species, elemental Hg (Hg^°), reactive gaseous Hg (RGM), and fine particulate Hg (Hg-PM_2.5) were measured at a coastal site near Weeks Bay, Alabama from April to August, 2005 and January to May, 2006. Mean concentrations of the species were 1.6 ± 0.3 ng m⁻³, 4.0 ± 7.5 pg m⁻³ and 2.7 ± 3.4 pg m⁻³, respectively. A strong diel pattern was observed for RGM (midday maximum concentrations were up to 92.7 pg m⁻³), but not for Hg^° or Hg-PM_2.5. Elevated RGM concentrations (>25 pg m⁻³) in April and May of 2005 correlated with elevated average daytime O₃ concentrations (>55 ppbv) and high light intensity (>500 W m⁻²). These conditions generally corresponded with mixed continental-Gulf and exclusively continental air mass trajectories. Generally lower, but still elevated, RGM peaks observed in August, 2005 and January–March, 2006 correlated significantly (p < 0.05) with peaks in SO₂ concentration and corresponded to periods of high light intensity and lower average daytime O₃ concentrations. During these times air masses were dominated by trajectories that originated over the continent. Elevated RGM concentrations likely resulted from photochemical oxidation of Hg^° by atmospheric oxidants. This process may have been enhanced in and by the near-shore environment relative to inland sites. The marine boundary layer itself was not found to be a significant source of RGM.     

Joint determination of K decay constants and Ar∗/K for the Fish Canyon sanidine standard, and improved accuracy for Ar/Ar geochronology

⁴⁰Ar/³⁹Ar and K-Ar geochronology have long suffered from large systematic errors arising from imprecise K and Ar isotopic data for standards and imprecisely determined decay constants for the branched decay of ⁴⁰K by electron capture and β⁻ emission. This study presents a statistical optimization approach allowing constraints from ⁴⁰K activity data, K-Ar isotopic data, and pairs of ²³⁸U-²⁰⁶Pb and ⁴⁰Ar/³⁹Ar data for rigorously selected rocks to be used as inputs for estimating the partial decay constants (λ_ε and λ_β) of ⁴⁰K and the ⁴⁰Ar∗/⁴⁰K ratio (κ_FCs) of the widely used Fish Canyon sanidine (FCs) standard. This yields values of κ_FCs = (1.6418 ± 0.0045) × 10⁻³, λ_ε = (0.5755 ± 0.0016) × 10⁻¹⁰ a⁻¹ and λ_β = (4.9737 ± 0.0093) × 10⁻¹⁰ a⁻¹. These results improve uncertainties in the decay constants by a factor of >4 relative to values derived from activity data alone. Uncertainties in these variables determined by our approach are moderately to highly correlated (cov(κ_FCs, λ_ε) = 7.1889 × 10⁻¹⁹, cov(κ_FCs, λ_β) = −7.1390 × 10⁻¹⁹, cov(λ_ε, λ_β) = −3.4497 × 10⁻²⁶) and one must take account of the covariances in error propagation by either linear or Monte Carlo methods. ⁴⁰Ar/³⁹Ar age errors estimated from these results are significantly reduced relative to previous calibrations. Also, age errors are smaller for a comparable level of isotopic measurement precision than those produced by the ²³⁸U/²⁰⁶Pb system, because the ⁴⁰Ar/³⁹Ar system is now jointly calibrated by both the ⁴⁰K and ²³⁸U decay constants, and because λ_ε(⁴⁰K) < λ(²³⁸U). Based on this new calibration, the age of the widely used Fish Canyon sanidine standard is 28.305 ± 0.036 Ma. The increased accuracy of ⁴⁰Ar/³⁹Ar ages is now adequate to provide meaningful validation of high-precision U/Pb or astronomical tuning ages in cases where closed system behavior of K and Ar can be established.     

Solvation processes in steam: Ab initio calculations of ion-solvent structures and clustering equilibria

Reports of the high ion content of steam and low-density supercritical fluids date back to the work of Carlon [Carlon H. R. (1980) Ion content of air humidified by boiling water.J. Appl.Phys.51, 171-173], who invoked ion and neutral-water clustering as mechanism to explain why ions partition into the low-density aqueous phase. Mass spectrometric, vibrational spectroscopic measurements and quantum chemical calculations have refined this concept by proposing strongly bound ion-solvent aggregates and water clusters such as Eigen- and Zundel-type proton clusters H₃O⁺·(H₂O)_m and the more weakly bound water oligomers (H₂O)_m. The extent to which these clusters affect fluid chemistry is determined by their abundance, however, little is known regarding the stability of such moieties in natural low-density high-temperature fluids. Here we report results from quantum chemical calculations using chemical-accuracy multi-level G3 (Curtiss-Pople) and CBS-Q theory (Peterson) to address this question. In particular, we have investigated the cluster structures and clustering equilibria for the ions and H₃S⁺·(H₂O)_m(H₂S)_n, where m ? 6 and n ? 4, at 300-1000 K and 1 bar as well as under vapor-liquid equilibrium conditions between 300 and 646 K. We find that incremental hydration enthalpies and entropies derived from van’t Hoff analyses for the attachment of H₂O and H₂S onto H₃O⁺, and H₃S⁺ are in excellent agreement with experimental values and that the addition of water to all three ions is energetically more favorable than solvation by H₂S. As clusters grow in size, the energetic trends of cluster hydration begin to reflect those for bulk H₂O liquids, i.e. calculated hydration enthalpies and entropies approach values characteristic of the condensation of bulk water (ΔH^o = −44.0 kJ mol⁻¹, ΔS^o = −118.8 J K mol⁻¹). Water and hydrogen sulfide cluster calculations at higher temperatures indicate that a significant fraction of H₃O⁺, and H₃S⁺ ions exists as solvated moieties.     

Streamwater chemistry in three contrasting monolithologic Czech catchments

Hydrochemical patterns resulting from differing bedrock geochemistry were ascertained by concurrent streamwater sampling in three small catchments, each underlain by geochemically contrasting silicate rock types in the western Czech Republic, Central Europe in 2001–2010. The catchments are situated 5–7 km apart in the Slavkov Forest and are occupied by Norway spruce (Picea abies) plantations. They have similar altitude, area, topography, mean annual air temperature, and atmospheric deposition fluxes. The amount of base cations oxides (Ca + Mg + Na + K) is markedly different among the three studied rocks (leucogranite 8%, amphibolite 22%, serpentinite 36%). The leucogranite contains a very small amount of MgO, while the serpentinite contains extremely large amounts of MgO. The amphibolite contains an intermediate amount of MgO and elevated CaO. The Lysina, on leucogranite, exhibited very small concentrations of Mg (median 0.4 mg L⁻¹) in streamwater; Pluh?v Bor, on serpentinite, contained extremely high concentrations of streamwater Mg (18 mg L⁻¹). Streamwater in the Na Zeleném catchment, on amphibolite, contained an intermediate amount of Mg and an elevated Ca. Very low pH (4.2), negative alkalinity (−60 μeq L⁻¹) and high inorganic monomeric Al concentrations (0.3 mg L⁻¹) were found in the stream draining leucogranite. Serpentinite streamwater exhibited the highest pH (7.6), alkalinity (+940 μeq L⁻¹) and Ni concentrations (100 μg L⁻¹). Aquatic chemistry reflected the composition of the underlying rocks within the studied catchments. Contrasting streamwater compositions of the studied catchments were generated according to the MAGIC model simulations mainly by differences in chemical weathering rates of base cations (65 meq m⁻² a⁻¹ at Lysina, 198 meq m⁻² a⁻¹ at Na Zeleném, and 241 meq m⁻² a⁻¹ at Pluh?v Bor).     

Determination of fluid/melt partition coefficients by LA-ICPMS analysis of co-existing fluid and silicate melt inclusions: Controls on element partitioning

Analyses of co-existing silicate melt and fluid inclusions, entrapped in quartz crystals in volatile saturated magmatic systems, allowed direct quantitative determination of fluid/melt partition coefficients. Investigations of various granitic systems (peralkaline to peraluminous in composition, log fO₂ = NNO−1.7 to NNO+4.5) exsolving fluids with various chlorinities (1-14 mol/kg) allowed us to assess the effect of these variables on the fluid/melt partition coefficients (D). Partition coefficients for Pb, Zn, Ag and Fe show a nearly linear increase with the chlorinity of these fluid (D_Pb ∼ 6 ∗ m_Cl, D_Zn ∼ 8 ∗ m_Cl, D_Ag ∼ 4 ∗ m_Cl, D_Fe ∼ 1.4 ∗ m_Cl, where m_Cl is the molinity of Cl). This suggests that these metals are dissolved primarily as Cl-complexes and neither oxygen fugacity nor the composition of the melt affects significantly their fluid/melt partitioning. By contrast, partition coefficients for Mo, B, As, Sb and Bi are highest in low salinity (1-2 mol/kg Cl) fluids with maximum values of D_Mo ∼ 20, D_B ∼ 15, D_As ∼ 13, D_Sb ∼ 8, D_Bi ∼ 15 indicating dissolution as non-chloride (e.g., hydroxy) complexes. Fluid/melt partition coefficients of copper are highly variable, but highest between vapor like fluids and silicate melt (D_Cu ? 2700), indicating an important role for ligands other than Cl. Partition coefficients for W generally increase with increasing chlorinity, but are exceptionally low in some of the studied brines which may indicate an effect of other parameters. Fluid/melt partition coefficients of Sn show a high variability but likely increase with the chlorinity of the fluid (D_Sn = 0.3-42, D_W = 0.8-60), and decrease with decreasing oxygen fugacity or melt peraluminosity.