Modelling the lava dome extruded at Soufrière Hills Volcano, Montserrat, August 2005–May 2006. Part I: Dome shape and internal structure

Lava domes comprise core, carapace, and clastic talus components. They can grow endogenously by inflation of a core and/or exogenously with the extrusion of shear bounded lobes and whaleback lobes at the surface. Internal structure is paramount in determining the extent to which lava dome growth evolves stably, or conversely the propensity for collapse. The more core lava that exists within a dome, in both relative and absolute terms, the more explosive energy is available, both for large pyroclastic flows following collapse and in particular for lateral blast events following very rapid removal of lateral support to the dome. Knowledge of the location of the core lava within the dome is also relevant for hazard assessment purposes. A spreading toe, or lobe of core lava, over a talus substrate may be both relatively unstable and likely to accelerate to more violent activity during the early phases of a retrogressive collapse. Soufrière Hills Volcano, Montserrat has been erupting since 1995 and has produced numerous lava domes that have undergone repeated collapse events. We consider one continuous dome growth period, from August 2005 to May 2006 that resulted in a dome collapse event on 20th May 2006. The collapse event lasted 3 h, removing the whole dome plus dome remnants from a previous growth period in an unusually violent and rapid collapse event. We use an axisymmetrical computational Finite Element Method model for the growth and evolution of a lava dome. Our model comprises evolving core, carapace and talus components based on axisymmetrical endogenous dome growth, which permits us to model the interface between talus and core. Despite explicitly only modelling axisymmetrical endogenous dome growth our core–talus model simulates many of the observed growth characteristics of the 2005–2006 SHV lava dome well. Further, it is possible for our simulations to replicate large-scale exogenous characteristics when a considerable volume of talus has accumulated around the lower flanks of the dome. Model results suggest that dome core can override talus within a growing dome, potentially generating a region of significant weakness and a potential locus for collapse initiation.

Modelling the lava dome extruded at Soufrière Hills Volcano, Montserrat, August 2005–May 2006. Part II: Rockfall activity and talus deformation

During many lava dome-forming eruptions, persistent rockfalls and the concurrent development of a substantial talus apron around the foot of the dome are important aspects of the observed activity. An improved understanding of internal dome structure, including the shape and internal boundaries of the talus apron, is critical for determining when a lava dome is poised for a major collapse and how this collapse might ensue. We consider a period of lava dome growth at the Soufrière Hills Volcano, Montserrat, from August 2005 to May 2006, during which a 100 × 106 m3 lava dome developed that culminated in a major dome-collapse event on 20 May 2006. We use an axi-symmetrical Finite Element Method model to simulate the growth and evolution of the lava dome, including the development of the talus apron. We first test the generic behaviour of this continuum model, which has core lava and carapace/talus components. Our model describes the generation rate of talus, including its spatial and temporal variation, as well as its post-generation deformation, which is important for an improved understanding of the internal configuration and structure of the dome. We then use our model to simulate the 2005 to 2006 Soufrière Hills dome growth using measured dome volumes and extrusion rates to drive the model and generate the evolving configuration of the dome core and carapace/talus domains. The evolution of the model is compared with the observed rockfall seismicity using event counts and seismic energy parameters, which are used here as a measure of rockfall intensity and hence a first-order proxy for volumes. The range of model-derived volume increments of talus aggraded to the talus slope per recorded rockfall event, approximately 3 × 103–13 × 103 m3 per rockfall, is high with respect to estimates based on observed events. From this, it is inferred that some of the volumetric growth of the talus apron (perhaps up to 60–70%) might have occurred in the form of aseismic deformation of the talus, forced by an internal, laterally spreading core. Talus apron growth by this mechanism has not previously been identified, and this suggests that the core, hosting hot gas-rich lava, could have a greater lateral extent than previously considered.

Rapid topographic change measured by high-resolution satellite radar at Soufriere Hills Volcano, Montserrat, 2008-2010

High-resolution satellite radar observations of erupting volcanoes can yield valuable information on rapidly changing deposits and geomorphology. Using the TerraSAR-X (TSX) radar with a spatial resolution of about 2 m and a repeat interval of 11-days, we show how a variety of techniques were used to record some of the eruptive history of the Soufriere Hills Volcano, Montserrat between July 2008 and February 2010. After a 15-month pause in lava dome growth, a vulcanian explosion occurred on 28 July 2008 whose vent was hidden by dense cloud. We were able to show the civil authorities using TSX change difference images that this explosion had not disrupted the dome sufficient to warrant continued evacuation. Change difference images also proved to be valuable in mapping new pyroclastic flow deposits: the valley-occupying block-and-ash component tending to increase backscatter and the marginal surge deposits reducing it, with the pattern reversing after the event. By comparing east- and west-looking images acquired 12 hours apart, the deposition of some individual pyroclastic flows can be inferred from change differences. Some of the narrow upper sections of valleys draining the volcano received many tens of metres of rockfall and pyroclastic flow deposits over periods of a few weeks. By measuring the changing shadows cast by these valleys in TSX images the changing depth of infill by deposits could be estimated. In addition to using the amplitude data from the radar images we also used their phase information within the InSAR technique to calculate the topography during a period of no surface activity. This enabled areas of transient topography, crucial for directing future flows, to be captured.

Stress field control of eruption dynamics at a rift volcano: Nyamuragira, D.R.Congo

Using the record of 30 flank eruptions over the last 110 years at Nyamuragira, we have tested the relationship between the eruption dynamics and the local stress field. There are two groups of eruptions based on their duration (< 80days >) that are also clustered in space and time. We find that the eruptions fed by dykes parallel to the East African Rift Valley have longer durations (and larger volumes) than those eruptions fed by dykes with other orientations. This is compatible with a model for compressible magma transported through an elastic-walled dyke in a differential stress field from an over-pressured reservoir (Woods et al., 2006). The observed pattern of eruptive fissures is consistent with a local stress field modified by a northwest-trending, right lateral slip fault that is part of the northern transfer zone of the Kivu Basin rift segment. We have also re-tested with new data the stochastic eruption models for Nyamuragira of Burt et al. (1994). The time-predictable, pressure-threshold model remains the best fit and is consistent with the typically observed declining rate of sulphur dioxide emission during the first few days of eruption with lava emission from a depressurising, closed, crustal reservoir. The 2.4-fold increase in long-term eruption rate that occurred after 1977 is confirmed in the new analysis. Since that change, the record has been dominated by short-duration eruptions fed by dykes perpendicular to the Rift. We suggest that the intrusion of a major dyke during the 1977 volcano-tectonic event at neighbouring Nyiragongo volcano inhibited subsequent dyke formation on the southern flanks of Nyamuragira and this may also have resulted in more dykes reaching the surface elsewhere. Thus that sudden change in output was a result of a changed stress field that forced more of the deep magma supply to the surface. Another volcano-tectonic event in 2002 may also have changed the magma output rate at Nyamuragira.

Pulsatory andesite lava flow at Bagana Volcano

Using a time series of TerraSAR-X spaceborne radar images we have measured the pulsatory motion of an andesite lava flow over a 14-month period at Bagana volcano, Papua New Guinea. Between October 2010 and December 2011, lava flowed continuously down the western flank of the volcano forming a 3 km-long blocky lava flow with a channel, levees, overflows and branches. We captured four successive pulses of lava advancing down the channel system, the first such behaviour of an andesite flow to be recorded using radar. Each pulse had a volume of the order of 107 m3 emplaced over many weeks. The average extrusion rate estimated from the radar data was 0.92 ± 0.35 m3 s-1 , and varied between 0.3 and 1.8 m3 s-1, with higher rates occurring earlier in each pulse. This, together with observations of sulphur dioxide emissions, explosions and incandescence suggest a variable supply rate of magma through Bagana’s conduit as the most likely source of the pulsatory behaviour.

Co-eruptive and inter-eruptive surface deformation measured by satellite radar interferometry at Nyamuragira volcano, D.R. Congo, 1996 to 2010.

Interferometric Synthetic Aperture Radar (InSAR) measurements of surface deformation at Nyamuragira Volcano between 1996 and 2010 reveal a variety of co-eruptive and inter-eruptive signals. During 7 of the 8 eruptions in this period deformation was measured that is consistent with the emplacement of shallow near-vertical dykes feeding the eruptive fissures and associated with a NNW-trending fissure zone that traverses the summit caldera. Between eruptions the caldera and the summit part of this fissure zone subsided gradually (b3–5 cm/year). We also find evidence of post-eruption subsidence around the sites of the main vents of some flank eruptions (2002, 2004, 2006, and 2010). In the 6 months prior to the 2010 eruption a10-km wide zone centred on the caldera inflated by 1–2 cm. The low magnitude of this signal suggests that the presumed magma reservoir at 3–8 km depth contains highly compressible magma with little stored elastic strain energy. To the north of the caldera the fissure zone splits into WNW and NE branches around a zone that has a distinct InSAR signal. We interpret this zone to represent an elevated, 'stable' block of basement rocks buried by lavas within the Rift Zone.

Coupled subdaily and multiweek cycles during the lava dome erupt of Soufriere Hills Volcano, Montserrat

Observations of volcanoes extruding andesitic lava to produce lava domes often reveal cyclic behaviour. At Soufriere Hills Volcano, Montserrat, cycles with sub-daily and multi-week periods have been recognised on many occasions. These two types of cycle have been modelled separately as stick-slip magma flow at the junction between a dyke and an overlying cylindrical conduit (Costa et al. 2012), and as the filling and discharge of magma through the elastic-walled dyke (Costa et al., 2007a) respectively. Here, we couple these two models to simulate the behaviour over a period of well-observed multi-week cycles, with accompanying sub-daily cycles, from 13 May to 21 September 1997. The coupled model captures well the asymmetrical first-order behaviour: the first 40% of the multi-week cycle consists of high rates of lava extrusion during short period/high amplitude sub-daily cycles as the dyke reservoir discharges itself. The remainder of the cycle involves increasing pressurization as more magma is stored, and extrusion rate falls, followed by a gradual increase in the period of the sub-daily cycles.

Measuring large topographic change with InSAR: lava thicknesses, extrusion rate and subsidence rate at Santiaguito volcano, Guatemala

Lava flows can produce changes in topography on the order of 10s-100s of metres. A knowledge of the resulting volume change provides evidence about the dynamics of an eruption. We present a method to measure topographic changes from the differential InSAR phase delays caused by the height differences between the current topography and a Digital Elevation Model (DEM). This does not require a pre-event SAR image, so it does not rely on interferometric phase remaining coherent during eruption and emplacement. Synthetic tests predicts that we can estimate lava thickness of as little as �9 m, given a minimum of 5 interferograms with suitably large orbital baseine separations. In the case of continuous motion, such as lava flow subsidence, we invert interferometric phase simultaneously for topographic change and displacement. We demonstrate the method using data from Santiaguito volcano, Guatemala, and measure increases in lava thickness of up to 140 m between 2000 and 2009, largely associated with activity between 2000 and 2005. We find a mean extrusion rate of 0.43 +/- 0.06 m3/s, which lies within the error bounds of the longer term extrusion rate between 1922-2000. The thickest and youngest parts of the flow deposit were shown to be subsiding at an average rate of �-6 cm/yr. This is the first time that flow thickness and subsidence have been measured simultaneously. We expect this method to be suitable for measurment of landslides and other mass flow deposits as well as lava flows.

Observations of the eruption of the Sarychev volcano and simulations using the HadGEM2 climate model

In June 2009 the Sarychev volcano located in the Kuril Islands to the northeast of Japan erupted explosively, injecting ash and an estimated 1.2 ± 0.2 Tg of sulfur dioxide into the upper troposphere and lower stratosphere, making it arguably one of the 10 largest stratospheric injections in the last 50 years. During the period immediately after the eruption, we show that the sulfur dioxide (SO2) cloud was clearly detected by retrievals developed for the Infrared Atmospheric Sounding Interferometer (IASI) satellite instrument and that the resultant stratospheric sulfate aerosol was detected by the Optical Spectrograph and Infrared Imaging System (OSIRIS) limb sounder and CALIPSO lidar. Additional surface‐based instrumentation allows assessment of the impact of the eruption on the stratospheric aerosol optical depth. We use a nudged version of the HadGEM2 climate model to investigate how well this state‐of‐the‐science climate model can replicate the distributions of SO2 and sulfate aerosol. The model simulations and OSIRIS measurements suggest that in the Northern Hemisphere the stratospheric aerosol optical depth was enhanced by around a factor of 3 (0.01 at 550 nm), with resultant impacts upon the radiation budget. The simulations indicate that, in the Northern Hemisphere for July 2009, the magnitude of the mean radiative impact from the volcanic aerosols is more than 60% of the direct radiative forcing of all anthropogenic aerosols put together. While the cooling induced by the eruption will likely not be detectable in the observational record, the combination of modeling and measurements would provide an ideal framework for simulating future larger volcanic eruptions.

Sensitivity of OMI SO2 measurements to variable eruptive behaviour at Soufrière Hills Volcano, Montserrat

Since 2004, the satellite-borne Ozone Mapping Instrument (OMI) has observed sulphur dioxide (SO2) plumes during both quiescence and effusive eruptive activity at Soufrière Hills Volcano, Montserrat. On average, OMI detected a SO2 plume 4-6 times more frequently during effusive periods than during quiescence in the 2008-2010 period. The increased ability of OMI to detect SO2 during eruptive periods is mainly due to an increase in plume altitude rather than a higher SO2 emission rate. Three styles of eruptive activity cause thermal lofting of gases (Vulcanian explosions; pyroclastic flows; a hot lava dome) and the resultant plume altitudes are estimated from observations and models. Most lofting plumes from Soufrière Hills are derived from hot domes and pyroclastic flows. Although Vulcanian explosions produced the largest plumes, some produced only negligible SO2 signals detected by OMI. OMI is most valuable for monitoring purposes at this volcano during periods of lava dome growth and during explosive activity.

Fire History at the Southwestern Great Plains Margin, Capulin Volcano National Monument

This study documents historic fire events at Capulin Volcano National Monument over the last four centuries using dendrochronologically dated fire scars at two sites: the lower volcano lava flows (the Boca) and the adjacent canyon slopes (Morrow Ranch). The mean fire interval (MFI) was 12 years at the Boca site (before 1890) and 5.4 years (1600-1750) and 19.1 years (1751-1890) at the Morrow Ranch site. Data from the Boca and Morrow Ranch sites combined with the extremely pyrogenic landscape position of the volcano slopes indicate that the volcano slopes likely burned more frequently (e.g., MFI <5 yr). Around 1750, the fire regime appeared to transition to longer fire intervals, greater temporal synchrony among fire-scarred trees, and a higher proportion of trees scarred in fire years. Temporal variability in the fire regime at Capulin Volcano may reflect changes in human populations, climate, and land use.

Yellowstone Super-Volcano: Evalutaion, Potential Threats, and Possible effects on Nebraska Citizens Health and Prosperity

Abstract Yellowstone National Park is located over a hot spot under the North American tectonic plate and holds a potentially explosive super-volcano that has the ability to cause deadly consequences on the North American continent. After an eruption the surrounding region would see the greatest devastation, covered by pyroclastic deposits and thick ash fall exterminating most all life and destroying all structures in its path. In landscapes of greater distance from the event the consequences will be less dramatic yet still substantial. Records of previous eruption data from the Yellowstone super-volcano show that the ash fall out from the eruption can cover areas as large as one million square kilometers and could leave Nebraska covered in ash up to 10 centimeters thick. This would cause destruction of agriculture, extensive damage to structures, decreased temperatures, and potential respiratory hazards. The effects of volcanic ash on the human respiratory system have been shown to cause acute symptoms from heavy exposure. Symptoms include nasal irritation, throat irritation, coughing, and if preexisting conditions are present some can develop bronchial symptoms, which can last for a few days. People with bronchitis and asthma are shown to experience airway irritation and uncomfortable breathing. In most occurrences, exposure of volcanic ash is too short to cause long-term health hazards. Wearing facial protection can alleviate much of the symptoms. Most of the long-term ramifications of the eruption will be from the atmospheric changes caused from disruption of solar radiation, which will affect much of the global population. The most pertinent concerns for Nebraska citizens are from the accumulation of ash deposits over the landscape and the climatic perturbations. Potential mitigation procedures are essential to prepare our essentially unaware population of the threat that they may soon face if the volcano continues on its eruption cycle.

Ages (U-Pb SHRIMP and LA ICPMS) and stratigraphic evolution of the Neoproterozoic volcano-sedimentary successions from the extensional Camaqua Basin, Southern Brazil

During the Ediacaran, southern Brazil was the site of multiple episodes of volcanism and sedimentation, which are best preserved in the 3000 km(2) Camaqua Basin. The interlayered sedimentary and volcanic rocks record tectonic events and paleoenvironmental changes in a more than 10 km-thick succession. In this contribution, we report new U-Pb and Sm-Nd geochronological constraints for the 605 to 580 Ma Born Jardim Group, the 570 Ma Acampamento Velho Formation, and a newly-recognized 544 Ma volcanism. Depositional patterns of these units reveal the transition from a restricted, fault-bounded basin into a wide, shallow basin. The expansion of the basin and diminished subsidence rates are demonstrated by increasing areal distribution and compressed isopachs and increasing onlap of sediments onto the basement to the west. The Sm-Nd isotopic composition of the volcanic rocks indicates mixed sources, including crustal rocks from the adjacent basement. Both Neoproterozoic and Paleoproterozoic sources are indicated for the western part of the basin, whereas only the older Paleoproterozoic signature can be discerned in the eastern part of the basin. (C) 2011 International Association for Gondwana Research. Published by Elsevier B.V. All rights reserved.

The submarine volcano eruption off El Hierro Island: Effects on the scattering migrant biota and the evolution of the pelagic communities

[EN] The submarine volcano eruption off El Hierro Island (Canary Islands) on 10 October 2011 promoted dramatic perturbation of the water column leading to changes in the distribution of pelagic fauna. To study the response of the scattering biota, we combined acoustic data with hydrographic profiles and concurrent sea surface turbidity indexes from satellite imagery. We also monitored changes in the plankton and nekton communities through the eruptive and post-eruptive phases. Decrease of oxygen, acidification, rising temperature and deposition of chemicals in shallow waters resulted in a reduction of epipelagic stocks and a disruption of diel vertical migration (nocturnal ascent) of mesopelagic organisms. Furthermore, decreased light levels at depth caused by extinction in the volcanic plume resulted in a significant shallowing of the deep acoustic scattering layer. Once the eruption ceased, the distribution and abundances of the pelagic biota returned to baseline levels. There was no evidence of a volcano-induced bloom in the plankton community.