Geology ReportsSearch

USGS · 70065872

Long-range hazard assessment of volcanic ash dispersal for a Plinian eruptive scenario at Popocatépetl volcano (Mexico): implications for civil aviation safety

Abstract

Popocatépetl is one of Mexico’s most active volcanoes threatening a densely populated area that includes Mexico City with more than 20 million inhabitants. The destructive potential of this volcano is demonstrated by its Late Pleistocene–Holocene eruptive activity, which has been characterized by recurrent Plinian eruptions of large magnitude, the last two of which destroyed human settlements in pre-Hispanic times. Popocatépetl’s reawakening in 1994 produced a crisis that culminated with the evacuation of two villages on the northeastern flank of the volcano. Shortly after, a monitoring system and a civil protection contingency plan based on a hazard zone map were implemented. The current volcanic hazards map considers the potential occurrence of different volcanic phenomena, including pyroclastic density currents and lahars. However, no quantitative assessment of the tephra hazard, especially related to atmospheric dispersal, has been performed. The presence of airborne volcanic ash at low and jet-cruise atmospheric levels compromises the safety of aircraft operations and forces re-routing of aircraft to prevent encounters with volcanic ash clouds. Given the high number of important airports in the surroundings of Popocatépetl volcano and considering the potential threat posed to civil aviation in Mexico and adjacent regions in case of a Plinian eruption, a hazard assessment for tephra dispersal is required. In this work, we present the first probabilistic tephra dispersal hazard assessment for Popocatépetl volcano. We compute probabilistic hazard maps for critical thresholds of airborne ash concentrations at different flight levels, corresponding to the situation defined in Europe during 2010, and still under discussion. Tephra dispersal mode is performed using the FALL3D numerical model. Probabilistic hazard maps are built for a Plinian eruptive scenario defined on the basis of geological field data for the “Ochre Pumice” Plinian eruption (4965 14 C yr BP). FALL3D model input eruptive parameters are constrained through an inversion method carried out with the semi-analytical HAZMAP model and are varied by sampling them using probability density functions. We analyze the influence of seasonal variations on ash dispersal and estimate the average persistence of critical ash concentrations at relevant locations and airports. This study assesses the impact that a Plinian eruption similar to the Ochre Pumice eruption would have on the main airports of Mexico and adjacent areas. The hazard maps presented here can support long-term planning that would help minimize the impacts of such an eruption on civil aviation.

Explore related subjects

90° N90° S · 180° W ← longitude → 180° E
Source-reported bounding extent: 15° to 30° latitude; -120° to -80° longitude. This indicates report coverage, not an exact sampling location. View area on OpenStreetMap.

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

Rosanna Bonasia, Chirara Scaini, Lucia Capra, Manuel Nathenson, Claus Siebe, Lilia Arana-Salinas, Arnau Folch. 2013-12-15. Long-range hazard assessment of volcanic ash dispersal for a Plinian eruptive scenario at Popocatépetl volcano (Mexico): implications for civil aviation safety. https://doi.org/10.1007/s00445-013-0789-z

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

Examples of eruption response teams from the Alaska Volcano Observatory

During times of eruption response, volcano observatories need to organize themselves differently than during normal operations. The number of formal operational roles grows to ensure that critical responsibilities are covered, including management of all activities at the observatory as well as increased staffing to ensure proper surveillance of data and issuance of timely notices and warnings. The scope and approach differ for each observatory and, in fact, for different eruptions. The Alaska Volcano Observatory (AVO) maintains an extensive monitoring program and issues forecasts and warnings about volcanic eruptions and unrest in Alaska. Since 2000, AVO has used formal roles to perform these duties and has implemented a variety of team approaches to respond to larger eruptions. For the Augustine (2006) and Redoubt (2009) eruptions, management scaled from 3 up to ~ 8 people in a command team to cover all aspects of the response. During these and other eruptions (e.g., Okmok, 2008), an operations room was staffed continuously to cover real-time responsibilities pertaining to monitoring and issuing alerts while the command team focused on overall management. More recently, such as for Bogoslof in 2016–2017 and Shishaldin in 2019 and 2023, AVO used a virtual real-time response team to handle warning tasks and variations of an Observatory Volcanic Event Response Team (OVERT; Moran et al. 2024) to manage overall observatory response activities. In 2025, AVO employed a formal OVERT for the first time to oversee its response to unrest at Mount Spurr. Frequent implementation, nimble scaling, constant evaluation, flexibility, and good communication make the team approach effective. We present examples of several response teams used over the last 25 years, and lessons learned from them, in the hope that these will be helpful to other observatories facing crisis responses. These examples may also allow stakeholders and the public to better understand how observatories work.

Alaska

Timescales of cumulate mobilization and mixing for the 1868 A.D. eruption of Mauna Loa, Island of Hawai‘i

The deadly 1868 A.D. eruption of Mauna Loa’s lower Southwest Rift Zone (Island of Hawai‘i) included a M7.9 earthquake and associated tsunami and landslides, demonstrating the severe hazards posed by Earth’s largest active subaerial volcano. To better understand the relationship between intense seismic activity, dike emplacement, magma storage, transport histories, and mobilization of olivine cumulates at Mauna Loa, we examine compositional zoning of olivine in the 1868 lava flows. Samples range from basalt (< 10% olivine) to picrite (30–40% olivine). The olivine cargo is heterogeneous (Fo 78.2–89.2 ; forsterite = [Mg/(Mg + Fe) × 100]) but dominated by ~ Fo 89 cores that lie above the Fe-Mg equilibrium field of host glasses. Crystal rims < Fo 80 are due to post-eruptive modification in slow cooling lava flows. Minor element compositions fall within the range of other Mauna Loa olivine erupted in the past 200 years. Olivine crystals exhibit both normal and complex Fo zoning patterns that yield timescales of diffusive re-equilibration that range from 3 to 258 days, with 72% of crystals recording 71 days or less. These timescales correspond to magmatic priming of the summit reservoir system ~ 2 months prior to the eruption and the M7.9 earthquake likely facilitated the transport of the crystal-rich summit-derived magmas downrift shortly prior to eruption. If the recently proposed faster Fe-Mg diffusion coefficient is used, timescales instead range from < 1 day to 25 days, with most recording 1 week or less. In this scenario, most of the olivine zoning would have to have been generated after the M7.9 earthquake perturbed the system.

Hawaii

Geochemistry of the 2022 Mauna Loa eruption: A comparison with earlier historical summit reservoir eruptions, with implications for magma supply and recharge

On November 28th, 2022, following a record historical repose period of 38 years, Mauna Loa erupted about 145 × 10 6 m 3 of lava and tephra over a 15-day period. The eruption was confined to the summit caldera region and the upper Northeast Rift Zone and is remarkably homogeneous in composition in both time and space. In these respects, it is typical of prior shallow summit reservoir magma bodies, recently estimated to be at a depth of around 1–2 km beneath the caldera. In contrast with these earlier magma bodies, which typically contain 6.7–7.1% MgO and are perched at the low-MgO end of olivine-control trends, the 2022 lava and tephra are more evolved with 6.24 + / − 0.03% MgO. This implies a temperature difference of around 11 °C with the prior 1984 magma. The simplest explanation is that over 38 years, cooling and crystallization of the remaining 1984 magma body has significantly exceeded magma recharge, giving rise to the evolved 2022 magma. The problem with this model is that we know from a variety of geophysical observations that in those 38 years, Mauna Loa has been erratically inflating, with heightened periods since around 2000 attributed to magma recharge. To reconcile these differences, we suggest instead that the 1984 magma cooled and crystallized much more extensively, from 1166 °C to around 1106 °C, co-crystallizing plagioclase, clinopyroxene, pigeonite and subsequently enstatite instead of pigeonite. At this point, the residual 1984 magma would have an MgO content around 4.2% and been about 50% solidified. Subsequent recharge and mixing by dominantly reservoir magmas, derived from a deeper 3–4-km intermediate magma reservoir, eventually produced the 2022 magma.

Hawaii