Universidad de Nariño

Ecuaciones Diferenciales
Proyecto final de Ecuaciones Diferenciales

Presentado por:

Álvarez Álvarez Richard Leonel

Caicedo Muñoz Andrea Alejandra
Gelpud Delgado Yenifer Damaris

Presentado a:

Eduardo Ibargüen Mondragón

May 2025
1 Resumen
Este proyecto de investigación propone un modelo matemático simplificado basado en el sistema masa-
resorte-subamortiguador (MRA) para analizar la respuesta lateral (cizallamiento) de una porción idealizada
de suelo en la zona de influencia del Volcán Galeras ante la excitación generada por terremotos volcano-
tectónicos (VTs). El objetivo principal es explorar cómo las propiedades elásticas y de amortiguamiento del
suelo podrı́an influir en la amplitud y frecuencia del desplazamiento lateral durante estos eventos sı́smicos,
con la potencial aplicación futura en el desarrollo de métodos de detección temprana. Se realizará una re-
visión bibliográfica sobre las caracterı́sticas de los VTs en el Galeras y las propiedades geotécnicas tı́picas de
suelos volcánicos. Posteriormente, se implementará un modelo MRA con parámetros estimados para simular
la respuesta a una excitación sinusoidal simplificada, representando la componente horizontal de las ondas
sı́smicas. Los resultados se analizarán en términos de amplitud de desplazamiento y frecuencia de resonancia.

2 Palabras claves
• Sistema Masa-Resorte-Amortiguador .
• Sismos Volcanotectónicos Volcán Galeras.
• Desplazamiento Lateral.
• Modelado Simplificado.
• Respuesta Sı́smica.
• Detección Temprana.

3 Planteamiento del problema

El Volcán Galeras, se localiza en el sur de Colombia, en el departamento de Nariño, y alcanza una altitud
máxima de 4276 metros sobre el nivel del mar. Su ubicación geográfica aproximada es 113′ N y 7721′ W .
En un radio de 15 km alrededor del cráter se encuentran varios municipios, incluida la ciudad de San Juan de
Pasto, que albergan en conjunto una población cercana a los 400.000 habitantes, conformando una zona de
alta densidad poblacional dentro del área de influencia directa del volcán. La actividad eruptiva del Galeras
en los últimos 500 años se ha caracterizado principalmente por eventos de tipo explosivo (Calvache, 1995;
Torres, 2010), lo que, sumado a su ubicación en una región sı́smicamente activa, representa una amenaza con-
siderable para las comunidades cercanas. Además, uno de los eventos eruptivos más significativos del Volcán
Galeras en los últimos 500 años ocurrió en 1936. Esta erupción fue de tal magnitud que pudo sentirse a más
de 15 kilómetros de distancia del cráter y estuvo acompañada por la expulsión de bloques volcánicos que
alcanzaron distancias superiores a los 3 kilómetros, posteriormente, se produjo un flujo piroclástico que de-
scendió por la ladera norte del volcán, generando un evento de alto impacto geológico y ambiental (Flórez et
al,2023). De acuerdo con el Servicio Geológico Colombiano, en caso de una erupción del Volcán Galeras, los
cinco municipios más cercanos al cráter ubicados a una distancia máxima de aproximadamente 14 kilómetros
presentan el mayor nivel de afectación proyectado, debido a su cercanı́a y exposición directa a los principales
eventos volcánicos (Servicio Geológico Colombiano, 2015).

La historia demuestra que incluso volcanes aparentemente inactivos o con actividad leve pueden causar
catástrofes de gran magnitud. Un ejemplo claro es la erupción del Volcán Vesubio en el año 79 d.C., que de-
struyó completamente las ciudades de Pompeya y Herculano. Se estima que entre 15.000 y 20.000 personas
vivı́an en la zona y que aproximadamente 5.000 fallecieron debido a quemaduras, inhalación de gases tóxicos
y asfixia por cenizas. El evento provocó una columna eruptiva que alcanzó la atmósfera superior, generando
cambios climáticos que llegaron a sentirse incluso en regiones distantes como Egipto y Siria. Investigaciones
recientes han permitido reconstruir los efectos de este desastre, incluyendo la devastación total en un radio de
más de 20 km desde el cráter (Lima, 2018). Igualmente la erupción del volcán Tambora en 1815, en Indone-
sia, es uno de los eventos volcánicos más poderosos registrados en la historia. La enorme cantidad de dióxido
de azufre expulsado a la atmósfera alteró significativamente el clima global, provocando que el año 1816 se
conociera como “el año sin verano”. Las temperaturas promedio descendieron entre 1 y 2,5 °C en regiones
como Nueva Inglaterra y Europa Occidental, afectando seriamente la producción agrı́cola. En Inglaterra, se
registró el junio más frı́o con una temperatura media de solo 12,9 °C, mientras que en Francia y Alemania
las cosechas fracasaron, agravando las condiciones de pobreza y hambre, e incluso impulsando migraciones
masivas hacia América y el este europeo (Flórez et al,2023).

En Colombia un antecedente trágico y representativo del riesgo volcánico en Colombia es la erupción del
Volcán Nevado del Ruiz ocurrida el 13 de noviembre de 1985, luego de 69 años de aparente inactividad. A
pesar de que las autoridades realizaron acciones preventivas y emitieron alertas, muchos habitantes de la ciu-
dad de Armero ubicada a 50 km del cráter minimizaron el riesgo, considerando que las emisiones de cenizas
y gases eran eventos comunes. No obstante, la erupción provocó el desprendimiento de material volcánico y
lahares (flujos de lodo), como consecuencia del derretimiento parcial del glaciar del volcán. Este fenómeno
alcanzó un radio de hasta 80 km, causando la muerte de aproximadamente 20.000 personas, la desaparición
de miles más y la destrucción total de Armero. Además, se registraron otras pérdidas humanas en municipios
como Villamarı́a y Chinchiná, y daños estructurales e impactos ambientales en localidades cercanas. Ac-
tualmente, el volcán sigue siendo una amenaza activa, y se estima que una futura erupción podrı́a dispersar
cenizas hasta más de 100 km, alcanzando incluso a Bogotá (Perdomo, 2021; Páez, 2020).

En este contexto, surge la necesidad de desarrollar un modelo matemático simplificado que permita es-
tudiar la respuesta lateral (cizallamiento) del suelo frente a perturbaciones sı́smicas generadas por sismos
volcanotectónicos. Este tipo de análisis es fundamental, ya que las caracterı́sticas geotécnicas de los suelos
volcánicos, como su elasticidad y capacidad de amortiguamiento, pueden amplificar o atenuar los efectos
sı́smicos sobre las edificaciones y estructuras crı́ticas. A pesar de que el Galeras se encuentra bajo monitoreo
permanente, no existe una estimación precisa sobre el comportamiento del suelo durante eventos sı́smicos
locales, lo cual representa una limitación significativa para la toma de decisiones en términos de diseño
sismo-resistente y gestión del riesgo.

Por lo tanto, este proyecto propone la implementación de un modelo masa-resorte-amortiguador (MRA)

como herramienta de análisis preliminar para simular la respuesta del suelo ante excitaciones sı́smicas hor-
izontales, permitiendo estimar frecuencias de resonancia y amplitudes de desplazamiento que podrı́an con-
tribuir al desarrollo de estrategias de mitigación, diseño estructural y detección temprana.

4 Justificacion del proyecto

En Colombia, gran parte de la población habita en regiones ubicadas sobre zonas de convergencia tectónica,
lo que implica una exposición significativa a amenazas sı́smicas de moderada a alta intensidad. No obstante,
dicha amenaza solo se traduce en un riesgo real cuando existen edificaciones vulnerables. De hecho, los
daños estructurales o el colapso de construcciones durante o después de un sismo suelen ser las principales
causas de pérdidas humanas, afectaciones económicas y conflictos sociales. Por ello, es fundamental que las
edificaciones se diseñen para resistir terremotos severos sin colapsar, y que puedan soportar sismos frecuentes
sin comprometer su integridad estructural ni la de sus componentes no estructurales (Oviedo & Duque, 2006)

El Volcán Galeras, clasificado como un estratovolcán andesı́tico, se localiza en el suroeste de Colombia,

en el departamento de Nariño, con una altitud máxima de 4276 metros sobre el nivel del mar. Su ubicación
geográfica precisa es aproximadamente 113′ N y 7721′ W . En un radio de 15 km alrededor del cráter se en-
cuentran varios municipios, incluida la ciudad de San Juan de Pasto, que en conjunto reúnen una población
cercana a los 400.000 habitantes. Esta área conforma la zona de influencia directa del volcán. Durante los
últimos 500 años, la actividad volcánica de Galeras ha estado dominada principalmente por erupciones ex-
plosivas (Calvache, 1995; Torres, 2010). Debido a su constante actividad, el Volcán Galeras fue seleccionado
en 1990 como parte del programa Volcanes de la Década, una iniciativa promovida por la Asociación Inter-
nacional de Vulcanologı́a y Quı́mica del Interior de la Tierra, cuyo objetivo es priorizar la investigación y
gestión del riesgo en volcanes con alta actividad y poblaciones cercanas (Mapas de Amenazas Volcánicas,
2020).

De esta manera el Volcán Galeras, ubicado cerca de la ciudad de Pasto, representa una amenaza volcánica
significativa para la región. La ocurrencia de sismos volcanotectónicos es un indicador importante de la
actividad interna del volcán y del potencial de erupción. Un modelo simplificado que explore la respuesta
lateral del suelo a estos eventos podrı́a contribuir a una mejor comprensión de la interacción entre la fuente
sı́smica y el medio de propagación local. Aunque limitado en su representación de la complejidad real, este
enfoque inicial podrı́a sentar las bases para investigaciones más avanzadas y para la eventual implementación
de sistemas de monitoreo que aprovechen las caracterı́sticas especı́ficas de la respuesta del suelo a los VTs.

5 Objetivos del proyecto

5.1 Objetivo general
Desarrollar un modelo matemático simplificado basado en el sistema masa-resorte-sub amortiguador para
analizar la respuesta lateral del suelo ante una excitación sinusoidal que simule la componente horizontal de
los terremotos volcanotectónicos en la zona del Volcán Galeras.
5.2 Objetivo especifico
• Identificar las propiedades fı́sicas del suelo en la zona del Volcán Galeras a partir de fuentes bibli-
ográficas.
• Formular un modelo masa-resorte-sub amortiguador (MRA) que represente la respuesta lateral del
suelo ante excitación sı́smica.
• Simular la respuesta del modelo MRA usando un caso ficticio con parámetros representativos.

6 Marco Conceptual
El Volcán Galeras, ubicado en el departamento de Nariño (Colombia), es uno de los sistemas volcánicos más
activos del paı́s, caracterizado por una compleja interacción entre procesos magmáticos y tectónicos. El estu-
dio de Vergara, Sánchez, Cadena y Márquez (2023) analiza 120 sismos volcanotectónicos (VTs) registrados
entre 2005 y 2018, proporcionando datos esenciales para modelar la respuesta dinámica del suelo mediante
un sistema masa-resorte-amortiguador (MRA). Estos eventos sı́smicos, con magnitudes entre 0,24 y 3,34 (p.
4), se clasifican en dos categorı́as según su profundidad y contenido frecuencial: los VT-A (profundos, ¿2
km, frecuencias 5-15 Hz) y los VT-B (somero, 1-2 km, frecuencias 1-5 Hz) (Vergara et. 2023, p. 3). La
distribución espacial de los hipocentros revela una correlación con estructuras tectónicas como la falla de
Buesaco (orientación SW-NE), lo que evidencia el control estructural en la generación de sismicidad.
El modelo MRA representa una porción homogénea de suelo volcánico mediante tres parámetros funda-
mentales: la masa equivalente (m), que cuantifica la inercia del suelo; la rigidez al corte (k), asociada a la
resistencia elástica del medio; y el amortiguamiento (c), que modela la disipación de energı́a. Para el Volcán
Galeras, estos parámetros se estiman a partir de estudios geotécnicos en suelos volcánicos no consolidados
(Torres et al., 2020, citado en Vergara et al., 2023):

• Masa (m): 1500-2500 kg/m², derivada de densidades tı́picas (ρ = 1500 − 2500 kg/m³) para un espesor
de capa de 1 m.
• Rigidez (k): 10-50 MN/m, equivalente a módulos de corte (G) de 10-50 MPa en suelos piroclásticos.
√
• Amortiguamiento (c): 5000-50 000 Ns/m, calculado como c = 2ζ km, donde ζ = 5% − 20% es el
amortiguamiento crı́tico.

La excitación sı́smica se modela como una fuerza armónica F (t) = F0 sin(ωt), donde F0 depende de
la aceleración del suelo (a = 0, 1 − 1, 0 m/s² para magnitudes M < 3, 5) y ω corresponde a la frecuencia
angular de los VTs (ω = 2πf , con f = 1 − 15 Hz). La ecuación diferencial del sistema es:

d2 x dx
m 2
+c + kx = F0 sin(ωt), (1)
dt dt
donde x(t) representa el desplazamiento lateral del suelo. La solución de esta ecuación combina una re-
spuesta transitoria (que decae exponencialmente debido al amortiguamiento) y una respuesta estacionaria
dominada por la frecuencia de excitación. La amplitud máxima de la respuesta estacionaria está dada por:
 F0 /m
A= p , (2)
(ωn − ω 2 )2 + (2ζωn ω)2
2
p
donde ωn = k/m es la frecuencia natural del sistema. En el Volcán Galeras, ωn oscila entre 2-10 Hz
(Vergara et.2023, p. 12), lo que genera riesgo de resonancia cuando ω ≈ ωn , especialmente en suelos con
bajo amortiguamiento (ζ ≈ 5%). Este fenómeno amplifica desplazamientos laterales hasta ≥ 1 cm incluso
para sismos pequeños (M < 3), según simulaciones teóricas.
El análisis de desajustes acumulados realizado por ( Vergara et.2023, p. 6-9) identifica variaciones en el
estado de esfuerzos a profundidades de 3-6 km y a lo largo del segmento NE de la falla de Buesaco. Estas
variaciones, asociadas a heterogeneidades litológicas y actividad magmática, subrayan la importancia de
incorporar cambios espaciotemporales en k y c para mejorar la precisión del modelo. Por ejemplo, intrusiones
magmáticas superficiales podrı́an reducir k en zonas especı́ficas, aumentando la susceptibilidad a resonancia.
Si bien el modelo MRA simplifica la complejidad real del suelo volcánico, proporciona una base cuan-
titativa para diseñar sistemas de monitoreo que detecten desplazamientos crı́ticos (xmax ) como indicadores
tempranos de inestabilidad. Futuros desarrollos podrı́an integrar registros sı́smicos reales del Observatorio
Vulcanológico y Sismológico de Pasto (OVSP-SGC) para validar y calibrar los parámetros del modelo.

7 Metodologia
7.1 Fase de revision bibliografica
• Revisión de bases de datos cientı́ficas (Scopus, Web of Science) y reportes del SGC sobre sismos VTs
en el Galeras.
• Análisis de estudios geotécnicos en suelos volcánicos de la región de Pasto para estimar parámetros del
modelo (módulo de corte, densidad, coeficiente de amortiguamiento).

7.2 Fase de recoleccion de datos

• Definición de rangos para masa equivalente (1500 − 2500kg/m³), rigidez al corte (10-50 MPa) y amor-
tiguamiento (5-20%).
• Selección de frecuencias de excitación (1-15 Hz) y amplitudes (0.1-1 m/s²) basadas en registros históricos
de VTs.

7.3 Fase de experimentacion del modelo

• Ecuación diferencial del sistema: Se utilizó la formulación clásica del sistema masa-resorte-amortiguador
(MRA) para desplazamiento lateral:

d2 x dx
m 2
+c + kx = F0 sin(ωt) (3)
dt dt
donde:
• m: Masa equivalente por unidad de área (kg/m²)
• c: Coeficiente de amortiguamiento viscoso (Ns/m)

• k: Rigidez al corte del suelo (N/m)

• F0 : Amplitud de la fuerza sı́smica equivalente (N)
• ω: Frecuencia angular de excitación (rad/s)
• Parámetros de entrada: Se asignaron valores basados en la revisión bibliográfica:

– Masa (m): 1500-2500 kg/m³ (densidad del suelo volcánico)

– Rigidez (k): 10-50 MPa (módulo de corte para suelos piroclásticos)
√
– Amortiguamiento (c): 5-20% del crı́tico (ccr = 2 km)
– Fuerza (F0 ): 0.1-1 m/s² (aceleración sı́smica reportada para VTs)
– Frecuencias (f ): 1-15 Hz (ω = 2πf )

7.4 Desarrollo analı́tico para el caso subamortiguado

• Ecuación diferencial del sistema La ecuación del sistema masa-resorte-amortiguador con excitación
sı́smica es:
mẍ + cẋ + kx = F0 sin(ωt) (4)

Dividiendo por m y definiendo parámetros adimensionales:

c
α= (Coeficiente de amortiguamiento) (5)
2m
r
k
ωn = (Frecuencia natural no amortiguada) (6)
m
p
ωd = ωn2 − α2 (Frecuencia natural amortiguada) (7)

• Solución homogénea La ecuación homogénea asociada es:

ẍ + 2αẋ + ωn2 x = 0 (8)

Con solución general:

xh (t) = e−αt [C1 cos(ωd t) + C2 sin(ωd t)] (9)
• Solución particular por variación de parámetros Asumimos una solución particular de la forma:

xp (t) = u1 (t)e−αt cos(ωd t) + u2 (t)e−αt sin(ωd t) (10)

El sistema para determinar u′1 y u′2 es:

u′1 e−αt cos(ωd t) + u′2 e−αt sin(ωd t) = 0 (11)

−u′1 e−αt [α cos(ωd t) + ωd sin(ωd t)]
F0
+u′2 e−αt [−α sin(ωd t) + ωd cos(ωd t)] = sin(ωt) (12)
m

• Cálculo del Wronskiano El determinante del sistema (Wronskiano) es:

e−αt cos(ωd t) e−αt sin(ωd t)

W = d −αt d −αt
dt [e cos(ωd t)] dt [e sin(ωd t)]
= ωd e−2αt (13)

• Solución para u′1 y u′2 Resolviendo el sistema (11)-(12):

F0
u′1 = − sin(ωt) sin(ωd t)eαt (14)
mωd
F0
u′2 = sin(ωt) cos(ωd t)eαt (15)
mωd

• Integración de u1 y u2 Usando identidades trigonométricas para sin A sin B y sin A cos B:

 (α+i(ω+ωd ))t
e(α+i(ω−ωd ))t

F0 e
u1 (t) = − + C3 (16)
2mωd α + i(ω + ωd ) α + i(ω − ωd )
 (α+i(ω+ωd ))t
e(α+i(ω−ωd ))t

F0 e
u2 (t) = + + C4 (17)
2mωd α + i(ω + ωd ) α + i(ω − ωd )

• Solución general completa Combinando (9), (16) y (17):

F0 /m
x(t) = e−αt [C1 cos(ωd t) + C2 sin(ωd t)] + p sin(ωt − ϕ) (18)
(ωn2 − ω 2 )2 + (2αω)2

donde el ángulo de fase ϕ es:  

2αω
ϕ = arctan (19)
2
ωn − ω 2
 • Amplitud máxima y resonancia La amplitud del estado estacionario es:

F0 /m
A= p (20)
(ωn − ω 2 )2 + (2αω)2
2

La condición de resonancia ocurre cuando:

s  2
α
ωres = ωn 1−2 (21)
ωn

Con amplitud máxima:

F
Amax = p 0 (22)
2mαωn 1 − (α/ωn )2

• Implicaciones para suelos volcánicos

– Bajos valores de α (tı́picos en suelos volcánicos) generan altas amplitudes (Amax ∝ 1/α)
– La frecuencia de resonancia ωres se aproxima a ωn cuando α → 0
– Suelos blandos (k ↓⇒ ωn ↓) son más susceptibles a frecuencias bajas de VTs

8 Resultados
Parámetros del modelo
Los parámetros empleados fueron seleccionados con base en datos caracterı́sticos de suelos volcánicos y
registros de sismos volcánicos documentados en la literatura técnica y cientı́fica:

• Masa equivalente del suelo: m = 2000 kg,

• Rigidez: k = 30 MPa = 30 × 106 N/m,
• Amortiguamiento: c = 10% del valor crı́tico, con coeficiente α = 0.0401,

• Frecuencia natural del sistema: ωn = 122.47 rad/s,

• Frecuencia de excitación sı́smica: f = 5 Hz, correspondiente a ω = 31.42 rad/s,
• Amplitud de aceleración sı́smica: asismo = 0.5 m/s2 ,
• Tasa de decaimiento de la excitación: β = 1.5,

• Fase inicial de la excitación: ϕ = π/4.

Ecuación de movimiento y componentes de la respuesta
La solución general propuesta para el desplazamiento del sistema es:
F0
x(t) = e−αt [C1 cos(ωd t) + C2 sin(ωd t)] + p m
e−βt sin(ωt − ϕ),
(ωn2 − ω 2 )2 + (2αω)2
la cual incorpora dos contribuciones fundamentales:

• Respuesta libre amortiguada: representada por el primer término, modela el comportamiento natu-
ral del suelo en ausencia de fuerza externa, con disipación progresiva de energı́a determinada por el
coeficiente α,
• Respuesta forzada amortiguada con envolvente temporal: el segundo término describe la acción
sı́smica externa modulada por una envolvente exponencial e−βt , que simula el carácter transitorio de
un sismo volcánico.

Justificación del ajuste y validez del modelo

En la formulación original, la excitación forzada es armónica y permanente, con una expresión del tipo:
F0
x(t) = e−αt [C1 cos(ωd t) + C2 sin(ωd t)] + p m
sin(ωt − ϕ).
(ωn2 − ω 2 )2 + (2αω)2
 Sin embargo, esta no refleja adecuadamente la naturaleza real de los VTs, cuya intensidad disminuye
progresivamente en el tiempo. Por ello, se incorpora el factor e−βt en la parte forzada, dando lugar a una
excitación sı́smica de amplitud decreciente. Este ajuste no altera la estructura fundamental de la solución
teórica, pero permite representar de forma más realista la atenuación gradual de un evento sı́smico.

Análisis cualitativo de la respuesta

La gráfica obtenida a partir del modelo exhibe un comportamiento oscilatorio amortiguado, cuya amplitud
decrece suavemente en el tiempo. Esta evolución refleja:

• La disipación progresiva de la energı́a mecánica del sistema, como resultado del amortiguamiento
interno (α),
• La atenuación natural de la excitación sı́smica (β), en concordancia con la duración finita de los eventos
volcánicos,
• Una dinámica consistente con observaciones experimentales y estudios sobre la respuesta de suelos
piroclásticos a vibraciones inducidas por actividad volcánica.
References
Vergara, M. et al. (2023). Análisis de mecanismos focales de sismos volcanotectónicos en el volcán
Galeras, Colombia. Boletı́n Geológico, 50(2).
https://doi.org/10.32685/0120-1425/bol.geol.50.2.2023.700

Cadena Ibarra, O. E. (2021). Modelos de fuente de sismicidad LP para la actividad del volcán
Galeras 2004–2010 (Colombia) [Tesis doctoral, Universidad Nacional de Colombia]. Repositorio
Institucional UNAL.

Oviedo, Juan Andrés, & Duque, Marı́a del Pilar. (2006). SISTEMAS DE CONTROL DE RE-
SPUESTA SÍSMICA EN EDIFICACIONES. Revista EIA, (6), 105-120. Retrieved May 25, 2025, from
http://www.scielo.org.co/scielo.phpscript=sci arttext&pid=S1794-12372006000200010&lng=en&tlng=es.

Calvache Velasco, Marta. The Geological Evolution of Galeras Volcanic Complex. D. Thesis.
Arizona State University. Arizona State. 1995. 185.p.

Torres Corredor, Roberto. (2010). Locating long period events and tremor beneath Galeras vol-
cano using seismic amplitudes stochastic approach. M. S. thesis. Building Research Institute. Tokyo.

Mapa de amenaza Volcán Galeras . (2020). Obtenido de Servicio Geológico Colombiano el 23 de

noviembre de 2020:
https://www2.sgc.gov.co/sgc/volcanes/VolcanGaleras/Paginas/Mapa-deamenaza.aspx

Servicio Geológico Colombiano. (2015). Mapa de amenaza volcánica del Volcán Galeras. Di-
rección de Geoamenazas.

Flórez-Molina, M. T., Parra-Sánchez, L. N., Calderón-Gutiérrez, J. A. (2023). Evidencias de

eventos extremos de precipitación y de explosiones volcánicas entre 1455-2022 a. d. en el humedal La
Bramadora, Antioquia, Colombia. Boletı́n de Ciencias de la Tierra, 53, 63–88.

Perdomo, J. (2021). ¿Qué ocurrio en Armero? con la Ex-Funcionaria de la Defensa Civil Colom-
biana Janneth Perdomo (E. Páez, Entrevistador)

Páez, W. (2020). Riesgos y efectos ante una ameza volcanicá del Volcán Galeras con el Capitan
(B). William R. Páez Sánchez -Ex.Comandante de Bomberos Bogotá - Intructor OFDA ”CPI-EDAN” -
Intructor de Primera Categoria Defensa Civil Colombiana. (E. Páez, Entrevistador)