Sensor Technologies and Methods for Geoinformatics and Remote Sensing
View this Special IssueResearch Article  Open Access
Mapping Pyroclastic Flow Inundation Using Radar and Optical Satellite Images and Lahar Modeling
Abstract
Sinabung volcano, located above the Sumatra subduction of the IndoAustralian plate under the Eurasian plate, became active in 2010 after about 400 years of quiescence. We use ALOS/PALSAR interferometric synthetic aperture radar (InSAR) images to measure surface deformation from February 2007 to January 2011. We model the observed preeruption inflation and coeruption deflation using Mogi and prolate spheroid sources to infer volume changes of the magma chamber. We interpret that the inflation was due to magma accumulation in a shallow reservoir beneath Mount Sinabung and attribute the deflation due to magma withdrawal from the shallow reservoir during the eruption as well as thermoelastic compaction of erupted material. The pyroclastic flow extent during the eruption is then derived from the LAHARZ model based on the coeruption volume from InSAR modeling and compared to that derived from the Landsat 7 Enhanced Thematic Mapper Plus (ETM+) image. The pyroclastic flow inundation extents between the two different methods agree at about 86%, suggesting the capability of mapping pyroclastic flow inundation by combing radar and optical imagery as well as flow modeling.
1. Introduction
Volcanic pyroclastic flows consist of a large amount of destructive mass of very hot ash, lava fragments, and gases ejected explosively from a volcano and typically flow downslope at great speed. Volcanic pyroclastic flows pose great threats to people and property downstream from volcanoes. In order to accurately estimate the extents that are likely to be impacted by pyroclastic flows, computer models, such as LAHARZ [1], have been developed on the basis of statistical analyses of flowpath geometries and mathematical analyses of flow physics. LAHARZ [1] has been applied at many volcanoes worldwide.
Mount Sinabung is a stratovolcano in the northwestern part of the Sumatra Island, Indonesia (Figure 1). Formed during the Pleistocene to Holocene, the volcano consists of andesitic and dacitic lavas. The Sumatra Island, located on the boundary between the Australian Plate to the south and the Eurasian Plate to the north, was formed by the melting of subducted rocks under enormous pressure and its subsequent eruptions in volcanoes. Sumatra and Java, Indonesia, include a long line of volcanoes along the boundary between the two plates and have produced a large number of earthquakes and volcanic eruptions. Located within 50 km southeast of Sinabung, Sumatra Island’s Toba volcano had the largest eruption with a volcanic explosivity index (VEI) of 8 in the geological record during last 25 million years [2].
Mount Sinabung has a 2460 m high conical structure with four overlapping craters about 60 to 300 m in diameter from north to south, respectively, and is mostly covered with dense vegetation exclusively at the top of the volcano (Figure 2). Sinabung volcano had its last historical activity about 400 years ago, and the 2010 eruption started on August 29 and lasted to the middle of September 2010. After the first reawakened volcanic eruption in 2010, Mount Sinabung erupted again in 2013, 2014, and 2016.
SAR interferometry (InSAR) is an efficient method for measuring surface deformation over large areas [3–7], and it has been applied worldwide to study volcanoes [8–13]. Chaussard and Amelung [10] processed 13 Lband ALOS/PALSAR images acquired before 2010 and showed progressive inflation at a rate of 2.2 cm/yr using smallbaseline subset (SBAS) InSAR technique [14–16]. They also measured the surface displacement before the eruption and during the eruption and estimated the depth of magma storage (0.9 km) using volcanic source modeling (i.e., Mogi and Yang) [17–21]. However, these previous InSAR studies did not calculate the substantive damages of the volcanic activities but merely observed the volume change of the magma chamber according to the arithmetical deformation of the surface.
In this study, we have carried out a study of pyroclastic flow mapping using satellite imagery from both radar and optical sensors as well as the LAHARZ program (Figure 3). We improve the preliminary results [22–24] by laying out details on how to derive the pyroclastic flow extent using maximum likelihood classification method based on Landsat 7 images where artifacts due to scanline corrector have been rectified. We generate the coeruption InSAR image by stacking several interferograms in order to improve the signaltonoise ratio (SNR) of the coeruption deformation map, based on which we derive the extrusion volume associated with the 2010 eruption at Sinabung from geophysical models. The goal is to generate a highSNR deformation map for an accurate estimate of the eruption volume and improve the earlier results of timeseries deformation measurements where noise was apparent in some individual InSAR images [22–24]. We derive the volume change associated with Sinabung’s volcanic eruption in 2010 using mechanical modeling, and the calculated volume change of the magma chamber is then applied to a pyroclastic flow modeling algorithm based on LAHARZ [1, 25–27] to obtain the hazard zone of pyroclastic flow deposit. We describe in detail how we model the pyroclastic inundation map using the LAHARZ model and discuss quantitatively the comparison of pyroclastic flow simulation maps and pyroclastic flow map derived from the Landsat 7 ETM+ image of Sinabung volcano (Table 1). Finally, we highlight that the integration of satellite image analysis and pyroclastic flow simulation using the LAHARZ program can accomplish volcano monitoring and pyroclastic flow inundation hazard mapping of active volcanoes and discuss the limitation of the proposed method.
(a)  
 
(b)  

2. Processing of Optical Image
A pyroclastic flow inundation map can be generated by field surveys from the crater rim to the furthest extent of the pyroclastic flow after a volcanic eruption. However, a field survey in an active volcano can be dangerous due to the exposure to hazardous gases and abrupt activities. In contrast, remote sensing techniques are a useful tool for generating pyroclastic flow deposit maps, providing a safe, cost effective alternative to field mapping.
In this study, we aim to generate a pyroclastic flow deposit map using a Landsat 7 ETM+ image acquired on 30 Jul. 2012 after the 2010 Sinabung eruption (Table 1). This image represents the land cover after the 2010 eruptions with an absence of atmospheric effects (Figure 4(a)). However, the Landsat 7 ETM+ data suffer from an artifact due to the failure of the scanline corrector (SLC) in the ETM+ instrument. The Landsat 7 data can be restored by a gapfilling and filtering technique described by Lee et al. [28, 29] (Figure 4(b)). Finally, we use the gapfilled Landsat 7 ETM+ data (Figure 4(b)) to derive landcover types through supervised maximum likelihood classification (Figure 4(c)). We establish six classes, including farm land (orange), forest (green), city (gray), shadow (dark gray), water (blue), and pyroclastic flow deposit area (red). Based on the classification maps, the pyroclastic flow deposits from the 2010 eruption were mainly contained over the southern and eastern upper flanks of the volcano (red in Figure 4(c)). The pyroclastic flow inundation area after the 2010 eruption of Sinabung volcano is estimated as 0.872 km^{2} (red in Figure 4(d)).
(a)
(b)
(c)
(d)
3. Processing of Radar Images
InSAR processing is carried out to characterize the surface deformation of the Sinabung volcano. Using 20 ALOS/PALSAR images from 20 Feb. 2007 to 16 Jan. 2011 (Table 1), 40 interferograms are generated. The best coherence was constrained on the volcano summit. Outside the summit, coherence was maintained only at many scattered patches of rocky and lessvegetated areas.
The trend of surface deformation shows an uplift pattern until the end of August 2010 but a deflation pattern from the eruption event to the end of the dataset, corresponding to the changes in magma volume inside the reservoir. Therefore, we divide our results into two different periods of inflation and deflation separated by the onset of the 2010 eruption. The timeseries SBAS technique is applied to measure the surface deformation of the Sinabung volcano at each period [16]. Using this method, we can minimize artificial atmospheric effects [30, 31] and accurately measure surface deformation over a specific time period. Figure 5(a) displays the stacked image showing inflation of magma volume change until the 2010 eruption, and Figure 5(b) represents the stacked image that shows the deflation pattern of surface deformation during the 2010 eruption.
(a)
(b)
The two stacked images of surface deformation (Figures 6(a) and 6(b)) have opposite fringe tendency, with a maximum 1.2 cm/yr of preeruption inflation (Figure 6(c)) and 7 cm/yr of coeruption deflation (Figure 6(d)).
4. Volume Change Model of Magma Chamber
The Mogi model, based on a point source embedded in an elastic homogeneous halfspace [32], is applied to estimate the location and volume change of the mean surface deformation rate maps (Figures 6(a) and 6(b)) before and during the eruption. The bestfit Mogi source models of the observed interferograms (Figures 7(a) and 7(b)) are represented by synthetic interferograms in Figures 6(c) and 6(d). Residual interferograms (Figures 7(g) and 7(h)) can be calculated from the difference between the observed and synthetic interferograms. The depth of source is estimated to be ~0.3 km below sea level (BSL) before eruption and ~0.6 km BSL during the eruption (Table 2). This Mogi model estimates an increase in volume before the eruption as 1.9 × 10^{−6} km^{3}/yr (Figure 7). During the deflation of the coeruption period, the volume decrease is estimated to be −2.7 × 10^{−5} km^{3}/yr (Figure 7). As the Mogi model is essentially a point source, we further generate a bestfit model for using a prolate spheroid [8] for the preeruption and coeruption periods. The prolate spheroid source model (Figure 7(e)) is calculated to be at a depth of 1.3 km BSL beneath the center of the caldera before the eruption and ~1.1 km BSL during the eruption. The RMSE between the observations and model predictions from the prolate spheroid model is less than the corresponding one from the Mogi source (Table 2). Results from the prolate spheroid model suggest that the source of coeruption deflation is a steeply dipping elongate spheroid with major and minor axes of 2.7 km and 1.2 km, respectively, and a loss in pressure of −62.2 MPa/yr (Table 3). The preeruption source axes are estimated to be 3.0 km and 2.4 km with an increased pressure of 4.3 MPa/yr (Table 3).


5. LAHARZ Processing
The rationale for the delineation of laharinundation hazard zones [33, 34] is derived from the scaling analysis of lahar paths and statistical analysis using the LAHARZ program [35–37]. These analyses generate semiempirical equations that predict the inundated valley crosssectional area (1) and the planimetric area (2) as a function of lahar volume with proportionality coefficients c and C using volcanic event data of 27 paths of lahars at nine volcanoes. where is the crosssectional area, is the planimetric area, and are proportionality coefficients, and is the volume of the lahar. These calibrated predictive equations involve all the information needed to calculate and plot inundation on topographic maps through an empirical formula. Statistical methods from the LAHARZ program result in experience coefficients, and , and ((3) and (4)).
PF inundation maps are generated in this study using a statistically constrained simulation model to estimate potential hazard areas [25] by adapting methodology from the LAHARZ program [35]. The prediction equations for crosssectional (5) and planimetric (6) areas of PF hazards for Sinabung volcano are calibrated using data from dozens of events at several volcanoes. The equations are expressed with empirical proportionality coefficients ( to 0.1, to 40) resulting from the adaptation of LAHARZ program statistical analyses [35].
The PF inundation map for the 2010 Sinabung eruption is created in this study using the PF modeling algorithm from the modified LAHARZ program with empirical coefficient values of 0.1 for C and 40 for c. Several values for the volume of the pyroclastic flow in (5) and (6) are tested, including the volume change estimate calculated using the Mogi model with deformation observed in the coeruption interferogram (Table 2).
The inundation maps for the PF (Figure 8) estimated from the LAHARZ program can be compared to the results of the supervised classification created using the posteruptive Landsat 7 ETM+ image from 30 Jul. 2012. The modified LAHARZ program estimates are calculated using the four different values for the PF volume (red, brown, orange, and yellow colors).
(a)
(b)
Among the four different values for volume tested using the modified LAHARZ method, the volume of −2.7 × 10^{−5} km^{3} (red color) agrees relatively well with the result from the supervised classification delineation for the inundation limits of the PF deposits from the Landsat image. Region B (volume of PF simulation) in Table 4 is also included within approximately 84% of the PF inundation area delineated using the Landsat image. Region A represents a too small PF simulation area compared with the Landsat image result. Regions C and D, however, have common areas of only 78% and 67%, respectively, between the PF simulations and pyroclastic inundation extents based on the Landsat image (Figure 8).

6. Discussion
InSAR images over Sinabung indicate volcanic inflation before the 2010 eruption and deflation during the eruption. Modeling the observed deformation using Mogi and prolate spheroid sources has allowed us to estimate the magma source location, excessive pressure, the shape of the source, and volume changes during the preeruption and coeruption periods. The model results suggest there is a magma storage underneath Sinabung. The reservoir is located ~1 km BSL (~3 km beneath the peak of Sinabung). The magma reservoir is best characterized with a prolate spheroid with the vertical axis of ~2.7 km and horizontal axis of 1.2 km.
The location of the shallow magma reservoir matches the observations of seismicity recorded by the local seismic network. The Indonesia’s Centre for Volcanology and Geological Hazard Mitigation (CVGHM) installed four seismic stations (KWR, SKN, SKM, and MRD; green stars in Figure 9(a)) around the crater area of Sinabung after the first eruption on 29 Aug. 2010. These stations recorded earthquakes that occurred on September 6 and 7 (red and blue circles in Figure 9(a)). Figure 9(b) displays hypocenters and the depths of earthquakes that occurred before and after another eruption on 7 Sept. 2010 [38]. A magma reservoir at ~1 km BSL under the peak of Sinabung seems to be consistent with the fact that most earthquakes are located shallower than 1 km BSL (Figure 9). Such a magma reservoir is responsible for the observed inflation before the 2010 eruption and the deflation associated with the 2010 eruption.
(a)
(b)
Pyroclastic inundation hazard maps are generated in this study based on a modified LAHARZ program (PF simulation) using the volume change calculated by the Mogi model from the deformation results observed using stacked interferograms and then are compared with the inundation map from the supervised classification method with Landsat imagery. When we generate these pyroclastic inundation hazard maps, the volume change associated with the deflation of the magma chamber after the eruption was found to be the betterfit value (approximately 86% overlap) based on the comparison with the pyroclastic inundation maps from the Landsat images and the other PF simulations. This study has demonstrated that a volcano hazard inundation map can be generated using the PF model [25] by inputting parameters estimated from models of the magma chamber based on InSARderived surface deformation measurements. However, our method has some limitations. First, the eruption volume is calculated from the volume change by modeling the volcanic deflation from InSAR. Due to the vascularity, the eruption volume can be different from the deflation volume [8]. Hence, the volume used for LAHARZ modeling represents the lower bound of the volume of eruptive material. Second, we have ignored other eruptive products such as volcanic ash, which can lower our estimate for the volume used in LAHARZ modeling. Third, we have also neglected other thermodynamic effects on the magma chamber in the derivation of a deflation volcano from InSARderived deformation measurements [8].
7. Conclusions
In this study, we firstly measure the mean surface deformation using stacked interferograms by ALOS/PALSAR data before and after the 2010 eruption of Sinabung volcano on the Island of Sumatra, Indonesia. The mean surface deformation rate maps are used in the Mogi and spheroid models to estimate the depth, volume changes, and dimensions of the magma source before and during the 2010 eruption. The estimated depth of the magma source is around 1 km BSL; the volume change during the coeruption is calculated. The change in volume estimated from the Mogi model is used to generate a PF inundation hazard map using the modified LAHARZ program. We finally have verified the pyroclastic flow inundation area using the inundation map from the supervised classification method based on Landsat 7 ETM+ imagery. The bestfit result is when the volume for the PF simulation is based on the volume change from modeling the coeruption InSAR deformation map. We conclude that the combination of satellite image analysis and PF simulation using the modified LAHARZ program can accomplish volcano monitoring and pyroclastic flow inundation hazard mapping of active volcanoes.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This research was supported by the National Research Foundation of Korea (NRF) grant funded by the Korea government (MSIP) (no. 2017R1A2B4003258); the Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Science, ICT & Future Planning (2015M1A3A3A02013416); and the Korea Meteorological Administration Research and Development Program under Grant KMIPA (20153071). The SAR data used in this study are copyrighted ©JAXA/METI and were provided by Alaska Satellite Facility.
References
 S.P. Schilling, “LAHARZ: GIS programs for automated delineation of lahar hazard zones,” Tech. Rep., U.S. Geological Survey Openfile Report 98638, Reston, Virginia, 1998. View at: Google Scholar
 Global Volcanism Program, Smithsonian, http://volcano.si.edu/.
 D. Massonnet and K.L. Feigl, “Radar interferometry and its application to changes in the Earth’s surface,” Reviews of Geophysics, vol. 36, no. 4, pp. 441–500, 1998. View at: Publisher Site  Google Scholar
 E. Rignot, “Changes in West Antarctic ice stream dynamics observed with ALOS PALSAR data,” Geophysical Research Letters, vol. 35, no. 12, 2008. View at: Publisher Site  Google Scholar
 A.H.M. Ng, H.C. Chang, L. Ge, C. Rizos, and M. Omura, “Assessment of radar interferometry performance for ground subsidence monitoring due to underground mining,” Earth, Planets and Space, vol. 61, no. 6, pp. 733–745, 2009. View at: Publisher Site  Google Scholar
 Z. Lu, J.W. Kim, H. Lee et al., “Helmand River hydrologic studies using ALOS PALSAR InSAR and ENVISAT altimetry,” Marine Geodesy, vol. 32, no. 3, pp. 320–333, 2009. View at: Publisher Site  Google Scholar
 H.S. Jung, Z. Lu, and L. Zhang, “Feasibility of alongtrack displacement measurement from Sentinel1 interferometric wideswath mode,” IEEE Transactions on Geoscience and Remote Sensing, vol. 51, no. 1, pp. 573–578, 2013. View at: Publisher Site  Google Scholar
 M.E. Pritchard and M. Simons, “Surveying volcanic arcs with satellite radar interferometry: the central Andes, Kamchatka, and beyond,” GSA Today, vol. 14, no. 8, pp. 4–10, 2004. View at: Publisher Site  Google Scholar
 Z. Lu and D. Dzurisin, InSAR Imaging of Aleutian Volcanoes: Monitoring a Volcanic Arc from Space, SpringerVerlag, Berlin Heidelberg, 2014. View at: Publisher Site
 E. Chaussard and F. Amelung, “Monitoring the ups and downs of Sumatra and Java with DInsar timeseries,” in American Geophysical Union, Fall Meeting, San Francisco, CA, USA, 2010, abstract G23C0838. View at: Google Scholar
 S. K. Ebmeier, J. Biggs, T. A. Mather, G. Wadge, and F. Amelung, “Steady downslope movement on the western flank of Arenal volcano, Costa Rica,” Geochemistry, Geophysics, Geosystems, vol. 11, no. 12, 2010. View at: Publisher Site  Google Scholar
 B. Philibosian and M. Simons, “A survey of volcanic deformation on Java using ALOS PALSAR interferometric time series,” Geochemistry, Geophysics, Geosystems, vol. 12, no. 11, 2011. View at: Publisher Site  Google Scholar
 R. Furuta, “Case study of small scale surface deformation detection using ALOS PALSAR differential interferometry,” in Proceedings of FRINGE 2009 Workshop, Frascati, Italy, 2009. View at: Google Scholar
 E. Chaussard and F. Amelung, “Precursory inflation of shallow magma reservoirs at west Sunda volcanoes detected by InSAR,” Geophysical Research Letters, vol. 39, no. 21, 2012. View at: Publisher Site  Google Scholar
 C.W. Lee, Z. Lu, H.S. Jung, J.S. Won, and D. Dzurisin, “Surface deformation of Augustine Volcano, 19922005, from multipleinterferogram processing using a refined Small Baseline Subset (SBAS) Interferometric Synthetic Aperture Radar (InSAR) approach,” in The 2006 Eruption of Augustine Volcano, Alaska, pp. 453–465, U.S. Geological Survey Professional Paper 1769, Alaska, USA, 2011. View at: Google Scholar
 C.W. Lee, Z. Lu, and H.S. Jung, “Simulation of timeseries surface deformation to validate a multiinterferogram InSAR processing technique,” International Journal of Remote Sensing, vol. 33, no. 22, pp. 7075–7087, 2012. View at: Publisher Site  Google Scholar
 Z. Lu, T. Masterlark, D. Dzurisin, R. Rykhus, and C. Wicks Jr., “Magma supply dynamics at Westdahl volcano, Alaska, modeled from satellite radar interferometry,” Journal of Geophysical Research: Solid Earth, vol. 108, no. B7, 2003. View at: Publisher Site  Google Scholar
 Z. Lu, T. Masterlark, and D. Dzurisin, “Interferometric synthetic aperture radar study of Okmok volcano, Alaska, 1992–2003: magma supply dynamics and postemplacement lava flow deformation,” Journal of Geophysical Research: Solid Earth, vol. 110, no. B2, 2005. View at: Publisher Site  Google Scholar
 T. Masterlark, “Magma intrusion and deformation predictions: sensitivities to the Mogi assumptions,” Journal of Geophysical Research: Solid Earth, vol. 112, no. B6, 2007. View at: Publisher Site  Google Scholar
 H. S. Jung, Z. Lu, J. S. Won, M. P. Poland, and A. Miklius, “Mapping threedimensional surface deformation by combining multipleaperture interferometry and conventional interferometry: application to the June 2007 eruption of Kilauea volcano, Hawaii,” IEEE Geoscience and Remote Sensing Letters, vol. 8, no. 1, pp. 34–38, 2011. View at: Publisher Site  Google Scholar
 C.W. Lee, “Baekdusan volcano timeseries analysis from 1992 to 1998 using multiinterferogram InSAR processing,” Terrestrial, Atmospheric and Oceanic Sciences, vol. 25, no. 6, p. 743, 2014. View at: Publisher Site  Google Scholar
 C. W. Lee, Z. Lu, J. W. Kim, and S. K. Lee, “Volcanic activity analysis of Mt. Sinabung in Indonesia using InSAR and GIS techniques,” in 2015 IEEE International Geoscience and Remote Sensing Symposium (IGARSS), pp. 4793–4796, Milan, Italy, 2015. View at: Publisher Site  Google Scholar
 C. W. Lee, S. K. Lee, Z. Lu, and J. W. Kim, “Volcanic activity analysis of Mt. Sinabung in Indonesia using remote sensing and GIS techniques,” in NearSurface Asia Pacific Conference, pp. 150–153, Waikoloa, HI, USA, 2015. View at: Publisher Site  Google Scholar
 C. W. Lee, Z. Lu, and J. W. Kim, “Monitoring Mount Sinabung in Indonesia using multitemporal InSAR,” Korean Journal of Remote Sensing, vol. 33, no. 1, pp. 37–46, 2017. View at: Publisher Site  Google Scholar
 C. Widiwijayanti, B. Voight, D. Hidayat, and S.P. Schilling, “Objective rapid delineation of areas at risk from blockandash pyroclastic flows and surges,” Bulletin of Volcanology, vol. 71, no. 6, pp. 687–703, 2009. View at: Publisher Site  Google Scholar
 B.E. Hubbard, M.F. Sheridan, G. CarrascoNunez, R. DiazCastellon, and S.R. Rodriguez, “Comparative lahar hazard mapping at Volcan Citlaltépetl, Mexico using SRTM, ASTER and DTED1 digital topographic data,” Journal of Volcanology and Geothermal Research, vol. 160, no. 12, pp. 99–124, 2007. View at: Publisher Site  Google Scholar
 R.M. Iverson, S.P. Schilling, and J.W. Vallance, “Objective delineation of laharinundation hazard zones,” Geological Society of America Bulletin, vol. 110, no. 8, pp. 972–984, 1998. View at: Publisher Site  Google Scholar
 C.W. Lee, M. Cho, and Y.J. Choi, “Satellite imagery correction method and system,” Korea Patent 101303730, 2013. View at: Google Scholar
 C.W. Lee, Z. Lu, J.S. Won, H.S. Jung, and D. Dzurisin, “Dynamic deformation of Seguam Island, Alaska, 1992–2008, from multiinterferogram InSAR processing,” Journal of Volcanology and Geothermal Research, vol. 260, pp. 43–51, 2013. View at: Publisher Site  Google Scholar
 G. Nico, R. Tome, J. Catalao, and P. M. A. Miranda, “On the use of the WRF model to mitigate tropospheric phase delay effects in SAR interferograms,” IEEE Transactions on Geoscience and Remote Sensing, vol. 49, no. 12, pp. 4970–4976, 2011. View at: Publisher Site  Google Scholar
 J. Catalao, G. Nico, R. Hanssen, and C. Catita, “Merging GPS and atmospherically corrected InSAR data to map 3D terrain displacement velocity,” IEEE Transactions on Geoscience and Remote Sensing, vol. 49, no. 6, pp. 2354–2360, 2011. View at: Publisher Site  Google Scholar
 K. Mogi, “Relations between the eruptions of various volcanoes and the deformations of the ground surfaces around them,” Bulletin of the Earthquake Research Institute, vol. 36, pp. 99–134, 1958. View at: Google Scholar
 R.J. Blong, “Volcanic hazards risk assessment,” in Monitoring and Mitigation of Volcanic Hazards, pp. 675–698, Springer, Berlin, Heidelberg, 1996. View at: Publisher Site  Google Scholar
 A.R. Darnell, J.C. Phillips, J. Barclay, R.A. Herd, A.A. Lovett, and P.D. Cole, “Developing a simplified geographical information system approach to dilute lahar modelling for rapid hazard assessment,” Bulletin of Volcanology, vol. 75, no. 4, 2013. View at: Publisher Site  Google Scholar
 E. Chaussard, F. Amelung, and Y. Aoki, “Characterization of open and closed volcanic systems in Indonesia and Mexico using InSAR time series,” Journal of Geophysical Research: Solid Earth, vol. 118, no. 8, pp. 3957–3969, 2013. View at: Publisher Site  Google Scholar
 C.S. Magirl, P.G. Griffiths, and R.H. Webb, “Analyzing debris flows with the statistically calibrated empirical model LAHARZ in southeastern Arizona, USA,” Geomorphology, vol. 119, no. 12, pp. 111–124, 2010. View at: Publisher Site  Google Scholar
 D. O. Dorta, G. Toyos, C. Oppenheimer, M. T. Pareschi, R. Sulpizio, and G. Zanchetta, “Empirical modelling of the May 1998 small debris flows in Sarno (Italy) using LAHARZ,” Natural Hazards, vol. 40, no. 2, pp. 381–396, 2007. View at: Publisher Site  Google Scholar
 M. Iguchi, K. Ishihara, Surono, and M. Hendrasto, “Learn from 2010 Eruptions at Merapi and Sinabung Volcanoes in Indonesia,” Disaster Prevention Research Institute Annuals, vol. 54, no. B, pp. 185–194, 2011. View at: Google Scholar
Copyright
Copyright © 2018 ChangWook Lee et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.