skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Data Assimilation for Tsunami Forecast With Ship‐Borne GNSS Data in the Cascadia Subduction Zone
Abstract An efficient and cost‐effective near‐field tsunami warning system is crucial for coastal communities. The existing tsunami forecasting system is based on offshore Deep‐Ocean Assessment and Reporting of Tsunamis and Global Navigation Satellite System (GNSS) buoys which are not affordable for many countries. A potential cost‐effective solution is to utilize position data from ships traveling in coastal and offshore regions. In this study, we examine the feasibility of using ship‐borne GNSS data in tsunami forecasting. We carry out synthetic experiments by applying a data assimilation (DA) method with ship position (elevation and velocity) data. Our findings show that the DA method can recover the reference model with high accuracy if a dense network of ship elevation data is used. However, the use of ship velocity data alone is unable to recover the reference model. In addition, we carried out sensitivity studies of the DA method to the ship spatial distribution. We find that a 20 km gap between the ships works well in terms of accuracy and computational time for the example source model that we explored. The highest accuracy is obtained when data from a sufficient number of ships traveling in and around the tsunami source area are available.  more » « less
Award ID(s):
1855090
PAR ID:
10369043
Author(s) / Creator(s):
 ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Earth and Space Science
Volume:
8
Issue:
3
ISSN:
2333-5084
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Precise point positioning (PPP) of ships using Global Navigation Satellite System (GNSS) data reveals the precise movements of marine vessels. This method may quantify anomalies in sea surface height with implications for oceanographic monitoring, exploration, and tsunami warning. The GNSS PPP data from theR/V Sikuliaq, a research ship of the University of Alaska Fairbanks, were processed to detect a small local tsunami generated by the Lowell Point landslide, which occurred near Seward, Alaska, on 8 May 2022 (UTC). The GNSS receiver aboard theR/V Sikuliaqrecorded the waves generated by the landslide, with a maximum wave amplitude of 6 cm and wave periods between 40 and 50 s. These results are consistent with simulations of the landslide event. 
    more » « less
  2. Tsunamigenic megathrust earthquakes along the Cascadia subduction zone present a major hazard concern. We can better prepare to model the earthquake source in a rapid manner by imbuing fault geometry constraints based on prior knowledge and by evaluating the capabilities of using existing GNSS sensors. Near-field GNSS waveforms have shown promise in providing rapid coarse finite-fault model approximations of the earthquake rupture that can improve tsunami modeling and response time. In this study, we explore the performance of GNSS derived finite-fault inversions and tsunami forecasting predictions in Cascadia that highlights the impact and potential of geodetic techniques and data in operational earthquake and tsunami monitoring. We utilized 1300 Cascadia earthquake simulations (FakeQuakes) that provide realistic (M7.5-9.3) rupture scenarios to assess how feasibly finite-fault models can be obtained in a rapid earthquake early warning and tsunami response context. A series of fault models with rectangular dislocation patches spanning the Cascadia megathrust area is added to the GFAST inversion algorithm to calculate slip for each earthquake scenario. Another method used to constrain the finite-fault geometry is from the GNSS-derived CMT fault plane solution. For the Cascadia region, we show that fault discretization using two rectangular segments approximating the megathrust portion of the subduction zone leads to improvements in modeling magnitude, fault slip, tsunami amplitude, and inundation. In relation to tsunami forecasting capabilities, we compare coastal amplitude predictions spanning from Vancouver Island (Canada) to Northern California (USA). Generally, the coastal amplitudes derived using fault parameters from the CMT solutions show an overestimation bias compared to amplitudes derived from the fixed slab model. We also see improved prediction values of the run-up height and maximum amplitude at 10 tide gauge stations using the fixed slab model as well. 
    more » « less
  3. We investigate the potential of using Global Navigation Satellite System (GNSS) observations to directly forecast full tsunami waveforms in real time. We train convolutional neural networks to use less than 9 min of GNSS data to forecast the full tsunami waveforms over 6 hr at select locations, and obtain accurate forecasts on a test data set. Our training and test data consists of synthetic earthquakes and associated GNSS data generated for the Cascadia Subduction Zone using the MudPy software, and corresponding tsunami waveforms in Puget Sound computed using GeoClaw. We use the same suite of synthetic earthquakes and waveforms as in earlier work where tsunami waveforms were used for forecasting, and provide a comparison. We also explore varying the number of GNSS stations, their locations, and their observation durations. 
    more » « less
  4. The paper discusses an intelligent vision-based control solution for autonomous tracking and landing of Vertical Take-Off and Landing (VTOL) capable Unmanned Aerial Vehicles (UAVs) on ships without utilizing GPS signal. The central idea involves automating the Navy helicopter ship landing procedure where the pilot utilizes the ship as the visual reference for long-range tracking; however, refers to a standardized visual cue installed on most Navy ships called the ”horizon bar” for the final approach and landing phases. This idea is implemented using a uniquely designed nonlinear controller integrated with machine vision. The vision system utilizes machine learning based object detection for long-range ship tracking and classical computer vision for the estimation of aircraft relative position and orientation utilizing the horizon bar during the final approach and landing phases. The nonlinear controller operates based on the information estimated by the vision system and has demonstrated robust tracking performance even in the presence of uncertainties. The developed autonomous ship landing system was implemented on a quad-rotor UAV equipped with an onboard camera, and approach and landing were successfully demonstrated on a moving deck, which imitates realistic ship deck motions. Extensive simulations and flight tests were conducted to demonstrate vertical landing safety, tracking capability, and landing accuracy. The video of the real-world experiments and demonstrations is available at this URL. 
    more » « less
  5. SUMMARY Tsunami generation by offshore earthquakes is a problem of scientific interest and practical relevance, and one that requires numerical modelling for data interpretation and hazard assessment. Most numerical models utilize two-step methods with one-way coupling between separate earthquake and tsunami models, based on approximations that might limit the applicability and accuracy of the resulting solution. In particular, standard methods focus exclusively on tsunami wave modelling, neglecting larger amplitude ocean acoustic and seismic waves that are superimposed on tsunami waves in the source region. In this study, we compare four earthquake-tsunami modelling methods. We identify dimensionless parameters to quantitatively approximate dominant wave modes in the earthquake-tsunami source region, highlighting how the method assumptions affect the results and discuss which methods are appropriate for various applications such as interpretation of data from offshore instruments in the source region. Most methods couple a 3-D solid earth model, which provides the seismic wavefield or at least the static elastic displacements, with a 2-D depth-averaged shallow water tsunami model. Assuming the ocean is incompressible and tsunami propagation is negligible over the earthquake duration leads to the instantaneous source method, which equates the static earthquake seafloor uplift with the initial tsunami sea surface height. For longer duration earthquakes, it is appropriate to follow the time-dependent source method, which uses time-dependent earthquake seafloor velocity as a forcing term in the tsunami mass balance. Neither method captures ocean acoustic or seismic waves, motivating more advanced methods that capture the full wavefield. The superposition method of Saito et al. solves the 3-D elastic and acoustic equations to model the seismic wavefield and response of a compressible ocean without gravity. Then, changes in sea surface height from the zero-gravity solution are used as a forcing term in a separate tsunami simulation, typically run with a shallow water solver. A superposition of the earthquake and tsunami solutions provides an approximation to the complete wavefield. This method is algorithmically a two-step method. The complete wavefield is captured in the fully coupled method, which utilizes a coupled solid Earth and compressible ocean model with gravity. The fully coupled method, recently incorporated into the 3-D open-source code SeisSol, simultaneously solves earthquake rupture, seismic waves and ocean response (including gravity). We show that the superposition method emerges as an approximation to the fully coupled method subject to often well-justified assumptions. Furthermore, using the fully coupled method, we examine how the source spectrum and ocean depth influence the expression of oceanic Rayleigh waves. Understanding the range of validity of each method, as well as its computational expense, facilitates the selection of modelling methods for the accurate assessment of earthquake and tsunami hazards and the interpretation of data from offshore instruments. 
    more » « less