**Marcelo Cando-Jácome and Antonio Martínez-Graña \***

Geology Department, External Geodynamics Area, Faculty of Sciences, University of Salamanca, Plaza Merced s/n, 37008 Salamanca, Spain; id00709713@usal.es

**\*** Correspondence: amgranna@usal.es; Tel.: +34-923-294496; Fax: +34-923-923294514

Received: 2 February 2019; Accepted: 25 March 2019; Published: 26 March 2019

**Abstract:** On 3 June 2018, a strong eruption of the Fuego volcano in Guatemala produced a dense cloud of 10-km-high volcanic ash and destructive pyroclastic flows that caused nearly 200 deaths and huge economic losses in the region. Subsequently, due to heavy rains, destructive secondary lahars were produced, which were not plotted on the hazard maps using the LAHAR Z software. In this work we propose to complement the mapping of this type of lahars using remote-sensing (Differential Interferometry, DINSAR) in Sentinel images 1A and 2A, to locate areas of deformation of the relief on the flanks of the volcano, areas that are possibly origin of these lahars. To determine the trajectory of the lahars, parameters and morphological indices were analyzed with the software System for Automated Geoscientific Analysis (SAGA). The parameters and morphological indices used were the accumulation of flow (FCC), the topographic wetness index (TWI), the length-magnitude factor of the slope (LS). Finally, a slope stability analysis was performed using the Shallow Landslide Susceptibility software (SHALSTAB) based on the Mohr–Coulomb theory and its parameters: internal soil saturation degree and effective precipitation, parameters required to destabilize a hillside. In this case, the application of this complementary methodology provided a more accurate response of the areas destroyed by primary and secondary lahars in the vicinity of the volcano.

**Keywords:** volcano deformation; lahars hazard; magma accumulation; pyroclastic flows; ash plumes
