Reconstruction of the non-linear wave at a buoy from shoreline data and applications to the tsunami inverse problem for piece-wise sloping bathymetry

Basic Information

• Paper ID: 2510.14177
• Title: Reconstruction of the non-linear wave at a buoy from shoreline data and applications to the tsunami inverse problem for piece-wise sloping bathymetry
• Authors: Oleksandr Bobrovnikov¹, Madison Jones², Shriya Prasanna³, Josiah Smith², Alexei Rybkin¹, Efim Pelinovsky⁴'⁵
• Affiliations: ¹University of Alaska Fairbanks, ²University of Colorado Boulder, ³University of Washington, ⁴HSE University, ⁵Institute of Applied Physics
• Classification: math.AP (Applied Mathematics Analysis)
• Publication Date: October 17, 2025
• Paper Link: https://arxiv.org/abs/2510.14177

Abstract

This paper addresses an important inverse problem: given tsunami wave runup data at the shoreline, can we recover its initial shape? The research is conducted within the framework of nonlinear shallow water equations, a model widely used in tsunami propagation and inundation studies. The authors consider finite-slope bathymetry and demonstrate that boundary conditions (water surface elevation and velocity) at a virtual buoy location can be recovered from shoreline data. Furthermore, they discuss coupling shallow water equations with Boussinesq equations in more complex piece-wise sloping bathymetry to recover initial conditions and introduce dispersive effects into the model.

Research Background and Motivation

Importance of the Problem

1. Disaster Threat: Major tsunami disasters such as the 1964 Alaska earthquake, 2004 Indian Ocean tsunami, and 2011 Great East Japan Earthquake have caused enormous casualties and economic losses, highlighting the importance of tsunami early warning and disaster mitigation.
2. Early Warning System Requirements: With the development of DART (Deep-ocean Assessment and Reporting of Tsunamis) systems, theoretical frameworks are needed to guide buoy placement and improve the accuracy and timeliness of tsunami warnings.
3. Mathematical Modeling Challenges: Tsunami propagation involves complex nonlinear effects, particularly in nearshore regions, requiring accurate mathematical models to describe wave behavior.

Limitations of Existing Methods

Existing tsunami inversion methods primarily include four approaches:

1. Methods combining seismic and waveform data
2. Reconstruction methods using waveform data alone
3. Inversion methods utilizing sediment distribution
4. Inversion methods based on inundation data

However, these methods have limitations when dealing with finite-slope bathymetry, particularly in addressing the combined effects of nonlinear and dispersive effects.

Research Motivation

The research motivation is to extend the authors' previous work on infinite-slope bathymetry to develop inversion theory applicable to more realistic finite-slope bathymetry and provide a theoretical framework for guiding buoy placement.

Core Contributions

1. Theoretical Extension: Extends tsunami inversion theory from infinite-slope to finite-slope bathymetry, providing complete analytical solutions.
2. Algorithm Development: Proposes a complete algorithmic procedure for recovering boundary conditions at a virtual buoy from shoreline runup data.
3. Model Coupling: Develops a method for coupling nonlinear shallow water equations (NSWE) with Boussinesq equations in piece-wise sloping bathymetry, simultaneously accounting for nonlinear and dispersive effects.
4. Numerical Verification: Validates the method's effectiveness through multiple physically realistic scenarios, including different types of initial waveforms.
5. Computational Efficiency: Provides an efficient algorithm implementation with O(N³) complexity, supporting rapid tsunami characteristic calculations.

Methodology Details

Problem Definition

Input: Shoreline runup data R(t) Output: Water surface elevation η(L,t) and velocity u(L,t) at a virtual buoy Constraints: Non-breaking wave condition, validity of Carrier-Greenspan transformation

Core Mathematical Framework

1. Nonlinear Shallow Water Equation System

In dimensionless units, system (1) is:

∂tη + u∂x(x + η) + (x + η)∂xu = 0  (mass conservation)
∂tu + u∂xu + ∂xη = 0              (momentum conservation)

2. Carrier-Greenspan Transformation

The key linearizing transformation (5):

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

This transformation converts the nonlinear system (1) into a linear system (6):

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

3. Forward Problem Solution

By introducing new variables ρ² = σ/σL and function θ(ρ,τ) = ψ(σ,τ) - ψb(τ), employing Fourier decomposition:

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

where J₀ is the Bessel function of the first kind, and jn is the n-th zero of J₀.

4. Core Inversion Algorithm Formula

Shoreline equation (20):

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

Inversion formula (24):

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

Complete Inversion Algorithm Procedure

1. Data Transformation: Convert R(t) to ψsh(τ) using CGT
2. Boundary Recovery: Calculate ψb(τ) from ψsh(τ) using formula (24)
3. Field Reconstruction: Compute ψ(σ,τ) and φ(σ,τ) throughout the domain using formulas (15,16)
4. Inverse Transformation: Perform inverse CGT to recover physical boundary conditions u(L,t) and η(L,t)

Technical Innovations

1. Analytical Solution Method: Unlike numerical methods, provides completely analytical inversion solutions with high computational efficiency and controllable precision.
2. Laplace Transform Technique: Cleverly employs Laplace transform and inverse fast Fourier transform (iFFT) to enhance computational efficiency.
3. Model Coupling Strategy: Innovatively couples NSWE with Boussinesq equations to simultaneously handle nonlinear and dispersive effects.

Experimental Setup

Numerical Verification Scenarios

The authors designed four different initial displacement scenarios for verification:

1. Bimodal Waveform: η(x,0) = 0.005 sech(x-16) - 0.003e^(-(x-13)²)
2. Soliton Waveform: η(x,0) = 0.005 cosh⁻²(2x-6)
3. Gaussian Waveform: η(x,0) = 0.005e^(-(x-7)²)
4. Multi-peak Waveform: η(x,0) = 0.005e^(-2(x-6)²) + 0.003e^(-(x-10)²)

Evaluation Metrics

• Accuracy Comparison: Comparing recovered ψb(τ) with exact value ψᵉb(τ)
• Boundary Condition Recovery: Verifying recovery accuracy of u(L,t) and η(L,t)
• Computational Efficiency: Measuring algorithm runtime and complexity

Implementation Details

• Boundary Setting: L = 1 (dimensionless units)
• Data Points: N = 1500 data points
• Truncation Parameter: Bessel function zeros truncated to n = 500
• Computing Platform: Google Colab environment

Experimental Results

Main Results

1. Inversion Accuracy Verification

From Figures 4 and 5, for all four test scenarios, the recovered boundary conditions ψb(τ) and φb(τ) show excellent agreement with exact solutions, validating the algorithm's effectiveness.

2. Computational Performance

• Runtime: Computation of 1500 data points takes 76 seconds on Google Colab
• Complexity Verification: Time complexity conforms to theoretical O(N³) analysis
• Fitting Result: Least squares fitting slope of time versus data points is 2.82, close to theoretical value of 3

3. Physical Boundary Condition Recovery

Figure 7 displays recovered physical boundary conditions u(L,t) and η(L,t) for four scenarios, demonstrating the algorithm's capability to accurately reconstruct wave field information at the virtual buoy.

Boussinesq Equation Coupling Experiments

1. Bi-soliton Solution Verification

Using Hirota bi-soliton solution as boundary conditions, successfully verifies the coupling method between NSWE and Boussinesq equations. Experimental parameters:

• q₁ = 0.1, q₂ = √0.1
• Various combinations of time parameters t₁, t₂

2. Initial Condition Reconstruction

Figure 13 shows initial condition reconstruction results at tsunami event time (tevent = -100s):

• Exact peak positions: x₁ = 119.55, x₂ = 123.63
• Recovered peak positions: x₁ = 119.50, x₂ = 123.64
• Minimal error, validating method accuracy

Linear Shallow Water Equation Comparison

Figure 16 displays comparison of initial states between Boussinesq model and linear shallow water equation model under identical boundary conditions, highlighting the impact of dispersive effects.

Related Work

Tsunami Inversion Research Directions

1. Seismic-Waveform Joint Inversion: Yokota et al. (2011), Yue et al. (2014)
2. Pure Waveform Data Inversion: Ho et al. (2019), Yamanaka and Tanioka (2024)
3. Sediment Inversion: Tang and Weiss (2015), Ioki and Tanioka (2016)
4. Inundation Data Inversion: Lee et al. (2021), Piatanesi et al. (1996)

Numerical Method Development

• Finite Difference/Finite Element: Kounadis and Dougalis (2020), Xing (2017)
• Physics-Informed Neural Networks: Brecht et al. (2025)
• Piece-wise Bathymetry Modeling: Synolakis (1987), Varsoliwala and Singh (2021)

Relative Advantages of This Work

Compared to existing work, this paper provides completely analytical inversion methods, avoiding cumulative errors from numerical methods, and systematically addresses the inversion problem for finite-slope bathymetry for the first time.

Conclusions and Discussion

Main Conclusions

1. Theoretical Completeness: Successfully extends tsunami inversion theory to finite-slope bathymetry, providing complete analytical solutions.
2. Practical Value: The algorithm can guide optimal placement of tsunami warning buoys such as DART systems, improving warning efficiency.
3. Model Integration: Successfully achieves coupling of NSWE with Boussinesq equations, providing new pathways for tsunami modeling in complex bathymetry.

Limitations

1. Non-breaking Wave Assumption: The algorithm applies only to non-breaking wave cases; the method fails when the Jacobian determinant becomes zero.
2. One-dimensional Model Limitation: Currently considers only one-dimensional cases, while actual tsunami propagation is a three-dimensional phenomenon.
3. Bathymetry Simplification: Although piece-wise bathymetry is considered, it remains highly simplified compared to actual seafloor topography.
4. Dispersive Effect Range: The applicability of Boussinesq equations remains controversial in certain cases.

Future Directions

1. Multi-dimensional Extension: Develop two-dimensional and three-dimensional inversion theory
2. Complex Bathymetry: Address inversion problems for multi-segment sloping bathymetry
3. Breaking Wave Criteria: Establish clear acceptability criteria for shoreline data
4. Practical Applications: Combine with actual tsunami observation data for validation

In-depth Evaluation

Strengths

1. Mathematical Rigor: Based on rigorous mathematical theory, providing complete analytical solutions that avoid error accumulation from numerical methods.
2. Computational Efficiency: The O(N³) complexity is acceptable for tsunami warning real-time requirements, with particularly evident computational speed advantages for one-dimensional models.
3. Strong Innovation: Systematically solves the tsunami inversion problem for finite-slope bathymetry for the first time and innovatively achieves coupling of different equations.
4. Sufficient Verification: Demonstrates method stability and accuracy through numerical verification across multiple waveform scenarios.

Shortcomings

1. Limited Applicability: The non-breaking wave assumption restricts method application in strongly nonlinear cases.
2. Simplified Bathymetry Model: Actual seafloor bathymetry is far more complex than piece-wise slopes; practical application effectiveness requires verification.
3. Lack of Real Data Validation: Uses only synthetic data for verification, lacking comparison with actual tsunami observation data.
4. Boundary Condition Compatibility: While compatibility conditions are mentioned, their practical impact requires more detailed discussion.

Impact

1. Theoretical Contribution: Provides important mathematical foundations for tsunami inversion theory, particularly in analytical methods.
2. Practical Value: Offers theoretical guidance for buoy placement in tsunami warning systems, with significant disaster mitigation importance.
3. Methodological Significance: The innovative application of Carrier-Greenspan transformation provides insights for other nonlinear wave problems.

Applicable Scenarios

1. Tsunami Warning System Design: Guides optimal placement of monitoring equipment such as DART buoys
2. Historical Tsunami Reconstruction: Reconstructs historical events based on paleotsunamis sediment data
3. Theoretical Research: Provides mathematical tools for tsunami propagation mechanism research
4. Engineering Applications: Reference for coastal protection engineering design

References

The paper cites abundant relevant literature, primarily including:

• Carrier and Greenspan (1958): Original work on CGT transformation
• Rybkin et al. (2021, 2023, 2024, 2025): Series of works by the author team
• Synolakis (1987): Classical work on piece-wise bathymetry modeling
• Hirota (1973): Important contribution to soliton solutions of Boussinesq equations

Overall Assessment: This is a high-quality mathematical physics paper that makes important contributions to tsunami inversion theory. The method has rigorous mathematical foundations and good computational efficiency, with significant theoretical and practical value for tsunami warning and disaster mitigation. Although certain limitations exist, it provides a solid foundation for further development in this field.