Dual chamber-conduit models of non-linear dynamics behaviour at Soufrière Hills Volcano, Montserrat
|
|
- Jemimah Copeland
- 5 years ago
- Views:
Transcription
1 Dual chamber-conduit models of non-linear dynamics behaviour at Soufrière Hills Volcano, Montserrat Oleg Melnik 1, Antonio Costa 2 1 Institute of Mechanics, Moscow State University, Moscow, Russia, and Earth Science Department, University of Bristol, Bristol, UK. 2 Environmental Systems Science Centre, University of Reading, Reading, UK, and Istituto Nazionale di Geofisica e Vulcanologia, Naples, Italy. To be submitted to the special volume The Eruption of Soufriere Hills Volcano, Montserrat from 2000 to 2010, Editors: G. Wadge, R. Robertson, B. Voight, Memoir of the Geological Society of London. 1. Introduction Lava dome eruptions are typically characterized by large fluctuations in discharge rate showing cyclic patterns of seismicity, ground deformation and volcanic activity. Time-scales associated with these fluctuations range from hours-days (short-term cycles) to years-decades (long-term cycles) with intermediate of 5 to 7 weeks as, for example, observed on the Soufrière Hills Volcano (SHV), Montserrat. Concerning short-term cycles, cyclic patterns of seismicity, ground deformation and volcanic activity have been documented at Mount Pinatubo, Philippines, in 1991 and Soufrière Hills Volcano, Montserrat, British West Indies, in At SHV, periodicity in seismicity and tilt ranged from ~3 to 30 hours (Voight et al., 1999; Odbert and Wadge, 2010). Several models have been proposed to explain the observed cyclicity in the last decade. In both Denlinger and Hoblitt (1999) and Wylie et al. (1999) models the lower part of the conduit acts like a capacitor that allows magma to be stored temporally in order to release it during the intense phase of the eruption. Costa et al.
2 (2012) generalized Denlinger and Hoblitt (1999) model to a more realistic system made of a compressible magma flowing through an elastic dyke evolving to a shallower cylindrical conduit. Lensky et al. (2008) explained short-term cyclic behaviour at the SHV as a result of gas diffusion into growing bubbles and filtration through the bubble network in stagnated magma column before the critical overpressure is reached followed by magma motion, depressurisation and stagnation of the plug at the top of the conduit. Regarding long-term cycles Costa et al. (2007a;b) considered an elastic dyke evolving to a cylinder toward the surface coupled with a shallow (~5 km) magma chamber. They showed that degassing-induced crystallization within the conduit coupled with the wallrock elasticity play a major control on process for the long-term cyclicity during lava dome building eruptions. In particular, they showed that there is a regime where the period of pulsations is controlled by the elasticity of the dyke (~weeks to months) and a regime where the period is controlled by the volume of the magma chamber (~years). Intermediate regimes are possible. The geometry proposed by Costa et al. (2007a;b) is also consistent with inversion of strain data for the March 2004 SHV explosion that highlight the presence of a spherical pressure source and a shallow dyke (Linde et al., 2010, Hautmann et al., 2009). Most of the existing conduit flow models for lava dome building eruptions assume a system consisting of a magma chamber (a large volume of magma located in elastic wallrocks) connected to the surface through a conduit. Magma chamber is assumed to be fed from below with a constant rate with a fresh magma. In few models (Dirksen et al, 2006; Maeda, 2000) the influx rate into the magma chamber is specified as a function of time. Modern geophysical data recorded during lava dome building eruptions indicate the presence of interconnected multiple magma storage regions. Concerning the SHV system different geometries have been proposed on the basis of both mineralogical and geodetic observations. Elsworth et al. (2008) proposed a model with two stacked magma reservoirs connected by vertical conduits (the shallow chamber was centered at 6 km and the deep at 12 km below sea level). Surface efflux and GPS station velocities were inverted to calculate crustal magma transfer. Results from the inversion also indicated that the eruptive episodes deplete the lower reservoir only, and not the upper reservoir.
3 A vertically elongated chamber idealized as a prolate ellipsoid was proposed by Voight et al. (2010). This chamber was assumed stratified with upper parts compressible due to exsolved gas phases (magma sponge), and fed at the base by deep influx estimated of the order of 2 m 3 /s. The top of the ellipsoid was assumed at 5 km below sea level and its centroid near 10 km. Whether SHV system is better represented by a single prolate chamber or a two chamber model is still not resolved. However Foroozan et al. (2010) showed that the deformations of GPS stations located near or far from the volcano are better explained by a two chamber model. In a more recent paper Foroozan et al. (2011) simultaneously inverted geodetic data and measured magma efflux rates with the additional constraint of constant basal influx. In contrast to Elsworth et al. (2008) they used data from 6 to 13 stations (versus 4 previously) and a compressible magma phase (versus incompressible previously). Their results suggest a basal supply of 1.2 m 3 /s for the activity from Centroids of the chambers were estimated to depths of about 5 and 19 km. In contrast to Elsworth et al. (2008), Foroozan et al. (2011) found that the data are characterized by synchronous deflation then inflation of both shallow and deep chambers. In order to account for the large-strain effects and modulus softening that are unrepresented in the extrapolated 1-D seismic velocity data, Foroozan et al. (2011) decreased the shear modulus of the wallrock G by a factor of 10, estimated on the basis of 1-D seismic velocity profile (Shalev et al., 2010). Rather than this arbitrary approach, we will use, following Costa et al. (2009; 2011), an empirical relationship between static and dynamic values of the Young modulus E of rocks (Wang, 2000): E static (GPa) = E dynamic (GPa) On the other hand the analysis of data for the period from the borehole strainmeter network at SHV gives good constraints on where the pressure changes occur. These observations are also consistent with a system made of a deep reservoir coupled with a shallow chamber connected to the surface through a shallow dyke (Sacks et al., 2011). Most numerical models of these eruptions assume a single magma chamber fed from below with a constant or prescribed time dependent influx rate and connected to the Earth surface through a conduit. Here we present a development of the model of
4 extrusive eruption by Costa et al. (2007a,b) considering a system made of two magma chambers located in elastic rocks and connected by a dyke between each other. The deep chamber is fed from below with a constant influx rate. For sake of comparison with previous studies (e.g., Hautmann et al., 2009; 2010) we use locations and volumes of the magma chambers inferred for SHV, Montserrat, from ground deformation studies by Hautmann et al. (2010) and also the more recent estimations obtained integrating seismic tomography with numerical models of the shallow magma chamber by Paulatto et al. (2012). The estimations of both chamber locations and volumes contain large uncertainties that need to be systematically investigated but here we will focus on exploring the control of two magma chambers on the dynamics of lava dome extrusion only. In order to study the role of a dual magma chamber system, however, we explore a large range of chamber volumes and two different chamber shapes (spheroidal and prolate) for both deep and shallow magma chamber. The model shows cyclic behaviour with a period that depends on the intensity of the influx rate and the chamber connectivity (described as the horizontal extent of the dyke connecting the two chambers). For a weak connectivity the overpressure in the lower chamber stays nearly constant during the cycle and the influx of fresh magma into the shallow chamber is also nearly constant. For a strong connectivity between the chambers their overpressures increases or decreases during the cycle in a synchronous way. Influx into the shallow chamber stays close to the extrusion rate at the surface. 2. A model for magma transfer between multiple reservoirs and dykes We consider a system like that proposed by Hautmann et al. (2010) represented in Fig. 1. The system is composed by a deep magma chamber fed from below and connected to a shallow magma chamber through a dyke. The shallow magma chamber is connected to the surface via a volcanic conduit made of dyke that evolves to a cylinder approaching to the surface as described Costa et al. (2007a,b). The system of equations for the description of magma transport from the deep to the shallow chamber is:
5 p z = ρg 4μ a2 + b 2 V a 2 b 2 S t + SV z = 0; S = πab a = a 0 + ΔP [ 2E (1 2ν)a 0 + 2(1 ν)b 0 ] ΔP b = b0 + [ 2(1 ν ) a0 (1 2ν ) b0 ]. 2E Here μ denotes magma viscosity, ΔP = p p e (z) is a magmatic overpressure, p is (1) 141 the pressure inside the dyke, p e (z) = ρ r gz σ t denotes the lithostatic pressure (g is acceleration due to gravity) minus σ t, i.e. tensile stress along the dyke due to the presence of magma chamber and an extensional/compressional far-field stress. S is the cross-section area of the dyke which is assumed to have an elliptical shape, a and b are the semi-axis of the ellipse. Wallrocks are assumed to be purely elastic with Young modulus E and Poisson ratio ν. The Muskhelishvili solution (Muskhelishvili, 1963) is used for calculation of the relation between magmatic overpressure and the deformation of the dyke. Densities of magma ρ and wallrocks ρ r, in the deeper part of the magmatic system, are assumed to be constant, z is a vertical coordinate increasing upwards, V is the ascent velocity. In order to simplify the equation system, since typically a 0 >> b 0 relative changes of a are negligible, thus, we can assume a constant and the initial value of b equal to zero, i.e. b 0 = 0. Moreover, at typical depths of several kilometers, we can reasonably neglect the presence of the free gas phase and assume incompressible magma of constant viscosity. For CO 2 rich magma this assumption will be violated due to deep exsolution of dissolved volatiles, but for the case of SHV, even assuming the upper bound estimation of CO 2 content in the magma (600ppm wt%, Edmonds et al., 2012) is exsolved, magma density variation is negligible, i.e. from 2520 to 2600 kg/m 3 when pressure changes from 240 to 300 MPa (moreover for such CO 2 content at 300 MPa magma bulk
6 modulus is ~3GPa, only few times smaller than rock rigidity at that depth). Hence we have: z (p + ρgz) = 4μ V b 2 b t + bv z = 0 (2) b = ΔP E (1 ν)a 0 = αδp where α = (1 ν)a 0 parabolic equation: ΔP t E. The above system can be rewritten as a non-linear α 2 ΔP 4 α 2 (ΩΔP 3 ) = 0 (3) 16μ z 2 4μ z where Ω= [ z p e (z) + ρgz]. In order to solve eq. (3) we impose the following boundary conditions: d dt ΔP = Q Q ch, d ( in, d ) P V ch, d ch, d 0 with Q = π a bv = [ ΔP + pe( z) + ρgz] 1 πa b 0 being the volumetric outflux from the 4μ z deep magma chamber and Q in,d being the influx into deep magma chamber from a deeper magmatic source. Q in,d will be considered as a given constant parameter in this study, although it can be influenced by pressure variations in the magmatic source, local tectonic stress and thermal development of the feeding system.
7 Eq. (3) is solved together with the system of equations that describes the flow from the shallow magma chamber through a dyke shaped conduit that evolves to a cylinder in the vicinity of the surface. These equations are described in detail in Costa et al. (2007a,b) and are summarized in the Appendix. FIG. 1 HERE Fig. 1 Sketch of the investigated system. Modified after Hautmann et al. (2010). 3. Parametric study Here we present results of two separate studies. First, we analyze the effects of a farfield stress and a pressurized deep magma chamber on the steady state solution of the flow equation in the upper conduit. Second, assuming a neutral far-field stress we present results on the role of the connectivity between deep and shallow magma chambers during an extrusive eruption on the transient behaviour of the system. 3.1 Effect of far-field stress The presence of a deep reservoir affects the dynamics and the cyclicity of the shallow part of the system, but it also has an important effect in perturbing the stress field around it. A pressurized magma chamber under the effect of a far-field stress can modify the tensile stress along the axis of the conduit (see Costa et al., 2011). This stress is calculated using the general analytical solution by Gao (1996) obtained in the limit of a plane 2D geometry (for the limitations of this assumption see Costa et al., 2011). We explore the effect of far field stresses, typical of arc volcanism environment, from compressional (-10 MPa) to extensional environment (10 MPa) (Hautmann et al., 2010; Roman et al., 2006) on the steady state solution of the flow equation in the upper conduit. Fig. 2 reports a set of steady-state solutions for fixed deep reservoir overpressure of 10 MPa and shows the control of the far-field stress field on the system sketched in Fig. 1. A compressional far-field stress tends to move the steady state solution toward larger pressures, whereas an extensional far-field stress moves the solution to lower
8 pressures allowing larger flow rates at low pressures. The conduit cross-section area is larger in the case of extensional environment, thus the transition to the upper branch moves to higher values of discharge rate because ascent velocity for large is smaller for the same discharge rate and crystals grow more efficiently (see Fig. 12b in Melnik and Sparks, 2005 for more detailed explanations). FIG. 2 HERE Fig. 2. A set of steady state solutions: dependence of discharge rate on chamber pressure for three different far-field stresses: compressional (-10 MPa; dashed-dotted line), neutral (full line), and extensional (10 MPa; dotted line). The remaining parameters are set as in Costa et al. (2007b). 3.2 Effect of coupling the two reservoirs In this section we explore the role of the connectivity between deep and shallow magma chamber during an extrusive eruption. We consider only the effect due to the two pressurized magma chambers and will set 0. As initial condition we assume that both magma chambers have zero overpressure with respect to the lithostatic pressure, the flow rate is zero everywhere and the dyke between the two chambers has a negligible thickness. At time t = 0 influx of magma into the deep magma chamber starts with a constant intensity Q in,d. Pressure in the deep chamber starts to grow leading to magma rising and opening of the pre-existing dyke through the rocks between two chambers. When the dyke reaches the shallow chamber, influx of the fresh magma leads to pressure growth and initiation process of the dome extrusion. Figs. 3a,b show the evolution of the chamber overpressures normalized to their maximum values (upper plots) and the influx into the shallow magma chamber Q in,s and discharge rate at the surface Q out (lower plots), for two cases: i) weak connectivity between the chambers (when the value of a 0 in Eq. (2) is small, i.e. 100 m) and ii) strong connectivity (when the value of a 0 in Eq. (2) is small, i.e. a 0 = 600 m).
9 In the case of weak connectivity the overpressure in the deep chamber increases to a high value (~70 MPa) and stays nearly constant while the overpressure in the shallow chamber shows a pronounced cyclic activity. In this case influx into the shallow chamber is almost constant although its variation is shifted with respect Q out. In the case of strong connectivity both chambers inflate and deflate simultaneously. Influx rate into the shallow chamber increases when the pressure in the shallow magma chamber decreases. Q out follows changes in Q in,s. For the intermediate values of a 0 the system stays in between these two end member cases. The case of strong connectivity is consistent with reconstructions of Foroozan et al. (2011). Behaviour for the weak connectivity is not supported by the recent analysis and interpretation of geodetic observations (Foroozan et al., 2011). In this case high chamber overpressure suggests a horizontal extension of the dyke and a consequent improvement of connectivity between two magma chambers. FIG. 3 HERE Fig 3. Time evolution of the chamber overpressures normalized to their maximum values in the upper plots, and of influx into the shallow magma chamber Q in,s (thicker lines) and discharge rate at the surface Q out (thinner lines) in the lower plots, for two cases: (a) weak connectivity between the chambers (a 0 = 100 m) and (b) strong connectivity (a 0 = 600 m). Π denotes the normalized to its maximum overpressure value for the deep magma chamber (thicker lines) and Π denotes the normalized overpressure for the shallow magma chamber (thinner lines). In order to quantify the inflation/deflation phase shift between two chambers we introduce a variable characterizing the value of offset as:,,, where T is the period of discharge rate variation. In the case when the two chambers work in phase the value of the offset is equal to 0 while an offset of 1 corresponds to completely opposite directions of chambers evolution. Fig. 4 shows the dependence of the offset on the intensity of deep magma chamber influx Q in,d. Curves are shown for the values of a 0 = 100 (short dashed line), 200, 300, 400, 500 and 600 m (solid line).
10 The offset progressively decreases with an increase of the horizontal dyke extent a 0 and increases with Q in,d because for larger influx intensity the pressure in the deep chamber reaches higher values (see Fig. 5) and the dyke cross-section area increases. This improves the connectivity between the chambers. FIG. 4 HERE Fig. 4. Dependence of the offset on the intensity of deep magma chamber influx Q in,d for the magma chamber volumes suggested by Hautmann et al. (2010), i.e. V ch,s = 4 km 3, V ch,d = 32 km 3, and different a 0. Fig. 5 shows the maximum overpressure P d that is reached in the deep chamber (a) and the amplitude of overpressure variation during the cycle ΔP d (b) as a function of Q in,d for different a 0. Small values of a 0 result in unrealistically high overpressures and small amplitude of pressure variations. However in real situations as the overpressure increases above a rock strength it is expected that a fracture will propagate at the tips of the dyke increasing the dyke dimension and moving the system to lower more realistic pressures. Increase in Q in,d leads to improvement of the connectivity between the chambers and thus amplitude of pressure variation. FIG. 5 HERE Fig. 5. Maximum overpressure P d that is reached in the deep chamber (a) and the amplitude of overpressure variation during the cycle ΔP d (b) as a function of Q in,d for V ch,s = 4 km 3, V ch,d = 32 km 3, as in Hautmann et al. (2010), and different a 0. The time lag between the beginning of the influx into the deep chamber and the start of the eruption is shown in Fig. 6 as a function of Q in,d for different a 0. It varies from ~3 months to more than a year. Larger influx intensities and better chamber connectivity lead in shortening of this interval. Mattioli et al. (2010) show that ground deformation and extrusion of the magma at the surface are closely correlated in time. This supports a hypothesis of good chamber connectivity for the SHV. FIG. 6 HERE
11 Fig. 6. Time lag between the beginning of the influx into the deep chamber and the start of the eruption for V ch,s = 4 km 3, V ch,d = 32 km 3, as in Hautmann et al. (2010), and different a Period of discharge rate variations In the case of weak connectivity between the chambers (Fig. 7) for fixed values of the shallow and deep magma chamber volumes (e.g., V ch,s = 4 km 3 and V ch,d = 32 km 3 ) the period of pulsations reaches a minimum value for some intermediate value of Q in,d corresponding to an unstable branch of the steady-state solution similar to the case of a single magma chamber (see Fig. 7 in Melnik and Sparks, 2005). When the connectivity is strong the period monotonically increases with increase in the influx rate into the deep chamber. In order to evaluate the role of magma chamber volumes on the period of discharge rate variation we performed a series of runs varying aspect ratios and volumes of both shallow (Fig. 8) and deep magma chambers (Fig. 9). FIG. 7 HERE Fig. 7. Period of pulsations as function of the influx into the deep chamber for V ch,s = 4 km 3, V ch,d = 32 km 3, as in Hautmann et al. (2010), and different a 0. Informed by the volume estimates of Paulatto et al. (2012), we varied the volume of the shallow magma chamber between 4 and 30 km 3 for different values of the deep magma chamber (from 30 to 60 km 3 ) fixing Q in,d = 2 m 3 /s, a value that represents a realistic estimation for the deep magma influx (Voight et al., 2010; Odbert and Wadge, 2010). Keeping fixed the volume of the deep chamber to V ch,d = 32 km 3, for V ch,s 20 km 3 the period monotonically increase from 125 to 150 days, but above 20 km 3 there is a bifurcation and the solution becomes multi-periodic. The case of V ch,s = 30 km 3 results in a double period behaviour with periods for small amplitude oscillations of order of 150 days whereas the large peak of increase in discharge rate occurs every 400 days (Fig. 8). Similar double periodicity behaviour was discovered first in Costa et al. (2007). FIG. 8 HERE
12 Fig 8. Change in the period of discharge rate variation with respect to the volume of the shallow magma chamber. In this case the volume of the deep magma chamber is fixed at 32 km 3 as estimated by Hautmann et al. (2010). The variation of discharge rate within one period of pulsation is shown in insets. Fig. 9 shows the dependence of the period and the shape of the pulsations as a function of the deep magma chamber volume. The volume of the shallow magma chamber is fixed at 15 km 3 in accord to Paulatto et al., (2012). Two shapes of the shallow magma chamber (spherical and prolate) are studied comparing results for aspect ratio of 1 and 1.5 respectively. For small volumes of the deep magma chamber (V ch,d < 37 km 3 ) the period does not strongly depend on the shallow chamber geometry and pulsations have a unimodal shape pattern. For larger magma chambers the pattern during the cycle consists of peaks in discharge rate separated by intervals of periodic variation in discharge rate with much lower amplitude that characterize major part of the cycle. FIG. 9 HERE Fig 9. Change in the period of discharge rate variation with respect to the volume of the deep magma chamber. The volume of the shallow magma chamber is fixed at 15 km 3. The patterns of variation of discharge rate within one period of pulsation is shown in the insets. We have discovered that a prolate shape of the shallow magma chamber leads to much longer periods of pulsations due to stronger influence of the chamber pressure on the opening of the dyke that connects it to the surface. The opposite effect was observed for the aspect ratio of deep magma chamber. A more spherical deep magma chamber tends to produce fluctuations with longer periodicity. Implications to SHV Observations at SHV suggest long-term cyclic behaviour with a period of 2 to 3 years. This period, for example, can be obtained for a prolate shallow magma chamber with a volume of ~15 km 3 and a deep magma chamber of ~50 km 3. Fig. 10 represents the variation in extrusion rate with time for these parameters.
13 Although with these magma chamber values we can obtain reasonable periodicity the simulated patterns of the flow rate variations, alternations between phases of high and low discharge rate, are not in agreement with observations for the SHV. Typical behaviour of the volcano involves repetitive periods of extrusion separated with repose periods with nearly no activity. Moreover, during the first three cycles the periods of extrusion were longer (~years) with respect the fourth and fifth cycle that show extrusion activity having a duration of months similar to the simulations showed in Fig. 10. The behaviour of the first three cycles is better reproduced assuming a narrow shallow dyke, d s (a 0 ~150m), as shown in Fig. 11. It is worth noting that this cyclic behaviour is closer to observations, and a systematic parametric investigation involving also the shallow part of the system is an object of ongoing research. It is also likely that magma rheology can play a key role on the behaviour of the extrusion pattern. For example, we would expect that for a Bingham rheology discharge rate between the two major pulses is zero until a critical chamber overpressure is reached (Melnik and Sparks, 2005). Moreover a more realistic description of the effective viscosity of magma during lava dome eruptions should account for the coupling with energy loss, viscous dissipation, two-dimensional and stick-slip effects (e.g., Costa and Macedonio, 2005; Costa et al. 2007c; 2012). However, due to the limitations of the current transient numerical code, we cannot investigate effect of complex rheologies. FIG. 10 HERE Fig. 10. Variation of discharge rate with time (V ch,s = 15 km 3, V ch,d = 50 km 3 ) for a shallow dyke width as in Costa et al. (2007a). FIG. 11 HERE Fig. 11. Variation of discharge rate with time (V ch,s = 15 km 3, V ch,d = 50 km 3 ) for a shallow dyke of 300m width. For comparison variation of discharge rate with time for a single chamber geometry and a shallow dyke of 300m (V ch,s = 30 km 3 ) is also shown.
14 Discussion and conclusion The models show that presence of two magma chambers connected by an elastic dyke significantly influences dynamics of magma extrusion in comparison with a case of a single magma chamber fed with a constant influx rate. The dynamics is strongly influenced by the amount of connectivity between the chambers measured as a length of large semi-axis of the dyke a 0. When a 0 is small (weak connectivity) the pressure in the deep chamber reaches high values and the influx into the shallow chamber is close to a constant value. This case, in terms of surface extrusion, is similar to the case of a single magma chamber but geodetic signal in this case will be different. At the initial stage of the eruption deep source will inflate producing wide spread deformation pattern. Later the pressure in the deep reservoir will come to a high constant value and ground deformation will be controlled by changes in pressure in the shallow chamber. In the case of strong connectivity (large a 0 ) initial ground deformation pattern will be controlled by an inflation of the deep magma chamber. Later in the eruptions both pressure sources work in phase. For the intermediate a 0 there is a transition between two end-members. It is clear that assumption of constant value of a 0 greatly oversimplifies the physics. In case of large overpressures, stress at dyke tips will exceed fracture toughness of the rocks and the dyke will expand horizontally. The maximum dyke extent will be controlled by horizontal chamber dimensions (Costa et al., 2011). When deep magma chamber deflates overpressure in the dyke will decrease and, as flow rate decreases, magma at dyke tips can solidify leading to a decrease in a 0 (Kavanagh and Sparks, 2011). Another strong simplification of the model comes from assumed Newtonian rheology of the magma. Although the system produces periodic behaviour the shape of discharge rate evolution during the cycle is different from observed at the SHV. Measurements of the dome volume suggest relatively slow initiation of the eruption with progressively increasing or nearly constant extrusion rate at later stages of the cycle followed by a rapid drop in discharge rate leading to stagnation of the eruption. Cycles produced by our model show rapid increase in discharge rate with more gradual reduction in eruption intensity. In the simulations with a shallow dyke of about 500m, as in Costa et al. (2007a), discharge rate remains always larger than zero
15 leading to continue extrusion of magma, whereas simulations using a narrower shallow dyke (~300m) show a better agreement. Finally, in order to understand if simulations can give insights about the issue single versus dual magma chamber system, we carried out also some simulations assuming only the presence of the shallow chamber. Results show that for a shallow dyke of ~300m we can have 2-3 periodicity assuming a chamber volume of about 30 km 3, implying that model results alone do not allow us to prefer a dual rather than a single chamber geometry. Predictions of the model must be verified against measurements of ground deformation, dome volume, gas flux and other geophysical data that provides constraints on the possible scenarios of the eruption. Only a robust approach that is based on an integrated system of geophysical observations will be able to give strong constraints to the model and give some insight into the geometry and dynamics of the system Acknowledgments. O.M. acknowledges the support from the ERC grant no (VOLDIES) and the Russian Foundation for Basic Research ( ). A.C. was supported by NERC research grant reference NE/H019928/1. We thank R.S.J. Sparks and G. Wadge for very useful comments on an earlier version of the manuscript
16 Appendix Governing equations in the shallow conduit (elliptical dyke to cylinder). The ascent of magma is simulated by 1-D transport equations. We assume that the conduit has an elliptical cross-section with area S = ab, with a and b major and minor semi-axes, respectively. Their vertical variation occurs at length-scales that are much larger than the dyke width. The set of cross-section averaged equations is presented below S t (Sρ ) + 1 m S x (Sρ V ) = G G m mc ph 1 S t (Sρ ) + 1 mc S x (Sρ V ) = G mc mc 1 S t (Sρ ) + 1 ph S x (Sρ V ) = G ph ph 1 S t (Sρ ) + 1 d S x (Sρ V ) = J d 1 S t (Sρ ) + 1 g S x (Sρ V ) = J g g (A1) (A2) (A3) (A4) (A5) Here t denotes time, x the vertical coordinate, ρ m, ρ ph, ρ mc, ρ d and ρ g are the densities of melt, phenocrysts, microlites, dissolved gas and exsolved gas respectively, and V and V g are the velocities of magma and gas, respectively. G ph, G mc represent the mass transfer rate due to crystallization of phenocrysts and microlites, respectively, and J
17 the mass transfer rate due to gas exsolution. Eq. (A1) represents the mass conservation for the melt phase, eqs. (A2) and (A3) are the conservation equations for microlites and phenocrysts respectively, eqs. (A4) and (A5) represent the conservation of the dissolved gas and of the exsolved gas respectively. Two momentum equations for the mixture as a whole (A6) and for the free gas phase (A7) describe the magma motion: p x = ρg 4μ a2 + b 2 a 2 b 2 V (A6) 501 V g V = k μ g p x (A7) Here p is the pressure, ρ the bulk density of magma, g the gravity acceleration, μ is the magma viscosity, k is the magma permeability and μ g is the gas viscosity. The main source of the magma temperature change is the release of the latent heat of crystallization (A8) S t (SρC mt ) + 1 S x (SρC mvt ) = L * (G mc + G ph ) (A8) Here C m is the bulk specific heat of magma, T is the bulk flow-averaged temperature, L * is latent heat of crystallization (viscous heating and heat flux to the host rocks are neglected). For details of equation of state and mass fluxes related to gas exsolution and crystallization see Costa et al. (2007a,b) For parameterizations of magma permeability k and magma viscosity μ we use:
18 k = k( α)= k 0 α j (A9) μ = μ m ( c,t)θ( β)η α,ca ( ) (A10) where k depends on bubble volume fraction of bubbles α, μ depends on water content c, temperature T, crystal content β, bubble fraction and capillary number Ca (Costa et al., 2007a,b). Variations of dyke semi-axis, a and b, with pressure are characterized by eq. (1). In order to get a smooth transition from the dyke at depth to a cylindrical conduit the value of a 0 is parameterized as: 524 x L a ( x) = A arctan + A T wt (A11) Here L T and w T are the height and the vertical extent of the cylinder to ellipse transition zone and constants A 1 and A 2 are calculated to satisfy conditions a 0 (L) = R and a 0 (0) = a 0, where R is the radius of the cylindrical part of the conduit and a 0 is the length of major dyke semi-axis at the inlet of the dyke. The value of b 0 is calculated in order to conserve the cross-section area of the unpressurized dyke, although it can also be specified independently References Costa, A., Macedonio G. (2005), Viscous heating in fluids with temperaturedependent viscosity: triggering of secondary flows, J. Fluid Mech., 540, Costa A., Melnik O., Sparks R.S.J. (2007a), Controls of conduit geometry and wallrock elasticity on lava dome eruptions. Earth Planet. Sci. Lett. 260: doi: /j.epsl
19 Costa A., Melnik O., Sparks R.S.J., Voight B. (2007b), The control of magma flow in dykes on cyclic lava dome extrusion. Geophys. Res. Lett. 34: L doi:1029/2006gl Costa, A., Melnik, O., Vedeneeva, E., (2007c), Thermal effects during magma ascent in conduits. J. Geophys. Res., 112, doi: /2007jb Costa, A., Sparks, R.S.J. Macedonio, G., Melnik, O. (2009), Effects of wall-rock elasticity on magma flow in dykes during explosive eruptions. Earth Planet. Sci. Lett. 288, , doi: /j.epsl Costa A., Gottsmann J., Melnik O., Sparks R.S.J. (2011), A stress-controlled mechanism for the intensity of very large magnitude explosive eruptions, Earth Planet. Sci. Lett., doi: /j.epsl Costa A., Wadge G., Melnik O. (2012), Cyclic extrusion of a lava dome based on a stick-slip mechanism. Earth Planet. Sci. Lett., in press Denlinger R.P., Hoblitt R.P. (1999), Cyclic eruptive behaviour of silicic volcanoes. Geology 27 (5): Dirksen O., Humphreys M.C.S., Pletchov P., Melnik O., Demyanchuk Y., Sparks R.S.J., Mahony S. (2006), The dome-forming eruption of Shiveluch Volcano, Kamchatka: Observation, petrological investigation and numerical modelling. J. Volcanol. Geotherm. Res. 155: doi: /j. jvolgeores Edmonds, M., M. Humphreys, E. Hauri, R.A. Herd, G. Wadge, H. Rawson, R. Ledden, M. Plail, J. Barclay, A. Aiuppa, T. Christopher, G. Giudice, R. Guida (2012), Pre-eruptive vapour and its role in controlling eruption style and longevity at Soufrière Hills Volcano, submitted to this issue. Elsworth, D., G. Mattioli, J. Taron, B. Voight, and R. Herd (2008), Implications of magma transfer between multiple reservoirs on eruption cycling, Science, 322(5899), , doi: /science
20 Foroozan, R., D. Elsworth, B. Voight, and G. Mattioli (2010), Dual reservoir structure at Soufriere Hills Volcano inferred from continuous GPS observations and heterogeneous elastic modeling, Geophys. Res. Lett., 37, L00E12, doi: /2010gl Foroozan, R., D. Elsworth, B. Voight, and G. S. Mattioli (2011), Magmatic metering controls the stopping and restarting of eruptions, Geophys. Res. Lett., 38, L05306, doi: /2010gl Hautmann S., Gottsmann J., Sparks R.S.J., Costa A., Melnik O., Voight B. (2009), Modelling ground deformation response to oscillating overpressure in a dyke conduit at Soufriere Hills Volcano, Montserrat, Tectonophysics, Vol. 471: 87-95, doi: /j.tecto Hautmann, S., Gottsmann J., Sparks R.S.J., Mattioli G.S., Sacks I.S., Strutt M.H. (2010), Effect of mechanical heterogeneity in arc crust on volcano deformation with application to Soufrière Hills Volcano, Montserrat, West Indies, J. Geophys. Res., 115, B09203, doi: /2009jb Kavanagh, J.L., R.S.J. Sparks (2011), Insights of dyke emplacement mechanics from detailed 3D dyke thickness datasets J. Geol. Soc., 168, , doi: / Lensky, N.G., Sparks. R.S.J., Navon. O. and Lyakhovsky, V. (2008), Cyclic activity at Soufrière Hills volcano, Montserrat. In Lane, S.J. & Gilbert, J.S., ((eds) Fluid Motions in Volcanic Conduits: A Source of Seismic and Acoustic Signals. Geological Society, London, Special Publications, 307, Linde, A. T., et al. (2010), Vulcanian explosion at Soufrière Hills Volcano, Montserrat on March 2004 as revealed by strain data, Geophys. Res. Lett., 37, L00E07, doi: / 2009GL Maeda I. (2000), Nonlinear visco-elastic volcanic model and its application to the recent eruption of Mt. Unzen. J. Volcanol. Geotherm. Res. 95: 35 47
21 Mattioli, G. S., R. A. Herd, M. H. Strutt, G. Ryan, C. Widiwijayanti, and B. Voight (2010), Long term surface deformation of Soufrière Hills Volcano, Montserrat from GPS geodesy: Inferences from simple elastic inverse models, Geophys. Res. Lett., 37, L00E13, doi: /2009gl Melnik, O., and R.S.J. Sparks, (2005), Controls on conduit magma flow dynamics during lava dome building eruptions, J. Geophys. Res., 110, B02209, doi: /2004jb Muskhelishvili, N. (1963), Some Basic Problems in the Mathematical Theory of Elasticity, Noordhof, Leiden, The Netherlands. Odbert, H.M., G. Wadge, (2009), Time series analysis of lava flux, J. Volcanol. Geotherm. Res., 188 (4), Paulatto, M., C. Annen, T.J. Henstock, E. Kiddle, T.A. Minshull, R.S.J. Sparks, B. Voight (2012). Magma chamber properties from integrated seismic tomography and thermal modelling at Montserrat, submitted to Geochem. Geophys. Geosyst. Roman, D. C., J. Neuberg, and R. R. Luckett (2006), Assessing the likelihood of volcanic eruption through analysis of volcanotectonic earthquake fault-plane solutions, Earth Planet. Sci. Lett., 248, , doi: /j.epsl Sacks, S., A. Linde, D. Hidayat (2011), Change in the Soufriere Hills Volcano magma system, Scientific Conference Soufriere Hills Volcano 15 years on, Montserrat, West Indies, 4-8 April Shalev, E., et al. (2010), Three dimensional seismic velocity tomography of Montserrat from the SEA CALIPSO offshore/onshore experiment, Geophys. Res. Lett., 37, L00E17, doi: /2010gl Voight, B., et al. (1999), Magma flow instability and cyclic activity at Soufriere Hills Volcano, Montserrat, Science, 283, Voight, B., C. Widiwijayanti, G. Mattioli, D. Elsworth, D. Hidayat, and M. Strutt (2010), Magma-sponge hypothesis and stratovolcanoes: Case for a compressible
22 reservoir and quasi-steady deep influx at Soufriere Hills Volcano, Montserrat, Geophys. Res. Lett., 37, L00E05, doi: /2009gl Wang, Z., (2000), Dynamic versus static elastic properties of reservoir rocks, SEG Books, Vol. 19, Pags Wylie J.J., Voight B., Whitehead J.A. (1999) Instability of magma flow from volatiledependent viscosity. Science 285:
23 FIGURE CAPTIONS Fig. 1 Sketch of the investigated system. Modified after Hautmann et al. (2010). Fig. 2. A set of steady state solutions: dependence of discharge rate on chamber pressure for three different far-field stresses: compressional (-10 MPa; dasheddotted line), neutral (full line), and extensional (10 MPa; dotted line). The remaining parameters are set as in Costa et al. (2007b). Fig 3. Time evolution of the chamber overpressures normalized to their maximum values in the upper plots, and of influx into the shallow magma chamber Q in,s (thicker lines) and discharge rate at the surface Q out (thinner lines) in the lower plots, for two cases: (a) weak connectivity between the chambers (a 0 = 100 m) and (b) strong connectivity (a 0 = 600 m). Π denotes the normalized to its maximum value overpressure for the deep magma chamber (thicker lines) and Π denotes the normalized overpressure for the shallow magma chamber (thinner lines). Fig. 4. Dependence of the offset on the intensity of deep magma chamber influx Q in,d for the magma chamber volumes suggested by Hautmann et al. (2010), i.e. V ch,s = 4 km 3, V ch,d = 32 km 3, and different a 0. Fig. 5. Maximum overpressure P d that is reached in the deep chamber (a) and the amplitude of overpressure variation during the cycle ΔP d (b) as a function of Q in,d for V ch,s = 4 km 3, V ch,d = 32 km 3, as in Hautmann et al. (2010), and different a 0. Fig. 6. Time lag between the beginning of the influx into the deep chamber and the start of the eruption for V ch,s = 4 km 3, V ch,d = 32 km 3, as in Hautmann et al. (2010), and different a 0.
24 Fig. 7. Period of pulsations as function of the influx into the deep chamber for V ch,s = 4 km 3, V ch,d = 32 km 3, as in Hautmann et al. (2010), and different a 0. Fig 8. Change in the period of discharge rate variation with respect to the volume of the shallow magma chamber. In this case the volume of the deep magma chamber is fixed at 32 km 3 as estimated by Hautmann et al. (2010). The variation of discharge rate within one period of pulsation is shown in insets. Fig 9. Change in the period of discharge rate variation with respect to the volume of the deep magma chamber. The volume of the shallow magma chamber is fixed at 15 km 3. The patterns of variation of discharge rate within one period of pulsation is shown in the insets. Fig. 10. Variation of discharge rate with time (V ch,s = 15 km 3, V ch,d = 50 km 3 ). Fig. 11. Variation of discharge rate with time (V ch,s = 15 km 3, V ch,d = 50 km 3 ) for a shallow dyke of 300m width. For comparison variation of discharge rate with time for a single chamber geometry and a shallow dyke of 300m (V ch,s = 30 km 3 ) is also shown.
25 Q out cylindrical part Dyke ds (connecting shallow magma chamber with the shallow cylindrical conduit) Dyke d d (connecting deep with shallow reservoir) P Q in,s Chamber Cs (shallow magma chamber centered at ~6 km depth, spherical geometry, volume V ch,s ) P Chamber C d (deep magma reservoir, centered at ~13 km depth, prolate geometry, volume V ch,d ) Q in,d (deep supply)
26 Discharge Rate (m 3 /s) Pressure, P (MPa) s σ =-10MPa ff σ = 0 MPa ff σ =10 MPa ff ΔP Deep Chamber = 10 MPa Extensional Compressional
27 Π Π 6 6 Q ins, Q out (m 3 /s) 4 Π Π Q ins, Q out (m 3 /s) time (days) time (days) (a) (b)
28 offset Q in,d (m 3 /s )
29 Δ (a) (b)
30 t beg (days) Q in,d (m 3 /s )
31 period (days) Q in,d (m 3 /s)
32 period (days) V ch,s (km 3 )
33 Aspect ratio period (days) V ch,d (km 3 )
34 8 6 Q out (m 3 /s) time (days)
35 12 1 chamber 2 chambers 8 Q out (m 3 /s) time (days)
Volcanic Eruptions: Cyclicity During Lava Dome Growth
Volcanic Eruptions: Cyclicity During Lava Dome Growth Oleg Melnik a,b, R. Stephen J. Sparks b, Antonio Costa b,c, Alexei A. Barmin a a Institute of Mechanics, Moscow State University, Russia b Earth Science
More informationUniversity of Bristol - Explore Bristol Research
Sheldrake, T. E., Aspinall, W. P., Odbert, H. M., Wadge, G., Sparks, R. S. J. (2017). Understanding causality and uncertainty in volcanic observations:: An example of forecasting eruptive activity on Soufrière
More informationJournal Publications (current reporting period , not reported below)
Journal Publications (2005 to 2011) A number of papers have been published related to CALIPSO, SEA-CALIPSO, and the CALIPSO Facility, including many recently. Seven papers were presented at the Soufriere
More informationChapter 12 Geodetic imaging of magma migration at Soufrière Hills Volcano 1995 to 2008
Downloaded from http://mem.lyellcollection.org/ by guest on June 3, 2014 Geological Society, London, Memoirs Chapter 12 Geodetic imaging of magma migration at Soufrière Hills Volcano 1995 to 2008 Derek
More informationDegassing during quiescence as a trigger of magma ascent and volcanic eruptions
Degassing during quiescence as a trigger of magma ascent and volcanic eruptions Társilo Girona 1,*, Fidel Costa 1,2, Gerald Schubert 3 1 Earth Observatory of Singapore, Nanyang Technological University,
More informationA simple formula for calculating porosity of magma in volcanic conduits during dome-forming eruptions
Earth Planets Space, 62, 483 488, 2010 A simple formula for calculating porosity of magma in volcanic conduits during dome-forming eruptions Tomofumi Kozono 1 and Takehiro Koyaguchi 2 1 National Research
More informationERTH 456 / GEOL 556 Volcanology. Lecture 06: Conduits
1 / 28 ERTH 456 / GEOL 556 Volcanology Lecture 06: Conduits Ronni Grapenthin rg@nmt.edu MSEC 356, x5924 hours: TR 3-4PM or appt. September 12, 2016 2 / 28 How does magma get from source to surface? What
More informationPUBLICATIONS. Geochemistry, Geophysics, Geosystems
PUBLICATIONS Geochemistry, Geophysics, Geosystems RESEARCH ARTICLE Key Points: Strain data analysis from explosions/ degassing events at Soufriere Hills Volcano Pressure release deep within the magmatic
More informationUnique strainmeter observations of Vulcanian explosions, Soufrière Hills Volcano, Montserrat, July 2003
Click Here for Full Article GEOPHYSICAL RESEARCH LETTERS, VOL. 37,, doi:10.1029/2010gl042551, 2010 Unique strainmeter observations of Vulcanian explosions, Soufrière Hills Volcano, Montserrat, July 2003
More informationHigh-temperature fracture of magma
High-temperature fracture of magma Hugh Tuffen Peter Sammonds Rosanna Smith Harry Pinkerton Don Dingwell Jon Castro Cracks, Fractures and Faults in the Earth Thursday 19 th June 2008 Montserrat (Sparks
More informationEffect of mechanical heterogeneity in arc crust on volcano deformation with application to Soufrière Hills Volcano, Montserrat, West Indies
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 115,, doi:10.1029/2009jb006909, 2010 Effect of mechanical heterogeneity in arc crust on volcano deformation with application to Soufrière Hills Volcano, Montserrat,
More informationAdvanced Workshop on Evaluating, Monitoring and Communicating Volcanic and Seismic Hazards in East Africa.
2053-11 Advanced Workshop on Evaluating, Monitoring and Communicating Volcanic and Seismic Hazards in East Africa 17-28 August 2009 Seismic monitoring on volcanoes in a multi-disciplinary context Jürgen
More informationGeophysical Research Letters
Geophysical Research Letters 16 FEBRUARY 2006 Volume 33 Number 3 American Geophysical Union Sudden pressure increase in magma reservoir triggered by gigantic lava dome collapse on Montserrat History of
More informationThis document is downloaded from DR-NTU, Nanyang Technological University Library, Singapore.
Title Author(s) Citation Explosion dynamics from strainmeter and microbarometer observations, Soufrière Hills Volcano, Montserrat : 2008 2009. Chardot, L.; Voight, Barry.; Foroozan, R.; Sacks, S.; Linde,
More informationControls on conduit magma flow dynamics during lava dome building eruptions
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 110,, doi:10.1029/2004jb003183, 2005 Controls on conduit magma flow dynamics during lava dome building eruptions O. Melnik 1 and R. S. J. Sparks Research Centre for
More informationInfluence of conduit flow mechanics on magma rheology and the growth style of lava domes
Influence of conduit flow mechanics on magma rheology and the growth style of lava domes Taha Husain, Derek Elsworth Energy and Mineral Engineering, G3 Center and EMS Energy Institute, The Pennsylvania
More informationVolcano Seismicity and Tremor. Geodetic + Seismic
Volcano Seismicity and Tremor Seismic Imaging Geodetic + Seismic Model based joint inversion Geodetic Monitoring How is magma stored in the crust? geometry, volume and physical state of crustal melts.
More informationIncorporating seismic observations into 2D conduit flow modeling
Journal of Volcanology and Geothermal Research 152 (2006) 331 346 www.elsevier.com/locate/jvolgeores Incorporating seismic observations into 2D conduit flow modeling L. Collier, J. Neuberg * Department
More informationOn the feasibility of magma fracture within volcanic conduits: Constraints from earthquake data and empirical modelling of magma viscosity
GEOPHYSICAL RESEARCH LETTERS, VOL. 38,, doi:10.1029/2011gl049297, 2011 On the feasibility of magma fracture within volcanic conduits: Constraints from earthquake data and empirical modelling of magma viscosity
More informationTopics. Magma Ascent and Emplacement. Magma Generation. Magma Rise. Energy Sources. Instabilities. How does magma ascend? How do dikes form?
Magma Ascent and Emplacement Reading: Encyclopedia of Volcanoes: Physical Properties of Magmas (pp. 171-190) Magma Chambers (pp. 191-206) Plumbing Systems (pp. 219-236) Magma ascent at shallow levels (pp.237-249)
More informationThis document is downloaded from DR-NTU, Nanyang Technological University Library, Singapore.
This document is downloaded from DR-NTU, Nanyang Technological University Library, Singapore. Title Magma degassing and intermittent lava dome growth. Author(s) Taisne, B.; Jaupart, C. Citation Taisne,
More informationVolcanic Eruptions and Hydrovolcanism
Find the Face Volcanic Eruptions and Hydrovolcanism Ocean Ridges Continental Rifts Subduction Zones: Continental Oceanic Back Arc Basins Hot Spots Plumes Cinder Cones Composite Volcanoes Shield VolcanoesCinder
More informationAbout Some Features of a Magma Flow Structure at Explosive Volcano Eruptions
About Some Features of a Magma Flow Structure at Explosive Volcano Eruptions V. Kedrinskiy 1 Introduction The cyclic character of magma ejections is one of the basic aspects in the research field of the
More informationSupplementary material
1 GSA DATA REPOSITORY 2011128 Melnik et al. Supplementary material 1. CSD acquisition: The groundmass plagioclase CSD (Fig. 1a and DR1d) was generated from element maps of 1024 by 800 points over an area
More informationVolcano inflation prior to an eruption: Numerical simulations based on a 1-D magma flow model in an open conduit
Earth Planets Space, 65, 1477 1489, 2013 Volcano inflation prior to an eruption: Numerical simulations based on a 1-D magma flow model in an open conduit Ryohei Kawaguchi, Takeshi Nishimura, and Haruo
More informationGravity Tectonics Volcanism Atmosphere Water Winds Chemistry. Planetary Surfaces
Gravity Tectonics Volcanism Atmosphere Water Winds Chemistry Planetary Surfaces Gravity & Rotation Polar flattening caused by rotation is the largest deviation from a sphere for a planet sized object (as
More informationThe Afar Rift Consortium
The Afar Rift Consortium University of Bristol University of Oxford University of Cambridge British Geological Survey University of Rochester University of Auckland, NZ Opportunistic programme following
More informationViscosity of magmas containing highly deformable bubbles
Journal of Volcanology and Geothermal Research 105 (2001) 19±24 www.elsevier.nl/locate/jvolgeores Viscosity of magmas containing highly deformable bubbles M. Manga a, *, M. Loewenberg b a Department of
More informationDegassing processes and recent activity at Volcán de Colima. Universidad de Colima, Mexico Corresponding author
Degassing processes and recent activity at Volcán de Colima Nick Varley * & Gabriel Reyes Dávila Universidad de Colima, Mexico Corresponding author email: nick@ucol.mx Volcán de Colima is currently in
More informationFracture of andesite in the brittle and brittleductile transition regimes
Fracture of andesite in the brittle and brittleductile transition regimes R. Smith 1,*, P.R. Sammonds 1, C.R.J. Kilburn 1, H. Tuffen 2,1 1. Department of Earth Sciences, UCL, Gower Street, London, WC1E
More informationInfluence of pre-existing volcanic edifice geometry on caldera formation
1 2 Author manuscript, published in "Geophysical Research Letters 38 (2011) L11305" DOI : 101029/2011gl047900 Influence of pre-existing volcanic edifice geometry on caldera formation V Pinel, 1 V Pinel,
More informationWhy do volcanoes (only sometimes) erupt explosively?
Why do volcanoes (only sometimes) erupt explosively? 2004-2008, effusive 1980, explosive Michael Manga Gonnermann and Manga, Annual Reviews of Fluids Mechanics, 2007 Why do volcanoes erupt explosively?
More informationThree dimensional seismic velocity tomography of Montserrat from the SEA CALIPSO offshore/onshore experiment
Click Here for Full Article GEOPHYSICAL RESEARCH LETTERS, VOL. 37,, doi:10.1029/2010gl042498, 2010 Three dimensional seismic velocity tomography of Montserrat from the SEA CALIPSO offshore/onshore experiment
More informationPhysics based models of ground deformation and extrusion rate at effusively erupting volcanoes
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 116,, doi:10.1029/2010jb007939, 2011 Physics based models of ground deformation and extrusion rate at effusively erupting volcanoes Kyle Anderson 1 and Paul Segall
More informationDynamic analysis. 1. Force and stress
Dynamic analysis 1. Force and stress Dynamics is the part of structural geology that involves energy, force, stress, and strength. It's very important to distinguish dynamic concepts from kinematic ones.
More informationvolcanic tremor and Low frequency earthquakes at mt. vesuvius M. La Rocca 1, D. Galluzzo 2 1
volcanic tremor and Low frequency earthquakes at mt. vesuvius M. La Rocca 1, D. Galluzzo 2 1 Università della Calabria, Cosenza, Italy 2 Istituto Nazionale di Geofisica e Vulcanologia Osservatorio Vesuviano,
More informationHow does magma reach the surface?
How does magma reach the surface? 2004-2008, effusive 1980, explosive Michael Manga Why do volcanoes (only sometimes) erupt explosively? 2004-2008, effusive 1980, explosive Michael Manga Gonnermann and
More informationEarthquake and Volcano Clustering at Mono Basin (California)
Excerpt from the Proceedings of the COMSOL Conference 2010 Paris Earthquake and Volcano Clustering at Mono Basin (California) D. La Marra *,1, A. Manconi 2,3 and M. Battaglia 1 1 Dept of Earth Sciences,
More informationJournal of Volcanology and Geothermal Research
Journal of Volcanology and Geothermal Research 285 (2014) 100 117 Contents lists available at ScienceDirect Journal of Volcanology and Geothermal Research journal homepage: www.elsevier.com/locate/jvolgeores
More informationPre-explosive conduit conditions of the 1997 Vulcanian explosions at Soufrière Hills Volcano, Montserrat: II. Overpressure and depth distributions
Pre-explosive conduit conditions of the 1997 Vulcanian explosions at Soufrière Hills Volcano, Montserrat: II. Overpressure and depth distributions Alain Burgisser, Laurent Arbaret, Timothy H. Druitt, Thomas
More informationAnalytical model of surface uplift above axisymmetric flat-lying magma intrusions: Implications for sill emplacement and geodesy
Author manuscript, published in "Journal of Volcanology and Geothermal Research 253 (2013) 114-130" DOI : 10.1016/j.jvolgeores.2012.12.006 Analytical model of surface uplift above axisymmetric flat-lying
More informationNumerical Simulation of magma plumbing system associated with the eruption at the Showa crater of Sakurajima inferred from ground deformation
Numerical Simulation of magma plumbing system associated with the eruption at the Showa crater of Sakurajima inferred from ground deformation Soma Minami 1, Masato Iguchi 2, Hitoshi Mikada 3, Tada-nori
More informationThermomechanical controls on magma supply and volcanic deformation: application to Aira caldera, Japan
SUPPLEMETARY MATERIAL FOR: Thermomechanical controls on magma supply and volcanic deformation: application to Aira caldera, Japan James Hickey 1*, Joachim Gottsmann 1, Haruhisa akamichi 2 & Masato Iguchi
More informationSurface changes caused by erosion and sedimentation were treated by solving: (2)
GSA DATA REPOSITORY 214279 GUY SIMPSON Model with dynamic faulting and surface processes The model used for the simulations reported in Figures 1-3 of the main text is based on two dimensional (plane strain)
More informationUsing impact exsolution to link the Chicxulub collision and Deccan volcanism
Using impact exsolution to link the Chicxulub collision and Deccan volcanism Kevin G. Harrison 1 1 Geosciences Department, 100 West College Street, PO Box 810, Denison University, Granville, OH 43023-0810
More informationChapter 5. Numerical Simulations of Volcanic Eruptions
Chapter 5 Numerical Simulations of Volcanic Eruptions Chapter Contents Chapter Contents 5.1. Previous Considerations 5.1.1. Objectives 5.1.. Geometry and Boundary Conditions 5.1.3. A Relation between W
More informationIntroduction to Volcanic Seismology
Introduction to Volcanic Seismology Second edition Vyacheslav M. Zobin Observatorio Vulcanolo'gico, Universidad de Colima, Colima, Col., Mexico ELSEVIER AMSTERDAM BOSTON HEIDELBERG LONDON * NEW YORK OXFORD
More informationViscous dissipation and temperature dependent shear thinning rheology
Viscous dissipation and temperature dependent shear thinning rheology M. Filippucci 1, A. Tallarico 1, M. Dragoni 2 1 Dipartimento di Scienze della Terra e Geoambientali, Università degli Studi di Bari
More informationThe role of bubble formation in volcanic eruption
The role of bubble formation in volcanic eruption Eyal Goldmann Division of Natural Sciences, El Camino College, Torrance, CA, 90506 Prepared for Geology 1 at El Camino College Fall 2009 1 1. Introduction
More informationVolcano-tectonic interaction at Soufrière Hills volcano, Montserrat (W.I.), constrained by dynamic gravity data
IOP Conference Series: Earth and Environmental Science Volcano-tectonic interaction at Soufrière Hills volcano, Montserrat (W.I.), constrained by dynamic gravity data To cite this article: Stefanie Hautmann
More informationVolcanic Plumes. JOHN WILEY & SONS Chichester New York Weinheim Brisbane Singapore Toronto
Volcanic Plumes R. S. J. SPARKS University of Bristol, UK M. I. BURSIK State University of New York, USA S. N. CAREY University of Rhode Island, USA J. S. GILBERT Lancaster University, UK L. S. GLAZE NASA/Goddard
More informationModeling seismic wave propagation during fluid injection in a fractured network: Effects of pore fluid pressure on time-lapse seismic signatures
Modeling seismic wave propagation during fluid injection in a fractured network: Effects of pore fluid pressure on time-lapse seismic signatures ENRU LIU, SERAFEIM VLASTOS, and XIANG-YANG LI, Edinburgh
More informationSupplementary material - Physical constraints for effective. magma-water interaction along volcanic conduits during silicic
GSA Data Repository 2018319 https://doi.org/10.1130/g45065.1. Supplementary material - Physical constraints for effective magma-water interaction along volcanic conduits during silicic explosive eruptions
More informationDurham Research Online
Durham Research Online Deposited in DRO: 31 March 2015 Version of attached le: Accepted Version Peer-review status of attached le: Peer-reviewed Citation for published item: Humphreys, Madeleine C.S. and
More informationEarthquake and Volcano Deformation
Earthquake and Volcano Deformation Paul Segall Stanford University Draft Copy September, 2005 Last Updated Sept, 2008 COPYRIGHT NOTICE: To be published by Princeton University Press and copyrighted, c
More informationGEODETIC EVIDENCE FOR SLOW INFLATION OF THE SANTORINI CALDERA
Proceedings, 11 th FIG Symposium on Deformation Measurements, Santorini, Greece, 2003. GEODETIC EVIDENCE FOR SLOW INFLATION OF THE SANTORINI CALDERA Stathis Stiros 1, Aris Chasapis 2 and Villy Kontogianni
More informationUniversity of Twente. Faculty of Mathematical Sciences. Two-dimensional magma-repository interactions. Memorandum No O. Bokhove.
Faculty of Mathematical Sciences t University of Twente The Netherlands P.O. Box 17 75 AE Enschede The Netherlands Phone: +31-53-893 Fax: +31-53-89311 Email: memo@math.utwente.nl www.math.utwente.nl/publications
More informationWhat is the Relationship between Pressure & Volume Change in a Magma Chamber and Surface Deformation at Active Volcanoes?
SSAC-pv2007.QE522.PL1.1 What is the Relationship between Pressure & Volume Change in a Magma Chamber and Surface Deformation at Active Volcanoes? What factors control the magnitude of surface deformation?
More informationAn Investigation on the Effects of Different Stress Regimes on the Magnitude Distribution of Induced Seismic Events
An Investigation on the Effects of Different Stress Regimes on the Magnitude Distribution of Induced Seismic Events Afshin Amini, Erik Eberhardt Geological Engineering, University of British Columbia,
More informationTWO-DIMENSIONAL MAGMA FLOW *
Iranian Journal of Science & Technology, Transaction A, Vol. 34, No. A2 Printed in the Islamic Republic of Iran, 2010 Shiraz University TWO-DIMENSIONAL MAGMA FLOW * A. MEHMOOD 1** AND A. ALI 2 1 Department
More informationFRICTIONAL HEATING DURING AN EARTHQUAKE. Kyle Withers Qian Yao
FRICTIONAL HEATING DURING AN EARTHQUAKE Kyle Withers Qian Yao Temperature Change Along Fault Mode II (plain strain) crack rupturing bilaterally at a constant speed v r Idealize earthquake ruptures as shear
More informationEAS 116 Earthquakes and Volcanoes
EAS 116 Earthquakes and Volcanoes J. Haase Forecasting Volcanic Eruptions Assessment of Volcanic Hazard Is that volcano active? Mount Lassen: 12000 BP and 1915 Santorini, IT: 180,000 BP, 70,000 BP, 21000
More informationSource parameters II. Stress drop determination Energy balance Seismic energy and seismic efficiency The heat flow paradox Apparent stress drop
Source parameters II Stress drop determination Energy balance Seismic energy and seismic efficiency The heat flow paradox Apparent stress drop Source parameters II: use of empirical Green function for
More informationVolcanoes in Compressional Settings (a seismological perspective)
Volcanoes in Compressional Settings (a seismological perspective) Diana C. Roman Department of Terrestrial Magnetism Carnegie Institution for Science December 11, 2016 AGU 2016 GeoPRISMS Mini-Workshop
More informationThe influence of short wavelength variations in viscosity on subduction dynamics
1 Introduction Deformation within the earth, driven by mantle convection due primarily to cooling and subduction of oceanic lithosphere, is expressed at every length scale in various geophysical observations.
More informationConstitution of Magmas. Magmas. Gas Law. Composition. Atomic Structure of Magma. Structural Model. PV = nrt H 2 O + O -2 = 2(OH) -
Constitution of Magmas Magmas Best, Ch. 8 Hot molten rock T = 700-1200 degrees C Composed of ions or complexes Phase Homogeneous Separable part of the system With an interface Composition Most components
More informationPRESSURE PERTURBATIONS IN TWO PHASE GEOTHERMAL RESERVOIR ASSOCIATED WITH SEISMISITY
PROCEEDINGS, Twenty-Eighth Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, January 27-29, 2003 SGP-TR-173 PRESSURE PERTURBATIONS IN TWO PHASE GEOTHERMAL RESERVOIR
More informationCOMPOSITION and PHYSICAL PROPERTIES GENERAL SUBJECTS. GEODESY and GRAVITY
COMPOSITION and PHYSICAL PROPERTIES Composition and structure of the continental crust Composition and structure of the core Composition and structure of the mantle Composition and structure of the oceanic
More informationOriginally published as:
Originally published as: Pascal, K., Neuberg, J., Rivalta, E. (2014): On precisely modelling surface deformation due to interacting magma chambers and dykes. - Geophysical Journal International, 196, 1,
More informationTransition in eruption style during the 2011 eruption of Shinmoe-dake, in the Kirishima volcanic group: Implications from a steady conduit flow model
Earth Planets Space, 65, 645 655, 2013 Transition in eruption style during the 2011 eruption of Shinmoe-dake, in the Kirishima volcanic group: Implications from a steady conduit flow model Ryo Tanaka and
More informationEVALUATING HEAT FLOW AS A TOOL FOR ASSESSING GEOTHERMAL RESOURCES
PROCEEDINGS, Thirtieth Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, January 31-February 2, 2005 SGP-TR-176 EVALUATING HEAT FLOW AS A TOOL FOR ASSESSING GEOTHERMAL
More informationGS388 Handout: Radial density distribution via the Adams-Williamson equation 1
GS388 Handout: Radial density distribution via the Adams-Williamson equation 1 TABLE OF CONTENTS ADIABATIC COMPRESSION: THE ADAMS WILLIAMSON EQUATION...1 EFFECT OF NON-ADIABATIC TEMPERATURE GRADIENT...3
More informationMEASUREMENT OF SURFACE DEFORMATION ASSOCIATED WITH HYDROFRACTURE
Proceedings, 11 th FIG Symposium on Deformation Measurements, Santorini, Greece, 2003. MEASUREMENT OF SURFACE DEFORMATION ASSOCIATED WITH HYDROFRACTURE Bassam Saleh 1, Balqies Sadoun 1 and Pierre Antoine
More informationDike induced Faults. Thomas Heinig. TU Bergakademie Freiberg, B. v. Cotta-Straße 2
Dike induced Faults Thomas Heinig TU Bergakademie Freiberg, B. v. Cotta-Straße 2 Abstract. The preferred region of dike induced faults can be found in rifts. They develop because of the extensional forces
More informationClock Advance in Earthquake Triggering of Volcanic Eruptions
Clock Advance in Earthquake Triggering of Volcanic Eruptions Mark Bebbington 1 and Warner Marzocchi 2 1 Volcanic Risk Solutions, Massey University, Palmerston North, New Zealand 2 Istituto Nazionale di
More informationFundamentals of Fluid Dynamics: Elementary Viscous Flow
Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research
More informationAnnual Project Report Project Title: Collaborative Research: The CALIPSO Project: Imaging the Magma Chamber on Montserrat
Annual Project Report Project Title: Collaborative Research: The CALIPSO Project: Imaging the Magma Chamber on Montserrat Official Reporting Period: 09/2010 to 08/2011 PI: G. Mattioli Original Awardee:
More informationVolcanic gas emissions from Soufrière Hills Volcano, Montserrat , with implications for mafic magma supply and degassing
Click Here for Full Article GEOPHYSICAL RESEARCH LETTERS, VOL. 37, L00E04, doi:10.1029/2009gl041325, 2010 Volcanic gas emissions from Soufrière Hills Volcano, Montserrat 1995 2009, with implications for
More informationTensor character of pore pressure/stress coupling in reservoir depletion and injection
Tensor character of pore pressure/stress coupling in reservoir depletion and injection Müller, B., Altmann, J.B., Müller, T.M., Weißhardt, A., Shapiro, S., Schilling, F.R., Heidbach, O. Geophysical Institute
More information3D simulations of an injection test done into an unsaturated porous and fractured limestone
3D simulations of an injection test done into an unsaturated porous and fractured limestone A. Thoraval *, Y. Guglielmi, F. Cappa INERIS, Ecole des Mines de Nancy, FRANCE *Corresponding author: Ecole des
More informationLateral extrusion and tectonic escape in Ilan Plain of northeastern Taiwan
Lateral extrusion and tectonic escape in Ilan Plain of northeastern Taiwan Angelier, J., Chang, T.Y., Hu, J.C., Chang, C.P., Siame, L., Lee, J.C., Deffontaines, B., Chu, H.T, Lu, C.Y., Does extrusion occur
More informationAnalytical Stress Modeling for Mine related Microseismicity
Analytical Stress Modeling for Mine related Microseismicity Himanshu Barthwal and Mirko van der Baan University of Alberta Summary Microseismicity is recorded in an underground mine by a network of 7 boreholes
More informationSupplementary information on the West African margin
Huismans and Beaumont 1 Data repository Supplementary information on the West African margin Interpreted seismic cross-sections of the north Angolan to south Gabon west African passive margins 1-3, including
More informationDynamics of two-phase conduit flow of high-viscosity gas-saturated magma: large variations of sustained explosive eruption intensity
Bull Volcanol 2000) 62 : 153±170 Springer-Verlag 2000 RESEARCH ARTICLE O. Melnik Dynamics of two-phase conduit flow of high-viscosity gas-saturated magma: large variations of sustained explosive eruption
More informationfriction friction a-b slow fast increases during sliding
µ increases during sliding faster sliding --> stronger fault --> slows sliding leads to stable slip: no earthquakes can start velocity-strengthening friction slow fast µ velocity-strengthening friction
More informationAnalytical and Numerical Investigations on the Vertical Seismic Site Response
Analytical and Numerical Investigations on the Vertical Seismic Site Response Bo Han, Lidija Zdravković, Stavroula Kontoe Department of Civil and Environmental Engineering, Imperial College, London SW7
More informationEarthquakes & Volcanoes
Earthquakes & Volcanoes Geology - the study of solid Earth, the rocks of which it is composed, and the processes by which they change geo = Earth; ology = study of Earth s Layers Plate Tectonics - the
More informationAnswer ALL questions in Section A, and TWO questions from Section B.
UNIVERSITY OF EAST ANGLIA School of Environmental Sciences Main Series Undergraduate Examination 2012-2013 Candidate s no.: GEODYNAMICS: EARTH S ENGINE ENV-2A43 Time allowed: 2 hours. Answer ALL questions
More informationEvolution of the 2007 earthquake swarm, Tanzania: Envisat and ALOS interferometry, ground observations and elastic modeling
Evolution of the 2007 earthquake swarm, Tanzania: Envisat and ALOS interferometry, ground observations and elastic modeling Gidon Baer, Yariv Hamiel, Gadi Shamir, Ran Nof Geological Survey of Israel East
More informationMIGRATING SWARMS OF BRITTLE-FAILURE EARTHQUAKES IN THE LOWER CRUST BENEATH MAMMOTH MOUNTAIN, CALIFORNIA
MIGRATING SWARMS OF BRITTLE-FAILURE EARTHQUAKES IN THE LOWER CRUST BENEATH MAMMOTH MOUNTAIN, CALIFORNIA David Shelly and David Hill GRL, October 2011 Contents Tectonic Setting Long Valley Caldera Mammoth
More informationEART162: PLANETARY INTERIORS
EART162: PLANETARY INTERIORS Francis Nimmo Last Week Global gravity variations arise due to MoI difference (J 2 ) We can also determine C, the moment of inertia, either by observation (precession) or by
More information7.2.1 Seismic waves. Waves in a mass- spring system
7..1 Seismic waves Waves in a mass- spring system Acoustic waves in a liquid or gas Seismic waves in a solid Surface waves Wavefronts, rays and geometrical attenuation Amplitude and energy Waves in a mass-
More informationFrictional rheologies have a wide range of applications in engineering
A liquid-crystal model for friction C. H. A. Cheng, L. H. Kellogg, S. Shkoller, and D. L. Turcotte Departments of Mathematics and Geology, University of California, Davis, CA 95616 ; Contributed by D.
More informationHeidy Mader, Michael Burton, Margherita Polacci. A combined physico-chemical model for passive degassing from Stromboli volcano.
Heidy Mader, Michael Burton, Margherita Polacci A combined physico-chemical model for passive degassing from Stromboli volcano. 2 Degassing at Stromboli Degassing from basaltic volcanoes is the most common
More informationAppendix: Nomenclature
Appendix: Nomenclature 209 Appendix: Nomenclature Chapter 2: Roman letters A amplitude of surface wave m A j (complex) amplitude of tidal component [various] a right ascension rad or a magnitude of tidal
More informationA possible model for initiation of ULF oscillation in magma
ANNALS OF GEOPHYSICS, VOL. 47, N. 1, February 004 A possible model for initiation of ULF oscillation in magma Yuri A. Kopytenko and Lidia V. Nikitina St. Petersburg Department of Institute of Terrestrial
More informationEART193 Planetary Capstone. Francis Nimmo
EART193 Planetary Capstone Francis Nimmo Last Time silicate melting How and why are melts generated? Increase in mantle potential temperature; or Reduction in solidus temperature (e.g. water); or Thinning
More information3D Finite Element Modeling of fault-slip triggering caused by porepressure
3D Finite Element Modeling of fault-slip triggering caused by porepressure changes Arsalan Sattari and David W. Eaton Department of Geoscience, University of Calgary Suary We present a 3D model using a
More informationWET EXPLOSIVE ERUPTIONS. Hawaii Photograph: Dorian Weisel
WET EXPLOSIVE ERUPTIONS Hawaii Photograph: Dorian Weisel WET EXPLOSIVE ERUPTIONS mechanisms hot magma/ hot rock + water pyroclasts + steam rapid expansion of gas fragmentation of magma + wall rock external
More informationExercise: concepts from chapter 8
Reading: Fundamentals of Structural Geology, Ch 8 1) The following exercises explore elementary concepts associated with a linear elastic material that is isotropic and homogeneous with respect to elastic
More information