Ocean surface wave patterns on TerraSAR-X images and follow-on applications

The TerraSAR-X (TS-X) satellite provides images of wave patterns at the ocean surface which are of great scientific value. With the relatively low orbit of the platform, SAR specific non-linear imaging effects of the moving sea surface play a minor role. As a result, TS-X is able to image ocean waves down to around ~30 m length. Data on the ocean wavelength and wave direction are inferred directly from the image spectra, i.e. without inversion techniques. In the Stripmap mode, TS-X image sequences may extend in length up to several hundred kilometres. Detecting the spatial variation of sea state parameters is the basis for different follow-on applications. We present selected use cases including the validation of sea state forecast models and a recent study on wave properties in the marginal ice zone. Another potential application included in our contribution, is the estimation of coastal bathymetry from ocean wavelength variations

Investigating High-Resolution Spatial Wave Patterns on the Canadian Beaufort Shelf Using Synthetic Aperture Radar Imagery at Herschel Island, Qikiqtaruk, Yukon, Canada

The Arctic is experiencing the greatest increase in air temperature on Earth. This significant climatic change is leading to a significant positive trend of increasing wave heights and greater coastal erosion. This in turn effects local economies and ecosystems. Increasing wave energy is one of the main drivers of this alarming trend. However, the data on spatial and temporal patterns of wave heights in the Arctic are either coarse, interpolated or limited to point measurements. The aim of this study is to overcome this shortcoming by using remote sensing data. In this study, the Synthetic Aperture Radar (SAR) satellite TerraSAR-X (TS-X) and TanDEM-X (TD-X) imagery are used to obtain sea state information with a high spatial resolution in Arctic nearshore waters in the Canadian Beaufort Sea. From the entire archive of the TS-X/TD-X StripMap mode with coverage around 30 km × 50 km acquired between 2009 and 2020 around Herschel Island, Qikiqtaruk (HIQ), all the ice-free scenes were processed. The resulting dataset of 175 collocated scenes was used to map the significant wave height (����) and to link spatial and temporal patterns to local coastal processes. Sea state parameters are estimated in raster format with a 600 m step using the empirical algorithm CWAVE_EX. The statistics of the ���� were aggregated according to spatial variability, seasonality and wind conditions. The results show that the spatial wave climate is clearly related to the dominant wind regime and seasonality. For instance, the aggregation of all the scenes recorded in July between 2009 and 2020 results in an average of 0.82 m ����, while in October the average ���� is almost 0.40 m higher. The analysis by wind direction shows that fetch length and wind speed are likely the most important variables influencing the spatial variability. A larger fetch under NW conditions results in a mean wave height of 0.92 m, while waves generated under ESE conditions are lower at 0.81 m on average

Sea State from High Resolution Satellite-borne Synthetic Aperture Radar Imagery

The Sea Sate Processor (SSP) was developed for fully automatic processing of high-resolution Synthetic Aperture Radar (SAR) data from TerraSAR-X (TS-X) satellites and implemented into the processing chain for Near Real Time (NRT) services in the DLR Ground Station "Neustrelitz". The NRT chain was organised and tested to provide the processed data to the German Weather Service (DWD) in order to validate the new coastal forecast model CWAM (Coastal WAve Model) in the German Bight of the North Sea with 900 m horizontal resolution. The NRT test-runs, wherein the processed TS-X data were transferred to DWD and then incorporated into forecast products reach the best performance about 10 min for delivery of processed TS-X data to DWD server after scene acquisition. To do this, a new empirical algorithm XWAVE_C (C = coastal) for estimation of significant wave height from X-band satellite-borne SAR data has been designed for coastal applications. The algorithm is based on the spectral analysis of subscenes and the empirical model function yields an estimation of integrated sea state parameters directly from SAR image spectra without transformation into wave spectra. To provide the raster coverage analysis, the SSP intends three steps of recognising and removing the influence of non-sea-state-produced signals in the Wadden Sea areas such as ships, buoys, dry sandbars as well as nonlinear SAR image distortions produced by e.g. short and breaking waves. For the validation, more than 150 TS-X StripMap scene sequences with a coverage of ~30 km × 300 km across the German Bight since 2013 were analysed and compared with in situ Buoy measurements from 6 different locations. On this basis, the SSP autonomous processing of TS-X Stripmap images has been confirmed to have a high accuracy with an error RMSE = 25 cm for the total significant wave height

Comparing wavelengths simulated by the coastal wave model CWAM and TerraSAR-X satellite data

The accuracy of the high resolution coastal wave forecast model CWAM is validated on the basis of sea state information from satellite images of TerraSAR-X (TS-X). Employing 2-dimensional Fast Fourier Transformation, image spectra are derived from TS-X and locally varying patterns of the peak wavelengths are provided. Subsequently, wavelength comparisons are performed between a typical set of TS-X scenes acquired in December 2013 over the German Bight and the model hindcasts. This results mostly in reasonable agreement. Potential wave modelling errors are discussed as well

Multiparametric Sea State from Spaceborne Synthetic Aperture Radar for Near Real Time Services

Spaceborne synthetic aperture radar (SAR) is a powerful tool for monitoring seas. The ability to work independently of sun illumination, cloud coverage and atmospheric conditions, as well as the capability of delivering spatial information, makes SAR one of the most perceptive instruments. The newest methods for processing SAR data with increased precision allow sea state fields to be estimated with local variabilities. For large areas in oceans where no in-situ measurements and only forecast predictions are available, this information is indispensable for global shipping and over human activity. Due to newest developments, the derived meteo-marine parameters can be transferred to weather services and to a ship’s bridge several minutes after acquisition, where the ship route can be optimized. The study presents a method and application for estimating series of integrated sea state parameters from satellite-borne SAR, allow processing of data from different satellites and modes in near real time (NRT). The developed Sea State Processor (SSP) estimates total significant wave height Hs, dominant and secondary swell and windsea wave heights, first, and second moment wave periods, mean wave period and period of wind sea. The algorithm was applied for the Sentinel-1 (S1) C-band Interferometric Wide Swath Mode (IW), Extra Wide (EW) and Wave Mode (WM) Level-1 (L1) products and also extended to the Xband TerraSAR-X (TSX) StripMap (SM) mode. The scenes are processed in raster and result in continuous sea state fields with the exception of S1 WV. Each 20 km × 20 km WV imagette, acquired every 100 km along the orbit, presents averaged values for each sea state parameter. The SSP was tuned and validated using two independent global wave models WAVEWATCH-3 (NOAA) and CMEMS (Copernicus) and NDBC buoys. The accuracy of Hs reaches an RMSE of 0.25 m by comparison with models (S1 WV); comparisons to NDBC worldwide buoys result into an RMSE of 0.3 m. Due to implemented parallelization, a fine rater step for scene processing can be practical applied: for example, S1 IW scene with coverage of 200 km × 250 km can be processed using raster step of 1 km (corresponds to ~50.000 subscenes) during minutes. The DLR Ground Station “Neustrelitz” applies SSP as part of a near real-time demonstrator service that involves a fully automated daily provision of surface wind and sea state parameters estimates from S1 IW for the North and Baltic Sea. All results and the presented methods are novel and provide a wide field for applications and implementations in prediction systems

Meteo-Marine Parameters from Sentinel-1 SAR Imagery: Towards Near Real-Time Services for the Baltic Sea

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