# A new airborne system for simultaneous high-resolution ocean vector current and wind mapping: first demonstration of the SeaSTAR mission concept in the macrotidal Iroise Sea

David Lewis Mccann, Adrien C. H. Martin, Karlus Macedo, Ruben Carrasco Alvarez, Jochen Horstmann, Louis Marié, José Márquez-Martínez, Marcos Portabella, Adriano Meta, Christine Gommenginger, et al.

<figure><img src="image_1.png" /></figure>

## To cite this version:

David Lewis Mccann, Adrien C. H. Martin, Karlus Macedo, Ruben Carrasco Alvarez, Jochen Horstmann, et al. A new airborne system for simultaneous high-resolution ocean vector current and wind mapping: first demonstration of the SeaSTAR mission concept in the macrotidal Iroise Sea. Ocean Science, 2024, 20 (5), pp.1109-1122. 10.5194/os-20-1109-2024. hal-04731661

HAL Id: hal-04731661

https://hal.science/hal-04731661v1

Submitted on 13 Oct 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

<figure><img src="image_4.png" /></figure>

# A new airborne system for simultaneous high-resolution ocean vector current and wind mapping: first demonstration of the SeaSTAR mission concept in the macrotidal Iroise Sea

David L. McCann<fnref n="1" />, Adrien C. H. Martin<fnref n="1" />,<fnref n="2" />, Karlus A. C. de Macedo<fnref n="3" />, Ruben Carrasco Alvarez<fnref n="4" />, Jochen Horstmann<fnref n="4" />, Louis Marié<fnref n="5" />, José Márquez-Martínez<fnref n="6" />, Marcos Portabella<fnref n="7" />, Adriano Meta<fnref n="3" />, Christine Gommenginger<fnref n="1" />, Petronilo Martin-Iglesias<fnref n="8" />, and Tania Casal<fnref n="8" />

**Correspondence:** David L. McCann (david.mccann@noc.ac.uk)

Received: 21 December 2023 – Discussion started: 4 January 2024

Revised: 19 July 2024 – Accepted: 26 July 2024 – Published: 9 September 2024

**Abstract.** Coastal seas, shelf seas and marginal ice zones are dominated by small-scale ocean surface dynamic processes that play a vital role in the transport and exchange of climate-relevant properties such as carbon, heat, water and nutrients between land, ocean, ice and atmosphere. Mounting evidence indicates that ocean scales below 10 km have far-ranging impacts on air–sea interactions, lateral ocean dispersion, vertical stratification, ocean carbon cycling and marine productivity – governing exchanges across key interfaces of the Earth system, the global ocean, and atmosphere circulation and climate. Yet, these processes remain poorly observed at the fine spatial and temporal scales necessary to resolve them. The Ocean Surface Current Airborne Radar (OSCAR) is a new airborne instrument with the capacity to inform these questions by mapping vectorial fields of total ocean surface currents and winds at high resolution over a wide swath. Developed for the European Space Agency (ESA), OSCAR is the airborne demonstrator of the satellite mission concept SeaSTAR, which aims to map total surface current and ocean wind vectors with unprecedented accuracy, spatial resolution and temporal revisit across all coastal seas, shelf seas and marginal ice zones. Like SeaSTAR, OSCAR is an active

microwave synthetic aperture radar along-track interferometer (SAR-ATI) with optimal three-azimuth sensing enabled by unique highly squinted beams. In May 2022, OSCAR was flown over the Iroise Sea, France, in its first scientific campaign as part of the ESA-funded SEASTARex project. The campaign successfully demonstrated the capabilities of OSCAR to produce high-resolution 2D images of total surface current vectors and near-surface ocean vector winds, simultaneously, in a highly dynamic, macrotidal coastal environment. OSCAR current and wind vectors show excellent agreement with ground-based X-band-radar-derived surface currents, numerical model outputs and NovaSAR-1 satellite SAR imagery, with root mean square differences from the X-band radar better than $0.2 \text{ m s}^{-1}$ for currents at 200 m resolution. These results are the first demonstration of simultaneous retrieval of total current and wind vectors from a high-squint three-look SAR-ATI instrument and the first geophysical validation of the OSCAR and SeaSTAR observing principle. OSCAR presents a remarkable new ocean observing capability to support the study of small-scale ocean dynamics and air–sea interactions across the Earth’s coastal, shelf and polar seas.

<footnote n="1">¹National Oceanography Centre, Liverpool, L55DA, UK</footnote>

<footnote n="2">²NOVELTIS, Labège, 31670, France</footnote>

<footnote n="3">³MetaSensing BV, Noordwijk, 2201 DK, the Netherlands</footnote>

<footnote n="4">⁴Department of Radar Hydrography, Helmholtz-Zentrum Hereon, Geesthacht, 21502, Germany</footnote>

<footnote n="5">⁵Laboratoire d'Océanographie Physique et Spatiale (LOPS), Ifremer, Plouzané, 29280, France</footnote>

<footnote n="6">⁶Radarmetrics SL, Santander, 39005, Spain</footnote>

<footnote n="7">⁷Physical and Technological Oceanography, Institut de Ciències del Mar, Barcelona, 08003, Spain</footnote>

<footnote n="8">⁸European Space Research and Technology Centre, European Space Agency, Noordwijk, 2201 AZ, the Netherlands</footnote>
# 1 Introduction

The Earth's coastal and shelf seas and marginal ice zones are dominated by dynamic surface processes that exist on much shorter spatial and temporal scales than in the open ocean (Schulz-Stellenfleth and Stanev, 2016; Kozlov et al., 2019). These small-scale dynamics (order of 1–10 km) play a key role in the exchange and transport of essential ocean and climate variables, biogeochemical processes, vertical mixing, and air–sea interactions (Martin and Richards, 2001; Levy et al., 2012; Sasaki et al., 2014); however global observations of dynamics at these scales are rare (McWilliams, 2016). Ocean surface currents and winds and their resultant air–sea interactions show increased complexity in coastal seas (Bricheno et al., 2013), shelf seas and in the vicinity of sea ice. The limited spatial resolution of wind predictions, model forcing and validation data has direct impacts on the ability to predict storms (Mass et al., 2002; Maskell, 2012) and assess wind power resources in developing areas (Ruiz et al., 2022). Wind observations are particularly lacking in coastal seas where satellite scatterometry is unable to make useful measurements within ~10 km of the coast (Martin et al., 2018). Simultaneous measurements of total current and wind vectors at fine scale are needed to improve understanding and predictive capability of coastal and shelf sea processes (Samelson, 2019; Villas Boas et al., 2019; Hauser et al., 2023). With ocean prediction systems aspiring to resolve small-scale processes operationally across the Earth's coastal, shelf and polar seas by the end of the decade (Holt et al., 2017), high-resolution mapping of these key physical properties will become increasingly important to validate and improve the numerical models used for short-term forecasting and climate change projections. There is therefore a clear and present need for synoptic measurements of currents and winds to meet the scientific, economic and societal challenges of the future.

SeaSTAR is a satellite mission concept submitted to the European Space Agency (ESA) Earth Explorer programme to study these fast-changing, small-scale ocean surface dynamics across the world's coastal oceans, shelf seas and marginal ice zones (Gommenginger et al., 2019). SeaSTAR is a synthetic aperture radar along-track interferometer (SAR-ATI) applied in three highly squinted directions in azimuth. SAR-ATI measures the Doppler shift in ocean microwave backscatter signals from pairs of complex SAR images separated by a short time lag. The measured Doppler shift (interferogram) relates to the motion of the ocean surface in the radar line of sight, which contains the total ocean surface current and the motion of microwave scatterers on the ocean surface (known as the wave Doppler or wind-wave-induced artefact surface velocity, WASV; Martin et al., 2016). The total surface current vector (TSCV), corresponding to the effective mass transport, contains contributions from surface wind drift, Ekman currents, geostrophic currents and Stokes drift but is dominated by tidal currents in many coastal environments by an order of magnitude. TSCV

is derived from the surface motion directly measured by the Doppler shift minus the WASV.

SAR-ATI has a long history, from the first broadside-only ATI measurements of ocean surface velocity by Goldstein and Zebker (1987) through to the two-look Wavemill concept (Buck, 2005) and the SeaSTAR three-look concept (Martin et al., 2018; Gommenginger et al., 2019). To date, no three-look SAR-ATI system has ever been flown in space. SeaSTAR represents a ground-breaking ocean observing system that would apply the cutting-edge performance and spatial resolution of SAR-ATI to deliver simultaneous fine-scale images of total surface current and wind vectors to inform these vital but poorly observed ocean processes.

The Ocean Surface Current Airborne Radar (OSCAR) instrument is a unique three-look SAR-ATI airborne system operating at Ku-band (13.5 GHz). Built and operated by MetaSensing BV in Italy and the Netherlands, OSCAR was developed for ESA to demonstrate the three-look SAR-ATI principle and the ability to simultaneously observe total surface current vectors (TSCVs) and ocean surface vector winds (OSVWs) in two dimensions, at high accuracy and fine resolution, in a single pass from an airborne platform (Trampuz et al., 2018). Within ESA Earth Explorer 11 Phase 0 activities, the SEASTARex project carried out OSCAR airborne campaigns to increase the scientific readiness level of the SeaSTAR mission concept and act as the first scientific demonstration and validation of the OSCAR instrument and the high-squint three-look SAR-ATI capability. Involving airborne OSCAR acquisitions, in situ observations, ground-based remote sensing, numerical modelling and earth observation, SEASTARex drew together expertise from nine international scientific institutions and specialists in the field of oceanography, airborne- and ground-based remote sensing, engineering, and earth observation sciences. This paper presents the activities and results from the scientific airborne campaign of the OSCAR instrument over the Iroise Sea in May 2022, representing the first instance of simultaneous retrieval of total surface current vectors and ocean surface vector winds based on experimental measurements from a high-squint three-look airborne SAR-ATI instrument.

# 2 OSCAR airborne system and data

## 2.1 OSCAR instrument

OSCAR is a unique high-squint three-look SAR-ATI instrument operating at Ku-band (13.5 GHz), featuring two pairs of interferometric antenna squinted at ±45° ("fore" and "aft") from the aircraft broadside and a conventional "zero-Doppler" antenna pointing to broadside ("middle"). All antennas pointed to the port (left) side of the aircraft and transmitted and received in vertical polarization ("VV"). The instrument is mounted on a three-axis gimbal (Fig. 1a) with a mounted inertial measurement unit (IMU) paired with a
high-precision global navigation satellite system (GNSS) receiver to stabilize the instrument pointing relative to the aircraft and resolve the pointing of the beams to a high precision. The antennas and gimbal are mounted in a purpose-built radome specifically designed for interferometry at Ku-band and Ka-band frequencies. OSCAR was installed on a Piper PA-31 Navajo airframe owned and flown by MetaSensing BV (Fig. 1b).

## 2.2 Data processing: simultaneous current and wind vector retrieval

The simultaneous retrieval of total surface current vector (TSCV) and ocean surface vector wind (OSVW) was performed using the SeaSTAR project in Python (Martin et al., 2023), which is based on the simultaneous wind-current retrieval method of Martin et al. (2018), adapted for the OSCAR three-look configuration. At the centre of the simultaneous wind-current retrieval method is the minimization of the cost function $J$, defined here as

$$ J(\mathbf{u}_{10}, c) = \frac{1}{N_S + N_D} \sum_{i=1}^{N_S} \left( \frac{\text{KuMod}(\mathbf{u}_{10}, \chi_i, \theta_i, p_i) - \sigma_{\text{obs},i}^0}{\Delta\sigma_i^0} \right)^2 + \frac{1}{N_S + N_D} \times \sum_{i=1}^{N_D} \left( \frac{\text{KuDop}(\mathbf{u}_{10}, \chi_i, \theta_i, p_i) + c_{\parallel,i} - \text{RSV}_{\text{obs},i}}{\Delta\text{RSV}_i} \right)^2 \quad (1) $$

where $i$ is the beam index in a given azimuth direction (fore, middle, aft), $\sigma_{\text{obs},i}^0$ is the observed normalized radar cross section (NRCS) in VV polarization for the beam index $i$, $\text{RSV}_{\text{obs},i}$ is the measured radial surface velocity for the beam index $i$ (in this case only fore and aft), $\mathbf{u}_{10}$ is the stress-equivalent OSVW at 10 m height (de Kloe et al., 2017), $c$ is TSCV, $N_S$ is the total number of observations for NRCS, $N_D$ is the total number of observations for RSV, KuMod is the predicted NRCS using the chosen geophysical model function (GMF) for Ku-band, KuDop is the predicted RSV using the chosen GMF for Ku-band WASV, $\chi$ is the azimuth look direction, $\theta$ is the incidence angle from nadir, $p$ is the radiometric polarization (in this case VV) and $c_{\parallel,i}$ is the component of TSCV along the azimuth look direction for beam index $i$.

It should be noted that aside from the adaption for three simultaneous looks, Eq. (1) differs from Martin et al. (2018), specifically in the definition of the wind vector. In Martin et al. (2018), $\mathbf{u}_{10} = \mathbf{u}_{\text{ERW}}$ is defined as the Earth relative wind (the same as for the Numerical Weather Prediction (NWP) product), whereas in this paper $\mathbf{u}_{10} = \mathbf{u}_{\text{OSVW}}$ is the stress-equivalent ocean surface wind at 10 m height, which is equivalent to the difference between the Earth relative wind and the surface current; i.e. $\mathbf{u}_{\text{OSVW}} = \mathbf{u}_{\text{ERW}} - c$. Despite $\mathbf{u}_{\text{ERW}}$ being the wind used in numerical weather prediction, the GMFs for $\sigma^0$ and Doppler shift used are at first order functions of $\mathbf{u}_{\text{OSVW}}$ and not $\mathbf{u}_{\text{ERW}}$, hence this slight modification.

For the results presented in this paper, the KuDop GMF used for Eq. (1) is an adapted version of the Envisat ASAR C-band (5.4 GHz) empirically derived C-DOP GMF (Mouche et al., 2012). The RSV derived from this model has proven to vary little with radar frequency and to be applicable over a wide range with examples in 9.8 GHz X-band (Martin et al., 2016; Elyouncha et al., 2024) up to 35.7 GHz Ka-band (Yurovsky et al., 2019; Polverari et al., 2022). The range of incidence angle for C-DOP is here extended with up to 62° (empirically derived only up to 44°). These assumptions are discussed in Sect. 5. Here C-DOP has been adapted in two ways: firstly, for Ku-band use via a frequency shift in the Doppler calculation and, secondly, extending the model’s applicable upper range of incidence angle from 44 to 62°. For the KuMod GMF the NSCAT-4DS model (Wang et al., 2017) is used.

The parameters $\Delta\sigma^0 = k_p\sigma^0$ (where $k_p$ is the radiometric resolution) and $\Delta\text{RSV}$ represent the uncertainties in the measurements. Here $k_p$ was estimated as a combination of contributions from instrument noise and geophysical noise (Mejia et al., 1999; Portabella and Stoffelen, 2006; Anderson et al., 2012). Instrument noise was estimated via a robust estimator of distribution (normalized interquartile range) of measured $\sigma^0$ for each look direction. Geophysical noise was estimated via comparison of measured $\sigma^0$ with the predicted NRCS using NSCAT-4DS. $\Delta\text{RSV}$ was estimated using open-ocean OSCAR data from the SEASTARex campaign and ground truth data from an acoustic Doppler current profiler (ADCP) mooring at 48.256° N, 5.249° W to estimate the observed WASV and compare it to the predicted WASV from Mouche et al. (2012). The estimated noise parameters used in this study were $k_p = 20\%$ and $\Delta\text{RSV} = 0.2 \text{ m s}^{-1}$.

The cost function $J$ is a unitless function of four unknown variables ($\mathbf{u}_{10}, c$). Minimizing the cost function finds the values of TSCV ($c$) and OSVW ($\mathbf{u}_{10}$) that best reduce the quadratic differences between the measured observables ($\sigma_{\text{obs},i}^0$, $\text{RSV}_{\text{obs},i}$) and the predicted quantities (KuMod, KuDop). As in scatterometry, the minimization returns up to four solutions (Portabella et al., 2002; Martin et al., 2018), leading to a well-known ambiguity problem and the requirement for an ambiguity removal procedure. A usual procedure to remove these ambiguities is to use additional geophysical information, e.g. from wind forecasts (Portabella and Stoffelen, 2004). In the case of the work presented within this paper, a simple ambiguity removal method was implemented, selecting the solution for $J$ closest in $\mathbf{u}_{10}$ space to wind vectors taken from the Meteo France operational Application of Research to Operations at Mesoscale (AROME) atmospheric wind model (Seity et al., 2011).

## 3 Iroise Sea airborne campaign

The SEASTARex airborne campaign was conducted between 17–26 May 2022 over the Iroise Sea, west of Brest, France,
<figure><img src="image_5.png" /><figcaption>**Figure 1.** (a) The OSCAR three-look Ku-band SAR instrument and three-axis stabilization gimbal; the two fore and aft ATI antenna pairs can be seen to the left and right of the boom with the scatterometer middle antenna in the centre. (b) The OSCAR instrument within the radome pod mounted to MetaSensing's Piper Navajo airframe. (c) The study area of the Iroise Sea, France; bathymetry is in metres below mean sea level from EMODnet. (d) The island of Ushant with the positions of the La Jument lighthouse and the over-island flight tracks marked.</figcaption></figure>

from the home airport of Morlaix, Brittany. Macrotidal and relatively shallow, the Iroise Sea is dominated by strong tidal currents and prevailing Atlantic swell interacting with complicated coastline morphology and bathymetry of the area (Muller et al., 2009). The island of Ushant (Ouessant in French) and its coastal waters, at the westerly end of the Molène archipelago in the Iroise Sea, experience some of the fastest tidal flows on the north European shelf, often exceeding $3 \text{ m s}^{-1}$ (Sentchev et al., 2013), and frequent westerly and south-westerly Atlantic storms. Figure 1c shows the location of the study area and the bathymetry (in metres below mean sea level) from the European Marine Observation and Data Network (EMODNet) harmonized digital terrain model.

The airborne scientific campaign consisted of 4 flight days of repeat acquisitions of OSCAR data over three main areas of interest: a site around the island of Ushant with highly dynamic, macrotidal and hydrodynamically heterogeneous currents; a site south of the La Jument lighthouse (Fig. 1d), with temporally and spatially homogenous conditions and deployed in situ buoy measurements (48.256°N, 5.249°W);

and a long, open-ocean flight further south to coincide with measurements from the ASCAT satellite scatterometer. This paper focuses on the results from the OSCAR flights over the dynamic tidal race around Ushant (flight tracks from Fig. 1d) on 17 and 22 May 2022.

## 3.1 OSCAR data acquisitions

The OSCAR flights were scheduled to occur at an altitude 3000 m above sea level for all days; however a low cloud ceiling on 17 May forced the aircraft to acquire data at an altitude of 1950 m. Acquisitions on 22 May were obtained at the intended altitude of 3050 m. All flights were performed at a mean ground speed of $80 \text{ m s}^{-1}$. The flights over Ushant were scheduled to capture peak tidal flows, with the aircraft passing over the island at 09:32 UTC on 17 May and 05:48 UTC on 22 May. OSCAR was configured to generate SAR imagery at 8 m pixel resolution in a 5 km swath, with incidence angles varying between 22–69° for the fore and aft squinted channels and between 16–61° for the middle channel. Level-0 processing (e.g. SAR focusing) and ra-
diometric and interferometric calibration of the OSCAR ac-
quisitions were performed by MetaSensing BV and Radar-
metrics SL. Radiometric calibration was performed via tar-
geted flights over corner reflectors of a known radar cross
section, with these data being recorded on each flight day be-
fore the scientific acquisitions took place. Additional flights
over the corner reflectors were performed at the end of each
flight to check that calibration parameters had not changed.
Additional residual calibration of NRCS was performed us-
ing OSCAR data from the open-ocean flights and computing
the incidence-angle-dependent bias with respect to NSCAT-
4DS. Interferometric calibration was performed using over-
land OSCAR data to assess the recorded Doppler velocities
of static land reflections (which should theoretically be zero).
Level-1 single-look complex Doppler images were spatially
smoothed using a 7-pixel (56 m) windowed mean and down-
sampled to 200 m ground resolution (using the mean value
in each 200 m cell) for input to the simultaneous retrieval.
This resolution was chosen as a trade-off between capturing
sub-mesoscale hydrodynamic features and saving computa-
tion time for the cost function. All campaign data are clas-
sified according to their stage in the processing chain, with
Level-1 (L1) data corresponding to SAR-processed data in
instrument (i.e. local) coordinates, Level-1c (L1c) data cor-
responding to calibrated data on a shared grid and Level-2
(L2) data corresponding to retrieved geophysical parameters
in a global coordinate system.

## 3.2 X-band marine radar

As part of the SEASTARex project, a marine radar was in-
stalled on the La Jument lighthouse overlooking the tidal race
to the south-west of Ushant, which is an ideal location for
observing the extreme tidal dynamics of the area (Filipot et
al., 2019). The coherent-on-receive X-band (9.3 GHz) radar
(Horstmann et al., 2021) was mounted on the lighthouse at
a height of 48 m above mean sea level and covered a ra-
dius between 52 and 3150 m, corresponding to grazing an-
gles between 40 and 0.78° (equal to incidence angles of 50
and 89.22° respectively). During the campaign the radar was
operated in its rotational mode, acquiring radar backscatter
intensity and radial Doppler velocity maps with a repetition
rate of approximately 0.5 Hz. All radar data were collected at
a pulse length of 0.5 ns and with a pulse repetition frequency
of 2 kHz at VV polarization with a 2.3 m slotted waveguide
antenna, resulting in an azimuthal resolution of 1.2°. Radar
data were acquired at 20 MHz, resulting in a range reso-
lution of 7.5 m. The radar was operated on 17 May 2022
between 06:41–12:00 UTC, collecting six 10 min video se-
quences every hour, which were used to compute surface
current fields. One of these sequences was coincident with
an OSCAR flight overpass of the radar at 09:38 UTC.

The surface current maps are derived from the wave sig-
nal within a 10 min marine radar image sequence, which is
then transformed from the spatial–temporal domain to the

wave number-frequency domain by a 3D fast Fourier trans-
form (FFT). Within this 3D radar image spectrum, the wave
signal is located on an inverted cone – the so-called disper-
sion shell defined by the linear dispersion relation of surface
waves. In the case of deep water with respect to the wave-
length, this dispersion shell is solely dependent on the sur-
face current (Senet at al., 2001; Huang et al., 2016). The 2D
fields of current vectors are determined by an algorithm that
searches for the current that maximizes the energy on the
dispersion shell using a brute-force optimization algorithm
(Streßer et al., 2017) and considering wave lengths between
15–125 m and wave periods between 4–20 s. The spatial win-
dow used for one individual current measurement spans an
area of 500 × 500 m, and the individual windows overlap
by 200 m, which corresponds to the pixel ground resolution
(posting) in the resulting 2D current map.

## 3.3 Numerical models

Data from two numerical models were used: the MARS2D
depth-averaged hydrodynamic model (Lazure and Dumas,
2008), run operationally by the Laboratory for Ocean
Physics and Satellite remote sensing (LOPS), and the Meteo
France operational AROME wind model (Seity et al., 2011).
MARS2D depth-averaged currents were provided at 15 min
intervals at a ground resolution of 250 m. Hourly AROME
wind forecasts were provided for an altitude of 10 m (u₁₀) at
a ground resolution of 0.025° (approx. 2 km).

## 3.4 Spaceborne SAR imagery

NovaSAR-1 is a UK-funded technology demonstration satel-
lite owned and operated by Surrey Satellite Technology Ltd.
(SSTL) that delivers high-resolution S-band SAR imagery.
NovaSAR-1 was commissioned to acquire images over the
study site, resulting in an overpass at 10:30 UTC on 17 May
2022, coinciding with the OSCAR flight over Ushant with
only an hour difference. The NovaSAR-1 data were commis-
sioned and accessed thanks to the support of Martin Cohen
at Airbus Defence and Space Ltd.

# 4 Results

## 4.1 OSCAR results for ebb tide on 17 May 2022

Figure 2 shows the retrieved TSCV (Fig. 2a) and OSVW
(Fig. 2b) at 200 m ground resolution from a single OS-
CAR acquisition flying north (looking left) over the island
of Ushant during ebb tide at 09:38 UTC on 17 May 2022.
Sea state conditions at the time of acquisition consisted of
a moderate swell from the south-west with a significant
wave height of 2.5 m and a peak period of 11 s. Winds
were approximately 7 m s⁻¹ from the south-east. The vec-
tors in Fig. 2a and b are plotted with a sub-sampled post-
ing of 400 m for clarity. The trapezoidal shape of the OS-
CAR swath is due to the squinted look direction (±45° in azimuth) of the fore and aft ATI channels and their orthogonal combination to compute the L2 current and wind vectors. For comparison, Fig. 2c shows depth-averaged current vectors at 250 m ground resolution (vector posting at 500 m) from the MARS2D ocean circulation model at 09:30 UTC, and Fig. 2d shows predicted stress-equivalent wind vectors at 10 m above the sea surface ($u_{10}$) at 2 km ground resolution from the hourly Meteo France AROME operational wind forecast model at 09:00 UTC. The OSCAR L2 swath can be seen outlined in black, overlain on both model outputs.

Overall, very good visual agreement can be seen between OSCAR and the two models, in both the magnitude and the direction of both TSCV and OSVW. OSCAR features additional complexity in the observed current field compared with the model, with sharper gradients and a separate submesoscale filament close to land in the northern jet north of the island. To the south of the domain, the counter-rotating eddy predicted by the model data (~48.38°N) is also captured in the OSCAR TSCV, albeit with subtle spatial differences and lower magnitude compared to the observations. Importantly, OSCAR reveals tidal currents north of the island that are both more intense and further from land (due to leeward sheltering of the tidal flow). For OSVW, good similarity can be seen between the retrieved OSVW from OSCAR (Fig. 2b) and AROME model data (Fig. 2d), with finescale variations in the near-surface winds present in the OSCAR data, especially around the western tips of the island, which could originate from local orographic perturbation. The large-scale north-south gradient in wind speed present in the numerical model is also observed in the OSCAR data, albeit with much finer-scale detail due to the extreme difference in resolution.

Further evidence of the exceptional accuracy of OSCAR in correctly measuring complex hydrodynamic structures around Ushant is shown in Fig. 3 via comparison and validation against surface current vectors derived from the X-band radar on the La Jument lighthouse. Figure 3a shows surface currents derived from 10 min of sequential X-band intensity imagery acquired over a radius of 3.1 km around the lighthouse at 09:30 UTC on 17 May, shown with a pixel ground resolution of 200 m, a vector posting of 400 m and the coincident OSCAR L2 swath overlain. Figure 3b shows co-located and superimposed vectors from both sensors plotted at 400 m posting. The two instruments show very good agreement, with a root mean square difference (RMSD) between the data of 0.18 m s⁻¹ for current velocity and 5.27° for current direction (Fig. 3c). Both the OSCAR system and the X-band radar measure the southern tidal jet further west than it appears in the model (Fig. 2c) and a northward shift in the counter-rotating flow. This is reflected by the RMSD between MARS2D and the measurements from these two remote sensing systems, which report an RMSD of 0.61 m s⁻¹ and 6.5° between OSCAR and the model and 0.47 m s⁻¹ and

16.6° between the X-band and the model for current velocity and current direction respectively. The high level of agreement between the two independent remote sensing measurements demonstrates the value of these systems in revealing and resolving inaccuracies that can occur in models of high-energy coastal environments.

A more detailed examination of Fig. 3b indicates that the best agreement can be seen at the western edge of the OSCAR swath, which corresponds to higher incidence angles. The comparison around the position of the high tidal flow (at lower incidence angles) is better for TSCV magnitude than it is for direction, with a clear shift in measured current direction between OSCAR and the X-band radar data. The validation of TSCV components from OSCAR and the X-band radar (Fig. 3c) reveals good agreement in U (east component; $R^2$ = 0.7) and excellent agreement in V (north component; $R^2$ = 0.96) across a wide range of flow conditions. The combined RMSD between the two sensors for TSCV at 200 m ground resolution is found to be 0.19 m s⁻¹. The departure between the two sensors in negative component velocities is related to the different OSCAR TSCV directions observed in the southern tidal jet.

Further validation of the ability of OSCAR to accurately measure fine-scale hydrodynamic structures in such a dynamic coastal environment is seen in Fig. 4. Here OSCAR TSCV vectors from 09:38 UTC on 17 May are overlain over a geo-registered SAR image from the NovaSAR-1 satellite around Ushant. The satellite image was acquired at 10:30 UTC, less than an hour after the OSCAR track, but still within the same ebb-tidal flow regime. The NovaSAR-1 S-band $\sigma^0$ image displays variations in the normalized radar cross section ($\sigma^0$) linked to the modulation of the ocean surface roughness by strong horizontal shear in the flow, producing clear intensity gradients that are coincident with these hydrodynamic structures. Excellent spatial correlation can be observed between the OSCAR TSCV vectors and these gradients in the SAR image, with good correspondence in the positions of both northern and southern tidal jets around Ushant in the spaceborne SAR $\sigma^0$ and OSCAR data. The NovaSAR-1 image also displays lighter gradients that coincide with both the counter-rotating flow to the south-west of Ushant and the small flow filament that are observed in the OSCAR data but are not present in the MARS2D output. The comparison with NovaSAR-1 confirms the validity of OSCAR TSCV placing the intense horizontal shear of the ebb-tidal jet north of the island (~48.48°N) approximately 1 km northward of the position predicted with MARS2D (Fig. 2c). The same is true for the position of the southern tidal jet towards the lower extremity of the OSCAR swath. This is an important result, as the exact position of such strong gradients in tidal flows such as these can be challenging to accurately predict using numerical models.
<figure><img src="image_6.png" /><figcaption>**Figure 2.** Simultaneously retrieved L2 data products from a single OSCAR acquisition over Ushant at 09:38 UTC on 17 May 2022 during an ebb tide: (a) OSCAR total surface current vector, (b) OSCAR ocean surface vector wind, (c) depth-averaged current velocity from the MARS2D model at 09:30 UTC and (d) wind speed ($u_{10}$) from the AROME operational wind model at 09:00 UTC. The outlines of the OSCAR L2 swaths are shown in black on the two model outputs.</figcaption></figure>

## 4.2 OSCAR results for flood tide on 22 May 2022

Figure 5 shows TSCV and OSVW results from an OSCAR acquisition flying north (looking left) and covering part of the flood-tidal regime over Ushant on 22 May 2022 at 05:48 UTC. Sea state conditions at the time of acquisition consisted of a light swell from the south-west with a significant wave height of 1 m and a peak period of 9 s. The wind was generally from the north-east and around 5–6 m s⁻¹. Figure 5a and b show OSCAR-retrieved TSCV and OSVW respectively. Figure 5c shows MARS2D-simulated currents around Ushant at 05:45 UTC, and Fig. 5d shows forecast $u_{10}$ wind vectors at 06:00 UTC from the AROME model. Once again, OSCAR presents very good overall agreement with the two models, in both the magnitude and the direction of both TSCV and OSVW. Moreover, the results confirm the

ability of OSCAR to resolve not only high tidal currents but also current fields in quiescent areas, for example in the shallow coastal areas around western Ushant that are sheltered from the northward flood tide which is present in both the predicted and the observed TSCV fields. Similar to the ebb-tide case on 17 May (previous section), the OSCAR TSCV field shows a more northerly extent of the accelerated tidal flow within the imaged swath than in the predicted data, highlighting the stark differences that can be obtained between observations and numerical models in such dynamic coastal areas. The OSCAR OSVW (Fig. 5b) captures the same general north-south trend and variability seen in the forecast AROME wind data (Fig. 5d), with the reduction in wind speeds to the lee of the island clearly visible in both the OSCAR results and the predictions.
<figure><img src="image_7.png" /><figcaption>**Figure 3.** (a) X-band-derived surface current vectors from the La Jument lighthouse from 10 min of radar data at 09:30 UTC on 17 May 2022, (b) co-located OSCAR total surface current vectors at 09:38 UTC and X-band current vectors, and (c) U (east) and V (north) current component direct comparison between the two sensors.</figcaption></figure>

On 22 May, three OSCAR acquisitions were made over the same area to the west of Ushant within a period of 15 min, with two passes in a southerly direction and one travelling north. The objective was to demonstrate the consistency of the OSCAR L2 products between successive passes and in different orientations. Figure 6 shows the comparison of TSCV for the three tracks (shown as coloured vectors) co-located on the same grid and overlain on MARS2D-predicted currents (black vectors). The swath outlines are in orange for the northward track at 05:48 UTC and in blue for the southward tracks at 05:39 and 05:54 UTC (for clarity only one swath is outlined for the two southward tracks). The agreement between successive passes is excellent, with different TSCV vectors nearly always perfectly superposed and quasi-undistinguishable at most grid points. The TSCV median sample standard deviation between the three flights (with each TSCV vector treated as an independent measurement, i.e. *n* = 3) is 0.10 m s⁻¹ for current velocity and 1.18° for current direction. It is worth noting that the radar lines of sight relative to the wind were not particularly favourable in this instance. As the wind direction was from the north-east (Fig. 5d), with the antenna squint angles at 45° and the flight directions either north or south, the southward

flights at 05:39 and 05:54 UTC (blue and green) correspond to one squinted beam looking crosswind and the other looking upwind, whereas the northward flight at 05:48 UTC (orange) had one squinted beam looking crosswind and the other downwind. This is important, as the performance of the three-look concept is known to be impacted when one of the azimuth lines of sight is aligned with the wind (Stoffelen and Portabella, 2006; Martin et al., 2018). However, the consistency between the TSCV obtained with upwind (blue and green) and downwind (orange) lines of sight is seen to be very good at the points where the swaths from the northward and southward flights overlap.

## 5 Discussion, conclusions and future work

This paper presents the results of the first scientific airborne campaign of the OSCAR high-squint three-look SAR-ATI instrument and represents the first demonstration of simultaneous retrieval of TSCV and OSVW in a single pass using a SeaSTAR-type observing concept in a highly dynamic coastal marine environment. The results are extremely promising, with OSCAR and the SeaSTAR inver-
<figure><img src="image_8.png" /><figcaption>**Figure 4.** OSCAR total surface current vectors from 09:38 UTC on 17 May 2022, overlain on a NovaSAR-1 S-band σ⁰ image acquired at 10:30 UTC over Ushant, and the surrounding waters during a period of ebb-tidal flow. The NovaSAR-1 image is courtesy of SSTL and Airbus.</figcaption></figure>

sion algorithm reporting a high degree of accuracy and self-consistency over a large swath in highly inhomogeneous coastal conditions. The OSCAR instrument demonstrates the capability to provide quantitative information about fine-scale dynamics that are not correctly represented with numerical models and are difficult to observe with ground-based sensors alone. In situ observations of hydrodynamics in high-energy macrotidal environments are costly and challenging (Neil and Hashemi, 2018) and provide only limited spatial detail needed to understand frontal dynamics in areas of strong spatial gradients. High-frequency (HF) radars could provide synoptic information about ocean currents, but their coverage is patchy and limited to a few coastal regions in some industrialized countries. Whilst HF radar data were recorded during the campaign, coverage around the island of Ushant, the focus of this paper, is unfortunately lacking, especially in the area to the west of the island that sits in

the radio shadow of the HF ground stations on the French mainland. However, HF coverage of the open-ocean homogeneous site sampled in the campaign was good, and its analysis and comparison with OSCAR will be a subject of future work. The quality of the OSCAR data obtained from these four flights spanning 2 d is notably impressive. Comparisons of the three flights for 22 May, acquired within a short time over roughly the same ground track but with radically different radar viewing directions with respect to the wind, give confidence in the OSCAR and SeaSTAR current and wind vector retrieval. Comparisons with high-resolution data from X-band radar, NovaSAR-1 and models demonstrate that the OSCAR airborne instrument and the SeaSTAR mission concept clearly bring great additional capability to the research community and present great promise for increasing our understanding of small-scale processes in coastal and shelf seas and marginal ice zones.

It is noteworthy that the retrieval returns excellent results within this exceptionally heterogeneous coastal environment whilst using the open-ocean wind-dependent WASV correction by Mouche et al. (2012). Importantly, no specific adjustments were made to account for additional sea state dependency of the WASV on factors like fetch, wave age, wave breaking at fronts or wave-current interactions over shallow bathymetry. Our results align with prior suggestions from Martin and Gommenginger (2017) that indicate that the effects on SAR-ATI of strong current gradients and wave age may hold secondary importance in these highly dynamic areas. There is a common perception that the accuracy of available open-ocean GMFs may be insufficient to predict the WASV in areas of strong surface current gradients and may require mixed-polarization sensor capability to correctly address additional Doppler effects at fronts (Kudryavtsev et al., 2014; Martin et al., 2018). A surprising result of our work is the consistency of TSCV retrieved using the C-DOP GMF of Mouche et al. (2012) and its predicted WASV in such a dynamic coastal tidal environment. The C-DOP GMF captures the average behaviour of the WASV based on a large globally distributed dataset of Envisat ASAR Doppler measurements over the open ocean (i.e. hydrodynamically homogeneous, microtidal environments) and as such does not contain the WASV response to shallow-water dynamics, wave-current interactions, strong tidal shear, etc. Additionally, the C-DOP GMF has been adapted for this work with respect to frequency (from C-band to Ku-band use) and extended in its upper applicable range of incidence angle from 44 to 62°. Initial comparison tests applying the GMF of Yurovsky et al. (2019) showed minor differences with our extended version of C-DOP at these high incidence angles. A more in-depth analysis of the WASV and the effect of these different GMFs will be conducted in a future study. The quality of the OSCAR-retrieved TSCV in the complicated Iroise Sea environment suggests that the assumed sensitivities of the WASV to hydrodynamic processes are perhaps not as important as previously assumed. Future work will consider other GMFs
<figure><img src="image_9.png" /><figcaption>**Figure 5.** Simultaneously retrieved currents and winds from a single OSCAR acquisition over Ushant at 05:48 UTC on 22 May 2022 during a flood tide: (a) OSCAR total surface current vector, (b) OSCAR ocean surface vector wind, (c) depth-averaged current velocity from the MARS2D model at 05:45 UTC and (d) wind speed ($u_{10}$) from the AROME operational wind model at 06:00 UTC. The outlines of the OSCAR L2 swaths are shown in black on the two model outputs.</figcaption></figure>

to compute the WASV, including those that depend on ocean wave parameters, such as the ones developed by Yurovsky et al. (2019). TSCV retrieval performance is especially sensitive to the choice of WASV at low incidence angles (Martin et al., 2018), so a different GMF may further improve the OSCAR results in the near-range parts of the swath. The full wind and wave parametrized model for Yurovsky et al. (2019) has not been tested with OSCAR data at this time but will be the basis of a larger validation study in the future.

All results of this study were obtained using the simultaneous retrieval of current and wind vectors, which is the chosen baseline algorithm for OSCAR and SeaSTAR Level-2 inversion. However, other inversion approaches exist; for example, one can consider a sequential retrieval approach where ancillary wind information (for example from AROME) serves to correct the WASV directly in each line of sight, be-

fore recovering the current vector (as in Martin et al., 2017). In this instance, computation time is almost instantaneous compared to a typical ~30 min computation on a 20-core machine to apply simultaneous inversion to a single track at 200 m resolution. However, our tests (not shown) indicate that the quality of the TSCV is degraded when using the sequential approach. The degradation is visible in the poorer performance against X-band marine radar currents but also in the mislocalization of OSCAR-retrieved current gradients with the visible roughness gradients in NovaSAR-1 imagery. This is particularly noticeable at low incidence angles and may be linked to errors inherent to the WASV GMF mentioned previously.

TSCV is defined as effective mass transport, being composed of any movement of the sea surface that is not accounted for by the WASV, i.e. from the wind (in the case
<figure><img src="image_10.png" /><figcaption>**Figure 6.** Total surface current vector comparison between three successive OSCAR tracks on 22 May, co-located onto a common grid and overlain on MARS2D current vectors from 05:30 UTC (black arrows). OSCAR swath outlines are depicted in their corresponding vector colour (blue, 05:39 UTC; orange, 05:48 UTC; green, 05:54 UTC). The tracks from 05:39 and 05:54 UTC share the same swath outline.</figcaption></figure>

of the GMF of Mouche et al., 2012, or wind and waves in the case of Yurovsky et al., 2019). In this sense, TSCV contains contributions from a wide variety of physical processes that may be isolated depending on the analysis or secondary information used. In the case of the work presented in this paper, the TSCV is completely dominated by the contribution from the tide, being orders of magnitude greater than e.g. Stokes drift or Ekman currents. Being a surface imaging system, OSCAR and the SeaSTAR methodology also capture only the ocean motion that is manifest at the surface. The currents sensed by the X-band radar used for comparison are barotropic, as are the currents produced by the MARS2D numerical model (which is depth averaged). Unravelling baroclinic contributions to the Doppler signal sensed at the surface will require a different experimental design in an area that is not dominated by tidal currents (e.g. the open ocean).

The comparison between OSCAR and ground-based X-band radar data is excellent, especially when considering the differences in frequency and imaging mechanisms between the two systems. Some consideration should be given to these differences and how they may account for a proportion of the scatter seen in Fig. 3c. For the X-band radar (e.g. using the methods of Streßer et al., 2017), the resulting

derived currents represent a weighted mean over the upper ocean, where the greatest weight is assigned to the surface, and the effective depth of the current relates to the maximum ocean wavelengths imaged by the radar and considered in the current fit. For OSCAR, the TSCV derives from Doppler shifts in backscatter at moderate, non-grazing angles (incidence angle as steep as 22°) where there are combined effects from small-scale surface roughness and longer ocean wavelengths, particularly wind waves around 10 m (Chapron et al., 2005). The agreement between OSCAR and the X-band-derived surface currents in Fig. 3b is seen to be better in one area (the northerly vectors in the centre left of the swath) than others (the bottom right of the swath). Aside from differences in imaging and retrieval methods between the two data, another source of the difference may be the combination of the antenna look direction relative to the surface current and errors in the WASV correction at low incidence angles in the near range of the swath. At these low incidence angles the WASV is more sensitive to wind (Mouche et al., 2012), which in the simultaneous retrieval also has its own sources of error, e.g. from the calibration of NRCS and its effect on the output of the NSCAT-4DS GMF. Combined with the fact that the current flows in broadly the same direction as the one of the look azimuths (for the aft beam) and perpendicular to another (the fore beam), errors in the WASV could be magnified more than in other areas of the swath. This is a complicated interplay of calibration factors, SAR interferometry and GMFs that warrants further investigation. Further investigations are also needed to better understand differences in TSCV measured by different sensors, including HF radars and other ocean current sensors like ADCP, and how to accommodate these differences when validating new sensors like OSCAR.

This paper presents only a preliminary validation of the OSCAR data to demonstrate the innovative and promising spatial mapping capabilities of this new system in the extreme macrotidal coastal environment close to Ushant. However, more comprehensive validation of the OSCAR L2 results is necessary to quantitatively compare the OSCAR data against established ocean measurements from ADCP and HF radar also collected during the SEASTARex project. For this, the focus of the validation must shift south, to a second instrumented site set up during SEASTAREx to validate OSCAR in a geophysically homogenous site to the south of Ushant (centred on 48.256°N, 5.249°W). Those larger datasets and comparisons will be the object of a separate publication. More recently, OSCAR flew in its second scientific campaign in May 2023 in the north-west Mediterranean Sea (north of the island of Menorca, Spain), which is known for its strong sub-mesoscale ocean dynamics. The OSCAR-Med airborne flights were timed to coincide with overpasses of the Surface Water and Ocean Topography (SWOT) mission during its 1d fast-repeat Cal/Val phase and in situ measurements from the BioSWOT-Med oceanographic ship campaign (Doglioli and Gregori, 2023). The results of the
OSCAR-Med campaign will be the subject of a future study. Finally, given the scientific objectives of the SeaSTAR mission concept to measure TSCV and OSVW in marginal ice zones, the team is currently exploring opportunities to fly a new OSCAR scientific campaign close to the sea ice edge.

*Code availability.* The repository for the SeaSTAR processing software used for the work in this study is available at https://doi.org/10.5281/zenodo.10026593 (Martin et al., 2023).

*Data availability.* OSCAR data from the Iroise Sea campaign are available from ESA's Earth Online database (https://doi.org/10.57780/esa-633ce94, ESA, 2024) and can also be provided by the corresponding author upon request.

*Author contributions.* ACHM, CG, JH, LM, JMM, MP and AM planned the campaign; TC and PMI contributed to the campaign planning and realization; AM managed the recording of the airborne data; DLM, ACHM, KM and RC developed the processing software; KM, JMM and MP performed the OSCAR Level-1 data processing and calibration; ACHM and DLM performed the OSCAR Level-2 data processing; RC, JH and LM recorded and processed the ground truth data; DLM wrote the paper; RC and JH provided text input for the paper; and DLM, ACHM and CG reviewed and edited the paper.

*Competing interests.* The contact author has declared that none of the authors has any competing interests.

*Disclaimer.* Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. While Copernicus Publications makes every effort to include appropriate place names, the final responsibility lies with the authors.

*Acknowledgements.* The authors would like to thank Christian Trampuz (formerly at MetaSensing BV), Hugo Kerhuel and Thomas Jamne for the preparation and execution of the OSCAR Iroise Sea campaign in May 2022. The authors would also like to thank Nicholas Gebert from ESA, who was the initial technical officer for the OSCAR system development; Rui Duarte and the team at France Énergies Marines (FEM), who provided essential access to the La Jument lighthouse and support during the installation and operation of the X-band radar; Thiago Luiz at MetaSensing BV for his work in the SAR processing of OSCAR data; and Guissepe Greco at ICM-CSIC, Spain, and Wenming Lin at NUIST, China, for their help with the radiometric calibration of OSCAR data.

Many thanks to Clive Neil at NOC for his essential assistance in setting up the SeaSTAR project software repository.

*Financial support.* This research has been supported by the European Space Agency (grant nos. 400017623/22/NL/IA and 4000116401/16/NL/BJ).

*Review statement.* This paper was edited by Ismael Hernández-Carrasco and reviewed by Alejandro Orfila and one anonymous referee.

## References

Anderson, C., Bonekamp, H., Duff, C., Figa-Saldana, J., and Wilson, J. J. W.: Analysis of ASCAT Ocean Backscatter Measurement Noise, IEEE T. Geosci. Remote, 50, 2449–2457, https://doi.org/10.1109/TGRS.2012.2190739, 2012.

Bricheno, L. M., Soret, A., Wolf, J., Jorba, O., and Baldasano, J. M.: Effect of high-resolution meteorological forcing on nearshore wave and current model performance, J. Atmos. Ocean. Tech., 30, 1021–1037, https://doi.org/10.1175/jtech-d-12-00087.1, 2013.

Buck, C.: An extension to the wide swath ocean altimeter concept, in: Proceedings of the 2005 IEEE International Geoscience and Remote Sensing Symposium, 29 July 2005 Seoul, Korea(South), IGARSS '05, vol. 8, IEEE, N. Y., 5436–5439, https://doi.org/10.1109/IGARSS.2005.1525970, 2005.

Chapron, B., Collard, F., and Ardhuin, F.: Direct measurements of ocean surface velocity from space: Interpretation and validation, J. Geophys. Res., 110, https://doi.org:/10.1029/2004JC002809, 2005.

de Kloe, J., Stoffelen, A., and Verhoef, A.: Improved Use of Scatterometer Measurements by Using Stress-Equivalent Reference Winds, IEEE J-STARS, 5, 2340–2347, https://doi.org/10.1109/JSTARS.2017.2685242, 2017.

Doglioli, A. and Gregori, G.: BioSWOT-Med cruise, RV L'Atalante, Ifremer, https://doi.org/10.17600/18002392, 2023.

Elyouncha, A., Broström, G., and Johnsen, H.: Synergistic utilization of spaceborne SAR observations for monitoring the Baltic Sea flow through the Danish straits, ESS Open Archive, https://doi.org/10.22541/essoar.171466079.98038905/v1, 2024.

European Space Agency (ESA): 2022 SEASTARex, [data set], https://doi.org/10.57780/esa-633ce94, 2024.

Filipot, J.-F., Guimaraes, P., Leckler, F., Hortsmann, J., Carrasco, R., Leroy, E., Fady, N., Accensi, M., Prevosto, M., Duarte, R., Roeber, V., Benetazzo, A., Raoult, C., Franzetti, M., Varing, A., and Le Dantec, N.: La Jument lighthouse: a real-scale laboratory for the study of giant waves and their loading on marine structures, Philos. T. Roy. Soc. A, 377, 20190008, https://doi.org/10.1098/rsta.2019.0008, 2019.

Goldstein, R. and Zebker, H.: Interferometric radar measurement of ocean surface currents, Nature, 328, 707–709, https://doi.org/10.1038/328707a0, 1987.

Gommenginger, C., Chapron, B., Hogg, A., Buckingham, C., Fox-Kemper, B., Eriksson, L., Soulat, F., Ubelmann, C., Ocampo-Torres, F., Nardelli, B. B., Griffin, D., Lopez-Dekker, P., Knudsen, P., Andersen, O., Stenseng, L., Stapleton, N., Perrie, W., Violante-Carvalho, N., Schulz-Stellenfleth, J., Woolf, D., Isern-Fontanet, J., Ardhuin, F., Klein, P., Mouche, A., Pascual, A., Capet, X., Hauser, D., Stoffelen, A., Morrow, R., Aouf, L.,
Breivik, Ø., Fu, L.-L., Johannessen, J. A., Aksenov, Y., Bricheno, L., Hirschi, J., Martin, A. C. H., Martin, A. P., Nurser, G., Polton, J., Wolf, J., Johnsen, H., Soloviev, A., Jacobs, G. A., Collard, F., Groom, S., Kudryavtsev, V., Wilkin, J., Navarro, V., Babanin, A., Martin, M., Siddorn, J., Saulter, A., Rippeth, T., Emery, B., Maximenko, N., Romeiser, R., Graber, H., Az-carate, A. A., Hughes, C. W., Vandemark, D., da Silva, J., van Leeuwen P. J., Naveira-Garabato, A., Gemmrich, J., Ma-hadevan, A., Marquez, J., Munro, Y., Doody, S., and Bur-bidge, G.: SEASTAR: A Mission to Study Ocean Subme-soscale Dynamics and Small-Scale Atmosphere-Ocean Pro-cesses in Coastal, Shelf and Polar Seas, Front. Mar. Sci., 6, 457, https://doi.org/10.3389/fmars.2019.00457, 2019.

Hauser, D., Abdalla, S., Ardhuin, F., Bidlot, J.-F., Bourassa, M., Cotton, D., Gommenginger, C., Evers-King, H., Johnsen, H., Knaff, J., Lavender, S., Mouche, A., Reul, N., Sampson, C., Steele, E. C. C., and Stoffelen, A.: Satellite Remote Sensing of Surface Winds, Waves, and Currents: Where are we Now?, Surv. Geophys., 44, 1357–1446, https://doi.org/10.1007/s10712-023-09771-2, 2023.

Holt, J., Hyder, P., Ashworth, M., Harle, J., Hewitt, H. T., Liu, H., New, A. L., Pickles, S., Porter, A., Popova, E., Allen, J. I., Sid-dorn, J., and Wood, R.: Prospects for improving the representa-tion of coastal and shelf seas in global ocean models, Geosci. Model Dev., 10, 499–523, https://doi.org/10.5194/gmd-10-499-2017, 2017.

Horstmann, J., Bódewadt, J., Carrasco, R., Cysewski, M., See-mann, J., and Stresser, M.: A coherent on receive x-band marine radar for ocean observations, Sensors, 21, 7828, https://doi.org/10.3390/s21237828, 2021.

Huang, W., Carrasco, R., Chengxi, S., Gill, E. W., and Horstmann, J.: Surface current measurements using X-band marine radar with vertical polarization, IEEE T. Geosci. Remote, 54, 2988–2996, https://doi.org/10.1109/TGRS.2015.2509781, 2016.

Kudryavtsev, V., Kozlov, I., Chapron, B., and Johann-essen, J.: Quad-polarization SAR features of ocean currents, J. Geophys. Res.-Oceans, 119, 6046–6065, https://doi.org/10.1002/2014JC010173, 2014.

Lazure, P. and Dumas, F.: An external–internal mode cou-pling for a 3D hydrodynamical model for applications at regional scale (MARS), Adv. Water Resour., 31, 233–250, https://doi.org/10.1016/j.adwwatres.2007.06.010, 2008.

Lévy, M., Ferrari, R., Franks, P. J., Martin, A. P., and Rivière, P.: Bringing physics to life at the submesoscale, Geophys. Res. Lett., 39, L14602, https://doi.org/10.1029/2012GL052756, 2012.

Martin, A. and Richards, K.: Mechanisms for vertical nutri-ent transport within a North Atlantic mesoscale eddy, Deep-Sea Res. Pt. II, 48, 757–773, https://doi.org/10.1016/S0967-0645(00)00096-5, 2001.

Martin, A. C. H. and Gommenginger, C.: Towards wide-swath high-resolution mapping of total ocean surface current vectors from space: Airborne proof-of-concept and validation, Remote Sens. Environ., 197, 58–71, https://doi.org/10.1016/j.rse.2017.05.020, 2017.

Martin, A. C. H., Gommenginger, C., Marquez, J., Doody, S., Navarro, V., and Buck, C.: Wind-Wave induced velocity in ATI SAR Ocean Surface Currents: First experimental evidence from an airborne campaign, J. Geophys. Res.-Oceans, 121, 1640–1653, https://doi.org/10.1002/2015JC011459, 2016.

Martin, A. C. H., Gommenginger, C., and Quilfen, Y.: Simulta-neous ocean surface current and wind vectors retrieval with squinted SAR interferometry: Geophysical inversion and per-formance assessment, Remote Sens. Environ., 216, 798–808, https://doi.org/10.1016/j.rse.2018.06.013, 2018.

Martin, A. C. H., McCann, D. L., Gommenginger, C., Macedo, K. A. C., and Le Merle, E.: Seastar project, Zenodo [code], https://doi.org/10.5281/zenodo.10026593, 2023.

Maskell, J.: Modelling storm surges in the Irish and Celtic seas us-ing a finite element model (TELEMAC), PhD dissertation, Uni-versity of Liverpool, 280 pp., 2012.

Mass, C. F., Owens, D., Westrick, K., and Colle, B. A.: Does increasing horizontal resolution pro-duce more skillful forecasts?, B. Am. Meteorol. Soc., 83, 406–430, https://doi.org/10.1175/1520-0477(2002)083<0407:DIHRPM>2.3.CO;2, 2002.

McWilliams, J. C.: Submesoscale currents in the ocean, Proc. Math. Phys. Eng. Sci., 472, 20160117, https://doi.org/10.1098/rspa.2016.0117, 2016.

Mejia, C., Badran, F., Bentamy, A., Crepon, M., Thiria, S., and Tran, N.: Determination of the geophysical model function of NSCAT and its corresponding variance by the use of neu-ral networks, J. Geophys. Res.-Oceans, 104, 11539–11556, https://doi.org/10.1029/1998JC900118, 1999.

Mouche, A. C. H., Collard, F., Chapron, B., Dagestad, K-F., Guitton, G., Johannessen, J., Kerabol, V., and Hansen, M. W.: On the use of Doppler shift for sea surface wind re-trieval from SAR, IEEE T. Geosci. Remote, 50, 2901–2909, https://doi.org/10.1109/TGRS.2011.2174998, 2012.

Muller, H., Blanke, B., Dumas, F., Lekien, F., and Ma-riette, V.: Estimating the Lagrangian residual circula-tion in the Iroise Sea, J. Marine Syst., 78, S17–S36, https://doi.org/10.1016/j.jmarsys.2009.01.008, 2009.

Neil, S. P. and Hashemi, M. R.: Chapter 8 – Ocean Mod-elling for Resource Characterization, in: Fundamentals of Ocean Renewable Energy, Academic Press, 193–235, https://doi.org/10.1016/B978-0-12-810448-4.00008-2, 2018.

Polverari, F., Wineteer, A., Rodríguez, E., Perkovic-Martin, D., Siqueira, P., Farrar, J. T., Adam, M., Closa Tarrés, M., and Edson, J. B.: A Ka-Band Wind Geophysical Model Function Using Doppler Scatterometer Measurements from the Air-Sea Interaction Tower Experiment, Remote Sens.-Basel, 14, 2067, https://doi.org/10.3390/rs14092067, 2022.

Portabella, M. and Stoffelen, A.: A probabilistic approach for sea-winds data assimilation, Q. J. Roy. Meteor. Soc., 130, 127–152, https://doi.org/10.1256/qj.02.205, 2004.

Portabella, M. and Stoffelen, A.: Scatterometer Backscatter Uncer-tainty Due to Wind Variability, IEEE T. Geosci. Remote, 44, 3356–3362, https://doi.org/10.1109/TGRS.2006.877952, 2006.

Portabella, M., Stoffelen, J., and Johannessen, A.: Toward an optimal inversion method for synthetic aperture radar wind retrieval, J. Geophys. Res.-Oceans, 107, 3086, https://doi.org/10.1029/2001JC000925, 2002.

Ruiz, S. A. G., Barriga, J. E. C., and Martínez, J. A.: Assessment and validation of wind power potential at convection-permitting resolution for the Caribbean region of Colombia, Energy, 244, 123127, https://doi.org/10.1016/j.energy.2022.123127, 2022.
Samelson, R. M.: Challenges and opportunities in coastal prediction, Eos T. Am. Geophys. Un., 100, https://doi.org/10.1029/2019EO113841, 2019.

Sasaki, H., Klein, P., and Sasai, Y.: Impact of oceanic-scale interactions on the seasonal modulation of ocean dynamics by the atmosphere, Nat. Commun., 5, 5636, https://doi.org/10.1038/ncomms6636, 2014.

Seity, Y., Brousseau, P., Malardel, S., Hello, G., Bénard, P., Bouttier, F., Lac, C., and Masson, V.: The AROME-France Convective-Scale Operational Model, Mon. Weather Rev., 139, 976–991, https://doi.org/10.1175/2010MWR3425.1, 2011.

Senet, C., Seemann, J., and Ziemer, F.: The near-surface current velocity determined from image sequences of the sea surface, IEEE T. Geosci. Remote, 39, 492–505, https://doi.org/10.1109/36.911108, 2001.

Sentchev, A., Forget, P., Barbin, Y., and Yaremchuk, M.: Surface circulation in the Iroise Sea (W. Brittany) from high resolution HF radar mapping, J. Marine Syst., 109–110, S153–S168, https://doi.org/10.1016/j.jmarsys.2011.11.024, 2013.

Schulz-Stellenfleth, J. and Stanev, E. V.: Analysis of the upscaling problem – A case study for the barotropic dynamics in the North Sea and the German Bight, Ocean Model., 100, 109–124, https://doi.org/10.1016/j.ocemod.2016.02.002, 2016

Stoffelen, A. and Portabella, M.: On Bayesian scatterometer wind inversion, IEEE T. Geosci. Remote, 44, 1523–1533, https://doi.org/10.1109/TGRS.2005.862502, 2006.

Streßer, M., Carrasco, R., and Horstmann, J.: Video-based estimation of surface currents using a low-cost quadcopter, IEEE Geosci. Remote S., 14, 2027–2031, https://doi.org/10.1109/LGRS.2017.2749120, 2017.

Trampuz, C., Gebert, N., Placidi, S., Izzy Hendricks, L., Speziali, F., Navarro, V., Martin, A. C. H., Gommenginger, C., Suess, M., and Meta, A.: OSCAR – The airborne interferometric and scatterometric radar instrument for accurate sea current and wind retrieval, EUSAR 2018, 12th European Conference on Synthetic Aperture Radar, Aachen, Germany, 4–7 June 2018, 1–16, 2018.

Villas Bôas, A. B., Ardhuin, F., Ayet, A., Bourassa, M. A., Brandt, P., Chapron, B., Cornuelle, B. D., Farrar, J. T., Fewings, M. R., Fox-Kemper, B., Gille, S. T., Gommenginger, C., Heimbach, P., Hell, M. C., Li, Q., Mazloff, M. R., Merrifield, S. T., Mouche, A., Rio, M. H., Rodriguez, E., Shutler, J. D., Subramanian, A. C., Terrill, E. J., Tsamados, M., Ubelmann, C., and van Sebille, E.: Integrated observations of global surface winds, currents, and waves: requirements and challenges for the next decade, Front. Mar. Sci., 6, 425, https://doi.org/10.3389/fmars.2019.00425, 2019.

Wang, Z., Stoffelen, A., Fois, F., Verhoef, A., Zhao, C., Lin, M., and Chen, G.: SST Dependence of Ku- and C-Band Backscatter Measurements, IEEE J-STARS, 10, 2135–2146, https://doi.org/10.1109/JSTARS.2016.2600749, 2017.

Yurovsky, Y., Kudryavtsev, V., Grodsky, S., and Chapron, B.: Sea surface Ka-band Doppler measurements: Analysis and model development, Remote Sens.-Basel, 11, 839, https://doi.org/10.3390/rs11070839, 2019.