Probability maps of pyroclastic density currents (PDCs) by using a modified approach of the energy cone model and SRTM 30 m elevation data.
Citation of the paper describing the branching strategy adopted here:
Aravena, A., Cioni, R., Bevilacqua, A., de' Michieli Vitturi, M., Esposti Ongaro, T., and Neri, A. (2020). Tree‐branching based enhancement of kinetic energy models for reproducing channelization processes of pyroclastic density currents. Journal of Geophysical Research: Solid Earth, 125, e2019JB019271.
Citation of the paper describing calibration strategies and user manuals:
Aravena, A., Bevilacqua, A., de' Michieli Vitturi, M., Esposti Ongaro, T., Neri, A., and Cioni, R. (2022). Calibration strategies of PDC kinetic energy models and their application to the construction of hazard maps. Bulletin of Volcanology, 84, 29.
Citation of the paper using this model systematically to study the effect of topography of a set of 50 volcanoes:
Aravena, A., Roche, O (2022). Influence of the topography of stratovolcanoes on the propagation and channelization of dense pyroclastic density currents analyzed through numerical simulations. Bulletin of Volcanology, 84, 67.
Application to the San Salvador and Nejapa-Chiltepe volcanic complexes:
Aravena, A., Bevilacqua, A., Neri, A., Gabellini, P., Ferrés, D., Escobar, D., Aiuppa, A., & Cioni, R. (2023). Scenario-based probabilistic hazard assessment for explosive events at the San Salvador Volcanic complex, El Salvador. Journal of Volcanology and Geothermal Research, 107809.
Bevilacqua, A., Aravena, A., Neri, A., Gutierrez, E., Escobar, D., Schliz, M., Aiuppa, A., and Cioni, R. (2021). Thematic vent opening probability maps and hazard assessment of small-scale pyroclastic density currents in the San Salvador volcanic complex (El Salvador) and Nejapa-Chiltepe volcanic complex (Nicaragua). Natural Hazards and Earth System Sciences, 21, 1639-1665.
Usage normal version: python3 ECMapProb.py
Input parameter: See the file input_data.py for the definition of the input parameters.
Usage user-friendly version: python3 Friendly_ECMapProb.py
This code needs:
- Installation of matplotlib (https://matplotlib.org/), numpy (http://www.numpy.org/), elevation (https://pypi.org/project/elevation/), Pillow (https://pillow.readthedocs.io/), utm (https://pypi.org/project/utm/), and tifffile (https://pypi.org/project/tifffile/).