The solubility of rhenium in silicate melts: Implications for the geochemical properties of rhenium at high temperatures

The solubility of rhenium (Re) in a haplobasaltic melt (anorthite-diopside eutectic composition) has been experimentally determined using the mechanically assisted equilibration technique at 1400°C as a function of oxygen fugacity (10−12 < fO2 ≤ 10−7 bar), imposed by CO-CO2 gas mixtures. Samples were analysed by laser ablation-inductively coupled plasma-mass spectrometry (LA-ICP-MS). This is a true microanalytical technique, which allows small-scale sample heterogeneity to be detected, while providing a limit of detection of 2 ppb Re. Time-resolved LA-ICP-MS spectra revealed the presence of suboptically sized micronuggets of Re in all samples, which, because they are present at the 0.5 to 10 ppm level, dominate the true solubilities of Re (<1 ppm at the conditions of the experiment) in bulk analyses of the samples. Nevertheless, the micronuggets could be filtered out from the time-resolved spectra to reveal accurate values of the true Re solubility. A number of time series of samples were taken at constant fO2 to demonstrate that the solubilities converge to a constant value. In addition, solubilities were measured after increasing and decreasing the imposed fO2. The results show that Re dissolves in the silicate melt as ReO2 (Re4+) and ReO3 (Re6+) species, with the latter predominating at typical terrestrial upper-mantle oxygen fugacities. The total solubility of Re is described by the following expression (fO2 in bars): [Re/ppb] = 9.7(±1.9) × 109 (fO2) + 4.2 (±0.3) × 1014 (fO2)1.5Assuming an activity coefficient for Re in Fe-rich metal of 1, this gives a value of DRemet/sil of 5 × 1010 at log fO2 = IW-2, appropriate for metal-silicate partitioning in an homogenously accreting Earth. Thus, Re is indeed very highly siderophile, and the mantle’s abundance cannot be explained by homogenous accretion

Viscoelasticity of crystal- and bubble-bearing rhyolite melts

The effect of non-deformable inclusions on the frequency-dependent rheology of a rhyolite melt plus crystals has been investigated using a sinusoidal torsion deformation device for measurements of shear viscosity and modulus in the frequency range of 5 mHz to 20 Hz at temperatures of 750–1050°C. The relaxed shear viscosity and unrelaxed shear modulus of rhyolite magma (rhyolite melt plus crystals plus bubbles) decreases with increasing bubble content and increases with the addition of crystals. At a crystal concentration of about 45% a relaxed value of the shear viscosity is not attainable. The presence of rigid inclusions results in an imaginary component of the shear modulus that becomes more symmetrical and shifted to the low-frequency—high-temperature range with respect to that for a crystal-free melt. The slope of log(Q−1) (internal friction) as a function of the dimensionless variable log(ωτ), is unaffected in the low-temperature—high-frequency range of crystals, with Q−1 ≈ 1/(ωτ)0.5 (the same as for bubble- and crystal-free rhyolite). For the present type of suspension, the internal friction is practically constant and independent of log(ωτ) in the high-temperature—low-frequency limit (ωτ 1). The shape of the Cole-Cole diagram becomes symmetrical and can be described as a Caputo body with parameter γ ≈ 0.45, whereas for bubble-bearing and inclusion-free rhyolite melts the shape of diagram relates to the β-relaxation exponent with β ≈ 0.5. The present work demonstrates that magma may or may not follow a power-law rheology depending on the relative volume proportion between crystals and bubbles

Permeability of compacting porous lavas

The highly transient nature of outgassing commonly observed at volcanoes is in part controlled by the permeability of lava domes and shallow conduits. Lava domes generally consist of a porous outer carapace surrounding a denser lava core with internal shear zones of variable porosity. Here we examine densification using uniaxial compression experiments on variably crystalline and porous rhyolitic dome lavas from the Taupo Volcanic Zone. Experiments were conducted at 900°C and an applied stress of 3MPa to 60% strain, while monitoring acoustic emissions to track cracking. The evolution of the porous network was assessed via X-ray computed tomography, He-pycnometry, and relative gas permeability. High starting connected porosities led to low apparent viscosities and high strain rates, initially accompanied by abundant acoustic emissions. As compaction ensued, the lavas evolved; apparent viscosity increased and strain rate decreased due to strain hardening of the suspensions. Permeability fluctuations resulted from the interplay between viscous flow and brittle failure. Where phenocrysts were abundant, cracks had limited spatial extent, and pore closure decreased axial and radial permeability proportionally, maintaining the initial anisotropy. In crystal-poor lavas, axial cracks had a more profound effect, and permeability anisotropy switched to favor axial flow. Irrespective of porosity, both crystalline samples compacted to a threshold minimum porosity of 17–19%, whereas the crystal-poor sample did not achieve its compaction limit. This indicates that unconfined loading of porous dome lavas does not necessarily form an impermeable plug and may be hindered, in part by the presence of crystals

Vesiculation and Quenching During Surtseyan Eruptions at Hunga Tonga-Hunga Ha'apai Volcano, Tonga

Surtseyan eruptions are shallow to emergent subaqueous explosive eruptions that owe much of their characteristic behavior to the interaction of magma with water. The difference in thermal properties between water and air affects the cooling and postfragmentation vesiculation processes in magma erupted into the water column. Here we study the vesiculation and cooling processes during the 2009 and 2014–2015 Surtseyan eruptions of Hunga Tonga‐Hunga Ha'apai volcano by combining 2‐D and 3‐D vesicle‐scale analyses of lapilli and bombs and numerical thermal modeling. Most of the lapilli and bombs show gradual textural variations from rim to core. The vesicle connectivity in the lapilli and bombs increases with vesicularity from fully isolated to completely connected and also increases from rim to core in transitional clasts. We interpret the gradual textural variations and the connectivity‐vesicularity relationships as the result of postfragmentation bubble growth and coalescence interrupted at different stages by quenching in water. The measured vesicle size distributions are bimodal with a population of small and large vesicles. We interpret this bimodality as the result of two nucleation events, one prefragmentation with the nucleation and growth of large bubbles and one postfragmentation with nucleation of small vesicles. We link the thermal model with the textural variations in the clasts—showing a dependence on particle size, Leidenfrost effect, and initial melt temperature. In particular, the cooling profiles in the bombs are consistent with the gradual textural variations from rim to core in the clasts, likely caused by variations in time available for vesiculation before quenching

Lava flow rheology: A comparison of morphological and petrological methods

In planetary sciences, the emplacement of lava flows is commonly modelled using a single rheological parameter (apparent viscosity or apparent yield strength) calculated from morphological dimensions using Jeffreys’ and Hulme’s equations. The rheological parameter is then typically further interpreted in terms of the nature and chemical composition of the lava (e.g., mafic or felsic). Without the possibility of direct sampling of the erupted material, the validity of this approach has remained largely untested. In modern volcanology, the complex rheological behaviour of lavas is measured and modelled as a function of chemical composition of the liquid phase, fractions of crystals and bubbles, temperature and strain rate. Here, we test the planetary approach using a terrestrial basaltic lava flow from the Western Volcanic Zone in Iceland. The geometric parameters required to employ Jeffreys’ and Hulme’s equations are accurately estimated from high-resolution HRSC-AX Digital Elevation Models. Samples collected along the lava flow are used to constrain a detailed model of the transient rheology as a function of cooling, crystallisation, and compositional evolution of the residual melt during emplacement. We observe that the viscosity derived from the morphology corresponds to the value estimated when significant crystallisation inhibits viscous deformation, causing the flow to halt. As a consequence, the inferred viscosity is highly dependent on the details of the crystallisation sequence and crystal shapes, and as such, is neither uniquely nor simply related to the bulk chemical composition of the erupted material. This conclusion, drawn for a mafic lava flow where crystallisation is the primary process responsible for the increase of the viscosity during emplacement, should apply to most of martian, lunar, or mercurian volcanic landforms, which are dominated by basaltic compositions. However, it may not apply to felsic lavas where vitrification resulting from degassing and cooling may ultimately cause lava flows to halt