When Does Magma Break

Geophysical signals arriving at the Earth’s surface originate from a source mechanism at depth but are not necessarily directly observable. Therefore, well-posed experiments can provide insights into source mechanics and, importantly, the parameters required to model aspects of the sources of unrest signals. In this Chapter we detail one such example of how experimental laboratory work has improved our understanding of unrest signals. We focus on the failure of single- and multi-phase magmas, demonstrating that the liquid viscosity, and therefore the temperature and volatile content of a magma of a given composition, is the limiting parameter in determining whether a magma will ascend viscously or whether it can fracture during ascent. This critical threshold is characterized by a Deborah number, the ratio of the timescale of relaxation to the timescale of local flow. We show that for single-phase magmatic liquids and for vigorously vesiculating magmas, a local Deborah number of \( 10^{ - 2} \) is the limit above which mixed viscoelastic behaviour including fracture propagation can be expected, and a Deborah number of \( 1 \) is the limit above which magma is dominantly elastic and responds in a brittle manner to applied stresses. These thresholds can be understood in terms of the onset and peak of the Debye relaxation process for viscoelastic liquids. The apparent validity of a Maxwell model permits us to predict the maximum stress that can be supported by a volcanic liquid deforming in the high Deborah number range. We use these constraints to provide a map of timescales on which we contour dominant system responses from viscous to purely brittle; valid for all magmatic liquids. Finally, we explore the scaling necessary to extend these conceptual insights to crystal- and bubble-bearing magmas valid under specific conditions. The competing timescales of deformation and relaxation in magma are relevant to unrest source mechanisms that originate from magma deformation, such as long-period seismic signals that are used to predict eruption timing.


Introduction
There are a wide variety of observable unrest signals at active volcanoes. The key unrest signals are those that are diagnostic of a new regime of behaviour or of how the system may evolve in the future. One such family of events are the low-frequency earthquakes at volcanoes, which are thought to directly represent magma movement at moderate to shallow depths (Chouet et al. 1994;Chouet 1996;Neuberg et al. 2006;Salvage and Neuberg 2016). Such events are powerful tools particularly because they are closely associated with surface activity (Miller et al. 1998). Accelerating event rates of low frequency earthquakes can be used to retrospectively predict an eruption or dome collapse event with reasonable accuracy (Salvage and Neuberg 2016). The utility of the low-frequency event type has been solidified since physical hypotheses for the source mechanism have been put forward (Neuberg et al. 2006) and tested in scaled laboratory expeirments (Benson et al. 2008;Tuffen et al. 2008). The phenomenological observation of the use of certain signal types over others is useful and has led to accurate retrospective eruption forecasts (Voight 1988;Chouet 1996). However, knowledge of the physics of the source leads to a more diverse range of tools being deployed to understand the evolution of volcanic systems into the future. In the context of low-frequency earthquakes at volcanoes, once a source mechanism is identified (Tuffen et al. 2003;Neuberg et al. 2006), numerical models of magma ascent in the conduit beneath volcanoes can be used to reproduce the depths of the source from first principles (Thomas and Neuberg 2012) and could eventually make forecasts of other observables at the surface that would be consistent with impending eruption.
In this Chapter we explore the source mechanism of low-frequency earthquakes at volcanoes from a physical perspective using a compilation of data from scaled laboratory experiments. We use these datasets to demonstrate that the critical threshold for fracture propagation in a viscoelastic fluid such as magma is universal. We propose that this potentially simplifies the inputs required for effective modelling of source mechanics of low-frequency events. More than this, we hope that this provides a good example of how laboratory work can provide valuable insight into the physical feasibility of models proposed to explain unrest signals at volcanoes.
2 Scaling the Viscous-to-Brittle Transition in Magmas

Single-Phase Magmatic Liquids
Newtonian viscous liquids deform and flow under applied stresses such that the rate of shear strain _ c resulting from an applied shear stress s is proportional to the viscosity l via s ¼ l_ c. However, as large shear stresses are applied to viscous liquids, they can exhibit viscoelastic behaviour, which manifests as an apparent non-Newtonian relationship between s and _ c and can result in fracture propagation locally in the liquid. Of interest here is the transition from simple, apparently Newtonian viscous behaviour to complex viscoelastic behaviour.
Maxwell described the simplest viscoelastic model in which a liquid has a characteristic time required to relax an applied shear stress, termed the Maxwell relaxation timescale k r , which is given by k r ¼ l=G 1 where G 1 is the shear modulus in the purely elastic regime. The Maxwell model of viscoelasticity permits us to predict that if a shear strain is imposed on a viscoelastic liquid, the resultant shear stress will rise rapidly to a peak value s i and will relax over time t according to s t ð Þ ¼ s i expðÀt=k r Þ. This simple concept has proved powerful in describing the transition between Newtonian and viscoelastic behaviour in volcanic liquids (Dingwell 1995(Dingwell , 1996. Indeed, this model is invoked to explain the transition between volcanic liquids and volcanic glasses on cooling (e.g. Stevenson et al. 1995), the fragmentation of fluid and bubbly magma undergoing decompression (Alidibirov and Dingwell 1996;Kameda and Kuribara 2008) and the conditions under which extensive shear fracture networks can form at or near volcanic conduit margins (Gonnermann and Manga 2003;Tuffen et al. 2003;Kendrick et al. 2014;Hornby et al. 2015). In all cases, it is useful to define a Deborah number De, which is the dimensionless ratio between the Maxwell relaxation time and the timescale of deformation k. For simple (viscometric) shearing flow, the latter timescale is 1=_ c and thus Our hypothesis is that a Deborah number of unity separates the regime between viscous ðDe ( 1) and elastic (De ) 1Þ responses to shear stresses, or equivalently, between a system that can relax shear stresses efficiently from one in which shear stresses can accumulate elastically. In two different experimental types, both Webb and Dingwell (1990a) and Cordonnier et al. (2012a, b) showed that in fact the first evidence of behaviour that is not purely viscous occurs at De ! 10 À2 and not exactly at unity. This observation will be discussed later.
The power of Eq. 1 is that simple parameters that can be measured in the laboratory can be used to predict where this transitional point of De ¼ 1 will be. This critical Deborah number can be termed De 0 : In what follows, we give examples of how l and G 1 have been found for a range of volcanic liquid compositions at magmatic temperatures and volatile contents before exploring the validity and applicability of the scaling in Eq. 1. The dimensionless nature of Eq. 1 means that it can be assessed for any system so long as k can be defined. At the core of the experimental toolkit is the determination of the fundamental quantities that are required to model volcanic processes involved in unrest. The most variable physical quantity in volcanic systems is the viscosity of the liquid phase. Using a range of techniques, this viscosity can be determined with prodigious accuracy and significant effort has been expended in mapping the full range of composition, temperature (see Giordano et al. 2008) and volatile content (Hess and Dingwell 1996) relevant to shallow magmatic settings. Multicomponent models have been proposed such that the viscosity of any composition of silicate volcanic liquid on Earth can now be predicted as a function of temperature including the range of conditions between shallow magma storage and the surface (Hess and Dingwell 1996;Giordano et al. 2008). This provides parameterization of lðTÞ for use in Eq. 1. The most commonly used form for lðTÞ is the non-Arrhenian Vogel-Fulcher-Tammann equation of the form where A; B; and C are coefficients that are experimentally determined and then parameterized as a function of families of oxides in the liquid structure (Giordano et al. 2008) or as a function of the dissolved water content (for calc-alkaline rhyolites: Hess and Dingwell 1996), to give two examples.
To give examples of the range of viscosities of interest in volcanic scenarios, we give end members in Fig. 1 for a calc-alkaline rhyolite using the model from Hess and Dingwell (1996), and for the basaltic liquid composition provided in Zhang et al. (1991) using the model from Giordano et al. (2008), both contoured for a range of water contents.
Unlike viscosity, the shear modulus in the elastic regime does not vary significantly in siliciate liquids and is not strongly dependent on composition or temperature. Indeed, a short survey of the values for G 1 in silicate glasses and liquids at a range of temperature and a huge range of composition, provided by Dingwell and Webb (1989) shows that G 1 ¼ 10 10AE0:5 Pa. Here we use this range in order to fully parameterize De (by predicting k r Þ as a function of temperature (via Eq. 2) for any silicate liquid in the shallow crust.

Extensions to Multiphase Magmas
Except in rare circumstances such as obsidian-forming eruptions, volcanic liquids are not often erupted without some proportion of suspended pore space (either as isolated bubbles or connected networks) and rigid crystals. In either case, the utility of Eq. 1 requires additional attention. We posit that for the liquid phase between the pores or the crystals, De given by Eq. 1 holds. However, we acknowledge that estimation or measurement of the rate of shear strain locally between pores or crystals would be difficult (Deubelbeiss et al. 2011). Therefore, a more robust criterion for the viscous-to-brittle transition in crystal-or bubble-bearing magma would require that we scale the bulk strain rate _ c b on the system for the effect of the suspended load.
For crystals, we show some first-order scaling attempts for their effect on the critical threshold for the onset of brittle behaviour in multiphase magma. The simplest view of the local flow of liquid in crystal-bearing magma under constant shear stress is that the rate of shear strain between the crystals should scale approximately with _ c ¼ where / x is the suspended crystal volume fraction and / m is a jamming fraction above which no more crystals can be added to the flowing system. This scaling would imply that the Deborah number for a crystal-bearing magma De x would be where / m is a function of crystal shape and roughness (Mueller et al. 2010;Mader et al. 2013). A critical value of De x , termed De 0 x ; is 10 À2 ; consistent with the limiting De 0 at / ¼ 0.
A new definition of the bulk failure criterion De 0 can now be made, which decreases as / ! / m so that De 0 ¼ De Fig. 2 we demonstrate how this concept predicts a linear relationship between De x and / x =/ m , which in turn shows that lower bulk strain rates are required to induce brittle behaviour in crystal-bearing magmas.
This hypothesis was found to hold for simple two-phase systems by Cordonnier et al. (2012a), but where those authors required that l be replaced by the suspension viscosity of the whole system, which is at odds with the scaling of strain rate for local liquid effects only. Nonetheless, this approach described their data satisfactorily and the conceptual insight that crystals locally increase the liquid rate of shear strain relative to the bulk value is robust, with the implication is that the whole suspension will begin fracturing at lower bulk rates of strain relative to a single-phase liquid of the same composition as crystals are added. It may be that additional second-order effects are important at high / x , which are not accounted for in this simple analysis.
For the case of porous magmas, there are two considerations: (1) The growth of bubbles exerts a rate of shear strain in the liquid concentrated at the bubble walls, which is broadly independent of any bulk shearing deformation, and (2) like for the crystal case, the presence of bubbles changes the partitioning of the bulk rate of shear strain in the liquid between the bubbles. In the case of scenario (1) (Hess and Dingwell 1996) for 0.01-1 wt% water, and a typical basaltic composition (composition from Zhang et al. 1991) calculated using the multicomponent viscosity model (Giordano et al. 2008) for 0-1 wt% water. The relaxation timescale is calculated assuming a composition independent value for G 1 of 10 10 Pa (Dingwell and Webb 1989) pressure in the bubble, relative to the hydrostatic pressure (Sparks 1978). This can be augmented for the diffusion-controlled gradient of viscosity in the immediate liquid shell around a growing bubble (Prousevitch et al. 1993;Lensky et al. 2001). The component of the rate of shear strain tangential to the bubble wall _ c h can then be computed throughout bubble growth (Ichihara et al. 2002). If we were then to input _ c h into our computation of De, we would predict under which conditions the system would meet the criterion of De 0 ¼ 10 À2 locally at the bubble rim. These conditions would be best cast in terms of the critical rate of pressure or temperature change, or bulk initial volatile content that would allow the bubble to grow sufficiently fast that fractures could be propagated at the bubble wall.
In case (2), the effect of bubbles on the bulk Deborah number is less clear. We speculate that at De ) 1, the stress required to fracture a bubbly system will be analogous to that required to fracture a vesicular glass. In this case, models for the effect of spherical cavities on the stress required for fracturing are valid and they predict that the stress is reduced significantly as the bulk gas volume fraction increases (Sammis and Ashby 1986). This has been confirmed in the high Deborah number regime for porous liquids analogous to volcanic systems (Vasseur et al. 2013) but remains untested in the low Deborah number regime where bubble deformation may be important.

Apparent Non-newtonian Effects
In single-phase liquids at Deborah numbers below the critical threshold at which fracturing is observed, there is evidence that during steady state shearing flow, there is a non-Newtonian relationship between applied shear stress and resultant shear rate of strain (Simmons et al. 1982;Dingwell and Webb 1989;Webb and Dingwell 1990a;Cordonnier et al. 2012b). In the onset of this non-Newtonian, the onset of this non-Newtonian behaviour appears to be well described by the value above which viscous dissipation of heat is active on the system length scale of interest (Costa and Macedonio 2005;Hess et al. 2008). This can be scaled by the Brinkman number Br where U g and U l are the gain and loss power densities, respectively, q is the areal heat flux out of the system and k is the thermal conductivity of  Fig. 2 The critical Deborah number required for fracturing De 0 is reduced when crystals are suspended in magma. To a first-order, the De 0 value is reduced proportional to De 0 ¼ De 0 x ð1 À / x =/ m Þ the material. U g ¼ l_ c 2 represents the amount of energy produced in a system of a given volume by viscous dissipation of heat, while U l ¼ kq represents the energy lost due to diffusive thermal equilibration. When Br ) 1, heat is efficiently produced and inefficiently lost from the system, resulting in a bulk temperature increase in the liquid. This would manifest itself as an apparent shear thinning rheology if the temperature increase where not locally accounted for, and would be most likely to be operative at high viscosities and high shear strain rates. We note that unlike the Deborah number, the Brinkman number is scale dependent (as q depends on the area available for heat transfer out of the system) and so should be assessed for each system scale separately.

The Universal Breaking Timescales of Volcanic Liquids
Laboratory data related to the viscous-to-brittle transition in magmas has been collected in a variety of geometries. Using single-phase liquids, there are two dominant geometries: (1) thin fibers of silicate liquid of basaltic, andesitic, phonolitic and rhyolitic composition were stretched under constant load in tension until the fibers snapped in a singular fracture event (Webb and Dingwell 1990b), and (2) cylinders of synthetic borosilicate liquid were compressed under constant load and the bulk temporal evolution of the rate of axial shortening was determined as viscous (relaxed) if it continuously increased to a steady or near-steady value, or brittle (unrelaxed) if the axial shortening rate jumped due to fracturing events (Cordonnier et al. 2012b). In another type of experiment, analogue vesicular liquids were decompressed at different rates from pressure (Kameda and Kuribara 2008;Kameda et al. 2013). In this type of decompression experiment, if the dominant response was viscous (relaxed), then the sample was seen to grow due to decompression-driven bubble growth, and if the dominant response was brittle (unrelaxed), then smooth sample inflation was punctuated by visible fractures opening. In these decompression experiments there additionally was a violent rupture mode in which the fracturing was pervasive and shattered the sample in a vigorous fragmentation event. These data are selected here because the liquids used have a known relaxation time under the conditions used in the experiments and care was taken by the authors who originated the work to investigate the shear rates of strain local to the bubbles (discussed above). In Fig. 3 we map these experimental results as a ratio of the deformation timescale to the relaxation timescale, which permits us to contour the plot for critical Deborah numbers. We colour-code the data according to the bulk mode of response of the sample to the deformation using green to represent purely viscous relaxed behaviour, orange to represent brittle unrelaxed behaviour and red to represent the complete violent rupture of the sample. These data suggest that the viscous relaxed behaviour transitions to unrelaxed brittle behaviour at a Deborah number of 10 À2 , and are consistent across a huge range of experimental conditions and across both the single-phase compression and vesicular decompression experiment types. Furthermore, we see that a Deborah number of unity consistently separates the experiments for which the bulk response was unrelaxed and brittle from those for which the response was violent rupture, fragmentation or complete failure. It appears that these two thresholds, De ¼ 10 À2 and De ¼ 1 are universal even when comparing analogue room temperature liquids with high temperature silicate liquids deformed in a variety of ways. This lends power to the scaling provided by the dimensionless Deborah number and implies that we need only to define the liquid viscosity and a characteristic rate of deformation in order to predict whether a system will flow viscously or rupture violently. For example, the working viscosity of the analogue fluid used in Kameda and Kuribara (2008) is 10 0 l 10 10 Pa.s, which extends to much lower used in Cordonnier et al. (2012b), and yet the scaling with the Deborah number holds across a broad range simply because the deformation timescale was also smaller in the former example.
In the viscous field, Cordonnier et al. (2012b) additionally recorded whether samples hosted a measurable temperature increase due to viscous dissipation of heat during the experiment, or not. In Fig. 3, we plot those with a measureable temperature increase as filled symbols and those without this feature as unfilled symbols. To explain this, we plot the threshold dimensionless Brinkman number of unity Br ¼ 1, using an estimation of the loss power density for the furnace and sample size used U l ¼ 10 4:5 W.m À2 (Cordonnier et al. 2012b). We note that the Brinkman number curve consistently divides the regimes of purely isothermal experiments and those with measurable heat gain. On this map, the position of Br ¼ 1 is non-unique and depends on sample size, such that on the scale of a volcanic conduit, for example, Br ¼ 1 would occur at much higher deformation timescales for the same relaxation timescale (the curve would shift up in Fig. 3). This implies that in the natural case, viscous dissipation of heat may be far more important than shown here for the sample length scale (Costa and Macedonio 2005;Mastin 2005;Costa et al. 2007). Nevertheless, the Deborah number limits discussed appear to be universal and, importantly, are scale-independent.
To gain helpful physical insight into why 10 À2 De 10 0 is the transitional window between a purely viscous and a purely elastic response of a liquid to a deformation, we provide low-strain, high frequency oscillatory rheological measurement data for similar liquids (Fig. 4).
Here, rods of single-phase liquid are subjected to a low-amplitude oscillatory strain with a forcing frequency x at a range of temperatures similar to those used in experiments presented in Fig. 2. Here, x is normalized with k r , yielding a dimensionless frequency or, equivalently, a Webb & Dingwell, 1990Cordonnier et al., 2012Kameda et al., 2008Kameda et al., 2013 Fig. 3 A map of behaviour arising from deformation experiments on natural (Webb and Dingwell 1990b) and synthetic (Cordonnier et al. 2012b) silicate liquids and on analogues of magmatic liquids undergoing vesiculation (Kameda and Kuribara 2008;Kameda et al. 2013). The experimental results were obtained by applying a characteristic timescale of deformation k of a liquid with a characteristic timescale of stress relaxation k r . Marked are ratios k r =k; equivalent to Deborah numbers De, of 10 À2 and 1, which separate experiments with a purely viscous response from those with an unrelaxed brittle response or a violent rupture response, respectively. Additionally marked is the Brinkman number Br of 1, which separates the boundary between isothermal viscous behaviour (above the line) and viscous behaviour in which the material heats up due to viscous dissipation as heat (below the line) frequency-domain version of the Deborah number. The metric that we choose to track is the imaginary component of the complex elastic modulus G 00 normalized by G 1 . When this value G 00 =G 1 is close to zero in the low Deborah number limit, the system is dominated by liquid behaviour and when it rises from zero at increasing Deborah number, there is an increasing component of the response to the forced oscillation that is elastic. We present collated data for a sodium disilicate synthetic composition (Webb 1991), rhyolitic compositions (Webb 1992;Bagdassarov and Dingwell 1993b) and a synthetic composition used as an analogue for type calc-alkaline rhyolite systems (Bagdassarov and Dingwell 1993a). This imaginary component can be described by the generalized Debye model for viscoelasticity for systems with a single characteristic relaxation time as G 00 ¼ xk r = xk r ð Þ 2 þ 1 h i . The data reproduce the broad shape of the Debye model, albeit with an overprediction of G 00 around xk r $ 1 for the rhyolitic liquids, which is poorly understood. Nevertheless, we point out that the first onset of a measurable elastic component to the response of the liquid to deformation occurs at De ¼ xk r % 10 À2 . Similarly, the point above which the majority of the response to deformation is elastic (the peak of G 00 ) occurs at De ¼ xk r % 1 (Fig. 3). We use this observation to validate the two thresholds found in Fig. 3. This implies that fractures can propagate in silicate liquids when there is even a small component of elastic behaviour (De ! 10 À2 ) and that those fractures can propagate vigorously when the elastic behaviour dominates over the viscous behaviour (De ! 1). An interpretation might also be that 10 À2 De 1 is the range in which fracture propagation is competing with viscous relaxation of stress and therefore the fractures are unlikely to be long, sharp-tipped or pervasive. And that at De ! 1, stress dissipation by fracture Normalized imaginary component, Na 2 Si 2 O 5 ; Webb, 1991 Rhyolite;Webb, 1992 Rhyolite;Bagdassarov & Dingwell, 1993 HPG8; 10 − 3 10 − 1 10 1 10 3 10 5 10 7 Deborah number (frequency domain), De Fig. 4 The response of synthetic (Webb 1991;Bagdassarov and Dingwell 1993a) and natural (Webb 1992;Bagdassarov and Dingwell 1993b)  propagation can be localized onto longer, sharp and pervasive fracture networks. The fact that a viscous limit to the expansion of vesicular magma also scales with our Deborah number criterion is tantalizing. Microphysically, in this case it is not fracture of a deforming homogeneous liquid, but fractures propagating at bubble walls as the local rate of strain in the liquid induced by expansion of the bubble meets the Deborah number criterion for fracture propagation (c.f. Kameda and Kuribara 2008). For complete rupture of the vesicular materialwhich is a fragmentation event or violent rupture in Fig. 3-the fractures propagating from the bubbles must interact. Presumably in the range 10 À2 De\1, fractures propagate but the fracture tips are blunted during competing viscous relaxation such that brittle behaviour can be observed but is not catastrophic to the system. Then, as the local strain rate at the bubble wall increases further and De approaches 1, the fractures propagate in a dominantly elastic medium and can span the inter-bubble distances, interacting to produce violent rupture. Therefore, it is clear that a scaling of the Deborah number concept to bubble wall dynamics would provide a fragmentation criterion for viscoelastic vesiculating magma (c.f. Namiki and Manga 2005;Koyaguchi et al. 2008;Namiki and Manga 2008).

Laboratory-Scale Unrest Signals
When magma breaks in a laboratory experiment, acoustic emissions-packets of acoustic energy -are released and can be recorded (Benson et al. 2008;Lavallee et al. 2008). These signals appear to represent large total released amounts of energy when k is short (high _ cÞ, compared to when k is long (low _ cÞ ). This supports our posit that fracture networks are more likely to be pervasive and large when De is large than when De is small. For single phase silicate liquids deformed at high De, Tuffen et al. (2008) showed that the experimental acoustic event frequency range and sample fracture length scale scaled with natural constraints of in-conduit volcanic fracture systems and natural frequencies of low-frequency volcano-seismicity, indicating that these viscoelastic fracturing events are indeed the likely source mechanism for low-frequency earthquakes at volcanoes. This confirmed the conclusion of Neuberg et al. (2006), who showed that these low-frequency events were most likely to be associated with repetitive fracturing of magma at a given depth in the conduit during ascent. Neuberg et al. (2006) further modelled magma ascent in confined geometry and showed that a De $ 1 is met at a depth of 830 m below the surface using parameters typical for recent eruptions at Soufriere Hills volcano and the magma thereof. Thomas and Neuberg (2012) predicted a deeper source of 1500 m using the same model approach but by invoking a conduit restriction (see Chap. 9), consistent with inversions for low-frequency sources using seismic data. Therefore, viscoelastic magma fracturing in the high Deborah number regime appears to be a consistent model for the source mechanics of low-frequency volcano seismicity, often used for eruption forecasting. Vasseur et al. (2015) showed that the forecastability of full sample rupture scales with the heterogeneity of the system-cast most simply as a porosity (Fig. 5). The implication is that the less vesicular the magma undergoing deformation, the less likely that accurate forecasts based on accelerated precursory signals can be made, with up to *120% error on the timing of the rupture event for single-phase homogeneous liquids in the high De limit. It may be that the magma vesicularity, pore-network structure (Vasseur et al. 2017), crystallinity and textural anisotropy, play key roles in determining forecasting success based on low-frequency earthquakes. More experimental work is clearly required in this area, along with more rigorous scaling between acoustic and seismic events that originate from magma failure.
An implication of the models explored here, encapsulated by Eq. 1, is that the peak stress supported by a liquid r m can be predicted as a function of De by r m ¼ G 1 De. Once De ! 10 À2 , however, this linear relationship appears to be invalid as fractures can form in the liquid. It is perhaps significant that the average strength of glass in the high De regime is $ 10 8 Pa (Simmons et al. 1982;Vasseur et al. 2013 Error on magma failure forecast Sintered glass beads; Vasseur et al., 2015;2017 Welded volcanic debris; Vasseur et al., 2017 Fig. 5 In the high-De regime, the error on a prediction of failure times scales with the porosity of the material such that low porosity magmas are unpredictable and high porosity magmas are predictable (Vasseur et al. 2015).
The failure times are recorded as a large stress drop during uniaxial loading, and the approach to failure is monitored using acoustic emissions generated by pre-failure micro-fracturing events 10 − 6 10 − 4 10 − 2 10 0 10 2 10 4 10 6 Deborah number, De σ g G ∞ Fig. 6 The peak stress that can be supported by a shearing liquid in the low De limit is given by r m ¼ G 1 De: And in the high De limit, this value appears to asymptotically approach r m ¼ 10 À2 G 1 which is equivalent to the strength of glassy materials in that same regime (Simmons et al. 1982) with a rupture threshold of 10 À2 G 1 . These relationships are explored in Fig. 6.
Here we summarize just one example of how targeted laboratory experiments and dimensional analysis have demonstrated the most likely source mechanism of unrest signals commonly used to monitor magma movement and predict impending eruptions. Clearly, with scaling arguments for the applicability of an experimental set up-such as we show here with the Deborah number analysis-laboratory-based work can provide new insights into unrest mechanisms that facilitate more accurate forward-modelling of geophysical signals.