Resumen del informe Se presenta un análisis exploratorio de los eventos sísmicos con magnitud igual o superior a 4.5 registrados entre 2006 y 2026 por el catálogo del Servicio Geológico de los Estados Unidos (USGS) en el norte de Suramérica y el sur de Centroamérica. El informe caracteriza los sismos según su magnitud y profundidad, su distribución espacial, su comportamiento temporal, y cierra con conclusiones y recomendaciones orientadas a la gestión del riesgo sísmico.
Los sismos son liberaciones súbitas de energía acumulada en la corteza terrestre, generalmente asociadas al movimiento relativo entre placas tectónicas. Aunque no es posible predecir con exactitud cuándo ni dónde ocurrirá el próximo evento, el análisis sistemático de catálogos históricos de sismicidad permite identificar zonas de mayor actividad, rangos de magnitud y profundidad más frecuentes, y patrones temporales que resultan clave para la planificación del riesgo, el diseño de normas de construcción sismorresistente y la preparación de la comunidad.
Los datos utilizados en este informe provienen del catálogo sísmico del USGS (United States Geological Survey), una de las agencias de referencia mundial en monitoreo sismológico. El USGS integra información de redes sismográficas globales y regionales para reportar, de forma casi inmediata, la ubicación (latitud, longitud y profundidad), la hora y la magnitud de los sismos detectados en todo el planeta.
La base de datos analizada se concentra principalmente en Colombia, Ecuador, Venezuela, Perú y Panamá, una región ubicada en el límite de convergencia entre varias placas tectónicas: la placa de Nazca, que subduce por debajo de la placa Suramericana a lo largo de la costa Pacífica (generando la llamada zona de subducción andina), y la placa del Caribe, cuyo límite con la placa Suramericana atraviesa el norte de Venezuela, Panamá y el noroccidente de Colombia. Esta combinación de límites de placa convierte a la región en una de las zonas con mayor actividad sísmica del continente, con capacidad de generar tanto sismos superficiales como sismos de profundidad intermedia asociados a la placa que subduce.
Caracterizar el comportamiento de los sismos de magnitud ≥ 4.5 registrados en la región entre 2006 y 2026, en términos de su magnitud, profundidad, ubicación geográfica y distribución en el tiempo, con el fin de identificar posibles patrones espaciales o temporales útiles para la gestión del riesgo sísmico.
ggplot(sismos, aes(x = mag)) +
geom_histogram(binwidth = 0.2, fill = "#2166ac", color = "#08306b", alpha = 0.9) +
labs(title = "Distribución de la magnitud de los sismos",
x = "Magnitud", y = "Frecuencia") +
theme_minimal(base_size = 13) +
theme(plot.title = element_text(color = "#08306b", face = "bold"))ggplot(sismos, aes(x = mag)) +
geom_density(fill = "#6baed6", color = "#08306b", alpha = 0.7) +
labs(title = "Densidad de la magnitud de los sismos",
x = "Magnitud", y = "Densidad") +
theme_minimal(base_size = 13) +
theme(plot.title = element_text(color = "#08306b", face = "bold"))| Overall (N=1173) |
|
|---|---|
| mag | |
| Mean (SD) | 4.78 (0.397) |
| Median [Min, Max] | 4.60 [4.50, 7.50] |
La magnitud de los sismos registrados oscila entre 4.5 y 7.5, con una fuerte concentración en el rango de sismos ligeros a moderados. Esto es esperable: el catálogo filtra sismos de magnitud igual o superior a 4.5, y en general la ocurrencia de sismos disminuye rápidamente a medida que aumenta la magnitud (relación de Gutenberg-Richter), por lo que los eventos mayores (magnitud ≥ 7) son, en comparación, mucho menos frecuentes.
Usando los rangos estándar de clasificación sísmica se agrupan los eventos así:
sismos$categoria <- cut(
sismos$mag,
breaks = c(4, 4.9, 5.9, 6.9, 7.9),
labels = c("ligero", "moderado", "fuerte", "mayor")
)
table1(~ categoria, data = sismos)| Overall (N=1173) |
|
|---|---|
| categoria | |
| ligero | 940 (80.1%) |
| moderado | 206 (17.6%) |
| fuerte | 20 (1.7%) |
| mayor | 7 (0.6%) |
ggplot(sismos, aes(x = categoria, fill = categoria)) +
geom_bar(color = "#08306b") +
scale_fill_manual(values = c(
"ligero" = "#c6dbef",
"moderado" = "#6baed6",
"fuerte" = "#2166ac",
"mayor" = "#08306b"
)) +
labs(title = "Número de sismos por categoría de magnitud",
x = "Categoría", y = "Número de sismos") +
theme_minimal(base_size = 13) +
theme(plot.title = element_text(color = "#08306b", face = "bold"),
legend.position = "none")La profundidad a la que ocurre un sismo es tan relevante como su magnitud: sismos de igual magnitud pueden tener efectos muy distintos en superficie según qué tan profundo se originen. Los sismos superficiales (menos de ~70 km) suelen ser los más destructivos para la infraestructura, mientras que los de profundidad intermedia o profunda (asociados aquí a la subducción de la placa de Nazca) suelen sentirse de forma más atenuada en superficie, aunque pueden percibirse en un área más extensa.
ggplot(sismos, aes(x = depth)) +
geom_histogram(binwidth = 20, fill = "#2166ac", color = "#08306b", alpha = 0.9) +
labs(title = "Distribución de la profundidad de los sismos",
x = "Profundidad (km)", y = "Frecuencia") +
theme_minimal(base_size = 13) +
theme(plot.title = element_text(color = "#08306b", face = "bold"))| Overall (N=1173) |
|
|---|---|
| depth | |
| Mean (SD) | 76.1 (63.8) |
| Median [Min, Max] | 51.3 [1.50, 437] |
ggplot(sismos, aes(x = depth, y = mag)) +
geom_point(color = "#2166ac", alpha = 0.5) +
geom_smooth(method = "loess", color = "#08306b", se = FALSE) +
labs(title = "Relación entre profundidad y magnitud",
x = "Profundidad (km)", y = "Magnitud") +
theme_minimal(base_size = 13) +
theme(plot.title = element_text(color = "#08306b", face = "bold"))La mayoría de los sismos registrados son relativamente superficiales a intermedios, con una cola de eventos que superan los 200 km de profundidad, consistentes con la actividad de la placa de Nazca en subducción. No se observa una relación fuerte y sistemática entre profundidad y magnitud: hay sismos moderados tanto en superficie como en profundidad, lo que confirma que ambas variables aportan información complementaria e independiente sobre el evento.
| Overall (N=1173) |
|
|---|---|
| pais | |
| Colombia | 671 (57.2%) |
| Colombia-Ecuador border region | 2 (0.2%) |
| Ecuador | 232 (19.8%) |
| near the coast of Ecuador | 2 (0.2%) |
| near the coast of Venezuela | 1 (0.1%) |
| near the north coast of Colombia | 1 (0.1%) |
| near the west coast of Colombia | 3 (0.3%) |
| northern Colombia | 8 (0.7%) |
| northern Peru | 4 (0.3%) |
| Panama | 56 (4.8%) |
| Panama-Colombia border region | 1 (0.1%) |
| Peru | 89 (7.6%) |
| Peru-Ecuador border region | 2 (0.2%) |
| south of Panama | 2 (0.2%) |
| Venezuela | 99 (8.4%) |
leaflet(sismos) %>%
addProviderTiles(providers$CartoDB.Positron) %>%
addCircleMarkers(
lng = ~longitude, lat = ~latitude,
radius = ~ifelse(mag < 5, 4, ifelse(mag < 6, 6, ifelse(mag < 7, 8, 11))),
color = "#08306b",
fillColor = "#2166ac",
fillOpacity = 0.6,
stroke = TRUE, weight = 1,
popup = ~paste0("<b>", place, "</b><br>Magnitud: ", mag, "<br>Profundidad: ", round(depth,1), " km")
)Dado el elevado número de eventos, un mapa de agrupamiento (clustering) facilita identificar visualmente las zonas de mayor concentración de sismos:
leaflet(sismos) %>%
addProviderTiles(providers$CartoDB.Positron) %>%
addCircleMarkers(
lng = ~longitude, lat = ~latitude,
color = "#08306b", fillColor = "#2166ac", fillOpacity = 0.6,
clusterOptions = markerClusterOptions()
)La distribución espacial muestra una alta concentración de eventos en Colombia (particularmente en el suroccidente, cerca de la costa Pacífica y la frontera con Ecuador) y en Ecuador, seguidos por Venezuela, Perú y Panamá. Este patrón coincide con el trazado de los límites de placa descritos en la introducción: la mayor densidad de sismos se ubica a lo largo del frente de subducción de la placa de Nazca (occidente de Colombia y Ecuador) y del límite con la placa del Caribe (norte de Venezuela y Panamá), más que distribuirse de forma aleatoria en el territorio.
| Overall (N=1173) |
|
|---|---|
| mes | |
| Jan | 108 (9.2%) |
| Feb | 118 (10.1%) |
| Mar | 108 (9.2%) |
| Apr | 79 (6.7%) |
| May | 104 (8.9%) |
| Jun | 86 (7.3%) |
| Jul | 86 (7.3%) |
| Aug | 95 (8.1%) |
| Sep | 97 (8.3%) |
| Oct | 95 (8.1%) |
| Nov | 96 (8.2%) |
| Dec | 101 (8.6%) |
ggplot(sismos, aes(x = mes)) +
geom_bar(fill = "#2166ac", color = "#08306b") +
labs(title = "Número de sismos por mes del año",
x = "Mes", y = "Número de sismos") +
theme_minimal(base_size = 13) +
theme(plot.title = element_text(color = "#08306b", face = "bold"))| Overall (N=1173) |
|
|---|---|
| as.factor(anio) | |
| 2006 | 13 (1.1%) |
| 2007 | 52 (4.4%) |
| 2008 | 37 (3.2%) |
| 2009 | 46 (3.9%) |
| 2010 | 42 (3.6%) |
| 2011 | 47 (4.0%) |
| 2012 | 36 (3.1%) |
| 2013 | 45 (3.8%) |
| 2014 | 82 (7.0%) |
| 2015 | 58 (4.9%) |
| 2016 | 48 (4.1%) |
| 2017 | 55 (4.7%) |
| 2018 | 72 (6.1%) |
| 2019 | 72 (6.1%) |
| 2020 | 74 (6.3%) |
| 2021 | 68 (5.8%) |
| 2022 | 62 (5.3%) |
| 2023 | 88 (7.5%) |
| 2024 | 64 (5.5%) |
| 2025 | 64 (5.5%) |
| 2026 | 48 (4.1%) |
ggplot(sismos, aes(x = as.factor(anio))) +
geom_bar(fill = "#2166ac", color = "#08306b") +
labs(title = "Número de sismos por año",
x = "Año", y = "Número de sismos") +
theme_minimal(base_size = 13) +
theme(plot.title = element_text(color = "#08306b", face = "bold"),
axis.text.x = element_text(angle = 45, hjust = 1))ggplot(sismos, aes(x = time, y = mag)) +
geom_point(color = "#2166ac", alpha = 0.5) +
geom_smooth(method = "loess", color = "#08306b", se = FALSE) +
labs(title = "Magnitud de los sismos a lo largo del tiempo",
x = "Fecha", y = "Magnitud") +
theme_minimal(base_size = 13) +
theme(plot.title = element_text(color = "#08306b", face = "bold"))El número de sismos por mes se distribuye de forma relativamente uniforme a lo largo del año, sin un patrón estacional marcado, lo cual es consistente con el origen tectónico (y no climático) de estos eventos. Por año, el número de sismos registrados se mantiene, en general, dentro de un rango comparable (con variaciones naturales), sin evidenciar una tendencia sostenida de aumento o disminución en la frecuencia de eventos de magnitud ≥ 4.5; los años recientes muestran valores en un rango similar a los del pasado. El año más reciente incluido en la base corresponde a un periodo parcial (aún no ha finalizado a la fecha de corte de los datos), por lo que su conteo no es directamente comparable con el de un año completo.