Satellite Altimeter Data Cuts North Sea Wave Model Errors by Twenty Percent

Every day, a fleet of satellites sweeps across the North Sea, bouncing radar pulses off the ocean surface and measuring the height of the waves below. Researchers in the Netherlands have now shown that feeding these measurements into a regional wave model through a sophisticated statistical technique called the Deterministic Ensemble Kalman Filter, or DEnKF, can substantially sharpen the accuracy of wave predictions. In a three-month experiment covering the winter of 2021 to 2022, a period that included the severe storms Corrie and Malik, the assimilation system reduced the error in predicted significant wave height by just over twenty percent at every one of the twenty-four independent validation buoy stations used in the study. The findings, published in Ocean Dynamics, mark an important step toward bringing ensemble-based data assimilation, long a staple of global wave forecasting, into the smaller and more challenging domain of shelf seas.

The work was carried out by C.W.E. de Korte, M. Verlaan, A. W. Heemink and B. Backeberg, affiliated with Delft University of Technology and the research institute Deltares. Their starting point was a persistent problem: third-generation wave models such as SWAN, the Simulating WAves Nearshore model used in the study, are highly reliable in the open ocean but continue to struggle in coastal shelf seas. Wave-current interactions, uncertain wind forcing, shallow-water effects and imperfect parametrizations of the physical source terms all introduce errors that are difficult to eliminate by calibration alone. Data assimilation offers a different route, blending real observations with the model’s own physics to nudge the simulated ocean state closer to reality without rewriting the underlying equations.

What sets this study apart from most earlier wave data assimilation efforts is the choice of state variable. Conventional operational schemes, such as Optimal Interpolation and three-dimensional variational methods, typically apply corrections only to significant wave height and then scale those corrections back onto the full wave spectrum using simplifying assumptions about how wave energy is distributed across frequencies and directions. The Dutch team instead placed the complete directional wave energy spectrum in the model state. Each ensemble member carried the spectrum across 32 frequency bands and 36 directional bins at every point of a 567-cell grid, producing a state vector of more than 650,000 elements per member. Because the ensemble evolves under the full SWAN physics, the corrections the filter produces are automatically physically consistent with the model, and integral parameters such as mean wave period adjust themselves without any ad hoc scaling.

The DEnKF itself is a deterministic variant of the classic Ensemble Kalman Filter. Rather than perturbing observations with random noise to propagate uncertainty, it updates the ensemble mean and the ensemble anomalies separately, avoiding the sampling errors that stochastic perturbations introduce. This is particularly valuable for small ensembles, and the team settled on 64 members after previous synthetic twin experiments showed the error statistics fully converged at that size. Uncertainty was injected into the system through the wind forcing, treated as the control variable, using a first-order autoregressive noise model with a spatial Gaussian correlation structure. Parameters were derived from the difference between HARMONIE wind analyses and forecasts, giving a standard deviation of two metres per second, a decorrelation timescale of fifteen hours and a spatial decorrelation length of 500 kilometres.

The observations came from seven nadir satellite altimeters: CFOSAT, Haiyang-2B, Cryosat-2, Jason-3, the two Sentinel-3 satellites, and Saral/AltiKa, all retrieved from the Copernicus Marine Environment Monitoring Service. Over the three-month window the satellites contributed 713 tracks over the North Sea, an average of about eight passes per day, with a mean interval of roughly three hours between passes but gaps stretching to nearly fifteen hours. Tracks were sub-sampled every 120 kilometres to avoid overloading individual grid cells, and a coastal mask excluded measurements within 50 kilometres of shore, where altimeter retrievals are known to be unreliable. Observation errors were assumed to be uncorrelated with a standard deviation of 0.2 metres. Hamill localisation experiments comparing the standard EnKF with the DEnKF across localisation radii of 100 to 500 kilometres showed the DEnKF with a 200-kilometre radius performed best, and that configuration became the final set-up.

Validation against the North Sea’s dense network of independent wave buoys delivered strikingly consistent results. Significant wave height errors dropped by a mean of 20.5 percent, from 0.39 metres in the free-running coarse model to 0.31 metres, a performance essentially matching the much finer SWAN-DCSM benchmark model run at roughly 3.6-kilometre resolution. Mean wave period improved at 21 of 23 stations with a ten percent reduction in root mean square error, while the peak period improved at 13 of 16 stations by about five percent. Wind speed showed modest improvements at some stations, though the researchers caution that the station anemometers, corrected to ten-metre equivalent heights assuming a neutral wind profile, carry their own uncertainties over the frequently non-neutral marine boundary layer. Not every parameter benefited: swell wave height degraded slightly on average, and the low-frequency inverse moment period and mean wave direction, each measured at only a handful of stations, also worsened marginally.

Spectral analysis explained the pattern. In unimodal sea states dominated by wind-driven waves, the assimilation corrected the entire wave spectrum in a way that closely matched buoy observations, as demonstrated during a storm peak on 20 January 2022 at the offshore station A121, where the analysis tracks the measured spectra hour by hour. But in mixed sea states where wind-sea and swell are clearly separated, the ensemble spread remained concentrated in the mid and high frequencies, because wind perturbations barely touch an independently propagating swell field and the altimeters measure only total significant wave height. Detailed examination of the largest swell errors revealed two distinct mechanisms: during short-fetch, rapidly rotating local wind conditions, the wind-based error covariances failed to represent the spatial scales at which swell actually varied between neighbouring stations, while a second error type, premature swell arrival at coastal stations, proved to be a systematic bias of the coarse-resolution model rather than a failure of the assimilation itself.

One of the most practically important findings concerns timing. The researchers binned all validation samples by the number of hours elapsed since the last satellite pass and found that prediction errors rose steadily with the length of the gap, particularly for stations in the open central North Sea. Coastal stations, whose errors are dominated by shallow-water processes rather than wind-driven corrections, were less sensitive to the satellite schedule. Because the orbital geometry of the contributing satellites fixes the timing of the gaps, these gaps recur with the tidal cycle, a phase-locking effect the authors flag as deserving further study. The message for forecasters is clear: the temporal density of observations matters, and merging additional data sources could deliver substantial gains.

The authors are candid about the limitations. The coarse 0.5-degree grid, chosen so that a 64-member ensemble of full spectra could be run at all, degrades accuracy near the coast, where resolution, missing triad interactions and a simplified setup all take their toll. Running a high-resolution model within the ensemble framework would require major advances in computing power and memory. Still, the path forward is mapped out: refining the wind noise model, assimilating additional integral wave parameters or even full spectra, incorporating continuous buoy measurements, and adding satellite SAR observations from missions such as Sentinel-1, SWOT and Sentinel-6. Beyond operational forecasting, the team points to wave reanalyses for risk and climate studies, and to the growing demand for high-quality training data for machine-learning wave models. For a shelf sea as busy and economically vital as the North Sea, better wave information from the satellites already overhead is a prize worth the computation.

Subject of Research: Ensemble-based assimilation of satellite altimeter wave measurements in a regional North Sea wave model

Article Title: Exploring the added value of assimilating satellite altimeter measurements in a North Sea wave model using the Deterministic Ensemble Kalman Filter

Article References: de Korte, C., Verlaan, M., Heemink, A. W., & Backeberg, B. (2026). Exploring the added value of assimilating satellite altimeter measurements in a North Sea wave model using the Deterministic Ensemble Kalman Filter. Ocean Dynamics, 76(10), Article 102. https://doi.org/10.1007/s10236-026-01858-9

Image Credits: AI Generated

DOI: 10.1007/s10236-026-01858-9

Keywords: wave data assimilation, Deterministic Ensemble Kalman Filter, satellite altimeter, SWAN wave model, North Sea, significant wave height, wave spectrum, ensemble forecasting, swell, coastal shelf seas, Ocean Dynamics, CMEMS