Probability maps of pyroclastic density currents (PDCs) by using a modified approach of the box 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.
Application to the San Salvador volcanic complex:
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.
Application to the Aso caldera:
Bevilacqua, A., Aravena, A., Aspinall, W., Costa, A., Mahony, S., Neri, A., Sparks, R., & Hill, B. (2022). Assessing minimum pyroclastic density current mass to impact critical infrastructures: example from Aso caldera (Japan). Natural Hazards and Earth System Sciences, 22(10), 3329-3348.
Usage normal version: python3 BoxMapProb.py
Input parameter: See the file input_data.py for the definition of the input parameters.
Usage user-friendly version: python3 Friendly_BoxMapProb.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/).