Granular flow rheology; the key to understanding the exceptional mobility of pyroclastic density currents
Lead Research Organisation:
University of Edinburgh
Department Name: Sch of Geosciences
Abstract
Pyroclastic density currents (PDCs) are hot avalanches of volcanic rock, pumice, ash and gas that descend the flanks of volcanoes. They can destroy and bury 100's km2 of terrain. Their high temperatures, inherent mobility and unpredictable nature render them one of the most hazardous volcanic phenomena. Since 1600AD, pyroclastic flows have resulted in over 90,000 deaths, 33% of all volcanic fatalities recorded, making them the single biggest cause of death at volcanoes.
Forecasting the flow paths and the extent of inundation by pyroclastic density currents at a given volcano depends on our understanding of (i) the flow mechanisms involved (ii) developing models that can faithfully capture the dynamic nature of those flows and accurately simulate past events, and (iii) applying those models probabilistically, so that all possible future scenarios at a given volcano can be considered in order to generate probabilistic hazard maps. Here we will tackle (i) and (ii), but our track history in (iii) demonstrates our longer-term intention. The rationale for this research therefore stems from both a strong end-user defined need, as well as motivation to advance the science of these complex multiphase (particle and gas) natural flows. The aim of this research is therefore to improve the capability of forecasting pyroclastic density current inundation zones around volcanoes by making breakthroughs in understanding the interplay between flow behaviour and how the rheological nature of the flow changes as it propagates.
During flow, pyroclastic density currents progressively develop regions that vary in their physical nature and flow mechanisms. Typically, the flows develop high particle concentrations at the base, with frictional or collisional contacts between the particles. An overriding ash cloud develops above this, where particle concentration is low and most particles are supported by turbulent convection of hot gases. As the flows propagate over topography, these upper and lower regions respond differently to changes in slope and valley confinement. Acceleration, deceleration and spreading of the upper and lower units occur at different points, and flow separation can be induced. The propensity for these upper ash clouds to separate from the parent basal flow and travel in unexpected directions often results in lethal consequences.
This research will focus on understanding the rheological variations in the basal granular flow and will consider how it may, in turn, modulate mass flux into the overriding ash cloud. We will test the hypothesis that variations in the basal undercurrent rheology, in part induced by topography, result in pore fluid pressure fluctuations that feed the generation and separation of upper turbulent ash clouds from their parent undercurrents. We will achieve this by integrating data obtained from complementary field, geomorphological, experimental and computational studies, in particular utilising cutting-edge modelling tools developed for engineering applications. We will build on important new advances in the understanding of industrial granular flows to characterise how flow rheology varies (through time and space), and what controls those variations.
Our results will form the basis for a new constitutive rheology description, providing a fundamental step forward by allowing advance from flow-averaged rheology laws currently employed in flow simulation tools used for hazard quantification. Extensions of this work, in particular the application of the new generation simulation tools will produce hazard maps that have lower associated uncertainties. Using methods we have already developed for probabilistic hazard mapping, we will quantify that degree of improvement. The project is timely and will benefit from synergy with a major Edinburgh-based initiative on industrial granular flows, as well as ongoing research by project partners.
Forecasting the flow paths and the extent of inundation by pyroclastic density currents at a given volcano depends on our understanding of (i) the flow mechanisms involved (ii) developing models that can faithfully capture the dynamic nature of those flows and accurately simulate past events, and (iii) applying those models probabilistically, so that all possible future scenarios at a given volcano can be considered in order to generate probabilistic hazard maps. Here we will tackle (i) and (ii), but our track history in (iii) demonstrates our longer-term intention. The rationale for this research therefore stems from both a strong end-user defined need, as well as motivation to advance the science of these complex multiphase (particle and gas) natural flows. The aim of this research is therefore to improve the capability of forecasting pyroclastic density current inundation zones around volcanoes by making breakthroughs in understanding the interplay between flow behaviour and how the rheological nature of the flow changes as it propagates.
During flow, pyroclastic density currents progressively develop regions that vary in their physical nature and flow mechanisms. Typically, the flows develop high particle concentrations at the base, with frictional or collisional contacts between the particles. An overriding ash cloud develops above this, where particle concentration is low and most particles are supported by turbulent convection of hot gases. As the flows propagate over topography, these upper and lower regions respond differently to changes in slope and valley confinement. Acceleration, deceleration and spreading of the upper and lower units occur at different points, and flow separation can be induced. The propensity for these upper ash clouds to separate from the parent basal flow and travel in unexpected directions often results in lethal consequences.
This research will focus on understanding the rheological variations in the basal granular flow and will consider how it may, in turn, modulate mass flux into the overriding ash cloud. We will test the hypothesis that variations in the basal undercurrent rheology, in part induced by topography, result in pore fluid pressure fluctuations that feed the generation and separation of upper turbulent ash clouds from their parent undercurrents. We will achieve this by integrating data obtained from complementary field, geomorphological, experimental and computational studies, in particular utilising cutting-edge modelling tools developed for engineering applications. We will build on important new advances in the understanding of industrial granular flows to characterise how flow rheology varies (through time and space), and what controls those variations.
Our results will form the basis for a new constitutive rheology description, providing a fundamental step forward by allowing advance from flow-averaged rheology laws currently employed in flow simulation tools used for hazard quantification. Extensions of this work, in particular the application of the new generation simulation tools will produce hazard maps that have lower associated uncertainties. Using methods we have already developed for probabilistic hazard mapping, we will quantify that degree of improvement. The project is timely and will benefit from synergy with a major Edinburgh-based initiative on industrial granular flows, as well as ongoing research by project partners.
Planned Impact
We identify 4 principal socio-economic beneficiaries of this research programme: volcanic hazard practitioners, policy makers, industry, and the general public in the UK.
1. Volcanic hazard practitioners: Volcanic hazard practitioners are individuals at government agencies around the world responsible for undertaking volcano monitoring and hazard assessments. This work will achieve impact through facilitating substantial improvements of the most widely used pyroclastic density current simulation tool used to undertake hazard assessments, TITAN2D. TITAN2D, is a shallow water depth averaged model, which simulates granular flows over natural terrains. We will work with project partners and TITAN2D developers at the University of Buffalo to translate our new variable rheology constitutive description into the TITAN2D code. Once a new rheology substitution for TITAN2D is formulated, we will work with collaborators at geological surveys and participants at Hazard Assessment workshops to implement and showcase the resultant new tool for hazard assessment.
2. Policy makers: Decision makers utilize hazard assessments and information about the state of unrest of a particular volcano to make hazard management decisions for short term hazard mitigation as well as long term planning. Through facilitating major improvements to robust volcanic hazard assessment, this project offers a pathway to reducing risk. As well as addressing NERCs increasing resilience to natural hazard impact area, this research also follows the UN International Strategy for Disaster Reduction guidelines for natural hazards and risk management by addressing the volcanic phenomena with the highest impact on lives, and potentially livelihoods and property assets. We will engage directly with UNISDR Science and Technical Advisory group and contribute to their case study documents. These activities will be facilitated by the BGS who have a world-leading track record in research impact.
3. Industry: This project involves interdisciplinary research, which can lead to innovative approaches and results affecting diverse academic fields as well as industry. Free-surface gravity currents are one of the most common particle transport mechanisms in the handling of industrial bulk solids (powders and granular materials), across diverse sectors such as the chemical, pharmaceutical, food and personal care industries. It is of mutual benefit to interact with industry and we plan to participate in the Annual Industry Club events to be organised through Co-I Sun's EPSRC project and host a joint NERC/EPSRC workshop in 2019 to promote synergy, inviting both academic and industrial partners and stakeholders.
4. General public in the UK: Increasing awareness of natural geophysical processes will help attract students into Earth science related degrees. In addition, Edinburgh is home to two excellent museums that cover geoscience topics, with which we plan to collaborate: the National Museum of Scotland (free entry, 1.9 million visitors/year), and Dynamic Earth (330K visitors/year). We will develop a mobile, interactive gravity current exhibit based on a recently constructed mini flume tank coupled with interactive computer-based activities and we will tour the exhibit around dedicated science outreach and open days throughout the duration of the project and beyond.
1. Volcanic hazard practitioners: Volcanic hazard practitioners are individuals at government agencies around the world responsible for undertaking volcano monitoring and hazard assessments. This work will achieve impact through facilitating substantial improvements of the most widely used pyroclastic density current simulation tool used to undertake hazard assessments, TITAN2D. TITAN2D, is a shallow water depth averaged model, which simulates granular flows over natural terrains. We will work with project partners and TITAN2D developers at the University of Buffalo to translate our new variable rheology constitutive description into the TITAN2D code. Once a new rheology substitution for TITAN2D is formulated, we will work with collaborators at geological surveys and participants at Hazard Assessment workshops to implement and showcase the resultant new tool for hazard assessment.
2. Policy makers: Decision makers utilize hazard assessments and information about the state of unrest of a particular volcano to make hazard management decisions for short term hazard mitigation as well as long term planning. Through facilitating major improvements to robust volcanic hazard assessment, this project offers a pathway to reducing risk. As well as addressing NERCs increasing resilience to natural hazard impact area, this research also follows the UN International Strategy for Disaster Reduction guidelines for natural hazards and risk management by addressing the volcanic phenomena with the highest impact on lives, and potentially livelihoods and property assets. We will engage directly with UNISDR Science and Technical Advisory group and contribute to their case study documents. These activities will be facilitated by the BGS who have a world-leading track record in research impact.
3. Industry: This project involves interdisciplinary research, which can lead to innovative approaches and results affecting diverse academic fields as well as industry. Free-surface gravity currents are one of the most common particle transport mechanisms in the handling of industrial bulk solids (powders and granular materials), across diverse sectors such as the chemical, pharmaceutical, food and personal care industries. It is of mutual benefit to interact with industry and we plan to participate in the Annual Industry Club events to be organised through Co-I Sun's EPSRC project and host a joint NERC/EPSRC workshop in 2019 to promote synergy, inviting both academic and industrial partners and stakeholders.
4. General public in the UK: Increasing awareness of natural geophysical processes will help attract students into Earth science related degrees. In addition, Edinburgh is home to two excellent museums that cover geoscience topics, with which we plan to collaborate: the National Museum of Scotland (free entry, 1.9 million visitors/year), and Dynamic Earth (330K visitors/year). We will develop a mobile, interactive gravity current exhibit based on a recently constructed mini flume tank coupled with interactive computer-based activities and we will tour the exhibit around dedicated science outreach and open days throughout the duration of the project and beyond.
Organisations
- University of Edinburgh (Lead Research Organisation)
- National Institute for Geophysics and Volcanology (INGV) (Collaboration)
- Resilience Development Initiative (Collaboration)
- University at Buffalo, State University of New York (Project Partner)
- United States Geological Survey (Project Partner)
- National Institute of Geophysics and Volcanology (Project Partner)
- Heriot-Watt University (Project Partner)
Publications
Barker T
(2021)
Exact solutions for steady granular flow in vertical chutes and pipes
in Journal of Fluid Mechanics
Langmuir, A.
(2019)
La historia se repite
Najuch T
(2019)
Analysis of two partially-saturated-cell methods for lattice Boltzmann simulation of granular suspension rheology
in Computers & Fluids
Zhu C
(2020)
Solid-like and liquid-like granular flows on inclined surfaces under vibration - Implications for earthquake-induced landslides
in Computers and Geotechnics
Zhu C
(2023)
A granular energy-controlled boundary condition for discrete element simulations of granular flows on erodible surfaces
in Computers and Geotechnics
Description | Multihazards in Indonesia |
Amount | £50,000 (GBP) |
Organisation | Natural Environment Research Council |
Sector | Public |
Country | United Kingdom |
Start | 04/2022 |
End | 03/2024 |
Description | Quantitative characterization of the rheology and transport-sedimentation mechanisms of debris flows by integrating deposits analyses, large-scale experiments and numerical modeling |
Amount | £98,476 (GBP) |
Funding ID | NAF\R2\180833 |
Organisation | The Royal Society |
Sector | Charity/Non Profit |
Country | United Kingdom |
Start | 12/2018 |
End | 12/2021 |
Description | University if Edinburgh - Global Impact Accelerator Account; "A new generation of Volcanic Hazard Maps that are useful, usable, and used." |
Amount | £33,315 (GBP) |
Organisation | University of Edinburgh |
Sector | Academic/University |
Country | United Kingdom |
Start | 11/2018 |
End | 03/2019 |
Title | Dataset of pyroclastic density current deposit characteristics |
Description | More than 200 publications were analysed to extract deposit, in particular grainsize information. The resultant dataset contains nearly 300 distributions associated with events at 20 volcanoes. |
Type Of Material | Database/Collection of data |
Year Produced | 2021 |
Provided To Others? | No |
Impact | Informed discussion on best practices on publishing grainsize information for use in further studies |
Title | IMEX_SfloW2D_v2 |
Description | The model simulates transient volcanic flows in 2D and can be applied to reproduce laboratory experiments and to simulate natural flows over topography. Model development led by colleagues at INGV Pisa, with contributions from this project. |
Type Of Material | Computer model/algorithm |
Year Produced | 2022 |
Provided To Others? | Yes |
Impact | Model applied in community cross-validation efforts |
URL | https://github.com/demichie/IMEX_SfloW2D_v2 |
Title | Physical properties and rheological characteristics of volcanic materials. |
Description | This database contains data collected during systematic powder rheometry laboratory measurements of granular material sampled from real pyroclastic density currents (PDC) deposits. By means of the compressibility, shear and wall friction tests performed using the FT4 powder rheometer (Freeman Technology), material bulk density, angle of internal friction (yield curve), wall friction angle and the flowability of the natural samples have been determined. These parameters have been evaluated for both bulk samples of PDC matrix material as well as specific grainsize sub-fractions. |
Type Of Material | Database/Collection of data |
Year Produced | 2019 |
Provided To Others? | Yes |
Impact | Preliminary results show that the parameters characterizing the analysed granular materials (e.g. the flowability) are strongly controlled by the material properties, and in particular grainsize and, more importantly, the content of fine ash. Additionally, the wall-friction angle is strongly controlled by the roughness of the wall, and in particular by the interplay between the roughness thickness and the characteristic grainsize of the material. |
Title | Updated dataset: Physical properties and rheological characteristics of volcanic materials |
Description | This is the updated version of the dataset first presented last year. The update includes new data. This database contains data collected during systematic powder rheometry laboratory measurements of granular material sampled from real pyroclastic density currents (PDC) deposits. By means of the compressibility, shear and wall friction tests performed using the FT4 powder rheometer (Freeman Technology), material bulk density, angle of internal friction (yield curve), wall friction angle and the flowability of the natural samples have been determined. These parameters have been evaluated for both bulk samples of PDC matrix material as well as specific grainsize sub-fractions. |
Type Of Material | Database/Collection of data |
Year Produced | 2020 |
Provided To Others? | Yes |
Impact | Preliminary results show that the parameters characterizing the analysed granular materials (e.g. the flowability) are strongly controlled by the material properties, and in particular grainsize and, more importantly, the content of fine ash. Additionally, the wall-friction angle is strongly controlled by the roughness of the wall, and in particular by the interplay between the roughness thickness and the characteristic grainsize of the material. We are now in the process of collecting new data of wall-friction angle with varying wall roughness |
Description | Co-organised a workshop on Understanding Geological Hazards to Support Disaster Risk Assessment in Indonesia |
Organisation | Resilience Development Initiative |
Country | Indonesia |
Sector | Academic/University |
PI Contribution | We organised a series of monthly meetings with RDI to plan an online workshop with several institutes from across Indonesia. The workshop covered several hazards, including landslide, volcanic, seismic and tsunami. We workshop resulted in the production of a collaborative report between BGS and RDI. |
Collaborator Contribution | RDI co-organised and hosted the online workshop, and co-wrote the related workshop report. |
Impact | A workshop report was produced |
Start Year | 2022 |
Description | Collaboration with INGV Pisa on numerical modelling of flows |
Organisation | National Institute for Geophysics and Volcanology (INGV) |
Country | Italy |
Sector | Public |
PI Contribution | Production of datasets and observations on pyroclastic density current deposits and their propagation, invitation to meetings and discussions, and co-development of numerical models |
Collaborator Contribution | Led production of open source numerical models of pyroclastic density currents |
Impact | Numerical model simulating pyroclastic density current propagation in 2D space |
Start Year | 2021 |
Description | Guest lecture for ITB international science seminar |
Form Of Engagement Activity | A talk or presentation |
Part Of Official Scheme? | No |
Geographic Reach | International |
Primary Audience | Postgraduate students |
Results and Impact | Gave an invited presentation to the International science seminar organised by ITB. Presentation attended by many universities across Indonesia through zoom, and was also broadcast on YouTube. Presentation focused on recent eruptive activity of Anak Krakatau, with the following discussion focusing on the physical process of explosive eruptions. |
Year(s) Of Engagement Activity | 2020 |
Description | Presentation of results at IAVCEI 2023 Conference Rotorua |
Form Of Engagement Activity | A talk or presentation |
Part Of Official Scheme? | No |
Geographic Reach | International |
Primary Audience | Professional Practitioners |
Results and Impact | Presented poster on dataset of physical parameters of pyroclastic density current (PDCs) deposits during IAVCEI 2023 meeting. This enabled discussion with others working on (PDCs), highly the dataset as a key new resource in the field. |
Year(s) Of Engagement Activity | 2023 |
URL | https://confer.eventsair.com/iavcei2023/programme |
Description | Presentation to NERC COMET Consortium |
Form Of Engagement Activity | A talk or presentation |
Part Of Official Scheme? | No |
Geographic Reach | National |
Primary Audience | Postgraduate students |
Results and Impact | Presented on recent research including that conducted on Anak Krakatau for this project for International Womens Day webinar |
Year(s) Of Engagement Activity | 2023 |
URL | https://youtu.be/rdqXuuslyMw |
Description | Workshop with Indonesian colleagues to identify areas for future work |
Form Of Engagement Activity | A formal working group, expert panel or dialogue |
Part Of Official Scheme? | No |
Geographic Reach | National |
Primary Audience | Policymakers/politicians |
Results and Impact | co-organised and hosted workshop with colleagues from RDI. Presentations on volcanic, landslide, earthquake and tsunami hazard, with discussion sessions for each topic to identify incountry needs, and potential for future collaborations |
Year(s) Of Engagement Activity | 2021 |