Doppler radar sounding of volcanic eruption dynamics at Mount Etna

Besides their common use in atmospheric studies, Doppler radars are promising tools for the active remote sensing of volcanic eruptions but were little applied to this field. We present the observations made with a mid-power UHF Doppler radar (Voldorad) during a 7-h Strombolian eruption at the SE crater of Mount Etna on 11–12 October 1998. Main characteristics of radar echoes are retrieved from analysis of Doppler spectra recorded in the two range gates on either side of the jet axis. From the geometry of the sounding, the contribution of uprising and falling ejecta to each Doppler spectrum can be discriminated. The temporal evolution of total power backscattered by uprising targets is quite similar to the temporal evolution of the volcanic tremor and closely reproduces the overall evolution of the eruption before, during and after its paroxysm. Moreover, during the sharp decrease of eruptive activity following the paroxysm, detailed analysis of video (from camera recording), radar and seismic measurements reveals that radar and video signals start to decrease simultaneously, approximately 2.5 min after the tremor decline. This delay is interpreted as the ascent time through a magma conduit of large gas slugs from a shallow source roughly estimated at about 500 m beneath the SE crater. Detailed analysis of eruptive processes has been also made with Voldorad operating in a high sampling rate mode. Signature of individual outburst is clearly identified on the half part of Doppler spectra corresponding to rising ejecta: temporal variations of the backscattered power exhibit quasi periodic undulations, whereas the maximum velocity measured on each spectrum displays a sharp peak at the onset of each outburst followed by a slow decay with time. Periodicity of power variations (between 3.8 and 5.5 s) is in agreement with the occurrence of explosions visually observed at the SE vent. Maximum vertical velocities of over 160 m s^–1 were measured during the paraoxysmal stage and the renewed activity. Finally, by using a simplified model simulating the radar echoes characteristics, we show that when Voldorad is operating in high sampling rate mode, the power and maximum velocity variations are directly related to the difference in size and velocity of particles crossing the antenna beam.Editorial responsibility: A. Woods     

Eruption dynamics of the SW crater of Stromboli volcano,Italy — An interdisciplinary approach

We report on the evaluation of in situ measurements of eruption velocities and relative mass flux, collected by Doppler radar, as well as acoustic and infrared data, that were recorded at Stromboli volcano. Doppler radar observations were made alongside thermal, acoustic and seismic measurements, to (1) further investigate the complex waveforms of the SW crater by combining infrasonic, infrared and Doppler radar measurements, (2) establish a relationship between infrared, acoustic and Doppler radar measurements and (3) verify that all instruments observe the same behavior in terms of relative mass eruption rate. We also explore the relationship between kinetic and acoustic energy released during an eruption. Comparing the different methods to each other we are able to show that the mass erupted can be estimated either from the total reflected Doppler radar energy or from infrared observations. However, neither thermal nor reflected energy can provide a value in terms of absolute mass in kg. The erupted masses of different eruptions can only be evaluated relative to each other. Using the combined three data sets we especially focus on the eruption dynamics of the SW crater of Stromboli, namely its fluctuations in eruption strength. These pulses in one eruptive event, so called ’pulsations’ dominate more than 40% of the eruptions of SW crater. Previous models that explain pulsations to be generated by multiple consecutively exploding bubbles are supported by our combined analysis of infrasound, thermal data, particle velocities, and the reflected energy.     

Strombolian explosive styles and source conditions: insights from thermal (FLIR) video

Forward Looking Infrared Radiometer (FLIR) cameras offer a unique view of explosive volcanism by providing an image of calibrated temperatures. In this study, 344 eruptive events at Stromboli volcano, Italy, were imaged in 2001–2004 with a FLIR camera operating at up to 30 Hz. The FLIR was effective at revealing both ash plumes and coarse ballistic scoria, and a wide range of eruption styles was recorded. Eruptions at Stromboli can generally be classified into two groups: Type 1 eruptions, which are dominated by coarse ballistic particles, and Type 2 eruptions, which consist of an optically-thick, ash-rich plume, with (Type 2a) or without (Type 2b) large numbers of ballistic particles. Furthermore, Type 2a plumes exhibited gas thrust velocities (>15 m s⁻¹) while Type 2b plumes were limited to buoyant velocities (<15 m s⁻¹) above the crater rim. A given vent would normally maintain a particular gross eruption style (Type 1 vs. 2) for days to weeks, indicating stability of the uppermost conduit on these timescales. Velocities at the crater rim had a range of 3–101 m s⁻¹, with an overall mean value of 24 m s⁻¹. Mean crater rim velocities by eruption style were: Type 1 = 34 m s⁻¹, Type 2a = 31 m s⁻¹, Type 2b = 7 m s⁻¹. Eruption durations had a range of 6–41 s, with a mean of 15 s, similar among eruption styles. The ash in Type 2 eruptions originates from either backfilled material (crater wall slumping or ejecta rollback) or rheological changes in the uppermost magma column. Type 2a and 2b behaviors are shown to be a function of the overpressure of the bursting slug. In general, our imaging data support a broadening of the current paradigm for strombolian behavior, incorporating an uppermost conduit that can be more variable than is commonly considered.     

Vent conditions for expected eruptions at Vesuvius

Determining consistent sets of vent conditions for next expected eruptions at Vesuvius is crucial for the simulation of the sub-aerial processes originating the volcanic hazard and the eruption impact. Here we refer to the expected eruptive scales and conditions defined in the frame of the EC Exploris project, and simulate the dynamics of magma ascent along the volcanic conduit for sub-steady phases of next eruptions characterized by intensities of the Violent Strombolian (VS), Sub-Plinian 2 (SP2), and Sub-Plinian 1 (SP1) scale. Sets of conditions for the simulations are determined on the basis of the bulk of knowledge on the past history of Vesuvius [Cioni, R., Bertagnini, A., Santacroce, R., Andronico, D., Explosive activity and eruption scenarios at Somma–Vesuvius (Italy): towards a new classification scheme. Journal of Volcanology and Geothermal Research, this issue.]. Volatile contents (H₂O and CO₂) are parameterized in order to account for the uncertainty in their expected amounts for a next eruption. In all cases the flow in the conduit is found to be choked, with velocities at the conduit exit or vent corresponding to the sonic velocity in the two-phase non-equilibrium magmatic mixture. Conduit diameters and vent mixture densities are found to display minimum overlapping between the different eruptive scales, while exit gas and particle velocities, as well as vent pressures, largely overlap. Vent diameters vary from as low as about 5 m for VS eruptions, to 35–55 m for the most violent SP1 eruption scale. Vent pressures can be as low as less than 1 MPa for the lowest volatile content employed of 2 wt.% H₂O and no CO₂, to 7–8 MPa for highest volatile contents of 5 wt.% H₂O and 2 wt.% CO₂ and large eruptive scales. Gas and particle velocities at the vent range from 100–250 m/s, with a tendency to decrease, and to increase the mechanical decoupling between the phases, with increasing eruptive scale. Except for velocities, all relevant vent quantities are more sensitive to the volatile content of the discharged magma for the highest eruptive scales considered.     

The 1996 eruption at Gjálp,Vatnajökull ice cap,Iceland: efficiency of heat transfer,ice deformation and subglacial water pressure

The 13-day-long Gjálp eruption within the Vatnajökull ice cap in October 1996 provided important data on ice–volcano interaction in a thick temperate glacier. The eruption produced 0.8 km³ of mainly volcanic glass with a basaltic icelandite composition (equivalent to 0.45 km³ of magma). Ice thickness above the 6-km-long volcanic fissure was initially 550–750 m. The eruption was mainly subglacial forming a 150–500 m high ridge; only 2–4% of the volcanic material was erupted subaerially. Monitoring of the formation of ice cauldrons above the vents provided data on ice melting, heat flux and indirectly on eruption rate. The heat flux was 5–6×10⁵ W m^-2 in the first 4 days. This high heat flux can only be explained by fragmentation of magma into volcanic glass. The pattern of ice melting during and after the eruption indicates that the efficiency of instantaneous heat exchange between magma and ice at the eruption site was 50–60%. If this is characteristic for magma fragmentation in subglacial eruptions, volcanic material and meltwater will in most cases take up more space than the ice melted in the eruption. Water accumulation would therefore cause buildup of basal water pressure and lead to rapid release of the meltwater. Continuous drainage of meltwater is therefore the most likely scenario in subglacial eruptions under temperate glaciers. Deformation and fracturing of ice played a significant role in the eruption and modified the subglacial water pressure. It is found that water pressure at a vent under a subsiding cauldron is substantially less than it would be during static loading by the overlying ice, since the load is partly compensated for by shear forces in the rapidly deforming ice. In addition to intensive crevassing due to subsidence at Gjálp, a long and straight crevasse formed over the southernmost part of the volcanic fissure on the first day of the eruption. It is suggested that the feeder dyke may have overshot the bedrock–ice interface, caused high deformation rates and fractured the ice up to the surface. The crevasse later modified the flow of meltwater, explaining surface flow of water past the highest part of the edifice. The dominance of magma fragmentation in the Gjálp eruption suggests that initial ice thickness greater than 600–700 m is required if effusive eruption of pillow lava is to be the main style of activity, at least in similar eruptions of high initial magma discharge.Editorial responsibility: J. Donnelly-Nolan     

The Taurewa Eruptive Episode: evidence for climactic eruptions at Ruapehu volcano, New Zealand

 The ca. 10,500 years B.P. eruptions at Ruapehu volcano deposited 0.2–0.3 km³ of tephra on the flanks of Ruapehu and the surrounding ring plain and generated the only known pyroclastic flows from this volcano in the late Quaternary. Evidence of the eruptions is recorded in the stratigraphy of the volcanic ring plain and cone, where pyroclastic flow deposits and several lithologically similar tephra deposits are identified. These deposits are grouped into the newly defined Taurewa Formation and two members, Okupata Member (tephra-fall deposits) and Pourahu Member (pyroclastic flow deposits). These eruptions identify a brief (<ca. 2000-year) but explosive period of volcanism at Ruapehu, which we define as the Taurewa Eruptive Episode. This Episode represents the largest event within Ruapehu's ca. 22,500-year eruptive history and also marks its culmination in activity ca. 10,000 years B.P. Following this episode, Ruapehu volcano entered a ca. 8000-year period of relative quiescence. We propose that the episode began with the eruption of small-volume pyroclastic flows triggered by a magma-mingling event. Flows from this event travelled down valleys east and west of Ruapehu onto the upper volcanic ring plain, where their distal remnants are preserved. The genesis of these deposits is inferred from the remanent magnetisation of pumice and lithic clasts. We envisage contemporaneous eruption and emplacement of distal pumice-rich tephras and proximal welded tuff deposits. The potential for generation of pyroclastic flows during plinian eruptions at Ruapehu has not been previously considered in hazard assessments at this volcano. Recognition of these events in the volcanological record is thus an important new factor in future risk assessments and mitigation of volcanic risk at Tongariro Volcanic Centre. Received: 5 July 1998 / Accepted: 12 March 1999     

Characterization of tephra deposits with limited exposure: the example of the two largest explosive eruptions at Nisyros volcano (Greece)

Explosive eruptions associated with tephra deposits that are only exposed in proximal areas are difficult to characterize. In fact, the determination of physical parameters such as column height, mass eruption rate, erupted volume, and eruption duration is mainly based on empirical models and is therefore very sensitive to the quality of the field data collected. We have applied and compared different modeling approaches for the characterization of the two main tephra deposits, the Lower Pumice (LP) and Upper Pumice (UP) of Nisyros volcano, Greece, which are exposed only within 5 km of the probable vent. Isopach and isopleth maps were compiled for two possible vent locations (on the north and on the south rim of the caldera), and different models were applied to calculate the column height, the erupted volume, and the mass eruption rate. We found a column height of about 15 km above sea level and a mass eruption rate of about 2 × 10⁷ kg/s for both eruptions regardless of the vent location considered. In contrast, the associated wind velocity for both UP and LP varied between 0 and 20 m/s for the north and south vent, respectively. The derived erupted volume for the south vent (considered as the best vent location) ranges between 2 and 27 × 10⁸ m³ for the LP and between 1 and 5 × 10⁸ m³ for the UP based on the application of four different methods (integration of exponential fit based on one isopach line, integration of exponential and power-law fit based on two isopach lines, and an inversion technique combined with an advection–diffusion model). The eruption that produced the UP could be classified as subplinian. Discrepancies associated with different vent locations are smaller than the discrepancies associated with the use of different models for the determination of erupted mass, plume height, and mass eruption rate. Proximal outcrops are predominantly coarse grained with ≥90 wt% of the clasts ranging between −6ϕ and 0ϕ. The associated total grainsize distribution is considered to result from a combination of turbulent fallout from both the plume margins and the umbrella region, and as a result, it is fines-depleted. Given that primary deposit thickness observed on Nisyros for both LP and UP is between 1 and 8 m, if an event of similar scale were to happen again, it would have a significant impact on the entire island with major damage to infrastructure, agriculture, and tourism. Neighboring islands and the continent could also be significantly affected.     

Gravity and aeromagnetic constraints on the structure of the Woods Mountains volcanic center, southeastern California

 The Woods Mountain volcanic center is a well-exposed, mildly alkaline volcanic center that formed during the Miocene in southeastern California. Detailed geologic mapping and geochemical studies have distinguished three major volcanic phases: precaldera, caldera forming, and postcaldera. Geologic mapping indicates that caldera formation occurred incrementally during eruptions of three large ignimbrites and continued into a period of voluminous intracaldera lava-flow eruptions. Rhyolitic ignimbrites and lava flows within the caldera are associated with large amplitude, circular gravity, and magnetic minima that are among the most prominent gravity and magnetic anomalies in southeastern California. Analysis of a Bouguer gravity anomaly map, reduced-to-the-pole magnetic intensity map, and three-dimensional gravity and magnetic models indicates that there is a single, funnel- to bowl-shaped caldera approximately 4 km thick and approximately 10 km wide at the surface. This model is consistent with other siliceous, pyroclastic-filled calderas on continental crust, except that most siliceous volcanic centers associated with more than one eruption are characterized by more than one caldera. Received: 20 December 1997 / Accepted: 15 October 1998     

GOES satellite and field observations of the 1998 eruption of Volcan Cerro Azul, Galápagos Islands

 The 1998 eruption of Volcan Cerro Azul, Isla Isabela, Galápagos Islands, was observed in near real-time by the Geostationary Operational Environmental Satellite-8 (GOES-8) weather satellite. Due to the remote location of the eruption site, 3.9-μm radiance values derived from GOES band 2 provide the best timing of the start and termination of the eruption, which occurred on 15 Sept. and 21 Oct., respectively. Throughout the 36-day long eruption, a total of 1335 thermal infrared images were collected, of which 851 were cloud-free and permitted the thermal anomaly to be detected. A detailed chronology including 77 separate events was assembled from the GOES data and field observations. Numerous attributes of the eruption were observed from the GOES data, including the sizes and dispersal of seven eruption plumes and the occurrence and timing of intra-caldera effusive activity. The growth of a lava flow on the SE flank, the formation of smoke and volcanic haze from the flank vent, and burning of vegetation caused by lava flows entering vegetated areas were monitored both on the ground and with the satellite data. In most cases GOES images were processed as they were received every 30 min and were then distributed over the Internet within minutes of reception. These data provided timely high-temporal information to field parties as well as enabled the documentation of the eruption. The GOES observations of Cerro Azul serve as a further example of the way in which the remote sensing community and field volcanologists can collaborate during future eruptions, and permit the temporal and spatial resolution requirements for future satellites systems to be better defined. Received: 27 April 1999 / Accepted: 21 February 2000     

Weather radar observations of the Hekla 2000 eruption cloud,Iceland

The Hekla eruption cloud on 26–27 February 2000 was the first volcanic cloud to be continuously and completely monitored advecting above Iceland, using the C-band weather radar near the Keflavík international airport. Real-time radar observations of the onset, advection, and waning of the eruption cloud were studied using time series of PPI (plan-position indicator) radar images, including VMI normal, Echotop, and Cappi level 2 displays. The reflectivity of the entire volcanic cloud ranges from 0 to >60 dBz. The eruption column above the vent is essentially characterised by VMI normal and Cappi level 2 values, >30 dBz, due to the dominant influence of lapilli and ash (tephra) on the overall reflected signal. The cloud generated by the column was advected downwind to the north-northeast. It is characterised by values between 0 and 30 dBz, and the persistence of these reflections likely result from continuing water condensation and freezing on ash particles. Echotop radar images of the eruption onset document a rapid ascent of the plume head with a mean velocity of ~30 to 50 m s^–1, before it reached an altitude of ~11–12 km. The evolution of the reflected cloud was studied from the area change in pixels of its highly reflected portions, >30 dBz, and tied to recorded volcanic tremor amplitudes. The synchronous initial variation of both radar and seismic signals documents the abrupt increase in tephra emission and magma discharge rate from 18:20 to 19:00 UTC on 26 February. From 19:00 the >45 dBz and 30–45 dBz portions of the reflected cloud decrease and disappear at about 7 and 10.5 h, respectively, after the eruption began, indicating the end of the decaying explosive phase. The advection and extent of the reflected eruption cloud were compared with eyewitness accounts of tephra fall onset and the measured mass of tephra deposited on the ground during the first 12 h. Differences in the deposit map and volcanic cloud radar map are due to the fact that the greater part of the deposit originates by fallout off the column margins and from the base of the cloud followed by advection of falling particle in lower level winds.Editorial responsibility: P. Mouginis-Mark     

Volcanic activity in the Azores. Report for 1959–1964

During the period 1959 to 1964 recorded volcanic activity in the Azores islands (Atlantic Ocean) was manifested by fumarolic fields in S. Miguel, Terceira, Graciosa, Pico and Faial and by two submarine eruptions off Pico and S. Jorge. Fumarolic activity increased in the Faial Capelinhos vent (in eruption 1957–58) after December 1963; it is concluded that the phenomenon is related to the 1963 and 1964 eruptions in the S. Jorge Channel. During February 1964 an important seismic swarm of more than 500 earthquakes, with volcanic tremor, occurred in S. Jorge.

     

New insights into Late Pleistocene explosive volcanic activity and caldera formation on Ischia (southern Italy)

A new pyroclastic stratigraphy is presented for the island of Ischia, Italy, for the period ∼75–50 ka BP. The data indicate that this period bore witness to the largest eruptions recorded on the island and that it was considerably more volcanically active than previously thought. Numerous vents were probably active during this period. The deposits of at least 10 explosive phonolite to basaltic-trachyandesite eruptions are described and interpreted. They record a diverse range of explosive volcanic activity including voluminous fountain-fed ignimbrite eruptions, fallout from sustained eruption columns, block-and-ash flows, and phreatomagmatic eruptions. Previously unknown eruptions have been recognised for the first time on the island. Several of the eruptions produced pyroclastic density currents that covered the whole island as well as the neighbouring island of Procida and parts of the mainland. The morphology of Ischia was significantly different to that seen today, with edifices to the south and west and a submerged depression in the centre. The largest volcanic event, the Monte Epomeo Green Tuff (MEGT) resulted in caldera collapse across all or part of the island. It is shown to comprise at least two thick intracaldera ignimbrite flow-units, separated by volcaniclastic sediments that were deposited during a pause in the eruption. Extracaldera deposits of the MEGT include a pumice fall deposit emplaced during the opening phases of the eruption, a widespread lithic lag breccia outcropping across much of Ischia and Procida, and a distal ignimbrite in south-west Campi Flegrei. During this period the style and magnitude of volcanism was dictated by the dynamics of a large differentiated magma chamber, which was partially destroyed during the MEGT eruption. This contrasts with the small-volume Holocene and historical effusive and explosive activity on Ischia, the timing and distribution of which has been controlled by the resurgence of the Monte Epomeo block. The new data contribute to a clearer understanding of the long-term volcanic and magmatic evolution of Ischia.     

Cyclic flank-vent and central-vent eruption patterns

 Many basaltic and andesitic polygenetic volcanoes have cyclic eruptive activity that alternates between a phase dominated by flank eruptions and a phase dominated by eruptions from a central vent. This paper proposes the use of time-series diagrams of eruption sites on each polygenetic volcano and intrusion distances of dikes to evaluate volcano growth, to qualitatively reconstruct the stress history within the volcano, and to predict the next eruption site. In these diagrams the position of an eruption site is represented by the distance from the center of the volcano and the clockwise azimuth from north. Time-series diagrams of Mauna Loa, Kilauea, Kliuchevskoi, Etna, Sakurajima, Fuji, Izu-Oshima, and Hekla volcanoes indicate that fissure eruption sites of these volcanoes migrated toward the center of the volcano linearly, radially, or spirally with damped oscillation, occasionally forming a hierarchy in convergence-related features. At Krafla, terminations of dikes also migrated toward the center of the volcano with time. Eruption sites of Piton de la Fournaise did not converge but oscillated around the center. After the convergence of eruption sites with time, the central eruption phase is started. The intrusion sequence of dikes is modeled, applying crack interaction theory. Variation in convergence patterns is governed by the regional stress and the magma supply. Under the condition that a balance between regional extension and magma supply is maintained, the central vent convergence time during the flank eruption phase is 1–10 years, whereas the flank vent recurrence time during the central eruption phase is greater than 100 years owing to an inferred decrease in magma supply. Under the condition that magma supply prevails over regional extension, the central vent convergence time increases, whereas the flank vent recurrence time decreases owing to inferred stress relaxation. Earthquakes of M≥6 near a volcano during the flank eruption phase extend the central vent convergence time. Earthquakes during the central eruption phase promote recurrence of flank eruptions. Asymmetric distribution of eruption sites around the flanks of a volcano can be caused by local stress sources such as an adjacent volcano. Received: 18 March 1996 / Accepted: 14 January 1997     

Development of unconfined historic lava flow fields in Tenerife: implications for the mitigation of risk from a future eruption

Historic and recent (last 2,000?years) eruptions on the active volcanic island of Tenerife have been predominantly effusive, indicating that this is the most probable type of activity to be expected in the near future. In the past, lava flow invasion caused major damage on the island, and as the population and infrastructure have increased dramatically since the last eruption, lava flows are the most important short-term volcanic risk on Tenerife. Hence, an understanding of lava flow behaviour is vital to manage risks from lava flows and minimise future losses on the island. This paper focuses on the lava flows from the historic eruptions in Tenerife, providing new data on the volumes emitted, advance rates and the timing of the emplacement of flows. The studies show three main stages in the development of unconfined flow fields: the first stage, corresponding to the fast advance of the initial fronts during the first 24?C36?h of eruption (reaching calculated velocities of up to 1.1?m/s); the second stage, in which fronts stagnate; and a third stage, in which secondary lava flows develop from breakouts 4?C7?days after the initial eruption and farther extend the flow field (velocities of up to 0.02?m/s have been calculated for this stage). The breakouts identified originated at sites both proximal and distal to the vent and, in both cases, caused damage through lengthening and widening the original flow field. Hence, the probability of damage from lavas to land and property is highest during stages 1 and 3, and this should be accounted for when planning the response to a future effusive eruption. Tenerife??s lava flows display a similar behaviour to that of lava flows on volcanoes characterised by basaltic effusive activity (such as Etna or Kilauea), indicating the possibility of applying forecasting models developed at those frequently active volcanoes to Tenerife.     

Tsunamis of volcanic origin: Summary of causes,with particular reference to Krakatoa, 1883

Known tsunamis of volcanic origin are reviewed and classified according to their causes. Earthquakes accompanying eruptions (excluding tectonic events which apparently triggered eruptions), pyroclastic flows, and submarine explosions have each accounted for about 20% of cases. Ten causes of volcanic tsunamis are discussed. From the risk point of view, those due to landslides are particularly dangerous. Eruptions at calderas are more likely to generate tsunamis than eruptions elsewhere. Of those killed directly by volcanic eruptions, nearly a quarter have died as a result of tsunamis. By transfer of energy to sea waves, a violent eruption, which would be comparatively harmless on land, extends greatly the radius over which destruction occurs. Krakatoa, 1883, is the only eruption sequence for which sufficient data exist for a detailed study of tsunamis. The times at which air and water waves generated by this sequence were recorded have been reread, and new origin times have been calculated and compared with observations made at the time. Origin times of successive pairs of air and water waves agree closely, except in some cases in which the tsunami arrived up to 15 minutes early, thus giving an apparent origin time 15 minutes before that of the corresponding air wave. This is explained by postulating that these tsunamis did not originate at the focus of the explosions, but at distances along the path towards the tide gauge, equivalent to those which would be covered by a tsunami in the time interval observed. The calculated point at which the largest recorded tsunami originated coincides with the outer edge of a bank of volcanic debris laid down during the eruption. This is interpreted as part of an unwelded ignimbrite deposit, the violent emplacement of which, within a minute or so of the explosion, generated the tsunami. A satisfactory correlation is established between explosions and deposits laid down by the eruptions, as described from a geological section close to the source vent. An outline is given of a proposed numerical index to define tsunamigenic potential at a given volcano. Such an index could be used to calculate the expected amplitudes of tsunamis at particular places in the vicinity, and hence could serve as a basis for tsunami risk contingency planning.     

Volcanism in Iceland in historical time: Volcano types,eruption styles and eruptive history

The large-scale volcanic lineaments in Iceland are an axial zone, which is delineated by the Reykjanes, West and North Volcanic Zones (RVZ, WVZ, NVZ) and the East Volcanic Zone (EVZ), which is growing in length by propagation to the southwest through pre-existing crust. These zones are connected across central Iceland by the Mid-Iceland Belt (MIB). Other volcanically active areas are the two intraplate belts of Öræfajökull (ÖVB) and Snæfellsnes (SVB). The principal structure of the volcanic zones are the 30 volcanic systems, where 12 are comprised of a fissure swarm and a central volcano, 7 of a central volcano, 9 of a fissure swarm and a central domain, and 2 are typified by a central domain alone.Volcanism in Iceland is unusually diverse for an oceanic island because of special geological and climatological circumstances. It features nearly all volcano types and eruption styles known on Earth. The first order grouping of volcanoes is in accordance with recurrence of eruptions on the same vent system and is divided into central volcanoes (polygenetic) and basalt volcanoes (monogenetic). The basalt volcanoes are categorized further in accordance with vent geometry (circular or linear), type of vent accumulation, characteristic style of eruption and volcanic environment (i.e. subaerial, subglacial, submarine).Eruptions are broadly grouped into effusive eruptions where >95% of the erupted magma is lava, explosive eruptions if >95% of the erupted magma is tephra (volume calculated as dense rock equivalent, DRE), and mixed eruptions if the ratio of lava to tephra occupy the range in between these two end-members. Although basaltic volcanism dominates, the activity in historical time (i.e. last 11 centuries) features expulsion of basalt, andesite, dacite and rhyolite magmas that have produced effusive eruptions of Hawaiian and flood lava magnitudes, mixed eruptions featuring phases of Strombolian to Plinian intensities, and explosive phreatomagmatic and magmatic eruptions spanning almost the entire intensity scale; from Surtseyan to Phreatoplinian in case of “wet” eruptions and Strombolian to Plinian in terms of “dry” eruptions. In historical time the magma volume extruded by individual eruptions ranges from ∼1 m³ to ∼20 km³ DRE, reflecting variable magma compositions, effusion rates and eruption durations.All together 205 eruptive events have been identified in historical time by detailed mapping and dating of events along with extensive research on documentation of eruptions in historical chronicles. Of these 205 events, 192 represent individual eruptions and 13 are classified as “Fires”, which include two or more eruptions defining an episode of volcanic activity that lasts for months to years. Of the 159 eruptions verified by identification of their products 124 are explosive, effusive eruptions are 14 and mixed eruptions are 21. Eruptions listed as reported-only are 33. Eight of the Fires are predominantly effusive and the remaining five include explosive activity that produced extensive tephra layers. The record indicates an average of 20–25 eruptions per century in Iceland, but eruption frequency has varied on time scale of decades. An apparent stepwise increase in eruption frequency is observed over the last 1100 years that reflects improved documentation of eruptive events with time. About 80% of the verified eruptions took place on the EVZ where the four most active volcanic systems (Grímsvötn, Bárdarbunga–Veidivötn, Hekla and Katla) are located and 9%, 5%, 1% and 0.5% on the RVZ–WVZ, NVZ, ÖVB, and SVB, respectively. Source volcano for ∼4.5% of the eruptions is not known.Magma productivity over 1100 years equals about 87 km³ DRE with basaltic magma accounting for about 79% and intermediate and acid magma accounting for 16% and 5%, respectively. Productivity is by far highest on the EVZ where 71 km³ (∼82%) were erupted, with three flood lava eruptions accounting for more than one half of that volume. RVZ–WVZ accounts for 13% of the magma and the NWZ and the intraplate belts for 2.5% each. Collectively the axial zone (RVZ, WVZ, NVZ) has only erupted 15–16% of total magma volume in the last 1130 years.     

An improved age framework for late Quaternary silicic eruptions in northern Central America

 Five new stepwise-heating ⁴⁰Ar/³⁹Ar ages and one new high-sensitivity ¹⁴C date of ash-fall and ash-flow deposits from late Quaternary silicic volcanoes in northern Central America document the eruption rates and frequencies of five major rhyodacite and rhyolite calderas (Atitlán, Amatitlán, Ayarza, Coatepeque, and Ilopango) located north of the basalt, andesite, and dacite stratovolcanoes of the Central American volcanic front. These deposits form extensive time-stratigraphic horizons that intercalate regionally, and knowledge of dates and stratigraphy provides a valuable framework for age determinations of more localized volcanic and nonvolcanic events. The new data, especially when integrated with previous stratigraphic and dating work, show that all five calderas erupted several times in the past 200 ka and, despite a lack of historic activity, should be considered as active centers that could produce highly explosive eruptions again. Because of their locations near the highly vulnerable economic hearts of Guatemala and El Salvador, the risks of eruptions from these calderas should be carefully considered along with risks of major earthquakes and volcanic front volcanoes, which are much more frequent but inflict less severe and extensive damage. This investigation also includes some examples of dating efforts that failed to produce reasonable results. Received: 15 May 1998 / Accepted: 18 January 1999     

Magma-hydrothermal system interaction inferred from volcanic gas measurements obtained during 2003–2008 at Meakandake volcano,Hokkaido, Japan

Phreatic eruptions occurred at the Meakandake volcano in 1988, 1996, 1998, 2006, and 2008. We conducted geochemical surveillance that included measurements of temperature, SO₂ emission rates, and volcanic gas composition from 2003 to 2008 at the Nakamachineshiri (NM), Northwest (NW), and Akanuma (AK) fumarolic areas, and the 96–1 vent, where historical eruptions had occurred. The elemental compositions of the gases discharged from the different areas are similar compared with the large variations observed in volcanic gases discharged from subduction zones. All the gases showed high apparent equilibrium temperatures, suggesting that all these gases originated from a common magmatic gas. The gases discharged from each area also exhibited different characteristics, which are probably the results of differences in the conditions of meteoric water mixing, quenching of chemical reactions, and vapor-liquid separation. The highest apparent equilibrium temperatures (about 500°C) were observed in the case of NW fumarolic gases, despite the low outlet temperature of about 100°C at these fumaroles. Since the NW fumaroles were formed as a result of the 2006 phreatic eruption, the high-temperature gas supply to the NW fumarole suggests that the phreatic eruption was caused by the ascent of high-temperature magmatic gases. The temperatures, compositions, and emission rates of the NM and 96–1 gases did not show any appreciable change after the 2006 eruption, indicating that each fumarolic system had a separate magmatic-hydrothermal system. The temperatures, compositions, and emission rates of the NM fumarolic gases were apparently constant, and these fumaroles are inferred to be formed by the evaporation of a hydrothermal system with a constant temperature of about 300°C. The 96–1 gas compositions showed large changes during continuous temperature decrease from 390° to 190°C occurred from 2003 to 2008, but the sulfur gas emission rates were almost constant at about four tons/day. At the 96–1 vent, the SO₂/H₂S ratio decreased, while the H₂/H₂O ratio remained almost constant; this was probably caused by the rock-buffer controlled chemical reaction during the temperature decrease.     

The dynamics and thermodynamics of large ash flows

 Ash flow deposits, containing up to 1000 km³ of material, have been produced by some of the largest volcanic eruptions known. Ash flows propagate several tens of kilometres from their source vents, produce extensive blankets of ash and are able to surmount topographic barriers hundreds of metres high. We present and test a new model of the motion of such flows as they propagate over a near horizontal surface from a collapsing fountain above a volcanic vent. The model predicts that for a given eruption rate, either a slow (10–100 m/s) and deep (1000–3000 m) subcritical flow or a fast (100–200 m/s) and shallow (500–1000 m) supercritical flow may develop. Subcritical ash flows propagate with a nearly constant volume flux, whereas supercritical flows entrain air and become progressively more voluminous. The run-out distance of such ash flows is controlled largely by the mass of air mixed into the collapsing fountain, the degree of fragmentation and the associated rate of loss of material into an underlying concentrated depositional system, and the mass eruption rate. However, in supercritical flows, the continued entrainment of air exerts a further important control on the flow evolution. Model predictions show that the run-out distance decreases with the mass of air entrained into the flow. Also, the mass of ash which may ascend from the flow into a buoyant coignimbrite cloud increases as more air is entrained into the flow. As a result, supercritical ash flows typically have shorter runout distances and more ash is elutriated into the associated coignimbrite eruption columns. We also show that one-dimensional, channellized ash flows typically propagate further than their radially spreading counterparts. As a Plinian eruption proceeds, the erupted mass flux often increases, leading to column collapse and the formation of pumiceous ash flows. Near the critical conditions for eruption column collapse, the flows are shed from high fountains which entrain large quantities of air per unit mass. Our model suggests that this will lead to relatively short ash flows with much of the erupted material being elutriated into the coignimbrite column. However, if the mass flux subseqently increases, then less air per unit mass is entrained into the collapsing fountain, and progressively larger flows, which propagate further from the vent, will develop. Our model is consistent with observations of a number of pyroclastic flow deposits, including the 1912 eruption of Katmai and the 1991 eruption of Pinatubo. The model suggests that many extensive flow sheets were emplaced from eruptions with mass fluxes of 10⁹–10¹⁰ kg/s over periods of 10³–10⁵ s, and that some indicators of flow "mobility" may need to be reinterpreted. Furthermore, in accordance with observations, the model predicts that the coignimbrite eruption columns produced from such ash flows rose between 20 and 40 km. Received: 25 August 1995 / Accepted: 3 April 1996     

Apoyo caldera, Nicaragua: A major quaternary silicic eruptive center

Apoyo caldera, near Granada, Nicaragua, was formed by two phases of collapse following explosive eruptions of dacite pumice about 23,000 yr B.P. The caldera sits atop an older volcanic center consisting of lava flows, domes, and ignimbrite (ash-flow tuff). The earliest lavas erupted were compositionally homogeneous basalt flows, which were later intruded by small andesite and dacite flows along a well defined set of N—S-trending regional faults. Collapse of the roof of the magma chamber occurred along near-vertical ring faults during two widely separated eruptions. Field evidence suggests that the climactic eruption sequence opened with a powerful plinian blast, followed by eruption column collapse, which generated a complex sequence of pyroclastic surge and ignimbrite deposits and initiated caldera collapse. A period of quiescence was marked by the eruption of scoria-bearing tuff from the nearby Masaya caldera and the development of a soil horizon. Violent plinian eruptions then resumed from a vent located within the caldera. A second phase of caldera collapse followed, accompanied by the effusion of late-stage andesitic lavas, indicating the presence of an underlying zoned magma chamber. Detailed isopach and isopleth maps of the plinian deposits indicate moderate to great column heights and muzzle velocities compared to other eruptions of similar volume. Mapping of the Apoyo airfall and ignimbrite deposits gives a volume of 17.2 km³ within the 1-mm isopach. Crystal concentration studies show that the true erupted volume was 30.5 km³ (10.7 km³ Dense Rock Equivalent), approximately the volume necessary to fill the caldera. A vent area located in the northeast quadrant of the present caldera lake is deduced for all the silicic pyroclastic eruptions. This vent area is controlled by N—S-trending precaldera faults related to left-lateral motion along the adjacent volcanic segment break. Fractional crystallization of calc-alkaline basaltic magma was the primary differentiation process which led to the intermediate to silicic products erupted at Apoyo. Prior to caldera collapse, highly atypical tholeiitic magmas resembling low-K, high-Ca oceanic ridge basalts were erupted along tension faults peripheral to the magma chamber. The injection of tholeiitic magmas may have contributed to the paroxysmal caldera-forming eruptions.