Messi a confronto i modelli di colonna eruttiva sviluppati dalla comunità vulcanologica internazionale per migliorare le stime del rischio per l’aviazione e dei pericoli associati alla ricaduta di lapilli e ceneri. I risultati dello studio, coordinato dall’INGV e dall’Università di Tokyo, sono stati pubblicati su un volume speciale di Journal of Volcanology and Geothermal Research
Sebbene la prima descrizione di una colonna eruttiva risalga a Plinio il giovane nella lettera a Tacito dove descriveva l'eruzione del Vesuvio del 79 d.C., la comprensione accurata della dinamica delle colonne eruttive, prodotte durante le eruzioni vulcaniche esplosive, rappresenta ancora una delle più importanti sfide della vulcanologia moderna, anche per le implicazioni che ha sul rischio per l’aviazione, legato alle emissioni di ceneri e, in generale, sui pericoli associati alla ricaduta di lapilli e ceneri.
A fare un confronto dei modelli di colonna eruttiva sviluppati dalle comunità vulcanologiche internazionali, uno studio coordinato dall’Istituto Nazionale di Geofisica e Vulcanologia (INGV) – sezione di Bologna e dall’Università di Tokyo. I risultati di questo lavoro di comparazione, promosso dalla Commissione di Tephra Hazard Modelling dell’International Association of Volcanology and Chemistry of Earth Interior (IAVCEI), sono stati pubblicati nel volume speciale 326 del Journal of Volcanology and Geothermal Research (
“Le dinamiche delle colonne eruttive”, spiega Antonio Costa, ricercatore dell’INGV e coordinatore dello studio, “sono influenzate significativamente dalle complesse interazioni tra di esse e l'atmosfera circostante. La comprensione di queste dinamiche è fondamentale per la stima del flusso di massa del materiale frammentato (tefra) rilasciato durante le eruzioni, quantità cruciale per i modelli di dispersione di tefra usati per la stima del rischio per l’aviazione e i pericoli associati alla ricaduta di lapilli e ceneri”.
Per queste ragioni negli ultimi decenni sono stati sviluppati diversi modelli di colonna eruttiva che vanno da semplici modelli empirici ai più recenti modelli basati sulla fluidodinamica computazionale.
“Lo studio confronta i risultati delle parametrizzazioni empiriche (modelli 0D), utilizzate per stimare il flusso di massa della colonna dalla stima della sua altezza, con quelli delle simulazioni di modelli numerici unidimensionali (1D) e tridimensionali (3D) sotto diverse condizioni vulcanologiche e meteorologiche, con l’obiettivo di valutare la robustezza dei modelli e gli aspetti che richiedono invece miglioramenti e ricerca futura. In particolare, lo studio ha coinvolto quattro modelli 0D, nove modelli 1D e quattro modelli 3D”, prosegue Costa.
“Uno dei risultati inaspettati dello studio è che le colonne eruttive “deboli” (ovvero caratterizzate da flussi di massa bassi) e “forti” (caratterizzate da flussi di massa alti), indipendentemente dalle condizioni di vento, sono governate da processi di inglobamento dell'aria molto diversi, che non sempre sono ben catturati dalle parametrizzazioni attualmente in uso, evidenziando quindi la necessità di ulteriori ricerche in questo campo”, continua il ricercatore.
Inoltre, il confronto con i modelli 3D ha mostrato chiaramente l'inadeguatezza delle attuali semplificazioni usate dai modelli 1D nel descrivere la dinamica delle colonne con flussi di massa più alti.
“Questi risultati hanno forti implicazioni pratiche in quanto diversi modelli confrontati nello studio sono utilizzati comunemente per stimare i flussi di massa durante eruzioni esplosive e la stima di questi flussi è un’informazione cruciale per i modelli di dispersione di cenere usati per le previsioni del trasporto delle nubi vulcaniche per la mitigazione del rischio legato all'effetto delle ceneri sugli aerei”, conclude Costa.
Abstract
This study compares and evaluates one-dimensional (1D) and three-dimensional (3D) numerical models of volcanic eruption columns in a set of different inter-comparison exercises. The exercises were designed as a blind test in which a set of common input parameters was given for two reference eruptions, representing a strong and a weak eruption column under different meteorological conditions. Comparing the results of the different models allows us to evaluate their capabilities and target areas for future improvement. Despite their different formulations, the 1D and 3D models provide reasonably consistent predictions of some of the key global descriptors of the volcanic plumes. Variability in plume height, estimated from the standard deviation of model predictions, is within ~20% for the weak plume and ~10% for the strong plume. Predictions of neutral buoyancy level are also in reasonably good agreement among the different models, with a standard deviation ranging from 9 to 19% (the latter for the weak plume in a windy atmosphere). Overall, these discrepancies are in the range of observational uncertainty of column height. However, there are important differences amongst models in terms of local properties along the plume axis, particularly for the strong plume. Our analysis suggests that the simplified treatment of entrainment in 1D models is adequate to resolve the general behaviour of the weak plume. However, it is inadequate to capture complex features of the strong plume, such as large vortices, partial column collapse, or gravitational fountaining that strongly enhance entrainment in the lower atmosphere. We conclude that there is a need to more accurately quantify entrainment rates, improve the representation of plume radius, and incorporate the effects of column instability in future versions of 1D volcanic plume models.
Fig. 1 da De Michieli Vitturi et al. (2016) JVGR, pag. 77-91, Vol. 326:
A. Aerial view showing Shinmoe-dake volcano peak erupting between Miyazaki and Kagoshima prefectures on January 27, 2011 (REUTERS/Kyodo)
B. The June 12, 1991 eruption column from Mount Pinatubo taken from the east side of Clark Air Base. U.S. Geological Survey Photograph taken by Richard P. Hoblitt
Fig. 2 Cover del volume speciale 326 del Journal of Volcanology and Geothermal Research