Explosive Volcanic Eruptions Research Articles

Chronometry and Speedometry of Magmatic Processes using Chemical Diffusion in Olivine, Plagioclase and Pyroxenes

The magmatic processes that fuel volcanism, crustal growth, ore formation and discharge of volcanic gases and aerosols to the atmosphere occur across a range of timescales, from millions of years to just a few seconds. For example, the production of new oceanic crust at mid-ocean ridges is a near-continuous process that can operate in any one ocean basin on timescales of more than 100 m.y. However, the driving force for such processes is the spreading of the ocean plates that happens on a cm/yr timescale. At the other end of the spectrum, explosive volcanic eruptions involve the ascent and fragmentation of magma at velocities of the order of 100 m/s such that the journey from a magma chamber to an ash cloud may take place in a matter of minutes. In this case the driving force is the rapid expansion of magmatic gas in response to changes in pressure. At intermediate timescales magmatic processes may give rise to hydrothermal ore deposits on timescales of less than a million years for an individual deposit, while growth of giant granite batholiths may require piecemeal assembly of magma batches on timescales of a few million years. Although each of these processes has a characteristic, time-averaged timescale on which it operates, this is typically the end result of one or more natural processes that operate on much shorter timescales. For example, mid-ocean ridges do not extrude magma continuously onto the ocean floor, mineralising fluids do not discharge continuously through the shallow crust, and granitic magmas do not dribble continuously into evolving batholithic chambers. In some cases it is the long-term timescales that are important, for example the spreading rate of ocean basins, in others it is the short-term timescales that are important, for example the episodic growth of lava domes at active volcanoes. Although …

Paleo-environmental record of polycyclic aromatic hydrocarbons and polychlorobiphenyls at the peripheral site GV7 in Victoria Land (East Antarctica)

In this paper we investigated the presence of Polycyclic Aromatic Hydrocarbons and Polychlorobiphenyls in a 50-m deep snow/firn core collected at the peripheral site GV7 in East Antarctica during the 2013–2014 XXIX Italian expedition. The concentration depth profile was obtained on the basis of the total concentration of fourteen PAHs and seven PCBs individually determined by gas chromatography triple quadrupole mass spectrometry. Both classes of pollutants showed synchronized concentration vs time profile throughout the whole period of time covered by the snow/firn core (1892–2012). A correlation between major explosive volcanic eruptions and the concentration maxima of the pollutants was found. PAH maximum (9 ng/L) was about twice the background level (5 ng/L). PCBs showed a similar but more limited trend with barely visible volcanic maxima. This concurrence highlights the contribution of the major explosive volcanic events to the global contamination level for PAHs, as expected, but also for PCBs whose industrial production and use began in 1930. Excluding the maximum values, PAHs and PCBs showed an increase in the period 1956–1986: PCBs from about 0.05 to 0.21 ng/l (400% increase), and PAHs from about 3.5 to 7.8 ng/l (100% increase). Finally, in the last decade (2000–2010) the trend of these pollutants was different: (i) PCBs constantly decrease (from 0.15 ng/L to 0.10 ng/L), thanks to the implemented restriction on their production and on their use only in closed systems in many countries; (ii) PAHs remains practically constant around 6.5 ng/L.

Estimating building vulnerability to volcanic ash fall for insurance and other purposes

Volcanic ash falls are one of the most widespread and frequent volcanic hazards, and are produced by all explosive volcanic eruptions. Ash falls are arguably the most disruptive volcanic hazard because of their ability to affect large areas and to impact a wide range of assets, even at relatively small thicknesses. From an insurance perspective, the most valuable insured assets are buildings. Ash fall vulnerability curves or functions, which relate the magnitude of ash fall to likely damage, are the most developed for buildings, although there have been important recent advances for agriculture and infrastructure. In this paper, we focus on existing vulnerability functions developed for volcanic ash fall impact on buildings, and apply them to a hypothetical building portfolio impacted by a modern-day Tambora 1815 eruption scenario. We compare and contrast the different developed functions and discuss some of the issues surrounding estimation of potential building damage following a volcanic eruption. We found substantial variability in the different vulnerability estimates, which contribute to large uncertainties when estimating potential building damage and loss. Given the lack of detailed and published studies of building damage resulting from ash fall this is not surprising, although it also appears to be the case for other natural hazards for which there are far more empirical damage data. Notwithstanding the potential limitations of some empirical data in constraining vulnerability functions, efforts are required to improve our estimates of building damage under ash fall loading through the collection of damage data, experimental testing and perhaps theoretical failure analysis. For insurance purposes, the current building typologies provided for use with vulnerability functions are too detailed to map to the relatively limited information on building types that is typically available to insurers. Thus, efforts to provide vulnerability functions that can be used where only limited information is available regarding building types would also be valuable, both for insurers and for at-risk areas that have not been subject to detailed building vulnerability surveys.

Ancient Turks and Climate Disasters in Central Asia (534–679 AD)

This paper is the first to analyse the effect climate disasters (AD 536, 581, 626, 679) had on the historical processes in the east of Central Asia in the Early Middle Ages. Climate disasters are abrupt, severe, and sometimes protracted periods of cooling and drought that produce a negative effect on the economy of a nomadic society, the reason for this being large explosive volcanic eruptions. Cooling and drought brought famine to the nomads and caused illnesses. In addition to that, all the natural catastrophes were accompanied by plague pandemics among the eastern Turks. The Turkic Empire can be split into several chronological periods characterised by significant events that changed the course of history of the nomadic state, namely, AD 534–545 — the rise of the Turkic Empire; AD 581–583 — the division of the Turkic Empire into the Western and the Eastern Empires; AD 627–630 — the Eastern Turks are conquered by China; AD 679–687 — the second rise of the Eastern Turkic Empire. The research carried out shows that there is a certain connection between the important historical events and climate disasters in the history of the Turkic Empire (AD 534–745), this factor making it possible for the author to come forth with the proposition that the climatic factor did have an effect on the historical processes in the east of Central Asia, especially on the territories that had a nomadic economy.

Estimating the frequency of volcanic ash clouds over northern Europe

Fine ash produced during explosive volcanic eruptions can be dispersed over a vast area, where it poses a threat to aviation, human health and infrastructure. Here, we focus on northern Europe, which lies in the principal transport direction for volcanic ash from Iceland, one of the most active volcanic regions in the world. We interrogate existing and newly produced geological and written records of past ash fallout over northern Europe in the last 1000 years and estimate the mean return (repose) interval of a volcanic ash cloud over the region to be 44±7 years. We compare tephra records from mainland northern Europe, Great Britain, Ireland and the Faroe Islands, with records of proximal Icelandic volcanism and suggest that an Icelandic eruption with a Volcanic Explosivity Index rating (VEI) ≥ 4 and a silicic magma composition presents the greatest risk of producing volcanic ash that can reach northern Europe. None of the ash clouds in the European record which have a known source eruption are linked to a source eruption with VEI < 4. Our results suggest that ash clouds are more common over northern Europe than previously proposed and indicate the continued threat of ash deposition across northern Europe from eruptions of both Icelandic and North American volcanoes.

Regional variability in the frequency and magnitude of large explosive volcanic eruptions

Quantifying the frequency at which volcanic eruptions of different size occurs is impor- tant for hazard assessment. Volcanic records can be used to estimate the recurrence rate of large-magnitude eruptions (magnitude ≥4), but recording biases that impact data complete- ness complicate analysis. To overcome these biases, we conceptualize the volcanic record as a series of individual and unique time series associated by a common behavior. Thus, we approach issues of completeness on a volcano-by-volcano basis and use a hierarchical Bayes- ian approach to characterize the common frequency-magnitude (f-M) behavior for different groups of volcanoes. We identify variations in the f-M relationship between different volcano types and between different volcanic arcs. By accounting for systematic under-recording in the volcanic record, we also calculate the global recurrence rates for large-magnitude erup- tions during the Holocene, which are similar to previous estimates. However, higher recur- rence rates for smaller-magnitude events are observed, which is a result of our adjustments for data completeness. Quantifying how the f-M relationship varies between different groups of volcanoes provides an opportunity to understand how the tectonic setting in uences f-M behavior, which is important to quantify long-term regional volcanic hazard.

Pumice attrition in an air-jet

We present the results from a series of jet-attrition experiments performed using a standard ASTM device (ASTM D5757-00) on naturally occurring ash-sized (<2mm) pumice, a product of explosive volcanic eruption comprising highly porous silicate glass. We investigate the effect of both feed grain size and attrition duration on the production of fines. We utilize a wet methodology for fines collection to ensure recovery of the total grain size distribution for each experimental run. The experiments convert a restricted size range of pumice particles to a bimodal population of parent and daughter particles. The bimodal distribution develops even after short (~15min) attrition times. With increased attrition time, the volume of daughter particles increases and the mode migrates to finer grain sizes. Jet attrition efficiency depends heavily on the particle size of the feed; our data show little attrition for a feed of 500μm vs. highly efficient attrition for a 250μm feed. Our rates of attrition for pumice are extremely high compared to rates recovered from experiments on limestone pellets. Fines production data are well modeled by:mfinesmbed0=0.2911−e−0.312twhere m0bed is the initial mass of particles in the bed, t is in hours, and the two adjustable coefficients dictate the long time limiting behaviour (0.291) and the rate at which the limit is reached (−0.312). This functional form provides more realistic limits in time while preserving a zero intercept and defining a plateau for long residence times.

Sequential development of tidal ravinement surfaces in macro- to hypertidal estuaries with high volcaniclastic input: the Miocene Puerto Madryn Formation (Patagonia, Argentina)

The late Miocene beds of the Puerto Madryn Formation (Provincia del Chubut, Argentina) are formed by shallow marine and estuarine sediments. The latter include several tidal-channel infills well exposed on the cliffy coast of the Peninsula Valdes. The Bahia Punta Fosil and Cerro Olazabal paleochannels are end members of these tidal channels and show a fining-upward infilling starting with intraformational channel lag conglomerates above deeply erosional surfaces interpreted as fluvial ravinement surfaces (the erosion surface formed in the purely fluvial or the fluvially dominated part of the estuary, where erosion is driven by fluvial processes). These are overlain and eventually truncated (and suppressed) by the tidal ravinement surface (TRS), in turn covered with high-energy, bioclastic conglomerates mostly formed in the “tidally dominated/fluvially influenced” part of an estuary. Above, large straight or arcuate point bars with alternatively sandy/muddy seasonal beds and varying trace and body fossil contents were deposited from the freshwater fluvially dominated to saline-water tidally dominated part of the estuary. The upper channel infill is formed by cross-bedded sands with mud drapes and seaward-directed paleocurrents, together with barren, volcaniclastic sandy to muddy heterolithic seasonal rhythmites, both deposited in the fluvially dominated part of the estuary. Volcanic ash driven by the rivers after large explosive volcanic eruptions on land resulted in sedimentation rates as high as 0.9 m per year, preserving (through burial) the morphology of tidal channels and TRSs. The channel deposits were formed in a tide-dominated, macrotidal to hypertidal open estuary with well-developed TRSs resulting from strong tidal currents deeply scouring into the transgressive filling of the channels and eventually cutting the fluvial ravinement surface. The TRSs extended upstream to the inner part of the estuary during long periods of low sedimentation rates, extended channel migration and sediment bypass, interrupted by transient, high volcaniclastic input. The tidal channels of the Puerto Madryn Formation constitute a unique example of estuary sedimentation with pulsed sediment supply in a macrotidal to hypertidal estuary.

Merging field mapping and numerical simulation to interpret the lithofacies variations from unsteady pyroclastic density currents on uneven terrain: The case of La Fossa di Vulcano (Aeolian Islands, Italy)

In order to obtain results from computer simulations of explosive volcanic eruptions, one either needs a statistical approach to test a wide range of initial and boundary conditions, or needs using a well-constrained field case study via stratigraphy. Here we followed the second approach, using data obtained from field mapping of the Grotta dei Palizzi 2 pyroclastic deposits (Vulcano Island, Italy) as input for numerical modeling. This case study deals with impulsive phreatomagmatic explosions of La Fossa Cone that generated ash-rich pyroclastic density currents, interacting with the topographic high of the La Fossa Caldera rim. One of the simplifications in dealing with well-sorted ash (one particle size in the model) is to highlight the topographic effects on the same pyroclastic material in an unsteady current. We demonstrate that by merging field data with 3D numerical simulation results it is possible to see key details of the dynamical current-terrain interaction, and to interpret the lithofacies variations of the associated deposits as a function of topography-induced sedimentation (settling) rate. Results suggest that a value of the sedimentation rate lower than 5kg/m2s at the bed load can still be sheared by the overlying current, producing tractional structures (laminae) in the deposits. Instead, a sedimentation rate higher than that threshold can preclude the formation of tractional structures, producing thicker massive deposits. We think that the approach used in this study could be applied to other case studies (both for active and ancient volcanoes) to confirm or refine such threshold value of the sedimentation rate, which is to be considered as an upper value as for the limitations of the numerical model.

Rapid laccolith intrusion driven by explosive volcanic eruption.

Magmatic intrusions and volcanic eruptions are intimately related phenomena. Shallow magma intrusion builds subsurface reservoirs that are drained by volcanic eruptions. Thus, the long-held view is that intrusions must precede and feed eruptions. Here we show that explosive eruptions can also cause magma intrusion. We provide an account of a rapidly emplaced laccolith during the 2011 rhyolite eruption of Cordón Caulle, Chile. Remote sensing indicates that an intrusion began after eruption onset and caused severe (>200 m) uplift over 1 month. Digital terrain models resolve a laccolith-shaped body ∼0.8 km3. Deformation and conduit flow models indicate laccolith depths of only ∼20–200 m and overpressures (∼1–10 MPa) that likely stemmed from conduit blockage. Our results show that explosive eruptions may rapidly force significant quantities of magma in the crust to build laccoliths. These iconic intrusions can thus be interpreted as eruptive features that pose unique and previously unrecognized volcanic hazards.

Impact of Stratospheric Volcanic Aerosols on Age-of-Air and Transport of Long-Lived Species

The radiative perturbation associated to stratospheric aerosols from major explosive volcanic eruptions may induce significant changes in stratospheric dynamics. The aerosol heating rates warm up the lower stratosphere and cause a westerly wind anomaly, with additional tropical upwelling. Large scale transport of stratospheric trace species may be perturbed as a consequence of this intensified Brewer–Dobson circulation. The radiatively forced changes of the stratospheric circulation during the first two years after the eruption of Mt. Pinatubo (June 1991) may help explain the observed trend decline of long-lived greenhouse gases at surface stations (approximately −8 and −0.4 ppbv/year for CH4 and N2O, respectively). This decline is partly driven by the increased mid- to high-latitude downward flux at the tropopause and also by an increased isolation of the tropical pipe in the vertical layer near the tropopause, with reduced horizontal eddy mixing. Results from a climate-chemistry coupled model are shown for both long-lived trace species and the stratospheric age-of-air. The latter results to be younger by approximately 0.5 year at 30 hPa for 3–4 years after the June 1991 Pinatubo eruption, as a result of the volcanic aerosols radiative perturbation and is consistent with independent estimates based on long time series of in situ profile measurements of SF6 and CO2. Younger age of air is also calculated after Agung, El Chichón and Ruiz eruptions, as well as negative anomalies of the N2O growth rate at the extratropical tropopause layer. This type of analysis is made comparing the results of two ensembles of model simulations (1960–2005), one including stratospheric volcanic aerosols and their radiative interactions and a reference case where the volcanic aerosols do not interact with solar and planetary radiation.

The first tephra evidence for a Late Glacial explosive volcanic eruption in the Arxan-Chaihe volcanic field (ACVF), northeast China

Abstract A 5 mm thick tephra layer has been identified in the lacustrine sediments of Moon Lake in the Arxan-Chaihe volcanic field (ACVF) in Greater Khingan Mountains (NE China). The visible tephra layer is clearly revealed as a distinct peak in magnetic susceptibility measurements. The tephra layer consists mainly of brown vesicular glass shards and minor amounts of plagioclase, olivine and clinopyroxene. Major and minor element analysis has been carried out on the glass shards and plagioclase minerals. Glass shards show low concentrations of K 2 O, similar to the eruptive products derived from post-Miocene volcanoes of the ACVF. The plagioclase phenocrysts in both lava and tephra from ACVF, and in the tephra recorded in Moon Lake are labradorites. During the Late Pleistocene to Holocene, there were also extensive explosive eruptions in the nearby Nuominhe volcanic field (NVF). Volcanic rocks from the ACVF are easily distinguished from those derived from the NVF, having distinctly different K 2 O concentrations. This compositional variation is likely the result of different magmatic processes operating in the ACVF and NVF. Radiocarbon dating on organic materials from the lacustrine sediments dates the tephra layer to ca. 14,200 cal yrs BP, which implies that it was generated by a previously unknown Late Pleistocene explosive eruption in the ACVF. These results, for the first time, give a direct tephra record in this area, and suggest that identification of further tephra and/or cryptotephra in local sedimentary basins such as crater lakes of scoria cones and maars will be significant for dating the Late Pleistocene to Holocene volcanic eruptions and will help to establish a detailed record of the volcanic activity in the ACVF. The newly discovered tephra layer also provides a dated tephrochronological marker layer, which will in future studies provide a means to synchronise local sedimentary records of the climatically variable Late Glacial.

3‐D high‐speed imaging of volcanic bomb trajectory in basaltic explosive eruptions

AbstractImaging, in general, and high speed imaging in particular are important emerging tools for the study of explosive volcanic eruptions. However, traditional 2‐D video observations cannot measure volcanic ejecta motion toward and away from the camera, strongly hindering our capability to fully determine crucial hazard‐related parameters such as explosion directionality and pyroclasts' absolute velocity. In this paper, we use up to three synchronized high‐speed cameras to reconstruct pyroclasts trajectories in three dimensions. Classical stereographic techniques are adapted to overcome the difficult observation conditions of active volcanic vents, including the large number of overlapping pyroclasts which may change shape in flight, variable lighting and clouding conditions, and lack of direct access to the target. In particular, we use a laser rangefinder to measure the geometry of the filming setup and manually track pyroclasts on the videos. This method reduces uncertainties to 10° in azimuth and dip angle of the pyroclasts, and down to 20% in the absolute velocity estimation. We demonstrate the potential of this approach by three examples: the development of an explosion at Stromboli, a bubble burst at Halema'uma'u lava lake, and an in‐flight collision between two bombs at Stromboli.

On the mechanism of explosive eruption of mount erebus volcano: the dynamics of the rupture structure in a cavitating layer

This paper presents the results of an experimental simulation of rupture development in heavily cavitating magma melt flow in volcanic conduits and its effect on the structure of explosive volcanic eruptions. The dynamics of the state of a layer of distilled water (similar in the density of cavitation nuclei to magma melt) under shock-wave loading was studied. The experiments were performed using electromagnetic hydrodynamic shock tubes (EM HST) with maximum capacitor bank energy of up to 100 J and 5 kJ. It was found that the topology of the rupture formed on the membrane surface did not change during its development. Empirical estimates were obtained for the proportion of the capacitor bank energy expended in the development of the rupture and the characteristic time of its existence. The study revealed a number of fundamentally new physical effects in the cavity dynamics in a cavitating medium: a cavitation “boundary layer” is formed on the surface of the quasi-empty rupture, which is transformed into a cluster of high energy density upon closure of the flow.

Inferring the Effects of Compositional Boundary Layers on Crystal Nucleation, Growth Textures, and Mineral Chemistry in Natural Volcanic Tephras through Submicron-Resolution Imaging

Crystal nucleation and growth are first order processes captured in volcanic rocks and record important information about the rates of magmatic processes and chemical evolution of magmas during their ascent and eruption. We have studied glass-rich andesitic tephras from the Central Plateau of the Southern Taupo Volcanic Zone by electron- and ion-microbeam imaging techniques to investigate down to sub-micrometre scale the potential effects of compositional boundary layers (CBLs) of melt around crystals on the nucleation and growth of mineral phases and the chemistry of crystal growth zones. We find that CBLs may influence the types of mineral phases nucleating and growing, and growth textures such as the development of swallowtails. The chemistry of the CBLs also has the capacity to trigger intermittent overgrowths of nanometre-scale bands of different phases in rapidly growing crystals, resulting in what we refer to as cryptic phase zoning. The existence of cryptic phase zoning has implications for the interpretation of microprobe compositional data, and the resulting inferences made on the conditions of magmatic evolution. Identification of cryptic phase zoning may in future lead to more accurate thermobarometric estimates and thus geospeedometric constraints. In future, a more quantitative characterization of CBL formation and its effects on crystal nucleation and growth may contribute to a better understanding of melt rheology and magma ascent processes at the onset of explosive volcanic eruptions, and will likely be of benefit to hazard mitigation efforts.

A new interpolation method to model thickness, isopachs, extent, and volume of tephra fall deposits

Tephra thickness distribution is the primary piece of information used to reconstruct the histories of past explosive volcanic eruptions. We present a method for modeling tephra thickness with less subjectivity than is the case with hand-drawn isopachs, the current, most frequently used method. The algorithm separates the thickness of a tephra fall deposit into a trend and local variations and models them separately using segmented linear regression and ordinary kriging. The distance to the source vent and downwind distance are used to characterize the trend model. The algorithm is applied to thickness datasets for the Fogo Member A and North Mono Bed 1 tephras. Simulations on subsets of data and cross-validation are implemented to test the effectiveness of the algorithm in the construction of the trend model and the model of local variations. The results indicate that model isopach maps and volume estimations are consistent with previous studies and point to some inconsistencies in hand-drawn maps and their interpretation. The most striking feature noticed in hand-drawn mapping is a lack of adherence to the data in drawing isopachs locally. Since the model assumes a stable wind field, divergences from the predicted decrease in thickness with distance are readily noticed. Hence, wind direction, although weak in the case of Fogo A, was not unidirectional during deposition. A combination of the isopach algorithm with a new, data transformation can be used to estimate the extent of fall deposits. A limitation of the algorithm is that one must estimate “by hand” the wind direction based on the thickness data.

Localized field-aligned currents and 4-min TEC and ground magnetic oscillations during the 2015 eruption of Chile’s Calbuco volcano

The Calbuco volcano in southern Chile erupted on April 22, 2015. About 2 h after the first eruption, a Swarm satellite passed above the volcano and observed enhancement of small-amplitude (~0.5 nT) magnetic fluctuations with wave-packet structure which extends 15° in latitude. Similar wave packet is seen at the geomagnetic conjugate point of the volcano. Just after the eruption, geomagnetic fluctuations with the spectral peaks around the vertical acoustic resonance periods, 215 and 260 s, were also observed at Huancayo Geomagnetic Observatory located on the magnetic equator. Besides these observations, around 4-min, i.e., 175, 205 and 260 s, oscillations of total electron content (TEC) were observed at global positioning system stations near the volcano. The horizontal propagation velocity and the spatial scale of the TEC oscillation are estimated to be 720 m/s and 1600 km, respectively. These observations strongly suggest that the atmospheric waves induced by explosive volcanic eruption generate TEC variation and electric currents. The Swarm observation may be explained as a manifestation of their magnetic effects observed in the topside ionosphere.

Magma storage in a strike-slip caldera.

Silicic calderas form during explosive volcanic eruptions when magma withdrawal triggers collapse along bounding faults. The nature of specific interactions between magmatism and tectonism in caldera-forming systems is, however, unclear. Regional stress patterns may control the location and geometry of magma reservoirs, which in turn may control the spatial and temporal development of faults. Here we provide new insight into strike-slip volcano-tectonic relations by analysing Bouguer gravity data from Ilopango caldera, El Salvador, which has a long history of catastrophic explosive eruptions. The observed low gravity beneath the caldera is aligned along the principal horizontal stress orientations of the El Salvador Fault Zone. Data inversion shows that the causative low-density structure extends to ca. 6 km depth, which we interpret as a shallow plumbing system comprising a fractured hydrothermal reservoir overlying a magmatic reservoir with vol% exsolved vapour. Fault-controlled localization of magma constrains potential vent locations for future eruptions.

The effects of dynamics on the triboelectrification of volcanic ash

AbstractLightning is often observed during explosive volcanic eruptions, and the charging processes associated with these displays have been attributed to several mechanisms. In this work we delineate a set of experiments designed to quantify silicate‐based triboelectric charging in the volcanic context. Using natural samples from three different volcanoes, we show that the rate of triboelectrification in a fluidized bed depends on the energy input into the granular system. Experiments are conducted employing nonintrusive electrostatic sensors, ensuring that all charge exchange arises solely from particle‐particle collisions. At higher fluidization energies, particles undergo more frequent and energetic collisions, facilitating the transfer of charge. This finding implies that triboelectric charging could help promote charging in regions of the eruptive system that contain numerous particle‐particle collisions such as the conduit and gas thrust regions. Our experiments also suggest that surface charge density is capped, at least in part, by atmospheric conditions, specifically the breakdown characteristics of the gas.

Stratospheric Aerosols from Major Volcanic Eruptions: A Composition-Climate Model Study of the Aerosol Cloud Dispersal and e-folding Time

Large explosive volcanic eruptions are capable of injecting considerable amounts of particles and sulfur gases above the tropopause, causing large increases in stratospheric aerosols. Five major volcanic eruptions after 1960 (i.e., Agung, St. Helens, El Chichón, Nevado del Ruiz and Pinatubo) have been considered in a numerical study conducted with a composition-climate coupled model including an aerosol microphysics code for aerosol formation and growth. Model results are compared between an ensemble of numerical simulations including volcanic aerosols and their radiative effects (VE) and a reference simulations ensemble (REF) with no radiative impact of the volcanic aerosols. Differences of VE-REF show enhanced diabatic heating rates; increased stratospheric temperatures and mean zonal westerly winds; increased planetary wave amplitude; and tropical upwelling. The impact on stratospheric upwelling is found to be larger when the volcanically perturbed stratospheric aerosol is confined to the tropics, as tends to be the case for eruptions which were followed by several months with easterly shear of the quasi-biennial oscillation (QBO), e.g., the Pinatubo case. Compared to an eruption followed by a period of westerly QBO, such easterly QBO eruptions are quite different, with meridional transport to mid- and high-latitudes occurring later, and at higher altitude, with a consequent decrease in cross-tropopause removal from the stratosphere, and therefore longer decay timescale. Comparing the model-calculated e-folding time of the volcanic aerosol mass during the first year after the eruptions, an increase is found from 8.1 and 10.3 months for El Chichón and Agung (QBO westerly shear), to 14.6 and 30.7 months for Pinatubo and Ruiz (QBO easterly shear). The corresponding e-folding time of the global-mean radiative flux changes goes from 9.1 and 8.0 months for El Chichón and Agung, to 28.7 and 24.5 months for Pinatubo and Ruiz.