Sismicidad histórica

Uno de los problemas que se presenta en la evaluación de la amenaza sísmica en una región determinada es la falta de información que permita en un largo periodo, por un lado, determinar el sismo característico o evento máximo posible y, por otro lado, calcular el periodo de retorno o tiempo promedio en el cual podría presentarse dicho evento. Para resolver este problema, se requiere realizar estudios confiables de sismicidad histórica que conduzcan al mejor conocimiento de los terremotos que han ocurrido, al ampliar la ventana de observación sismológica al reciente histórico del país; entiéndase para Colombia, a partir de la conquista española a comienzos del siglo XVI.

Los estudios históricos tienen la finalidad de realizar la búsqueda de diversos tipos de información (observaciones de campo, relatos detectados en fuentes escritas como manuscritos, libros y prensa o, incluso, instrumental) sobre los daños y efectos causados por un terremoto en un área específica. Dependiendo de la historia sísmica que tenga la región donde se pretenda llevar a cabo un estudio sobre los sismos ocurridos en tiempos lejanos y del tipo de información a consultar, según ^{Gouin (1994)}, ^{Audemard (1998}, 2003, 2005) y ^{Salcedo-Hurtado (1999)}, los periodos de análisis para los estudios históricos sobre terremotos pueden ser divididos en tres tipos:

De esta manera, los datos históricos sobre terremotos aportan información sobre los daños ocasionados en diversas poblaciones, permitiendo obtener los datos cualitativos que son interpretados a la luz de las escalas macrosísmicas para evaluar el grado de intensidad en cada pueblo o localidad, conocido como centro o punto de intensidad; esto posibilita la elaboración de mapas de intensidades a partir de los cuales se trazarían, cuando se tienen puntos suficientes, las curvas de isosistas, que corresponden a curvas que encierran zonas de igual intensidad.

La investigación histórica, en especial la que corresponde al periodo histórico propiamente dicho, se desarrolla mediante la lectura sistemática de distintas fuentes, permitiendo ampliar el conocimiento sobre eventos símicos conocidos (investigación intensivista) o encontrar nuevos sismos no incluidos en los catálogos sismológicos (investigación extensivista) (^{Rodríguez de la Torre, 1993}). En estos casos, es muy importante realizar búsqueda en diversos tipos de espacios como archivos históricos, archivos locales (municipales o de ciudades y notariales), archivos de iglesias y monasterios, archivos de dominios aristocráticos.

La información que se recoja en estos archivos tiene la ventaja de ser un material que corresponde a relatos de testigos oculares, y puede dar cuenta sobre daños en edificaciones y diferentes tipos de efectos en diversas áreas. Sin embargo, también tiene la desventaja de que estos archivos frecuentemente están en desorden: los investigadores pueden encontrar información solamente estudiando el archivo completo o por casualidad; igualmente, la información puede ser muy limitada y a veces dudosa (^{Salcedo-Hurtado, 1999}), por lo cual se requiere de cierta experticia para la interpretación y para la adecuada asignación de la intensidad sísmica. En este sentido, Rodríguez y ^{Audemard (2003)} dan cuenta de un rosario de limitaciones que tiene esta fuente documental particular, con base en casos venezolanos estudiados por esta aproximación metodológica.

De acuerdo con ^{Salcedo-Hurtado (1999)}, un estudio histórico sobre los terremotos requiere de la búsqueda de fuentes primarias, que permitan tener registros detallados tomados de fuentes confiables para que sean interpretadas con fines sismológicos. ^{Stucchi (1994)} señala que, para poder ser utilizados con fines de interpretaciones sismológicas, los registros históricos necesitan ser ubicados en el espacio y el tiempo, además, interpretados en términos de la intensidad macrosísmica; tal operación produce las llamadas intensidades puntuales (o centros de intensidad), las cuales pueden ser asumidas como observaciones macrosísmicas. Dichos datos de intensidad alimentan los catálogos sísmicos de una región.

Clasificación y valoración de las fuentes bibliográficas

Un aspecto fundamental para la evaluación de intensidades macrosísmicas de terremotos históricos concierne a la valoración de las fuentes bibliográficas. ^{Izquierdo-Álvarez (1999)} considera que al valorar las fuentes de información se debe tener claro su tipo y finalidad (si se trata de un informe para pedir ayudas económicas o de un diario personal), fiabilidad, grado de originalidad (si la información es de primera mano o proviene de otra fuente); se debe ser cuidadoso para que la información que se extraiga de cada fuente sea lo más exacta posible, tratando de perder la menor cantidad de datos. Es importante tener mucho cuidado con la localización espacio-temporal, puesto que errores de esta clase pueden dar lugar a duplicaciones de eventos, a atribuir efectos a terremotos o lugares equivocados, entre otros aspectos.

Intensidad y mapa de intensidad sísmica

Existen diversos tipos de parámetros mediante los cuales se puede valorar la fuerza de un terremoto, entre ellos parámetros instrumentales y no instrumentales. Entre los instrumentales se encuentra la magnitud, que corresponde a una escala cuantitativa que inicialmente Charles Richter relacionó con la energía liberada por un sismo, que es independiente del lugar de observación (^{Lay y Wallace, 1995}; ^{Bolt, 2003}). Por su parte, entre los parámetros no instrumentales más representativos está la intensidad sísmica, representada como I.

La intensidad sísmica se entiende como una unidad de medida adimensional que clasifica la severidad o las manifestaciones de las vibraciones sísmicas en la superficie terrestre; se mide en grados y se especifica con un número romano previamente determinado por una escala específica. La medición de la intensidad sísmica no requiere de instrumentos para ser evaluada; sin embargo, se hace a través de ciertos sensores usados históricamente, que pueden ser divididos en cuatro grupos (^{Grünthal, 2009}):

El hecho de que la intensidad no sea medida instrumentalmente le da un carácter subjetivo a su valoración, la cual depende de cómo el observador interprete la información sobre los daños y efectos causados por terremotos, descritos acorde con los sensores indicados.

Existen diversas escalas para valorar el nivel de intensidad de un terremoto. En su mayoría, estas escalas estiman doce grados del I al XII, en los cuales I corresponde a vibraciones débiles que solo son notadas por los instrumentos, sin que sean sentidas por las personas; XII hace referencia a fuertes catástrofes, donde ninguna construcción se sostiene, hay grandes cambios en el relieve y muchas víctimas; desde el grado II hasta el XI se consideran como efectos intermedios. Entonces, se acostumbra designar los grados de intensidad con numeración que va desde I a XII, como es el caso de la Escala Macrosísmica Europea, que corresponde a una actualización de la escala MSK-64, realizada por la Comisión Europea en el año 1998 (^{Grünthal, 2009}).

La evaluación de la intensidad sísmica de un terremoto muestra que es mayor en la zona epicentral (área de la superficie terrestre donde se observa la mayor cantidad y más severos daños causados por el terremoto) y disminuye con la distancia a media que el observador se aleja de la fuente sísmica. La intensidad evaluada en las vecindades del epicentro se conoce como la intensidad epicentral y se presenta como I₀. A partir de este criterio, se pueden establecer modelos de atenuación de la intensidad, los cuales son fundamentales para la valoración de la amenaza y el riesgo sísmico.

La asignación del nivel de intensidad de un terremoto en diferentes sitios sobre la superficie terrestre conduce a la realización de los denominados mapas de intensidad. La elaboración de un mapa de intensidades conlleva a elegir sitios georreferenciados sobre la superficie de la Tierra, que han sido afectados por las ondas sísmicas generadas por un terremoto fuerte, y asociarles un valor de intensidad conforme a una escala macrosísmica elegida previamente.

Mapas de isosistas

Un mapa de isosistas consiste en el trazado de curvas que encierran áreas de igual nivel de intensidad sísmica en una región que ha sido afectada durante la propagación de las ondas sísmicas generadas por un terremoto. A medida que las ondas sísmicas se alejan del epicentro, se espera un decrecimiento de su intensidad y, dependiendo del tipo de terreno que atraviesen, su propagación se ve favorecida o retardada: hay mejor propagación en la dirección que transcurre la estructura de la fuente generadora del sismo, generalmente una falla geológica, y mayor atenuación en sentido perpendicular, razón por la cual las líneas de isosistas suelen tener forma de elipse relativamente deformada (^{Musson y Cecić, 2012}). En algunos casos estas curvas se simplifican y tienden a dibujarse de manera circular; el caso es que la realización de las isosistas implica el trazo de isolíneas que encierran iguales valores de intensidad.

Generalmente, el epicentro macrosísmico de un terremoto está rodeado por la isosista de mayor valor, la cual bordea la zona pleistosísmica respondiendo, con alguna aproximación, a la forma del foco, que aparece como la última proyección en la superficie de la tierra. Dependiendo de la rapidez con que se atenúen las ondas sísmicas, en la medida que aumenta la distancia epicentral, se configura la distribución y forma de las isosistas subsiguientes, las cuales reflejan claramente los campos de emisión, transmisión y recepción de un terremoto.

Las curvas de isosistas han sido utilizadas con diversos propósitos en el campo de la sismología y la ingeniería y para establecer las áreas de mayor intervención en la mitigación y atención de los efectos causados por fuertes terremotos. Por ejemplo, ^{Grases y Rodríguez (2001)} establecen correlaciones que permiten calcular la magnitud de un sismo histórico (no-instrumental) a partir del tamaño o superficie de las áreas de una intensidad dada. De igual manera, como lo muestran ^{Salcedo-Hurtado et al. (2021a)}, los mapas de isosistas también son usados para calcular diversos parámetros focales de los terremotos.

Información instrumental

La primera fuente instrumental que da cuenta de este evento sísmico es el catálogo elaborado por ^{Gutenberg y Richter (1954)}, que lo registra a las 23:41 UTC, localizado en el océano Pacífico en las coordenadas 3,0ºN de latitud y 78,0ºO de longitud, a 170 km de profundidad. Por su parte, ^{Engdahl y Villaseñor (2002)} establecen que el sismo tuvo su epicentro cerca del municipio de Tuluá, en el departamento del Valle del Cauca, con coordenadas de 4,02ºN de latitud y 76,07ºO de longitud, con profundidad superficial, convencionalmente establecida en 33 km. En el catálogo de ^{CERESIS (1985)} se transcriben los parámetros instrumentales inicialmente reportados por Gutenberg y Richter. El proyecto Global Earthquake Model (^{ISG-GEM, 2015}) se definen las coordenadas del epicentro en 4,279°N y 75,331°O, y se calcula una magnitud Mw 6,3.

Información y estudios macrosísmicos

El primer reporte macrosísmico de este terremoto aparece en el informe presentado por el Grupo de Estudio de Riesgo Sísmico del Suroccidente Colombiano (^{GERSCO, 1987}); estos datos son retomados y ampliados por ^{Espinosa-Baquero (2003)}. El trabajo de ^{Meyer (1990)} también retoma los datos macrosísmicos de GERSCO y propone que el epicentro se localiza al norte de la ciudad de Cali, con origen en la zona de Wadatti-Benioff (entiéndase sismo de subducción), cambiando la localización previamente dada por ^{Engdahl y Villaseñor (2002)}. Por su parte, ^{Mendoza et al. (2004)} reportan nuevos datos macrosísmicos que indican que los daños causados por este evento se concentran principalmente en poblaciones del departamento del Valle del Cauca. Las primeras hipótesis de intensidad macrosísmica son aportadas por ^{Espinosa-Baquero (2003)}. Otros autores como ^{INGEOMINAS (2005a}, 2010), ^{Salcedo-Hurtado et al. (2007)} y ^{Martínez (2014)} han aportado nuevos datos para mejorar el conocimiento de los parámetros sismológicos de este evento sísmico. Cabe señalar que ninguna de las versiones del libro Historia de los terremotos en Colombia del padre Jesús Emilio Ramírez (^{Ramírez, 1975}) incluye información sobre los efectos causados por este evento.

En la Figura 2 se muestra la secuencia temporal y la naturaleza de los estudios anteriormente realizados sobre el terremoto del 07 de junio de 1925, ocurrido en el suroccidente colombiano. Por su parte, la Tabla 1 presenta las intensidades asignadas por cada uno de los estudios anteriores, acorde con la información recopilada. Como se aprecia, por causa de este terremoto resultaron afectados diversos municipios de varios departamentos de las regiones del Pacífico, Eje Cafetero, centro y norte de Colombia; sin embargo, los mayores efectos se presentaron en la ciudad de Cali y en otras poblaciones del departamento del Valle del Cauca.

Se puede apreciar que en el caso del estudio de ^{Espinosa-Baquero (2003)}, se asigna intensidad de VIII a las poblaciones de Cali, La Cumbre, Restrepo y Yotoco. Los demás autores asignan un valor máximo de intensidad entre VII-VIII para la ciudad de Cali. Puede verse que no todas las poblaciones mencionadas en la Tabla 1 tienen asignaciones de intensidad, lo cual muestra un vacío de información macrosísmica que es complementado con la búsqueda de nueva información realizada en el presente estudio, como se muestra más adelante.

Figura 2 Estado del arte del terremoto del 07 de junio de 1925 (modificado de ^{Salcedo-Hurtado et al., 2007}). 

La Figura 3 presenta, de manera gráfica, la comparación de la asignación de las intensidades macrosísmicas dada por los estudios anteriores en diversas poblaciones del territorio colombiano, conforme se discrimina en la Tabla 1.

Figura 3 Intensidades macrosísmicas asignadas por estudios anteriores a diversas poblaciones afectadas por el terremoto del 07 de junio de 1925. 

Para establecer la relación entre los diversos valores de intensidad asignados en los estudios anteriores, se aplicó el coeficiente de correlación lineal de Pearson, expresado como (1):

donde σ _XY es la covarianza muestral de las variables (X, Y), σ _X es la desviación estándar muestral de la variable X y σ _Y , la desviación estándar muestral de la variable Y.

El coeficiente de correlación es adimensional y varía en el intervalo [−1, +1], donde −1 expresa una correlación máxima inversa, mientras +1 representa una correlación máxima directa; un valor de correlación 0 indica que no existe correlación alguna entre las variables consideradas. Este procedimiento ha sido usado por otros autores (^{Salcedo-Hurtado y Castaño-Castaño, 2011}; ^{Leal-Guzmán et al., 2018}). La Tabla 2 muestra el resultado de la correlación entre las asignaciones de intensidad para el terremoto del 07 de junio de 1925, dadas por los estudios macrosísmicos anteriores.

Como se puede apreciar, independientemente del tipo de escala macrosísmica utilizada en cada estudio, se nota muy baja correlación entre los valores de intensidad asignados; la mayor correlación se presenta en las asignaciones dadas por ^{INGEOMINAS (2010)} y ^{Martínez (2014)}, con un coeficiente del orden de 0,82. La correlación que existe entre la evaluación hecha por ^{Espinosa-Baquero (2003)} con las demás fuentes es muy baja. El estudio de ^{Espinosa-Baquero (2003)} se toma como referencia debido a que es el que presenta la mayor cantidad de poblaciones evaluadas. Sin embargo, la evaluación de las intensidades en dicho estudio se realizó con base en la escala de Mercalli modificada, que, en los últimos años, en Colombia, ha entrado en desuso.

Así, la reevaluación de las intensidades macrosísmicas de este terremoto se hace necesaria, por un lado, debido a que los estudios macrosísmicos anteriores no unifican la información y tampoco realizan la evaluación de intensidades para todas las poblaciones afectadas. Por otro lado, por la disparidad o falta de correlación de las asignaciones anteriores, lo que muestra gran incertidumbre entre ellas, como se indica en la Tabla 2.

Estos aspectos marcan la necesidad de hacer nuevas interpretaciones de la información y reevaluación de la intensidad macrosísmica para el terremoto del 07 de junio de 1925, preferiblemente usando una única escala de intensidad, en este caso la EMS-98, cuya aplicación, como lo muestran ^{Salcedo-Hurtado y Castaño-Castaño (2011)}, ^{Salcedo-Hurtado y Gómez- Capera (2013)} y ^{Gómez-Capera et al. (2014)}, presenta ventajas frente a resultados obtenidos con otras escalas macrosísmicas tradicionalmente usadas en Colombia.

El presente estudio aporta nuevos datos macrosísmicos: por un lado, reporta daños y efectos ocasionados a varias poblaciones no señaladas en los estudios anteriores y, por el otro, propone nueva hipótesis de intensidad, en cada una de las poblaciones de diferentes departamentos afectadas por este terremoto, dejando puertas abiertas para que en estudios macrosísmicos posteriores se evalúen los parámetros sismológicos como el epicentro, la profundidad y la magnitud, entre otros (^{Blake, 1941}; ^{Shebalin, 1968}, 1971; ^{Bakun y Wentworth, 1997}; ^{Salcedo-Hurtado et al., 2021a}).

Fuentes de información histórica

En este estudio, la búsqueda de información se ha llevado a cabo en archivos, bibliotecas, hemerotecas, librerías, centros de documentación de diferentes ciudades. Acorde con la información contenida, los documentos obtenidos y revisados fueron categorizados como fuentes primarias y secundarias, principalmente.

Las fuentes primarias corresponden a documentos originales como libros, revistas, periódicos, diarios, documentos oficiales de instituciones públicas, informes técnicos y de investigación, que contienen información primaria o contemporánea con la fecha del evento y que narran los hechos observados por testigos directos. Por su parte, las fuentes secundarias hacen alusión a publicaciones como libros y artículos, que dan cuenta del evento como producto del análisis, reorganización de datos e información contenida en otros trabajos primarios originales.

Captura y almacenamiento de la información

Para capturar y almacenar la información se utilizó una ficha bibliográfica (^{Salcedo-Hurtado, 1999}), donde se documentan los datos históricos sobre los efectos o daños ocasionados por el terremoto en cuestión. En esta ficha se hace la transcripción de la información encontrada en cada fuente bibliográfica. Esto se hace describiendo tanto la fuente donde se encontró la noticia como las características del evento sísmico, dentro de las cuales se detallan los daños o efectos tal como aparece en la noticia, sin modificar o corregir ningún apartado de la misma. Dado que la noticia sísmica aparece en diversas fuentes bibliográficas, se requiere sintetizar la información sobre los datos encontrados. Para ello, se elaboró una segunda ficha denominada “Ficha Resumen de Efectos del Terremoto” (FRET), en la cual se consignan los diferentes efectos con los cuales se evalúa la intensidad sísmica (^{Salcedo- Hurtado et al., 2007}).

Figura 4 Esquema de procedimientos para la actualización del conjunto de datos macrosísmicos (modificado de ^{Guidoboni y Stucchi, 1993}). 

Asignación de intensidades

La evaluación y asignación de intensidades se realiza mediante una escala macrosísmica que permita procesar la información sobre los efectos ocasionados por un terremoto. Para este estudio se empleó la Escala Macrosísmica Europea 1998 (EMS-98) que clasifica la intensidad del terremoto en 12 niveles, desde I (no sentido) hasta XII (^{Grünthal, 2009}).

La valoración de intensidad realizada en el presente estudio introduce un factor de confiabilidad de la información y certeza en la hipótesis de la intensidad asignada, lo cual se hace mediante un factor de calidad Q, descrito como: A - el sismo fue sentido y la información es suficiente y veraz para dar un valor apropiado de intensidad; B - el sismo fue sentido, pero la información no es suficiente para dar un valor de intensidad confiable; C - la información es insuficiente y dudosa (^{Salcedo-Hurtado, 1999}, 2002).

Elaboración del mapa de intensidades e isosistas

Una vez asignada la intensida a cada uno de los lugares señalados por las diversas fuentes bibliográficas ⸻noticias sísmicas encontradas sobre el evento bajo estudio donde se haya sentido el sismo o este haya causado algún tipo de efecto⸻, se procede a la elaboración del respectivo mapa de intensidades que indica el valor de intensidad asignada a cada población, acorde con la escala macrosísmica usada, en este caso EMS-98.

Las poblaciones afectadas han sido georreferenciadas acorde con la base cartográfica oficial para Colombia; con ello se realiza el trazado de las líneas de isosistas siguiendo una figura elíptica (^{Musson y Cecić, 2012}) y empleando sistemas de información geográfica, de tal manera que cada población y su respectivo valor o punto de intensidad quede ubicado dentro de un área que represente un valor de intensidad, encerrada por su correspondiente isolínea.

Es posible que algún punto de intensidad quede ubicado fuera de su respectiva área de intensidad y aparezca dentro de una isosista que encierre otra área con valor de intensidad diferente; por ejemplo, un punto de intensidad V puede quedar ubicado dentro del área de una isosista de intensidad IV o VI, esto se debe a las variaciones inducidas por efectos locales y por el correcto trazado de las líneas de isosistas, que buscan guardar la forma de una elipse (^{Musson y Cecić, 2012}). Es preciso aclarar que el espaciado de las isolíneas no necesariamente es homogéneo en todo el campo macrosísmico del respectivo terremoto.

Nueva investigación histórica

Como resultado de los procesos metodológicos descritos anteriormente, la reevaluación de la intensidad macrosísmica del terremoto del 07 de junio de 1925 permitió establecer los siguientes productos: en primer lugar, se clasificaron las fuentes bibliográficas encontradas que contienen información sobre el evento señalado; luego, se discriminaron los datos macrosísmicos descritos en cada fuente de información, clasificando los tipos y modos de daños y efectos causados en cada población afectada o en aquellas donde se reportó que el sismo al menos fue sentido; por último, se asignó o reevaluó el valor de intensidad en cada población, de acuerdo con la descripción de la Escala Macrosísmica Europea de 1998 (^{Grünthal, 2009}), y se elaboraron los mapas de intensidad e isosistas.

Clasificación de las fuentes de información

Discriminación de daños y efectos causados por el terremoto

A partir de la información extraída de cada una de las fuentes, se discriminan los daños y efectos en cada población mencionada, los cuales, según la escala EMS-98, se categorizan desde débiles hasta totalmente devastadores, en cuanto afecten personas, objetos, construcciones y el entorno natural (efectos sismogeológicos). Estos últimos corresponden a efectos complejos influenciados con frecuencia por distintos factores como la estabilidad inherente de una ladera, el nivel freático, etc. (^{Grünthal, 2009}).

Principales efectos en personas, objetos y construcciones: fueron diversos los efectos reportados por los habitantes de las poblaciones donde el sismo del 07 de junio de 1925 fue sentido o produjo daños. La información al respecto fue extraída de las diferentes fuentes bibliográficas descritas anteriormente. Es importante tener en cuenta que para la descripción de los efectos se hizo la transcripción exacta de cómo la población percibió y relató el acontecimiento, esto con el fin de realizar la asignación de la intensidad con la Escala Macrosísmica Europea (EMS-98). Los efectos en objetos se manifiestan en el sentido de que diversos elementos observados por las personas, se balancean, caen y, en muchos casos, se rompen a causa de la vibración sísmica. Al igual que los efectos en personas, para la determinación de intensidades, es necesario mostrar la información puntual de los efectos en objetos ocasionados por el evento sísmico. Por último, se tiene que las afectaciones en construcciones pueden impactar seriamente en el estado y la calidad de las viviendas, que pueden llegar a sufrir colapso parcial o incluso total de las mismas. Esto provocaría más consecuencias e incremento de la intensidad macrosísmica. En el Anexo 3 se muestra, de manera resumida, los efectos principales causados en las personas, objetos y construcciones en algunas de las poblaciones donde las diversas fuentes consultadas indican algún tipo de afectación provocada por el terremoto del 07 de junio de 1925. Se nota que los mayores efectos fueron reportados en la ciudad de Cali.

Efectos en el entorno natural: a pesar de que el evento sísmico del 07 de junio de 1925 ocasionó varios daños en la ciudad de Cali, no hubo reportes de efectos en el medio ambiente, tales como deslizamientos o eyección de lodo de manera súbita; por el contrario, el comportamiento de los suelos aparentemente fue estable, a excepción de la población de Piendamó en el departamento del Cauca, donde se reportaron efectos de este tipo (Tabla 3).

Reevaluación de intensidades macrosísmicas

De acuerdo con la información recopilada y la respectiva clasificación y discriminación de los efectos y daños ocasionados, se reevaluó la intensidad macrosísmica en cada una de las poblaciones donde el sismo del 07 de junio de 1925 fue sentido o dejó algún tipo de afectación; la Tabla 4 muestra el resultado de este proceso.

En la Figura 5 se presenta un diagrama de barras que indica el nivel de intensidad asignada a cada una de las poblaciones donde se evaluó la respectiva información macrosísmica. En esta figura se aprecia que, en relación con los estudios macrosísmicos anteriores sobre el terremoto del 07 de junio de 1925 (ver Figura 3), el presente trabajo complementa la información asignando la intensidad a cada una de las poblaciones afectadas según la información de las diferentes fuentes bibliográficas que se consultaron.

Los valores de intensidad presentados en la Tabla 4 y la Figura 5 muestran que a la ciudad de Santiago de Cali se le ha asignado el mayor nivel de intensidad, clasificada con daños severos (VIII), debido a que se dieron afectaciones serias e incluso colapsos parciales en muchas de las edificaciones presentes para la época del evento, lo que causó pánico extremo entre los habitantes de la ciudad.

Obsérvese que la información utilizada, aunque para algunas poblaciones puede ser poca, ha permitido dar asignaciones de intensidad de buena calidad, a pesar de que algunos valores pueden ser intermedios, como es el caso de las intensidades asignadas a los municipios de La Cumbre y La Unión en el departamento del Valle del Cauca, donde se reportan valores de VI-VII y V-VI, respectivamente. Debe quedar claro que estos no corresponden a valores decimales de 6,5 y 5,5, esto significa que la información reportada, por ejemplo, para el caso del municipio de La Cumbre, contiene algunos elementos que permiten asignar intensidad de VI, pero además indica datos del nivel VII que no son suficientes para que este sea el valor final de la intensidad en dicha población.

Figura 5 Histograma de asignación de intensidades macrosísmicas en diversas poblaciones afectadas por el terremoto del 07 de junio de 1925 en el suroccidente colombiano. 

Mapa de intensidad mejorado

En la Figura 6 se presenta el mapa de puntos de intensidad macrosísmica, donde se distribuyen las intensidades asignadas, en la escala EMS-98, a cada una de las 42 poblaciones que, de acuerdo con los datos reportados por las diferentes fuentes de información, han sido afectadas por dicho evento sísmico (ver Tabla 4). Puede observarse que estos puntos geográficamente se distribuyen principalmente en los departamentos del Valle del Cauca, Tolima, Cundinamarca, Chocó, Quindío y Cauca; también se tienen puntos de intensidad en los departamentos de Risaralda, Caldas, Antioquia y Boyacá. El Valle del Cauca es el departamento con mayor número de poblaciones afectadas, específicamente 13 municipios con reportes de algún tipo de afectación.

La distribución geográfica de estos puntos de intensidad marca una ligera tendencia SSO-NNE, donde los mayores valores de intensidad se asignan a las poblaciones del departamento del Valle del Cauca, con intensidades desde V hasta VIII, y el valor máximo se asigna a la ciudad de Cali, a partir de la cual se aprecia el decaimiento de la intensidad tanto hacia el NNE como hacia el SSO, marcando mayor elongación hacia el NNE. Los puntos de menor intensidad mostrados en la Tabla 4, asignados a las ciudades de Barranquilla y Cúcuta con valores de II y III, respectivamente, por razones de la escala del mapa, no se muestran en la Figura 6.

Mapa de isosistas y epicentro macrosísmico

Las isosistas, en esencia, son la representación gráfica del campo macrosísmico de un terremoto dado a partir de la irradiación de la energía sísmica desde el foco. En el presente trabajo se trazan las curvas de isosistas correspondientes al campo macrosísmico del terremoto del 07 de junio de 1925 ocurrido en el suroccidente colombiano (Figura 7). En este mapa se nota que la isolínea de mayor valor corresponde a la zona afectada con intensidad VIII, que encierra la ciudad de Santiago de Cali; las demás isolíneas varían en intensidades entre IV y VII en la escala EMS-98. Las isolíneas dibujadas tienen forma elíptica mostrando que su eje mayor tiene una tendencia SSO-NNE, dirección en la cual se supone la mejor propagación de la energía elástica. Mientras que el eje menor aparece en la dirección ONO-ESE, en la cual se tiene mayor atenuación de la energía elástica.

En este caso, se evalúa el epicentro macrosísmico del terremoto del 07 de junio de 1925, el cual preliminarmente se ha ubicado en las coordenadas: 3,45° norte y 76,49° oeste, en cercanías de la ciudad de Cali, donde se evaluó la intensidad máxima.

Figura 6 Mapa de centros de intensidad macrosísmica del terremoto del 07 de junio de 1925. 

Varios autores y metodologías hacen uso de los mapas de isosistas para evaluar diversos parámetros macrosísmicos de los terremotos, como es el caso de leyes de atenuación, profundidad focal, magnitud, entre otros (^{Barbano et al., 1979}; ^{Bottari et al., 1982a}, 1982b, 1986; ^{Shebalin, 1971}, 1974, 1991; ^{Grases y Rodríguez, 2001}; ^{Marín-Arias et al., 2006}, 2009; ^{Salcedo-Hurtado et al., 2021a}).