La necesidad y la oportunidad de estudiar el impacto económico de los tsunamis en las costas españolas

En los últimos 20 años los tsunamis han sido responsables de cientos de miles de víctimas mortales directas y más de 11 millones de desplazados (1), lo que convierte a este tipo de fenómenos en los más mortíferos de entre los peligros naturales de la historia de la humanidad. No todas las zonas del mundo son propensas a sufrir tsunamis, solo lo son las regiones próximas a una masa de agua importante, o conectadas a ella por medio de algún accidente geológico; y entre las zonas donde puede generarse un tsunami, la probabilidad de que ocurra es diferente en cada una de ellas. Más de la mitad de los registros de tsunamis (2) tienen como punto de origen el Océano Pacífico (figura 1), aunque es más acertado decir que de todos los registros de tsunamis el 60% aproximadamente se originan en el denominado Anillo de Fuego, la zona próxima a las costas de Asia y América, donde las placas oceánicas colisionan con las continentales y generan una intensa actividad sísmica y volcánica. Llama la atención que, si bien en el Océano Atlántico (incluyendo el Mar Caribe) se acumula un 10% de los eventos, en el Mediterráneo la cifra es del orden del 13%. Es llamativo porque, puesto en perspectiva, el Mediterráneo es unas 70 veces más pequeño que el Océano Pacífico en términos de área, por lo que, por unidad de superficie, el Mar Mediterráneo es una zona en apariencia más tsunamigénica que el Anillo de Fuego. Se puede plantear como hipótesis que en el Mediterráneo el registro de tsunamis es más completo que en otras partes del globo, por razones relacionadas con los asentamientos humanos, su nivel de desarrollo y su capacidad para registrarlos. También se puede plantear como hipótesis que ni siquiera el catálogo de eventos en el Mediterráneo sea completo, ya que es creíble que hayan ocurrido tsunamis que hayan pasado inadvertidos o cuya huella geológica aún no se haya descubierto. A pesar de estas dos hipótesis de no-completitud del registro, se pueden obtener algunas conclusiones preliminares de gran interés. Por ejemplo, y con efecto sobre las costas de la Península Ibérica, la probabilidad de que ocurra un tsunami en el Mediterráneo es del 14% cada año, mientras que en el Atlántico es del 5% cada año, o tal vez sea aún más frecuente si se confirmaran las hipótesis de no-completitud. La IOC-UNESCO (3) es más contundente en su forma de expresar la probabilidad, si cabe, pues señala que para un intervalo promedio de 30 años, es casi seguro que se produzca un tsunami significativo en el Mediterráneo.

Se estima que en las zonas costeras habita más del 40% de la población mundial (4), pero como se ha visto en eventos de tsunami recientes (Sumatra, 2004; Japón, 2011; Indonesia, 2018), los afectados no son solo los habitantes de la zona, sino también los turistas allí donde esta actividad tiene gran desarrollo asociado a la costa. Tanto es así, que en el tercer aniversario del Día Mundial para la Concienciación sobre los Tsunamis, se ha puesto parte del foco en esta particularidad del riesgo de tsunamis. 

En España (5) se estimó para el año 2017 que en las zonas de playa en particular, las pernoctaciones alcanzaron su máximo en unas 40.000 personas al día. Tomando como referencia el tsunami de Lisboa de 1755, se puede hacer un cálculo aproximado del número de españoles potencialmente víctimas de un tsunami equivalente a aquél que pudiera ocurrir hoy. Así, en 1755, en España, perecieron unas 1.214 personas (6). Por aquel entonces la población era de unos 9 a 10 millones de habitantes (7) y dado que el crecimiento poblacional hasta la actualidad ha sido cercano a x5 (o incluso x20 en zonas costeras), la cifra orientativa sería entre 5.000 y 24.000 personas. Es decir, que en términos de víctimas mortales, la magnitud de referencia para un evento similar al de 1755 podría presuponerse teóricamente en hasta unas 60.000 personas, incluyendo locales y turistas. En términos de afectados, la cifra ha de ser mucho mayor. En todo caso, esta aproximación debe ser revisada y trabajada para ofrecer una cifra apoyada en los métodos más avanzados, y usando el mayor número de datos y de la mejor calidad que resulte posible.

Para entender qué es un tsunami hay que empezar por el concepto más sencillo, el de una ola. En un estanque profundo de aguas tranquilas, si se lanza una piedra se forman en el agua unas ondas u olas. La distancia de “una ola” a otra suele ser bastante pequeña en este ejemplo, del orden de unos pocos centímetros, y la altura de las olas es aún menor. A la distancia entre dos crestas sucesivas la llamamos longitud de onda; al tiempo que tardan dos crestas en pasar por un mismo punto lo llamamos periodo; por último, a la mitad de la altura entre la cresta y el valle, lo llamamos amplitud (figura 2). Si se pudiera seguir con un rotulador a las moléculas de agua en esas olas, cerca de la superficie se dibujaría un círculo y, a medida que se va profundizando, se dibujarían círculos de menor radio cada vez. Hay una profundidad a partir de la cual el agua del estanque ya no se movería, esta profundidad es igual a la mitad de la longitud de onda. El ejercicio teórico del rotulador y las moléculas de agua serviría también para observar otra cosa: el desplazamiento de una ola ideal no transporta agua. Por eso, llamamos a este tipo de fenómenos perturbaciones oscilatorias. Si dos personas sujetan una cuerda, cada una por un extremo, y la agitan formando ondas, se comprobará que, en efecto, no hay transporte de masa, pues ninguno gana ni pierde cuerda, sólo la agitan cada uno por su extremo. Es decir, sólo hay transporte de energía. La energía de esas olas se describe mejor por medio de su longitud de onda que a través de su altura, pues a mayor longitud de onda, el volumen perturbado será mayor.

En el mar, las olas están provocadas principalmente por el viento y por el efecto de la atracción gravitatoria de los cuerpos cercanos, el sistema Sol-Tierra-Luna. En números redondos, las olas del Océano Atlántico tienen una longitud de onda de unos 90 metros y, en el Pacífico, de unos 300 m. Pues bien, un tsunami, a diferencia de otros tipos de perturbaciones oscilatorias en masas de agua, tiene una longitud de onda del orden de las decenas de kilómetros. Mientras que una ola común tiene un periodo del orden de 15 segundos, un tsunami tiene periodos típicos del orden de 10 a 30 minutos (8) , o incluso más. Dado que la profundidad media de los océanos es de unos 3.000 metros, un tsunami afecta o moviliza toda la columna de agua del océano, o casi toda. 

A medida que una ola se acerca a la costa se sobre-eleva, rompe y después progresa como una lámina hacia el interior (figura 2) aunque todo ello está en función, además de las características de la ola, de la geometría de la costa, incluyendo la parte emergida y la sumergida. El hecho de que las olas se deformen y rompan proviene de la imposibilidad del agua de seguir un patrón circular cuando la profundidad es menor que la mitad de longitud de onda. Tanto es así, que llega un momento en el que la oscilación se puede simplificar a un movimiento paralelo al fondo, en cizalla. Un tsunami, por este motivo, es más parecido a una zona de swash (o de esparcimiento) que a una ola propiamente dicha, y su avance tierra adentro sigue una oscilación (más o menos) paralela al relieve (si ignoramos la turbulencia). Cuando predomina el movimiento del agua hacia tierra, se dice que el tsunami está en avance (o en run-up en inglés) y cuando predomina el movimiento del agua hacia el mar se dice que el tsunami está en retroceso, en resaca (o en backwash en inglés). Mientras que en avance el proceso es decelerado (en contra de la gravedad y sobre una superficie no inundada y, por tanto, de alta fricción), en resaca el proceso es acelerado (a favor de la gravedad y moviéndose sobre una superficie inundada); por lo que un tsunami en resaca es igual o más peligroso que en avance.

Mar adentro, en zonas profundas, la altura de un tsunami es muy pequeña, de unos pocos centímetros, y se desplaza muy rápido, comparable a la velocidad de un vuelo comercial (8) (unos 900 km/h). Sin embargo, a medida que se acerca a la costa, por este efecto de no poder desarrollar el movimiento oscilatorio, la altura del tsunami se va haciendo cada vez mayor. La perspectiva de un tsunami desde tierra es parecida a un aumento súbito del nivel del mar, que progresa hacia tierra a una velocidad parecida a la de un caballo a galope (8) (casi 40 km/h). En definitiva, un tsunami puede provocar inundaciones costeras cuyo poder destructivo radica tanto en la velocidad con la que se desplaza el agua como en el volumen de agua y el número de veces que se repite, ya que no es “una ola” sino “un tren de olas” o una sucesión de olas. A esto hay que añadirle que a medida que el flujo avanza se va enriqueciendo en sedimentos (rocas, arena) y todo tipo de objetos (coches, residuos, troncos y restos vegetales) que se añaden a su capacidad para causar daños.

A estas alturas llega la pregunta de qué provoca un tsunami. En la base de datos de la NOAA (2) ya se apuntan algunos de los procesos que con mayor asiduidad dan lugar a tsunamis, que en orden de frecuencia serían: los terremotos, los movimientos de ladera y las erupciones volcánicas. Pero en realidad, cualquier fenómeno que tenga capacidad para desplazar un gran volumen de agua en poco tiempo será capaz de generar un tsunami.

Los terremotos son un fenómeno muy habitual, de hecho, están ocurriendo constantemente, pero sólo detectamos los que están en el rango de medición de los instrumentos. En promedio, al año se registran en el mundo más de un millón de terremotos de magnitud Mw (9) superior a 2. A medida que la magnitud es mayor, la frecuencia decrece rápidamente. Cada año se registran, de media, unos 20 terremotos 7 Mw o mayores (figura 3). 
 

Para que un terremoto provoque un tsunami, es necesario que se manifieste deformación con componente vertical en la superficie del fondo de la masa de agua (figura 4). Los terremotos en zonas de fallas convergentes o divergentes (mecanismo focal de falla inversa o normal) son más propensos a dar lugar a tsunamis que aquellos que se producen en zonas de fallas transformantes (o transcurrentes), ya que un mecanismo de falla inverso o normal implica una deformación mayoritaria en la componente vertical, mientras que las fallas transformantes implican deformación sobre todo en la componente horizontal. No obstante, como la deformación de las rocas por un terremoto es finita, los terremotos profundos (con hipocentro a más de 100 km) no pueden dar lugar a tsunamis, aun siendo de gran magnitud, pues la deformación no llega a ser significativa en el fondo de la masa de agua.

Dado que la profundidad es un aspecto tal vez más crítico que el mecanismo focal, un mapa de la distribución de terremotos por profundidad permite dar una idea preliminar de las zonas más propensas a generar tsunamis desencadenados por terremotos, que delimitan bastante bien los contornos de las placas tectónicas. Según el catálogo del Centro Sismológico Internacional, la gran mayoría de terremotos se producen a una profundidad inferior a 50 km (figura 5) y éstos ya sí podrían ser generadores de tsunamis si su magnitud fuera importante. No obstante, la mayoría de los terremotos más superficiales son de mecanismo predominantemente transformante. Por este motivo los tsunamis de origen sísmico no son tan frecuentes como los terremotos superficiales.

En principio, la idea de que para generar un tsunami es necesario deformación de la corteza con componente vertical suena bastante obvia. Un experimento casero rudimentario demuestra que una variación vertical brusca en el fondo de un cubo con agua genera olas en la superficie. Sin embargo, un terremoto no sólo puede deformar la corteza por un salto vertical, también puede desplazar el fondo oceánico una falla transformante si el relieve es lo bastante complejo y abrupto como para que un frente vertical se vea desplazado. Es la primera de las dos hipótesis que a día de hoy se barajan como posible causa del reciente tsunami de Indonesia (13) del 28 de septiembre de 2018. La segunda hipótesis que se plantea para la generación de este tsunami (con unas 2.000 víctimas y unos 5.000 desaparecidos) es que el terremoto haya desencadenado un deslizamiento submarino y éste, a su vez, haya generado el tsunami13. Dos casos claros de estudio de deformación de la corteza con salto vertical son el tsunami de Sumatra del 26 de diciembre de 2004 (con varios centenares de miles de víctimas) y el terremoto de Japón del 11 de marzo de 2011 (con decenas de miles de víctimas). Ambos con mecanismo focal de falla inversa. Este tipo de ruptura, la inversa o compresiva, es la que puede dar lugar a tsunamis con mayor facilidad, ya que la energía que pueden acumular las rocas en compresión es mucho mayor que la energía que acumulan en tracción y, por tanto, a igualdad de otros factores, la extensión de la zona de deformación suele ser mucho mayor que en los mecanismos normales o transformantes.

Durante el encuentro internacional de expertos en fuentes tsunamigénicas, celebrado en la Universidad de Málaga los días 6 y 7 de noviembre de 2017 (14), se identificaron tres zonas probables y capaces de generar tsunamis con potencial impacto significativo en las costas españolas. Se trata del área del golfo de Cádiz sensu lato (donde se originó el terremoto y tsunami de Lisboa de 1755), el Mar de Alborán (donde hay abundante actividad sísmica) y el margen norte de África (de donde se originó el primer registro de tsunami sismogénico el 21 de mayo de 2003 con daños económicos, sobre todo en las Baleares, pagados por el Consorcio de Compensación de Seguros por un monto de unos 350.000 euros (15).

Los movimientos de ladera ocurren a modo de corolario de toda una secuencia de procesos geológicos. La fuerza que actúa en última instancia es la gravedad, ocasionando que una parte de una ladera se movilice pendiente abajo hasta una nueva situación de equilibrio. En laderas emergidas a veces es suficiente un cambio de humedad, el tránsito de un animal, el viento, la lluvia o un terremoto. La relación entre la magnitud de los terremotos y el área afectada por movimientos de ladera (AAML) es conocida desde hace décadas (16) y se está progresando hacia conceptos más depurados como, por ejemplo, la relación entre la aceleración pico del suelo y el AAML (17). Sea como fuere, que una masa de roca caiga en un cuerpo de agua puede provocar una gran ola. De hecho, el mayor tsunami jamás registrado en términos de altura alcanzada (el 9 de julio de 1958) fue provocado por un movimiento de ladera que terminó en la bahía de Lituya, en Alaska, generando una ola cuyo remonte (run-up) fue de más de 520 metros (18). La simulación de este tipo de mega-tsunamis representa un reto para la los modelos matemáticos y la única simulación numérica realizada con éxito, hasta el momento, en la geometría real de la Bahía de Lituya, corresponde al grupo de investigación en Ecuaciones Diferenciales, Análisis Numérico y Aplicaciones de la Universidad de Málaga (grupo EDANYA (19)). Otro ejemplo imposible de obviar en toda introducción a los tsunamis es el de la presa de Vajont, en Italia, el 9 de octubre de 1963. La ladera que deslizó junto al embalse de Vajont, entró en la masa de agua a unos 100 km/h, generando una ola que sobrepasó la represa (de 262 m de alto) liberando 30 millones de metros cúbicos cauce abajo (20). Se cobró cerca de 2.000 víctimas.

A diferencia de los tsunamis provocados por terremotos, los provocados por movimientos de ladera afectan a una extensión menor y se consideran, por tanto, tsunamis locales. Este criterio puede ser utilizado para discernir, al menos en parte, el origen de un tsunami cuando su fuente es submarina. Sin embargo, dado que la correlación entre terremotos y movimientos del terreno en zonas emergidas es conocida, cabe pensar que exista también una correlación entre los terremotos y los deslizamientos submarinos, por lo que la génesis de muchos tsunamis puede ser una intrincada combinación de ambos efectos. Otros fenómenos de naturaleza co-sísmica, como la licuefacción de suelos, podrían también constituir desencadenantes de movimientos de ladera, incluidos los submarinos.

En realidad, no fue hasta el evento de Papúa Nueva Guinea del 17 de julio de 1998 (un terremoto y deslizamiento submarino) cuando se comenzó a considerar seriamente el efecto de los deslizamientos submarinos como causa significativa de tsunamis (21). Este tsunami se cobró más de 2.000 víctimas. Como consecuencia, se desarrollaron campañas de investigación batimétricas para describir mejor el fondo marino y los depósitos del deslizamiento, mejorando desde entonces el reconocimiento y la cartografía de otros eventos similares. De hecho, la hipótesis de un gran deslizamiento asociado al flanco volcánico de la isla de La Palma, en el archipiélago canario, donde existen registros de multitud de mega-deslizamientos (22), tuvo una gran repercusión mediática. Sin embargo, el estudio de este tipo de deslizamientos apunta a que su frecuencia es del orden de las decenas de miles de años, relacionados con los procesos volcánicos de desarrollo de los edificios volcánicos.

En España (y en otras regiones volcánicas próximas) la probabilidad de que se desarrolle un mega-deslizamiento es muy pequeña, por lo que tiene más interés desde el punto de vista científico que desde la perspectiva de planificación. Al menos, así se concluyó durante la reunión internacional de expertos de 2017 celebrada en la Universidad de Málaga (14). Durante este encuentro también se concluyó que la posible ocurrencia de deslizamientos más pequeños y próximos a núcleos de población sí podría generar impactos localmente significativos. No obstante, las laderas de los volcanes italianos Etna y Stromboli llevan décadas siendo monitorizadas con el doble propósito de vigilancia volcánica y de mega deslizamiento vinculado también a la ocurrencia de tsunamis (23).

Los procesos volcánicos incluyen muchos tipos de fenómenos diferentes y la mayoría de ellos implican o pueden implicar una movilización de un gran volumen de material. Entre estos procesos están los deslizamientos, las explosiones o los flujos de piroclastos. El evento volcánico y ulterior tsunami más conocido tal vez sea el de la isla de Krakatoa, que el 27 de agosto de 1883 culminaba su erupción con una explosión generalizada del edificio volcánico. Gran parte de la isla desapareció como consecuencia de esa explosión. 296 ciudades fueron destruidas y perecieron unas 36.000 personas, más del 90% como consecuencia del tsunami (24). La energía que se calcula que liberó esta explosión se compara con unas 10.000 bombas atómicas como la arrojada sobre Hiroshima (10). En comparación con los tsunamis provocados por terremotos, los provocados por volcanes representan del orden del 5% en los registros actuales (2). No obstante, de nuevo se puede plantear la no-completitud de los registros, pues es muy complicado realizar observaciones de este tipo de fenómenos. Existen unos 60 registros de tsunamis provocados por volcanes (2) durante el siglo XX. Algunos autores (24) reconocen hasta ocho mecanismos por los que un volcán puede provocar un tsunami, a saber: explosiones submarinas; flujos piroclásticos; lahares; terremotos de origen volcánico; movimientos de ladera; colapso de frentes de lava en la zona costera; colapso por subsidencia de calderas volcánicas; y ondas expansivas o de choque por explosiones en áreas emergidas. De ellos, parece ser que el más frecuente sería el provocado por explosiones submarinas.

En España, las Islas Canarias son la única región volcánicamente activa, pero no existe hasta la fecha ningún estudio que permita establecer un vínculo creíble (o útil a efectos de planificación) entre la probabilidad de erupción y la probabilidad de que ésta generase un tsunami; al menos más allá del colapso de los flancos mencionado en el apartado de movimientos de ladera (caso de la Isla de La Palma). Por otra parte, otros volcanes que podrían generar tsunamis (aunque es discutible que llegaran a afectar a las costas españolas) son el Etna, las Azores, el Stromboli (que el 30 de diciembre de 2002 generó un evento de este tipo (23) u otros aún más lejanos. 

En el archipiélago Balear, las inundaciones costeras provocadas por fenómenos atmosféricos son relativamente comunes y por ello, reciben un nombre propio: rissagas. En la prensa y en la literatura científica se han rebautizado como meteotsunamis, un término que se está popularizando. La última rissaga que afectó a las costas españolas ocurrió el pasado 15 de julio de 2018 y golpeó sobre todo al puerto de Alcudia. Hasta mediados de 1930, se pensaba que este tipo de oscilaciones extraordinarias (de mayor rango que las marejadas) respondían a configuraciones astronómicas, al igual que las mareas. Surge entonces un trabajo (25) que correlaciona los cambios de presión atmosférica con la generación de estos ascensos o descensos del nivel del mar, y que se ha ido depurando después ampliamente en la literatura científica. El fenómeno está actualmente descrito como una respuesta de amplificación de las oscilaciones en el mar y las oscilaciones atmosféricas (26), que genera perturbaciones que, a escala temporal y espacial, son mucho más cercanas a las de un tsunami que a las de otro tipo de fenómenos costeros con potencial para causar daños (como tormentas u ondas estacionarias o seiches (27)). Es decir, se trata de un acoplamiento de una oscilación de fuerte gradiente barométrico en la atmósfera con las ondas en la superficie del mar de origen gravitatorio o provocado por fuertes vientos.

No es posible abordar los tsunamis sin, al menos, hacer una mención a este particular desencadenante. En primer lugar, hay que tener en consideración que de recibir el impacto de un meteorito, lo más probable es que ocurra en alguno de nuestros océanos, pues ocupan el 70% de la superficie planetaria. De manera similar a un desprendimiento de roca, el impacto de un meteorito sobre una masa de agua puede provocar un tsunami, con la diferencia nada desdeñable de uno respecto a otro de la hiper-velocidad con la que llegan los meteoritos a la Tierra. Si los tsunamis de origen volcánico plantean dificultades en su observación, este tipo de eventos es aún más difícil de describir o caracterizar; de hecho, aún no hay ningún registro histórico fehaciente del tándem meteorito-tsunami, pero se presupone que el evento de impacto meteorítico al que se atribuye la extinción masiva de hace unos 65 millones de años, pudo haber generado un gran tsunami. La simulación de este tipo de eventos se suele hacer de manera similar al de una explosión puntual (28), lo que permite aproximar mejor las ondas de choque que provocaría la llegada de un meteorito a hipervelocidad. Por ahora, estos ejercicios tienen interés científico, pero haría falta mejorar nuestra capacidad de detección de cuerpos celestes para poder dar una cifra creíble sobre la probabilidad de este tipo de impactos.

Todo estudio de peligros naturales se hace con un único objetivo: minimizar los impactos indeseados de su ocurrencia (figura 6). Dado que la predicción de muchos eventos naturales no es posible, se recurre a diferentes técnicas o estrategias para su evaluación, que pueden incluir el uso de datos históricos o geológicos con traslado, tal cual, a la situación actual o mediante herramientas de generación de escenarios con estimación de probabilidad de ocurrencia, lo que implica necesariamente el uso de técnicas de simulación.

A la hora de evaluar el impacto de un peligro natural es necesario establecer algunas hipótesis de trabajo y, en particular, resulta útil asociar la probabilidad del desencadenante a la probabilidad de la ocurrencia del peligro natural. Con ello se desprecia la variabilidad en los resultados debidos a cambios en otras variables y, por ello, se requiere una actualización continua de este tipo de estudios. Dado que la cantidad de desencadenantes de un proceso natural y, en particular de los tsunamis, puede ser muy variada, cualquier aproximación centrada en un único desencadenante infravalora el peligro real de la ocurrencia del fenómeno. No obstante, con los datos disponibles hoy en día se puede aproximar bastante si se centra el caso de estudio en la causa más probable (o más probada), que para los tsunamis son los terremotos. Otra forma de abordar un peligro natural es simular un evento previo conocido, partiendo de la hipótesis de que “si ha ocurrido en el pasado, podría volver a ocurrir en el futuro”. En sentido estricto esta hipótesis es buena si los procesos son cíclicos (como la evaporación y precipitación de agua) pero para procesos no-cíclicos (como los movimientos de ladera o los terremotos) la hipótesis pierde validez en favor de una aproximación más probabilista. Incluso en casos de inundación por precipitación extrema, si bien el desencadenante es un proceso cíclico, no lo es su resultado, la inundación, pues los cambios morfológicos del cauce y los cambios que ocurren a nivel de cuenca dan como resultado inundaciones diferentes para precipitaciones idénticas en distintos momentos. Es por ello que la validez de este tipo de estudios se suele acotar a un futuro próximo, que para las inundaciones fluviales es de seis años según la Directiva Europea de Inundaciones (29).

La dificultad de trabajar con terremotos como desencadenantes de tsunami, es decir, con epicentro en una masa de agua, radica en que sus características no se suelen estimar con el mismo nivel de incertidumbre que los terremotos continentales. Esto se debe a la distribución de los sismógrafos, ubicados en su inmensa mayoría en tierra por razones obvias del coste de mantenimiento y a que se puede caracterizar mejor un terremoto con epicentro en terreno emergido por medio de estudios complementarios (geológicos e incluso históricos en proximidad con poblaciones). La falta de estudios de geología marina enfocados a la caracterización de las fuentes sísmicas tsunamigénicas es por tanto un reto que debe abordarse para reducir la incertidumbre. Entre tanto, muchas de las fuentes se caracterizan de acuerdo al estado de conocimiento actual, que se centra sobre todo en los registros sísmicos instrumentales y en la información derivada de las campañas oceanográficas, que a estos efectos aportan información de tipo batimétrico y geofísico. Con esta información se delinean y caracterizan las propiedades de las fuentes y se calcula su probabilidad de ruptura asociada a la magnitud y a la profundidad epicentral.

Dado que la mayor incertidumbre en el estudio de tsunamis es su fuente y debido a que es muy caro el estudio de las fuentes sismo-tsunamigénicas marinas (que implica trabajos oceanográficos intensivos), reducir las incertidumbres que sí se pueden controlar adquiere especial relevancia. Por ello, a nivel mundial existen iniciativas de benchmarking, una especie de competición de validación de casos de estudio, tanto reales como de laboratorio, para poner a punto los modelos de cálculo. Estos ensayos los lidera y dirige la NOAA en su Programa Nacional de Mitigación del Peligro de Tsunamis (más conocido por sus siglas, NTHMP del inglés National Tsunami Hazard Mitigation Program) que establece los requisitos o estándares a cumplir.

La fuente tsunamigénica, cuando es de origen sísmico, se suele aproximar siguiendo el modelo de Okada (30), que se apoya en la teoría de elasticidad para dar lugar a la deformación cortical que, en primera instancia, da lugar al tsunami. 

A partir de este modelo se deriva la deformación cortical que da lugar a la “condición inicial”, es decir, a la situación de deformación de la superficie libre del mar que genera el tsunami. De aquí en adelante queda simular la evolución de la onda inicial, primero en mar abierto para, finalmente, recrear su impacto en la costa.

Simular flujos es, posiblemente, uno de los mayores retos de la ciencia moderna, a pesar de que las ecuaciones que describen su movimiento son de principios del siglo XIX. A estas ecuaciones se las conoce como “ecuaciones de Navier-Stokes” (ENS) en honor a sus padres, el francés Claude-Louis Henri Navier y el irlandés George Gabriel Stokes. Se derivan de tres de los más importantes pilares de la ciencia, la Ley de Conservación de la Masa (Ley de Lomonosov-Lavoisier) la Ley de Conservación del Momento (Segunda Ley de Newton) y la Ley de Conservación de la Energía (Primera Ley de la Termodinámica). Las ENS no tienen solución unívoca conocida a día de hoy (salvo para casos particulares y sencillos) y existe toda una rama dedicada a su aproximación mediante métodos numéricos con la ayuda de computadores, conocida como CFD (Computer Fluid Dynamics, o dinámica de fluidos computacional) (31). Una de las aproximaciones (simplificaciones) más socorrida para simular agua (y otros fluidos newtonianos) es utilizar las ecuaciones “no lineales de aguas someras” (en inglés Non-Linear Shallow Water Equations, NLSWE) o Ecuaciones de Saint-Venant (en honor a su desarrollador, el francés Jean Claude Barré de Saint-Venant).

Para el caso de flujos encauzados (como el agua de un canal) se puede simplificar el movimiento del agua considerando que, para una sección en particular del cauce, el flujo sigue un único vector de dirección (simulación uni-dimensional). Para el caso de los tsunamis, la aproximación puede ser unidimensional (en perfiles) pero se desprecian los efectos relacionados con la interacción de las ondas durante su aproximación a la costa y debidas a su geometría, así como los efectos anidados de las reflexiones y difracciones de las ondas. Con ello se elimina toda posibilidad de evaluar efectos de amplificación costeros, que son críticos para explicar muchas circunstancias de un tsunami. A pesar de ello, la aproximación 1D es muy apropiada para abordar casos de estudio generalizados (o donde no sea posible realizar otro tipo de aproximaciones) siempre que se realice por especialistas de primer orden aportando resultados confiables (32), que pueden llegar a ser contrastados con modelos físicos reducidos a escala. Las aproximaciones bi-dimensionales consideran que, para un punto cualquiera del dominio de cálculo, el agua puede desplazarse en cualquier dirección de un plano. La diferencia en términos computacionales de un modelo 1D a un modelo 2D es de varios órdenes de magnitud de diferencia, tanto en demanda de datos como en tiempo de cálculo. En un ordenador convencional, una simulación puede tardar semanas en completarse dependiendo, sobre todo, del tamaño y resolución del dominio de cálculo. La simplificación a un modelo 3D, incrementa a su vez la complejidad otro tanto, y a día de hoy sólo se utiliza para casos de mucho detalle (en ningún caso a escala general o de cuenca) y, sobre todo, en el ámbito industrial (turbinas, motores, pilares de puentes, presas) a un coste computacional muy alto. 

La literatura científica está repleta de ensayos de simulación de tsunamis, y no ha sido hasta tiempo reciente que la computación empieza a abandonar, de forma lenta y progresiva, el cálculo en CPU (Central Processing Unit, o unidad central de procesamiento) para ejecutarse en GPU (Graphics Processing Unit, o unidad de procesamiento gráfico) ganando con ello velocidades de cálculo en el orden de 100-200 veces más rápidas (debido a que la arquitectura de las GPUs permite realizar miles  de operaciones sencillas simultáneamente, lo que hace posible la alta paralelización de determinados algoritmos de cálculo).

Uno de los modelos de simulación de tsunamis más avanzados, en cuanto a algoritmos precisos y computacionalmente eficientes, ha sido desarrollado en España por el grupo EDANYA (33) de la Universidad de Málaga. Este equipo de matemáticos ha estudiado e implementado decenas de esquemas numéricos, comparado sus resultados y tiempos de cálculo hasta seleccionar el código óptimo para la simulación numérica de tsunamis en el contexto de los Sistemas de Alerta Temprana de Tsunamis. Se trata del modelo de simulación Tsunami-HySEA, englobado en la familia de códigos HySEA (34) (Hyperbolic Systems and Efficent Algorithms o sistemas hiperbólicos y algoritmos eficientes), desarrollado en lenguaje CUDA para GPUs NVIDIA. La eficiencia de cálculo para resolver la propagación ha permitido reducir los tiempos de cálculo desde los días (o semanas) de los medios convencionales, a sólo algunos minutos e incluso, en la actualidad, a segundos con el uso de múltiples tarjetas NVIDIA de última generación (Teslas V100). Tanto es así, que este código ya ha sido integrado como módulo de cálculo en Sistema de Alerta Temprana en varios países, incluido el nuestro, aunque por cuestiones diversas, en España seguimos dependiendo de matrices de decisión y bases de datos pre-calculadas. Este sistema, desarrollado en la Universidad de Málaga, recibió el pasado mes de marzo el prestigioso galardón NVIDIA Global Impact Award, otorgado por primera vez en su historia a un organismo de investigación no-estadounidense y, además, ha superado los estándares requeridos por el NTHMP (35).

Pero para esta parte del estudio de los tsunamis, la propagación, solo se tiene en cuenta la deformación del fondo y su efecto en mar abierto. Al llegar a tierra, comienza la fase de cálculo de la inundación.

Aunque simular la inundación en tierra con modelos bidimensionales implica el uso de las mismas ecuaciones que la simulación en mar abierto, la demanda computacional es mucho mayor. Esto es debido a múltiples factores, entre los que se encuentran la divergencia y convergencia de flujos por la presencia de obstáculos en el terreno; a que se puede “secar” y “mojar” varias veces un mismo punto; pero sobre todo, a que la resolución necesaria para trabajar en la costa debe ser mucho mayor que en mar abierto. Por ejemplo, pasar de una resolución de 400x400 metros a 50x50 metros implica tardar unas 500 veces más; es decir, de 5 minutos a 2.560 minutos (casi dos días). 

Cuando se habla de simulación de inundaciones, en realidad lo que se genera no es “una simulación”. Lo normal es que se hagan múltiples simulaciones, de manera que aquellos errores detectables en una simulación sean solventados en la siguiente y así sucesivamente hasta que las incertidumbres no sean significativas para los objetivos de cálculo. No es infrecuente que el proceso pueda implicar decenas de ensayos previos. Por si fuera poco, generalmente se simulan distintos escenarios y se repiten todas las simulaciones “finales” para comprobar la estabilidad del resultado, aplicando variaciones a distintos parámetros (análisis de sensibilidad).

El modelo Tsunami-HySEA utiliza una estrategia de mallas estructuradas anidadas sobre el mismo esquema de cálculo, de manera que la de menor resolución (nodos de cálculo más alejados o píxeles más grandes) realiza una computación más ágil y sirve para aportar la condición de contorno a la siguiente malla de mayor resolución y así sucesivamente según convenga. El mismo código numérico realiza ambas fases de la simulación: la propagación en mar abierto y la inundación costera, adaptando los algoritmos numéricos implementados en cada caso de manera adecuada. Esto representa una gran ventaja a la hora de su utilización en la práctica, ya que evita pasos intermedios de procesado de datos para la adaptación del modelo de propagación al de inundación.

Para poder construir un modelo de simulación de tsunamis de origen sísmico hacen falta muchos datos. Por un lado, es necesario contar con un catálogo de eventos sísmicos y con un catálogo de fallas activas. Luego, hace falta conocer los parámetros que definen la fuente sísmica, la falla, que se simplifica a un plano horizontal (figura 7). Por otro lado, se requiere de información del relieve emergido susceptible de ser inundado (topografía) e información del relieve sumergido (batimetría), así como un mapa de las distintas rugosidades al flujo en la parte emergida.

Para el caso de estudio de España, el catálogo de eventos sísmicos y la topografía los proporciona el Instituto Geográfico Nacional de forma libre y para todo el público vía web. En cuanto a las fallas activas, el Instituto Geológico y Minero de España proporciona, en acceso abierto, la Base de Datos sobre Fallas Cuaternarias Activas (QAFI, Quaternary Active Faults of Iberia), así como otras fuentes de información necesarias para comprender el marco tectónico regional (el mapa GEODE –geología digital continua de España y productos derivados o asociados). En cuanto a las fuentes sísmicas tsunamigénicas, aún se está trabajando para consensuar los parámetros. En cuanto a la batimetría, el Instituto Geológico y Minero de España pertenece a la red EMODnet (European Marine Observation and Data Network (36) que genera coberturas batimétricas continuas lo más homogéneas y coherentes posible, utilizando datos de diversas campañas oceanográficas.  La última versión disponible es de septiembre de 2018, con un píxel de 250x250 metros. Además, es necesario contar con batimetrías de la zona próxima a la costa a mayor resolución, que pueden obtenerse de distintos proveedores, como el Instituto Hidrográfico de la Marina, el Ministerio para la Transición Ecológica  y proyectos de investigación oceanográfica (como el proyecto FAUCES (37), con socios de sobrado prestigio, como el Instituto Geológico y Minero de España, el Instituto de Ciencias del Mar [Centro Superior de Investigaciones Científicas], el Instituto Español de Oceanografía, la Universidad de Salamanca y la Comisión de Geología Marina de la Sociedad Geológica de España).

Uno de los cuellos de botella de toda simulación es la elección cuidadosa de las variables a obtener. A mayor cantidad de variables que se pida devolver al modelo, más se ralentiza el proceso de cálculo, pues se puede llegar a necesitar mucho tiempo en la escritura de grandes volúmenes de datos. Además, se requiere una gran inversión de tiempo en su post-procesado y presentación, de modo que resulte fácil de interpretar incluso para el público más experto. En todo caso, típicamente son objeto de interés las siguientes tres variables de la fase de propagación (figura 8 y figura 9): 
 

• las alturas de ola en su proximidad a la costa, pues aportan un valor de referencia clave para  estimar la “intensidad” esperable de un tsunami a lo largo de una costa; 
• el tiempo de llegada de la primera ola (positiva o negativa) o de la ola más significativa (representa la ventana máxima disponible para toda acción de autoprotección y alerta temprana); 
• registros mareográficos virtuales en puntos singulares (donde se puede observar la oscilación completa del tsunami a lo largo del tiempo). 

En cuanto a la parte de inundación, típicamente se usan las siguientes variables (figura 10):
 

• la altura máxima de remonte del tsunami (o run-up, que proporciona la información sobre el límite afectado tierra adentro en términos de altitud sobre el nivel del mar que alcanza el tsunami a lo largo de la costa);
• la profundidad máxima del agua en cada punto inundado (lo que incluye implícitamente la extensión máxima o alcance máximo);
• la velocidad máxima del agua en cada punto inundado (si bien es cierto que es menos común que la profundidad porque el tratamiento es muy exigente; para empezar, implica dos valores por píxel: módulo en dirección este-oeste y módulo en dirección norte-sur y para continuar, su variabilidad espacial es muy grande);
• el flujo específico del agua en cada punto (que integra el concepto clásico de volumen de agua movilizado por unidad de tiempo en cada punto, con o sin indicación de dirección). 

Todos estos parámetros se pueden obtener distribuidos en el tiempo, por ejemplo, al objeto de generar vídeos o material más detallado, pero dado que el volumen de los datos así expresado es muy grande, su procesado es muy lento. 

Si se parte del concepto de tsunami como el de una inundación, se podría dar por válido el esquema clásico que correlaciona velocidad de flujo y profundidad (o calado, figura 11) (40) con la probabilidad binaria de daño. Pero los tsunamis plantean muchos retos adicionales. Si se tratara de un único flujo, como en el caso de una inundación, el esquema sería perfectamente aplicable, pero como ya se apuntaba en apartados anteriores, un tsunami tiene dos flujos como mínimo, el de avance y el de resaca. Si en avance un tsunami causa daños, los daños que causa el tsunami en resaca se añaden a los anteriores. Incluso esta hipótesis sería válida para el intervalo del espacio de alcance máximo de la ola más significativa (que no tiene porqué ser la primera), pero en los intervalos donde alcanza la segunda ola más significativa, ya hay más de cuatro flujos. Por su parte, esta simplificación sería válida para tsunamis de origen lejano, pero para los tsunamis provocados por la causa más frecuente, los terremotos, y en caso de terremotos próximos a la costa, es posible esperar un cierto nivel de daños como consecuencia de la aceleración del suelo provocada por el propio terremoto. Tal y como se ha visto en el caso del reciente tsunami de Indonesia de 2018, los daños de un tsunami se pueden añadir a toda la cadena de efectos directos e indirectos de un terremoto: colapsos parciales o totales de infraestructuras (incapaces de soportar las cargas adicionales de la oscilación del suelo), licuefacción de suelos, movimientos de ladera (deslizamientos, desprendimientos), agrietamientos del suelo, y el sinfín de causas concatenadas (incendios, polución). 

La aproximación cualitativa para la evaluación del riesgo por tsunamis más utilizada (41) y que se pudo contrastar (42) tras el terremoto de Sumatra de 2004, es objeto de multitud de variantes para su aplicación en distintos casos de estudio. Se trata del modelo PTVA-3 (Papathoma Tsuami Vulnerability Assessment, o evaluación de la vulnerabilidad por tsunamis de Papathoma –el autor), que se expresa en forma del índice RVI (Relative Vulnerability Index, o índice de vulnerabilidad relativo). A pesar de que el autor original habla de vulnerabilidad, incorpora parámetros de peligrosidad (calado, orientación de la edificación respecto al flujo), por lo que en el mejor de los casos se trata de una vulnerabilidad específica, o incluso de riesgo relativo en el marco de los conceptos europeos de Zonas de Riesgos Naturales (43), (44), que se apoya en los conceptos de la UNISDR (45).

Cabe señalar que para el caso de estudio en España, habría que revisar la validez de este modelo, y en particular para su integración en la aproximación general desarrollada por el IGME para el Consorcio de Compensación de Seguros (46) (CCS) a efectos de cálculo cuantitativo de pérdidas económicas. La interpretación cuantitativa de otros trabajos similares podría ser particularmente útil en este sentido, como la Guía para la reducción de la vulnerabilidad de los edificios frente a las inundaciones, del CCS (47).

Del Área de Riesgos Geológicos del Instituto Geológico y Minero de España surgió a principios de los años 2000 una iniciativa para desarrollar estudios de riesgos geológicos que fue paralizada por falta de presupuesto, el Plan PRIGEO. Se buscaron apoyos y oportunidades y se alentaron iniciativas técnicas y científicas para el desarrollo de una herramienta que los técnicos y científicos estimamos “imprescindible” para un país moderno: cartografía de peligros naturales. Los principales peligros estaban considerados: las inundaciones, los terremotos, los movimientos de ladera y el volcanismo; y se hicieron varios esfuerzos por incluir otros peligros que por entonces eran menos conocidos fuera del sector técnico, como los tsunamis o las arcillas expansivas.

El primer peligro para el que Europa consiguió cuantificar y visualizar mediante una herramienta cartográfica fue el de inundación y ha sido solo gracias a la promulgación de la Directiva Europea de Inundaciones (29) que entró en vigor el año 2007 y a su ulterior implementación en las regulaciones nacionales. Menos de un año después, el IGME publicaba un manual técnico para la elaboración de mapas de inundaciones (48); y cuatro años después el MAPAMA publicaba la guía para elaborar el Sistema Nacional de Cartografía de Zonas Inundables (49)'(50).

Todo ello fue fruto de un importante esfuerzo previo de diferentes grupos de trabajo en los que el IGME y la Dirección General de Protección Civil y Emergencias (DGPCYE) pidieron, con insistencia, que se incluyeran los tsunamis como parte del análisis de los riesgos de inundación costeros. Pero los tsunamis quedaron explícitamente excluidos de la Directiva Europea de Inundaciones, de su transposición y relegados por completo. El tsunami de Sumatra de 2004 y el de Japón de 2011 dispararon el interés mundial por estos eventos y en la administración española creció el interés para abordar esta cuestión.

En 2015, se dio un importante paso en favor del estudio de los tsunamis en España: la DGPCYE publicaba la Directriz Básica correspondiente (51), y dos años después en 2017, se presentaba el estudio cualitativo de las áreas de interés (52), que destaca las costas Mediterráneas Peninsulares, el Golfo de Cádiz, Canarias y Baleares como áreas prioritarias frente a las costas gallegas y cántabras para abordar futuros trabajos.