A method for estimating meteo-marine parameters from satellite Synthetic Aperture Radar (SAR) data, intended for near-real-time (NRT) service over the Baltic Sea, is presented and validated. Total significant wave height data are retrieved with an empirical function CWAVE_S1-IW, which combines spectral analysis of Sentinel-1A/B Interferometric Wide swath (IW) subscenes with wind data derived with common C-Band Geophysical Model Functions (GMFs). In total, 15 Sentinel-1A/B scenes (116 acquisitions) over the Baltic Sea were processed for comparison with off-shore sea state measurements (52 collocations) and coastal wind measurements (357 colocations). Sentinel-1 wave height was spatially compared with WAM wave model results (Copernicus Marine Environment Monitoring Service (CMEMS). The comparison of SAR-derived wave heights shows good agreement with measured wave heights correlation r of 0.88 and with WAM model (r = 0.85). The wind speed estimated from SAR images yields good agreement with in situ data (r = 0.91). The study demonstrates that the wave retrievals from Sentinel-1 IW data provide valuable information for operational and statistical monitoring of wave conditions in the Baltic Sea. The data is valuable for model validation and interpretation in regions where, and during periods when, in situ measurements are missing. The Sentinel-1 A/B wave retrievals provide more detailed information about spatial variability of the wave field in the coastal zone compared to in situ measurements, altimetry wave products and model forecast. Thus, SAR data enables estimation of storm locations and areal coverage. Methods shown in the study are implemented in NRT service in German Aerospace Center’s (DLR) ground station Neustrelitz.

Multi-mission satellite altimetry (e.g., ERS, Envisat, TOPEX/Poseidon, Jason) data have enabled a synoptic-scale view of ocean variations in past decades. Since 2016, the Sentinel-3 mission has provided better spatial and temporal sampling compared to its predecessors. The Sentinel-3 Ku/C Radar Altimeter (SRAL) is one of the synthetic aperture radar altimeters (SAR Altimeter) which is more precise for coastal and lake observations. The article studies the performance of the Sentinel-3 Level-2 sea level altimetry products in the coastal areas of the Baltic Sea and on two lakes of Estonia. The Sentinel-3 data were compared with (i) collocated Global Navigation Satellite System (GNSS) ship measurements, (ii) the Estonian geoid model (EST-GEOID2017) together with sea-level anomaly corrections from the tide gauges, and (iii) collocated buoy measurements. The comparisons were carried out along seven Sentinel-3A/B tracks across the Baltic Sea and Estonian lakes in 2019. In addition, the Copernicus Marine Environment Monitoring Service (CMEMS) Level-3 sea-level products and the Nucleus for European Modelling of the Ocean (NEMO) reanalysis outcomes were compared with measurements from Estonia’s 21 tide gauges and the buoy deployed offshore. Our results showed that the uncertainty of the Sentinel-3 Level-2 altimetry product was below decimetre level for the seacoast and the selected lakes of Estonia. Results from CMEMS Level-3 altimetry products showed a correlation of 0.83 (RMSE 0.18 m) and 0.91 (RMSE 0.27 m) when compared against the tide gauge measurements and the NEMO model, respectively. The overall performance of the altimetry products was very good, except in the immediate vicinity of the coastline and for the lakes, where the accuracy was nearly three times lower than for the open sea, but still acceptably good.

In this work, remote sensing synthetic aperture radar (SAR) data from X-band TerraSAR-X and TanDEM-X (TS-X and TD-X) satellites have been used to adopt the algorithms for estimating sea state parameters in the specific condition of the Baltic Sea with archipelago islands and where short steep sea state dominates. Since the moving targets can be defocused and shifted in SAR images, sea state consisting of short windsea waves with strong local orbital velocities and wave breaking needs additional effort for accurate estimation of the total significant wave height that consists of swell and windsea parts. The XWAVE_C algorithm, developed for the North Sea, where the long swell waves coming from the Atlantic Ocean are present during storms, was further enhanced for the short steep windsea which dominates under ordinary storm conditions in the Baltics. For the empirical XWAVE_C model function, based on the spectral analysis of subscenes as well as on local wind information, an additional term was incorporated for assessment the minimal windsea significant wave height by applying JONSWAP wave spectra. A term to compensate spectral distortions triggered by windsea waves moving in SAR flight direction has also been introduced. In total, 95 TS-X/TD-X StripMap scenes between 2012 and 2017 were acquired in Eastern Baltic Sea, processed and analysed. The wave height results from SAR images were compared with collocated in situ data from 11 available buoys. The analysed data include both high and low windsea conditions. The comparison of SAR-derived wave heights with measured wave heights shows high agreement with a correlation coefficient r of 0.88. The wind speed, estimated from SAR images, was compared to measurements from 14 collocated in situ stations, yielding a high agreement with an r value of 0.90. This article is focused on the algorithm developments; however, it is also the first study of sea state retrieval in the Baltic Sea using high-resolution satellite-based techniques. The results show the local variability in the wave fields connected to atmospheric features. The observed local wave height can increase by 1–2 m in kilometre-size cells that are accompanied by wind gusts. The developed algorithms are installed in the German Aerospace Center’s (DLR) ground station Neustrelitz and can also be used in near-real-time.

A method for reconstructing weekly Baltic gridded sea levels was developed and tested. This method uses input data from tide gauge and altimetry observations. The reconstruction is based on sea level empirical orthogonal function (EOF) modes, calculated as spatiotemporal statistics from daily model reanalysis results available from the Copernicus Marine Service for the 1993–2021 period. In the semi-enclosed, tideless Baltic Sea, the four leading EOF modes cover 99% of the sea level variance. Several experiments with different combinations of input data were carried out. This method was validated against coastal tide gauges and altimetry data. The best reconstruction was obtained when both the tide gauges and altimetry data were used as inputs. An assessment of the centered root-mean-square difference (cRMSD) of the reconstruction versus the tide gauges revealed a value of 0.05 m, and a result of 0.10 m was revealed versus altimetry. The average coefficient of determination (R²) was determined to be 0.93 for the tide gauges and 0.82 for the altimetry data. In the cases where only one type of input data was used, the reconstruction worsened with respect to other data sources. The reconstruction method demonstrated its usefulness for the reconstruction of coastal sea levels in unsampled locations and the calculation of changes in sea volume.

This article proposes to apply long-short-term memory (LSTM) deep learning models to transform Sentinel-1 A/B interferometric wide (IW) swath image data into the wave density spectrum. Although spectral wave estimation methods for synthetic aperture radar data have been developed, similar approaches for coastal areas have not received enough attention. Partially, this is caused by the lack of high-resolution wave-mode data, as well as the nature of wind waves that have more complicated backscattering mechanisms compared to the swell waves for which the aforementioned methods were developed. The application of the LSTM model has allowed the transformation of the Sentinel-1 A/B IW one-dimensional image spectrum into wave density spectra. The best results in the test dataset led to the mean Pearson's correlation coefficient 0.85 for the comparison of spectra and spectra. The result was achieved with the LSTM model using $VV$ and $VH$ polarization spectra fed into the model independently. Experiments with LSTM neural networks that classify images into wave spectra with the Baltic Sea dataset demonstrated promising results in cases where empirical methods were previously considered.

Circularly polarized X-band coastal marine radar data have been used to develop the algorithm for estimating significant wave height (HS) in the Tallinn Bay in the Gulf of Finland directly from radar images. Since sea state in the Tallinn Bay is mainly dominated by slight (WMO-3) windsea, the traditional methods where backscatter intensity variance spectrum is transferred to wave spectrum do not resolve wave height retrievals with sufficient accuracy. In contrast, an empirical method which uses only image spectrum and its parameters has proven to be applicable under these conditions. A wave height retrieval algorithm was developed for the short steep wind sea using image spectrum parameters in addition to Grey level co-occurrence matrix (GLCM) statistics of the radar signal intensity. In total, 1678 collocation pairs from October 18, 2016 to November 14, 2016 were used in the algorithm tuning process. The HS results from radar images were collocated with in situ data from three buoys representing variable meteo-marine conditions. The comparison of radar-derived HS with measured HS shows high agreement with a correlation coefficient r of 0.78 (RMSE--0.23 m) for tuning data set. The method validation with independent data sets from January and June showed high correlation values of 0.82 (RMSE--0.26 m) and 0.89 (RMSE--0.25 m) correspondingly. In case of HS > 0.5 m, the validation resulted in higher correlation (over 0.93) and lower RMSE (from 0.15 to 0.21 m). The spatial variability of wave height from radar imagery is demonstrated based on the commonly occurring north-western storm.