Super-Volcanoes and Dangerous Bubbles

Gas bubbles tend to accumulate in the upper layers of magma reservoirs, which are only a few kilometers beneath the earth’s surface, building up pressure that can then be abruptly liberated by eruption. These bubbles mainly contain water vapor but also sulfur.

Sulfur-rich eruptions

“Such volcanic eruptions can be extremely powerful and spew an enormous amount of ash and sulfur to the surface,” says Andrea Parmigiani, a post-doc in the Institute of Geochemistry and Petrology at ETH Zurich. “We’ve known for some time that gas bubbles play a major role in such events, but we had only been able to speculate on how they accumulate in magma reservoirs.”

The 1815 Mount Tambora eruption. The red areas are maps of the thickness of volcanic ashfall.

Credit: Wikimedia Commons

Together with other scientists from ETH Zurich and Georgia Institute of Technology (Georgia Tech), the researchers studied the behavior of bubbles with a computer model.

The scientists used theoretical calculations and laboratory experiments to examine in particular how bubbles in crystal-rich and crystal-poor layers of magma reservoirs move buoyantly upward. In many volcanic systems, the magma reservoir consists mainly of two zones: an upper layer consisting of viscous melt with almost no crystals, and a lower layer rich in crystals, but still containing pore space.

Super bubbles meander through a maze

When Andrea Parmigiani, Christian Huber and Olivier Bachmann started this project, they thought that the bubbles, as they moved upwards through crystal-rich areas of the magma reservoirs, would dramatically slow down, while they would go faster in the crystal-poor zones.

“Instead, we found that, under volatile-rich conditions, they would ascend much faster in the crystal-rich zones, and accumulate in the melt-rich portions above” says Parmigiani.

Parmigiani explains this as follows: when the proportion of bubbles in the pore space of the crystal-rich layers increases, small individual bubbles coalesce into finger-like channels, displacing the existing highly viscous melt. These finger-like channels allow for a higher vertical gas velocity. The bubbles, however, have to fill at least 10 to 15 % of the pore space.

“If the vapor phase cannot form these channels, individual bubbles are mechanically trapped,” says the earth scientist. As these finger-like channels reach the boundary of the crystal-poor melt, individual, more spherical bubbles detach, and continue their ascent towards the surface. However, the more bubble, the more reduce their migration velocity is.

Tambora volcano on Indonesia’s Sumbawa Island was the site of the world’s largest historical eruption in April 1815. This NASA Landsat mosaic shows the 6-km-wide caldera truncating the 2850-m-high summit of the massive volcano. Pyroclastic flows during the 1815 eruption reached the sea on all sides of the 60-km-wide volcanic peninsula, and the ejection of large amounts of tephra caused world-wide temperature declines in 1815 and 1816.

Credit: NASA

This is because each bubble creates a return flow of viscous melt around it. When an adjacent bubble feels this return flow, it is slowed down. This process was demonstrated in a laboratory experiment conducted by Parmigiani’s colleagues Salah Faroughi and Christian Huber at Georgia Tech, using water bubbles in a viscous silicone solution.

“Through this mechanism, a large number of gas bubbles can accumulate in the crystal-poor melt under the roof of the magma reservoir. This eventually leads to over-pressurization of the reservoir,” says lead author Parmigiani. And because the bubbles also contain sulfur, this also accumulates, explaining why such a volcano might emit more sulfur than expected based on its composition.

What this means for the explosivity of a given volcano is still unclear. “This study focuses primarily on understanding the basic principles of gas flow in magma reservoirs; a direct application to prediction of volcanic behavior remains a question for the future,” says the researcher, adding that existing computer models do not depict the entire magma reservoir, but only a tiny part of it: roughly a square of a few cubic centimeter with a clear boundary between the crystal-poor and crystal-rich layers.

To calculate this small volume, Parmigiani used high-performance computers such as the Euler Cluster at ETH Zurich and a supercomputer at the Swiss National Supercomputing Centre in Lugano.

For the software, the researcher had access to the open-source library Palabos, which he continues to develop in collaboration with researchers from University of Geneva. “This software is particularly suitable for this type of simulation,” says the physicist.

Contacts and sources: