Reconstrucción de la onda no lineal en una boya a partir de datos de costa y aplicaciones al problema inverso de tsunami para batimetría de pendiente por tramos

Información Básica

• ID del Artículo: 2510.14177
• Título: Reconstruction of the non-linear wave at a buoy from shoreline data and applications to the tsunami inverse problem for piece-wise sloping bathymetry
• Autores: Oleksandr Bobrovnikov¹, Madison Jones², Shriya Prasanna³, Josiah Smith², Alexei Rybkin¹, Efim Pelinovsky⁴'⁵
• Instituciones de los Autores: ¹Universidad de Alaska Fairbanks, ²Universidad de Colorado Boulder, ³Universidad de Washington, ⁴Universidad HSE, ⁵Instituto de Física Aplicada
• Clasificación: math.AP (Análisis Matemático Aplicado)
• Fecha de Publicación: 17 de octubre de 2025
• Enlace del Artículo: https://arxiv.org/abs/2510.14177

Resumen

Este artículo aborda un importante problema inverso: ¿es posible recuperar la forma inicial de una onda de tsunami dada la información de su ascenso en la costa? La investigación se realiza dentro del marco de las ecuaciones no lineales de aguas someras, un modelo ampliamente utilizado en estudios de propagación y inundación de tsunamis. Los autores consideran una batimetría de fondo marino con pendiente finita y demuestran que es posible recuperar las condiciones de frontera en una boya virtual (desplazamiento del nivel del agua y velocidad) a partir de datos de la línea de costa. Además, se discute el acoplamiento de las ecuaciones de aguas someras con las ecuaciones de Boussinesq en batimetrías de pendiente por tramos más complejas, para recuperar las condiciones iniciales e introducir efectos dispersivos en el modelo.

Antecedentes de Investigación y Motivación

Importancia del Problema

1. Amenaza de Desastres: El terremoto de Alaska de 1964, el tsunami del Océano Índico de 2004 y el Gran Terremoto de Japón Oriental de 2011 causaron enormes pérdidas de vidas humanas y económicas, destacando la importancia de la alerta temprana de tsunamis y la mitigación de desastres.
2. Necesidad de Sistemas de Alerta: Con el desarrollo del sistema DART (Evaluación Profunda del Océano y Reporte de Tsunamis), existe la necesidad de un marco teórico que guíe el despliegue de boyas para mejorar la precisión y oportunidad de las alertas de tsunami.
3. Desafíos de Modelado Matemático: La propagación de tsunamis implica efectos no lineales complejos, particularmente en regiones cercanas a la costa, requiriendo modelos matemáticos precisos para describir el comportamiento de las ondas.

Limitaciones de Métodos Existentes

Los métodos de inversión de tsunamis existentes incluyen principalmente cuatro enfoques:

1. Métodos que combinan datos sísmicos y de forma de onda
2. Métodos de reconstrucción utilizando solo datos de forma de onda
3. Métodos de inversión utilizando distribución de sedimentos
4. Métodos de inversión basados en datos de inundación

Sin embargo, estos métodos tienen limitaciones al tratar batimetrías de pendiente finita, particularmente en considerar los efectos combinados de no linealidad y dispersión.

Motivación de la Investigación

La motivación de este trabajo es extender investigaciones previas de los autores en batimetrías de pendiente infinita, desarrollar una teoría de inversión aplicable a batimetrías de pendiente finita más realistas, y proporcionar un marco teórico para guiar el despliegue de boyas.

Contribuciones Principales

1. Extensión Teórica: Ampliar la teoría de inversión de tsunamis desde fondos marinos de pendiente infinita a batimetrías de pendiente finita, proporcionando una solución analítica completa.
2. Desarrollo de Algoritmos: Proponer un procedimiento algorítmico completo para recuperar las condiciones de frontera en una boya virtual a partir de datos de ascenso en la costa.
3. Acoplamiento de Modelos: Desarrollar un método para acoplar las ecuaciones no lineales de aguas someras (NSWE) con las ecuaciones de Boussinesq en batimetrías de pendiente por tramos, considerando simultáneamente efectos no lineales y dispersivos.
4. Verificación Numérica: Validar la efectividad del método a través de múltiples escenarios físicamente realistas, incluyendo diferentes tipos de formas de onda iniciales.
5. Eficiencia Computacional: Proporcionar una implementación de algoritmo eficiente con complejidad O(N³), que permite cálculos rápidos de características de tsunamis.

Explicación Detallada de Métodos

Definición de la Tarea

Entrada: Datos de ascenso en la costa R(t) Salida: Desplazamiento del nivel del agua η(L,t) y velocidad u(L,t) en la boya virtual Restricciones: Condición de onda no rompiente, validez de la transformación de Carrier-Greenspan

Marco Matemático Principal

1. Sistema de Ecuaciones No Lineales de Aguas Someras

En unidades adimensionales, el sistema (1) es:

∂tη + u∂x(x + η) + (x + η)∂xu = 0  (Conservación de masa)
∂tu + u∂xu + ∂xη = 0              (Conservación de momento)

2. Transformación de Carrier-Greenspan

La transformación de linealización clave (5):

φ(σ,τ) = u(x,t)
σ = x + η(x,t)
ψ(σ,τ) = η(x,t) + u²(x,t)/2
τ = t - u(x,t)

Esta transformación convierte el sistema no lineal (1) en el sistema lineal (6):

∂τψ + σ∂σφ + φ = 0
∂τφ + ∂σψ = 0

3. Solución del Problema Directo

Mediante la introducción de nuevas variables ρ² = σ/σL y la función θ(ρ,τ) = ψ(σ,τ) - ψb(τ), se utiliza descomposición de Fourier:

θ(τ,ρ) = Σ cn(τ)J₀(jnρ)

donde J₀ es la función de Bessel de primera clase y jn es el n-ésimo cero de J₀.

4. Fórmulas Principales del Algoritmo de Inversión

Ecuación de costa (20):

ψsh(τ) = Σ (bn/√an) ∫₀τ sin(√an(τ-ξ))ψ''b(ξ)dξ + ψb(τ)

Fórmula de inversión (24):

ψb(τ) = L⁻¹[Lψsh/(s²Σ bn/(an+s²) + 1)⁻¹]

Procedimiento Completo del Algoritmo de Inversión

1. Conversión de Datos: Utilizar CGT para convertir R(t) en ψsh(τ)
2. Recuperación de Frontera: Calcular ψb(τ) a partir de ψsh(τ) mediante la fórmula (24)
3. Reconstrucción de Campo: Utilizar las fórmulas (15,16) para calcular ψ(σ,τ) y φ(σ,τ) en todo el dominio
4. Transformación Inversa: Ejecutar la CGT inversa para recuperar las condiciones de frontera físicas u(L,t) y η(L,t)

Puntos de Innovación Técnica

1. Método de Solución Analítica: A diferencia de los métodos numéricos, proporciona una solución de inversión completamente analítica con alta eficiencia computacional y precisión controlable.
2. Técnica de Transformada de Laplace: Aplicación ingeniosa de la transformada de Laplace y la transformada rápida de Fourier inversa (iFFT) para mejorar la eficiencia computacional.
3. Estrategia de Acoplamiento de Modelos: Acoplamiento innovador de NSWE con ecuaciones de Boussinesq, manejando simultáneamente efectos no lineales y dispersivos.

Configuración Experimental

Escenarios de Verificación Numérica

Los autores diseñaron cuatro escenarios de desplazamiento inicial diferente para verificación:

1. Forma de Onda Bimodal: η(x,0) = 0.005 sech(x-16) - 0.003e^(-(x-13)²)
2. Forma de Onda de Solitón: η(x,0) = 0.005 cosh⁻²(2x-6)
3. Forma de Onda Gaussiana: η(x,0) = 0.005e^(-(x-7)²)
4. Forma de Onda Multimodal: η(x,0) = 0.005e^(-2(x-6)²) + 0.003e^(-(x-10)²)

Indicadores de Evaluación

• Comparación de Precisión: Comparar el ψb(τ) recuperado con el valor exacto ψᵉb(τ)
• Recuperación de Condiciones de Frontera: Verificar la precisión de recuperación de u(L,t) y η(L,t)
• Eficiencia Computacional: Medir el tiempo de ejecución del algoritmo y la complejidad

Detalles de Implementación

• Configuración de Frontera: L = 1 (unidades adimensionales)
• Número de Puntos de Datos: N = 1500 puntos de datos
• Parámetro de Truncamiento: Ceros de función de Bessel truncados en n = 500
• Plataforma de Cálculo: Entorno Google Colab

Resultados Experimentales

Resultados Principales

1. Verificación de Precisión de Inversión

A partir de las figuras 4 y 5, se puede observar que para todos los cuatro escenarios de prueba, la consistencia entre las condiciones de frontera recuperadas ψb(τ) y φb(τ) y la solución exacta es excelente, validando la efectividad del algoritmo.

2. Rendimiento Computacional

• Tiempo de Ejecución: El cálculo de 1500 puntos de datos requirió 76 segundos en Google Colab
• Verificación de Complejidad: La complejidad de tiempo se ajusta al análisis teórico de O(N³)
• Resultado de Ajuste: La pendiente del ajuste de mínimos cuadrados entre tiempo y número de puntos de datos es 2.82, cercana al valor teórico de 3

3. Recuperación de Condiciones de Frontera Física

La figura 7 muestra las condiciones de frontera física recuperadas u(L,t) y η(L,t) en los cuatro escenarios, demostrando que el algoritmo puede reconstruir con precisión la información del campo de ondas en la boya virtual.

Experimentos de Acoplamiento de Ecuaciones de Boussinesq

1. Verificación de Solución de Doble Solitón de Hirota

Utilizando la solución de doble solitón de Hirota como condición de frontera, se verificó exitosamente el método de acoplamiento entre NSWE y ecuaciones de Boussinesq. Parámetros experimentales:

• q₁ = 0.1, q₂ = √0.1
• Diferentes combinaciones de parámetros de tiempo t₁, t₂

2. Reconstrucción de Condiciones Iniciales

La figura 13 muestra los resultados de reconstrucción de condiciones iniciales en el momento del evento de tsunami (tevent = -100s):

• Posición de pico exacta: x₁ = 119.55, x₂ = 123.63
• Posición de pico recuperada: x₁ = 119.50, x₂ = 123.64
• Error mínimo, validando la precisión del método

Comparación con Ecuaciones Lineales de Aguas Someras

La figura 16 muestra la comparación entre el modelo de Boussinesq y el modelo lineal de aguas someras bajo las mismas condiciones de frontera, destacando el impacto de los efectos dispersivos.

Trabajos Relacionados

Direcciones de Investigación en Inversión de Tsunamis

1. Inversión Conjunta Sísmica-Forma de Onda: Yokota et al. (2011), Yue et al. (2014)
2. Inversión de Datos Puros de Forma de Onda: Ho et al. (2019), Yamanaka and Tanioka (2024)
3. Inversión de Sedimentos: Tang and Weiss (2015), Ioki and Tanioka (2016)
4. Inversión de Datos de Inundación: Lee et al. (2021), Piatanesi et al. (1996)

Desarrollo de Métodos Numéricos

• Diferencias Finitas/Elementos Finitos: Kounadis and Dougalis (2020), Xing (2017)
• Redes Neuronales Informadas por Física: Brecht et al. (2025)
• Modelado de Topografía por Tramos: Synolakis (1987), Varsoliwala and Singh (2021)

Ventajas Relativas de Este Trabajo

En comparación con trabajos existentes, este artículo proporciona un método de inversión completamente analítico, evitando errores acumulativos de métodos numéricos, y es el primero en abordar sistemáticamente el problema de inversión para topografías de pendiente finita.

Conclusiones y Discusión

Conclusiones Principales

1. Completitud Teórica: Extensión exitosa de la teoría de inversión de tsunamis a batimetrías de fondo marino con pendiente finita, proporcionando una solución analítica completa.
2. Valor Práctico: El algoritmo puede utilizarse para guiar la colocación óptima de boyas de sistemas de alerta de tsunami como DART, mejorando la eficiencia de alerta.
3. Integración de Modelos: Acoplamiento exitoso de NSWE con ecuaciones de Boussinesq, proporcionando nuevas vías para el modelado de tsunamis en topografías más complejas.

Limitaciones

1. Suposición de Onda No Rompiente: El algoritmo solo es aplicable en casos de ondas no rompientes; el método falla cuando el determinante jacobiano es cero.
2. Limitación de Modelo Unidimensional: Actualmente solo se considera el caso unidimensional, mientras que la propagación real de tsunamis es un fenómeno tridimensional.
3. Simplificación de Topografía: Aunque se considera topografía por tramos, sigue siendo mucho más simplificada que la topografía real del fondo marino.
4. Rango de Efectos Dispersivos: La aplicabilidad de las ecuaciones de Boussinesq sigue siendo controvertida en algunos casos.

Direcciones Futuras

1. Extensión Multidimensional: Desarrollar teoría de inversión bidimensional y tridimensional
2. Topografía Compleja: Abordar problemas de inversión con batimetrías de pendiente múltiple por tramos
3. Criterios de Onda Rompiente: Establecer criterios claros de aceptabilidad para datos de costa
4. Aplicación Práctica: Validación con datos reales de observación de tsunamis

Evaluación Profunda

Fortalezas

1. Rigor Matemático: Basado en teoría matemática rigurosa, proporciona soluciones analíticas completas, evitando la acumulación de errores de métodos numéricos.
2. Eficiencia Computacional: La complejidad O(N³) es aceptable bajo los requisitos de tiempo real de alerta de tsunami, con ventajas evidentes de velocidad computacional para modelos unidimensionales.
3. Fuerte Innovación: Primera solución sistemática del problema de inversión de tsunami para topografías de pendiente finita, con acoplamiento innovador de diferentes ecuaciones.
4. Verificación Suficiente: Verificación numérica a través de múltiples escenarios de forma de onda, demostrando la estabilidad y precisión del método.

Deficiencias

1. Limitación de Rango de Aplicabilidad: La suposición de onda no rompiente limita la aplicación del método en casos altamente no lineales.
2. Simplificación del Modelo de Topografía: La topografía real del fondo marino es mucho más compleja que la topografía por tramos, y la efectividad de la aplicación práctica del método requiere verificación.
3. Falta de Verificación con Datos Reales: Solo se utiliza datos sintéticos para verificación, faltando comparación con datos reales de observación de tsunamis.
4. Compatibilidad de Condiciones de Frontera: Aunque se mencionan condiciones de compatibilidad, no se discuten en detalle sus impactos en aplicaciones prácticas.

Impacto

1. Contribución Teórica: Proporciona una base matemática importante para la teoría de inversión de tsunamis, particularmente en métodos analíticos.
2. Valor Práctico: Proporciona orientación teórica para el despliegue de boyas en sistemas de alerta de tsunami, con importante significado en mitigación de desastres.
3. Significado Metodológico: La aplicación innovadora de la transformación de Carrier-Greenspan proporciona referencias para otros problemas de ondas no lineales.

Escenarios de Aplicación

1. Diseño de Sistemas de Alerta de Tsunami: Guiar la colocación óptima de dispositivos de monitoreo como boyas DART
2. Reconstrucción de Tsunamis Históricos: Reconstruir eventos históricos basándose en datos de sedimentos de tsunamis antiguos
3. Investigación Teórica: Proporcionar herramientas matemáticas para investigación de mecanismos de propagación de tsunamis
4. Aplicación de Ingeniería: Referencia de diseño para ingeniería de protección costera

Referencias

El artículo cita una amplia literatura relacionada, incluyendo principalmente:

• Carrier and Greenspan (1958): Trabajo original de la transformación CGT
• Rybkin et al. (2021, 2023, 2024, 2025): Serie de trabajos del equipo de autores
• Synolakis (1987): Trabajo clásico en modelado de topografía por tramos
• Hirota (1973): Contribución importante a soluciones de solitones de ecuaciones de Boussinesq

Evaluación General: Este es un artículo de alta calidad en física matemática que realiza contribuciones importantes a la teoría de inversión de tsunamis. El método tiene una base matemática rigurosa y buena eficiencia computacional, con importante valor teórico y práctico para alerta de tsunami y mitigación de desastres. Aunque tiene algunas limitaciones, proporciona una base sólida para desarrollo futuro en este campo.