Showing papers in "Bulletin of Volcanology in 2013"

Origin and ascent history of unusually crystal-rich alkaline basaltic magmas from the western Pannonian Basin

Abstract: The last eruptions of the monogenetic Bakony-Balaton Highland Volcanic Field (western Pannonian Basin, Hungary) produced unusually crystal- and xenolith-rich alkaline basalts which are unique among the alkaline basalts of the Carpathian–Pannonian Region. Similar alkaline basalts are only rarely known in other volcanic fields of the world. These special basaltic magmas fed the eruptions of two closely located volcanic centres: the Bondoro-hegy and the Fuzes-to scoria cone. Their uncommon enrichment in diverse crystals produced unique rock textures and modified original magma compositions (13.1–14.2 wt.% MgO, 459–657 ppm Cr, and 455–564 ppm Ni contents). Detailed mineral-scale textural and chemical analyses revealed that the Bondoro-hegy and Fuzes-to alkaline basaltic magmas have a complex ascent history, and that most of their minerals (∼30 vol.% of the rocks) represent foreign crystals derived from different levels of the underlying lithosphere. The most abundant xenocrysts, olivine, orthopyroxene, clinopyroxene, and spinel, were incorporated from different regions and rock types of the subcontinental lithospheric mantle. Megacrysts of clinopyroxene and spinel could have originated from pegmatitic veins/sills which probably represent magmas crystallized near the crust–mantle boundary. Green clinopyroxene xenocrysts could have been derived from lower crustal mafic granulites. Minerals that crystallized in situ from the alkaline basaltic melts (olivine with Cr-spinel inclusions, clinopyroxene, plagioclase, and Fe–Ti oxides) are only represented by microphenocrysts and overgrowths on the foreign crystals. The vast amount of peridotitic (most common) and mafic granulitic materials indicates a highly effective interaction between the ascending magmas and wall rocks at lithospheric mantle and lower crustal levels. However, fragments from the middle and upper crust are absent from the studied basalts, suggesting a change in the style (and possibly rate) of magma ascent in the crust. These xenocryst- and xenolith-rich basalts yield divers tools for estimating magma ascent rate that is important for hazard forecasting in monogenetic volcanic fields. According to the estimated ascent rates, the Bondoro-hegy and Fuzes-to alkaline basaltic magmas could have reached the surface within hours to few days, similarly to the estimates for other eruptive centres in the Pannonian Basin which were fed by “normal” (crystal and xenoliths poor) alkaline basalts.

The potential failure of Monte Nuovo at Ischia Island (Southern Italy): numerical assessment of a likely induced tsunami and its effects on a densely inhabited area

Abstract: Ischia is the emergent top of a large volcanic complex that rises more than 1,000 m above the sea floor, at the north-western end of the Gulf of Naples. Caldera resurgence in the central part of the island has resulted in the formation of differentially displaced blocks, among which Mt. Epomeo (787 m a.s.l.) is the most uplifted. Deformation and slope instability have been recognised as common features induced by a block resurgence mechanism that causes uplift and favours gravitational loading and flank failure. The Monte Nuovo block, a topographic high on the north-western flank of Mt. Epomeo, has recently been interpreted as a block affected by deep-seated gravitational slope deformation. This block may undergo a catastrophic failure in the case of renewal of magmatic activity. This paper investigates the potential failure of the Monte Nuovo block as a rockslide-debris avalanche, the consequent tsunami generation and wave propagation, and discusses the catastrophic effects of such an event. Mobilization-prone volume has been estimated at about 160·106 m3 and would move from a maximum elevation of 400 m a.s.l. The landslide itself would sweep away a densely populated territory as large as 3.5 km2. The highest waves generated by the tsunami, on which this paper is mainly focussed, would hit the northern and western shores of Ischia. However, the high coast would prevent inundation and limit devastation to beaches, harbours and surrounding areas. Most of the tsunami energy would head towards the north-east, hitting the Campania coast. Severe inundation would affect an area of up to 20 km2 around the mouth of the Volturno river, including the urban area of Castel Volturno. In contrast, less energy would travel towards the south, and the Gulf of Naples would be perturbed by long persisting waves of limited damaging potential.

Storage and eruption of near-liquidus rhyolite magma at Cordon Caulle, Chile

Abstract: The last three eruptions at the Cordon Caulle volcanic complex, Chile, have been strikingly similar in that they have started with relatively short pre-eruptive warning and produced chemically homogeneous rhyolite to rhyodacite magma with glassy to aphyric texture. These characteristics collectively call for an understanding of the storage conditions leading to the rise and extraction of crystal-poor silicic magma from volcanoes. We have analyzed and experimentally reproduced the mineral assemblage and glass chemistry in rhyolite magma produced in the most recent eruption of Cordon Caulle, and we use these to infer magma storage and ascent conditions. Fe–Ti oxide mineral geothermometry suggests that the rhyolite was stored at ∼870–920 °C. At these temperatures, the phenocryst assemblage (plag∼An37 > cpx + opx > mag + ilm) can be reproduced under H2O-saturated conditions of between 100 and 50 MPa, corresponding to crustal depths between about 2.5 and 5.0 km. The shallow and relatively hot magma storage conditions have implications for the rapid onset, degassing efficiency, and progression from explosive to mixed pyroclastic-effusive eruption style at Cordon Caulle.

Plume height, volume, and classification of explosive volcanic eruptions based on the Weibull function

Abstract: The Weibull distribution between volume and square root of isopach area has been recently introduced for determining volume of tephra deposits, which is crucial to the assessment of the magnitude and hazards of explosive volcanoes. We show how the decay of the size of the largest lithics with the square root of isopleth area can also be well described using a Weibull function and how plume height correlates strongly with corresponding Weibull parameters. Variations of median grain size (Mdϕ) values with square root of area of the associated contours can be, similarly, well fitted with a Weibull function. Weibull parameters, derived for both the thinning of tephra deposits and the decrease of grain size (both maximum lithic diameter and Mdϕ), with a proxy for the distance from vent (e.g., square root of isoline areas) can be combined to classify the style of explosive volcanic eruptions. Accounting for the uncertainty in the derivation of eruptive parameters (e.g., plume height and volume of tephra deposits) is crucial to any classification of eruptive style and hazard assessment. Considering a typical uncertainty of 20 % for the determination of plume height, a new eruption classification scheme based on selected Weibull parameters is proposed. Ultraplinian, Plinian, Subplinian, and small–moderate explosive eruptions are defined on the ground of plume height and mass eruption rate. Overall, the Weibull fitting represents a versatile and reliable strategy for the estimation of both the volume of tephra deposits and the height of volcanic plumes and for the classification of eruptive style. Nonetheless, due to the typically large uncertainties (mainly due to availability of data, compilation of isopach and isopleth maps, and discrepancies from empirical best fits), plume height, volume, and magnitude of explosive eruptions cannot be considered as absolute values, regardless of the technique used. It is important that various empirical and analytical methods are applied in order to assess such an uncertainty.

Review of eruptive activity at Tianchi volcano, Changbaishan, northeast China: implications for possible future eruptions

Abstract: One of the largest explosive eruptions in the past several thousand years occurred at Tianchi volcano, also known as Changbaishan, on the China–North Korea border. This historically active polygenetic central volcano consists of three parts: a lower basaltic shield, an upper trachytic composite cone, and young comendite ash flows. The Millennium Eruption occurred between 938 and 946 ad, and was preceded by two smaller and chemically different rhyolitic pumice deposits. There has been at least one additional, small eruption in the last three centuries. From 2002 to 2005, seismicity, deformation, and the helium and hydrogen gas contents of spring waters all increased markedly, causing regional concern. We attribute this event to magma recharge or volatile exhalation or both at depth, followed by two episodes of addition of magmatic fluids into the overlying aquifer without a phreatic eruption. The estimated present magma accumulation rate is too low by itself to account for the 2002–2005 unrest. The most serious volcanic hazards are ash eruption and flows, and lahars. The available geological information and volcano monitoring data provide a baseline for comprehensive assessment of future episodes of unrest and possible eruptive activity.

European summer temperature response to annually dated volcanic eruptions over the past nine centuries

Abstract: The drop in temperature following large volcanic eruptions has been identified as an important component of natural climate variability. However, due to the limited number of large eruptions that occurred during the period of instrumen- tal observations, the precise amplitude of post-volcanic cooling is not well constrained. Here we present new evidence on summer temperature cooling over Europe in years following volcanic eruptions. We compile and analyze an updated net- work of tree-ring maximum latewood density chronologies, spanning the past nine centuries, and compare cooling sig- natures in this network with exceptionally long instru- mental station records and state-of-the-art general circu- lation models. Results indicate post-volcanic June-August cooling is strongest in Northern Europe 2 years after an eruption (−0.52±0.05 °C), whereas in Central Europe the temperature response is smaller and occurs 1 year after an eruption (−0.18±0.07 °C). We validate these estimates by comparisonwith the shorter instrumental network and evaluate the statistical significance of post-volcanic summer tempera- ture cooling in the context of natural climate variability over the past nine centuries. Finding no significant post-volcanic temperature cooling lasting longer than 2 years, our results question the ability of large eruptions to initiate long-term temperature changes through feedback mechanisms in the climate system. We discuss the implications of these findings with respect to the response seen in general circulation models and emphasize the importance of considering well- documented, annually dated eruptions when assessing the

Determination of the largest clast sizes of tephra deposits for the characterization of explosive eruptions: a study of the IAVCEI commission on tephra hazard modelling

Abstract: The distribution of clasts deposited around a volcano during an explosive eruption typically contoured by isopleth maps provides important insights into the associated plume height, wind speed and eruptive style. Nonetheless, a wide range of strategies exists to determine the largest clasts, which can lead to very different results with obvious implications for the characterization of eruptive behaviour of active volcanoes. The IAVCEI Commission on Tephra Hazard Modelling has carried out a dedicated exercise to assess the influence of various strategies on the determination of the largest clasts. Suggestions on the selection of sampling area, collection strategy, choice of clast typologies and clast characterization (i.e. axis measurement and averaging technique) are given, mostly based on a thorough investigation of two outcrops of a Plinian tephra deposit from Cotopaxi volcano (Ecuador) located at different distances from the vent. These include: (1) sampling on a flat paleotopography far from significant slopes to minimize remobilization effects; (2) sampling on specified-horizontal-area sections (with the statistically representative sampling area depending on the outcrop grain size and lithic content); (3) clast characterization based on the geometric mean of its three orthogonal axes with the approximation of the minimum ellipsoid (lithic fragments are better than pumice clasts when present); and (4) use of the method of the 50th percentile of a sample of 20 clasts as the best way to assess the largest clasts. It is also suggested that all data collected for the construction of isopleth maps be made available to the community through the use of a standardized data collection template, to assess the applicability of the new proposed strategy on a large number of deposits and to build a large dataset for the future development and refinement of dispersal models.

Compositionally zoned crystals and real-time degassing data reveal changes in magma transfer dynamics during the 2006 summit eruptive episodes of Mt. Etna

Abstract: One of the major objectives of volcanology remains relating variations in surface monitoring signals to the magmatic processes at depth that cause these variations. We present a method that enables compositional and temporal information stored in zoning of minerals (olivine in this case) to be linked to observations of real-time degassing data. The integrated record may reveal details of the dynamics of gradual evolution of a plumbing system during eruption. We illustrate our approach using the 2006 summit eruptive episodes of Mt. Etna. We find that the history tracked by olivine crystals, and hence, most likely the magma pathways within the shallow plumbing system of Mt. Etna, differed considerably between the July and October eruptions. The compositional and temporal record preserved in the olivine zoning patterns reveal two mafic recharge events within months of each other (June and September 2006), and each of these magma supplies may have triggered the initiation of different eruptive cycles (July 14–24 and August 31–December 14). Correlation of these observations with gas monitoring data shows that the systematic rise of the CO2/SO2 gas values is associated with the gradual (pre- and syn-eruptive) supply of batches of gas-rich mafic magma into segments of Etna’s shallow plumbing system, where mixing with pre-existing and more evolved magma occurred.

Acute sedimentation response to rainfall following the explosive phase of the 2008-2009 eruption of Chaitén volcano, Chile

Abstract: The 10-day explosive phase at the start of the 2008–2009 eruption of Chaiten volcano in southern Chile (42.83°S, 72.65°W) blanketed the steep, rain-forest-cloaked, 77-km2 Chaiten River drainage basin with 3 to >100 cm of tephra; predominantly fine to extremely fine rhyolitic ash fell during the latter half of the explosive phase. Rain falling on this ash blanket within days of cessation of major explosive activity generated a hyperconcentrated-flow lahar, followed closely by a complex, multi-day, muddy flood (streamflow bordering on dilute hyperconcentrated flow). Sediment mobilized in this lahar-flood event filled the Chaiten River channel with up to 7 m of sediment, buried the town of Chaiten (10 km downstream of the volcano) in up to 3 m of sediment, and caused the lower 3 km of the channel to avulse through the town. Although neither the nature nor rate of the sedimentation response is unprecedented, they are unusual in several ways: (1) Nearly 70 percent of the aggradation (almost 5 m) in the 50–70-m-wide Chaiten River channel was caused by a lahar, triggered by an estimated 20 mm of rainfall over a span of about 24 h. An additional 2 m of aggradation occurred in the next 24–36 h. (2) Direct damage to the town was accomplished by the sediment-laden water-flood phase of the lahar-flood event, not the lahar phase. (3) The volume of sediment eroded from hillslopes and delivered to the Chaiten River channel was at least 3–8 × 106 m3—roughly 15–40 % of the minimum tephra volume that mantled the Chaiten River drainage basin. (4) The acute sedimentation response to rainfall appears to have been due to the thickness and fineness of the ash blanket (inhibiting infiltration of rain) and the steepness of the basin’s hillslopes. Other possible factors such as the prior formation of an ash crust, development of a hydrophobic surface layer, or large-scale destruction of rain-intercepting vegetation did not play a role.

The nature and formation of cristobalite at the Soufrière Hills volcano, Montserrat: implications for the petrology and stability of silicic lava domes

Abstract: Cristobalite is commonly found in the dome lava of silicic volcanoes but is not a primary magmatic phase; its presence indicates that the composition and micro-structure of dome lavas evolve during, and after, emplacement. Nine temporally and mineralogically diverse dome samples from the Soufriere Hills volcano (SHV), Montserrat, are analysed to provide the first detailed assessment of the nature and mode of cristobalite formation in a volcanic dome. The dome rocks contain up to 11 wt.% cristobalite, as defined by X-ray dif- fraction. Prismatic and platy forms of cristobalite, identified by scanning electron microscopy (SEM), are commonly found in pores and fractures, suggesting that they have precipitated from a vapour phase. Feathery crystallites and micro-crystals of cristobalite and quartz associated with vol- canic glass, identified using SEM-Raman, are interpreted to have formed by varying amounts of devitrification. We dis- cussmechanismsof silicatransport and cristobaliteformation, and their implications for petrological interpretations and domestability. We conclude: (1) that silica maybetransported in the vapour phase locally, or from one part of the magmatic system to another; (2) that the potential for transport of silica into the dome should not be neglected in petrological and geochemical studies because the addition of non-magmatic phases may affect whole rock composition; and (3) that the extent of cristobalite mineralisation in the dome at SHV is sufficient to reduce porosity—hence, permeability—and may impact on the mechanical strength of the dome rock, thereby potentially affecting dome stability.

Rhyolitic volcanism of the central Snake River Plain: a review

Abstract: The central Snake River Plain (CSRP) of southern Idaho and northern Nevada, USA, forms part of the Columbia River–Yellowstone large igneous province. Volcanic rocks of the province are compositionally bimodal (basalt–rhyolite), and the rhyolites produce a broadly time-transgressive record of a hotspot which is currently located under Yellowstone. Snake River Plain rhyolites represent hot (>850 °C), dry magmas and have field characteristics consistent with high emplacement temperatures. Individual ignimbrite sheets reach 1,000 km3 and exhibit little to no compositional zonation on a large scale but reveal considerable complexity on a crystal scale, particularly with regard to pyroxene compositions. Multiple pyroxene compositions may exist in a single ignimbrite which, along with multiple glass compositions in widely dispersed fallout tephra, suggests complex storage of rhyolite prior to eruption. Unlike most igneous rocks, the mineral cargo of the CSRP rhyolites exhibits little isotopic variability, with unimodal 87Sr/86Sr values returned from plagioclase grains inferred to represent the combination of strong crystal–melt coupling and rapid diffusional re-equilibriation. All the rhyolites within the CSRP have a characteristic low-δ18O signature; with >20,000 km3 of rhyolite exhibiting this depletion, the CSRP represents the largest low-δ18O province on Earth. The low-18O nature of the rhyolites requires assimilation of hydrothermally altered materials which may be from altered Eocene batholithic rocks or from down-dropped intra-caldera tuffs. The wide range of crustal assimilants, with highly variable radiogenic isotope characteristics, available in the CSRP is permissive of a variety of petrogenetic models based on radiogenic isotopic data.

Unbedded diatreme deposits reveal maar-diatreme-forming eruptive processes: Standing Rocks West, Hopi Buttes, Navajo Nation, USA

Abstract: Diatremes provide partial records of how dyke-delivered magma, periodically interacting with water, produces the largest known cylindrical conduit structures of any volcano type. We address how pre-eruptive country rock is disrupted and redistributed to form the diatreme structure during an eruption by establishing the internal architecture of the unbedded diatreme at Standing Rocks West, and the volumes and sources of wall-rock within the diatreme and in a complementary tephra ring at Teshim, both in the Hopi Buttes volcanic field. The unbedded diatreme is dominantly exposed as a massif comprising multiple subvertical columns, interpreted as ephemeral-conduit deposits of well-mixed, poorly sorted, composite and juvenile pyroclast-rich deposits that truncate marginal layered deposits, plus a peripheral heterolithic country rock breccia. Wall-rock clasts are scarce in the massif and were mainly sourced from ~75 to 245 m depth below the paleosurface. By contrast, the tephra ring has abundant wall-rock fragments, dominantly from the upper 75 m of the pre-eruptive sedimentary sequence. The diatreme structure is interpreted to have been formed by many small-volume explosions fed by a low flux of basaltic magma. Late stage explosive activity was rooted mostly within pyroclastic debris at shallow to intermediate depths within the diatreme structure. This resulted in damped and shifting intra-diatreme explosions and jets that facilitated gradual mixing, recycling and remobilization of debris in the diatreme, with incremental addition of juvenile material and possibly a local rise in the crater floor.

Evolution of the crustal magma plumbing system during the build-up to the 22-ka caldera-forming eruption of Santorini (Greece)

Abstract: The formation of shallow, caldera-sized reservoirs of crystal-poor silicic magma requires the generation of large volumes of silicic melt, followed by the segregation of that melt and its accumulation in the upper crust. The 21.8 ± 0.4-ka Cape Riva eruption of Santorini discharged >10 km3 of crystal-poor dacitic magma, along with 0.004 km3 year-1. Rapid ascent and accumulation of the Cape Riva dacite may have been caused by an increased flux of mantle-derived basalt into the crust, explaining the occurrence of hybrid andesites (formed by the mixing of olivine basalt and dacite in approximately equal proportions) in the Cape Riva and late Therasia products. Pressurisation of the upper crustal plumbing system by sustained, high-flux injection of dacite and basalt may have triggered the transition from prolonged, largely effusive activity to explosive eruption and caldera collapse.

Spherulites and lithophysae—200 years of investigation on high-temperature crystallization domains in silica-rich volcanic rocks

Abstract: High-temperature crystallization domains (HTCDs) including spherulites and lithophysae form during cooling of silica-rich lava and welded ignimbrites. Spherulites grow in silicate melts or hot glass and they display a radiating or microcrystalline texture, typically consisting of cristobalite, tridymite, and sanidine. Lithophysae are HTCDs comprising one or more cavities. This contribution reviews the research and discussions on HTCDs carried out over the last 200 years. The emphasis, here, is on lithophysae and summarizes current knowledge of their formation. A number of parameters influence the initiation and growth of lithophysae, as well as, their shapes and internal textures. The most likely cause of cavity formation is transient tensional stress that produces a mechanical opening and widening at the interface between the crystallization front and the host melt (e.g., where T > T g ). Cavity growth and expansion forced by rising vapor pressure is considered less important. In some cases, further growth of HTCD cavities results from vapor phase corrosion and brecciation.

Architecture and emplacement of flood basalt flow fields: case studies from the Columbia River Basalt Group, NW USA

Abstract: The physical features and morphologies of collections of lava bodies emplaced during single eruptions (known as flow fields) can be used to understand flood basalt emplacement mechanisms. Characteristics and internal features of lava lobes and whole flow field morphologies result from the forward propagation, radial spread, and cooling of individual lobes and are used as a tool to understand the architecture of extensive flood basalt lavas. The features of three flood basalt flow fields from the Columbia River Basalt Group are presented, including the Palouse Falls flow field, a small (8,890 km2, ∼190 km3) unit by common flood basalt proportions, and visualized in three dimensions. The architecture of the Palouse Falls flow field is compared to the complex Ginkgo and more extensive Sand Hollow flow fields to investigate the degree to which simple emplacement models represent the style, as well as the spatial and temporal developments, of flow fields. Evidence from each flow field supports emplacement by inflation as the predominant mechanism producing thick lobes. Inflation enables existing lobes to transmit lava to form new lobes, thus extending the advance and spread of lava flow fields. Minimum emplacement timescales calculated for each flow field are 19.3 years for Palouse Falls, 8.3 years for Ginkgo, and 16.9 years for Sand Hollow. Simple flow fields can be traced from vent to distal areas and an emplacement sequence visualized, but those with multiple-layered lobes present a degree of complexity that make lava pathways and emplacement sequences more difficult to identify.

Magma discharge variations during the 2011 eruptions of Shinmoe-dake volcano, Japan, revealed by geodetic and satellite observations

Abstract: We present precise geodetic and satellite observation-based estimations of the erupted volume and discharge rate of magma during the 2011 eruptions of Kirishima-Shinmoe-dake volcano, Japan. During these events, the type and intensity of eruption drastically changed within a week, with three major sub-Plinian eruptions on January 26 and 27, and a continuous lava extrusion from January 29 to 31. In response to each eruptive event, borehole-type tiltmeters detected deflation of a magma chamber caused by migration of magma to the surface. These measurements enabled us to estimate the geodetic volume change in the magma chamber caused by each eruptive event. Erupted volumes and discharge rates were constrained during lava extrusion using synthetic aperture radar satellite imaging of lava accumulation inside the summit crater. Combining the geodetic volume change and the volume of lava extrusion enabled the determination of the erupted volume and discharge rate during each sub-Plinian event. These precise estimates provide important information about magma storage conditions in magma chambers and eruption column dynamics, and indicate that the Shinmoe-dake eruptions occurred in a critical state between explosive and effusive eruption.

40 Ar/ 39 Ar dating, geochemistry, and isotopic analyses of the quaternary Chichinautzin volcanic field, south of Mexico City: implications for timing, eruption rate, and distribution of volcanism

Abstract: Monogenetic structures located at the southern and western ends of the Chichinautzin volcanic field (Trans-Mexican Volcanic Belt, Central Mexico) yield 40Ar/39Ar ages ranging from 1.2 Ma in the western portion of the field to 1.0–0.09 Ma in the southern portion, all of which are older than the <0.04 Ma age previously established for the entire volcanic field. These new ages indicate: (1) an eruption rate of 0.47 km3/kyr, which is much lower than the 11.7 km3/kyr previously estimated; (2) that the Chichinautzin magmatism coexisted with the Zempoala (0.7 Ma) and La Corona (1.0 Ma) polygenetic volcanoes on the southern edge of Las Cruces Volcanic Range (Trans-Mexican Volcanic Belt); and confirm (3) that the drainage system between the Mexico and Cuernavaca basins was closed during early Pleistocene forming the Texcoco Lake. Whole-rock chemistry and Sr, Nd, and Pb isotopic data indicate heterogeneous magmatism throughout the history of Chichinautzin activity that likely reflects variable degrees of slab and sediment contributions to the mantle wedge, fractional crystallization, and crustal assimilation. Even with the revised duration of volcanism within the Chichinautzin Volcanic Field, its eruption rate is higher than most other volcanic fields of the Trans-Mexican Volcanic Belt and is comparable only to the Tacambaro-Puruaran area in the Michoacan-Guanajuato Volcanic Field to the west. These variations in eruption rates among different volcanic fields may reflect a combination of variable subduction rates of the Rivera and Cocos plates along the Middle America Trench, as well as different distances from the trench, variations in the depth with respect to the subducted slab, or the upper plate characteristics.

The 1815 Tambora ash fall: implications for transport and deposition of distal ash on land and in the deep sea

Abstract: Tambora volcano lies on the Sanggar Peninsula of Sumbawa Island in the Indonesian archipelago. During the great 1815 explosive eruption, the majority of the erupted pyroclastic material was dispersed and subsequently deposited into the Indian Ocean and Java Sea. This study focuses on the grain size distribution of distal 1815 Tambora ash deposited in the deep sea compared to ash fallen on land. Grain size distribution is an important factor in assessing potential risks to aviation and human health, and provides additional information about the ash transport mechanisms within volcanic umbrella clouds. Grain size analysis was performed using high precision laser diffraction for a particle range of 0.2 μm–2 mm diameter. The results indicate that the deep-sea samples provide a smooth transition to the land samples in terms of grain size distributions despite the different depositional environments. Even the very fine ash fraction (<10 μm) is deposited in the deep sea, suggesting vertical density currents as a fast and effective means of transport to the seafloor. The measured grain size distribution is consistent with an improved atmospheric gravity current sedimentation model that takes into account the finite duration of an eruption. In this model, the eruption time and particle fall velocity are the critical parameters for assessing the ash component depositing while the cloud advances versus the ash component depositing once the eruption terminates. With the historical data on eruption duration (maximum 24 h) and volumetric flow rate of the umbrella cloud (∼1.5–2.5 × 1011 m3/s) as input to the improved model, and assuming a combination of 3 h Plinian phase and 21 h co-ignimbrite phase, it reduces the mean deviation of the predicted versus observed grain size distribution by more than half (∼9.4 % to ∼3.7 %) if both ash components are considered.

Unraveling the solidification path of a pahoehoe “cicirara” lava from Mount Etna volcano

Abstract: The solidified surface of a lava flow reflects the viscosity of its molten fraction and the crystal content during flow; crystal-poor basaltic lavas produce pahoehoe fields, whereas crystal-rich ones solidify with aa carapaces. At Mount Etna, volcano aa morphologies are very common, whereas pahoehoe lavas are rare. The latter are locally named “cicirara” due to the presence of centimeter-sized plagioclase phenocrysts much more abundant than in aa lavas. The phenocryst content of “cicirara” lavas contrasts with the low viscosity generally associated with pahoehoe morphology. Therefore, to reconcile the discrepancy between textural and volcanic observations, we have studied the most primitive pahoehoe “cicirara” lava sampled until now. Two samples at 0.5 and 1 m from the bottom of the 2-m thick lava flow were investigated on the basis of their mineral compositional variations and textural features, i.e., size frequency and crystal size distribution (CSD). Results coupled with rheological models indicate that only large phenocrysts of plagioclase (>1 mm) and clinopyroxene have grown before eruption. Thermobarometric models and petrological computations based on the composition of plagioclase and clinopyroxene phenocryst cores highlight that only a small amount (10–15 vol.%) of crystals equilibrated at 12 km of depth. Cumulative size frequency and CSD data also indicate that plagioclase and clinopyroxene phenocryst rims grew heterogeneously and coalesced around their cores at depths <1 km, before eruption. In this view, the “cicirara” lava was erupted with a low crystalline content that favoured the formation of its pahoehoe surface; however, crystals with a size <1 mm (~75 vol.%) solidified at post-eruptive conditions. Our findings underline that the emplacement of high-viscosity aa or low-viscosity pahoehoe lavas is driven by the degree of undercooling imposed by the volatile exsolution rate in the shallowest portion of the Etnean plumbing system. A slow magma ascent rate promotes significant intratelluric degassing and widespread nucleation; consequently, the viscosity of the suspension significantly increases leading to an aa morphology. In contrast, pahoehoe “cicirara” lavas are associated with a rapid rise to the surface of poorly degassed, undercooled magmas.

Eruption processes and deposit characteristics at the monogenetic Mt. Gambier Volcanic Complex, SE Australia: implications for alternating magmatic and phreatomagmatic activity

Abstract: The ∼5 ka Mt Gambier Volcanic Complex in the Newer Volcanics Province, Australia is an extremely complex monogenetic, volcanic system that preserves at least 14 eruption points aligned along a fissure system The complex stratigraphy can be subdivided into six main facies that record alternations between magmatic and phreatomagmatic eruption styles in a random manner The facies are (1) coherent to vesicular fragmental alkali basalt (effusive/Hawaiian spatter and lava flows); (2) massive scoriaceous fine lapilli with coarse ash (Strombolian fallout); (3) bedded scoriaceous fine lapilli tuff (violent Strombolian fallout); (4) thin–medium bedded, undulating very fine lapilli in coarse ash (dry phreatomagmatic surge-modified fallout); (5) palagonite-altered, cross-bedded, medium lapilli to fine ash (wet phreatomagmatic base surges); and (6) massive, palagonite-altered, very poorly sorted tuff breccia and lapilli tuff (phreato-Vulcanian pyroclastic flows) Since most deposits are lithified, to quantify the grain size distributions (GSDs), image analysis was performed The facies are distinct based on their GSDs and the fine ash to coarse+fine ash ratios These provide insights into the fragmentation intensities and water–magma interaction efficiencies for each facies The eruption chronology indicates a random spatial and temporal sequence of occurrence of eruption styles, except for a “magmatic horizon” of effusive activity occurring at both ends of the volcanic complex simultaneously The eruption foci are located along NW–SE trending lineaments, indicating that the complex was fed by multiple dykes following the subsurface structures related to the Tartwaup Fault System Possible factors causing vent migration along these dykes and changes in eruption styles include differences in magma ascent rates, viscosity, crystallinity, degassing and magma discharge rate, as well as hydrological parameters

Stress field change around the Mount Fuji volcano magma system caused by the Tohoku megathrust earthquake, Japan

Abstract: Crustal deformation by the M w 9.0 megathrust Tohoku earthquake causes the extension over a wide region of the Japanese mainland. In addition, a triggered M w 5.9 East Shizuoka earthquake on March 15 occurred beneath the south flank, just above the magma system of Mount Fuji. To access whether these earthquakes might trigger the eruption, we calculated the stress and pressure changes below Mount Fuji. Among the three plausible mechanisms of earthquake–volcano interactions, we calculate the static stress change around volcano using finite element method, based on the seismic fault models of Tohoku and East Shizuoka earthquakes. Both Japanese mainland and Mount Fuji region are modeled by seismic tomography result, and the topographic effect is also included. The differential stress given to Mount Fuji magma reservoir, which is assumed to be located to be in the hypocentral area of deep long period earthquakes at the depth of 15 km, is estimated to be the order of about 0.001–0.01 and 0.1–1 MPa at the boundary region between magma reservoir and surrounding medium. This pressure change is about 0.2 % of the lithostatic pressure (367.5 MPa at 15 km depth), but is enough to trigger an eruptions in case the magma is ready to erupt. For Mount Fuji, there is no evidence so far that these earthquakes and crustal deformations did reactivate the volcano, considering the seismicity of deep long period earthquakes.

Short time scales of magma-mixing processes prior to the 2011 eruption of Shinmoedake volcano, Kirishima volcanic group, Japan

Abstract: We estimated time scales of magma-mixing processes just prior to the 2011 sub-Plinian eruptions of Shinmoedake volcano to investigate the mechanisms of the triggering processes of these eruptions. The sequence of these eruptions serves as an ideal example to investigate eruption mechanisms because the available geophysical and petrological observations can be combined for interpretation of magmatic processes. The eruptive products were mainly phenocryst-rich (28 vol%) andesitic pumice (SiO2 57 wt%) with a small amount of more silicic pumice (SiO2 62–63 wt%) and banded pumice. These pumices were formed by mixing of low-temperature mushy silicic magma (dacite) and high-temperature mafic magma (basalt or basaltic andesite). We calculated the time scales on the basis of zoning analysis of magnetite phenocrysts and diffusion calculations, and we compared the derived time scales with those of volcanic inflation/deflation observations. The magnetite data revealed that a significant mixing process (mixing I) occurred 0.4 to 3 days before the eruptions (pre-eruptive mixing) and likely triggered the eruptions. This mixing process was not accompanied by significant crustal deformation, indicating that the process was not accompanied by a significant change in volume of the magma chamber. We propose magmatic overturn or melt accumulation within the magma chamber as a possible process. A subordinate mixing process (mixing II) also occurred only several hours before the eruptions, likely during magma ascent (syn-eruptive mixing). However, we interpret mafic injection to have begun more than several tens of days prior to mixing I, likely occurring with the beginning of the inflation (December 2009). The injection did not instantaneously cause an eruption but could have resulted in stable stratified magma layers to form a hybrid andesitic magma (mobile layer). This hybrid andesite then formed the main eruptive component of the 2011 eruptions of Shinmoedake.

Topographic and Stochastic Influences on Pahoehoe Lava Lobe Emplacement

Abstract: A detailed understanding of pāhoehoe emplacement is necessary for developing accurate models of flow field development, assessing hazards, and interpreting the significance of lava morphology on Earth and other planetary surfaces. Active pāhoehoe lobes on Kīlauea Volcano, Hawai'i, were examined on 21–26 February 2006 using oblique time series stereo-photogrammetry and differential global positioning system measurements. During this time, the local discharge rate for peripheral lava lobes was generally constant at 0.0061 ± 0.0019 m3/s, but the areal coverage rate of the lobes exhibited a periodic increase every 4.13 ± 0.64 min. This periodicity is attributed to the time required for the pressure within the liquid lava core to exceed the cooling-induced strength of its margins. The pāhoehoe flow advanced through a series of down-slope and cross-slope breakouts, which began as ∼0.2-m-thick units (i.e., toes) that coalesced and inflated to become approximately meter-thick lobes. The lobes were thickest above the lowest points of the initial topography and above shallow to reverse-facing slopes, defined relative to the local flow direction. The flow path was typically controlled by high-standing topography, with the zone directly adjacent to the final lobe margin having an average relief that was a few centimeters higher than the lava-inundated region. This suggests that toe-scale topography can, at least temporarily, exert strong controls on pāhoehoe flow paths by impeding the lateral spreading of the lobe. Observed cycles of enhanced areal spreading and inflated lobe morphology are also explored using a model that considers the statistical likelihood of sequential breakouts from active flow margins and the effects of topographic barriers.

The influence of thermal and cyclic stressing on the strength of rocks from Mount St. Helens, Washington

Abstract: Stratovolcanoes and lava domes are particularly susceptible to sector collapse resulting from wholesale rock failure as a consequence of decreasing rock strength. Here, we provide insights into the influence of thermal and cyclic stressing on the strength and mechanical properties of volcanic rocks. Specifically, this laboratory study examines the properties of samples from Mount St. Helens; chosen because its strength and stability have played a key role in its history, influencing the character of the infamous 1980 eruption. We find that thermal stressing exerts different effects on the strengths of different volcanic units; increasing the heterogeneity of rocks in situ. Increasing the uniaxial compressive stress generates cracking, the timing and magnitude of which was monitored via acoustic emission (AE) output during our experiments. AEs accelerated in the approach to failure, sometimes following the pattern predicted by the failure forecast method (Kilburn 2003). Crack damage during the experiments was tracked using the evolving static Young’s modulus and Poisson’s ratio, which represent the quasi-static deformation in volcanic edifices more accurately than dynamic elastic moduli which are usually implemented in volcanic models. Cyclic loading of these rocks resulted in a lower failure strength, confirming that volcanic rocks may be weakened by repeated inflation and deflation of the volcanic edifice. Additionally, volcanic rocks in this study undergo significant elastic hysteresis; in some instances, a material may fail at a stress lower than the peak stress which has previously been endured. Thus, a volcanic dome repeatedly inflated and deflated may progressively weaken, possibly inducing failure without necessarily exceeding earlier conditions.

Relating vesicle shapes in pyroclasts to eruption styles

Abstract: Vesicles in pyroclasts provide a direct record of conduit conditions during explosive volcanic eruptions. Although their numbers and sizes are used routinely to infer aspects of eruption dynamics, vesicle shape remains an underutilized parameter. We have quantified vesicle shapes in pyroclasts from fall deposits of seven explosive eruptions of different styles, using the dimensionless shape factor $\Omega $ , a measure of the degree of complexity of the bounding surface of an object. For each of the seven eruptions, we have also estimated the capillary number, Ca, from the magma expansion velocity through coupled diffusive bubble growth and conduit flow modeling. We find that $\Omega $ is smaller for eruptions with ${\rm{Ca}} \ll 1$ than for eruptions with Ca $\gg 1$ . Consistent with previous studies, we interpret these results as an expression of the relative importance of structural changes during magma decompression and bubble growth, such as coalescence and shape relaxation of bubbles by capillary stresses. Among the samples analyzed, Strombolian and Hawaiian fire-fountain eruptions have ${\rm{Ca}} \ll 1$ , in contrast to Vulcanian, Plinian, and ultraplinian eruptions. Interestingly, the basaltic Plinian eruptions of Tarawera volcano, New Zealand in 1886 and Mt. Etna, Italy in 122 BC, for which the cause of intense explosive activity has been controversial, are also characterized by ${\rm{Ca}} \gg 1$ and larger values of $\Omega $ than Strombolian and Hawaiian style (fire fountain) eruptions. We interpret this to be the consequence of syn-eruptive magma crystallization, resulting in high magma viscosity and reduced rates of bubble growth. Our model results indicate that during these basaltic Plinian eruptions, buildup of bubble overpressure resulted in brittle magma fragmentation.

Origin and evolution of the Deception Island caldera (South Shetland Islands, Antarctica)

Abstract: Deception Island has been interpreted variously as a classical ring fault caldera, as a tectonically controlled collapse caldera or as a tectonic depression. Review of previous studies combined with new fieldwork has allowed us to obtain a more precise model of the formation and internal structure of the Deception Island caldera. It formed as a result of the explosive eruption of basaltic-to-andesitic magmas, mostly as pyroclastic density currents representing in total a bulk volume of the order of 90 km3. Caldera collapse occurred rapidly along a polygonal structural network consisting of several pre-existing major normal faults. These faults, which originated as a result of regional tectonics, controlled pre- and post-caldera volcanism on the island. The formation of the caldera generated a very active geothermal system inside its depression, which is responsible for most of the present-day seismic activity and may also have a significant influence on the observed surface deformation. Our results do not support the hypothesis that there is a large but shallow, active magma chamber beneath the current caldera; instead we suggest that recent eruptions have been fed by small batches of deeper-sourced magmas. The intrusive remains of these eruptions and probably of other minor intrusions that have not reached the surface provide the main heat source that sustains the current geothermal system.

Impact of volcanism on the evolution of Lake Van I: evolution of explosive volcanism of Nemrut Volcano (eastern Anatolia) during the past >400,000 years

Abstract: The historically active Nemrut Volcano (2,948 m asl) (Eastern Anatolia), rising close to the western shore of huge alkaline Lake Van, has been the source of intense Plinian eruptions for >530,000 years (drilled lake sediments). About 40 widespread, newly recognized trachytic and less common rhyolitic fallout tephras and ca. 12 interbedded ignimbrites, sourced in Nemrut Volcano, are documented in stratigraphic traverses throughout an area of >6,000 km2 mostly west of Lake Van. Phenocrysts in the moderately peralkaline trachytes and rarer large-volume comenditic rhyolites comprise anorthoclase, hedenbergite-augite, fayalite and, especially in trachytic units, augite, minor aenigmatite, apatite and quartz, and rare chevkinite and zircon. Dacitic to rhyolitic tephras from nearby calcalkalic Suphan Volcano (4,058 m asl), locally interbedded with Nemrut tephras, are characterized by disequilibrium phenocryst assemblages (biotite, augitic clinopyroxene and hypersthene, minor olivine, common crystal clots and/or, in some deposits, amphibole). The magma volume (DRE) of the largest Nemrut tephra sheet (AP-1) described in detail may exceed 30 km3. Extreme facies and systematic compositional changes are documented in the ca. 30 ka Nemrut Formation (NF) deposits formed from one large and complex eruption (thick rhyolitic fallout overlain by ignimbrite, welded agglutinate, overbank surge deposits, and final more mafic fallout deposits). Common evidence of magma mixing in Nemrut ignimbrites reflects eruption from compositionally zoned magma reservoirs. Several young Cekmece Formation trachytes overlying ca. 30 ka old NF deposits and the late trachytes of the NF deposits show compositional affinities to tephra from Suphan Volcano possibly due to temporary influx of Suphan magmas into the Nemrut system following the evacuation of >10 km3 magma (DRE) during the caldera-forming NF eruption. Axes of large fallout fans are dominantly SW–NE but W–E in the younger sheets resembling the direction of the present dominant wind field. Growth of Nemrut volcanic edifice and its peripheral domes since before 0.5 Ma in the hinge area between the Van and Mus tectonic basins is likely to have been the major factor in isolating Lake Van basin thus initiating the origin and subsequent alkaline evolution of the lake. This alkalinity was later significantly controlled by climate forcing. Internal forcing mechanisms (volcanic and geodynamic) may also have contributed to major lake level changes in addition to climate forcing.

Biogeochemical processes involving dissolved CO 2 and CH 4 at Albano, Averno, and Monticchio meromictic volcanic lakes (Central-Southern Italy)

Abstract: This paper focuses on the chemical and isotopic features of dissolved gases (CH4 and CO2) from four meromictic lakes hosted in volcanic systems of Central–Southern Italy: Lake Albano (Alban Hills), Lake Averno (Phlegrean Fields), and Monticchio Grande and Piccolo lakes (Mt. Vulture). Deep waters in these lakes are characterized by the presence of a significant reservoir of extra-atmospheric dissolved gases mainly consisting of CH4 and CO2. The δ13C-CH4 and δD-CH4 values of dissolved gas samples from the maximum depths of the investigated lakes (from −66.8 to −55.6 ‰ V-PDB and from −279 to −195 ‰ V-SMOW, respectively) suggest that CH4 is mainly produced by microbial activity. The δ13C-CO2 values of Lake Grande, Lake Piccolo, and Lake Albano (ranging from −5.8 to −0.4 ‰ V-PDB) indicate a significant CO2 contribution from sublacustrine vents originating from (1) mantle degassing and (2) thermometamorphic reactions involving limestone, i.e., the same CO2 source feeding the regional thermal and cold CO2-rich fluid emissions. In contrast, the relatively low δ13C-CO2 values (from −13.4 to −8.2 ‰ V-PDB) of Lake Averno indicate a prevalent organic CO2. Chemical and isotopic compositions of dissolved CO2 and CH4 at different depths are mainly depending on (1) CO2 inputs from external sources (hydrothermal and/or anthropogenic); (2) CO2–CH4 isotopic exchange; and (3) methanogenic and methanotrophic activity. In the epilimnion, vertical water mixing, free oxygen availability, and photosynthesis cause the dramatic decrease of both CO2 and CH4 concentrations. In the hypolimnion, where the δ13C-CO2 values progressively increase with depth and the δ13C-CH4 values show an opposite trend, biogenic CO2 production from CH4 using different electron donor species, such as sulfate, tend to counteract the methanogenesis process whose efficiency achieves its climax at the water–bottom sediment interface. Theoretical values, calculated on the basis of δ13C-CO2 values, and measured δ13CTDIC values are not consistent, indicating that CO2 and the main carbon-bearing ion species (HCO3 −) are not in isotopic equilibrium, likely due to the fast kinetics of biochemical processes involving both CO2 and CH4. This study demonstrates that the vertical patterns of the CO2/CH4 ratio and of δ13C-CO2 and δ13C-CH4 are to be regarded as promising tools to detect perturbations, related to different causes, such as changes in the CO2 input from sublacustrine springs, that may affect aerobic and anaerobic layers of meromictic volcanic lakes.

An integrated field-numerical approach to assess slope stability hazards at volcanoes: the example of Pacaya, Guatemala

Abstract: Pacaya is an active stratovolcano located 30 km south of Guatemala City, Guatemala. A large (0.65 km3) sector collapse of the volcano occurred 0.6–1.6 ka B.P., producing a debris avalanche that traveled 25 km SW of the edifice. The current cone has since rebuilt within the scarp of this ancestral collapse. The structural setting of the volcano, along with two recent smaller-volume collapses in 1962 and 2010, suggests gravitational instability of this volcano. To assess Pacaya’s stability and potential for another large lateral collapse of the active cone, standard engineering methodologies for studying non-volcanic slopes were used to examine the SW flank of the edifice. A geomechanical model was developed based on the physical–mechanical material properties of Pacaya’s intact rocks and rock mass characteristics found through field observations and laboratory tests. Slope stability was analyzed in several scenarios with the Limit Equilibrium Method (LEM) and Finite Element Method (FEM), including static conditions (i.e., under gravity forces only), and considering the application of magma pressure and seismic force as triggering mechanisms for slope failure. Results show that the edifice remains stable under gravity alone; however, a large-scale collapse could be triggered by reasonable ranges of magma pressure (≥7.7 MPa if constant along a dyke) or peak ground acceleration (≥460 cm/s2). Results also suggest that a layer of pyroclastics beneath the edifice could have controlled the ancestral sector collapse. Structural analysis shows that a transtensional stress regime is causing a NW–SE orientation of aligned features at the surface, and may be a controlling mechanism for the direction of a future collapse. FEM results are concordant with those from LEM and reveal that maximum shear strain patterns within the edifice may account for long lava flows erupted from lower vent elevations.

The formation of columnar joints produced by cooling in basalt at Staffa, Scotland

Abstract: Columnar jointing in basaltic lava flows on the island of Staffa, NW Scotland, was studied using a combination of field mapping and measurement of column dimensions, sample petrology and measurements of plagioclase crystal size distributions (CSDs) interpreted using theoretical models of cooling. Four different lava flow units were measured, and column ordering was assessed using the hexagonality index and relative standard deviations of column side length, top area and internal angle. Upper and lower colonnades consist of dominantly 5-, 6- and 7-sided columns, with a hexagonality index value very similar to that of Giant’s Causeway and other basaltic columnar jointed localities. CSDs from samples at different heights within one colonnade were used to infer the propagation of the solidus isotherm, which was consistent with a convective cooling mechanism within the colonnade interior. Sample petrology and CSD measurements suggest that entablature can form both by the interaction of propagating joint sets and flooding of the flow surface by water, and the most widely exposed unit on Staffa shows evidence of both mechanisms operating on the same flow. Crystal size distribution measurements can provide a useful tool for field interpretation of lava flow cooling mechanisms.