Calculation of lava effusion rates from Landsat TM data

 We present a thermal model to calculate the total thermal flux for lava flowing in tubes, on the surface, or under shallow water. Once defined, we use the total thermal flux to estimate effusion rates for active flows at Kilauea, Hawaii, on two dates. Input parameters were derived from Landsat Thematic Mapper (TM), field and laboratory measurements. Using these parameters we obtain effusion rates of 1.76±0.57 and 0.78±0.27 m³ s^–1 on 23 July and 11 October 1991, respectively. These rates are corroborated by field measurements of 1.36±0.14 and 0.89±0.09 m³s^–1 for the same dates (Kauahikaua et al. 1996). Using weather satellite (AVHRR) data of lower spatial resolution, we obtain similar effusion rates for an additional 26 dates between the two TM-derived measurements. We assume that, although total effusion rates at the source declined over the period, the shut down of the ocean entry meant that effusion rates for the surface flows alone remained stable. Such synergetic use of remotely sensed data provides measurements that can (a) contribute to monitoring flow-field evolution, and (b) provide reliable numerical data for input into rheological and thermal models. We look forward to being able to produce estimates for effusion rates using data from high-spatial-resolution sensors in the earth observing system (EOS) era, such as Landsat 7, the hyperspectral imager, the advanced spaceborne thermal emission spectrometer, and the advanced land imager. Received: 25 July 1997 / Accepted: 26 February 1998     

Simulations of the 2004 lava flow at Etna volcano using the magflow cellular automata model

Since the mechanical properties of lava change over time, lava flows represent a challenge for physically based modeling. This change is ruled by a temperature field which needs to be modeled. MAGFLOW Cellular Automata (CA) model was developed for physically based simulations of lava flows in near real-time. We introduced an algorithm based on the Monte Carlo approach to solve the anisotropic problem. As transition rule of CA, a steady-state solution of Navier-Stokes equations was adopted in the case of isothermal laminar pressure-driven Bingham fluid. For the cooling mechanism, we consider only the radiative heat loss from the surface of the flow and the change of the temperature due to mixture of lavas between cells with different temperatures. The model was applied to reproduce a real lava flow that occurred during the 2004–2005 Etna eruption. The simulations were computed using three different empirical relationships between viscosity and temperature.     

New formulae for estimating lava flow volumes at Mt. Etna Volcano, Sicily

 A new data set of Etna lava flows erupted since 1868 has been compiled from eight topographic maps of the volcano published at intervals since then. Volumes of 59 flows or groups of flows were measured from topographic difference maps. Most of these volumes are likely to be considerably more accurate than those published previously. We cut the number of flow volumes down to 25 by selecting those examples for which the volume of an individual eruption could be derived with the highest accuracy. This refined data set was searched for high correlations between flow volume and more directly measurable parameters. Only two parameters showed a correlation coefficient of 70% or greater: planimetric flow area A (70%) and duration of the eruption D (79%). If only short duration (<18 days) flows were used, flow length cubed, L³, had a correlation coefficient of 98%. Using combinations of measured parameters, much more significant correlations with volume were found. Dh had a correlation coefficient of 90% (h is the hydrostatic head of magma above the vent), and  , 92% (where W is mean width and E is the degree of topographic enclosure), and a combination of the two , 97%. These latter formulae were used to derive volumes of all eruptions back to 1868 to compare with those from the complete data set. Values determined from the formulae were, on average, lower by 16% (Dh), 7% (, and 19% . Received: 30 November 1998 / Accepted: 20 June 1999     

Oblique photogrammetry with visible and thermal images of active lava flows

Digital images from hand-held cameras are increasingly being acquired for scientific purposes, particularly where non-contact measurement is required. However, they frequently consist of oblique views with significant camera-to-object depth variations and occlusions that complicate quantitative analyses. Here, we report the use of oblique photogrammetric techniques to determine ground-based thermal camera orientations (position and pointing direction), and to generate scene information for lava flows at Mount Etna, Sicily. Multiple images from a consumer grade digital SLR camera are used to construct a topographic model and reference associated ground-based thermal imagery. We present data collected during the 2004–2005 eruption and use the derived surface model to apply viewing distance corrections (to account for atmospheric attenuation) to the thermal images on a pixel-by-pixel basis. For viewing distances of ~100 to 400 m, the corrections result in systematic changes in emissive power of up to ±3% with respect to values calculated assuming a uniform average viewing distance across an image.     

Viscous Newtonian laminar flow in a rectangular channel: application to Etna lava flows

 We introduce a 3D model for near-vent channelized lava flows. We assume the lava to be an isothermal Newtonian liquid flowing in a rectangular channel down a constant slope. The flow velocity is calculated with an analytical steady-state solution of the Navier-Stokes equation. The surface velocity and the flow rate are calculated as functions of the flow thickness for different flow widths, and the results are compared with those of a 2D model. For typical Etna lava flow parameters, the influence of levees on the flow dynamics is significant when the flow width is less than 25 m. The model predicts the volume flow rate corresponding to the surface velocity, taking into account that both depend on flow thickness. The effusion rate is a critical parameter to evaluate lava flow hazard. We propose a model to calculate the effusion rate given the lava flow width, the topograhic slope, the lava density, the surface flow velocity, and either the lava viscosity or the flow thickness. Received: 20 January 1998 / Accepted: 8 January 1999     

The exceptional activity and growth of the Southeast Crater, Mount Etna (Italy), between 1996 and 2001

Between 1971 and 2001, the Southeast Crater was the most productive of the four summit craters of Mount Etna, with activity that can be compared, on a global scale, to the opening phases of the Pu‘u ‘Ō‘ō-Kūpaianaha eruption of Kīlauea volcano, Hawai‘i. The period of highest eruptive rate was between 1996 and 2001, when near-continuous activity occurred in five phases. These were characterized by a wide range of eruptive styles and intensities from quiet, non-explosive lava emission to brief, violent lava-fountaining episodes. Much of the cone growth occurred during these fountaining episodes, totaling 105 events. Many showed complex dynamics such as different eruptive styles at multiple vents, and resulted in the growth of minor edifices on the flanks of the Southeast Crater cone. Small pyroclastic flows were produced during some of the eruptive episodes, when oblique tephra jets showered the steep flanks of the cone with hot bombs and scoriae. Fluctuations in the eruptive style and eruption rates were controlled by a complex interplay between changes in the conduit geometry (including the growth of a shallow magma reservoir under the Southeast Crater), magma supply rates, and flank instability. During this period, volume calculations were made with the aid of GIS and image analysis of video footage obtained by a monitoring telecamera. Between 1996 and 2001, the bulk volume of the cone increased by ~36×10⁶ m³, giving a total (1971–2001) volume of ~72×10⁶ m³. At the same time, the cone gained ~105 m in height, reaching an elevation of about 3,300 m. The total DRE volume of the 1996–2001 products was ~90×10⁶m³. This mostly comprised lava flows (72×10⁶ m³) erupted at the summit and onto the flanks of the cone. These values indicate that the productivity of the Southeast Crater increased fourfold during 1996–2001 with respect to the previous 25 years, coinciding with a general increase in the eruptive output rates and eruption intensity at Etna. This phase of intense summit activity has been followed, since the summer of 2001, by a period of increased structural instability of the volcano, marked by a series of important flank eruptions.     

Kinematic significance of mingling-rolling structures in lava flows: a case study from Porri Volcano (Salina, Southern Tyrrhenian Sea)

 A basaltic andesite lava flow from Porri Volcano (Salina, Southern Tyrrhenian Sea) is composed of two different magmas. Magma A (51 vol.% of crystals) has a dacitic glass composition, and magma B (18 vol.% of crystals), a basaltic glass composition. Magma B is hosted in A and consists of sub-spherical enclaves and boudin-like, banding and rolling structures (RS). Four types of RS have been recognized: σ–type;δ–type; complex σ-δ–types and transitional structures between sub-spherical enclaves and rolling structures. An analysis of the RS has been performed in order to reconstruct the flow kinematics and the mechanism of flow emplacement. Rolling structures have been selected in three sites located at different distances from the vent. In all sites most RS show the same sense of shear. Kinematic analysis of RS allows the degree of flow non-coaxiality to be determined. The non-coaxiality is expressed by the kinematic vorticity number Wk, a measure of the ratio Sr between pure shear strain rate and simple shear strain rate. The values of Wk calculated from the measured shapes of microscopic RS increase with increasing distance from the vent, from approximately 0.5 to 0.9. Results of the structural analysis reveal that the RS formed during the early–intermediate stage of flow emplacement. They represent originally sub-spherical enclaves deformed at low shear strain. At higher strain, RS deformed to give boudin-like and stretched banding structures. Results of the kinematic analysis suggest that high viscosity lava flows are heterogeneous non-ideal shear flows in which the degree of non-coaxiality increases with the distance from the vent. In the vent area, deformation is intermediate between simple shear and pure shear. Farther from the vent, deformation approaches ideal simple shear. Lateral extension processes occur only in the near-vent zone, where they develop in response to the lateral push of magma extruded from the vent. Lateral shortening processes develop in the distal zone and record the gravity-driven movement of the lava. The lava flow advanced by two main mechanisms, lateral translation and rolling motion. Lateral translation equals rolling near the vent, while rolling motion prevailed in the distal zones. Received: 6 November 1997 / Accepted: 29 November 1997     

Pyroclastic density currents resulting from the interaction of basaltic magma with hydrothermally altered rock: an example from the 2006 summit eruptions of Mount Etna,Italy

After 16 months of quiescence, Mount Etna began to erupt again in mid-July 2006. The activity was concentrated at and around the Southeast Crater (SEC), one of the four craters on the summit of Etna, and eruptive activity continued intermittently for 5 months. During this period, numerous vents displayed a wide range of eruptive styles at different times. Virtually all explosive activities took place at vents at the summit of the SEC and on its flanks. Eruptive episodes, which lasted from 1 day to 2 weeks, became shorter and more violent with time. Volcanic activity at these vents was often accompanied by dramatic mass-wasting processes such as collapse of parts of the cone, highly unusual flowage processes involving both old rocks and fresh magmatic material, and magma–water interaction. The most dramatic events took place on 16 November, when numerous rockfalls and pyroclastic density currents (PDCs) were generated during the opening of a large fracture on the SE flank of the SEC cone. The largest PDCs were clearly triggered explosively, and there is evidence that much of the energy was generated during the interaction of intruding magma with wet rocks on the cone’s flanks. The most mobile PDCs traveled up to 1 km from their source. This previously unknown process on Etna may not be unique on this volcano and is likely to have taken place on other volcanoes. It represents a newly recognized hazard to those who visit and work in the vicinity of the summit of Etna.     

Geochemistry of gases and waters discharged by the mud volcanoes at Paternò, Mt. Etna (Italy)

 Approximately 20 km south of Mt. Etna craters, at the contact between volcanic and sedimentary formations, three mud volcanoes discharge CO₂-rich gases and Na–Cl brines. The compositions of gas and liquid phases indicate that they are fed by a hydrothermal system for which temperatures of 100–150  °C were estimated by means of both gas and solute geothermometry. The hydrothermal system may be associated with CO₂-rich groundwaters over a large area extending from the central part of Etna to the mud volcanoes. Numerous data on the He, CH₄, CO₂ composition of the gases of the three manifestations, sampled over the past 5 years, indicate clearly that variations are due to separation processes of a CO₂-rich gas phase from the liquid. The effects of these processes have to be taken into account in the interpretation of the monitoring data collected for the geochemical surveillance of Etna volcano. Received: 4 September 1995 / Accepted: 14 February 1996     

Towards an Automatic Monitoring System of Infrasonic Events at Mt. Etna: Strategies for Source Location and Modeling

Active volcanoes characterized by open conduit conditions generate sonic and infrasonic signals, whose investigation provides useful information for both monitoring purposes and studying the dynamics of explosive processes. In this work, we discuss the automatic procedures implemented for a real-time application to the data acquired by a permanent network of five infrasound stations running at Mt. Etna volcano. The infrasound signals at Mt. Etna consist in amplitude transients, called infrasound events. The adopted procedure uses a multi-algorithm approach for event detection, counting, characterization and location. It is designed for an efficient and accurate processing of infrasound records provided by single-site and array stations. Moreover, the source mechanism of these events can be investigated off-line or in near real-time by using three different models: (1) Strombolian bubble; (2) resonating conduit and (3) Helmholtz resonator. The infrasound waveforms allow us to choose the most suitable model, to get quantitative information about the source and to follow the time evolution of the source parameters.     

Complex effusive events at Kīlauea as documented by the GOES satellite and remote video cameras

 GOES provides thermal data for all of the Hawaiian volcanoes once every 15 min. We show how volcanic radiance time series produced from this data stream can be used as a simple measure of effusive activity. Two types of radiance trends in these time series can be used to monitor effusive activity: (a) Gradual variations in radiance reveal steady flow-field extension and tube development. (b) Discrete spikes correlate with short bursts of activity, such as lava fountaining or lava-lake overflows. We are confident that any effusive event covering more than 10,000 m² of ground in less than 60 min will be unambiguously detectable using this approach. We demonstrate this capability using GOES, video camera and ground-based observational data for the current eruption of Kīlauea volcano (Hawai'i). A GOES radiance time series was constructed from 3987 images between 19 June and 12 August 1997. This time series displayed 24 radiance spikes elevated more than two standard deviations above the mean; 19 of these are correlated with video-recorded short-burst effusive events. Less ambiguous events are interpreted, assessed and related to specific volcanic events by simultaneous use of permanently recording video camera data and ground-observer reports. The GOES radiance time series are automatically processed on data reception and made available in near-real-time, so such time series can contribute to three main monitoring functions: (a) automatically alerting major effusive events; (b) event confirmation and assessment; and (c) establishing effusive event chronology. Received: 12 January 1999 / Accepted: 13 July 1999     

The dynamics of a channel-fed lava flow on Pico Partido volcano, Lanzarote

A short length of channel on Pico Partido volcano, Lanzarote, provides us the opportunity to examine the dynamics of lava flowing in a channel that extends over a sudden break in slope. The 1–2-m-wide, 0.5–2-m-deep channel was built during the 1730–1736 eruptions on Lanzarote and exhibits a sinuous, well-formed channel over a steep (11° slope) 100-m-long proximal section. Over-flow units comprising smooth pahoehoe sheet flow, as well as evidence on the inner channel walls for multiple (at least 11) flow levels, attest to unsteady flow in the channel. In addition, superelevation is apparent at each of the six bends along the proximal channel section. Superelevation results from banking of the lava as it moves around the bend thus causing preferential construction of the outer bank. As a result, the channel profile at each bend is asymmetric with an outer bank that is higher than the inner bank. Analysis of superelevation indicates flow velocities of ~8 m s^–1. Our analysis of the superelevation features is based on an inertia-gravity balance, which we show is appropriate, even though the down-channel flow is in laminar flow. We use a viscosity-gravity balance model, together with the velocities calculated from superelevation, to obtain viscosities in the range 25–60 Pa s (assuming that the lava behaved as a Newtonian liquid). Estimated volume fluxes are in the range 7–12 m³ s^–1. An apparent down-flow increase in derived volume flux may have resulted from variable supply or bulking up of the flow due to vesiculation. Where the channel moves over a sharp break in slope and onto slopes of ~6°, the channel becomes less well defined and widens considerably. At the break of slope, an elongate ridge extends across the channel. We speculate that this ridge was formed as a result of a reduction in velocity immediately below the break of slope to allow deposition of entrained material or accretion of lava to the channel bed as a result of a change in flow regime or depth.     

CO2 output and δ13C(CO2) from Mount Etna as indicators of degassing of shallow asthenosphere

 An estimated average CO₂ output from Etna's summit craters in the range of 13±3 Mt/a has recently been determined from the measured SO₂ output and measured CO₂/SO₂ molar ratios. To this amount the CO₂ output emitted diffusely from the soil (≈ 1 Mt/a) and the amount of CO₂ dissolved in Etna's aquifers (≈ 0.25 Mt/a) must be added. Data on the solubility of CO₂ in Etnean magmas at high temperature and pressure allow the volume of magma involved in the release of such an amount of this gas to be estimated. This volume of magma (≈ 0.7 km³/a) is approximately 20 times greater than the volume of magma erupted annually during the period 1971–1995. On the basis of C-isotopic data of CO₂ collected in the Etna area and of new hypotheses on the source of Mediterranean magmas, significant contributions of CO₂ from non-magmatic sources to the total output from Etna are unlikely. Such large outputs of CO₂ and also of SO₂ from Etna could be due to an anomalously shallow asthenosphere beneath the volcano that allows a continuous escape of gases toward the surface, even without migration of magma. Received: 7 August 1996 / Accepted: 9 November 1996     

Observations of Mt. Etna volcanic ash plumes in 2006: An integrated approach from ground-based and polar satellite NOAA–AVHRR monitoring system

Mt. Etna, in Sicily (Italy), is one of the world's most frequent emitters of volcanic plumes. During the last ten years, Etna has produced copious tephra emission and fallout that have damaged the inhabited and cultivated areas on its slopes and created serious hazards to air traffic. Recurrent closures of the Catania International airport have often been necessary, causing great losses to the local economy. Recently, frequent episodes of ash emission, lasting from a few hours to days, occurred from July to December 2006, necessitating a look at additional monitoring techniques, such as remote sensing. The combination of a ground monitoring system with polar satellite data represents a novel approach to monitor Etna's eruptive activity, and makes Etna one of the few volcanoes for which this surveillance combination is routinely available.In this work, ash emission information derived from an integrated approach, based on comparing ground and NOAA–AVHRR polar satellite observations, is presented. This approach permits us to define the utility of real time satellite monitoring systems for both sporadic and continuous ash emissions. Using field data (visible observations, collection of tephra samples and accounts by local inhabitants), the duration and intensity of most of the tephra fallout events were evaluated in detail and, in some cases, the order of magnitude of the erupted volume was estimated. The ground data vs. satellite data comparison allowed us to define five different categories of Etna volcanic plumes according to their dimensions and plume height, taking into account wind intensity. Using frequent and good quality satellite data in real time, this classification scheme could prove helpful for investigations into a possible correlation between eruptive intensity and the presence and concentration of ash in the volcanic plume. The development and improvement of this approach may constitute a powerful warning system for Civil Protection, thus preventing unnecessary airport closures.     

Thick lava flows of Karisimbi Volcano, Rwanda: insights from SIR-C interferometric topography

 We use a digital elevation model (DEM) derived from interferometrically processed SIR-C radar data to estimate the thickness of massive trachyte lava flows on the east flank of Karisimbi Volcano, Rwanda. The flows are as long as 12 km and average 40–60 m (up to >140 m) in thickness. By calculating and subtracting a reference surface from the DEM, we derived a map of flow thickness, which we used to calculate the volume (up to 1 km³ for an individual flow, and 1.8 km³ for all the identified flows) and yield strength of several flows (23–124 kPa). Using the DEM we estimated apparent viscosity based on the spacing of large folds (1.2×10¹² to 5.5×10¹² Pa s for surface viscosity, and 7.5×10¹⁰ to 5.2×10¹¹ Pa s for interior viscosity, for a strain interval of 24 h). We use shaded-relief images of the DEM to map basic flow structures such as channels, shear zones, and surface folds, as well as flow boundaries. The flow thickness map also proves invaluable in mapping flows where flow boundaries are indistinct and poorly expressed in the radar backscatter and shaded-relief images. Received: 6 September 1997 / Accepted: 15 May 1998     

Thermal monitoring of Lascar Volcano, Chile, using infrared data from the along-track scanning radiometer: a 1992–1995 time series

 Lascar Volcano (22°22'S, 67°44'W) is the most active volcano of the central Andes of northern Chile. Activity since 1984 has been characterised by periods of lava dome growth and decay within the active crater, punctuated by explosive eruptions. We present herein a technique for monitoring the high-temperature activity within the active crater using frequent measurements of emitted shortwave infrared (SWIR) radiation made by the spaceborne along-track scanning radiometer (ATSR). The ATSR is an instrument of low spatial resolution (pixels 1 km across) that shares certain characteristics with the MODIS instrument, planned for use as a volcano monitoring tool in the NASA EOS Volcanology Project. We present a comprehensive time series of over 60 cloud- and plume-free nighttime ATSR observations for 1992–1995, a period during which Lascar experienced its largest historical eruption. Variations in short wavelength infrared flux relate directly to changes in high-temperature surfaces within the active crater. From these data, interpretations can be made that supplement published field reports and that can document the presence and status of the lava dome during periods where direct, ground-based, observations are lacking. Our data agree with less frequent information collected from sensors with high spatial resolution, such as the Landsat thematic mapper (Oppenheimer et al. 1993) and are consistent with field observations and models that relate subsidence of the dome to subsequent explosive eruptions (Matthews et al., 1997). Most obviously, Lascar's major April 1993 eruption follows a period in which the magnitude of emitted shortwave infrared radiation fell by 90%. At this time subsidence of the 1991–1992 lava dome was reported by field observers and this subsidence is believed to have impeded the escape of hot volatiles and ultimately triggered the eruption (Smithsonian Institution 1993a). Extrapolating beyond the period for which field observations of the summit are available, our data show that the vulcanian eruption of 20 July 1995 occurred after a period of gradual increase in short wavelength infrared flux throughout 1994 and a more rapid flux decline during 1995. We attribute this additional, otherwise undocumented, cycle of increasing and decreasing SWIR radiance as most likely representing variations in degassing through fumaroles contained within the summit crater. Alternatively, it may reflect a cycle of dome growth and decay. The explosive eruption of 17 December 1993 appears to have followed a similar, but shorter, variation in SWIR flux, and we conclude that large explosive eruptions are more likely when the 1.6-μm signal has fallen from a high to a low level. The ATSR instrument offers low-cost data at high temporal resolution. Despite the low spatial detail of the measurements, ATSR-type instruments can provide data that relate directly to the status of Lascar's lava dome and other high-temperature surfaces. We suggest that such data can therefore assist with predictions of eruptive behaviour, deduced from application of physical models of lava dome development at this and similar volcanoes. Received: 1 October 1996 / Accepted: 13 January 1997     

The changing morphology of an open lava channel on Mt. Etna

An open channel lava flow on Mt. Etna (Sicily) was observed during May 30–31, 2001. Data collected using a forward looking infrared (FLIR) thermal camera and a Minolta-Land Cyclops 300 thermal infrared thermometer showed that the bulk volume flux of lava flowing in the channel varied greatly over time. Cyclic changes in the channel's volumetric flow rate occurred over several hours, with cycle durations of 113–190 min, and discharges peaking at 0.7 m³ s⁻¹ and waning to 0.1 m³ s⁻¹. Each cycle was characterized by a relatively short, high-volume flux phase during which a pulse of lava, with a well-defined flow front, would propagate down-channel, followed by a period of waning flow during which volume flux lowered. Pulses involved lava moving at relatively high velocities (up to 0.29 m s⁻¹) and were related to some change in the flow conditions occurring up-channel, possibly at the vent. They implied either a change in the dense rock effusion rate at the source vent and/or cyclic-variation in the vesicle content of the lava changing its bulk volume flux. Pulses would generally overspill the channel to emplace pāhoehoe overflows. During periods of waning flow, velocities fell to 0.05 m s^–1. Blockages forming during such phases caused lava to back up. Occasionally backup resulted in overflows of slow moving ‘a‘ā that would advance a few tens of meters down the levee flank. Compound levees were thus a symptom of unsteady flow, where overflow levees were emplaced as relatively fast moving pāhoehoe sheets during pulses, and as slow-moving ‘a‘ā units during backup. Small, localized fluctuations in channel volume flux also occurred on timescales of minutes. Volumes of lava backed up behind blockages that formed at constrictions in the channel. Blockage collapse and/or enhanced flow under/around the blockage would then feed short-lived, wave-like, down-channel surges. Real fluctuations in channel volume flux, due to pulses and surges, can lead to significant errors in effusion rate calculations. Editorial responsibility: A. Woods     

Voluminous lava flows at Oldoinyo Lengai in 2006: chronology of events and insights into the shallow magmatic system

The largest natrocarbonatite lava flow eruption ever documented at Oldoinyo Lengai, NW Tanzania, occurred from March 25 to April 5, 2006, in two main phases. It was associated with hornito collapse, rapid extrusion of lava covering a third of the crater and emplacement of a 3-km long compound rubbly pahoehoe to blocky aa-like flow on the W flank. The eruption was followed by rapid enlargement of a pit crater. The erupted natrocarbonatite lava has high silica content (3% SiO₂). The eruption chronology is reconstructed from eyewitness and news media reports and Moderate Resolution Imaging Spectroradiometer (MODIS) satellite data, which provide the most reliable evidence to constrain the eruption’s onset and variations in activity. The eruption products were mapped in the field and the total erupted lava volume estimated at 9.2 ± 3.0 × 10⁵ m³. The event chronology and field evidence are consistent with vent construct instability causing magma mixing and rapid extrusion from shallow reservoirs. It provides new insights into and highlights the evolution of the shallow magmatic system at this unique natrocarbonatite volcano.     

Relationships between deformation and mixing processes in lava flows: a case study from Salina (Aeolian Islands,Tyrrhenian Sea)

 The massive unit of a lava flow from Porri volcano (Salina, Aeolian Islands) displays many unusual structures related to the physical interaction between two different magmas. The magma A represents approximately 80% of the exposed lava surface; it has a crystal content of 51 vol.% and a dacitic glass composition (SiO₂=63–64 wt.%). The magma B has a basaltic-andesite glass composition (SiO₂=54–55 wt.%) and a crystal content of approximately 18 vol.%. It occurs as pillow-like enclaves, banding, boudin-like and rolling structures which are hosted in magma A. Structural analysis suggests that banding and boudin-like structures are the result of the deformation of enclaves at different shear strain. The linear correlation between strain and stratigraphic height of the measured elements indicates a single mode of deformation. We deduce that the component B deformed according to a simple shear model. Glass analyses of the A–B boundary indicate that A and B liquids mix together at high shear strain, whereas only mingling occurs at low shear strain. This suggests that the amount of deformation (i.e. forced convection) plays an important role in the formation of hybrid magmas. High shear strain may induce stretching, shearing and rolling of fluids which promote both forced convection and dynamical diffusion processes. These processes allow mixing of magmas with large differences in their physical properties. Received: 15 July 1995 / Accepted: 30 May 1996     

The character of long-term eruptions: inferences from episodes 50–53 of the Pu'u 'Ō'ō-Kūpaianaha eruption of Kīlauea Volcano

 The Pu'u 'Ō'ō-Kūpaianaha eruption on the east rift zone of Kīlauea began in January 1983. The first 9 years of the eruption were divided between the Pu'u 'Ō'ō (1983–1986) and Kūpaianaha (1986–1992) vents, each characterized by regular, predictable patterns of activity that endured for years. In 1990 a series of pauses in the activity disturbed the equilibrium of the eruption, and in 1991, the output from Kūpaianaha steadily declined and a short-lived fissure eruption broke out between Kūpaianaha and Pu'u 'Ō'ō. In February 1992 the Kūpaianaha vent died, and, 10 days later, eruptive episode 50 began as a fissure opened on the uprift flank of the Pu'u 'Ō'ō cone. For the next year, the eruption was marked by instability as more vents opened on the flank of the cone and the activity was repeatedly interrupted by brief pauses in magma supply to the vents. Episodes 50–53 constructed a lava shield 60 m high and 1.3 km in diameter against the steep slope of the Pu'u 'Ō'ō cone. By 1993 the shield was pockmarked by collapse pits as vents and lava tubes downcut as much as 29 m through the thick deposit of scoria and spatter that veneered the cone. As the vents progressively lowered, the level of the Pu'u 'Ō'ō pond also dropped, demonstrating the hydraulic connection between the two. The downcutting helped to undermine the prominent Pu'u 'Ō'ō cone, which has diminished in size both by collapse, as a large pit crater formed over the conduit, and by burial of its flanks. Intervals of eruptive instability, such as that of 1991–1993, accelerate lateral expansion of the subaerial flow field both by producing widely spaced vents and by promoting surface flow activity as lava tubes collapse and become blocked during pauses. Received: 1 July 1997 / Accepted: 23 October 1997