Damage amplification during repetitive seismic waves in mechanically loaded rocks.

Cycles of stress build-up and release are inherent to tectonically active planets. Such stress oscillations impart strain and damage, prompting mechanically loaded rocks and materials to fail. Here, we investigate, under uniaxial conditions, damage accumulation and weakening caused by time-dependent creep (at 60, 65, and 70% of the rocks' expected failure stress) and repeating stress oscillations (of ± 2.5, 5.0 or 7.5% of the creep load), simulating earthquakes at a shaking frequency of ~ 1.3 Hz in volcanic rocks. The results show that stress oscillations impart more damage than constant loads, occasionally prompting sample failure. The magnitudes of the creep stresses and stress oscillations correlate with the mechanical responses of our porphyritic andesites, implicating progressive microcracking as the cause of permanent inelastic strain. Microstructural investigation reveals longer fractures and higher fracture density in the post-experimental rock. We deconvolve the inelastic strain signal caused by creep deformation to quantify the amount of damage imparted by each individual oscillation event, showing that the magnitude of strain is generally largest with the first few oscillations; in instances where pre-existing damage and/or the oscillations' amplitude favour the coalescence of micro-cracks towards system scale failure, the strain signal recorded shows a sharp increase as the number of oscillations increases, regardless of the creep condition. We conclude that repetitive stress oscillations during earthquakes can amplify the amount of damage in otherwise mechanically loaded materials, thus accentuating their weakening, a process that may affect natural or engineered structures. We specifically discuss volcanic scenarios without wholesale failure, where stress oscillations may generate damage, which could, for example, alter pore fluid pathways, modify stress distribution and affect future vulnerability to rupture and associated hazards.

Quantification of ash sedimentation dynamics through depolarisation imaging with AshCam.

Even modest ash-rich volcanic eruptions can severely impact a range of human activities, especially air travel. The dispersal of ash in these eruptions depends critically on aggregation and sedimentation processes - however these are difficult to quantify in volcanic plumes. Here, we image ash dynamics from mild explosive activity at Santiaguito Volcano, Guatemala, by measuring the depolarisation of scattered sunlight by non-spherical ash particles, allowing the dynamics of diffuse ash plumes to be investigated with high temporal resolution (>1 Hz). We measure the ash settling velocity downwind from the main plume, and compare it directly with ground sampled ash particles, finding good agreement with a sedimentation model based on particle size. Our new, cost-effective technique leverages existing technology, opening a new frontier of integrated ash visualisation and ground collection studies which could test models of ash coagulation and sedimentation, leading to improved ash dispersion forecasts. This will provide risk managers with improved data quality on ash location, reducing the economic and societal impacts of future ash-rich eruptions.

Facile production of nanocomposites of carbon nanotubes and polycaprolactone with high aspect ratios with potential applications in drug delivery.

The geometries and surface properties of nanocarriers greatly influence the interaction between nanomaterials and living cells. In this work we combine multiwalled carbon nanotubes (CNTs) with poly-ε-caprolactone (PCL) to produce non-spherical nanocomposites with high aspect ratios by using a facile emulsion solvent evaporation method. Particles were characterised by dynamic light scattering (DLS), scanning electron microscopy (SEM), atomic force microscopy (AFM) and asymmetric flow field flow fractionation (AF4). Different sizes and morphologies of nanoparticles were produced depending on the concentration of the sodium dodecyl sulphate (SDS), CNTs and PCL. Rod-like PCL-CNT nanostructures with low polydispersity were obtained with 1.5 mg mL-1 of SDS, 0.9 mg mL-1 of CNTs and 10 mg mL-1 PCL. AFM analysis revealed that the PCL and PCL-CNT nanocomposite had comparatively similar moduli of 770 and 560 MPa respectively, indicating that all the CNTs have been coated with at least 2 nm of PCL. Thermogravimetric analysis of the PCL-CNT nanocomposite indicated that they contained 9.6% CNTs by mass. The asymmetric flow field flow fractionation of the samples revealed that the PCL-CNT had larger hydrodynamic diameters than PCL alone. Finally, the drug loading properties of the nanocomposites were assessed using docetaxel as the active substance. The nanocomposites showed comparable entrapment efficiencies of docetaxel (89%) to the CNTs alone (95%) and the PCL nanoparticles alone (81%). This is a facile method for obtaining non-spherical nanocomposites that combines the properties of PCL and CNTs such as the high aspect ratio, modulus. The high drug entrapment efficiency of these nanocomposites may have promising applications in drug delivery.

Disclosing the temperature of columnar jointing in lavas.

Columnar joints form by cracking during cooling-induced contraction of lava, allowing hydrothermal fluid circulation. A lack of direct observations of their formation has led to ambiguity about the temperature window of jointing and its impact on fluid flow. Here we develop a novel thermo-mechanical experiment to disclose the temperature of columnar jointing in lavas. Using basalts from Eyjafjallajökull volcano (Iceland) we show that contraction during cooling induces stress build-up below the solidus temperature (980 °C), resulting in localised macroscopic failure between 890 and 840 °C. This temperature window for incipient columnar jointing is supported by modelling informed by mechanical testing and thermal expansivity measurements. We demonstrate that columnar jointing takes place well within the solid state of volcanic rocks, and is followed by a nonlinear increase in system permeability of <9 orders of magnitude during cooling. Columnar jointing may promote advective cooling in magmatic-hydrothermal environments and fluid loss during geothermal drilling and thermal stimulation.

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.

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.

Seismic and experimental insights into eruption precursors at Volcán de Colima.

We combine geophysical and experimental observations to interpret preeruptive unrest at Volcán de Colima in 1998. 17,893 volcanic earthquakes were detected between 1 October and 31 December 1998, including 504 clusters. Using seismic ambient noise interferometry, we observe a drop in velocity prior to the eruption linked to damage accumulation during magma ascent. This is supported by experimental observations where static stress causes a velocity decrease prior to failure. Furthermore, we observe acoustic emission clusters during the experiments, with lower porosity samples producing higher numbers of repeaters. This behavior introduces tensile failure as an additional viable mechanism for clusters during magma ascent. The findings suggest that preeruptive magma ascent may be monitored to variable degrees of accuracy via descriptions of damage accumulation and associated seismic velocity changes.

Conduit margin heating and deformation during the AD 1886 basaltic Plinian eruption at Tarawera volcano, New Zealand.

During explosive eruptions, a suspension of gas and pyroclasts rises rapidly within a conduit. Here, we have analysed textures preserved in the walls of a pyroclastic feeder dyke of the AD 1886 Tarawera basaltic Plinian fissure eruption. The samples examined consist of basaltic ash and scoria plastered onto a conduit wall of a coherent rhyolite dome and a welded rhyolitic dome breccia. We examine the textural evidence for the response of the wall material, built of ∼75 vol.% glass and ∼25 vol.% crystals (pore-free equivalent), to mass movement in the adjacent conduit. In the rhyolitic wall material, we quantify the orientation and aspect ratio of biotite crystals as strain markers of simple shear deformation, and interpret juxtaposed regions of vesiculation and vesicle collapse as evidence of conduit wall heating. Systematic changes occur close to the margin: (1) porosity is highly variable, with areas locally vesiculated or densified, (2) biotite crystals are oriented with their long axis parallel to the margin, (3) the biotites have greater aspect ratios close to the margin and (4) the biotite crystals are fractured. We interpret the biotite phenocryst deformation to result from crystal fracture, rotation and cleavage-parallel bookcase translation. These textural observations are inferred to indicate mechanical coupling between the hot gas-ash jet and the conduit wall and reheating of wall rock rhyolite. We couple these observations with a simple 1D conductive heating model to show what minimum temperature the conduit wall needs to reach in order to achieve a temperature above the glass transition throughout the texturally-defined deformed zone. We propose that conduit wall heating and resulting deformation influences conduit margin outgassing and may enhance the intensity of such large basaltic eruptions.

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.

Thermal vesiculation during volcanic eruptions.

Terrestrial volcanic eruptions are the consequence of magmas ascending to the surface of the Earth. This ascent is driven by buoyancy forces, which are enhanced by bubble nucleation and growth (vesiculation) that reduce the density of magma. The development of vesicularity also greatly reduces the 'strength' of magma, a material parameter controlling fragmentation and thus the explosive potential of the liquid rock. The development of vesicularity in magmas has until now been viewed (both thermodynamically and kinetically) in terms of the pressure dependence of the solubility of water in the magma, and its role in driving gas saturation, exsolution and expansion during decompression. In contrast, the possible effects of the well documented negative temperature dependence of solubility of water in magma has largely been ignored. Recently, petrological constraints have demonstrated that considerable heating of magma may indeed be a common result of the latent heat of crystallization as well as viscous and frictional heating in areas of strain localization. Here we present field and experimental observations of magma vesiculation and fragmentation resulting from heating (rather than decompression). Textural analysis of volcanic ash from Santiaguito volcano in Guatemala reveals the presence of chemically heterogeneous filaments hosting micrometre-scale vesicles. The textures mirror those developed by disequilibrium melting induced via rapid heating during fault friction experiments, demonstrating that friction can generate sufficient heat to induce melting and vesiculation of hydrated silicic magma. Consideration of the experimentally determined temperature and pressure dependence of water solubility in magma reveals that, for many ascent paths, exsolution may be more efficiently achieved by heating than by decompression. We conclude that the thermal path experienced by magma during ascent strongly controls degassing, vesiculation, magma strength and the effusive-explosive transition in volcanic eruptions.

Heterogeneity: The key to failure forecasting.

Elastic waves are generated when brittle materials are subjected to increasing strain. Their number and energy increase non-linearly, ending in a system-sized catastrophic failure event. Accelerating rates of geophysical signals (e.g., seismicity and deformation) preceding large-scale dynamic failure can serve as proxies for damage accumulation in the Failure Forecast Method (FFM). Here we test the hypothesis that the style and mechanisms of deformation, and the accuracy of the FFM, are both tightly controlled by the degree of microstructural heterogeneity of the material under stress. We generate a suite of synthetic samples with variable heterogeneity, controlled by the gas volume fraction. We experimentally demonstrate that the accuracy of failure prediction increases drastically with the degree of material heterogeneity. These results have significant implications in a broad range of material-based disciplines for which failure forecasting is of central importance. In particular, the FFM has been used with only variable success to forecast failure scenarios both in the field (volcanic eruptions and landslides) and in the laboratory (rock and magma failure). Our results show that this variability may be explained, and the reliability and accuracy of forecast quantified significantly improved, by accounting for material heterogeneity as a first-order control on forecasting power.

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.

Viscosity measurements of crystallizing andesite from Tungurahua volcano (Ecuador).

Viscosity has been determined during isothermal crystallization of an andesite from Tungurahua volcano (Ecuador). Viscosity was continuously recorded using the concentric cylinder method and employing a Pt-sheathed alumina spindle at 1 bar and from 1400°C to subliquidus temperatures to track rheological changes during crystallization. The disposable spindle was not extracted from the sample but rather left in the sample during quenching thus preserving an undisturbed textural configuration of the crystals. The inspection of products quenched during the crystallization process reveals evidence for heterogeneous crystal nucleation at the spindle and near the crucible wall, as well as crystal alignment in the flow field. At the end of the crystallization, defined when viscosity is constant, plagioclase is homogeneously distributed throughout the crucible (with the single exception of experiment performed at the lowest temperature). In this experiments, the crystallization kinetics appear to be strongly affected by the stirring conditions of the viscosity determinations. A TTT (Time-Temperature-Transformation) diagram illustrating the crystallization "nose" for this andesite under stirring conditions and at ambient pressure has been constructed. We further note that at a given crystal content and distribution, the high aspect ratio of the acicular plagioclase yields a shear-thinning rheology at crystal contents as low as 13 vol %, and that the relative viscosity is higher than predicted from existing viscosity models. These viscosity experiments hold the potential for delivering insights into the relative influences of the cooling path, undercooling, and deformation on crystallization kinetics and resultant crystal morphologies, as well as their impact on magmatic viscosity.

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.

Sedimentology and geomorphology of the deposits from the August 2006 pyroclastic density currents at Tungurahua volcano, Ecuador.

The deposits of the pyroclastic density currents from the August 2006 eruption of Tungurahua show three facies associations depending on the topographic setting: the massive, proximal cross-stratified, and distal cross-stratified facies. (1) The massive facies is confined to valleys on the slopes of the volcano. It contains clasts of >1 m diameter to fine ash material, is massive, and interpreted as deposited from dense pyroclastic flows. Its surface can exhibit lobes and levees covered with disk-shaped and vesicular large clasts. These fragile large clasts must have rafted at the surface of the flows all along the path in order to be preserved, and thus imply a sharp density boundary near the surface of these flows. (2) The proximal cross-stratified facies is exposed on valley overbanks on the upper part of the volcano and contains both massive coarse-grained layers and cross-stratified ash and lapilli bedsets. It is interpreted as deposited from (a) dense pyroclastic flows that overflowed the gentle ridges of valleys of the upper part of the volcano and (b) dilute pyroclastic density currents created from the dense flows by the entrainment of air on the steep upper flanks. (3) The distal cross-stratified facies outcrops as spatially limited, isolated, and wedge-shaped bodies of cross-stratified ash deposits located downstream of cliffs on valleys overbanks. It contains numerous aggrading dune bedforms, whose crest orientations reveal parental flow directions. A downstream decrease in the size of the dune bedforms, together with a downstream fining trend in the grain size distribution are observed on a 100-m scale. This facies is interpreted to have been deposited from dilute pyroclastic density currents with basal tractional boundary layers. We suggest that the parental flows were produced from the dense flows by entrainment of air at cliffs, and that these diluted currents might rapidly deposit through "pneumatic jumps". Three modes are present in the grain size distribution of all samples independently of the facies, which further supports the interpretation that all three facies derive from the same initial flows. This study emphasizes the influence of topography on small volume pyroclastic density currents, and the importance of flow transformation and flow-stripping processes.

Kimberlite ascent by assimilation-fuelled buoyancy.

Kimberlite magmas have the deepest origin of all terrestrial magmas and are exclusively associated with cratons. During ascent, they travel through about 150 kilometres of cratonic mantle lithosphere and entrain seemingly prohibitive loads (more than 25 per cent by volume) of mantle-derived xenoliths and xenocrysts (including diamond). Kimberlite magmas also reputedly have higher ascent rates than other xenolith-bearing magmas. Exsolution of dissolved volatiles (carbon dioxide and water) is thought to be essential to provide sufficient buoyancy for the rapid ascent of these dense, crystal-rich magmas. The cause and nature of such exsolution, however, remains elusive and is rarely specified. Here we use a series of high-temperature experiments to demonstrate a mechanism for the spontaneous, efficient and continuous production of this volatile phase. This mechanism requires parental melts of kimberlite to originate as carbonatite-like melts. In transit through the mantle lithosphere, these silica-undersaturated melts assimilate mantle minerals, especially orthopyroxene, driving the melt to more silicic compositions, and causing a marked drop in carbon dioxide solubility. The solubility drop manifests itself immediately in a continuous and vigorous exsolution of a fluid phase, thereby reducing magma density, increasing buoyancy, and driving the rapid and accelerating ascent of the increasingly kimberlitic magma. Our model provides an explanation for continuous ascent of magmas laden with high volumes of dense mantle cargo, an explanation for the chemical diversity of kimberlite, and a connection between kimberlites and cratons.