Crystal plasticity as an indicator of the viscous-brittle transition in magmas

Weighted Burgers Vector analysis of orientation fields from high-angular resolution electron backscatter diffraction

The Weighted Burgers Vector (WBV) method can extract information about dislocation types and densities present in distorted crystalline materials from electron backscatter diffraction (EBSD) maps, using no assumptions about which slip systems might be present. Furthermore, high-angular resolution EBSD (HR-EBSD) uses a cross-correlation procedure to increase the angular precision of EBSD measurements by an order of magnitude compared to conventional EBSD. However, the WBV technique has not previously been applied to HR-EBSD data and therefore it remains unclear as to which low-angle substructures can be reliably characterised by WBV analysis of conventional EBSD data and which require additional HR-EBSD processing. To establish some practical examples that can be used to guide future data-acquisition strategies, we compare the output of the WBV method when applied to conventional EBSD data and HR-EBSD data collected from the most common minerals in Earth's lower crust (plagioclase feldspar) and upper mantle (olivine). The results demonstrate that HR-EBSD and WBV processing are complementary techniques. The increase in angular precision achieved with HR-EBSD processing allows low-angle (on the order of 0.1°) structures, which are obscured by noise in conventional EBSD data, to be analyzed quantitatively using the WBV method. Combining the WBV and HR-EBSD methods increases the precision of calculated WBV directions, which is essential when using information about active slip systems to infer likely deformation mechanisms from naturally deformed microstructures. This increase in precision is particularly important for low-symmetry crystals, such as plagioclase, that have a wide range of available slip systems that vary in relative activity with changing pressure, temperature and differential stress. Because WBV directions are calculated using no assumptions about which slip systems may be present, combining this technique with HR-EBSD to refine the precision of lattice orientation gradients is ideal for investigating complex natural materials with unknown deformation histories.

Caldera-forming eruptions of mushy magma modulated by feedbacks between ascent rate, gas retention/loss and bubble/crystal framework interaction

Caldera-forming eruptions of mushy silicic magma are among the most catastrophic natural events on Earth. In such magmas, crystals form an interlocking framework when their content reaches critical thresholds, resulting in the dramatic increase in viscous resistance to flow. Here, we propose a new mechanism for the ascent of mushy magma based on microstructural observations of crystal-rich silicic pumices and lavas from the Central Andes and decompression experiments. Microstructural data include spherical vesicles and jigsaw-puzzle association of broken crystals in pumices, whereas there is limited breakage of crystals in lavas. These observations insinuate that shearing of magma during ascent was limited. Decompression experiments reveal contrasting interaction between growing gas bubbles and the crystal framework in crystal-rich magma. Under slow decompression typical of effusive eruptions, gas extraction is promoted, whereas under rapid decompression, bubbles are retained and the crystal framework collapses. This feedback between decompression rate, retention of gas bubbles, and integrity of the crystal framework leads to strong non-linearity between magma decompression rate and eruption explosivity. We extend these findings to caldera-forming eruptions of crystal-rich magma where large overpressures are induced by caldera-collapse, resulting in magma plug-flow, rapid decompression facilitated by shear-localization at conduit margins, and explosive eruption.