JDR Vol.17 No.5 pp. 805-817
doi: 10.20965/jdr.2022.p0805


Continuously Operable Simulator and Forecasting the Deposition of Volcanic Ash from Prolonged Eruptions at Sakurajima Volcano, Japan

Masato Iguchi*,†, Haruhisa Nakamichi*, Kosei Takishita*,**, and Alexandros P. Poulidis***

*Disaster Prevention Research Institute, Kyoto University
1722-19 Sakurajima-Yokoyama-cho, Kagoshima, Kagoshima 891-1419, Japan

Corresponding author

**Division of Earth and Planetary Sciences, Graduate School of Science, Kyoto University, Kyoto, Japan

***Institute of Environmental Physics, University of Bremen, Bremen, Germany

January 24, 2022
June 15, 2022
August 1, 2022
volcanic ash, meteorological-ash-dispersion, X-MP radar, disdrometer, Sakurajima volcano

At Sakurajima volcano, frequent Vulcanian eruptions have been seen at the summit crater of Minamidake since 1955. In addition to this eruption style, the eruptive activities of Strombolian type and prolonged ash emission also occur frequently. We studied the design of a simulator of advection-diffusion-fallout of volcanic ash emitted continuously. The time function of volcanic ash eruption rate is given by a linear combination of the volcanic tremor amplitude and the volume change of the pressure source obtained from the ground motion. The simulation is repeated using discretized values of the eruption rate time function at an iteration time interval of the simulation. The integrated value of the volcanic ash deposition on the ground obtained from each individual simulation is used to estimate the value of the ash fallout. The plume height is given by an empirical equation proportional to a quarter of the power of the eruption rate. Since the wind velocity field near the volcano is complicated by the influence of the volcanic topography, the predicted values published by meteorological organizations are made in high resolution by Weather Research and Forecasting (WRF) for the simulation. We confirmed that an individual simulation can be completed within a few minutes of iteration interval time, using the PUFF model as the Lagrangian method and FALL3D-8.0 as the Eulerian method on a general-purpose PC. Except during rainfall, the radar reflectivity, the count of particles per particle size, and fall velocity obtained by the disdrometers can be used for the quasi-real time matching of the plume height calculated from the eruption rate and the ash fall deposition rate obtained from the simulation.

Cite this article as:
M. Iguchi, H. Nakamichi, K. Takishita, and A. Poulidis, “Continuously Operable Simulator and Forecasting the Deposition of Volcanic Ash from Prolonged Eruptions at Sakurajima Volcano, Japan,” J. Disaster Res., Vol.17, No.5, pp. 805-817, 2022.
Data files:
  1. [1] T. Kobayashi, “Geology of Sakurajima Volcano: A review,” Bull. Volcanol. Soc. Japan, Vol.27, No.4, pp. 277-292, 1982 (in Japanese with English abstract).
  2. [2] J. Hickey, J. Gottsmann, H. Nakamichi, and M. Iguchi, “Thermomechanical controls on magma supply and volcanic deformation: Application to Aira caldera, Japan,” Sci. Rep., Vol.6, Article No.32691, doi: 10.1038/srep32691, 2016.
  3. [3] S. Scollo, A. Folch, and A. Costa, “A parametric and comparative study of different tephra fallout models,” J. Volcanol. Geotherm. Res., Vol.176, No.2, pp. 199-211, 2008.
  4. [4] C. Bonadonna et al., “Probabilistic modeling of tephra dispersal: Hazard assessment of a multiphase rhyolitic eruption at Tarawera, New Zealand,” J. Geophys. Res. Solid Earth, Vol.110, No.B3, Article No.B03203, 2005.
  5. [5] A. Folch and A. Felpeto, “A coupled model for dispersal of tephra during sustained explosive eruptions,” J. Volcanol. Geotherm. Res., Vol.145, Nos.3-4, pp. 337-349, 2005.
  6. [6] K. Mannen, “Total grain size distribution of a mafic subplinian tephra, TB-2, from the 1986 Izu-Oshima eruption, Japan: An estimation based on a theoretical model of tephra dispersal,” J. Volcanol. Geotherm. Res., Vol.155, Nos.1-2, pp. 1-17, 2006.
  7. [7] J. Fero, S. N. Carey, and J. T. Merrill, “Simulation of the 1980 eruption of Mount St. Helens using the ash-tracking model PUFF,” J. Volcanol. Geotherm. Res., Vol.175, No.3, pp. 355-366, 2008.
  8. [8] A. Folch, C. Cavazzoni, A. Costa, and G. Macedonio, “An automatic procedure to forecast tephra fallout,” J. Volcanol. Geotherm. Res., Vol.177, No.4, pp. 767-777, 2008.
  9. [9] M. Iguchi, “Method for real-time evaluation of discharge rate of volcanic ash – Case study on intermittent eruptions at the Sakurajima volcano, Japan –,” J. Disaster Res., Vol.11, No.1, pp. 4-14, 2016.
  10. [10] A. Folch, A. Costa, and G. Macedonio, “FALL3D: A computational model for transport and deposition of volcanic ash,” Comput. Geosci., Vol.35, No.6, pp. 1334-1342, 2009.
  11. [11] A. P. Poulidis, T. Takemi, and M. Iguchi, “Experimental high-resolution forecasting of volcanic ash hazard at Sakurajima, Japan,” J. Disaster Res., Vol.14, No.5, pp. 786-797, 2019.
  12. [12] H. L. Tanaka and K. Yamamoto, “Numerical simulation of volcanic plume dispersal from Usu volcano in Japan on 31 March 2000 using PUFF model,” Earth Planet. Space, Vol.54, No.7, pp. 743-752, 2002.
  13. [13] H. L. Tanaka and M. Iguchi, “Numerical simulations of volcanic ash plume dispersal for Sakura-jima using a real-time emission rate estimation,” J. Disaster Res., Vol.14, No.1, pp. 160-172, 2019.
  14. [14] A. P. Poulidis, I. A. Renfrew, and A. J. Matthews, “Thermally induced convective circulation and precipitation over an isolated volcano,” J. Atmos. Sci., Vol.73, No.4, pp. 1667-1686, 2016.
  15. [15] W. C. Skamarock et al., “A description of the advanced research WRF model version 4 (No.980 NCAR/TN-556+STR),” Technical Report, National Center for Atmospheric Research, doi: 10.5065/1dfh-6p97, 2019.
  16. [16] M. Maki et al., “Preliminary results of weather radar observations of Sakurajima volcanic smoke,” J. Disaster Res., Vol.11, No.1, pp. 15-30, 2016.
  17. [17] K. Takishita, A. P. Poulidis, and M. Iguchi, “In-situ measurement of tephra deposit load based on a disdrometer network at Sakurajima volcano, Japan,” J. Volcanol. Geotherm. Res., Vol.421, Article No.107442, 2022.
  18. [18] K. Mogi, “Relation between the eruptions of various volcanoes and the deformations of the ground surface around them,” Bull. Earthq. Res. Inst., Univ. Tokyo, Vol.36, pp. 99-134, 1958.
  19. [19] M. Iguchi et al., “Characteristics of volcanic activity at Sakurajima volcano’s Showa crater during the period 2006 to 2011,” Bull. Volcanol. Soc. Japan, Vol.58, No.1, pp. 115-135, 2013.
  20. [20] M. Iguchi et al., “Integrated monitoring of volcanic ash and forecasting at Sakurajima volcano, Japan,” J. Disaster Res., Vol.14, No.5, pp. 798-809, 2019.
  21. [21] L. G. Mastin et al., “A multidisciplinary effort to assign realistic source parameters to models of volcanic ash-cloud transport and dispersion during eruptions,” J. Volcanol. Geotherm. Res., Vol.186, Nos.1-2, pp. 10-21, 2009.
  22. [22] M. Iguchi, “Eruptive activity of Sakurajima volcano during the period from July 2011 to June 2012,” M. Iguchi (Ed.), “Report on Integrated volcano observation for the study on preparation process of eruption at Sakurajima volcano 2012,” pp. 53-61, Kyoto University, 2013, (in Japanese with English abstract) [accessed July 13, 2022]
  23. [23] A. P. Poulidis, T. Takemi, and M. Iguchi, “The effect of wind and atmospheric stability on the morphology of volcanic plumes from vulcanian eruptions,” J. Geophys. Res. Solid Earth. Vol.124, No.8, pp. 8013-8029, 2019.
  24. [24] A. T. Prata, L. Mingari, A. Folch, G. Macedonio, and A. Costa, “FALL3D-8.0: A computational model for atmospheric transport and deposition of particles, aerosols and radionuclides – Part 2: Model validation,” Geosci. Model Dev., Vol.14, No.1, pp. 409-436, 2021.
  25. [25] B. R. Morton, G. I. Taylor, and J. S. Turner, “Turbulent gravitational convection from maintained and instantaneous sources,” Proc. R. Soc. A, Vol.234, No.1196, pp. 1-23, 1956.
  26. [26] A. P. Poulidis, T. Takemi, M. Iguchi, and I. A. Renfrew, “Orographic effects on the transport and deposition of volcanic ash: A case study of Mount Sakurajima, Japan,” J. Geophys. Res. Atmos., Vol.122, No.17, pp. 9332-9350, 2017.
  27. [27] A. P. Poulidis, S. Biass, G. Bagheri, T. Takemi, and M. Iguchi, “Atmospheric vertical velocity – A crucial component in understanding proximal deposition of volcanic ash,” Earth Planet. Sci. Lett., Vol.566, Article No.116980, 2021.
  28. [28] A. P. Poulidis and T. Takemi, “A 1998–2013 climatology of Kyushu, Japan: Seasonal variations of stability and rainfall,” Int. J. Climatol., Vol.37, No.4, pp. 1843-1858, 2017.
  29. [29] M. C. Diaz Vecino et al., “Aerodynamic characteristics and genesis of aggregates at Sakurajima Volcano, Japan,” Sci. Rep., Vol.12, Article No.2044, 2022.
  30. [30] R. Gunn and G. D. Kinzer, “The terminal velocity of fall for water droplets in stagnant air,” J. Meteorol., Vol.6, No.4, pp. 243-248, 1949.
  31. [31] A. Folch et al., “FALL3D-8.0: A computational model for atmospheric transport and deposition of particles, aerosols and radionuclides – Part 1: Model physics and numerics,” Geosci. Model Dev., Vol.13, No.3, pp. 1431-1458, 2020.
  32. [32] M. Iguchi, “A vertical expansion source model for the mechanisms of earthquakes originated in the magma conduit of an andesitic volcano: Sakurajima, Japan,” Bull. Volcanol. Soc. Jpn., Vol.39, No.2, pp. 49-67, 1994.
  33. [33] M. Iguchi, T. Tameguri, J. Hirabayashi, and H. Nakamichi, “Forecasting volcanic eruption of Sakurajima volcano based on magma intrusion rate,” Bull. Volcanol. Soc. Jpn., Vol.64, No.2, pp. 33-51, 2019 (in Japanese with English abstract).

*This site is desgined based on HTML5 and CSS3 for modern browsers, e.g. Chrome, Firefox, Safari, Edge, Opera.

Last updated on Aug. 05, 2022