# Joint analysis of coastal altimetry and high-frequency (HF) radar data: observability of seasonal and mesoscale ocean dynamics in the Bay of Biscay

Ivan Manso-Narvarte, Ainhoa Caballero, Anna Rubio, Claire Dufau, Florence

Birol

► To cite this version:

Ivan Manso-Narvarte, Ainhoa Caballero, Anna Rubio, Claire Dufau, Florence Birol. Joint analysis of coastal altimetry and high-frequency (HF) radar data: observability of seasonal and mesoscale ocean dynamics in the Bay of Biscay. OCEAN SCIENCE, 2018, 14, pp.1265-1281. 10.5194/os-14-1265-2018. insu-03670734

HAL Id: insu-03670734

https://insu.hal.science/insu-03670734

Submitted on 17 May 2022

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_1.png" /></figure>

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

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

# Joint analysis of coastal altimetry and high-frequency (HF) radar data: observability of seasonal and mesoscale ocean dynamics in the Bay of Biscay

Ivan Manso-Narvarte<fnref n="1" />, Ainhoa Caballero<fnref n="1" />, Anna Rubio<fnref n="1" />, Claire Dufau<fnref n="2" />, and Florence Birol<fnref n="3" />

**Correspondence:** Ivan Manso-Narvarte (imanso@azti.es)

Received: 11 March 2018 – Discussion started: 26 March 2018

Revised: 20 July 2018 – Accepted: 23 August 2018 – Published: 16 October 2018

**Abstract.** Land-based coastal high-frequency (HF) radar systems provide operational measurements of coastal surface currents (within 1–3 m depth) with high spatial (300 m–10 km) and temporal (≤ 1 h) sampling resolutions, while the near-continuous altimetry missions provide information, from 1993 until today, on geostrophic currents in the global ocean with typical along-track and temporal sampling resolutions of > 7 km and > 9 days, respectively. During the last years, the altimetry community has made a step forward in improving these data in the coastal area, where the data present lower quality than in the open ocean. The combination of HF radar and altimetry measurements arises as a promising strategy to improve the continuous monitoring of the coastal area (e.g. by expanding the measurements made by HF radars to adjacent areas covered by the altimetry or by validating/confirming improvements brought by specific coastal algorithms or new altimeter missions). A first step towards this combination is the comparison of both data sets in overlapping areas.

In this study, a HF radar system and two Jason-2 satellite altimetry products with different processing are compared over the period from 1 January 2009 to 24 July 2015. The results provide an evaluation of the performance of different coastal altimetry data sets within the study area and a better understanding of the ocean variability contained in the HF radar and altimetry data sets. Both observing systems detect the main mesoscale processes within the study area (the Iberian Poleward Current and mesoscale eddies), and the highest correlations between radar and altimetry (up to 0.64) occur in the slope where the Iberian Poleward Current rep-

resents a significant part of the variability in the circulation. Besides, the use of an Ekman model, to add the wind-induced current component to the altimetry-derived geostrophic currents, increases the agreement between both data sets (increasing the correlation by around 10%).

## 1 Introduction

Ocean dynamics result from a combination of processes of different timescales and space scales. However, and mainly due to technical limitations, this complexity cannot be captured by the existing observational systems if each observing technique is analysed individually, since they are designed for resolving certain scales. Nowadays, there is a growing tendency to combine different observing systems for a more complete description and understanding of the ocean dynamics. Current observatories are designed to monitor, in an operational way, the ocean environment to support the human activities concentrated in the coast (Liu et al., 2015). In recent years, great effort has been focused on the development and improvement of these platforms. In the framework of European projects such as JERICO (2007–2013) and JERICO-NEXT (2014–ongoing, http://www.jerico-ri.eu, last access: 1 October 2018) progress has been made, on the one hand, regarding the improvement and creation of coastal observational networks and, on the other hand, regarding the applications of these observatories for addressing European marine policies. As an example, the latter project aims to inte-

<footnote n="1">¹AZTI-Marine Research, Herrera Kaia, Portualdea z/g, 20110 Pasaia, Spain</footnote>

<footnote n="2">²CLS Space Oceanography Division, 8–10 Rue Hermès, 31526 Ramonville-St.-Agne, France</footnote>

<footnote n="3">³LEGOS, 14 Ave. E. Belin, 31400 Toulouse, France</footnote>
grate emerging methodologies and technologies to proceed
towards the automated monitoring at a high temporal and
spatial resolution of wider areas in order to provide the best
possible data and products necessary for the implementation
of the European Marine Strategy Framework Directive. For
this purpose, there is ongoing research on assessing the in-
terconnection among physics, biogeochemistry, and biology
at different spatial and temporal scales. In addition to the
development of coastal observatories, there are global ini-
tiatives, such as the GLOBCURRENT project (2014–2017;
http://www.globcurrent.org/, last access: 1 October 2018),
aimed at advancing in the evaluation of the synergy of satel-
lite sensors and in situ data for the quantitative estimation of
ocean surface currents (e.g. Rio et al., 2014).

Among the different methodologies to retrieve surface cur-
rents, two are particularly interesting due to their high poten-
tial complementarity: satellite altimetry and land-based high-
frequency (HF) radar (HFR) systems. The former technique
consists in a constellation of altimeters onboard satellites
measuring the global sea level, with a revisit period greater
than a week and a track distance around tens of kilometres.
These continuous sea level series are today close to com-
pleting 25 years of data, resolving the ocean dynamics from
mesoscale to near-climate scale. HFRs are designed to mea-
sure the local ocean surface dynamics with a high time and
space resolution. However, altimetry and HFR do not capture
exactly the same dynamics. Altimetry detects surface cur-
rents that are in geostrophic equilibrium (by excluding the di-
rect response of the surface layer to the wind and then part of
the high-frequency variations), whereas HFRs measure sur-
face total currents, i.e. the geostrophic and ageostrophic com-
ponents (like wind-driven and inertial currents or the wave-
induced Stokes drift; e.g. Graber et al., 1997; Law, 2001;
Arduhuin et al., 2009).

Besides the effort made for collecting data from differ-
ent platforms, methods for combining these data are under
development. Recent studies focused on the evaluation of
the performance of altimetry using HFRs, concluding that
HFRs offer a way to improve the validation of altimetry prod-
ucts for coastal areas (Chavanne and Klein, 2010; Liu et al.,
2012; Pascual et al., 2015; Troupin et al., 2015; Roesler et
al., 2013). One of the most extended approaches found in the
literature to study the synergy between altimetry and HFR
data consists in the comparison of the total across-track cur-
rents in the along-track direction (e.g. Morrow et al., 2017;
Troupin et al., 2015; Pascual et al., 2015). The combination
of HFR and altimetry could help to potentiate their strengths
by, for example, expanding the spatial and temporal coverage
of the HFR systems or evaluating and correcting the altimet-
ric signal near the coast.

In this study, we focus on the south-eastern Bay of Bis-
cay (SE-BoB), which is characterized by the presence of
canyons (e.g. Capbreton canyon), by an abrupt change in the
orientation of the coast, and by a narrow shelf and slope. The
winter surface circulation in the SE-BoB is mainly related

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

**Figure 1.** Study area, observational systems and main character-
istics of the ocean circulation (figure modified from Rubio et al.,
2018). The winter IPC is represented by blue solid arrows, whereas
the blue hollow arrows show the mesoscale eddy regime (although
only anticyclonic arrows are represented, eddies of anticyclonic
and cyclonic polarity are observed in different locations along the
slope). The bold black lines delimit the HFR total-current footprint.
The black stars represent the HFR stations: Matxitxako (left) and
Higer (right). Jason-2 tracks 213 and 248 are represented by black
crosses and the part of the track used in this study is marked in red.
Grey lines: 1000, 3000, and 4000 m isobaths.

to the Iberian Poleward Current (IPC), which affects the up-
per 300 m of the water column. In winter, the IPC flows over
the slope, advecting warm surface waters (Le Cann and Ser-
pette, 2009; Charria et al., 2013) eastwards along the Span-
ish coast and northwards along the French coast (Fig. 1). In
summer, the flow is reversed, being 3 times weaker than in
winter (Solabarrieta et al., 2014). Overlaid to the density-
driven slope circulation, wind-induced currents are the main
drivers of the surface circulation in the area (e.g. Lazure,
1997; Solabarrieta et al., 2015). During autumn and winter,
south-westerly winds dominate and generate northward and
eastward drift over the shelf. The wind regime changes to
the NE during spring, when it causes sea currents to turn to-
ward the W–SW along the Spanish coast. The summer sit-
uation is similar to that of spring, but the weakness of the
winds and the greater variability in the direction of the gen-
eral drift make currents more variable (González et al., 2004;
Lazure, 1997; Solabarrieta et al., 2015). In addition to these
processes, mesoscale eddies in the SE-BoB are generated,
mainly during winter, by the interaction of the IPC with the
abrupt bathymetry (Pingree and Le Cann, 1992) (Fig. 1). The
combination of these processes makes the SE-BoB an area of interesting complexity.

The existence of a long historical time series of surface current fields from a long-range HFR system in the area provides an invaluable opportunity to explore the benefit of a combined analysis of satellite and land-based remote-sensing ocean currents. This HFR network (two sites; see Fig. 1) is part of the coastal observatory of the SE-BoB, also composed of a network of oceano-meteorological coastal stations and two slope buoys. The performance of this system and its potential for the study of ocean processes and of transport patterns in the area have already been demonstrated by previous work (e.g. Solabarrieta et al., 2015; Rubio et al., 2018). With regard to the usefulness of altimetry for describing ocean dynamics in the BoB, several studies have proven its suitability to study processes that go from mesoscale (Dussurget et al., 2011; Herbert et al., 2011; Caballero et al., 2008b, 2014, 2016) to climate scale (e.g. Pingree, 2005).

The main objectives of this study are first to obtain a diagnosis of the agreement of the surface currents measured by altimetry and HFR over the SE-BoB and, second, to evaluate the observability of certain mesoscale processes by both measuring systems.

## 2 Data and methods

### 2.1 Data

#### 2.1.1 HFR data

HFRs are remote-sensing instruments that send radio waves to the ocean surface and use the signal backscattered by the waves to infer the radial velocity of the surface current (toward or away from each HFR antenna). They can measure surface currents over wide areas with high spatial (300 m–10 km) and temporal (≤ 1 h) resolution. In this study, surface currents were obtained by means of two long-range HFR antennae. These antennae emit at a central frequency of 4.5 MHz and with an operational 30 kHz bandwidth. They are located at Matxitxako and Higer capes (Fig. 1) and have provided operational data since 2009 (with some interruptions mostly due to maintenance stops or malfunctioning related to severe atmospheric conditions). The averaged Doppler backscatter spectrum obtained from the received signal (in a 3 h window) is processed to obtain hourly radial currents using the MUSIC algorithm (Schmidt, 1986). The coverage of radial data is up to 150 km, and the range cell and angular resolution are set to 5 km and 5°, respectively. Radial data are quality controlled using advanced procedures based on velocity and variance thresholds, signal-to-noise ratios, and radial total coverage. Since the deployment of the HFR system, the receiving antenna pattern of the two HFR sites has been calibrated at least every 2 years. A more detailed description of the system and of the HFR data validation exercises are provided by Solabarrieta et al. (2014, 2015, 2016) and Rubio et al. (2011, 2018).

To obtain total currents gridded onto a 5 km resolution regular orthogonal mesh, a least mean square algorithm (spatial interpolation radius of 10 km) was applied by using the HFR_Progs Matlab package (https://cencalarchive.org/~cocmpmb/COCMP-wiki/index.php/Main_Page#HFR_Progs_Documentation, last access: 1 October 2018), based on Gurgel (1994) and Lipa and Barrick (1983). Then, using the same grid, radial velocities were processed with HFR_Progs to generate spatially gap-filled open-mode analysis (OMA) total currents (Kaplan and Lekien, 2007). In order to generate hourly total fields, 85 OMA modes, built setting a minimum spatial scale of 20 km, were used. A first analysis of the comparisons between HFR and altimetry showed that the results obtained using total currents generated by least mean squares and OMA were very similar. Thus, only results using OMA currents are presented in this work.

The typical spatial scales resolved by the HFRs depend on the resolution of the data and thus mainly on the operation frequency of the systems (Rubio et al., 2017). For the SE-BoB, the spatial scales resolved are typically of O(15–20) km.

#### 2.1.2 Altimetry data

The basic principle of the altimetry technology is to send a radar signal to the sea surface and then to measure the reflected return echo. The time needed for the signal to go and come back determines the distance between the altimeter and the sea surface (called the range). A physically based model (Brown, 1977) is adjusted to the resulting signal, called waveform, providing different parameters, including range. To reduce the measurement noise, the result is averaged and the final data rate is classically (and in our case) 1 Hz (i.e. one datum every ~7 km along the satellite track). By subtracting the range to the satellite orbit altitude (with respect to the reference ellipsoid) and by applying several corrections (e.g. tropospheric and ionospheric effects on the radar wave, sea surface bias), the sea surface height (SSH) is obtained (SSH = orbit – range – corrections). To retrieve the total geostrophic currents, in balance with the SSH gradients, the surface height must be referenced to the geoid; this height is called dynamic topography (DT). However, since geoids are not known with enough accuracy, only geostrophic anomaly currents can be accurately derived from sea level anomaly (SLA). The SLA is referenced to a temporal average (Le Traon et al., 2003) by subtracting a mean profile to the SSH (SLA = SSH – (SSH)). For more information about the SLA and the reference surfaces, the reader is referred to https://www.aviso.altimetry.fr/en/techniques/altimetry/principle/basic-principle.html (last access: 1 October 2018).
The two different along-track SLA time series used herein come from Jason-2's 248 and 213 tracks from cycle 18 to 259 and with a revisit period of ~10 days. Track 248 covers the HFR footprint area, whereas track 213 only crosses a small area at the NW of the HFR total-current coverage (Fig. 1); therefore, track 213 is only useful for the comparison with HFR radial data.

One of the altimetry data sets used in this study is the CTOH (Center for Topographic studies of the Ocean and Hydrosphere)-XTRACK product (http://ctoh.legosobs.mip.fr/products/coastal-products/coastal-products-1/sla-1hz, last access: 1 October 2018), which provides SLA data specifically processed for coastal areas. The filtered version of the product was used, meaning that a 40 km cut-off spatial Loess filter (Cleveland and Devlin, 1988) was also applied in order to reduce the remaining noise in the along-track SLA. For the same Jason-2 altimeter measurements, a differently processed data set was also used to assess possible discrepancies between the two altimetric products. This product is a preliminary version of the 2018 update of the reprocessed global ocean along-track level-3 data provided by the Copernicus Marine Environment Monitoring Service (CMEMS). As for the CTOH-XTRACK product, data filtered with a 65 km cut-off spatial Lanczos filter (Pujol et al., 2016) were used. Note that, on the one hand, the accuracy of altimetry data is lower in the 20–30 km coastal band, so it might be a source of differences between altimetry and HFR data, and that, on the other hand, the larger the oceanic signal (larger signal-to-noise ratio), the lower this effect will be observed to be. In the SE-BoB, part of the slope (characteristically narrow) is located within the 20–30 km coastal band; therefore, the altimetry data will be affected by this inaccuracy in this area.

### 2.1.3 Wind data and sea surface temperature images

Hourly wind data from the Weather Research and Forecasting model (WRF; http://mandeo.meteogalicia.es/thredds/catalogos/WRF_2D/catalog.html, last access: 1 October 2018) were provided by the meteorological agency of Galicia (MeteoGalicia). This model, with a native resolution of 12 km, reproduces the offshore wind fields of the SE-BoB (Ferrer et al., 2010) with reasonable accuracy. In this study, the WRF gridded fields were interpolated to the Jason-2 along-track points. Due to the cloudy weather in the SE-BoB, the most appropriate infrared sea surface temperature (SST) images were selected one by one from the AVHRR (Advanced Very High Resolution Radiometer) sensor series of 1 km resolution to process level-2 SST maps.

### 2.2 Methods

Since the time resolution of the altimetry (~10 days) is lower than that of the HFR (hourly), the HFR current data were filtered using a low-pass filter, based on a 10-day running average. The objective was to remove the high-

frequency (HF) signals contained in the data (see, for instance, Solabarrieta et al., 2014) and part of the ageostrophic signals to make the HFR data closer to the measurements from altimetry. Several tests were carried out to determine the sensibility of the results to the temporal window chosen for the filter. For that purpose, 2-, 5-, 10- and 15-day low-pass filter configurations were tested, and the resulting currents were compared to those derived from the SLA. The 10-day filtered HFR currents provided the highest correlation with the altimetry data. It is worth noting that the 10-day running average filtered out a significant part of the wind-induced currents, except for the low-frequency (LF) Ekman component.

Since from the altimetry data used here we can only obtain sea surface anomaly currents, as explained in Sect. 2.1.2, the comparisons with the radar data were carried out in terms of anomaly. In order to obtain the HFR anomaly currents (AC<sub>HFR</sub>), the temporal average of the HFR currents for the study period was subtracted to the series of low-pass HFR currents (Eq. 1). Then, to obtain SLA relative to the same period, the average SLA for the study period was subtracted to the SLA series (Eq. 2).

$$AC_{HFR} = V - \langle V \rangle_{t_1-t_2} \quad (1)$$

$$SLA' = SLA - \langle SLA \rangle_{t_1-t_2} \quad (2)$$

Hence, the time referenced SLA' and AC<sub>HFR</sub> were obtained as suggested in Pujol et al. (2016), where $t_1 = 1$ January 2009 and $t_2 = 24$ July 2015. Hereafter, when referring to AC<sub>HFR</sub> computed using Eq. (1), we will use AC<sub>HFR,R</sub> to refer to current anomalies computed directly from the radial HFR components and AC<sub>HFR,T</sub> to refer to current anomalies computed from total OMA HFR currents.

For the altimetry velocity, across-track geostrophic anomaly currents (AC<sub>G</sub>) were inferred by means of the finite-difference geostrophic velocity equation (Eq. 3):

$$AC_G = -\frac{g \Delta SLA'}{f \Delta x}, \quad (3)$$

where g is the acceleration of gravity, f is the Coriolis parameter, x is the along-track distance, and SLA' is the time referenced SLA. AC<sub>G</sub> was estimated along the altimeter track by a three-point central difference operator.

For the statistical comparison between AC<sub>G</sub> and AC<sub>HFR</sub>, two different strategies were used. The first and simplest approach is the comparison of HFR radial across-track currents with altimetry across-track currents at two given points (Sect. 2.2.1), which permits the direct use of radial HFR currents, that is, AC<sub>HFR,R</sub>. The second approach is the along-track comparison (Sect. 2.2.2) of total OMA HFR and altimetry across-track currents, which provides additional information on the spatial variability of the agreement between both data sets. Finally, in order to take into account the LF Ekman component that remains in the low-pass-filtered HFR data, a model for the computation of this component was used (Sect. 2.2.3).
<figure><img src="image_5.png" /><figcaption>Figure 2. Schematic view of the pointwise comparison and of the data used for this approximation. Jason-2's 213 and 248 tracks are depicted by orange lines. The HFR radial directions from both sites are represented by grey lines and the selected radial directions (from the Matxitxako site) for the pointwise comparison are plotted in red (the central radial orthogonal to the track) and in blue (the adjacent radials). Points E and W are those along both tracks, where each of the HFR radial directions are orthogonal to the Jason-2's track directions. At point E, the blue and red crosses show the selected points of HFR radial directions where radial currents are averaged in a first step to obtain the corresponding values on the track (green crosses). Then, the along-track green crosses are averaged to obtain the corresponding HFR current values at the point at issue (the same process is carried out for point W, but it is not represented in the figure). Grey lines: 200, 1000, and 2000 m isobaths.</figcaption></figure>

## 2.2.1 Pointwise comparison

This method, previously applied in Liu et al. (2012), consists in a direct comparison between HFR and altimetry data at a certain point, where one of the HFR radial directions (red lines in Fig. 2) crosses the altimeter track perpendicularly. This approximation allows us to directly use the radar radial currents, which are in the same direction as the across-track $AC_G$. This method was applied for Jason-2 altimeter tracks 213 and 248 shown in Fig. 1, by using the corresponding $AC_{HFR,R}$ measured by the Matxitxako antenna at two points: point E on track 248 and point W on track 213 (Fig. 2). Note that the radials from Higer station were not used, since they did not provide any orthogonal radial direction for track 248 and were too gappy for track 213 due to the large distance between this track and the antenna. Point E is located near the Capbreton canyon at 43.75° N, 2.05° W, 46.86 km from the coast and at a depth of 500 m (on the slope), whereas point W is located at 43.80° N, 3.58° W, 40.73 km from the coast and at a depth of 3000 m (on the abyssal plain). It is worth noting that by using directly radial

currents, additional errors that propagate in the combination of HFR radials into HFR OMA currents are avoided.

In order to make the computations more robust to the potential absence of HFR or altimetry data at points W and E, nearby points were considered to obtain the across-track currents (Fig. 2). With regard to $AC_{HFR,R}$, radial directions of ±5°, away from the orthogonal radial direction were also considered. For each of them, the points in a range of ±5 km from the track along the central and adjacent radials were selected. Then, the across-track $AC_{HFR,R}$ in our point was obtained by firstly averaging the values for each radial, so that only three values along the track were obtained. Finally, the three values were again averaged to get the corresponding across-track $AC_{HFR,R}$.

On the other hand, the $AC_G$ in the three along-track points considered for the HFR case were averaged to obtain the $AC_G$ in the central point (E and W). This permitted us to ensure a similar spatial smoothing for both data sets. It must be mentioned that at point W and for the CTOH altimetry product, the points were located between the along-track HFR points, so instead of considering three along-track points, two points on each side of the central point (W) were selected.

## 2.2.2 Along-track current comparison

In order to assess the variability in the comparison between HFR and altimetry from the coast to the open ocean, the comparison between the across-track $AC_{HFR,T}$ and $AC_G$ along track 248 was performed. For that purpose, $AC_{HFR,T}$ was interpolated into the along-track altimetry points, and it was rotated to the across-track direction. Then, for each along-track $AC_G$ point, the average with its four adjacent points was calculated. As in the previous case, this permitted us to ensure a similar spatial smoothing for both data sets. The sensitivity to the number of adjacent points considered was tested, and this approach was the one that provided the best adjustment to the HFR data.

## 2.2.3 Ekman currents

Ekman currents were estimated to evaluate what their contribution to LF currents in the area was, and how this component contributed to part of the differences observed between HFR and altimetry. Three different ways of calculating Ekman currents were tested to infer which one provided the best results in the comparisons: the rule of thumb that states that the surface currents are 3% of the wind velocity; Ekman equations for the surface (Ekman, 1905); and the model M1 proposed in Rio and Hernandez (2003). Ultimately, the M1 model offered the best results:

$$ \mathbf{u}_{ek} = \frac{b\tau}{\sqrt{f}} e^{i\theta}, \qquad (4) $$

where f is the Coriolis parameter, $\mathbf{u}_{ek}$ is the Ekman current vector (meridional and zonal components), b is the amplitude
<table><caption>**Table 1.** Statistics of different points for the study period. *r* is the correlation parameter, and RMSD is the root mean square difference of the velocity anomalies between HFR and altimetry for different altimetry products (i.e. CMEMS and CTOH). The mean and the standard deviation (SD) are also estimated for each data set. *E* and *W* are the points where HFR radial velocities are orthogonal to tracks 248 and 213, respectively (see Fig. 1). *E*<sub>T</sub> means that the parameters are estimated at point *E* but using the HFR current velocity anomaly fields rotated to across-track direction instead of using the HFR radial velocity anomalies. That HFR data have also been used for Max<sub>T</sub>, which is the track point where the correlation *r* is maximal (point 2.10° W, 43.82° N for CTOH and point 2.09° W, 43.80° N for CMEMS). RAD refers to HFR data, which really means radial anomaly current for *E*<sub>R</sub> and *W*<sub>R</sub> and across-track rotated OMA anomaly currents for *E*<sub>T</sub> and Max<sub>T</sub>. There is no data in Max<sub>T</sub> because they are different for each altimetry product.</caption><thead><tr><th rowspan="2"></th><th rowspan="2"></th><th colspan="2">r</th><th colspan="2">RMSD (cm s<sup>-1</sup>)</th><th colspan="3">Mean &plusmn; SD (cm s<sup>-1</sup>)</th></tr><tr><th>CMEMS</th><th>CTOH</th><th>CMEMS</th><th>CTOH</th><th>CMEMS</th><th>CTOH</th><th>RAD</th></tr></thead><tbody><tr><td rowspan="2">E<sub>R</sub></td><td>AC<sub>G</sub></td><td>0.53</td><td>0.48</td><td>7.4</td><td>7.5</td><td>0.3 &plusmn; 8.7</td><td>-0.1 &plusmn; 8.5</td><td rowspan="2">0.1 &plusmn; 5.9</td></tr><tr><td>AC<sub>G</sub> + AC<sub>E</sub></td><td>0.64</td><td>0.59</td><td>7.1</td><td>7.4</td><td>0.3 &plusmn; 9.7</td><td>-0. &plusmn; 9.4</td></tr><tr><td rowspan="2">W<sub>R</sub></td><td>AC<sub>G</sub></td><td>0.61</td><td>0.64</td><td>9.2</td><td>7.8</td><td>-0.2 &plusmn; 11.9</td><td>0.2 &plusmn; 9.9</td><td rowspan="2">0.0 &plusmn; 9.3</td></tr><tr><td>AC<sub>G</sub> + AC<sub>E</sub></td><td>0.67</td><td>0.70</td><td>8.9</td><td>7.4</td><td>-0.1 &plusmn; 12.8</td><td>0.2 &plusmn; 10.6</td></tr><tr><td rowspan="2">E<sub>T</sub></td><td>AC<sub>G</sub></td><td>0.56</td><td>0.53</td><td>7.1</td><td>6.5</td><td>-0.1 &plusmn; 8.2</td><td>-0.3 &plusmn; 7.6</td><td rowspan="2">0.2 &plusmn; 5.4</td></tr><tr><td>AC<sub>G</sub> + AC<sub>E</sub></td><td>0.65</td><td>0.62</td><td>7.1</td><td>6.4</td><td>-0.1 &plusmn; 9.3</td><td>-0.4 &plusmn; 8.6</td></tr><tr><td rowspan="2">Max<sub>T</sub></td><td>AC<sub>G</sub></td><td>0.60</td><td>0.55</td><td>6.5</td><td>6.1</td><td>-0.2 &plusmn; 7.9</td><td>-0.2 &plusmn; 7.2</td><td rowspan="2">-</td></tr><tr><td>AC<sub>G</sub> + AC<sub>E</sub></td><td>0.68</td><td>0.65</td><td>6.6</td><td>6.0</td><td>-0.3 &plusmn; 8.9</td><td>-0.3 &plusmn; 8.3</td></tr></tbody></table>

parameter, $\theta$ is the phase parameter, and $\tau$ is the wind stress vector (meridional and zonal components) that was obtained by means of the bulk–flux formula (Stewart, 2009):

$$ \tau = \rho_{\text{air}} C_d |\mathbf{w}| \mathbf{w}, \qquad (5) $$

where $\rho_{\text{air}}$ is the density of the air (1.22 kg m<sup>-3</sup>), $\mathbf{w}$ is the wind vector, and $C_d$ is the drag coefficient proposed by Large and Pond (1981). *b* and $\theta$ were acquired by adjusting the model (by a least square fit) to the real data in the study area. These parameter values were taken from Calabiero et al. (2008a), where from the fit of wind stress measurements and drifter-derived currents of the Bay of Biscay to Eq. (4), the following results were obtained: $b = 4.45 \times 10^{-3}$ m<sup>2</sup> kg<sup>-1</sup> s<sup>1/2</sup> and $\theta = -23.68^{\circ}$. Once the parameters were obtained, $\mathbf{u}_{\text{ek}}$ was estimated using the wind data series described in Sect. 2.1.3.

Ekman currents initially computed in the locations of the WRF model nodes were interpolated and rotated (from zonal and meridional directions to along-track and across-track directions). For the pointwise comparison, they were interpolated in *E*, in *W*, and at their adjacent points and rotated to obtain the across-track component. Then, they were averaged to obtain the across-track Ekman current velocity at each point (*E* and *W*). In the along-track current comparison, they were interpolated to the altimetry along-track points and, then, rotated to get the across-track component. Ekman currents were also 10-day running averaged to remove the HF Ekman signal. For the comparisons with $AC_{\text{HFR}}$, Ekman anomaly currents ($AC_E$) were obtained by subtracting the average value of the study period like in Eqs. (1) and (2). Then, they were added to $AC_G$, and this is indicated as $AC_G + AC_E$ in Table 1.

The results of all the comparisons described above are presented in terms of the correlation coefficient or correlation (*r*) (with a confidence level of 100% and 90% for the pointwise comparison and the along-track current comparison, respectively) and the root mean square difference (RMSD) between the across-track currents measured by each system. In addition, the mean and the standard deviation (SD) of such currents were also analysed. All these parameters were computed for the study period. Since the HFR radials are used at the pointwise comparison, the crossing points are called $E_R$ and $W_R$. In the along-track comparison, the point *E* is precisely one of those along-track points. Therefore, considering that the rotated $AC_{\text{HFR},T}$ are used in this comparison, the point is named $E_T$. The points along the track where *r* is maximal (point 2.10° W, 43.82° N for CTOH and point 2.09° W, 43.80° N for CMEMS) are also called $Max_T$. All these points and their statistics are displayed in Table 1, as well as in Figs. 3–5, and they are also discussed in the next section. Moreover, in Fig. 5, apart from the results of $E_T$ and $Max_T$, the results for all of track 248 are also shown.

Note that the HFR–altimetry comparisons were carried out for the CMEMS and CTOH data sets and that each comparison was also made with and without adding $AC_E$.

# 3 Results and discussion

## 3.1 Statistical results

Table 1 provides an overview of the statistical results of the HFR and altimetry data set comparison at points $E_R$ and $W_R$ for the pointwise method and at points $E_T$ and $Max_T$ for the along-track method. Despite the differences between the
<figure><img src="image_6.png" /><figcaption>**Figure 3.** Across-track $AC_G + AC_E$ and $AC_{HFR,R}$. (a) CTOH, CMEMS, and HFR data sets at point $E_R$. (b) CTOH-HFR and CMEMS-HFR residuals at point $E_R$. (c) CTOH, CMEMS, and HFR data sets at point $W_R$. (d) CTOH-HFR and CMEMS-HFR residuals at point $W_R$. Black arrows depict the slope current intensifications mentioned in the text.</figcaption></figure>

two measuring systems, the comparison between $AC_{HFR}$ and $AC_G$ ($AC_{HFR}$ and $AC_G + AC_E$) shows maximum correlations of 0.64 (0.70). Concerning the pointwise comparisons between $AC_G$ and $AC_{HFR,R}$ it can be observed that $r$ depends on the location as well as on the altimetry product considered. In any case, the addition of $AC_E$ increases $r$ by 6%–11% for both altimetry products: from 0.53 and 0.48

to 0.64 and 0.59 in $E_R$ and from 0.61 and 0.64 to 0.67 and 0.70 in $W_R$ for CMEMS and CTOH, respectively. The RMSD also decreases between 0.1 and 0.4 cm s⁻¹. Therefore, in general terms, the addition of $AC_E$ decreases the differences between $AC_{HFR,R}$ and $AC_G$. However, it adds variability: the SD for $AC_G + AC_E$ increases by 0.7–1.0 cm s⁻¹. It should be pointed out that the SD of $AC_{HFR,R}$ is lower than
the SD of $AC_G$ (where CMEMS's variability is higher than CTOH's), probably due to the filtering of HF signals applied to the radar or because the finite-difference operator that estimates the geostrophic velocity from SLA' increases the variability. The same effect was observed in Liu et al. (2012) and in Verron et al. (2018). Nevertheless, although the addition of $AC_E$ increases the variability even more, it improves the results in terms of correlation and RMSD.

Figure 3 displays the time series of $AC_{HFR,R}$ and $AC_G + AC_E$ at points $E_R$ and $W_R$. In general, there is global agreement between the three series ($AC_G$ (CTOH) + $AC_E$, $AC_G$ (CMEMS) + $AC_E$, $AC_{HFR,R}$), and all the data sets observe the same variability. Although in Fig. 3 it is difficult to detect any differences between locations and data sets, the lowest RMSD is observed in $E_R$ for CMEMS $AC_G$ (RMSD = 7.1 cm s⁻¹), while the RMSD rises up to 8.9 cm s⁻¹ in $W_R$ for CTOH. A larger variability in $W_R$ is shown in Table 1, which is also noticeable in Fig. 3a and c.

In terms of correlation, the results suggest greater agreement in $W_R$ (0.67 and 0.70 for CMEMS and CTOH, respectively) than in $E_R$ (0.64 and 0.59 for CMEMS and CTOH, respectively). This could be explained by the fact that point $W_R$ is located in a deeper area and is more influenced by the IPC that flows over the slope. The IPC shows lower variability at point $E_R$ than at point $W_R$, where the signal of the slope current is more persistent. These spatial differences agree with what was observed in the area by Rubio et al. (2009) from the analysis of the ADCP (Acoustic Doppler Current Profiler) time series of two buoys located over the slope in locations comparable to points *E* and *W*. These authors suggest that while there is a clear along-slope transport with intense mesoscale variability at Matxitxako buoy (slightly east from point *W*), at Donostia buoy (close to point *E*) the influence of the slope circulation was less significant. The lower current velocities and lower vertical coherence observed at Donostia buoy during winter could be linked to the complex bathymetry, which might force the IPC to flow over deeper grounds out of the point measured by the buoy. This could also explain why the addition of $AC_E$ increases more *r* in $E_R$ (by 11%) than in $W_R$ (by 6%), where the circulation has a stronger geostrophic component.

With regard to the performances of the two altimetry products, it must be highlighted that CTOH shows higher (lower) *r* (RMSD) in $W_R$, while CMEMS shows higher (lower) *r* (RMSD) in $E_R$. However, the differences are small and do not permit us to draw conclusions on their relative accuracy.

In $E_R$ ($W_R$) when *r* is higher for CMEMS (CTOH) than for CTOH (CMEMS), the RMSD is lower in the former. Therefore, the higher the correlation, the smaller the difference between $AC_G$ and $AC_{HFR,R}$. However, the difference in the RMSD at each point (i.e. between the RMSD value of each altimetry product with the radar) is not related to the difference in *r*.

Figure 3b and d also show the residuals between $AC_G$ and $AC_{HFR,R}$ for each altimetry product. It can be observed

that in $W_R$ (Fig. 3d), the residuals' amplitudes are larger for CMEMS than for CTOH, agreeing with the higher RMSD value observed for CMEMS $AC_G$. In $E_R$ the similar RMSD values observed for both altimetry products agree with the similar amplitudes of the residuals in Fig. 3b.

The IPC winter intensification is visible in all data sets, being stronger in $AC_G$ for both altimetry products. There are some remarkable intensifications, for instance in November 2014, when a strong peak is shown in all the series and where the $AC_{HFR,R}$ signal is higher in $W_R$ than in $E_R$ (see black arrows in Fig. 3a and c), approaching $AC_G$. There is another remarkable intensification in winter 2009, when the IPC shows a pronounced peak (see black arrows in Fig. 3a and c), especially in $E_R$, where CTOH is very similar to $AC_{HFR,R}$. This intensification is not so clear at point $W_R$, but it is still noticeable compared to the rest of the period. Equatorward slope current intensifications can be also observed through the whole period, as for instance in May 2011, where the peak is more prominent in $W_R$ than in $E_R$ (and even more for CMEMS) (see black arrows in Fig. 3a and c).

Since the presence of a stronger IPC signal is expected to improve the correlation between HFR and SLA data sets and the IPC shows marked seasonality, a monthly analysis has been carried out (Fig. 4). The monthly values of the statistical parameters shown in the figure, have been computed considering all the available data for that month during all the study period.

It can be observed that in terms of monthly mean currents, the three time series have the same tendency and that in general there is low discrepancy among them. No significant differences in terms of monthly patterns are observed among the two altimetry products. The winter poleward current intensification is evident from October to January with a maximum in November (ranging for all data sets from 7.5 to 13.4 cm s⁻¹ in $E_R$ and from 11.7 to 14.8 cm s⁻¹ in $W_R$). In $E_R$, the current intensification is still perceptible in February, whereas in $W_R$ there is a small increase in March in $AC_G$ but not in $AC_{HFR,R}$. From March to September the mean is maintained around zero, with a slight decrease reaching the minimum in September (ranging from -7.3 to -3.9 cm s⁻¹ in $E_R$ for all data sets and from -5.4 to -3.6 cm s⁻¹ in $W_R$) and showing an equator-ward weak mean slope current for that period. The addition of $AC_E$ slightly strengthens the intensity of the slope currents for both poleward and equatorward directions. This can be explained by the general wind patterns of the area, which are in agreement with the main local geostrophic regime, although winter south-westerlies are stronger than the summer north-easterlies (Herbert et al., 2011).

As has been mentioned previously, Figure 4 shows lower SD values for $AC_{HFR,R}$ along all the period. These values are slightly increased by the addition of $AC_E$, especially in winter when winds are stronger. At point $E_R$, the tendency of the SD is similar to that of the mean, being slightly higher in January and then stable (with small oscillations) until September.
<figure><img src="image_7.png" /><figcaption>**Figure 4.** Monthly statistical parameters of the comparison between $AC_G$ vs. $AC_{HFR,R}$ and $AC_G + AC_E$ vs. $AC_{HFR,R}$ at points $E_R$ and $W_R$ for both altimetry products: (a) CTOH $AC_G$ at point $W_R$, (b) CTOH $AC_G$ at point $E_R$, (c) CMEMS $AC_G$ at point $W_R$, and (d) CMEMS $AC_G$ at point $E_R$. The normalized RMSD is the quotient between the RMSD and the values of the corresponding $AC_{HFR,R}$.</figcaption></figure>
<figure><img src="image_8.png" /><figcaption>**Figure 5.** Mean and variance along track 248 of the different data sets – $AC_{HFR,T}$, $AC_G$ (CTOH) and $AC_G$ (CMEMS) – and correlation and RMSD between (a) $AC_{HFR,T}$ and $AC_G$ and (b) $AC_{HFR,T}$ and $AC_G + AC_E$. The grey-coloured area corresponds to the slope between the 200 and 1000 m isobaths. The black line is the location of the point *E*, and the first point of the track is around 11 km away from the coast. Correlation values are not plotted for confidence levels under 90% (i.e. at the points along the track with a distance from the first point beyond 120 km). Note that the first point is different for CTOH and CMEMS (which are around 10.5 and 13 km away from the coast, respectively) and that the spacing between the points is slightly different. Moreover, the first point in CTOH is removed because it is an outlier. The radar points are the same as CMEMS points.</figcaption></figure>

In the last 3 months of the year, it is increased. Therefore, there is a higher variability in late autumn and winter, probably due to the slope current intensification and the stronger winds.

At point $W_R$, the highest variability takes place in the first 4 months of the year. This increase is coherent with the intensification of the slope current and the development of an anticyclonic structure in March and especially in April near Torrelavega canyon (Caballero et al., 2014). Afterwards, the variability is practically maintained, with small oscillations and an increase in CMEMS data.

Regarding the RMSD, the patterns are similar to those of the SD of $AC_G$, but with variability. In general, the addition of $AC_E$ improves the results by a slight decrease in the RMSD. The major differences between $AC_G$ and $AC_G + AC_E$ time series are observed in February and April at point $E_R$. This can be clearly seen in the normalized RMSD,

which measures the value of the RMSD with respect to the values of $AC_{HFR,R}$.

The same statistical parameters computed along track 248 are shown in Fig. 5 in addition to the correlation in order to study the spatial variability in the comparison between $AC_{HFR,T}$ and $AC_G$. Temporal statistics considering all the study period for each point of the track are plotted as a function of the distance to the first point of the track.

$AC_{HFR,T}$ show mean values close to zero along all the track, with low variability. Currents are oriented poleward over the shelf and upper slope over grounds shallower than 1000 m (i.e. until the point where the track leaves the slope current area). From there on, the mean is oriented equatorward. The mean $AC_G$ is also close to zero; however, it shows larger variability, changing between positive and negative values along the track, with a lack of agreement between both altimetry products at some points. Nevertheless, it should
be pointed out that in any case the mean anomaly currents are very weak and that the high SD values compensate for changes in the mean values. The addition of $AC_E$ does not cause any spatial variation, and it barely changes the values.

With regard to the variability, it is higher close to the coast. For $AC_G$, it slowly decreases as it gets away from the first point of the track, until the 1000 m isobath of the slope is reached (where the grey area ends in the figure). From there on, it is kept almost constant with a slight local maximum at around 120 km from the first point. For $AC_{HFR,T}$, the variability decreases until the point $E_T$ is reached (where the Capbreton canyon area ends), and afterwards, it is also almost constant, with two maximums at around 60 and 105 km from the first point.

It can be once again observed that the addition of the $AC_E$ slightly increases the variability (Fig. 5), that the variability in $AC_{HFR}$ is lower than that of $AC_G$, and that CMEMS's variability is higher than CTOH's. Moreover, in Fig. 5, all these results are provided for all the track.

The highest *r* is observed between the 200 and 1000 m isobaths of the slope, where the IPC signal is the strongest, and thus the geostrophic component measured by the HFR is also stronger. At the points furthest from the coast, the *r* decreases. This can be linked to the absence of a strong and persistent geostrophic component and a higher signal-to-noise ratio for the HFR data (which increases as we get away from the antennae).

The maximum (minimum) values of *r* (RMSD) occur at around 37 and 45 km from the first point of the track for CMEMS and CTOH, respectively (Max$_T$ points). These values could be explained by the fact that those points are located in the middle of the slope, where the slope current is stronger and where they are out of the Capbreton canyon area. At the same time, in that area, the slope current direction is nearly orthogonal to the track, so that the across-track component is stronger. For CTOH, the values around the maximum are relatively high, that is why the maximum is not a prominent peak. The same happens for CMEMS, but with a sharper peak and higher value (see Table 1). The addition of $AC_E$ increases the *r* by 8%–10% as can be seen in the table (for $E_T$ and Max$_T$ points).

In general, the addition of $AC_E$ slightly increases (decreases) the *r* (RMSD), along all the track. This fact is not perceptible in the figure, where it seems that the addition of $AC_E$ does not make any difference, but it can be observed in the values of Table 1, except at point $E_T$ for CMEMS data where it does not change.

## 3.2 Observability of mesoscale processes in HFR and altimetry data sets

In order to provide a complementary insight into the synergies and differences between HFR and altimeter data, in this section the observability of different processes detected by HFR and altimetry is qualitatively analysed. Since the data

are spatially filtered (for $AC_G$) or time filtered (for $AC_{HFR}$ and $AC_E$), the detectable processes are mesoscale, seasonal, and interannual processes, such as mesoscale eddies and the IPC. Only CMEMS data are used for this analysis, since the statistical results are very similar for both altimetry products and CMEMS data have less data gaps in the period and study area, which is more suitable for monitoring ocean processes.

### 3.2.1 Observability of the IPC

Along-track values of SLA', across-track $AC_E$ and across-track $AC_{HFR,T}$ for all the study period are shown in Fig. 6. The highest SLA' values are observed during late autumn and winter, whereas the lowest ones are observed in spring and summer, especially from March to July (Fig. 6a). This is coherent with the contribution of the main driving factor of the seasonal SLA variability in the area, the steric effect. It was observed in Caballero et al. (2008b) that the maximum amplitude of sea level in the BoB is reached in October, whilst the minimum takes place in April. $AC_E$ shows a poleward seasonality with intensifications mainly in autumn and winter (usually from November to February) and weaker equator-ward currents in spring and summer (usually from March to October) (Fig. 6b). This fact agrees with the general wind pattern in the area. Along-track SLA' gradients indicate winter slope current intensifications (IPC events) mostly from November to January and from the coast to the 1000 m isobath, approximately. The poleward intensification in winter 2014/2015 is the most remarkable, as already described in Rubio et al. (2018). In spring and summer the gradients are weaker and even suggest equator-ward currents along all the track.

The IPC events are also detected by $AC_{HFR,T}$ in late autumn and winter, mainly in the points nearest to the coast between the 200 and 1000 m isobaths (Fig. 6c). In spring and summer, although there are also several poleward current pulses, they are weaker. During this period, equator-ward current pulses are also observed. From the joint analysis of SLA' and $AC_{HFR,T}$ data series, four main IPC events can be detected along the study period (Fig. 6a and c).

More details on these events are provided in Fig. 7, where four selected HFR total current field (obtained from OMA as explained in Sect. 2.1.1) snapshots are shown. Although each event is presented for a specific date, they last around 2–3 weeks (not shown), with the dates displayed in the figure being representative of all the period. Note that the SST maps do not show the same date as HFR snapshots and SLA' data, due to the limitations of this technique under cloudy conditions. For the four events, HFR total currents show a typical IPC spatial pattern, with poleward circulation along the slope intensifying between the 200 and 1000 m isobaths (Solabarrieta et al., 2015). The poleward patterns observed by the HFR agree with the $AC_G$ observed along the altimetry track, which shows that poleward currents intensified over the slope. For the four events, the SST images show that the current intensi-
<figure><img src="image_9.png" /><figcaption>**Figure 6.** Time evolution from 1 January 2009 to 1 July 2015 along track 248 (y axis shows distance to the first point of the track in kilometres) of (a) CMEMS SLA' (cm), (b) $AC_E$ (cm s$^{-1}$), and (c) $AC_{HFR}$ (cm s$^{-1}$). The distance from the first point of the track to the coast is around 6 km. The horizontal black lines delimit the slope area between the 200 and 1000 m isobaths. The black diamonds depict the IPC intensification signals whereas the black inverted triangles show the eddy events, all mentioned in the text.</figcaption></figure>

fications along the slope are related to an increase of 0.5–1 °C over time (not shown), which is the increase in temperature that is typically associated with the slope current intensification in the study area (Esnaola et al., 2013). The spatial extension of the warm water masses and the IPC along the French shelf/slope depends on the event, and it coincides with the area where the highest agreement between $AC_G$ and HFR total currents is observed. During the IPC event of November 2009, the warm water tongue is closer to the coast. In this event, the strongest agreement between $AC_G$ and HFR total currents is observed over the slope, while they disagree in the north-western area of the domain. Otherwise, during December 2010, the warm water extends over adjacent offshore areas, as well as the area presenting the strongest agreement between HFR and the altimeter. In three of the four events, $AC_G$ and HFR total currents show a lower fit over the shelf. These observations corroborate the results obtained in Sect. 3.1, where the best statistical results are obtained for the data pairs inside the slope area (Fig. 5).

The relationship between the IPC and the NAO (North Atlantic Oscillation) in the study area was described in Garcia-Soto et al. (2002) and Garcia-Soto (2004). They concluded that for strong IPC years, January water warmings (as a signal of the IPC) were related to negative NAO index values in the previous months (November, December). On the other hand, the eastern Atlantic (EA) is also considered to be a possible factor of the IPC intensification, with positive EA values related to current intensifications.

<table><caption>**Table 2.** NAO and EA indexes in the previous 2 months of the events. $m-2$ means 2 months before the event whereas $m-1$ means 1 month before.</caption><thead><tr><th rowspan="2">Event</th><th colspan="2">NAO</th><th colspan="2">EA</th></tr><tr><th>$m-2$</th><th>$m-1$</th><th>$m-2$</th><th>$m-1$</th></tr></thead><tbody><tr><td>November 2009</td><td>1.62</td><td>-0.61</td><td>0.59</td><td>0.96</td></tr><tr><td>January 2011</td><td>-1.84</td><td>-1.8</td><td>0.24</td><td>-0.49</td></tr><tr><td>January 2014</td><td>0.81</td><td>0.79</td><td>0.09</td><td>1.3</td></tr><tr><td>November 2014</td><td>1.72</td><td>-0.87</td><td>0.2</td><td>1.02</td></tr></tbody></table>

For the four events studied here, the relationship between the IPC intensification and those indexes is shown in Table 2. In general, the NAO (EA) indexes are negatively (positively) related to the IPC in strong intensification periods; however, this relationship does not always apply (see, for instance, the events of 2011 and 2014 where a negative EA index and a positive NAO index are observed, respectively). Moreover, the intensity of the currents is not related to the amplitude of the index, leading to the same conclusion as Le Cann and Serpette (2009) and Le Hènaff et al. (2011).
<figure><img src="image_10.png" /><figcaption>**Figure 7.** Snapshots showing four slope current intensification events observed by HFR, altimetry, and SST (see the dates of the events depicted in Fig. 6) in November 2009 (a), December 2010–January 2011 (b), January 2014 (c), and November 2014 (d). The small arrows depict the HFR current fields (not rotated) whereas the thick ones indicate the across-track $AC_G$. The black line shows the altimeter’s track. The colormap depicts the sea surface temperature (°C) with values referenced in each colour scale. Note that the dates corresponding to the data are not the same for SST or for HFR and $AC_G$ (specified in each panel’s title).</figcaption></figure>

### 3.2.2 Observability of mesoscale eddies

Figure 8 shows four examples of eddies detected by the HFR and the altimeter. Although the effect of the presence of mesoscale eddies has not been explored in terms of statistical results, there is a qualitative agreement between $AC_G$ and HFR total currents when eddies are observed in the area covered by the two measuring systems (even if the eddy core is not crossed by the altimeter track). This happens when either the eddies cross the track of the altimeter or when the size of the eddies, whose centre is located off the track, is large enough to be observable by the altimetry. Across-track $AC_G$ are generally in agreement with the HFR current fields, mainly in terms of current directions. For example, on 15 September 2012 the altimeter crosses a small anticyclone located to the north of the Capbreton canyon head. The maximum $SLA'$ and minimum $AC_G$ near the core of the eddy and the patterns of the $SLA'$ and $AC_G$ north and southwards of the core agree with the structure detected by the HFR. An eddy with a similar diameter and located near the same area, is observed on 15 February 2015. In this case the eddy is cyclonic, and though the HFR and altimetry currents in the area occupied by the structure agree with each other, this is not the case in the rest of the track. During winter, on

15 March 2014, a cyclone is detected by both data sets near the head of the Capbreton canyon, in addition to another cyclone in the north-western part of the domain. Finally, on 7 December 2014, an anticyclone, more energetic than all the former eddies, is observed over the Capbreton canyon. This anticyclone was analysed by Rubio et al. (2018), showing that it had an important role in the offshore transport of coastal waters. North of this eddy, the altimetry and the HFR detect a cyclonic circulation, but in this case, it is not clear from the HFR total current fields that the structure is an eddy.

## 4 Summary and conclusions

In this study, we have investigated the synergies and differences between land-based HFR and satellite altimetry, two remote-sensing techniques that provide measurements of the ocean surface currents at different temporal and spatial scales. A general agreement between HFR and altimetry was observed in the study area, with correlations ranging up to 0.7. The comparisons were carried out in terms of time anomaly of currents, following different approaches with radial and total OMA HFR data. In all cases the addition of
<figure><img src="image_11.png" /><figcaption>**Figure 8.** Four mesoscale eddies observed in the study area. The dots show the points of track 248 of the CMEMS database. SLA' (cm) values are indicated in the colour scale. Black arrows depict the HFR current fields. Red arrows correspond to across-track $AC_G$ derived from the SLA' values in the dots. Grey lines: 200, 1000, and 2000 m isobaths. Note that the scale of each kind of arrow is not the same.</figcaption></figure>

the LF Ekman component ($AC_E$) to the geostrophic component ($AC_G$) improved the results.

The best agreement between both data sets was observed in the slope area, mainly between 200 and 1000 m isobaths, where the surface circulation was dominated by a more energetic geostrophic component. In the coastal area, the agreement between both data sets was lower. $AC_G$ have higher variability than $AC_{HFR}$, which could be explained by the error propagation in the finite-difference operator that estimates the geostrophic velocity from SLA' or by a more effective filtering of HF signals in the radar data (where HF signals can be efficiently removed using a running average temporal filter).

In terms of monthly mean currents, north-eastward currents were observed in all data sets in late autumn and in winter, while weaker north-eastward and south-westward currents were observed in spring and summer. In the winter period, higher variability was also observed at points *E* and *W*, possibly related to a more energetic slope current regime. Additionally, high variability was observed at point *W* in March–April, probably linked to an anticyclonic structure near Torrelavega canyon.

Four IPC events were isolated and described further by means of additional SST data. From this analysis we conclude that during the IPC intensifications the qualitative agreement between $AC_G$ and HFR total currents is great and well related to the SST anomalies. In addition, a relationship between strong IPC events with negative (positive) NAO (EA) indexes in the previous months was detected for three of the four events. On the other hand, although the effect in terms of statistical results of the presence of mesoscale eddies was not explored, there is a qualitative agreement between $AC_G$ and HFR total currents when eddies are observed in the area.

The low correlation between HFR and altimetry observed in some areas and periods can be due to several factors. It is worth noting that both technologies are based on different physical approaches to measure currents, at different spatial and temporal scales, and work under different physical assumptions. Besides, the quality of the radar data is expected to decrease in the furthest points from the antennae and varies as a function of the angle formed by the radial current components used for total current estimations (affecting the along-track comparison). Altimetry also has its own lim-
itations and might have errors in the data editing procedure or in the corrections.

Future work should be oriented towards a better understanding of the relationship of the surface circulation and the dynamics of the subsurface layers by means of the combination of remote observations with data in the water column. Since the comparison near the shoreline is inconsistent, another future work line could be the investigation of the assumptions of geostrophic balance in the coastal area and the merging of altimetry and radar measurements to improve both products. In addition, further comparison with HFR data and higher-resolution coastal altimetry products would enable a better understanding of the differences between both observing systems.

*Data availability.* The CTOH-XTRACK altimetry product is available on the CTOH website (http://ctoh.legosobs-mip.fr/products/coastal-products/coastal-products-1/sla-1hz, last access: 8 October 2018).

The CMEMS along-track level-3 altimetry product is available on the CMEMS website (http://marine.copernicus.eu/services-portfolio/access-to-products/?option=com_csw&view=details&product_id=SEALEVEL_GLOPHY_L3_REPOBSERVATIONS_008_045, last access: 8 October 2018).

The SST data used were retrieved from the NERC Earth Observation Data Acquisition and Analysis Service website (http://www.neodaas.ac.uk/data, last access: 8 October 2018).

The HF radar data are available at http://www.euskoos.eus/en/radar-higer-en/ (last access: 8 October 2018) and at http://www.emodnet-physics.eu/map/platinfo/piradar.aspx?platformid=10274 (last access: 8 October 2018). The OMA HF radar data can be provided upon request to the authors.

*Author contributions.* IMN, AC, AR: contribution to the main structure and contents of all sections; drafting, review and final approval of the submitted version. In addition, IMN produced the figures. CD, FB: contribution to the main structure and contents; review of all sections and advice with regard to altimetry issues.

*Competing interests.* The authors declare that they have no conflict of interest.

*Special issue statement.* This article is part of the special issue "Coastal marine infrastructure in support of monitoring, science, and policy strategies". It is not associated with a conference.

*Acknowledgements.* This study has been supported by the JERICO-NEXT project, funded by the European Union's Horizon 2020 research and innovation programme under grant agreement no. 654410 and the COMBAT project supported by the second call of the Service Evolution of CMEMS. The work of Anna Rubio was partially supported by the LIFE-LEMA

project (LIFE15 ENV/ES/000252). This study has also been undertaken with the financial support of the Department of Environment, Regional Planning, Agriculture and Fisheries of the Basque Government (Marco Program). Ivan Manso-Narvarte was also supported by a PhD fellowship from the Department of Environment, Regional Planning, Agriculture and Fisheries of the Basque Government. The HFR system, whose data have been used herein, is owned by the Directorate of Emergency and Meteorological Services (DAEM) of the Basque Government. Altimetry data used in this study were developed, validated, and distributed partly by the CTOH/LEGOS, France, and partly by Collecte Localisation Satellites (CLS), France, as the Sea Level Thematic Assembly Center of the Copernicus Marine Environment Monitoring Service. The SST level-2 images were produced and distributed by the NERC Earth Observation Data Acquisition and Analysis Service (NEODAAS, http://www.neodaas.ac.uk/data, last access: 1 October 2018). This is contribution number 880 of the Marine Research Division of AZTI-Tecnalia.

Edited by: Stefania Sparnocchia

Reviewed by: two anonymous referees

## References

Ardhuin, F., Marie, L., Rascle, N., Forget, P., and Roland, A.: Observation and estimation of lagrangian, stokes and eulerian currents induced by wind and waves at the sea surface, J. Phys. Oceanogr., 39, 2820–2838, https://doi.org/10.1175/2009JPO4169.1, 2009.

Brown, G. S.: The average impulse response of a rough surface and its applications, IEEE J. Ocean. Eng., 25, 67–74, https://doi.org/10.1109/TAP.1977.1141536, 1977.

Caballero, A., Espino, M., Sagarminaga, Y., Ferrer, L., Uriarte, A., and Gonzalez, M.: Simulating the migration of drifters deployed in the Bay of Biscay, during the Prestige crisis, Mar. Pollut. Bull., 56, 475–482, https://doi.org/10.1016/j.marpolbul.2007.11.005, 2008a.

Caballero, A., Pascual, A., Dibarboure, G., and Espino, M.: Sea level and Eddy Kinetic Energy variability in the Bay of Biscay, inferred from satellite altimeter data, J. Marine Syst., 72, 116–134, 2008b.

Caballero, A., Ferrer, L., Rubio, A., Charria, G., Taylor, B. H., and Grima, N.: Monitoring of a quasi-stationary eddy in the Bay of Biscay by means of satellite, in situ and model results, Deep-Sea Res. Pt. II, 106, 23–37, https://doi.org/10.1016/J.DSR2.2013.09.029, 2014.

Caballero, A., Rubio, A., Ruiz, S., Le Cann, B., Testor, P., Mader, J., and Hernández, C.: South-eastern Bay of Biscay eddy-induced anomalies and their effect on chlorophyll distribution, J. Marine Syst., 162, 57–72, https://doi.org/10.1016/j.jmarsys.2016.04.001, 2016.

Charria, G., Lazure, P., Le Cann, B., Serpette, A., Reverdin, G., Louazel, S., Batifoulier, F., Dumas, F., Pichon, A., and Morel, Y.: Surface layer circulation derived from Lagrangian drifters in the Bay of Biscay, J. Marine Syst., 109–110, S60–S76, 2013.

Chavanne, C. P. and Klein, P.: Can oceanic submesoscale processes be observed with satellite altimetry?, Geophys. Res. Lett., 37, L22602, https://doi.org/10.1029/2010GL045057, 2010.
Cleveland, W. S. and Devlin, S. J.: Locally Weighted Regression: An Approach to Regression Analysis by Local Fitting, J. Am. Stat. Assoc., 83, 403, 1988.

Dussurget, R., Birol, F., Morrow, R., and De Mey, P.: Fine resolution altimetry data for a regional application in the Bay of Biscay, Mar. Geod., 34, 447–476, 2011.

Ekman, V. W.: On the influence of the Earth’s rotation on ocean currents, Ark. Mat. Aston. Fys., 2, 1–53, 1905.

Esnaola, G., Sáenz, J., Zorita, E., Fontán, A., Valencia, V., and Lazure, P.: Daily scale wintertime sea surface temperature and IPC-Navidad variability in the southern Bay of Biscay from 1981 to 2010, Ocean Sci., 9, 655–679, https://doi.org/10.5194/os-9-655-2013, 2013.

Ferrer, L., Liria, P., Bolaños, R., Balseiro, C., Carracedo, P., González-Marco, D., González, M., Fontán, A., Mader, J., and Hernández, C.: Reliability of coupled meteorological and wave models to estimate wave energy resource in the Bay of Biscay, in: Proceedings of the 3rd International Conference on Ocean Energy (ICOE), 6–8 October 2010, Bilbao, Spain, 6 pp., 2010.

Garcia-Soto, C.: “Prestige” oil spill and Navidad flow, J. Mar. Biol. Assoc. UK, 84, 297–300, 2004.

Garcia-Soto, C., Pingree, R. D., and Valdés, L.: Navidad development in the southern Bay of Biscay: Climate change and swoddy structure from remote sensing and in situ measurements, J. Geophys. Res.-Oceans, 107, 3118, https://doi.org/10.1029/2001JC001012, 2002.

González, M., Uriarte, A., Fontán, A., Mader, J., and Gyssels, P.: Marine dynamics, in: Oceanography and Marine Environment of the Basque Country, 1, Elsevier Oceanography Series, edited by: Borja, A. and Collins, M., Amsterdam, the Netherlands, 133–157, 2004.

Graber, H. C., Haus, B. K., Chapman, R. D., and Shay, L. K.: HF radar comparisons with moored estimates of current speed and direction: expected differences and implications, J. Geophys. Res., 102, 18749–18766, https://doi.org/10.1029/97JC01190, 1997.

Gurgel, K.-W.: Shipborne measurement of surface current fields by HF radar, L'Onde – Electr., 74, 54–59, 1994.

Herbert, G., Ayoub, N., Marsaleix, P., and Lyard, F.: Signature of the coastal circulation variability in altimetric data in the southern Bay of Biscay during winter and fall 2004, J. Marine Syst., 88, 139–158, 2011.

Kaplan, D. M. and Lekien, F.: Spatial interpolation and filtering of surface current data based on open-boundary modal analysis, J. Geophys. Res., 112, C12007, https://doi.org/10.1029/2006JC003984, 2007.

Large, W. and Pond, S.: Open ocean momentum flux measurements in moderate to strong winds, J. Phys. Oceanogr., 11, 324–336, 1981.

Law, K.: Measurements of Near Surface Ocean Currents Using HF Radar, PhD thesis, University of California, Santa Cruz, 2001.

Lazure, P.: La circulation des eaux dans le Golfe de Gascogne, in: 10émes encontres interregionales de l'AGLIA, Saint Jean de Luz, 83–88, 1997.

Le Cann, B. and Serpette, A.: Intense warm and saline upper ocean inflow in the southern Bay of Biscay in autumn–winter 2006–2007, Cont. Shelf. Res., 29, 1014–1025, 2009.

Le Hènaff, M., Roblou, L., and Bouffard, J.: Characterizing the Navidad current interannual variability using coastal altimetry, Ocean Dynam., 61, 425–437, 2011.

Le Traon, P. Y., Faugère, Y., Hernandez, F., Dorandeau, J., Mertz, F., and Ablain, M.: Can we merge GEOSAT Follow-On with TOPEX/Poseidon and ERS-2 for an improved description of the ocean circulation?, J. Atmos. Ocean. Tech., 20, 889–895, 2003.

Lipa, B. J. and Barrick, D. E.: Least-Squares Methods for the Extraction of Surface Currents from CODAR Crossed-Loop Data: Application at ARSLOE, IEEE J. Ocean. Eng., 8, 1–28, https://doi.org/10.1109/JOE.1983.1145578, 1983.

Liu, Y., Weisberg, R. H., Vignudelli, S., Roblou, L., and Merz, C. R.: Comparison of the X-TRACK altimetry estimated currents with moored ADCP and HF radar observations on the West Florida Shelf, Adv. Space Res., 50, 1085–1098, https://doi.org/10.1016/j.asr.2011.09.012, 2012.

Liu, Y., Kerkering, H., and Weisberg, R. H. (Eds.): Coastal Ocean Observing Systems, Elsevier, Academic Press, London, UK, 461 pp., 2015.

Morrow, R., Carret, A., Birol, F., Nino, F., Valladeau, G., Boy, F., Bachelier, C., and Zakardjian, B.: Observability of fine-scale ocean dynamics in the northwestern Mediterranean Sea, Ocean Sci., 13, 13–29, https://doi.org/10.5194/os-13-13-2017, 2017.

Pascual, A., Lana, A., Troupin, C., Ruiz, S., Faugère, Y., Escudier, R., and Tintoré, J.: Assessing SARAL/AltiKa Delayed-Time Data in the Coastal Zone: Comparisons with HF Radar Observations, Mar. Geod., 38, 260–276, https://doi.org/10.1080/01490419.2015.1019656, 2015.

Pingree, R. D.: North Atlantic and North Sea Climate Change: curl up, shut down, NAO and Ocean Colour, J. Mar. Biol. Ass. UK, 85, 1301–1315, 2005.

Pingree, R. D. and Le Cann, B.: Three anticyclonic Slope Water Oceanic eDDIES (SWODDIES) in the southern Bay of Biscay in 1990, Deep-Sea Res., 39, 1147–1175, 1992.

Pujol, M.-I., Faugère, Y., Taburet, G., Dupuy, S., Pelloquin, C., Ablain, M., and Picot, N.: DUACS DT2014: the new multi-mission altimeter data set reprocessed over 20 years, Ocean Sci., 12, 1067–1090, https://doi.org/10.5194/os-12-1067-2016, 2016.

Rio, M.-H. and Hernandez, F.: High-frequency response of wind-driven currents measured by drifters and altimetry over the world ocean, J. Geophys. Res., 108, 3283, https://doi.org/10.1029/2002JC001655, 2003.

Rio, M.-H., Mulet, S., and Picot, N.: Beyond GOCE for the ocean circulation estimate: Synergetic use of altimetry, gravimetry, and in situ data provides new insight into geostrophic and Ekman currents, Geophys. Res. Lett., 41, 8918–8925, https://doi.org/10.1002/2014GL061773, 2014.

Roesler, C. J., Emery, W. J., and Kim, S. Y.: Evaluating the use of high-frequency radar coastal currents to correct satellite altimetry, J. Geophys. Res.-Oceans, 118, 3240–3259, https://doi.org/10.1002/jgrc.20220, 2013.

Rubio, A., Blanke, B., Speich, S., Grima, N., and Roy, C.: Mesoscale eddy activity in the southern Benguela upwelling system from satellite altimetry and model data, Prog. Oceanogr., 83, 288–295, 2009.

Rubio, A., Reverdin, G., Fontán, A., González, M., and Mader, J.: Mapping near-inertial variability in the SE Bay of Biscay from HF radar data and two offshore moored buoys, Geophys. Res. Lett., 38, L19607, https://doi.org/10.1029/2011GL048783, 2011.
Rubio, A., Mader, J., Corgnati, L., Mantovani, C., Griffa, A., Novellino, A., Quentin, C., Wyatt, L., Schulz-Stellenfleth, J., Horstmann, J., Lorente, P., Zambianchi, E., Hartnett, M., Fernandes, C., Zervakis, V., Gorringe, P., Melet, A., and Puillat, I.: HF radar activity in European coastal seas: next steps towards a pan-European HF radar network, Front. Mar. Sci., 4, 8, https://doi.org/10.3389/fmars.2017.00008, 2017.

Rubio, A., Caballero, A., Orfila, A., Hernández-Carrasco, I., Ferrer, L., González, M., Solabarrieta, L., and Mader, J.: Eddy-induced cross-shelf export of high Chl-a coastal waters in the SE Bay of Biscay, Remote Sens. Environ., 205, 290–304, 2018.

Schmidt, R. O.: Multiple emitter location and signal parameter estimation, IEEE T. Antenn. Propag., 34, 276–280, 1986.

Solabarrieta, L., Rubio, A., Castanedo, S., Medina, R., Fontán, A., González, M., Fernández, V., Charria, G., and Hernández, C.: Surface water circulation patterns in the southeastern Bay of Biscay: new evidences from HF radar data, Cont. Shelf Res., 74, 60–76, 2014.

Solabarrieta, L., Rubio, A., Cárdenas, M., Castanedo, S., Esnaola, G., Méndez, F. J., Medina, R., and Ferrer, L.: Probabilistic relationships between wind and surface water circulation patterns in the SE Bay of Biscay, Ocean Dynam., 65, 1289–1303, 2015.

Solabarrieta, L., Frolov, S., Cook, M., Paduan, J., Rubio, A., González, M., Mader, J., and Charria, G.: Skill assessment of HF radar-derived products for Lagrangian simulations in the Bay of Biscay, J. Atmos. Ocean. Tech., 33, 2585–2597, https://doi.org/10.1175/jtech-d-16-0045.1, 2016.

Stewart, R. H.: Introduction to Physical Oceanography, University Press of Florida, 2009.

Troupin, C., Pascual, A., Valladeau, G., Pujol, I., Lana, A., Heslop, E., Ruiz, S., Torner, M., Picot, N., and Tintoré, J.: Illustration of the emerging capabilities of SARAL/AltiKa in the coastal zone using a multi-platform approach, Adv. Space Res., 55, 51–59, 2015.

Verron, J., Bonnefond, P., Aouf, L., Birol, F., Bhowmick, S. A., Calmant, S., Conchy, T., Crétaux, J.-F., Dibarbourne, G., Dubey, A. K., Faugère, Y., Guerreiro, K., Gupta, P. K., Hamon, M., Jebri, F., Kumar, R., Morrow, R., Pascual, A., Pujol, M.-I., Rémy, E., Rémy, F., Smith, W. H. F., Tournadre, J., and Vergara, O.: The Benefits of the Ka-Band as Evidenced from the SARAL/AltiKa Altimetric Mission: Scientific Applications, Remote Sens., 10, 163, https://doi.org/10.3390/rs10020163, 2018.