Fractional-Order Viscoelastic Modeling of Seismic P-, S-, and Rayleigh-Wave Attenuation: Marmara Seismotectonic Context, Fault Geometry, and a Reproducible Strong-Motion Framework

Bajo la superficie de la Tierra, el suelo no se comporta como un simple bloque sólido de roca. Cuando un terremoto ocurre, la energía se desplaza hacia afuera en forma de ondas, pero a medida que estas ondas atraviesan la corteza, pierden energía y cambian de forma. Este proceso, conocido como atenuación, es la razón por la cual un suelo que tiembla se siente diferente a una distancia de la fuente que justo al lado de la ruptura. Durante décadas, los científicos han intentado describir esta pérdida de energía utilizando las reglas estándar de la física, tratando la corteza terrestre como un resorte perfectamente elástico o como un fluido viscoso simple. Sin embargo, el mundo real suele ser más complejo de lo que permiten estos dos extremos. El suelo es una mezcla desordenada de roca, fluido y grietas que recuerda sus movimientos pasados de una manera que la física estándar tiene dificultades para capturar. Aquí es donde una rama de las matemáticas llamada cálculo fraccionario resulta útil. En lugar de utilizar números enteros para describir cómo reaccionan los materiales, el cálculo fraccionario permite valores intermedios, ofreciendo una forma de modelar la "memoria" de la roca con muchas menos suposiciones.

En un nuevo estudio centrado en el Mar de Mármara, una región de intensa actividad geológica justo al sur de Estambul, los investigadores Taylan Demir y Atakan Koçyiğit han construido un marco matemático más preciso para describir cómo se desvanecen las ondas sísmicas. El Mar de Mármara se asienta sobre la Falla de Anatolia del Norte, una enorme grieta en la corteza terrestre donde las placas tectónicas se rozan entre sí. Debido a que esta falla pasa directamente por debajo de una gran ciudad, comprender exactamente cómo viajan y se debilitan las ondas sísmicas en esta zona específica es crítico para la ingeniería de edificios más seguros y para predecir la fuerza de futuros temblores. Los investigadores se propusieron probar si un modelo fraccionario, que tiene en cuenta el comportamiento complejo y dependiente de la historia de la roca, podía explicar mejor el desvanecimiento de las ondas sísmicas que los modelos tradicionales. No intentaron predecir cuándo ocurriría el próximo gran terremoto; en su lugar, se centraron en cómo reconstruir con precisión el temblor que ocurre después de que se registra un evento.

El equipo comenzó mapeando el escenario geológico. Examinaron los segmentos específicos de la falla que corren bajo el mar, incluyendo las ramas de Tekirdağ, la Cuenca Central y Kumburgaz, y analizaron la historia de los terremotos en la región, como el gran sismo de İzmit de 1999 y el evento de Silivri de 2019, que fue menor. Este contexto geológico era esencial porque las matemáticas que estaban desarrollando debían estar fundamentadas en la forma y estructura reales de la zona de la falla. Luego introdujeron una nueva forma de describir el comportamiento de la roca. Imaginen la corteza terrestre no como un simple resorte que vuelve a su posición instantáneamente, sino como un material que se relaja lentamente con el tiempo, conservando la memoria de cómo fue estirado o comprimido. Los investigadores utilizaron una herramienta matemática llamada derivada de Caputo para describir este efecto de memoria. Esta herramienta les permitió crear un modelo que conecta las propiedades físicas de la roca con la forma en que las ondas pierden energía mientras viajan.

Utilizando este nuevo modelo, los investigadores derivaron cómo deberían comportarse tres tipos diferentes de ondas sísmicas: las ondas P, que son las más rápidas y llegan primero; las ondas S, que sacuden el suelo de lado a lado; y las ondas Rayleigh, que ruedan a lo largo de la superficie. Demostraron matemáticamente que su modelo siempre resulta en una pérdida de energía, lo cual es una condición necesaria para cualquier descripción realista de un terremoto. Si un modelo sugiriera que las ondas pueden ganar energía mientras viajan, sería físicamente imposible. Sus ecuaciones mostraron que el enfoque fraccionario produce una curva de pérdida de energía suave y continua que coincide con la forma en que se comportan las rocas reales en un amplio rango de frecuencias. Esto representa una mejora significativa respecto a los modelos antiguos, que a menudo tienen dificultades para ajustarse a los datos en diferentes frecuencias sin volverse excesivamente complicados.

Para probar si esta teoría realmente funciona, el equipo realizó un experimento controlado utilizando datos generados por computadora. Crearon un escenario de terremoto sintético donde conocían exactamente las reglas de cómo debían desvanecerse las ondas. Luego, pidieron a su modelo fraccionario que descifrara esas reglas basándose únicamente en los patrones de las ondas. Los resultados fueron sorprendentes. El modelo fraccionario recuperó los parámetros correctos con una precisión muy alta, reduciendo el error en sus predicciones en casi un 87 por ciento en comparación con un modelo tradicional que utilizaba reglas de números enteros. En el modelo tradicional, los investigadores tuvieron que forzar los números para que encajaran, lo que resultó en un ajuste deficiente con el comportamiento real de la onda. El modelo fraccionario, al permitir ese efecto de memoria "intermedio", capturó la física de la pérdida de la onda de una manera mucho más natural.

El estudio también analizó cómo este marco podría aplicarse a datos del mundo real de la región de Mármara. Los investigadores identificaron el terremoto de Silivri de 2019 como un candidato perfecto para una primera prueba real. Esbozaron un plan claro sobre cómo los ingenieros y sismólogos podrían descargar grabaciones reales de este evento y utilizar sus nuevas ecuaciones para reconstruir la intensidad del temblor. Esto permitiría calcular medidas precisas del movimiento del suelo, como la aceleración o velocidad máxima, que son los números que los arquitectos utilizan para diseñar edificios resistentes a los terremotos. Los autores fueron cuidadosos al declarar que esta herramienta es para analizar lo que ya ha sucedido, no para predecir el futuro. No nos dice cuándo ocurrirá la próxima gran ruptura, pero proporciona una lente mucho más nítida para comprender la violencia del temblor una vez que ocurre.

Al combinar una prueba matemática rigurosa con un mapa geológico claro y un plan para el análisis de datos reales, este artículo traslada la conversación de la teoría abstracta a la aplicación práctica. Demuestra que el enfoque fraccionario no es solo una curiosidad matemática, sino una forma viable y más precisa de describir la compleja realidad de las ondas sísmicas en el Mar de Mármara. El trabajo cierra la brecha entre la geología desordenada y estratificada de la zona de la falla y las ecuaciones limpias y predecibles necesarias para la seguridad en la ingeniería. Si bien la calibración completa frente a miles de registros de terremotos reales es una tarea para el futuro, este estudio ha sentado las bases necesarias, demostrando que el método funciona en la teoría y en la simulación, y está listo para ser probado contra el suelo real que se encuentra bajo Estambul.