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

Titel

Reconstructie van de niet-lineaire golf op een boei vanuit kustdata en toepassingen op het inverse probleem voor de tsunami bij een kust met een stuksgewijs hellende bodem.

Probleemstelling

Het artikel adresseert het inverse probleem voor tsunami's: gegeven de oploopdata (run-up) van een tsunami-golf aan de kust, kan men de initiële vorm van de golf (de initiële verplaatsing en snelheid) of de randvoorwaarden op een virtuele boei terugvinden?

Traditionele modellen gebruiken vaak lineaire benaderingen of beschouwen oneindig hellende bodems. Deze studie richt zich op een meer realistische scenario:

1. Een eindig hellend strand (finite sloping bathymetry) waarbij niet-lineaire effecten dominant zijn.
2. Een stuksgewijs hellende bodem (piece-wise sloping), bestaande uit een hellend gebied (Region 1) en een vlakke, semi-oneindige regio (Region 2) waar dispersie belangrijk wordt.
3. Het doel is om de golfveld-data op een virtuele boei (buiten de kust) te reconstrueren op basis van waarnemingen aan de kustlijn, en vervolgens de initiële condities van de tsunami te bepalen.

Methodologie

De auteurs hanteren een analytische en numerieke aanpak die is opgebouwd uit drie hoofdfasen:

1. Het Inverse Probleem voor de Niet-Lineaire Ondiepe Watervergelijkingen (NSWE)
Voor het hellende gebied (Region 1) worden de NSWE (Non-linear Shallow Water Equations) gebruikt. Omdat deze vergelijkingen niet-lineair zijn en een vrij oppervlak hebben, is een directe oplossing moeilijk.

• Carrier-Greenspan Transformatie (CGT): De auteurs gebruiken deze exacte transformatie om de niet-lineaire NSWE te lineariseren voor een specifiek type bodemprofiel (lineair hellend). Dit transformeert het probleem naar het hodograafvlak $(\sigma, \tau)$.
• Linearisatie en Oplossing: In het hodograafvlak wordt het probleem een lineaire partiële differentiaalvergelijking (een vergelijking van Bessel-type). De auteurs gebruiken een Fourier-decompositie (met Bessel-functies) om de oplossing voor de directe problemen (van randvoorwaarden naar kustoploop) af te leiden.
• Inverse Stap: Door de Laplace-transformatie toe te passen op de afgeleide "kustvergelijking" (shoreline equation), kunnen ze de relatie omkeren. Dit stelt hen in staat om de randvoorwaarden op de boei ($\psi_b, \phi_b$) te berekenen uit de gemeten kustoploop $R(t)$.
• Terugtransformatie: Via de inverse CGT worden de fysieke variabelen (waterverplaatsing $\eta$ en snelheid $u$) op de boei teruggevonden.

2. Numerieke Validatie
De auteurs hebben een "manufactured solution" (een kunstmatige oplossing) gegenereerd om het algoritme te testen:

• Ze simuleerden een voorwaartse probleem met verschillende initiële profielen (bijv. soliton-achtige vormen).
• Vervolgens gebruikten ze het inverse algoritme op de gegenereerde kustdata.
• Resultaat: De gereconstrueerde randvoorwaarden kwamen zeer nauwkeurig overeen met de oorspronkelijke data, wat de validiteit van de analytische methode bevestigt.

3. Koppeling met Dispersieve Modellen (Boussinesq en Lineaire SWE)
Om dispersie (golfspreiding) te modelleren in de diepere wateren (Region 2), koppelen ze het NSWE-model aan andere vergelijkingen:

• Boussinesq-vergelijking: Gebruikt voor het niet-lineaire, dispersieve regime. De auteurs gebruiken een heuristische aanpak waarbij ze aannemen dat de golf op de grens van de twee regio's kan worden beschreven als een som van solitons (meerdere solitons). Ze passen een "least squares"-fit toe op de gereconstrueerde data om de parameters van deze solitons (snelheid, amplitude) te bepalen en zo de initiële condities in het diepe water te reconstrueren.
• Lineaire Ondiepe Watervergelijkingen (SWE): Als alternatief wordt een model getest waarbij beide regio's door lineaire SWE worden bestuurd, vergelijkbaar met eerdere werken van Synolakis (1987).

Belangrijkste Bijdragen

1. Analytische Oplossing voor het Inverse Probleem: Voor het eerst wordt een volledig analytisch algoritme gepresenteerd om de initiële condities van een tsunami op een virtuele boei te reconstrueren, puur op basis van kustoploopdata, binnen het kader van de niet-lineaire ondiepe watervergelijkingen op een eindig hellend strand.
2. Overbrugging van Modellen: De studie introduceert een methode om niet-lineaire modellen (NSWE) in ondiep water te "stitchen" (koppelen) met dispersieve modellen (Boussinesq) in dieper water. Dit maakt een meer realistische reconstructie mogelijk die rekening houdt met zowel niet-lineariteit als dispersie.
3. Efficiëntie: Het algoritme is computatie-efficiënt. Hoewel de complexiteit $O(N^3)$ is (voornamelijk door de berekening van de Bessel-coëfficiënten), kan het gebruik van de inverse Fast Fourier Transform (iFFT) voor de Laplace-inversie de efficiëntie aanzienlijk verbeteren onder bepaalde fysieke aannames.
4. Toepassing op Boei-plaatsing: De resultaten bieden een theoretisch kader om de optimale plaatsing van waarnemingsboeien (zoals het DART-systeem) te bepalen, zodat deze de golfvorm vóór de kust nauwkeurig kunnen voorspellen.

Resultaten

• Validatie: Numerieke experimenten met verschillende initiële profielen (Gaussische, sech-kwadratische en gemengde vormen) tonen een uitstekende overeenkomst tussen de oorspronkelijke en gereconstrueerde randvoorwaarden op de boei.
• Reconstructie van Solitons: In het stuksgewijs hellende model slaagden de auteurs erin om de parameters van twee solitons in het diepe water nauwkeurig te reconstrueren uit de kustdata, zelfs met de aanwezigheid van reflecties. De foutmarges in de positie van de pieken waren verwaarloosbaar klein (bijv. exacte piek $x=119.55$ vs. gereconstrueerd $x=119.50$).
• Vergelijking Modellen: De vergelijking tussen het Boussinesq-model en het lineaire SWE-model toonde aan dat beide methoden bruikbaar zijn, maar dat het Boussinesq-model beter geschikt is voor het modelleren van dispersieve effecten die optreden in dieper water.

Betekenis en Toekomstperspectief

De studie is van groot belang voor de tsunami-mitigatie en vroegtijdige waarschuwing. Door de initiële golfvorm nauwkeurig te kunnen reconstrueren vanuit kustmetingen, kunnen waarschuwingssystemen sneller en accurater worden. Het biedt een wiskundige basis voor het optimaliseren van het netwerk van DART-buoy's.

Beperkingen en Toekomstig Werk:

• Het huidige model is één-dimensionaal (1D). Uitbreiding naar 2D is een logische volgende stap.
• De methode vereist dat de golf niet breekt (geen "gradient catastrophe"). De auteurs plannen onderzoek naar de voorwaarden waaronder de kustdata "toelaatbaar" is (d.w.z. leidt tot een niet-brekende golf).
• Verdere verfijning van de bodemprofielen (meerdere hellingen) en het integreren van resonantie-effecten (versterking van de oploop) worden als toekomstige uitdagingen geïdentificeerd.

Kortom, dit artikel levert een robuust wiskundig raamwerk dat de brug slaat tussen kustobservaties en de initiële tsunami-condities in de diepe oceaan, met name door de combinatie van niet-lineaire dynamica en dispersieve effecten.