The influence of viscous and latent heating on crystal-rich magma flow in a conduit

The flow dynamics of crystal-rich high-viscosity magma is likely to be strongly influenced by viscous and latent heat release. Viscous heating is observed to play an important role in the dynamics of fluids with temperature-dependent viscosities. The growth of microlite crystals and the accompanying release of latent heat should play a similar role in raising fluid temperatures. Earlier models of viscous heating in magmas have shown the potential for unstable (thermal runaway) flow as described by a Gruntfest number, using an Arrhenius temperature dependence for the viscosity, but have not considered crystal growth or latent heating. We present a theoretical model for magma flow in an axisymmetric conduit and consider both heating effects using Finite Element Method techniques. We consider a constant mass flux in a 1-D infinitesimal conduit segment with isothermal and adiabatic boundary conditions and Newtonian and non-Newtonian magma flow properties. We find that the growth of crystals acts to stabilize the flow field and make the magma less likely to experience a thermal runaway. The additional heating influences crystal growth and can counteract supercooling from degassing-induced crystallization and drive the residual melt composition back towards the liquidus temperature. We illustrate the models with results generated using parameters appropriate for the andesite lava dome-forming eruption at Soufriere Hills Volcano, Montserrat. These results emphasize the radial variability of the magma. Both viscous and latent heating effects are shown to be capable of playing a significant role in the eruption dynamics of Soufriere Hills Volcano. Latent heating is a factor in the top two kilometres of the conduit and may be responsible for relatively short-term (days) transients. Viscous heating is less restricted spatially, but because thermal runaway requires periods of hundreds of days to be achieved, the process is likely to be interrupted. Our models show that thermal evolution of the conduit walls could lead to an increase in the effective diameter of flow and an increase in flux at constant magma pressure.

Pb-210-Ra-226 and Ra-228-Th-232 systematics in young arc lavas: implications for magma degassing and ascent rates

New data show that island arc rocks have (Pb-210/Ra-226)(o) ratios which range from as low as 0.24 up to 2.88. In contrast, (Ra-22S/Th-232) appears always within error of I suggesting that the large Ra-226-excesses observed in arc rocks were generated more than 30 years ago. This places a maximum estimate on melt ascent velocities of around 4000 m/year and provides further confidence that the Ra-226 excesses reflect deep (source) processes rather than shallow level alteration or seawater contamination. Conversely, partial melting must have occurred more than 30 years prior to eruption. The Pb-210 deficits are most readily explained by protracted magma degassing. Using published numerical models, the data suggest that degassing occurred continuously for periods up to several decades just prior to eruption but no link with eruption periodicity was found. Longer periods are required if degassing is discontinuous, less than 100% efficient or if magma is recharged or stored after degassing. The long durations suggest much of this degassing occurs at depth with implications for the formation of hydrothermal and copper-porphyry systems. A suite of lavas erupted in 1985-1986 from Sangeang Api volcano in the Sunda arc are characterised by deficits of Pb-210 relative to Ra-226 from which 6-8 years of continuous Rn-222 degassing would be inferred from recent numerical models. These data also form a linear (Pb-210)/Pb-(Ra-226)/Pb array which might be interpreted as a 71-year isochron. However, the array passes through the origin suggesting displacement downwards from the equiline in response to degassing and so the slope of the array is inferred not to have any age significance. Simple modelling shows that the range of (Ra-226)/Pb ratios requires thousands of years to develop consistent with differentiation occurring in response to cooling at the base of the crust. Thus, degassing post-dated, and was not responsible for magma differentiation. The formation, migration and extraction of gas bubbles must be extremely efficient in mafic magma whereas the higher viscosity of more siliceous magmas retards the process and can lead to Pb-210 excesses. A possible negative correlation between (Pb-210/Ra-226)(o) and SO2 emission rate requires further testing but may have implications for future eruptions. (C) 2004 Elsevier B.V. All rights reserved.

Tracing pre-eruptive magma degassing using (Pb-210/Ra-226) disequilibria in the volcanic deposits of the 1980-1986 eruption of Mount St. Helens

Disequilibria between Pb-210 and Ra-226 can be used to trace magma degassing, because the intermediate nuclides, particularly Rn-222, are volatile. Products of the 1980-1986 eruptions of Mount St. Helens have been analysed for (Pb-210/Ra-226). Both excesses and deficits of Pb-210 are encountered suggesting rapid gas transfer. The time scale of diffuse, non-eruptive gas escape prior to 1980 as documented by Pb-210 deficits is on the order of a decade using the model developed by Gauthier and Condomines (Earth Planet. Sci. Lett. 172 (1999) 111-126) for a non-renewed magma chamber and efficient Rn removal. The time required to build-up Pb-210 excess is much shorter (months) as can be observed from steady increases of (Pb-210/Ra-226) with time during 1980-1982. The formation of Pb-210 excess requires both rapid gas transport through the magma and periodic blocking of gas escape routes. Superposed on this time trend is the natural variability of (Pb-210/Ra-226) in a single eruption caused by tapping magma from various depths. The two time scales of gas transport, to create both Pb-210 deficits and Pb-210 excesses, cannot be reconciled in a single event. Rather Pb-210 deficits are associated with pre-eruptive diffuse degassing, while Pb-210 excesses document the more vigorous degassing associated with eruption and recharge of the system. (c) 2006 Elsevier B.V. All rights reserved.

Pb-210-Ra-226 and Ra-228-Th-232 systematics in young arc lavas: implications for magma degassing and ascent rates

New data show that island arc rocks have (Pb-210/Ra-226)(o) ratios which range from as low as 0.24 up to 2.88. In contrast, (Ra-22S/Th-232) appears always within error of I suggesting that the large Ra-226-excesses observed in arc rocks were generated more than 30 years ago. This places a maximum estimate on melt ascent velocities of around 4000 m/year and provides further confidence that the Ra-226 excesses reflect deep (source) processes rather than shallow level alteration or seawater contamination. Conversely, partial melting must have occurred more than 30 years prior to eruption. The Pb-210 deficits are most readily explained by protracted magma degassing. Using published numerical models, the data suggest that degassing occurred continuously for periods up to several decades just prior to eruption but no link with eruption periodicity was found. Longer periods are required if degassing is discontinuous, less than 100% efficient or if magma is recharged or stored after degassing. The long durations suggest much of this degassing occurs at depth with implications for the formation of hydrothermal and copper-porphyry systems. A suite of lavas erupted in 1985-1986 from Sangeang Api volcano in the Sunda arc are characterised by deficits of Pb-210 relative to Ra-226 from which 6-8 years of continuous Rn-222 degassing would be inferred from recent numerical models. These data also form a linear (Pb-210)/Pb-(Ra-226)/Pb array which might be interpreted as a 71-year isochron. However, the array passes through the origin suggesting displacement downwards from the equiline in response to degassing and so the slope of the array is inferred not to have any age significance. Simple modelling shows that the range of (Ra-226)/Pb ratios requires thousands of years to develop consistent with differentiation occurring in response to cooling at the base of the crust. Thus, degassing post-dated, and was not responsible for magma differentiation. The formation, migration and extraction of gas bubbles must be extremely efficient in mafic magma whereas the higher viscosity of more siliceous magmas retards the process and can lead to Pb-210 excesses. A possible negative correlation between (Pb-210/Ra-226)(o) and SO2 emission rate requires further testing but may have implications for future eruptions. (C) 2004 Elsevier B.V. All rights reserved.

The magma budget of Volcan Arenal, Costa Rica from 1980 to 2004

Using topographic data collected by radar interferometry, stereo-photogrammetry, and field survey we have measured the changing surface of Volcan Arenal in Costa Rica over the period from 1980 to 2004. During this time this young volcano has mainly effused basaltic andesite lava, continuing the activity that began in 1968. Explosive products form only a few percent of the volumetric output. We have calculated digital elevation models for the years 1961, 1988 and 1997 and modified existing models for 2000 and 2004. From these we have estimated the volume of lava effused and coupled this with the data presented by an earlier study for 1968-1980. We find that a dense rock equivalent volume of 551 M m(3) was effused from 1968 to 2004. The dense rock equivalent effusion rate fell from about 2 m(3) s(-1) to about 0.1-0.2 m(3) s(-1) over the same period, with an average rate of about 0.5 m(3) s(-1). Between 1980 and 2004, the average effusion rate was 0.36 m(3) s(-1), a similar rate to that measured between 1974 and 1980. There have been two significant deviations from this long-term rate. The effusion rate increased from 1984 to 1991, at the same time as explosivity increased. After a period of moderate effusion rates in the 1990s, the rate fell to lower levels around 1999. (c) 2006 Elsevier B.V. All rights reserved.

Geochemical precursors to volcanic activity at Mount St. Helens, USA

The importance of the interplay between degassing and crystallization before and after the eruption of Mount St. Helens (Washington, USA) in 1980 is well established. Here, we show that degassing occurred over a period of decades to days before eruptions and that the manner of degassing, as deduced from geochemicai signatures within the magma, was characteristic of the eruptive style. Trace element (lithium) and short-lived radioactive isotope (lead-210 and radium-226) data show that ascending magma stalled within the conduit, leading to the accumulation of volatiles and the formation of lead-210 excesses, which signals the presence of degassing magma at depth.

Lava production at Soufriere Hills Volcano, Montserrat: 1995-2009

[1] We estimate that about 1 km3 of andesitic lava has been produced at Soufrière Hills Volcano, Montserrat from 1995 to 2009. There were three major episodes of extrusion, each lasting about 2 to 3.5 years and producing about 280 to 340 M m3 of lava, and one minor episode. Our estimates account for the dense rock equivalent volumetric contributions from the core and talus components of the lava dome, pyroclastic flow deposits and air-fall deposits. By 2005 at least two thirds of the erupted mass has already entered the sea. The average lava flux across the major extrusion episodes has been 3–5 m3s−1, with short-period (10–15 days) pulses up to 10–20 m3s−1. The first and third episodes of extrusion show similar flux histories suggesting similar behaviour of the system ten years apart. Waning flux towards the end of each episode may be caused by declining overpressure in the magma reservoir.

Database Engines for Geographical Information Systems

Our ability to identify, acquire, store, enquire on and analyse data is increasing as never before, especially in the GIS field. Technologies are becoming available to manage a wider variety of data and to make intelligent inferences on that data. The mainstream arrival of large-scale database engines is not far away. The experience of using the first such products tells us that they will radically change data management in the GIS field.

Ground deformation at Soufriere Hills Volcano, Montserrat during 1998-2000 measured by radar interferometry and GPS

We examine the motion of the ground surface on the Soufriere Hills Volcano, Montserrat between 1998 and 2000 using radar interferometry (InSAR). To minimise the effects of variable atmospheric water vapour on the InSAR measurements we use independently-derived measurements of the radar path delay from six continuous GPS receivers. The surfaces providing a measurable inter-ferometric signal are those on pyroclastic flow deposits, mainly emplaced in 1997. Three types of surface motion can be discriminated. Firstly, the surfaces of thick, valley-filling deposits subsided at rates of 150-120 mm/year in the year after emplacement to 50-30 mm/year two years later. This must be due to contraction and settling effects during cooling. The second type is the near-field motion localised within about one kilometre of the dome. Both subsidence and uplift events are seen and though the former could be due to surface gravitational effects, the latter may reflect shallow (< 1 km) pressurisation effects within the conduit/dome. Far-field motions of the surface away from the deeply buried valleys are interpreted as crustal strains. Because the flux of magma to the surface stopped from March 1998 to November 1999 and then resumed from November 1999 through 2000, we use InSAR data from these two periods to test the crustal strain behaviour of three models of magma supply: open, depleting and unbalanced. The InSAR observations of strain gradients of 75-80 mm/year/krn uplift during the period of quiescence on the western side of the volcano are consistent with an unbalanced model in which magma supply into a crustal magma chamber continues during quiescence, raising chamber pressure that is then released upon resumption of effusion. GPS motion vectors agree qualitatively with the InSAR displacements but are of smaller magnitude. The discrepancy may be due to inaccurate compensation for atmospheric delays in the InSAR data. (c) 2005 Elsevier B.V. All rights reserved.