Skip to article frontmatterSkip to article content
Site not loading correctly?

This may be due to an incorrect BASE_URL configuration. See the MyST Documentation for reference.

Casas en el borde de una falla, con ondas sísmicas

Este capítulo sigue el molde de siete puntos (Molde de capítulo de amenaza). A diferencia de la inundación, el sismo es una amenaza geofísica: no depende del clima y su evaluación descansa sobre la ingeniería sísmica y no sobre la hidrología. Por eso este capítulo usa más el Atlas de Riesgo de Colombia y la literatura nacional que los marcos CLIMAAX y CCDR-tools, que se centran en amenazas hidrometeorológicas GFDRR and World Bank, 2025.

1. Fenómeno

Los sismos son movimientos del suelo producidos por una liberación repentina de energía. Las fuerzas tectónicas y convectivas del interior de la Tierra acumulan esfuerzos de forma gradual en el contacto entre placas o en fallas de la corteza. Cuando se alcanza la resistencia de la roca, la energía se libera como ondas que se propagan en todas direcciones y hacen vibrar el terreno Unidad Nacional para la Gestión del Riesgo de Desastres, 2019.

Las ondas de cuerpo (P, de compresión, y S, de corte) son las responsables de las mayores afectaciones sobre lo construido. Las ondas de superficie (Love y Rayleigh) producen movimientos horizontales de corte y movimientos elípticos verticales en la superficie Unidad Nacional para la Gestión del Riesgo de Desastres, 2019. La intensidad se atenúa con la distancia a la fuente porque las ondas pierden energía al viajar, de modo que aceleración, velocidad y desplazamiento del terreno son distintos en cada punto Unidad Nacional para la Gestión del Riesgo de Desastres, 2019.

Marco tectónico de Colombia

Colombia se ubica en la convergencia de tres placas: Nazca, Caribe y Suramérica. Su movimiento relativo originó el relieve y la estructura de las cordilleras, y la convergencia entre Nazca y Suramérica también se absorbe en grandes sistemas de fallas intracontinentales, responsables de la sismicidad superficial Unidad Nacional para la Gestión del Riesgo de Desastres, 2019. Entre los sistemas de fallas más destacados están el de Romeral y el del Piedemonte Llanero, que se extienden desde Ecuador hasta el norte del país Unidad Nacional para la Gestión del Riesgo de Desastres, 2019.

Hay además fallas con poco peso en la sismicidad nacional pero muy importantes a escala local, que deben entrar en los estudios de microzonificación Unidad Nacional para la Gestión del Riesgo de Desastres, 2019.

Variable de intensidad

La amenaza se expresa con parámetros de movimiento fuerte: aceleración máxima del terreno (PGA), velocidad (PGV), desplazamiento (PGD) y ordenadas del espectro de respuesta. El PGA es el parámetro más empleado Unidad Nacional para la Gestión del Riesgo de Desastres, 2019. Para amenazas secundarias, como la licuefacción, se usan además variables geotécnicas (ver el punto 2).

Definición según los HIPs de UNDRR

Los perfiles de información de amenazas (Perfiles de información de amenazas (HIPs) de UNDRR) que corresponden a este capítulo son:

CódigoPerfil (clúster)Definición (traducción propia)
GH0101Terremoto (geológico, sísmico)Término que describe el deslizamiento repentino de una falla y la sacudida del suelo que produce la energía sísmica radiada durante ese evento. El deslizamiento puede deberse a cambios de esfuerzo en la Tierra o a actividad volcánica o magmática (USGS).
GH0307Licuefacción (geológico, falla del terreno)Pérdida de resistencia de sedimentos sueltos, saturados o casi saturados, en o cerca de la superficie, por sacudidas fuertes como las de un sismo, por cargas cíclicas o vibraciones, o por exceso de presión de poros (AGI; USGS).
GH0311Ruptura superficial y fisuramiento (geológico, otros)Desplazamientos localizados del terreno por esfuerzos de tensión, compresión y corte, sobre todo en sedimentos no consolidados; la ruptura superficial es la prolongación hacia arriba del deslizamiento de la falla en profundidad (Arizona Geological Survey; USGS).

Fuente: UNDRR & ISC, 2025.

Qué se incluye y por qué.

2. Amenaza

Por qué probabilista

Los terremotos pequeños (Mw ≤ 3,5) son mucho más frecuentes que los grandes (Mw ≥ 7,5), pero no es posible estimar cuándo, dónde ni con qué magnitud ocurrirá el próximo Unidad Nacional para la Gestión del Riesgo de Desastres, 2019. La alta incertidumbre motivó los análisis probabilistas de amenaza sísmica, que estiman a largo plazo las posibles localizaciones y frecuencias de los terremotos y propagan las incertidumbres Unidad Nacional para la Gestión del Riesgo de Desastres, 2019.

Etapas del análisis probabilista

El método fue establecido por Cornell (1968) y Esteva (1970) y se resume en cuatro pasos, según Reiter (1990) Unidad Nacional para la Gestión del Riesgo de Desastres, 2019:

  1. Caracterizar las fuentes que influyen en el sitio, por su geometría y por la distribución de probabilidad de inicio de la ruptura. Es usual modelarlas como planos y asignar probabilidad uniforme dentro de cada uno.

  2. Definir la sismicidad de cada fuente a partir del catálogo sísmico y sus características sismotectónicas.

  3. Estimar el movimiento fuerte en el sitio con funciones de atenuación que consideren su incertidumbre.

  4. Combinar las incertidumbres de localización, tamaño y atenuación para obtener las tasas anuales de excedencia de los parámetros de movimiento fuerte.

La tasa de excedencia de una intensidad en un punto es la suma ponderada de las probabilidades de excedencia sobre todos los escenarios que definen la amenaza, con las frecuencias anuales de los terremotos como pesos. Esos escenarios deben ser estocásticos, mutuamente excluyentes y colectivamente exhaustivos Unidad Nacional para la Gestión del Riesgo de Desastres, 2019. En notación habitual, para una intensidad aa y fuentes ii:

ν(a)  =  ∑i∫M0Mu∣dλi(M)dM∣ P(A>a∣M,Ri) dM\nu(a) \;=\; \sum_i \int_{M_0}^{M_u} \left|\frac{d\lambda_i(M)}{dM}\right| \, P\big(A>a \mid M, R_i\big)\, dM

donde λi(M)\lambda_i(M) es la tasa de excedencia de magnitudes de la fuente ii y P(A>a∣M,Ri)P(A>a\mid M,R_i) la probabilidad de que la intensidad supere aa dado un sismo de magnitud MM a la distancia RiR_i. La forma de la ecuación es la clásica de Cornell y Esteva. El Atlas la describe con palabras y la sustenta en el teorema de la probabilidad total Unidad Nacional para la Gestión del Riesgo de Desastres, 2019.

Sismicidad de las fuentes

Es común suponer que los terremotos de una fuente siguen un proceso de Poisson: el número de ocurrencias en una ventana de tiempo es independiente del de otra, es decir, el proceso no tiene memoria. La tasa de excedencia de magnitudes λ(M)\lambda(M) suele ser una versión modificada de la relación de Gutenberg y Richter, definida por tres parámetros que se estiman por máxima verosimilitud e información experta Unidad Nacional para la Gestión del Riesgo de Desastres, 2019:

ParámetroSignificado
β\betaPendiente de la curva de recurrencia de magnitudes en la zona de magnitudes bajas
MuM_uMagnitud máxima, según la máxima longitud de ruptura posible y las características morfotectónicas
λ0\lambda_0Número promedio de eventos por año con magnitud mayor que la umbral M0M_0

Ejemplo con datos reales: sismicidad con el catálogo del USGS

El cuaderno Sismicidad de Colombia con el catálogo del USGS descarga el catálogo del USGS para Colombia y su área de influencia (2000–2024, magnitud de 3,0 o más) y muestra cómo se obtienen los insumos de esta sección. Está adaptado del cuaderno de la Subdirección para el Conocimiento del Riesgo de la UNGRD, con licencia CC BY 4.0 (cuadernos-reproducibles-crc004). Con los datos descargados al ejecutarlo:

ResultadoValor
Eventos del catálogo4 201
Sismos superficiales (hasta 70 km) / intermedios (70 a 300 km)2 290 / 1 910
Magnitud de completitud (McM_c, máximo de la curvatura)4,3
bb por máxima verosimilitud (Aki-Utsu), sismos con M≥McM \ge M_c1,06 ± 0,02
bb por regresión sobre la curva acumulada0,93

Lo que enseña el ejemplo:

Los valores salen de una descarga del 2 de noviembre de 2026; el USGS puede revisar magnitudes, por lo que las cifras pueden variar un poco al repetirlo.

Atenuación del movimiento fuerte

La intensidad se modela como variable aleatoria, usualmente con distribución lognormal cuyos momentos entrega la función de atenuación. Estas funciones dependen de la magnitud, de la distancia fuente–sitio y del mecanismo focal de ruptura Unidad Nacional para la Gestión del Riesgo de Desastres, 2019. Existen varias definiciones de distancia (epicentral, hipocentral, a la superficie de falla) y la elección debe corresponder al modelo adoptado Unidad Nacional para la Gestión del Riesgo de Desastres, 2019.

Para el norte de Suramérica hay modelos regionalizados. El modelo NoSAm para sismos corticales regionaliza el modelo global ASK14 de NGA-West2, corrige su desajuste con los datos observados para magnitudes moderadas y distancias intermedias y grandes, y añade atenuación para trayectorias que cruzan el arco volcánico Arteta et al., 2023.

Amenaza sísmica nacional y local

La evaluación nacional que respalda la NSR-10 siguió cinco pasos: modelar las fuentes sismogénicas, asignar los eventos de los catálogos a cada fuente, caracterizar cada fuente con parámetros estadísticos, aplicar el modelo probabilista (software CRISIS V7.2) e interpretar los resultados para el diseño Unidad Nacional para la Gestión del Riesgo de Desastres, 2019. El Atlas muestra el mapa para un periodo de retorno de 475 años y un periodo estructural de 0,5 s, calculado en roca firme. Sirve como insumo para microzonificación y estudios de sitio, no como valor directo en el suelo de la ciudad Unidad Nacional para la Gestión del Riesgo de Desastres, 2019.

Amenazas secundarias

Un marco probabilista de amenaza de licuefacción para sitios con una capa superior no licuable sobre una capa licuable se aplicó al distrito de Aguablanca, en Cali, un área de susceptibilidad moderada a alta. Está pensado para regiones con pocos datos, porque aprovecha la información sísmica y geotécnica local Cruz et al., 2026. Los tsunamis de origen sísmico se tratan en Tsunami.

3. Exposición

La exposición sísmica es el inventario de edificaciones e infraestructura con su ubicación, tipología estructural, material, altura, edad y uso. Es la parte que más varía entre ciudades.

Un ejemplo colombiano de construcción del modelo de exposición para una provincia combinó información de la Fundación del Modelo Global de Terremotos (GEM), encuestas y el censo nacional. Se asignaron curvas de fragilidad y vulnerabilidad a los tipos de edificios de la región Feliciano et al., 2023. Ese enfoque es replicable en otras ciudades cuando no hay catastro con tipología estructural.

La edad es un atributo clave: una parte significativa del inventario nacional data de antes del primer código de diseño sísmico obligatorio (Decreto 1400 de 1984) Feliciano et al., 2023Unidad Nacional para la Gestión del Riesgo de Desastres, 2019.

4. Vulnerabilidad

La vulnerabilidad física sísmica se expresa con funciones de fragilidad (probabilidad de superar un estado de daño dada la intensidad) y funciones de vulnerabilidad (relación entre intensidad y pérdida esperada). El Atlas exige definir cada función por su valor esperado y su varianza, porque las pérdidas por una misma intensidad no son idénticas entre edificios similares Unidad Nacional para la Gestión del Riesgo de Desastres, 2019.

Hay estudios recientes de fragilidad y vulnerabilidad por tipología constructiva de Colombia: mampostería confinada Archbold et al., 2025 y mampostería no reforzada informal Reyes et al., 2025. La elección entre derivar curvas propias o usar curvas genéricas afecta el resultado: un estudio metodológico sobre varias ciudades con distintos niveles de amenaza examina los sesgos de los métodos de derivación de la vulnerabilidad en lugares con poca información Hoyos & Hernandez, 2022.

Para la dimensión social de la vulnerabilidad se ofrecen índices demográficos aplicados a la pérdida directa (ver Vulnerabilidad social y capacidades). El estudio de la provincia de Sabana Centro incorporó uno Feliciano et al., 2023.

5. Riesgo

El flujo es el del Atlas: catálogo de eventos estocásticos, exposición, vulnerabilidad y cálculo de pérdidas. Sus resultados son la curva de excedencia de pérdidas, la PAE y la PMP (Unidad Nacional para la Gestión del Riesgo de Desastres, 2019; ver Evaluación probabilista y pérdida anual esperada).

Dos modalidades se usan en la práctica:

ModalidadUsoEjemplo
ProbabilistaPAE, PMP, seguros, priorizaciónAtlas de Riesgo Unidad Nacional para la Gestión del Riesgo de Desastres, 2019
EscenariosPlanes de emergencia, comunicaciónDieciocho escenarios de sismo de 475 años de periodo de retorno simulados con la herramienta OpenQuake para Sabana Centro Feliciano et al., 2023

Incertidumbre epistémica. Combinar varios modelos de amenaza produce una incertidumbre epistémica que afecta las estimaciones probabilistas de riesgo sísmico Marulanda et al., 2021.

6. Clima

Los sismos no son sensibles al clima. CCDR-tools lo plantea así: el cambio climático no afecta de forma significativa las amenazas geofísicas, que incluyen terremotos, tsunamis y actividad volcánica, aunque sí puede afectar las hidrometeorológicas GFDRR and World Bank, 2025.

Para la amenaza, el Atlas la considera estacionaria y sus órdenes de magnitud normalmente no cambian en la historia. Aun así, los estudios deben actualizarse periódicamente porque hay más registros de movimiento fuerte y mejores funciones de atenuación Unidad Nacional para la Gestión del Riesgo de Desastres, 2019.

Lo que sí cambia con el tiempo es la exposición y, en menor medida, la vulnerabilidad (crecimiento urbano, reforzamiento, nueva normativa). Por eso los escenarios de este capítulo son de exposición, no climáticos. Esto es un criterio de método de este libro, no una afirmación de las fuentes citadas.

7. Datos y medidas

Datos disponibles

FuenteQué ofreceReferencia
Servicio Geológico Colombiano (SGC)Catálogo e información histórica de sismos, mapas de sismicidad, Red Sismológica Nacional y Red Nacional de AcelerógrafosUnidad Nacional para la Gestión del Riesgo de Desastres, 2019
Red Nacional53 estaciones sismológicas en zonas de alta actividad y 103 acelerógrafos (cifras de 2019)Unidad Nacional para la Gestión del Riesgo de Desastres, 2019
OSSOSismómetros en Valle del Cauca, Cauca y Risaralda, con reportes diarios y mensualesUnidad Nacional para la Gestión del Riesgo de Desastres, 2019
Redes urbanasLISA (Bogotá) y SisMan LISA (Manizales), acelerógrafos interconectadosUnidad Nacional para la Gestión del Riesgo de Desastres, 2019
AISMapas de amenaza, microzonificación y normas de la NSRUnidad Nacional para la Gestión del Riesgo de Desastres, 2019
GEMModelos de exposición y fragilidad de edificacionesFeliciano et al., 2023

Las cifras de estaciones son de 2019 y deben actualizarse al aplicarlas.

Normativa

El Código Colombiano de Construcciones Sismo Resistentes se expidió con el Decreto 1400 de 1984, después del terremoto de Popayán de 1983. Fue actualizado con el Decreto 33 de 1998 (NSR-98) y con el Decreto 926 de 2010 (NSR-10) Unidad Nacional para la Gestión del Riesgo de Desastres, 2019. La actualización NSR-22 se prepara en el comité AIS-300, que propone coeficientes sísmicos de diseño; sus autores critican el uso del movimiento fuerte orientado al riesgo (RTGM) por transmitir una falsa seguridad Bernal & Cardona, 2021.

Incertidumbres principales

Medidas según la Ley 1523

ProcesoMedidaFuente
ConocimientoMicrozonificación sísmica, redes de acelerógrafos, escenarios de riesgoUnidad Nacional para la Gestión del Riesgo de Desastres, 2019Feliciano et al., 2023
Reducción: prevenciónNormas de diseño y construcción sismo resistente, ordenamiento territorial según la amenazaUnidad Nacional para la Gestión del Riesgo de Desastres, 2019
Reducción: mitigaciónReforzamiento de edificaciones existentes; planes de mitigación para infraestructura escolarUnidad Nacional para la Gestión del Riesgo de Desastres, 2019Fernández et al., 2021
ManejoAtención de emergencias apoyada en las redes de acelerógrafosUnidad Nacional para la Gestión del Riesgo de Desastres, 2019

Síntesis de evidencia

BloqueEstudioDiseño y lugarResultado relevante
Fuente sísmicaChicangana-Montón et al. (2022)Sismicidad y sismotectónica del Sistema de Fallas de AlgecirasCuatro sismos históricos con M > 6,5 desde finales del siglo XVIII, uno de ellos el de 1967 (Mw 7,2)
Fuente sísmicaSalcedo-Hurtado et al. (2021); Di Giacomo & Sarabia (2021)Reevaluación macrosísmica del sismo de 1925 y de otros previos a la era digitalLos datos macrosísmicos permiten reubicar y redimensionar sismos históricos
Fuente sísmicaMuñoz Lopez et al. (2020); Jarrin et al. (2022)Calibración de la magnitud local en cinco zonas; movimiento de la placa de Nazca con GPSMejoran el catálogo y el marco tectónico
Fuente sísmicaGarcía-Arias et al. (2026)Aprendizaje automático para el Nido Sísmico de BucaramangaEl nido es uno de los tres volúmenes con más sismicidad intermedia del mundo y cercano a centros urbanos
Movimiento del sueloArteta et al. (2023)Modelo de movimiento del suelo para sismos corticales del norte de SuraméricaLos sismos corticales cercanos a ciudades contribuyen mucho a la amenaza; el sismo del Eje Cafetero (1999, Mw 6,2) causó daños del 1,9 % del PIB
Movimiento del sueloGómez-García et al. (2023); Londoño et al. (2022); Celorio et al. (2025)Espesor sismogénico cortical; atenuación en el Valle Medio del Magdalena; longitudes de atenuación urbanasInsumos para modelos de atenuación y de efectos de ciudad
Efectos de sitioMercado et al. (2024)Parámetros estimadores de efectos de sitio, norte de SuraméricaCompara estimadores simplificados de la amplificación local
VulnerabilidadChávez-García et al. (2021)Informe de daños de Armenia (1999) frente al número de pisos y edad de la construcción, con dos evaluaciones previas (1993 y 2004)El daño fue desproporcionado (aceleración pico de 580 gal) y se relaciona con la vulnerabilidad y con los efectos de sitio
VulnerabilidadFernández et al. (2022)Metodología simplificada por componentes para escuelasLos daños en escuelas se deben también a diseño y construcción deficientes
RiesgoAcevedo et al. (2020)Evaluación de riesgo sísmico y escenarios para edificaciones residenciales de Bogotá (8,0 millones), Medellín (2,5) y Cali (2,4), con metodología uniformePermite comparar las tres ciudades con los mismos modelos de amenaza, vulnerabilidad y exposición
RiesgoHoyos & Hernandez (2022); Marulanda et al. (2021)Sesgos al derivar vulnerabilidad en varias ciudades; incertidumbre epistémica en riesgo sísmicoLa derivación de vulnerabilidad se sesga según el nivel de amenaza; la incertidumbre epistémica es relevante en estimados probabilistas
Exposición y respuestaHurtado et al. (2026); Aroquipa & Hurtado (2022); Ordaz et al. (2023)Predicción de colapso con aprendizaje automático; resiliencia de edificios; aceleraciones de diseñoEnlazan daño observado, resiliencia y normativa
ComunicaciónAcevedo et al. (2025)Iniciativa de transferencia de conocimiento de riesgo sísmicoPara que el riesgo sirva a la toma de decisiones, quien lo genera debe transferirlo a los usuarios finales

Lectura conjunta. Los estudios revisados sugieren, con confianza media, que el riesgo sísmico urbano se calcula con un marco común pero cambia por tres factores que el análisis de cada ciudad debe resolver: las fuentes corticales cercanas, los efectos de sitio y los sesgos al derivar la vulnerabilidad.

Pendiente

Bibliografía recomendada

Lista de los 194 registros de la hoja de vigilancia científica (pestaña Vigilancia_CTeI, columna AMENAZA_ESPECIFICA = «Terremoto (sismo)»; «Licuefacción»), en formato APA 7 con metadatos de Crossref. Son lecturas recomendadas: no están verificadas contra el texto completo ni se usan como evidencia hasta que se citen en el capítulo. Las entradas sin DOI carecen de autores en la hoja.

References
  1. GFDRR and World Bank. (2025). Risk Data Library: Country Climate and Disaster Risk Screening Tools (CCDR-tools). https://gfdrr.github.io/CCDR-tools/home.html
  2. Unidad Nacional para la Gestión del Riesgo de Desastres. (2019). Atlas de Riesgo de Colombia: revelando los desastres latentes. UNGRD.
  3. UNDRR, & ISC. (2025). 2025 Hazard Information Profiles (HIPs) online reference. PreventionWeb, DRR glossary. https://www.preventionweb.net/drr-glossary/hips
  4. Arteta, C. A., Pajaro, C. A., Mercado, V., Montejo, J., Arcila, M., & Abrahamson, N. A. (2023). Ground-Motion Model (GMM) for Crustal Earthquakes in Northern South America (NoSAm Crustal GMM). Bulletin of the Seismological Society of America, 113(1), 186–203. 10.1785/0120220168
  5. Cruz, A., Karimzadeh, S., Chieffo, N., Sandoval, E., & Lourenço, P. B. (2026). Probabilistic liquefaction hazard assessment in urban centres with a non-liquefiable surface layer overlying a liquefiable stratum. Bulletin of Earthquake Engineering, 24(9), 6809–6841. 10.1007/s10518-026-02512-x
  6. Feliciano, D., Arroyo, O., Cabrera, T., Contreras, D., Valcárcel Torres, J. A., & Gómez Zapata, J. C. (2023). Seismic risk scenarios for the residential buildings in the Sabana Centro province in Colombia. Natural Hazards and Earth System Sciences, 23(5), 1863–1890. 10.5194/nhess-23-1863-2023
  7. Archbold, J., Carrillo, J., Orjuela Mejia, J., & Piedrahita, J. (2025). Seismic fragility assessment of confined masonry houses in Colombia. Earthquake Spectra, 41(5), 3586–3611. 10.1177/87552930251378782
  8. Reyes, J. C., Torres, E., Portilla, S., Riaño, A. C., & Paul Smith-Pardo, J. (2025). Analytical seismic fragility and vulnerability curves for Colombian informal unreinforced masonry buildings. Earthquake Spectra, 41(5), 3556–3585. 10.1177/87552930251360037
  9. Hoyos, M. C., & Hernandez, A. F. (2022). Seismic Risk Assessment of Multiple Cities: Biases in the Vulnerability Derivation Methods for Urban Areas With Different Hazard Levels. Frontiers in Earth Science, 10. 10.3389/feart.2022.910118
  10. Marulanda, M. C., de la Llera, J. C., Bernal, G. A., & Cardona, O. D. (2021). Epistemic uncertainty in probabilistic estimates of seismic risk resulting from multiple hazard models. Natural Hazards, 108(3), 3203–3227. 10.1007/s11069-021-04820-8
  11. Bernal, G., & Cardona, O. D. (2021). Actualización de los coeficientes sísmicos de diseño estructural para la Norma Colombiana de Construcción Sismo Resistente NSR-22. In Investigaciones en Gestión del Riesgo de Desastres en Colombia 2021. 10.5281/zenodo.21627086
  12. Fernández, R., Yamin, L., Reyes, J. C., Garcia, A., Fuentes, G., & Echeverry, J. (2021). Mitigación del riesgo sísmico de la infraestructura escolar. In Investigaciones en Gestión del Riesgo de Desastres en Colombia 2021. 10.5281/zenodo.21627106
  13. Chicangana-Montón, G., Bocanegra-Gómez, A., Pardo-Mayorga, J., Salcedo-Hurtado, E. de J., Gómez-Capera, A., & Vargas-Jiménez, C. A. (2022). Sismicidad y sismotectónica para el sector norte del ámbito del Sistema de Fallas de Algeciras, Cordillera Oriental, Colombia. Boletín de Geología, 44(1). 10.18273/revbol.v44n1-2022005
  14. Salcedo-Hurtado, E. de J., Blandón-Silva, J. C., & Ríos-García, L. (2021). Reevaluación de intensidades macrosísmicas del terremoto del 07 de junio de 1925 en Colombia. Boletín de Geología, 43(3). 10.18273/revbol.v43n3-2021002
  15. Di Giacomo, D., & Sarabia, A. M. (2021). Use of macroseismic and instrumental data to reassess earthquake locations: Examples from pre-digital earthquakes in Colombia. Journal of South American Earth Sciences, 111, 103467. 10.1016/j.jsames.2021.103467