Oxide nanolitisation-induced melt iron extraction causes viscosity jumps and enhanced explosivity in silicic magma.

Explosivity in erupting volcanoes is controlled by the degassing dynamics and the viscosity of the ascending magma in the conduit. Magma crystallisation enhances both heterogeneous bubble nucleation and increases in magma bulk viscosity. Nanolite crystallisation has been suggested to enhance such processes too, but in a noticeably higher extent. Yet the precise causes of the resultant strong viscosity increase remain unclear. Here we report experimental results for rapid nanolite crystallisation in natural silicic magma and the extent of the subsequent viscosity increase. Nanolite-free and nanolite-bearing rhyolite magmas were subjected to heat treatments, where magmas crystallised or re-crystallised oxide nanolites depending on their initial state, showing an increase of one order of magnitude as oxide nanolites formed. We thus demonstrate that oxide nanolites crystallisation increases magma bulk viscosity mainly by increasing the viscosity of its melt phase due to the chemical extraction of iron, whereas the physical effect of particle suspension is minor, almost negligible. Importantly, we further observe that this increase is sufficient for driving magma fragmentation depending on magma degassing and ascent dynamics.

Effect of water on the glass transition of a potassium-magnesium carbonate melt.

Calorimetric measurements of the glass transition temperatures (Tg) of hydrous carbonate melts are reported on a near-eutectic composition of 55 mol% K2CO3 - 45 mol% MgCO3 with up to 42 mol% bulk H2O dissolved in the carbonate melt. Hydrous melts were quenched from 750°C to transparent and crystal-free glasses and were subsequently analysed for water content before and after measuring Tg by high-sensitivity differential scanning calorimetry. The glass transition and limited fictive temperatures as a function of the water content were determined at 10 K/min cooling/heating rates resulting in Tg ranging from 245°C at nominally anhydrous conditions to 83°C in the presence of 42 mol% H2O in the glass. Through a generalized Gordon-Taylor analysis, the factors k (7.27), k0 (3.2) and the interaction parameter Ax (0.49) were derived. The limited fictive temperature of a hypothetically, zero water containing 55 mol% K2CO3 - 45 mol% MgCO3 glass is 232 ± 5°C (505 K). The high value of the interaction parameter A indicates strong specific molecular interactions between water and the carbonates in the glassy state and a large decrease in the excess enthalpy of mixing during the conversion of the glassy into the liquid state at the glass transition. This article is part of the theme issue 'Exploring the length scales, timescales and chemistry of challenging materials (Part 1)'.

Glass transition temperatures and crystallization kinetics of a synthetic, anhydrous, amorphous calcium-magnesium carbonate.

We report the first calorimetric observations of glass transition temperatures and crystallization rates of anhydrous, amorphous calcium-magnesium carbonate using fast scanning differential scanning calorimetry. Hydrous amorphous Ca0.95Mg0.05CO3 · 0.5H2O (ACMC) solid was precipitated from a MgCl2-NaHCO3 buffered solution, separated from the supernatant, and freeze-dried. An aliquot of the freeze-dried samples was additionally dried at 250°C for up to 6 h in a furnace and in a high-purity N2 atmosphere to produce anhydrous ACMC. The glass transition temperature of the anhydrous Ca0.95Mg0.05CO3 was determined by applying different heating rates (1000-6000 K s-1) and correcting for thermal lag to be 376°C and the relaxational heat capacity was determined to be Cp = 0.16 J/(g K). Additionally, the heating rate dependence of the temperature that is associated with the corrected crystallization peaks is used to determine the activation energy of crystallization to be 275 kJ mol-1. A high-resolution transmission electron microscopy study on the hydrous and anhydrous samples provided further constraints on their compositional and structural states. This article is part of the theme issue 'Exploring the length scales, timescales and chemistry of challenging materials (Part 1)'.

Experimental generation of fulgurite under realistic lightning discharge conditions.

Fulgurites have been documented in geological deposits from throughout Earth's history. They have also been assigned a potential role in prebiotic chemistry as a source of reactants. Fulgurites are generated in nature by cloud-to-ground lightning strikes. The unpredictability in space and time of the occurrence of lightning events has limited the investigation of both the mechanisms by, and the conditions under, which fulgurites form. A laboratory-based approach can mitigate these limitations. Here, we describe experimentally generated fulgurites generated from Laacher See volcanic ash. We employ a DC source with a trigger-pulse setup in a high voltage laboratory, whose capabilities enable experimental conditions that correspond closely to the electrical characteristics of natural lightning strikes. The experimentally generated fulgurites closely resemble naturally-occurring fulgurites in both state and texture. These experimental investigations yield a high reproducibility of the characteristic of fulgurites generated under well-constrained conditions, enabling some inferences to be made regarding the processes involved in the generation of fulgurites in nature. This work provides a basis for a systematic characterization of experimental fulgurites and the characteristic of lightning discharges.

Bubble rise in molten glasses and silicate melts during heating and cooling cycles.

The Hadamard-Rybczynski equation describes the steady-state buoyant rise velocity of an unconfined spherical bubble in a viscous liquid. This solution has been experimentally validated for the case where the liquid viscosity is held constant. Here, we extend this result for non-isothermal conditions, by developing a solution for bubble position in which we account for the time-dependent liquid viscosity, liquid and gas densities, and bubble radius. We validate this solution using experiments in which spherical bubbles are created in a molten silicate liquid by cutting gas cavities into glass sheets, which are stacked, then heated through the glass transition interval. The bubble-bearing liquid, which has a strongly temperature-dependent viscosity, is subjected to various heating and cooling programs such that the bubble rise velocity varies through the experiment. We find that our predictions match the final observed position of the bubble measured in blocks of cooled glass to within the experimental uncertainty, even after the application of a complex temperature-time pathway. We explore applications of this solution for industrial, artistic, and natural volcanological applied problems.

Combined effusive-explosive silicic volcanism straddles the multiphase viscous-to-brittle transition.

Magma is a viscoelastic fluid that can support fracture propagation when local shear stresses are high, or relax and flow when shear stresses are low. Here we present experiments to confirm this using synthetic and natural magmatic liquids across eruptive conditions and use Maxwell's linear viscoelasticity to parameterize our results and predict the maximum stresses that can be supported during flow. This model proves universal across a large range of liquid compositions, temperatures, crystallinity and rates of strain relevant to shallow crustal magma ascent. Our results predict that the 2008 Volcán Chaitén eruption resided in the viscous field at the onset of magma ascent, but transitioned to a mixed viscous-brittle regime during degassing, coincident with the observed combined effusive-explosive behaviour during dome extrusion. Taking a realistic maximum effusive ascent rate, we propose that silicic eruptions on Earth may straddle the viscous-to-brittle transition by the time they reach the surface.

Aggregation in particle rich environments: a textural study of examples from volcanic eruptions, meteorite impacts, and fluidized bed processing.

Aggregation is a common process occurring in many diverse particulate gas mixtures (e.g. those derived from explosive volcanic eruptions, meteorite impact events, and fluid bed processing). It results from the collision and sticking of particles suspended in turbulent gas/air. To date, there is no generalized model of the underlying physical processes. Here, we investigate aggregates from 18 natural deposits (16 volcanic deposits and two meteorite impact deposits) as well as aggregates produced experimentally via fluidized bed techniques. All aggregates were analyzed for their size, internal structuring, and constituent particle size distribution. Commonalities and differences between the aggregate types are then used to infer salient features of the aggregation process. Average core to rim ratios of internally structured aggregates (accretionary lapilli) is found to be similar for artificial and volcanic aggregates but up to an order of magnitude different than impact-related aggregates. Rim structures of artificial and volcanic aggregates appear to be physically similar (single, sub-spherical, regularly-shaped rims) whereas impact-related aggregates more often show multiple or irregularly shaped rims. The particle size distributions (PSDs) of all three aggregate types are similar (< 200 µm). This proves that in all three environments, aggregation occurs under broadly similar conditions despite the significant differences in source conditions (particle volume fraction, particle size distribution, particle composition, temperature), residence times, plume conditions (e.g., humidity and temperature), and dynamics of fallout and deposition. Impact-generated and volcanic aggregates share many similarities, and in some cases may be indistinguishable without their stratigraphic context.

Enhancement of eruption explosivity by heterogeneous bubble nucleation triggered by magma mingling.

We present new evidence that shows magma mingling can be a key process during highly explosive eruptions. Using fractal analysis of the size distribution of trachybasaltic fragments found on the inner walls of bubbles in trachytic pumices, we show that the more mafic component underwent fracturing during quenching against the trachyte. We propose a new mechanism for how this magmatic interaction at depth triggered rapid heterogeneous bubble nucleation and growth and could have enhanced eruption explosivity. We argue that the data support a further, and hitherto unreported contribution of magma mingling to highly explosive eruptions. This has implications for hazard assessment for those volcanoes in which evidence of magma mingling exists.

Wetting and Spreading of Molten Volcanic Ash in Jet Engines.

A major hazard to jet engines posed by volcanic ash is linked to the wetting and spreading of molten ash droplets on engine component surfaces. Here, using the sessile drop method, we study the evolution of the wettability and spreading of volcanic ash. We employ rapid temperature changes up to 1040-1450 °C, to replicate the heating conditions experienced by volcanic ash entering an operating jet engine. In this scenario, samples densify as particles coalesce under surface tension until they form a large system-sized droplet (containing remnant gas bubbles and crystals), which subsequently spreads on the surface. The data exhibit a transition from a heterogeneous to a homogeneous wetting regime above 1315 °C as crystals in the drops are dissolved in the melt. We infer that both viscosity and microstructural evolution are key controls on the attainment of equilibrium in the wetting of molten volcanic ash droplets.

Local geology controlled the feasibility of vitrifying Iron Age buildings.

During European prehistory, hilltop enclosures made from polydisperse particle-and-block stone walling were exposed to temperatures sufficient to partially melt the constituent stonework, leading to the preservation of glassy walls called 'vitrified forts'. During vitrification, the granular wall rocks partially melt, sinter viscously and densify, reducing inter-particle porosity. This process is strongly dependent on the solidus temperature, the particle sizes, the temperature-dependence of the viscosity of the evolving liquid phase, as well as the distribution and longevity of heat. Examination of the sintering behaviour of 45 European examples reveals that it is the raw building material that governs the vitrification efficiency. As Iron Age forts were commonly constructed from local stone, we conclude that local geology directly influenced the degree to which buildings were vitrified in the Iron Age. Additionally, we find that vitrification is accompanied by a bulk material strengthening of the aggregates of small sizes, and a partial weakening of larger blocks. We discuss these findings in the context of the debate surrounding the motive of the wall-builders. We conclude that if wall stability by bulk strengthening was the desired effect, then vitrification represents an Iron Age technology that failed to be effective in regions of refractory local geology.

Topological inversions in coalescing granular media control fluid-flow regimes.

Sintering-or coalescence-of viscous droplets is an essential process in many natural and industrial scenarios. Current physical models of the dynamics of sintering are limited by the lack of an explicit account of the evolution of microstructural geometry. Here, we use high-speed time-resolved x-ray tomography to image the evolving geometry of a sintering system of viscous droplets, and use lattice Boltzmann simulations of creeping fluid flow through the reconstructed pore space to determine its permeability. We identify and characterize a topological inversion, from spherical droplets in a continuous interstitial gas, to isolated bubbles in a continuous liquid. We find that the topological inversion is associated with a transition in permeability-porosity behavior, from Stokes permeability at high porosity, to percolation theory at low porosity. We use these findings to construct a unified physical description that reconciles previously incompatible models for the evolution of porosity and permeability during sintering.

Raman spectra of Martian glass analogues: A tool to approximate their chemical composition.

Raman spectrometers will form a key component of the analytical suite of future planetary rovers intended to investigate geological processes on Mars. In order to expand the applicability of these spectrometers and use them as analytical tools for the investigation of silicate glasses, a database correlating Raman spectra to glass composition is crucial. Here we investigate the effect of the chemical composition of reduced silicate glasses on their Raman spectra. A range of compositions was generated in a diffusion experiment between two distinct, iron-rich end-members (a basalt and a peralkaline rhyolite), which are representative of the anticipated compositions of Martian rocks. Our results show that for silica-poor (depolymerized) compositions the band intensity increases dramatically in the regions between 550-780 cm-1 and 820-980 cm-1. On the other hand, Raman spectra regions between 250-550 cm-1 and 1000-1250 cm-1 are well developed in silica-rich (highly polymerized) systems. Further, spectral intensity increases at ~965 cm-1 related to the high iron content of these glasses (~7-17 wt % of FeOtot). Based on the acquired Raman spectra and an ideal mixing equation between the two end-members we present an empirical parameterization that enables the estimation of the chemical compositions of silicate glasses within this range. The model is validated using external samples for which chemical composition and Raman spectra were characterized independently. Applications of this model range from microanalysis of dry and hydrous silicate glasses (e.g., melt inclusions) to in situ field investigations and studies under extreme conditions such as extraterrestrial (i.e., Mars) and submarine volcanic environments.

Volcanic ash melting under conditions relevant to ash turbine interactions.

The ingestion of volcanic ash by jet engines is widely recognized as a potentially fatal hazard for aircraft operation. The high temperatures (1,200-2,000 °C) typical of jet engines exacerbate the impact of ash by provoking its melting and sticking to turbine parts. Estimation of this potential hazard is complicated by the fact that chemical composition, which affects the temperature at which volcanic ash becomes liquid, can vary widely amongst volcanoes. Here, based on experiments, we parameterize ash behaviour and develop a model to predict melting and sticking conditions for its global compositional range. The results of our experiments confirm that the common use of sand or dust proxy is wholly inadequate for the prediction of the behaviour of volcanic ash, leading to overestimates of sticking temperature and thus severe underestimates of the thermal hazard. Our model can be used to assess the deposition probability of volcanic ash in jet engines.

Magma mixing induced by particle settling.

A time series of experiments at high temperature have been performed to investigate the influence of particle settling on magma mixing. A natural rhyolite glass was held above a natural basalt glass in a platinum crucible. After melting of the glasses at superliquidus temperatures, a platinum sphere was placed on the upper surface of the rhyolitic melt and sank into the experimental column (rhyolitic melt above basaltic melt). Upon falling through the rhyolitic-basaltic melt interface, the Pt sphere entrained a filament of rhyolitic melt in its further fall. The quenched products of the experiments were imaged using X-ray microCT methods. The images of our time series of experiments document the formation of a rhyolite filament as it is entrained into the underlying basalt by the falling platinum sphere. When the Pt particle reached the bottom of the crucible, the entrained rhyolitic filament started to ascend buoyantly up to the initial rhyolitic-basaltic interface. This generated a significant thickness increase of a comingled "melange" layer at the interface due to "liquid rope coiling" and piling up of the filament. As a consequence, the basalt/rhyolite interface was greatly enlarged and diffusive hybridisation greatly accelerated. Further, bubbles, originating at the interface, are observed to have risen into the overlying rhyolite dragging basalt filaments with them. Upon crossing the basalt/rhyolite interface, the bubbles have non-spherical shapes as they adapt to the differing surface tensions of basaltic and rhyolitic melts. Major element profiles, measured across the rhyolite filaments, exhibit asymmetrical shapes from the rhyolite into the basalt. Na and Ti reveal uphill diffusion from the rhyolite towards the interface in the filament cross sections. These results reveal the potential qualitative complexity of the mingling process between rhyolitic and basaltic magmas in the presence of sinking crystals. They imply that crystal-rich magma mingling may be expected to be accelerated with respect to crystal-poor systems. We urge the further fluid dynamic analysis of these phenomena, obtainable for the first time using detailed tomographic imaging.

Fault rheology beyond frictional melting.

During earthquakes, comminution and frictional heating both contribute to the dissipation of stored energy. With sufficient dissipative heating, melting processes can ensue, yielding the production of frictional melts or "pseudotachylytes." It is commonly assumed that the Newtonian viscosities of such melts control subsequent fault slip resistance. Rock melts, however, are viscoelastic bodies, and, at high strain rates, they exhibit evidence of a glass transition. Here, we present the results of high-velocity friction experiments on a well-characterized melt that demonstrate how slip in melt-bearing faults can be governed by brittle fragmentation phenomena encountered at the glass transition. Slip analysis using models that incorporate viscoelastic responses indicates that even in the presence of melt, slip persists in the solid state until sufficient heat is generated to reduce the viscosity and allow remobilization in the liquid state. Where a rock is present next to the melt, we note that wear of the crystalline wall rock by liquid fragmentation and agglutination also contributes to the brittle component of these experimentally generated pseudotachylytes. We conclude that in the case of pseudotachylyte generation during an earthquake, slip even beyond the onset of frictional melting is not controlled merely by viscosity but rather by an interplay of viscoelastic forces around the glass transition, which involves a response in the brittle/solid regime of these rock melts. We warn of the inadequacy of simple Newtonian viscous analyses and call for the application of more realistic rheological interpretation of pseudotachylyte-bearing fault systems in the evaluation and prediction of their slip dynamics.

Approximate chemical analysis of volcanic glasses using Raman spectroscopy.

The effect of chemical composition on the Raman spectra of a series of natural calcalkaline silicate glasses has been quantified by performing electron microprobe analyses and obtaining Raman spectra on glassy filaments (~450 µm) derived from a magma mingling experiment. The results provide a robust compositionally-dependent database for the Raman spectra of natural silicate glasses along the calcalkaline series. An empirical model based on both the acquired Raman spectra and an ideal mixing equation between calcalkaline basaltic and rhyolitic end-members is constructed enabling the estimation of the chemical composition and degree of polymerization of silicate glasses using Raman spectra. The model is relatively insensitive to acquisition conditions and has been validated using the MPI-DING geochemical standard glasses1 as well as further samples. The methods and model developed here offer several advantages compared with other analytical and spectroscopic methods such as infrared spectroscopy, X-ray fluorescence spectroscopy, electron and ion microprobe analyses, inasmuch as Raman spectroscopy can be performed with a high spatial resolution (1 µm2) without the need for any sample preparation as a nondestructive technique. This study represents an advance in efforts to provide the first database of Raman spectra for natural silicate glasses and yields a new approach for the treatment of Raman spectra, which allows us to extract approximate information about the chemical composition of natural silicate glasses using Raman spectroscopy. We anticipate its application in handheld in situ terrestrial field studies of silicate glasses under extreme conditions (e.g. extraterrestrial and submarine environments).

Volcanic sintering: Timescales of viscous densification and strength recovery.

[1] Sintering and densification are ubiquitous processes influencing the emplacement of both effusive and explosive products of volcanic eruptions. Here we sinter ash-size fragments of a synthetic National Institute of Standards and Technology viscosity standard glass at temperatures at which the resultant melt has a viscosity of ∼108-109 Pa.s at 1bar to assess sintering dynamics under near-surface volcanic conditions. We track the strength recovery via uniaxial compressive tests. We observe that volcanic ash sintering is dominantly time dependent, temperature dependent, and grain size dependent and may thus be interpreted to be controlled by melt viscosity and surface tension. Sintering evolves from particle agglutination to viscous pore collapse and is accompanied by a reduction in connected porosity and an increase in isolated pores. Sintering and densification result in a nonlinear increase in strength. Micromechanical modeling shows that the pore-emanated crack model explains the strength of porous lava as a function of pore fraction and size.

Synthesis, crystal structure, and properties of two modifications of MgB(12)C(2).

Single crystals of two modifications of the new magnesium boride carbide MgB(12)C(2) were synthesized from the elements in a metallic melt by using tantalum ampoules. Crystals were characterized by single-crystal X-ray diffraction and electron microprobe analysis (energy-dispersive (EDX) and wavelength-dispersive (WDX) X-ray spectroscopy). Orthorhombic MgB(12)C(2) is formed in a Cu/Mg melt at 1873 K. The crystal structure of o-MgB(12)C(2) (Imma, Z=4, a=5.6133(10), b=9.828(2), c=7.9329(15) A, 574 reflections, 42 variables, R(1)(F)=0.0208, wR(2)(I)=0.0540) consists of a hexagonal primitive array of B(12) icosahedra with Mg atoms and C(2) units in trigonal-prismatic voids. Each icosahedron has six exohedral B--B and six B--C bonds. Carbon is tetrahedrally coordinated by three boron atoms and one carbon atom with a remarkably long C--C distance of 1.727 A. Monoclinic MgB(12)C(2) is formed in an Al/Mg melt at 1573 K. The structure of m-MgB(12)C(2) (C2/c, Z=4, a=7.2736(11), b=8.7768(13), c=7.2817(11) A, beta=105.33(3) degrees , 1585 reflections, 71 variables, R(1)(F)=0.0228, wR(2)(I)=0.0610) may be described as a distorted cubic close arrangement of B(12) icosahedra. Tetrahedral voids are filled by C atoms and octahedral voids are occupied by Mg atoms. The icosahedra are interconnected by four exohedral B--B bonds to linear chains and by eight interstitial C atoms to form a three-dimensional covalent network. Both compounds fulfill the electron-counting rules of Wade and Longuet-Higgins.