Disappearance of a crater lake: implications for potential explosivity at Soufrière volcano,St Vincent,Lesser Antilles

Soufrière volcano in St Vincent, West Indies, is one of the most active volcanoes in the Eastern Caribbean with at least six eruptions since 1718 AD, the latest of which occurred in 1979. Prior to the 1979 eruption, the active crater hosted deep-water lakes during periods of repose, which were always replenished within a few years after the eruptions. In 1979, the crater was filled with 10⁸ m³ of fragmental material and, despite constant precipitation, has remained virtually dry ever since, with the exception of a small shallow pond. A resistivity survey was conducted in July 2006 to investigate groundwater occurrence in the crater. Results from the resistivity data inversion on several 2-D profiles show a shallow horizontal conductor across the crater floor, consistent with a water-saturated aquifer. They also show that the post-1979 pond, currently present in the crater lake is in fact an outcropping part of the groundwater water reservoir. The reservoir water table is ∼28 m above the pre-1979 lake level and reflects mass equilibrium in the system where constant seepage underground balances the meteoric recharge. We suggest that the groundwater body extends at depth to the bottom of the pre-1979 crater lake, either due to a significant structural discontinuity or because of a reduction of permeability at depth. The estimated maximum volume of water stored underground is 10–30 × 10⁶ m³ and energy considerations indicate that 2.4–7.3 × 10¹⁰ kg of magma would potentially be sufficient to vaporise the whole groundwater body. This amount of magma represents only 13–41% of the mass erupted during the last eruption in 1979 which was the smallest of the past 3 eruptions (1902, 1971–72, 1979). Since explosive phreatic or phreatomagmatic eruptions at Soufrière seem to be linked to magma-water interaction within confined space, the results from this survey suggests that phreatic or phreatomagmatic activity is a distinct possibility during future magma intrusion in the summit area, despite the apparent disappearance of water in the summit crater.     

Phenocryst-hosted melt inclusions record stalling of magma during ascent in the conduit and upper magma reservoir prior to vulcanian explosions, Soufrière Hills volcano, Montserrat, West Indies

The mechanics of explosive eruptions influence magma ascent pathways. Vulcanian explosions involve a stop–start mechanism that recurs on various timescales, evacuating the uppermost portions of the conduit. During the repose time between explosions, magma rises from depth and refills the conduit and stalls until the overpressure is sufficient to generate another explosion. We have analyzed major elements, Cl, S, H₂O, and CO₂ in plagioclase-hosted melt inclusions, sampled from pumice erupted during four vulcanian events at Soufrière Hills volcano, Montserrat, to determine melt compositions prior to eruption. Using Fourier transform infrared spectroscopy, we measured values up to 6.7 wt.% H₂O and 80 ppm CO₂. Of 42 melt inclusions, 81 % cluster between 2.8 and 5.4 wt.% H₂O (57 to 173 MPa or 2–7 km), suggesting lower conduit to upper magma reservoir conditions. We propose two models to explain the magmatic conditions prior to eruption. In Model 1, melt inclusions were trapped during crystal growth in magma that was stalled in the lower conduit to upper magma reservoir, and during trapping, the magma was undergoing closed-system degassing with up to 1 wt.% free vapor. This model can explain the melt inclusions with higher H₂O contents since these have sampled the upper parts of the magma reservoir. However, the model cannot explain the melt inclusions with lower H₂O because the timescale for plagioclase crystallization and melt inclusion entrapment is longer than the magma residence time in the conduit. In Model 2, melt inclusions were originally trapped at deeper levels of the magma chamber, but then lost hydrogen by diffusion through the plagioclase host during periodic stalling of the magma in the lower conduit system. In this second scenario, which we favor, the melt inclusions record re-equilibration depths within the lower conduit to upper magma reservoir.     

Seismicity of Ruapehu volcano,New Zealand, 1971–1996: a review

From 1971 until 1995, the style of seismicity at Ruapehu changed little, reflecting a period of relatively low eruptive activity and consequent long-term stability within the vent system. Volcanic earthquakes and volcanic tremor were both dominated by a frequency of about 2 Hz. Volcanic earthquakes accompanied all phreatic and phreatomagmatic eruptions, but not small hydrothermal eruptions that originated within Crater Lake. Furthermore, more than half of the M_L>3 volcanic earthquakes and changes in the reduced displacement of 2 Hz volcanic tremor by as much as a factor of 20 occurred without any accompanying eruptive activity. Three and 7 Hz volcanic tremor were also recorded, although never at lower-elevation seismometers. At times, this tremor was stronger at the summit seismometer than the 2 Hz tremor. Their source regions were independent of the 2 Hz source, and located at shallower depths. Volcano-tectonic earthquakes were generally unrelated to eruptive activity. The seismicity accompanying the 1995–1996 eruptive activity was significantly different from that of the period 1971 to 1995, and included volcanic tremor with a frequency of less than 1 Hz, simultaneous changes in the amplitude of the previously independent 2 Hz and 7 Hz volcanic tremor, and finally a change in the frequency content of volcanic earthquakes and volcanic tremor from 2 Hz to wideband. Path transmission effects play an important role in determining the characteristics of seismograms at Ruapehu. The presence of Crater Lake affects both the style of eruptions and the accompanying seismicity.     

Time-scales of recent Phlegrean Fields eruptions inferred from the application of a ‘diffusive fractionation’ model of trace elements

The variation of chemical element compositions in two pyroclastic sequences (Astroni 6 and Averno 2, Phlegrean Fields, Italy) is studied. Both sequences are compositionally zoned indicating a variability of melt compositions in the magma chamber prior to eruption. A clear dichotomy between the behaviour of major vs. trace elements is also observed in both sequences, with major elements displaying nearly linear inter-elemental trends and trace elements showing a variable scattered behaviour. Together with previous petrological investigations these observations are consistent with the hypothesis that magma mixing processes played a key role in the evolution of these two magmatic systems. Recently it has been suggested that mixing processes in igneous systems may strongly influence the mobility of trace elements inducing a ‘diffusive fractionation’ phenomenon, whose extent depends on the mixing time-scale. Here we merge information from 1) numerical simulations of magma mixing, and 2) magma mixing experiments (using as end-members natural compositions from Phlegrean Fields) to derive a relationship relating the degree of ‘diffusive fractionation’ to the mixing time-scales. Application of the ‘diffusive fractionation’ model to the two studied pyroclastic sequences allowed us to apply the relationship derived by numerical simulations and experiments to estimate the mixing time-scales for these two magmatic systems. Results indicate that mixing processes in Astroni 6 and Averno 2 systems lasted for approximately 2 and 9 days, respectively, prior to eruption.     

HCl emissions at Soufrière Hills Volcano,Montserrat, West Indies,during a second phase of dome building: November 1999 to October 2000

HCl:SO₂ mass ratios measured by open path Fourier transform spectroscopy (OP-FTIR) in the volcanic plume at Soufrière Hills Volcano, Montserrat, are presented for the second phase of dome building between November 1999 and November 2000. HCl:SO₂ mass ratios of greater than 1 and HCl emission rates of greater than 400 t day^–1 characterise periods of dome building for this volcano. The data suggest that chlorine partitions into a fluid phase as the magma decompresses and exsolves water during ascent. This is substantiated by a correlation between chlorine and water content in the melt (derived from the geochemical analysis of plagioclase melt inclusion and matrix glasses from phase I and II of dome growth). The matrix glass from the November 1999 and March 2000 domes indicate an open system degassing regime with a fluid-melt partition coefficient for chlorine of the order of 250–300. September 1997 glasses have higher chlorine contents and may indicate a switch to closed system degassing prior to explosive activity in September and October 1997. The OP-FTIR HCl time series suggests that HCl emission rate is strongly related to changes in eruption rate and we infer an emission rate of over 13.5 kt day^–1 HCl during a period of high extrusion rate in September 2000. A calculation of the HCl emission rate expected for varying extrusion rates from the open-system degassing model suggests a HCl emission rate of the order of 1–4 kt day^–1 is indicative of an extrusion rate of between 2 and 8 m³ s^–1. Monitoring of HCl at Soufrière Hills Volcano provide a proxy for extrusion rate, with changes in ratio between HCl and SO₂ occurring rapidly in the plume. Order of magnitude changes occur in HCl emission rates over the time-scale of hours to days, making these changes easy to detect during the day-to-day monitoring of the volcano. Mean water emission rates are calculated to range from 9–24 kt day^–1 during dome building activity, calculated from the predicted mass ratio of H₂O:HCl in the fluid at the surface and FTIR-derived HCl emission rates.     

The 1976–1982 Strombolian and phreatomagmatic eruptions of White Island,New Zealand: eruptive and depositional mechanisms at a ‘wet’ volcano

White Island is an active andesitic-dacitic composite volcano surrounded by sea, yet isolated from sea water by chemically sealed zones that confine a long-lived acidic hydrothermal system, within a thick sequence of fine-grained volcaniclastic sediment and ash. The rise of at least 10⁶ m³ of basic andesite magma to shallow levels and its interaction with the hydrothermal system resulted in the longest historical eruption sequence at White Island in 1976–1982. About 10⁷ m³ of mixed lithic and juvenile ejecta was erupted, accompanied by collapse to form two coalescing maar-like craters. Vent position within the craters changed 5 times during the eruption, but the vents were repeatedly re-established along a line linking pre-1976 vents. The eruption sequence consisted of seven alternating phases of phreatomagmatic and Strombolian volcanism. Strombolian eruptions were preceded and followed by mildly explosive degassing and production of incandescent, blocky juvenile ash from the margins of the magma body. Phreatomagmatic phases contained two styles of activity: (a) near-continuous emission of gas and ash and (b) discrete explosions followed by prolonged quiescence. The near-continuous activity reculted from streaming of magmatic volatiles and phreatic steam through open conduits, frittering juvennile shards from the margins of the magma and eroding loose lithic particles from the unconsolidated wall rock. The larger discrete explosions produced ballistic block aprons, downwind lobes of fall tephra, and cohesive wet surge deposits confined to the main crater. The key features of the larger explosions were their shallow focus, random occurrence and lack of precursors, and the thermal heterogeneity of the ejecta. This White Island eruption was unusual because of the low discharge rate of magma over an extended time period and because of the influence of a unique physical and hydrological setting. The low rate of magma rise led to very effective separation of magmatic volatiles and high fluxes of magmatic gas even during phreatic phases of the eruption. While true Strombolian phases did occur, more frequently the decoupled magmatic gas rose to interact with the conduit walls and hydrothermal system, producing phreatomagmatic eruptions. The form of these wet explosions was governed by a delicate balance between erosion and collapse of the weak conduit walls. If the walls were relatively stable, fine ash was slowly eroded and erupted in weak, near-continous phreatomagmatic events. When the walls were unstable, wall collapse triggered larger discrete phreatomagmatic explosions.     

Unusual types of degassing from melts of peripheral magma chambers of “dormant” El’brus volcano, Russia: Geochemical and mineralogical features

In addition to traditional degassing of the melt in the subsurface magma chambers of the “dormant” El’brus volcano, alsodegassing through pores and microcracks that occur in the top of magma chambers has also been detected. It is proven by studies of compactness, porosity, and permeability of the rocks. The speeds at which gases (H, He, H₂S, CO₂, F, and Cl) pass through gneiss and volcanic rocks were estimated. Magma chambers on the ground surface are expressed in stable thermal anomalies revealed by night-time thermal sounding from an NOAA satellite. The presence of magma chambers at depths of 2–12 km was proven by magnetotelluric sounding [Sobisevich et al., 2003] and gravity studies. In addition to occasional “columns” of bright-white fluorescence above the thermal anomalies, aerosol “clouds” and hydrogen flows were detected by lidar and hydrogen surveying [Alekseev et al., 2007, 2009]. Observation at the same sites detected steam outbursts occurring periodically, the snow-ice cover thaws and the smell of hydrogen sulfide is felt. Geochemical characteristics of degassing were studied by snow sampling from up to 1 m deep pits. They were taken within contours of the thermal anomalies, above active fault zones, in the sites of bright-white fluorescence “columns,” and on a new fumarole locality. It is shown that the degassing of melt was accompanied by the gas transporting many elements (Li, B, Si, P, S, Ca, Zn, Pb, Mo, Ba, W, Hg, Ag, U, Th, I, Au, and Pt) in a fine-grained state (a few microns or possibly nanometers) with an active participation of F and Cl. Native platinum, chalcopyrite, halite, sylvite, barite, gypsum, zircon, opal, chlorinated organics, etc. were for the first time discovered in the Mt. El’brus area using electron microscope studies of solid residue from dehydrated snow samples. “Hidden” ore mineralization genetically related to degassing of melts enriched in ore elements may be supposedly found in paleo- and present-day areas of volcanic activity.     

Lava penetrating water: the different behaviours of pāhoehoe and ‘a‘ā at the Nesjahraun, Tingvellir, Iceland

The Nesjahraun is a basaltic lava flow erupted from a subaerial fissure, extending NE along the Tingvellir graben from the Hengill central volcano that produced pāhoehoe lava followed by ‘a‘ā. The Nesjahraun entered Iceland’s largest lake, Tingvallavatn, along its southern shore during both phases of the eruption and exemplifies lava flowing into water in a lacustrine environment in the absence of powerful wave action. This study combines airborne light detection and ranging, sidescan sonar and Chirp seismic data with field observations to investigate the behaviour of the lava as it entered the water. Pāhoehoe sheet lava was formed during the early stages of the eruption. Along the shoreline, stacks of thin (5–20 cm thick), vesicular, flows rest upon and surround low (<5 m) piles of coarse, unconsolidated, variably oxidised spatter. Clefts within the lava run inland from the lake. These are 2–5 m wide, >2 m deep, ∼50 m long, spaced ∼50 m apart and have sub-horizontal striations on the walls. They likely represent channels or collapsed tubes along which lava was delivered into the water. A circular rootless cone, Eldborg, formed when water infiltrated a lava tube. Offshore from the pāhoehoe lavas, the gradient of the flow surface steepens, suggesting a change in flow regime and the development of a talus ramp. Later, the flow was focused into a channel of ‘a‘ā lava, ∼200–350 m wide. This split into individual flow lobes 20–50 m wide along the shore. ‘A‘ā clinker is exposed on the water’s edge, as well as glassy sand and gravel, which has been locally intruded by small (<1 m), irregularly shaped, lava bodies. The cores of the flow lobes contain coherent, but hackly fractured lava. Mounds consisting predominantly of scoria lapilli and the large paired half-cone of Grámelur were formed in phreatomagmatic explosions. The ‘a‘ā flow can be identified underwater over 1 km offshore, and the sidescan data suggest that the flow lobes remained coherent flowing down a gradient of <10°. The Nesjahraun demonstrates that, even in the absence of ocean waves, phreatomagmatic explosions are ubiquitous and that pāhoehoe flows are much more likely to break up on entering the water than ‘a‘ā flows, which, with a higher flux and shallow underlying surface gradient, can penetrate water and remain coherent over distances of at least 1 km.     

Seismic signals in geothermal areas of active volcanism: a case study from ‘La Fossa’, Vulcano (Italy)

Since the last eruption (1888–1890) volcanism at Vulcano, Aeolian Archipelago, southern Tyrrhenian Sea, has taken the form of persistent fumarolic activity. The gas-vapour phases of the geothermal systems are mainly discharged within two restricted areas about 1 km apart from each other, in the northern part of the island. These areas are La Fossa crater, and the beach fumaroles of the Baia di Levante. Fluids released at the two main fumarolic fields display quite different chemical and temperature characteristics, implying different origins. The local seismicity essentially takes the form of discrete shocks of shallow origin (depth1 km) at La Fossa, usually with energy < 10¹³ ergs. They are thought to be related to the uprise of pressurized hot gases and vapours discharged at the crater fumaroles. The present investigation points to the existence of two principal categories of seismic events (called M-shocks and N-shocks). These are short events (normally < 10 s). M-type shocks are thought to be due to resonance vibrations within the interior of the volcano, probably driven by the excitation of shock-waves within cavities deeply affected by deposition and alteration of self-sealant hydrothermal minerals. N-type events display features that resemble those of volcano-tectonic earthquakes, but have no recognizable S-phases. Here they are tentatively attributed to microfracturing of rocks which have been extensively hydrothermally altered. Results of the present study permit a preliminary conceptual model of the local shallow seismic processes in the framework of geochemical modelling of fumarolic activity and geological inferences from geothermal drilling.     

Formation of an ‘a’ā lava delta: insights from time-lapse multibeam bathymetry and direct observations during the Stromboli 2007 eruption

Five consecutive multibeam bathymetries collected before, during, and after the 2007 Stromboli eruption, combined with visual inspections, allowed us to document the morphological evolution of an ‘a’ā lava-fed delta and to reconstruct the main processes acting during its submarine emplacement. The 2007 Stromboli delta extended down to 600-m water depth and covered an area of 420?×?10³ m², with a maximum thickness of 65 m and a total estimated volume of ≈7?×?10⁶ m³, i.e., three times larger than its subaerial counterpart. The lava delta grew mainly through the emplacement of discrete lobes about 50–150 m in size. Lobes were fed from point sources along the paleoshoreline, and their subaqueous pathways seem to be mainly controlled by the submarine morphology, with flows mostly filling in depressions left by previous lobes. The main controlling factors on the lobe morphology and thickness are the effusion rates and the pre-eruption morphology, i.e., the geometry and gradients of the basal surface. Data also shows that sudden slope failure of portions of the submarine delta may occur simultaneously with accretion, implying that a significant part of the delta material can be transported to greater depths by submarine gravity flows. The present study is relevant for future monitoring and hazard assessment during the growth of active lava-fed deltas as well as for a better interpretation of ancient volcaniclastic successions inland.     

Testing the use of quartz ‘micro-hole’ photon-simulated luminescence for dating sediments from the central Lomonosov Ridge,Arctic Ocean

In the Arctic Ocean, direct dating methods are needed as an alternative to the radiocarbon (¹⁴C) method and to various indirect approaches for a longer stratigraphy. In past attempts to develop a luminescence sediment dating, the use of fine-silt (4–11 μm) mixture of quartz and feldspar grains from core tops has often produced large age overestimates by several ka. A recent application of micro-focused laser (‘micro-hole’) photon-stimulated luminescence (PSL) to medium-silt to fine-sand quartz grains (11–105 μm) from the core tops at the Alaska margin has been usefully accurate. To extend this approach to the central Arctic Ocean and to a larger grain size range, we applied micro-hole PSL dating to >11 μm quartz grains from core tops (0.5–2 cm horizon) from two sites on the central Lomonosov Ridge. We obtain a burial age estimate of ca. 2 ka for 11–62 μm grains at a multicore site 18 MC within a perched intra-ridge basin, in accord with ¹⁴C ages obtained on foraminifers. At nearby site 19 MC on the erosive ridge top, the micro-hole PSL dating of >90 μm quartz grains produces a burial age estimate of ∼ca. 25 ka, in accord with a foraminiferal ¹⁴C age of ca. 26 ka. However, the 11–90 μm grains from the same sample produce a much younger burial age estimate of ca. 9 ka. Thus, these two size fractions of quartz grains record different burial times and different deposition agents (icebergs vs. sea ice), providing insight into past sedimentary processes. Overall, our results confirm an earlier conclusion from micro-hole PSL dating study at the Alaska margin that medium to coarse silt fractions of quartz grains (11–90 μm or at least 62 μm) is the preferred material for direct dating of the last daylight exposure of detrital sediment in the Arctic Ocean.     

OSL dating of individual quartz ‘supergrains’ from the Ancient Middle Palaeolithic site of Cuesta de la Bajada,Spain

Dose saturation represents a fundamental limitation for obtaining finite optically stimulated luminescence (OSL) quartz ages over Middle Pleistocene timescales. Single-grain OSL studies typically reveal significant intra- and inter-sample variability in quartz dose saturation properties at the individual grain level. This enhanced variability may offer the potential to obtain extended-range chronologies exceeding the traditional upper age limits of multiple-grain OSL dating. However, there have been relatively few detailed assessments of single-grain OSL properties over high dose ranges. In this study we investigate extended-range single-grain OSL dating potential at Cuesta de la Bajada, one of the most important Ancient Middle Palaeolithic sites in the Iberian Peninsula. The quartz samples from this site exhibit exceptional dose saturation properties and contain significant populations of individual ‘supergrains’ with bright OSL signals, very high characteristic saturation dose (D₀) limits of 200 to >600 Gy, and dose-response curves that closely conform to single saturating exponential functions. Assessments of OSL signal composition and the ability to successfully recover a known laboratory dose of 470 Gy support the potential for obtaining high equivalent dose (D_e) values with reasonable fitting uncertainties. A series of novel sensitivity tests are used to assess potential biases in supergrain D_e estimation over high dose ranges related to (i) thermally transferred signal contributions, (ii) choice of dose-response curve fitting function, and (iii) insufficient dose saturation properties of individual grains. The single-grain D_e values obtained using the standard quality assurance criteria and novel sensitivity tests are consistent at 2σ, and are in agreement with quartz Ti-centre ages obtained on two of three paired electron spin resonance (ESR) samples collected from Cuesta de la Bajada. These comparisons support the suitability of our single-grain OSL results and suggest there may be good potential for using quartz supergrains to establish extended-range chronologies at some Middle Pleistocene sites. Comparisons with other single-grain OSL studies across the Iberian central plains suggest that favourable dose saturation properties may be influenced by regional-scale geological controls. The importance of undertaking single-grain OSL dating at Cuesta de la Bajada is also demonstrated by the results of synthetic aliquot experiments, which reveal multiple-grain age offsets of 60–170 ka when unsuitable grain types are inadvertently included in the final age calculations. The single-grain OSL results indicate that the main archaeological horizon at Cuesta de la Bajada (unit CB3) accumulated between 264 ± 22 and 293 ± 24 ka. These ages are consistent with the chronologies of other key early Middle Palaeolithic sites in the region (Ambrona AS6 and Gran Dolina TD10) and indicate that the Ancient Middle Palaeolithic techno-complex likely spread across Iberia from at least Marine Isotope Stage 9 onwards.     

Marine resources,biophysical processes,and environmental management of a tropical shelf seaway: Torres Strait,Australia–Introduction to the special issue

This special issue of Continental Shelf Research contains 20 papers giving research results produced as part of Australia's Torres Strait Co-operative Research Centre (CRC) Program, which was funded over a three-year period during 2003–2006. Marine biophysical, fisheries, socioeconomic-cultural and extension research in the Torres Strait region of northeastern Australia was carried out to meet three aims: 1) support the sustainable development of marine resources and minimize impacts of resource use in Torres Strait; 2) enhance the conservation of the marine environment and the social, cultural and economic well being of all stakeholders, particularly the Torres Strait peoples; and 3) contribute to effective policy formulation and management decision making. Subjects covered, including commercial and traditional fisheries management, impacts of anthropogenic sediment inputs on seagrass meadows and communication of science results to local communities, have broad applications to other similar environments.     

Modeling the Mogi magma source centre of the Santorini (Thera) volcano, Aegean Sea, Greece, 1994–1999, based on a numerical-topological approach

The aim of this study is the refinement of the dynamics of a recent (1994?C1999) minor, slow-inflation episode of the Santorini (Thera) volcano, famous for the Minoan (??3600 B.C.) eruption and the identification of the parameters of the magmatic source responsible for the inflation. Based on the Mogi source equations, on geodetic observations of base-line changes, on a topological, grid-search approach and on the reasonable assumption that the magma source remained practically stable in map view during the inflation period, we have been able to refine the location and depth (approximately 2.7 km) of the magma center. A tendency for increase of the magma pressure with time, roughly corresponding to a sphere with radius between 30 and 60 m, and a short deflation interval were also documented. The overall modeling was based on a topological method of inversion in two steps and for a selected 4-D grid. At a first step the system of Mogi-source equations was approximated by the intersection of the 4-D subspaces (defined by sets of grid points) each satisfying one observation equation on the basis of a grid-search procedure. At a second step, the best estimate of the Mogi source solution and its full variance-covariance matrix were defined using a common stochastic approach. The overall approach leads to a solution of a system of equations focusing on a 4-D space bounding significant minima in the misfits between model and observed values, and not on solutions focusing on single points, usually trapped in local minima. This study is important to understand a new phase of volcanic unrest since January 2011, while the proposed methodology, inspired from traditional navigation methods may be useful for other inversion problems leading to redundant systems of highly non-linear equations with n unknowns (i.e. topological solutions in the n-D space).