IP Library › Granted Patent US 12,650,304
Granted Patent B2
US 12,650,304 · App. 17/989,922 · Granted Jun 9, 2026

Systems and methods for generating a phase-resolved ocean wave forecasts with ensemble based data assimilation

Inventors: Yulin Pan (Ann Arbor, MI); Guangyao Wang (Tianjin, CN); Yuminghao Xiao (Wuhan, CN)
Assignee: THE REGENTS OF THE UNIVERSITY OF MICHIGAN
G01C13/004B63B43/02B63B51/00G01C21/203
View Patent ↗
Loading inventors, assignments & file history…
Monitor This Case
Get email alerts when status or documents change.
Order Certified Copies
Most orders are placed with the USPTO same day — all within 24 business hours.
Order via The Patent Place →
Pre-filled with this patent's details
Quick Facts
Patent No.
US 12,650,304
App. No.
17/989,922
Granted
Jun 9, 2026
Kind
B2
Abstract

Systems and methods for generating a phase-resolved ocean wave forecast with ensemble based data assimilation are disclosed. An example method includes receiving radar data corresponding to an ocean surface, and determining a surface elevation and a surface potential of a portion of the ocean surface. The example method also includes generating an ensemble of perturbed ocean surface data, and applying a phase-resolved nonlinear wave model to the ensemble of perturbed ocean surface data to generate a set of forecast ocean surface data. The example method also includes receiving a subsequent set of radar data corresponding to the ocean surface, and determining a subsequent surface elevation and surface potential of the portion of the ocean surface. The example method also includes combining, by applying an ensemble Kalman filter, the set of forecast ocean surface data with the subsequent surface elevation and surface potential to generate a phase-resolved ocean wave forecast.

Claims (63)

1 . A method for generating a phase-resolved ocean wave forecast with ensemble based data assimilation, the method comprising:

(a) receiving, at one or more processors, a set of radar data corresponding to an ocean surface;

(b) determining, by the one or more processors based on the set of radar data, a surface elevation and a particle velocity of a portion of the ocean surface;

(c) generating, by the one or more processors, an ensemble of perturbed ocean surface data based on the surface elevation and the particle velocity, wherein said ensemble is generated by adding random noise to the surface elevation;

(d) applying, by the one or more processors, a phase-resolved nonlinear wave model to the ensemble of perturbed ocean surface data to generate a set of forecast ocean surface data;

(e) receiving, at the one or more processors, a subsequent set of radar data corresponding to the ocean surface;

(f) determining, by the one or more processors based on the subsequent set of radar data, a subsequent surface elevation and a subsequent particle velocity of the portion of the ocean surface; and

(g) combining, by the one or more processors applying an ensemble Kalman filter, the set of forecast ocean surface data with the subsequent surface elevation and the subsequent particle velocity to generate a phase-resolved ocean wave forecast.

2 . The method of claim 1 , further comprising:

identifying, by the one or more processors, a rogue wave based on the phase-resolved ocean wave forecast; and

generating, by the one or more processors applying a modified route planner, a modified travel route for a ship to avoid the rogue wave.

3 . The method of claim 1 , wherein the method further comprises:

iteratively performing (e)-(g) for the portion of the ocean surface until the portion is not included in the subsequent set of radar data.

4 . The method of claim 1 , wherein the set of radar data includes radar data corresponding to a plurality of points across a region of the ocean surface, the portion of the ocean surface corresponds to a respective point of the plurality of points, and wherein the method further comprises:

performing (b)-(g) for each respective point of the plurality of points.

5 . The method of claim 1 , further comprising:

linearly inflating, by the one or more processors utilizing an adaptive inflation algorithm, the set of forecast ocean surface data.

6 . The method of claim 1 , further comprising:

identifying, by the one or more processors, a first predictable zone and a first unpredictable zone within the set of radar data;

identifying, by the one or more processors, a second predictable zone and a second unpredictable zone within the subsequent set of radar data; and

combining, by the one or more processors applying the ensemble Kalman filter, the set of forecast ocean surface data with the subsequent surface elevation and the subsequent particle velocity to generate the phase-resolved ocean wave forecast based on radar data corresponding to an overlap region between the first predictable zone and the second predictable zone.

7 . The method of claim 1 , wherein the phase-resolved nonlinear wave model comprises a high-order spectral (HOS) method.

8 . The method of claim 1 , wherein the ensemble of perturbed ocean surface data includes a plurality of surface elevations and a plurality of particle velocities, and the method further comprises:

generating, by the one or more processors utilizing a set of measurement error statistics, the ensemble of perturbed ocean surface data based on the surface elevation and the particle velocity.

9 . The method of claim 1 , further comprising:

causing, by the one or more processors, a user device to display the phase-resolved ocean wave forecast on a user interface of the user device for viewing by a user.

10 . A system for generating a phase-resolved ocean wave forecast with ensemble based data assimilation, the system comprising:

a memory storing a set of computer-readable instructions comprising at least a phase-resolved nonlinear wave model and an ensemble Kalman filter; and

a processor interfacing with the memory, and configured to execute the set of computer-readable instructions to cause the processor to:

(a) receive a set of radar data corresponding to an ocean surface;

(b) determine, based on the set of radar data, a surface elevation and a particle velocity of a portion of the ocean surface;

(c) generate an ensemble of perturbed ocean surface data based on the surface elevation and the particle velocity, wherein said ensemble is generated by adding random noise to the surface elevation;

(d) apply the phase-resolved nonlinear wave model to the ensemble of perturbed ocean surface data to generate a set of forecast ocean surface data;

(e) receive a subsequent set of radar data corresponding to the ocean surface;

(f) determine, based on the subsequent set of radar data, a subsequent surface elevation and a subsequent particle velocity of the portion of the ocean surface; and

(g) combine, by applying the ensemble Kalman filter, the set of forecast ocean surface data with the subsequent surface elevation and the subsequent particle velocity to generate a phase-resolved ocean wave forecast.

11 . The system of claim 10 , wherein the set of computer-readable instructions further cause the processor to:

identify a rogue wave based on the phase-resolved ocean wave forecast; and

apply a modified route planner to generate a modified travel route for a ship to avoid the rogue wave.

12 . The system of claim 10 , wherein the set of computer-readable instructions further cause the processor to:

iteratively perform (e)-(g) for the portion of the ocean surface until the portion is not included in the subsequent set of radar data.

13 . The system of claim 10 , wherein the set of radar data includes radar data corresponding to a plurality of points across a region of the ocean surface, the portion of the ocean surface corresponds to a respective point of the plurality of points, and wherein the set of computer-readable instructions further cause the processor to:

perform (b)-(g) for each respective point of the plurality of points.

14 . The system of claim 10 , wherein the set of radar data is collected by a plurality of ships and a plurality of buildings.

15 . The system of claim 10 , wherein the ensemble of perturbed ocean surface data includes a plurality of surface elevations and a plurality of particle velocities, and wherein the set of computer-readable instructions further cause the processor to:

generate, by utilizing a set of measurement error statistics, the ensemble of perturbed ocean surface data based on the surface elevation and the particle velocity.

16 . The system of claim 10 , further comprising a user interface, and wherein the set of computer-readable instructions further cause the processor to:

cause the user interface to display the phase-resolved ocean wave forecast for viewing by a user.

17 . A non-transitory computer-readable storage medium having stored thereon a set of instructions, executable by at least one processor, for generating a phase-resolved ocean wave forecast with ensemble based data assimilation, the instructions comprising:

(a) instructions for receiving a set of radar data corresponding to an ocean surface;

(b) instructions for determining, based on the set of radar data, a surface elevation and a particle velocity of a portion of the ocean surface;

(c) instructions for generating an ensemble of perturbed ocean surface data based on the surface elevation and the particle velocity, wherein said ensemble is generated by adding random noise to the surface elevation;

(d) instructions for applying a phase-resolved nonlinear wave model to the ensemble of perturbed ocean surface data to generate a set of forecast ocean surface data;

(e) instructions for receiving a subsequent set of radar data corresponding to the ocean surface;

(f) instructions for determining, based on the subsequent set of radar data, a subsequent surface elevation and a subsequent particle velocity of the portion of the ocean surface; and

(g) instructions for combining, by applying an ensemble Kalman filter, the set of forecast ocean surface data with the subsequent surface elevation and the subsequent particle velocity to generate a phase-resolved ocean wave forecast.

18 . The non-transitory computer-readable storage medium of claim 17 , wherein the instructions further comprise:

instructions for identifying a rogue wave based on the phase-resolved ocean wave forecast; and

instructions for applying a modified route planner to generate a modified travel route for a ship to avoid the rogue wave.

19 . The non-transitory computer-readable storage medium of claim 17 , wherein the instructions further comprise:

instructions for iteratively performing (e)-(g) for the portion of the ocean surface until the portion is not included in the subsequent set of radar data.

20 . The non-transitory computer-readable storage medium of claim 17 , wherein the set of radar data includes radar data corresponding to a plurality of points across a region of the ocean surface, the portion of the ocean surface corresponds to a respective point of the plurality of points, and wherein the instructions further comprise:

instructions for performing (b)-(g) for each respective point of the plurality of points.

Assignments (1)
ASSIGNMENT OF ASSIGNOR'S INTEREST Recorded Nov 19, 2022
From: PAN, YULIN; WANG, GUANGYAO; XIAO, YUMINGHAO
To: THE REGENTS OF THE UNIVERSITY OF MICHIGAN
Reel/Frame 061833/0345 →
Continuity (2)
Provisional Application 63264289 · Nov 18, 2021
Related Publication 20230152092A1 · May 18, 2023
References Cited (26)
US 5528493A · Potter · 1996 [cited by examiner]
US 5546084A · Hindman · 1996 [cited by examiner]
US 9223058B1 · Uy · 2015 [cited by examiner]
US 10677891B2 · Rudzinsky et al. · 2020 [cited by applicant]
US 20090303109A1 · Erkocevic-Pribic · 2009 [cited by examiner]
US 20120130569A1 · Huntsberger · 2012 [cited by examiner]
US 20170307748A1 · Johnson · 2017 [cited by examiner]
US 20180081028A1 · Kusters, Jr. · 2018 [cited by examiner]
US 20180081054A1 · Rudzinsky · 2018 [cited by examiner]
US 20190161152A1 · Kusters, Jr. · 2019 [cited by examiner]
US 20220155069A1 · Previsic · 2022 [cited by examiner]
US 20240300626A1 · Pan · 2024 [cited by examiner]
CN 105389468B · 2017 [cited by applicant]
WO WO2012117195A1 · 2012 [cited by applicant]
WO WO2020131187A2 · 2020 [cited by applicant]
WO WO2020180818A1 · 2020 [cited by examiner]
Xiao et al., “Time-Optimal Path Planning in an Evolving Ocean Wave Field based on Reachability Theory,” IEEE Transactions and Control Systems Technology, vol. 31 (2023). [cited by applicant]
Wang et al., “Phase-resolved ocean wave forecast with ensemble-based data assimilation,” Journal of Fluid Mechanics, (2021). [cited by applicant]
Waseda et al., “Directional Coherent Wave Group From an Assimilated Non-Linear Wavefield,” Front. Phys., vol. 9 (2021). [cited by applicant]
Yoon et al., “An Explicit Data Assimilation Scheme for a Nonlinear Wave Prediction Model Based on a Pseudo-Spectral Method,” IEEE Journal of Oceanic Engineering, (2016). [cited by applicant]
Fujimoto et al., “Ensemble-Based Variational Method for Nonlinear Inversion of Surface Gravity Waves,” Journal of Atmospheric and Oceanic Technology, vol. 37 (2020). [cited by applicant]
Serpoushan et al., “An Ensemble Kalman Filter Data Assimilation Scheme for Modeling the Wave Climate in Persian Gulf,” Proceedings of the ASME 2013 32nd International Conference on Ocean, Offshore and Arctic Engineering… [cited by applicant]
Zamani et al., Non-Linear Wave Data Assimilation with an ANN-Type Wind-Wave Model and Ensemble Kalman Filter (EnKF), Applied Mathematical Modeling, vol. 34 (2010). [cited by applicant]
Desmars et al., “Experimental and Numerical Assessment of Deterministic Nonlinear Ocean Waves Prediction Algorithms Using Non-Uniformly Sampled Wave Gauges,” Ocean Engineering, (2020). [cited by applicant]
Simpson et al., “Wave-by-Wave Forecasting via Assimilation of Marine Radar Data,” Journal of Atmospheric and Oceanic Technology, vol. 37 (2020). [cited by applicant]
Kollisch et al., Nonlinear real time prediction of ocean surface waves, Ocean Engineering, vol. 157 pp. 387-400 (2018). [cited by applicant]