El Campo Volcánico Monogenético Guamuez-Sibundoy, ubicado en el límite entre los departamentos de Nariño y Putumayo, está conformado por 19 centros de emisión identificados hasta la fecha, en cercanías a poblaciones que albergan en conjunto cerca de 27.000 personas. Con el objetivo de evaluar la probabilidad espacial de la formación de nuevos centros eruptivos en el sector y establecer escenarios potenciales de actividad eruptiva futura, articularmente en relación con flujos de lava, se utilizó una metodología probabilística de largo plazo. Esto involucró la recopilación de datos geológicos, cartográficos, topográficos y meteorológicos, integrados en una base de datos geoespacial desarrollada en ArcGIS 10.8©. Los datos más relevantes para la evaluación de la probabilidad espacial incluyeron fallas, centros de emisión, actividad sísmica y manifestaciones hidrotermales, los cuales fueron subclasificados en subconjuntos y se les asignaron valores de importancia relativa mediante el método de Estimación del Juicio de Expertos (EJE). Posteriormente, se obtuvieron las Funciones de Densidad de Probabilidad (PDFs) de cada conjunto de datos mediante el plugin QVAST, permitiendo obtener un primer mapa de susceptibilidad del campo volcánico. Este mapa fue actualizado teniendo en cuenta la reciente actividad sísmica ocurrida en agosto de 2024. La simulación de escenarios de amenaza para flujos de lava se llevó a cabo a partir de los puntos con mayor valor de probabilidad obtenidos en cada mapa, utilizando el plugin Q-LavHA, sobre un DEM de 12,5 m de resolución. Los parámetros físicos de los flujos considerados para las simulaciones incluyen espesor, longitud y factores como la topografía y el perfil del terreno. Los resultados indicaron que la zona con mayor probabilidad de apertura de nuevos centros de emisión se encuentra al este de la Laguna de La Cocha, en cercanía a la mayor concentración de volcanes. Con respecto a las simulaciones de flujos de lava realizadas en los dos puntos de ambos mapas, se observaron dos trayectorias hacia el norestesureste y hacia el noroeste respectivamente, cada uno abarcando un área máxima de 6,1 km2 y 3,8 km2. Estos resultados destacan la importancia de definir el área con mayor probabilidad, ya que a su vez, esto conllevará un impacto diferente en las zonas a partir de donde se realizan las simulaciones de escenarios eruptivos. Asimismo, estos resultados permiten pronosticar los posibles impactos de la actividad eruptiva en el área, proporcionando una base sólida para la gestión del riesgo y la planificación de medidas preventivas.