Fine-Scale Ocean Currents Derived From in situ Observations in Anticipation of the Upcoming SWOT Altimetric Mission

1. Introduction

Over the last few decades, remote sensing observations of sea surface height (SSH) have greatly increased our understanding of the global ocean large and mesoscale circulation (e.g., Chelton et al., 2011; Le Traon, 2013). Altimetric observations provide global coverage of the ocean surface on a daily basis. The spatial resolution of along-track altimetric observations is however coarse, with nominal wavelengths ranging from ~40 to ~110 km (Dufau et al., 2016). This results in a mean effective spatial resolution for gridded SSH maps of ~200 km wavelength for the global ocean at mid-latitudes, and ~130 km wavelength for the Mediterranean Sea (Ballarotta et al., 2019). These resolutions are insufficient to capture the entire range of mesoscale dynamics, especially in regional seas such as the Mediterranean, where the Rossby radius of deformation is small (~5–15 km) and mesoscale processes are characterized by shorter spatial scales than in other regions of the ocean (Chelton et al., 1998; Beuvier et al., 2012; Escudier et al., 2016; Barceló-Llull et al., 2019; Kurkin et al., 2020).

The new Surface Water and Ocean Topography (SWOT) satellite mission will be launched in 2022 and it is considered to be the next major breakthrough in satellite ocean observation (Morrow et al., 2019). The SWOT mission aims to provide SSH measurements in two dimensions along a wide-swath altimeter track with an expected effective resolution down to wavelengths of 15–30 km. This will allow, in some regions, observation of the full range of mesoscale dynamics, i.e., scales larger than the first Rossby radius of deformation (Fu and Ferrari, 2008; Fu and Ubelmann, 2014; Wang et al., 2019). During the fast-sampling phase after launch, SWOT will provide observations of SSH in a 1-day-repeat orbit in specific areas of the world ocean for instrumental calibration and validation (d'Ovidio et al., 2019). The region around the Balearic Islands (western Mediterranean Sea) is one of the selected areas for the SWOT fast-sampling phase, and will be the target of several in situ sampling campaigns aimed at validation of SWOT measurements, and evaluation of the subsurface biophysical activity and its interaction with the surface dynamics resolved by SWOT.

In preparation for these experiments, a multi-platform experiment called PRE-SWOT was conducted south of the Balearic Islands in 2018. The PRE-SWOT cruise objective was to collect in situ data from different observational platforms in order to explore the three-dimensional (3D) circulation at scales of 20 km wavelength that SWOT aims to resolve (Morrow et al., 2019) and evaluate the variability not captured by the current altimetric constellation. The scales of 20 km wavelength are within the range of fine-scales (1–100 km) as defined by the altimetric community (d'Ovidio et al., 2019; Morrow et al., 2019). The study of the oceanic features associated to these scales (such as fronts, meanders, eddies, and filaments) is important because they play a critical role in the distribution of heat, salt, gases, carbon, and nutrients in the global ocean (e.g., Lévy et al., 2001; Thomas, 2008; Mahadevan, 2016; Klein et al., 2019). Understanding the 3D dynamics of these features and their impact on the large scale ocean circulation and climate system is one of the major challenges for the next decade in physical oceanography (e.g., Young and Sikora, 2003; Kwon et al., 2010; Ma et al., 2016; Su et al., 2018; Bishop et al., 2020; Small et al., 2020). Integrated approaches, combining multi-platform in situ data, remote sensing observations and numerical modeling, constitute an innovative methodology for the evaluation and understanding of the 3D pathways associated with these structures. Some of the multi-platform experiments that have recently followed this approach are LatMix (Shcherbina et al., 2015), AlborEx (Pascual et al., 2017; Ruiz et al., 2019), LASER (D'Asaro et al., 2020), and CALYPSO (Mahadevan et al., 2020; Freilich and Mahadevan, 2021; Tarry et al., 2021).

One of the motivations to study the scales that SWOT aims to resolve is their importance in terms of coupled biophysical processes (Lévy et al., 2018). Vertical motions associated with these features drive the vertical exchange of tracers, including nutrients and passive marine organisms, between the surface layers and the abyss. However, direct measurements of vertical velocities (w) are difficult to obtain due to their small magnitudes when compared to horizontal velocities (u): ${O}$(w) ~ (10⁻³-${10}^{-4})·{O}$(u) ~ 10 m day⁻¹ for the mesoscale and ${O}$(w) ~ (${10}^{-2})·{O}$(u) ~ 100 m day⁻¹ for the submesoscale (Mahadevan and Tandon, 2006; Baschek and Farmer, 2010; D'Asaro et al., 2011, 2018). In consequence, indirect approaches have been applied to calculate the vertical velocity field from observational data using distinct forms of the so-called omega equation or inverse methods (e.g., Hoskins et al., 1978; Viúdez et al., 1996; Pallàs-Sanz and Viúdez, 2005; Thomas et al., 2010; Barceló-Llull et al., 2017). For flows with low Rossby numbers, the quasi-geostrophic (QG) omega equation has been widely used to compute the vertical velocity field from observations of temperature and salinity (e.g., Tintoré et al., 1991; Pollard and Regier, 1992; Allen and Smeed, 1996; Pascual et al., 2015; Barceló-Llull et al., 2016; Ruiz et al., 2019; Buongiorno Nardelli, 2020). Pietri et al. (2021) have recently analyzed the skills and limitations of the QG omega equation in several regions of the ocean with different dynamics. Using a fully eddy-resolving numerical simulation with a 1/16° horizontal resolution, they demonstrate that the QG omega equation provides satisfactory results for scales larger than ~10 km.

With the objective of improving the accuracy and reliability of altimetric observations, the European Union Copernicus Programme launched the Sentinel-3A satellite mission in 2016. The Sentinel-3A synthetic aperture radar altimeter (SRAL) is based on the synthetic aperture radar mode (SARM) principle proposed by Raney (1998). It has been demonstrated that, compared to conventional altimetry (or low-resolution mode), this technique significantly reduces noise levels in the measurements and increases the along-track spatial resolution (Boy et al., 2017; Heslop et al., 2017; Sánchez-Román et al., 2020). To focus on the analysis of (i) the fine-scale sea surface circulation over the region of study and (ii) the limitations of present-day nadir altimeters, we additionally compare observations from Sentinel-3A—together with in situ observations from an ocean glider and a ship-based hydrographic survey—with gridded altimetric data. This comparison complements previous studies reported in the area by using a multi-platform approach (e.g., Ruiz et al., 2009; Bouffard et al., 2010; Cotroneo et al., 2016; Heslop et al., 2017; Aulicino et al., 2018).

The objective of this study is to compare observations from present-day altimeters with measurements from different in situ instruments to explore the horizontal circulation at scales of 20 km that SWOT aims to resolve and evaluate the limitations of the current altimetric constellation. We also explore the 3D QG vertical velocity field associated with the fine-scale dynamics of the region of study. In addition, we analyse the relationship between the vertical velocity field and nutrient concentrations measured from water samples, and the chlorophyll-a signal observed from an ocean glider.

2. Data and Methods

2.1. The PRE-SWOT Cruise

The PRE-SWOT cruise was conducted aboard the R/V García del Cid between 5 and 17 May 2018 south of the Balearic Islands (western Mediterranean Sea). Observations from in situ platforms, including two gliders, twelve drifters, a Conductivity Temperature Depth (CTD) probe attached to a rosette system, a hull-mounted 75 kHz RDI Acoustic Doppler Current Profiler (ADCP), and water samples taken at different depths, were collected together with satellite data.

2.2. Rosette CTD Casts and Water Samples

The field experiment was divided into two legs. Leg 1 took place between 6 and 10 May 2018 (4.5 days) with CTD sampling over a 3D high-resolution regular grid in order to resolve the fine-scales expected from SWOT (Gómez-Navarro et al., 2020). During Leg 1 we made 63 rosette CTD casts at intervals of 10 km, resulting in a regular grid of 60 × 80 km (Figure 1B). The maximum depth of most of the CTD casts was 500 m (pink dots in Figure 1B), with 3 casts down to 1,000 m (blue dots in Figure 1B). The rosette carried an SBE 911plus CTD instrument that was additionally equipped with SeaPoint Fluorometer and Turbidity sensors, a SBE43 dissolved oxygen sensor, a WetLabs Transmissometer (CST-974DR), and a Biospherical irradiance radiometer. These were attached to the rosette system of 12 oceanographic 12-liter Niskin bottles. Water samples were collected with the rosette at different depths in some stations in order to calibrate salinity, dissolved oxygen and in vivo fluorescence (yellow stars in Figure 1B). Additionally, water samples for nutrients, pigment signatures, and phytoplankton were taken in selected stations, chosen from remote sensing ocean color maps in order to support an analysis of the main observed features (green dots in Figure 1B).

Leg 2 took place between 13 and 16 May 2018 and focused on completing four transects parallel to Sentinel-3A track 244 that crossed the study region on 13 May 2018 (Figure 1C). For each transect the water column was sampled from the surface to 500 m depth with 10 CTD casts 10 km apart; the transects were separated by 10 km. Water samples for nutrients, pigment signatures and phytoplankton determination were collected in all the rosette casts of the three transects closest to Sentinel-3A track 244 (green dots in Figure 1C). Water samples for calibration of sensors were collected in one cast out of three along these three transects (yellow stars in Figure 1C), while along the western transect only water samples for conductivity calibration were collected (blue stars in Figure 1C). A detailed description of the CTD data processing and calibration can be found in the PRE-SWOT cruise report (Barceló-Llull et al., 2018). We used the Thermodynamic Equations of Seawater (TEOS-10) functions to calculate potential temperature from in situ temperature, and the potential density anomaly (McDougall and Barker, 2011).

2.3. ADCP

Current velocities were continuously recorded using a hull-mounted 75 kHz RDI ADCP at a transit speed of ~8 knots during the two legs of the experiment. The ADCP provided raw data from the surface to ~800 m in bin sizes of 8 and 16 m. The raw data were quality controlled, corrected for heading misalignment and edited with the Common Oceanographic Data Access System (CODAS, Firing et al., 1995). Further details of the ADCP data processing can be found in the PRE-SWOT cruise report (Barceló-Llull et al., 2018).

2.4. Drifters

During Leg 1, twelve surface velocity program (SVP) drifters with a drogue centered at 15 m depth were deployed within the region of study. The aim was to sample frontal areas and study relative dispersion by horizontal stirring at the mesoscale and submesoscale. For that purpose, the drifters were arranged into triangles with a separation between drifters of 5 km (red dots and triangles in Figure 1B). Drifter deployment locations were chosen based on fronts that were visible in ocean color imagery. To remove the inertial signals, the drifter trajectories were filtered using a fifth-order Butterworth low-pass filter with a cutoff of 1.5 times the inertial period of the area, which corresponds to 19 hours (Essink et al., 2019). The zonal and meridional velocity components of each drifter were computed by finite differencing their trajectories.

2.5. Gliders

As part of the PRE-SWOT multi-platform experiment, two underwater gliders followed the Sentinel-3A track 244 (red line in Figures 1A,C): one Slocum from the Balearic Islands Coastal Observing and Forecasting System (SOCIB) and one SeaExplorer from the Mediterranean Institute of Oceanography (MIO). Both gliders were equipped with a pumped CTD sensor (Seabird's GPCTD). The Slocum glider was also equipped with an Anderaa optode for the measurements of dissolved oxygen concentration (O2) and a WetLabs ECO FLNTU for the measurements of chlorophyll-a (chl-a) fluorescence (λEx/λEm: 470/695 nm) and turbidity at 700 nm. The SeaExplorer was equipped with an O2 sensor (Seabird's SBE-43F), from which concentrations are derived, and with a WetLabs ECO FLBBCD for the measurements of chl-a fluorescence (λEx/λEm: 470/695 nm), backscattering at 700 nm (BB700) and humic-like fluorophore fluorescence (λEx/λEm: 370/460 nm). The manufacturer's calibrations were applied to all these sensors. Additionally, the SeaExplorer was equipped with two MiniFluo sensors capable of measuring fluorescent dissolved organic matter (DOM), including organic pollutants such as Polycyclic Aromatic Hydrocarbons (PAHs).

On 3 May 2018, 2 days before the PRE-SWOT cruise started, the Slocum and the SeaExplorer gliders were deployed south of Mallorca. The glider mission plan involved 2 phases. In the first phase, the gliders transited during 2 days between the launch waypoint and the survey start-waypoint (38.94°N, 2.97°E); this phase also served as a test and validation. During the second phase, that extended between 5 and 15 May 2018, the SOCIB glider followed the Sentinel-3A track 244 as closely as possible, whilst the MIO glider sailed 5 km east of the satellite track, in parallel with the SOCIB glider. Over 10 days, both gliders made a round-trip over their respective 100-km long tracks, and finished their missions back at the survey start-waypoint. Glider profiles sampled the water column from the sea surface to a maximum depth of 650 m for both units and were collected at a spatial resolution of about 8 km. Raw glider data from the SOCIB unit were processed with the SOCIB Glider Toolbox, which includes a thermal lag correction (Garau et al., 2011; Troupin et al., 2015). Raw SeaExplorer data were also processed with the SOCIB Glider Toolbox, but at that moment no thermal lag correction was available for the SeaExplorer. We used temperature and salinity observations from the Slocum glider to calculate the dynamic height and compare the scales resolved by the glider with the scales resolved by present altimetric products and CTD observations. To evaluate the impact of vertical velocities on the vertical distribution of tracers, we use chl-a concentrations derived from the SeaExplorer fluorescence measurements and equivalent of a mono-culture of phytoplankton Thalassiosira weissflogii (manufacturer's calibration, https://www.whoi.edu/fileserver.do?id=199125&pt=2&p=207009).

2.6. Satellite Data

The choice of location for sampling of oceanic features was guided and updated in real time by remote sensing imagery during the cruise. Ocean color (OC), SSH, and derived geostrophic currents from remote sensing were essential components of the design of the sampling strategy. The OC data were obtained from the NASA Ocean Color Level 1&2 Browser (https://oceancolor.gsfc.nasa.gov). Level-2 data acquired by both the Moderate Resolution Imaging Spectroradiometer (MODIS) sensor aboard the Aqua satellite, and the Visible and Infrared Imager/Radiometer Suite (VIIRS) sensor carried by the Suomi-NPP satellite were used to design the Leg 1 CTD sampling strategy and the drifter release location. MODIS Level-2 data have a spatial resolution of 1 km, while VIIRS Level-2 data have a spatial resolution of 750 m.

Before and during the cruise, SSH and derived geostrophic currents were obtained from the near-real-time gridded L4 multi-mission altimeter product for the Mediterranean Sea that is distributed by the Copernicus Marine Environment Monitoring Service (CMEMS). After the cruise ended, the reprocessed product had become available. Absolute Dynamic Topography (ADT) was also used to assess the oceanographic context of the experiment. The Mediterranean gridded Level-4 (L4) altimetric products have daily time resolution and 0.125° × 0.125° grid resolution. ADT is computed by adding the mean dynamic topography (MDT) from Rio et al. (2014) to the satellite sea level anomaly (SLA) observations.

We also used 1 Hz along-track level-3 (L3) ADT data collected by the Sentinel-3A satellite mission along the ground-track 244 (cycle 31). This cycle is associated with the passover of the satellite over the Algerian Basin on 13 May 2018 (red line in Figure 1A). We use the reprocessed product for the European Seas distributed by CMEMS with a spatial resolution of 7 km. Cycle 31 is close in time with the glider mission described above: the glider return transects along the Sentinel-3A ground-track were conducted between 9 and 14 May 2018, whilst the satellite overpassed the region on 13 May 2018. Thus, there is a temporal lag < 5 days between the first glider observations on the return transects along the Sentinel-3A ground-track and the satellite overpass. Note that this temporal lag is lower than the temporal lag of 8–11 days that Aulicino et al. (2018) had between observations in their inter-comparison of four SARAL-AltiKa cycles with glider transects in the same region of study. Only data captured between 37.00°N and 39.05°N along the ground-track were extracted in order to achieve close co-location with the glider measurements; this also reduces contamination of the altimetry data due to proximity to the coast (Aulicino et al., 2018). At the time of the satellite passover, the Slocum glider was located at 38.63°N, 2.80°E (cyan star in Figure 1A).

2.7. Optimal Interpolation

Potential temperature, practical salinity, potential density anomaly, and dissolved oxygen concentration (hereafter temperature, salinity, density, and oxygen, respectively) data from Leg 1 were interpolated onto a regular grid with a vertical resolution of 5 m and a horizontal resolution of 2 km. Processed ADCP velocity was vertically smoothed using a Loess filter with a half-power filter cutoff of 60 m and interpolated onto the same regular grid. Optimal interpolation was used for the horizontal interpolation (Bretherton et al., 1976). The data covariance was calculated using a 2D Gaussian function with semimajor and semiminor axes of Lx = Ly = 20 km (see section 3.1.5 for a comparison with the empirical correlation calculated from drifter and CTD observations). With this interpolation, wavelengths smaller than the covariance function correlation lengths (Lx and Ly) are filtered (Pallàs-Sanz et al., 2010). The mean fields are computed from a planar fit for temperature, salinity, density, and oxygen, and are assumed to be constant for ADCP velocity (Rudnick, 1996). The uncorrelated noise for the interpolation of the ADCP velocity is assumed to be 18% of the signal variance; this is estimated as the ratio of noise-to-signal variance considering an instrumental error with a standard deviation of 0.01 m s⁻¹ (Allen et al., 1996; Gomis et al., 2001). The uncorrelated noise for the interpolation of CTD observations is set to 3% of the signal variance (Barceló-Llull et al., 2017); this value is based on instrumental error and the variance of the field. Other authors have used a range of noise-to-signal ratios from 0.01% (Ruiz et al., 2019) to 5% (Rudnick, 1996), in our case we have tested this ratio for different instrumental errors and found that 3% ensures the filtering of noise due to instrumental errors. Note that the observations reconstructed through optimal interpolation are collected during 4.5 days and with this interpolation we assume that they are quasi-synoptic, i.e., a stationary representation of the ocean at the scales resolved of 20 km (e.g., Rudnick, 1996; Pascual et al., 2004; Barceló-Llull et al., 2017; Ruiz et al., 2019).

2.8. Dynamic Height and Geostrophic Velocity Calculation

Dynamic height (DH) and geostrophic velocity fields were computed from the optimally interpolated density data from Leg 1, assuming a reference level of 1,000 m depth (1,010 dbar), the maximum depth of the CTD measurements. DH was also calculated using the Slocum glider transect data. Temperature and salinity measurements collected during the return transect (between 9 and 14 May 2018) were used to compute DH at 30 m depth, as the glider subsurface inflection depth was 15 m and the spatial resolution of glider profiles was low in the shallower layers (Heslop et al., 2017). The reference level for the calculation was set to 650 dbar (maximum depth of the glider measurements). To compare glider and CTD observations from Leg 1, we also computed DH from the optimally interpolated density data, with an assumed reference level of 650 dbar.

Leg 2 CTD data obtained along the four transects parallel to the glider tracks (Figure 1C) were used to compute DH at 30 m depth. In Leg 2, CTD casts sampled the water column down to 500 m depth, and the reference level for the inference of the DH was set to 500 dbar. The DH calculated along the four CTD transects was compared to that computed from the Slocum glider during the return transect, also referenced to 500 dbar for consistency.

The DH calculated from glider data is compared to the ADT obtained from both the Sentinel-3A dataset and the L4 altimetric gridded product. Prior to the comparison, a low-pass Loess filter (Cleveland and Devlin, 1988) with a cutoff of 30 km was applied to the glider, CTD and Sentinel-3A datasets in order to remove measurement noise and small-scale variability (see e.g. Heslop et al., 2017). The efficiency of this approach for the comparison of DH and ADT, and the derived geostrophic velocities, has been demonstrated in previous studies in the same area (Heslop et al., 2017; Aulicino et al., 2018).

2.9. Finite-Size Lyapunov Exponents

Horizontal advection was analyzed using Lagrangian methods, comparing trajectories of in situ and numerical drifters. The trajectories of the numerical drifters were obtained by integrating the horizontal velocity fields derived from the altimetry, CTD, and ADCP observations with a fourth-order Runge-Kutta integrator and a time step of 3 h, assuming time varying altimetric currents and stationary CTD-derived and ADCP velocity fields (d'Ovidio et al., 2013) during the integration period of 6 days. The period of the numerical integration corresponds to the time needed to reach a predefined final separation between numerical drifters. Here we integrated the numerical drifters until they were separated by 45 km, which corresponds to the final separation between in situ drifters after the drifter distance increased exponentially (red lines in Figures 8E,F). The velocity fields were interpolated linearly in space (and in time for the altimetry-derived currents) to the numerical drifter position. The trajectories were used to compute finite-size Lyapunov exponents (FSLEs) (Aurell et al., 1997; Lacorata et al., 2001; d'Ovidio et al., 2004). For clarity, we recall here some basic concepts related to this technique. An FSLE is estimated by finding the time (τ) needed for initially close particles separated by an initial distance d₀ at time t₀ to reach a prescribed final separation d_f. The equation that defines an FSLE is the following:

An FSLE, therefore, measures the exponential rate of separation among particle trajectories. It can be computed forward or backward in time. The forward calculation studies the dynamics that drifters initialized nearby undergo in terms of their relative distance. On the other hand, the backward separation is used to estimate regions where water masses far away in the past are brought into contact by the circulation. The points of the line along which these confluence dynamics occur all have relatively high values in maps of FSLEs, and they form what is called Lagrangian fronts (Prants et al., 2014) and belong to the class of Lagrangian Coherent Structures (Haller and Yuan, 2000). Water parcels and arrays of drifters are stretched along Lagrangian fronts, which act as barriers to transport. More details on Lyapunov computation and its meaning can be found in d'Ovidio et al. (2004) and Lehahn et al. (2018).

Backward FSLEs were computed from the numerical trajectories in order to identify Lagrangian fronts (and hence transport barriers) present in the velocity fields derived from altimetry, CTD, and ADCP observations. This was done in order to test whether the trajectories of drifters released during the cruise were aligned along these fronts. For this application, the prescribed initial and final separation distances were optimized for better visualization of the fronts with values of 1 and 20 km, respectively, which are typical values used for this type of application (d'Ovidio et al., 2004; Hernández-Carrasco et al., 2011). This calculation was repeated at different depths (5, 50, 100, and 200 m) using the CTD-derived velocity fields to evaluate the vertical extension of the Lagrangian fronts.

Forward FSLEs were computed in order to compare the exponential rate of separation of numerical and in situ drifters. For a reliable comparison, the initial and final separations used to compute the Lyapunov exponents forward in time from the velocity fields are defined as 7 and 45 km, respectively, which correspond to the separations of in situ drifters for the period in which their relative separation is in an exponential regime (see Figures 8E,F).

2.10. Vertical Velocity Estimation

The quasi-geostrophic (QG) vertical velocity field was calculated through the QG omega equation (Hoskins et al., 1978; Tintoré et al., 1991; Buongiorno Nardelli et al., 2012; Pascual et al., 2015) using the optimally interpolated density field and the derived geostrophic velocity from Leg 1. The QG omega equation is valid for Rossby numbers < < 1 and is presented here in its Q vector form:

where

w is the vertical velocity, u_geo is the geostrophic velocity vector, ρ is density, ρ₀ is the mean density, N² is the Brunt–Väisälä frequency ($N^2=-\frac{g}{{\rho }_0}\frac{\partial \rho }{\partial z}$), f is the Coriolis parameter, considered constant and computed at the mean latitude, and g is gravity. In this implementation, N² only depends on depth and is estimated as the horizontal average of the Brunt–Väisälä frequency estimated at each grid point. The Q vector represents the deformation of the horizontal density gradient by the geostrophic velocity field. The forcing term of the QG omega equation is on the right-hand side of (2), while on the left-hand side an elliptic operator is applied to the vertical velocity. The QG omega equation is solved by applying an iterative relaxation method with Dirichlet boundary conditions (w = 0).

2.11. Oceanographic Context From Satellites and Sampling Region Determination

During the month before the PRE-SWOT experiment, ocean color images provided essential information at high resolution for the design of the experiment. However, observations without cloud-coverage in the area of study were limited to some specific days. Figures 2A–C show three examples of the best observations of ocean color preceding the PRE-SWOT experiment. On 17 April 2018, relatively high values of chl-a of the order of 0.3 mg m⁻³ were observed in an elongated meander at ~2.50°E, 38.00°N extending to the northeast (Figure 2A). Three days after, the chl-a signal was concentrated at ~3.20°E, 38.60°N and had values similar to those observed in the region southeast of Mallorca and south of Menorca (Figure 2B). Two days before the cruise experiment, this signal was displaced to ~2.90°E, 38.50°N; this map was used to determine the position of the CTD casts in Leg 1 and the drifter release locations (Figure 2C).

According to altimetric observations, the target sampling area (boxes in Figure 2C) was surrounded by three oceanographic features: a cyclonic eddy in the Algerian basin, a cyclonic eddy south of Menorca, and an anticyclonic eddy east of Ibiza (Figures 2D–F). The western part of the region with the chl-a signal observed at ~2.90°E, 38.50°N on 3 May 2018 (Figure 2C) was characterized by southwestern velocities of the order of 0.20 m s⁻¹ during the month before the cruise experiment. The persistence of the meandering flow and the presence of an associated front visible on ocean color maps determined the region selected for the cruise experiment within the domain that will be sampled by SWOT during the fast-sampling phase. Note that ocean color maps (Figures 2A–C) reveal the presence of filaments, eddies and fronts that are not detected in present-day altimetric observations (Figures 2D–F). The in situ sampling is then expected to provide information at high resolution which is out of reach of present-day altimeters but potentially detectable by future SWOT observations.

3. Results

3.1. Hydrodynamic Signature of a Meandering Front

3.1.1. Hydrological Characterization

Maps of salinity at different depths from the Leg 1 CTD sampling reveal the presence of a meandering front in the upper layers that separates saltier water at the northeastern edge of the domain from fresher water at the southeastern edge (Figures 3A,E). The signal of the front in salinity is apparent from the surface to ~150 m depth (not shown), having a difference of 0.8 between maximum and minimum salinity values at 5 m, the shallowest layer of CTD observations, that decreases with depth. The map of temperature at 5 m depth (Figure 3B) is characterized by a meridional distribution with warmer water at the western edge and colder water at the eastern edge of the sampled domain. Below 5 m and until ~150 m, the temperature signal resembles the front with colder water in the northeastern part of the domain and warmer water in the southeastern part (Figure 3F). The front detected between 5 and ~150 m depth is characterized by saltier and colder water on the northeastern side that is the signal of surface local Atlantic Water (local AW)—AW already modified by a long residence in the Mediterranean. On the southeastern side the front has fresher and warmer water that is the signature of recent Atlantic Water (recent AW) coming from the Strait of Gibraltar (Barceló-Llull et al., 2019). In this depth layer, density is dominated by salinity over thermal effects (Figures 3C,G). At 5 m depth, dissolved oxygen is maximum at the saltier side of the front, with concentration values of 7.54 mg l⁻¹ (Figure 3D). West of 3.00°E the concentration is minimum reaching values of 7.29 mg l⁻¹. The dissolved oxygen distribution at 50 m depth does not seem related to the front (Figure 3H).

3.1.2. Water Masses

To analyse the different water masses residing in the region of study, temperature-salinity diagrams of all the CTD profiles from Leg 1 are shown in Figure 4. A local maximum of temperature and salinity between 200 and 400 m depth reveals the presence of Levantine Intermediate water (LIW), and temperature values between 12.7°C and 12.9°C are the signal of Western Mediterranean Deep water (WMDW) (Pinot et al., 2002; Balbín et al., 2012; Barceló-Llull et al., 2019). Western Mediterranean Intermediate water (WMIW), associated with a relative minimum of temperature at intermediate layers, is not observed in the sampled region. At ~50 m depth some profiles of saltier and colder local AW have higher concentrations of oxygen than the other profiles (Figure 4A), which may be the signal of a local process of oxygenation. A maximum of dissolved oxygen is normally formed above the deep chlorophyll maximum as a result of photosynthetic activity by phytoplankton cells. This has been observed in the northwestern Mediterranean by Estrada et al. (1999), Segura-Noguera et al. (2016), and Lefevre (2019). Between 300 and 500 m depth, higher temperature and salinity values are associated with lower oxygen concentration (Figure 4B). At depths higher than 500 m, the oxygen concentration is almost constant.

3.1.3. Geostrophic and Ageostrophic Horizontal Velocities

Geostrophic velocities calculated from CTD observations (u_geo) at 20 m depth reveal a meandering flow coming from the center north of the domain and changing its direction to the east with a maximum speed of the order of 0.20 m s⁻¹ (Figure 5A, black arrows). The southeastern part of the domain is characterized by a northeastward flow reaching a maximum speed of 0.26 m s⁻¹, whilst the western part shows an anticyclonic circulation with maximum speeds of the order of 0.15 m s⁻¹. The ADCP velocity (u_adcp), that includes the geostrophic and ageostrophic components of the flow, at 20 m depth shows a similar circulation with differences in direction in the southeastern and northwestern regions of the domain (Figure 5A, cyan arrows). The geostrophic currents derived from altimetry (u_alt; Figure 5B) do not resolve the features that the u_geo and u_adcp maps show. This is due to the lower effective spatial resolution of altimetric gridded products in the Mediterranean (~130 km wavelength, Ballarotta et al., 2019) in comparison with CTD and ADCP observations (~20 km wavelength) that limits the representation of the circulation in the region. Even the general circulation pattern shown in the altimetric map does not correspond to the higher-resolution in situ velocity fields. Only in the western side of the domain, where a southwestward flow dominates the circulation, u_alt have similarities with u_geo and u_adcp. A visual comparison with the trajectories of the drifters released during the cruise also highlights the limitations of present-day altimetry (see drifter trajectories in Figure 5). The geostrophic velocity calculated from CTD data and the total ADCP horizontal velocity show good agreement with the trajectories followed by the drifters, while altimetric currents cannot explain their eastward displacement.

To compare quantitatively the drifter movement with u_alt, u_geo and u_adcp, we have estimated the meridional and zonal velocity components of each drifter from their trajectories between 5 and 12 May 2018. Table 1 shows the Pearson correlation coefficients obtained between each component of the drifter velocity and the corresponding component of the altimetric, geostrophic and ADCP velocities. The zonal component (u) of the drifter velocity has similar correlation coefficients with the three velocity fields and is maximum with the geostrophic velocity (0.90). The meridional component (v) has a correlation coefficient of 0.59 with the altimetric meridional component that increases with the ADCP (0.70) and geostrophic (0.75) velocities. Hence, among the three datasets, the geostrophic velocity field calculated from CTD data is the one showing the best agreement with the trajectories followed by drifters. This result suggests that the strongest component of the horizontal currents responsible for the drifter displacement is geostrophic, and hence potentially resolvable in higher-resolution altimetric maps like the ones that SWOT should provide.

Regarding the vertical distribution of horizontal velocities, zonal sections of the meridional component of the ADCP velocity (v_adcp) show good correspondence with the meridional component of the geostrophic velocity (v_geo) calculated from CTD observations (Figure 6). However, ADCP velocity has higher magnitude and vertical variability than the geostrophic velocity field. Vertical profiles of the root mean square difference (RMSD) show that the differences in magnitude are higher in the upper layers and decrease with depth, in correspondence with the decrease in magnitude of both velocities. The correlation coefficient vertical profiles reveal that both fields have maximum correlation in the upper ~150 m depth, where the front is located. The correlation decreases with depth, where the vertical variability of the ADCP velocity is more evident in comparison with the smoother geostrophic velocity. Dashed (dotted) lines in Figure 6 represent the RMSD and correlation coefficient vertical profiles between the meridional (zonal) components of both velocity fields considering all data in each depth layer. The zonal components have higher RMSD and lower correlation coefficient than the meridional components of the velocities. This translates in higher differences between the zonal velocities along meridional sections (not shown) than between the meridional velocities along the zonal transects shown in Figure 6. Hence, the effect of having more synoptic observations along meridional transects than along zonal transects does not imply smaller differences between ADCP and geostrophic velocity fields. This is due to the optimal interpolation used to reconstruct the observations, which considers all data in each depth layer. Note that these velocities are obtained from two independent instruments and following different procedures, hence, different sources of error may contribute to the small differences observed between both velocity fields. First, the geostrophic velocity calculated from CTD data may be underestimated as it is calculated assuming a level of no motion of 1,000 m, while the ADCP velocity does not have this constraint and also includes the barotropic component of the horizontal velocity, which may increase its magnitude. In addition, from 500 to 1,000 m depth the number of CTD casts is reduced and this could imply an additional smoothing to the interpolated fields. Also, the ADCP velocity field has smaller features on the southeastern and northwestern corners of the domain (Figure 5A) that could be due to (i) the higher resolution along meridional transects of ADCP data compared to the lower resolution of CTD stations, and (ii) potential artifacts of the optimal interpolation close to the boundaries that may affect the estimation of the derivatives involved in the computation of u_geo. Lastly, the ADCP velocity has more vertical variability (Figure 6) that could be the signal of errors coming from the instrument. In addition to the different sources of error, the ADCP velocity may include ageostrophic motions that are not represented in the geostrophic velocity field, such as the cyclostrophic component that could have a contribution in areas of strong curvature, internal waves and inertial motions (Pallàs-Sanz and Viúdez, 2005; Morrow et al., 2019). The small differences observed between ADCP and CTD-derived geostrophic velocities, with maximum RMSD of the order of 0.06 m s⁻¹ and correlation coefficients of ~0.9 in the upper layers, and the different sources of error that may introduce these differences, suggest that the geostrophic component of the horizontal flow may dominate the local dynamics of this region down to scales of 20 km.

The vertical relative vorticity field scaled by the planetary vorticity, or Rossby number, has been estimated from u_geo (Ro^g) and from u_adcp (Ro). Both fields have similar distributions and intensities with small differences due to the deviations observed between velocity fields (Figure 7). Good correspondence between both fields is observed below 20 m depth (not shown), with a higher decrease with depth of the Ro^g magnitude due to the smoother u_geo field (Figure 6). Ro has values ranging from −0.26 to 0.20 and Ro^g between −0.24 and 0.32 in the region of study and at the scales resolved of ~20 km, in accordance with a dominance of the geostrophic component of the flow in this region for the period analyzed. Note that with higher-resolution observations we would sample smaller-scale features and the Rossby number would potentially be higher. Hence, all our estimates are based on the scales we resolve of ~20 km.

3.1.4. Finite-Size Lyapunov Exponents

The backward Lyapunov exponents computed from u_alt, u_geo and u_adcp are shown in Figures 8A–C. Over the Lagrangian fronts we plot the trajectories of two drifters launched in situ on 8 May 2018. After their release, the two drifters moved away from their deployment locations, marking an increasing separation with time (Figure 8E). The behavior of the two drifters particularly underlines the classical presence of a so-called hyperbolic point in the velocity field (e.g., see Figure 3 in Lehahn et al., 2007): the drifters separate along an imaginary line, sometimes referred to as a Lagrangian front or as a Lagrangian Coherent Structure (Haller and Yuan, 2000; Prants et al., 2014). The FSLE maps estimated from u_geo and u_adcp (Figures 8B,C) confirm this suggestion (d'Ovidio et al., 2004), while the Lagrangian fronts derived from u_alt (Figure 8A) are inconsistent with the trajectories of the two drifters, even in qualitative terms. The eastern drifter in particular crosses an altimetry-derived Lagrangian front, which is supposed to behave as a transport barrier. The situation is different for the case of the CTD- and ADCP-derived Lyapunov exponents, which display a better agreement between Lagrangian fronts and drifter trajectories. Indeed, the trajectory of the two drifters is consistent with an expected separation along the attracting front underlined by CTD- and ADCP-derived Lyapunov exponents (Figures 8B,C). The Lagrangian front detected in the CTD-derived Lyapunov exponents appears to extend vertically and vanish approximately at 200 m depth (Figure 8D).

In order to compare the exponential rate of separation computed from the in situ drifter trajectories and the numerical drifters, we proceeded in the following steps. We firstly identified a temporal window over which the relative separation between the two drifters is approximately exponential. This temporal window extends from 10 to 50 h since 8 May 2018 at 10:00:00 (Figures 8E,F). We then used the initial and final separations corresponding to the extrema of this temporal window (7 and 45 km, respectively) to recompute the finite-size Lyapunov exponents from the velocity fields by simulating the trajectories of numerical drifters. The corresponding u_alt, u_geo and u_adcp finite-size Lyapunov exponents, recomputed in this area over the Lagrangian fronts, are shown in Table 2. The exponential rate of separation for the in situ drifters is 1.097 days⁻¹. Altimetry again has the largest mismatch with respect to the in situ drifters, grossly underestimating their rate of separation by more than 50%. CTD- and ADCP-derived exponential rates of separations are closer to the one of the in situ drifters.

3.1.5. Vertical Velocities and Impact on Biochemical Variability

The horizontal distribution of the QG vertical velocity field (w_QG) is similar throughout the upper 300 m depth and is characterized by an upwelling cell at the southern edge of the domain (3.05°E, 38.42°N) surrounded by two downwelling cells at 3.28°E, 38.48°N and 2.85°E, 38.53°N (Figure 9A). At 85 m depth w_QG has maximum values ranging between −7.7 and 7.9 m day⁻¹. However, a sensitivity test reveals the dependence of the extreme values of w_QG to the correlation scales used to interpolate the original density field (Table 3). An increase of the correlation scales from 20 to 25 km would induce a decrease on the extreme values of w_QG of about 50%. Note that applying smaller correlation scales (within the range allowed by the sampling) in the optimal interpolation enables the representation of smaller-scale variability in the interpolated fields, and this may include higher ageostrophic and divergent motions that would translate to larger vertical velocities, consistent with the results from the sensitivity test. A criterion used to define the correlation scales is the estimation of the empirical correlation (Gomis et al., 2001; Pascual et al., 2004). Figure 10 shows the correlation of the original density field at 100 m depth, and the correlation of the drifter velocities (which have higher resolution than the CTD casts). The Gaussian function that visually better resembles the empirical correlations is likely to be the one with a characteristic scale of 20 km, especially for the meridional component of the drifter velocity and for density. This analysis supports the correlation scales used in the optimal interpolation.

The spatial variability of nitrate and chl-a concentrations is consistent with the vertical velocity field calculated from CTD data (Figures 9B–E). An enhancement of nitrate concentration within the photic layer was detected at some stations located close to the upwelling cell: concentrations higher than ~5 μmol l⁻¹ at ~75 m are found in the 38.42°N and 38.50°N zonal sections at the stations located at ~2.90°E and ~3.00°E, and also in the 38.50°N zonal section at the station located at ~3.10°E (Figures 9C,D). This enhanced concentration of nitrate within the photic layer may indicate a recent input of nutrients from deeper layers induced by the upwelling cell centered at 3.05°E, 38.42°N (Figure 9A). A vertical section of chl-a concentration from the SeaExplorer glider return transect shows higher values between 100 and 180 m depth than in the surroundings at ~2.85°E, 38.50°N. This may be the signal of a chl-a subduction related to the western downwelling cell centered at 2.84°E, 38.53°N.

3.2. Scales Detected by the Slocum Glider Along the Sentinel-3A Track

3.2.1. Glider Data vs. Altimetry Data

Figure 11 shows the along-track ADT from the Sentinel-3A satellite mission and the DH computed from glider observations, together with the ADT from the gridded altimetric products interpolated onto the glider transect. The glider DH partially captures the signature of two positive anomalies located south of the Balearic Islands at ~38.00°N and ~38.50°N. These two features are located around 0.5 degrees apart in longitude, thus, the length scale resolved by the glider data in this region has a radius of ~13 km, the same order of magnitude of the Rossby radius of deformation in the western Mediterranean Sea (Escudier et al., 2016; Barceló-Llull et al., 2019). Not surprisingly, these two signatures are not captured by the gridded ADT interpolated onto the glider transect (brown line in Figure 11) due to both its lower spatial resolution and the correlation length scale applied to construct the gridded product (100 km). In consequence, the gridded ADT anomaly shows one single positive anomaly and does not resolve the variability at the glider scales. The Sentinel-3A ADT properly represents the signal of these two positive anomalies with some differences in amplitude and position with respect to the signals detected in glider data. The amplitude of the northern signature is higher in Sentinel-3A ADT anomaly while the opposite occurs in the southern feature. In addition, the maximum anomalies of both fields are slightly displaced and the length scale solved by the Sentinel-3A dataset decreases to ~8 km. The glider completed the transect in almost 6 days so these differences may be due to the time lag between the Sentinel-3A and glider observations. A vertical line in Figure 11 shows the location of the glider during the Sentinel-3A overpass. North and south of this point the glider sampling has a time difference varying from hours to days with respect to the Sentinel-3A sampling as we approach to the limits of the transect. This fact explains the discrepancies in height of up to 2.5 cm observed in the southern feature, and also the displacement of its maximum, captured by both platforms.

3.2.2. Glider and CTD Dynamic Height

To evaluate the contribution of the deeper layers on the estimation of DH, Figure 12A shows the DH computed using Leg 1 CTD data interpolated onto the glider track and assuming the reference level at 650 dbar (orange line) and 1,000 dbar (blue line). The mean values of DH are subtracted to compare the variability of both lines, which is almost identical with an absolute maximum around 38.60°N and minimum values at the bounds of the transect. Thus, the contribution of the water column between 650 and 1,000 m depth to the computed DH is a constant value that increases the absolute value of DH, but does not change the horizontal variability. As a conclusion, in the study area, the contribution of the deeper layers in the DH anomaly can be considered negligible.

The DH computed using the Slocum glider return data is compared to the DH calculated from Leg 2 CTD data along the four transects parallel to the glider track (Figures 12B,C). Overall, the DH shows weak features with magnitudes ranging between 2 and 3 cm from both platform observations. The spatiotemporal variability of the circulation in the region of study introduces the differences observed between the four CTD transects, as they are sampled in different locations and times. The fine-scale signal detected by both platforms shows an evolution between transect 4, sampled on 13-14 May 2018, and transect 1, sampled on 16 May 2018. The maximum of DH located at ~38.50°N in transect 4 and in the glider data is displaced to ~38.25°N in transects 2 and 1. The CTD DH more similar to the glider DH is the one computed along transect 3 (green line in Figure 12C). Both have similar values from 38.50°N toward the north, even though the data have been collected ~10 km apart. The reason of this similarity is that both platforms sampled the region almost simultaneously and they could capture the same signature. Note that even if transect 3 DH has similar values than glider DH north of 38.50°N, transect 4 has a maximum of DH located at the same latitude as the glider data. The DH calculated using transect 2 CTD data (yellow line in Figure 12C) shows a different signal than the signature observed by the glider due to the difference in time of both samplings, even though data were collected along the same track.

4. Discussion and Conclusions

We have analyzed the impact of the scales resolved by different observational platforms on the representation of the horizontal velocity field associated with a meandering front observed from the surface to ~150 m depth. We compared the surface geostrophic velocity field obtained from present-day altimetric products with (i) the geostrophic velocity calculated from CTD measurements, (ii) the total horizontal velocity measured by the ADCP and (iii) drifter trajectories. The CTD-derived geostrophic velocity and ADCP horizontal velocity have similar distributions with maximum root mean square differences between both fields of the order of 0.06 m s⁻¹ and correlation coefficients of ~0.9 in the upper layers (Figure 6). The differences between ADCP and geostrophic velocities could be attributed to different sources of error that arise from the different instruments and procedures used to obtain these fields: level of no motion used to calculate the geostrophic velocity, reduced number of CTD casts between 500 and 1,000 m depth, higher resolution of ADCP data along meridional transects, errors associated with the optimal interpolation, and presence of internal waves and inertial motions in the ADCP velocity field. The small differences observed between both fields, with maximum RMSD of the order of 0.06 m s⁻¹, and the different sources of error that may introduce these differences, suggest that the geostrophic component of the horizontal flow may dominate the local dynamics of this region down to scales of 20 km. The Rossby number computed from each velocity field has absolute values lower than 0.32 at the scales resolved by the sampling (20 km wavelength). The present-day resolution of altimetric products precludes estimation of the horizontal velocities obtained from the ADCP and CTD sampling at scales of ~20 km. This is corroborated by the dispersion of drifters released within the study region: while the CTD-derived and ADCP velocities are in good agreement with the drifter trajectories, the altimetric currents cannot justify their eastward displacement. Quantitatively, the CTD-derived geostrophic velocity has a higher correlation with the drifter velocities than the other fields. This suggests that the geostrophic velocity calculated from CTD observations, which are sampled with a resolution similar to the target SWOT resolution, correctly represents the actual currents that drive the drifters over the sampling domain. In addition, note that the CTD and ADCP observations reconstructed through optimal interpolation are collected during 4.5 days, and they are assumed to be quasi-synoptic to apply the interpolation (e.g., Rudnick, 1996; Pascual et al., 2004; Barceló-Llull et al., 2017; Ruiz et al., 2019). This approximation may introduce errors to the reconstructed fields and derived variables due to the lack of synopticity of the observations (Allen et al., 2001; Gomis et al., 2005). The correspondence detected between the trajectory of in situ drifters and the CTD-derived and ADCP velocity fields supports the approximation of quasi-synoptic observations in the PRE-SWOT experiment.

The Lagrangian analysis performed by comparing the relative separation of in situ drifters with the numerical drifters derived from the assumed stationary CTD and ADCP velocity fields is quite consistent in showing that the region of study has strong frontogenetic dynamics. Interestingly, these dynamics fall into a “dark spot” for altimetry due to its resolution limitations. More specifically, the FSLEs derived from the CTD and ADCP velocity fields are high (reaching maximum values of 0.72 and 0.78 days⁻¹, which is typical of boundary currents estimated from altimetry: see e.g., Figure 1 in Hernández-Carrasco et al., 2012), and approximately double that measured from altimetry in this region. The exponential rate of separation estimated for two in situ drifters is similar to the CTD- and ADCP-derived FSLE maximum values, while nadir altimetry FSLEs misrepresent this value. In addition, the position of the Lagrangian front detected in the ADCP- and CTD-derived FSLE maps is consistent with the in situ drifter trajectories, while the altimetry-derived Lagrangian front, that should represent a transport barrier, is crossed by the eastern in situ drifters. The results from the Lagrangian analysis support the approximation of quasi-synoptic observations considered to optimally interpolate the CTD and ADCP observations, and the assumption of stationary velocity fields used to integrate the numerical drifters. We consider that the approximation of stationarity is applicable in our PRE-SWOT Lagrangian analysis because (i) the duration of the integration is short (6 days or less) and (ii) we are not analyzing eddy retention but only the position and intensity of Lagrangian fronts (d'Ovidio et al., 2013). We are looking forward to future satellite missions like SWOT to have synoptic, time-varying maps of the circulation at these spatial scales. For Lagrangian studies, these future observations will allow us to explore the evolution in time of transport barriers, and see the implication of their dynamics in terms of biogeochemistry and biodiversity, possibly testing hypotheses that at the moment can only be analyzed on modeled experiments (e.g., Lévy et al., 2015).

The rapid relative separation of two drifters along the FSLE front indicates that horizontal stirring stretches any tracer present in the region along that front. The estimation of the backward FSLE along the Lagrangian front provides in fact the intensity of the convergent dynamics acting in the direction orthogonal to the front, and hence the rate at which the passive tracer gradients intensify. The presence of these frontogenetic dynamics induced by horizontal stirring is consistent with the presence of the density gradient, which is the source of the vertical velocities derived through the QG omega equation, as well as with the front visible in the satellite-derived chl-a images. Note that ADCP- and CTD-derived FSLEs provide similar results. This observation suggests again that the dynamics of this region are mostly in geostrophic balance. The Lagrangian front observed during our field experiment is an example of what we seek to resolve in the near future from high-resolution altimetric observations that SWOT is expected to provide, and that today is only accessible from in situ observations.

The vertical velocity field was calculated through integration of the QG omega equation using the CTD-derived geostrophic velocity and density fields. The w_QG distribution is characterized by an upwelling cell surrounded by two downwelling cells with a magnitude that is maximum at ~ 85 m depth and decreases downward and upward. The magnitude of w_QG depends on the correlation scales used to interpolate the original fields. Using correlation scales of 20 km, w_QG has maximum and minimum values of 7.9 m day⁻¹ and −7.7 m day⁻¹, respectively. With correlation scales larger than 30 km the magnitude decreases more than 60% for minimum values and 80% for maximum values. The Gaussian function that is most similar to the empirical correlations, estimated using the original density CTD data and both components of the drifter velocities, is that with a standard deviation of 20 km, supporting the choice of spatial scales used for the optimal interpolation. We have analyzed the spatial variability of nitrate and chl-a in relation to the vertical velocity field. An enhancement of nitrate concentration within the photic layer was detected at two stations located close to the upwelling cell. This may indicate a recent input of nutrients from deeper layers induced by the upwelling cell. Chl-a concentration derived from SeaExplorer glider measurements had high values at ~2.85°E, 38.50°N between 100 and 180 m depth. This may be the signal of a subduction of chl-a related to the western downwelling cell centered at 2.84°E, 38.53°N.

The DH estimated from glider observations shows the signature of two positive anomalies separated by 0.5 degree in longitude (Figure 11), thus, the length scale resolved by the glider data has a radius of ~13 km, the same order of magnitude of the Rossby radius of deformation in the western Mediterranean Sea (Escudier et al., 2016; Barceló-Llull et al., 2019). While the altimetric gridded ADT does not capture this variability due to its lower spatial resolution and the correlation length scale applied to construct the product, the Sentinel-3A ADT correctly represents the signal of these two positive anomalies, with some differences in amplitude and position with respect to the signals detected by the glider. The positive anomalies detected in the Sentinel-3A ADT and glider DH are slightly displaced, and the length scale resolved by the Sentinel-3A dataset decreases to ~8 km. These differences may be due to the time lag between the glider and Sentinel-3A observations.

A comparison between the CTD and glider observations reveals the temporal variability of the estimated DH, which changes over the four transects sampled in Leg 2. The CTD transect with a DH variability similar to that calculated from glider observations was the one that sampled the water column almost simultaneously with the glider, even though they were 10 km apart. The two CTD transects that had a larger time lag with respect to the glider observations have larger differences in DH variability. A sensitivity test of the DH estimated using CTD observations down to 650 and 1,000 m shows that the contribution of the deeper layers is an almost constant mean value and does not introduce significant horizontal variability. In the same region of study, Aulicino et al. (2018) previously obtained similar results from a DH sensitivity test calculated by assuming different reference levels between 200 and 900 m depth. This suggests that, when a quasi-synoptic sampling is needed, in the study region CTD casts do not need to sample the water column down to 1,000 m depth in order to estimate DH. This has the advantage that it is possible to sample small study regions more quickly.

We have used in situ observations from different observational platforms to demonstrate the limitations of present-day nadir altimetry to resolve the circulation in a region of the western Mediterranean Sea with relatively moderate dynamic activity (Figure 2). The altimetric observations are too coarse to capture the evolving circulation over the period of the cruise: some features are missing, others are misplaced, and most are significantly weaker than in the in situ observations. One important consequence is that dispersion processes are incorrectly represented, as evidenced by the drifter experiment. This misrepresentation of lateral stirring results in data suggesting much weaker fronts. As a consequence, the frontogenetic region sampled during the PRE-SWOT cruise, which has a clear signal in biochemistry, has almost no signal in altimetric maps. The implication is that, at present, we must use in situ instruments to characterize the physical signatures of the features we have surveyed. Interestingly, the component which is missing in the horizontal circulation captured by nadir altimetry appears to be mainly geostrophic down to scales of 20 km wavelength in the region of study and during the sampling period, and hence it should be resolvable by future altimetric satellites with higher-resolution such as SWOT. At smaller scales it is expected that it will be a significant departure from geostrophy and, therefore, other satellite missions measuring ocean surface currents are required (Gommenginger et al., 2019).

The integrated multi-platform and interdisciplinary sampling of the region south of the Balearic Islands has revealed interesting results about the potential scales that may be resolved by SWOT, and also the limitations of current nadir altimetry. Our study quantifies and shows that the geostrophic approximation at scales of 20 km in the region of study is a good approximation. We have also calculated the QG vertical velocity field and detected a correspondence with the spatial distribution of nitrate and chl-a concentrations. Through the multi-platform and interdisciplinary PRE-SWOT experiment in collaboration with the French campaign PROTEVSMED-SWOT, we have evaluated (i) the 3D dynamics of the study region, which has been previously scarcely sampled with only some recent studies using high-resolution glider data (Heslop et al., 2017; Aulicino et al., 2018), and (ii) the coupling between fine-scale physics and biological measurements (Tzortzis et al., 2021). With this joint experiment we have evaluated the dynamics of the western Mediterranean SWOT crossover and gained experience in multi-institution and multi-platform campaign coordination.

The Mediterranean Sea has lower first Rossby radius of deformation (~5–15 km, Beuvier et al., 2012; Escudier et al., 2016; Barceló-Llull et al., 2019; Kurkin et al., 2020) than other regions of the global ocean at mid-latitudes (~20–40 km, Chelton et al., 1998). This means that the spatial scales are smaller and nadir altimetry only captures a small part of the mesoscale dynamics. Hence, in this region SWOT is expected to provide a substantial improvement in the observation of these scales (d'Ovidio et al., 2019). In addition, the Mediterranean Sea is characterized by low tides and relatively weak internal wave dynamics that make this region interesting for interdisciplinary studies (d'Ovidio et al., 2019). Note that due to the particular characteristics of the Mediterranean Sea, the conclusions from this study are not extendable to other regions with different dynamics.

In future experiments, it may be interesting to explore the region northwest of the Balearic Islands, which has been documented to have higher dynamic activity than the region studied here and is under the SWOT path during the fast-sampling phase. In future work, we will use high-resolution regional ocean circulation models to conduct Observing System Simulation Experiments to (i) evaluate the quasi-synopticity of the observations collected in multi-platform experiments similar to PRE-SWOT, (ii) compare the optimal interpolation with other methods of reconstruction that can account for the temporal variability of the observations, and (iii) study optimal sampling strategies for multi-platform experiments aimed at SWOT validation.

Data Availability Statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: http://dx.doi.org/10.20350/digitalCSIC/8640.

Author Contributions

BB-L led the writing of the manuscript. AP designed, obtained the funding for and led the PRE-SWOT experiment. BB-L, AS-R, and EC post-processed the in situ observations. FD and GF conducted the Lagrangian analysis. All authors contributed to the analysis of the results and in the collection of the in situ observations.

Funding

The PRE-SWOT experiment was funded by the Spanish Research Agency and the European Regional Development Fund (AEI/FEDER, UE) under Grant Agreement (CTM2016-78607-P). BB-L was supported by a postdoctoral contract in the framework of the EuroSea project that has received funding from the European Union's Horizon 2020 research and innovation programme under Grant Agreement No. 862626. This study is also a contribution to the SWOT Science Team (TOSCA-CNES VERSO and TOSCA-CNES BIOSWOT-AdAC project proposals) and to the EuroSea project (Horizon 2020 research and innovation programme under Grant Agreement No. 862626).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher's Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

We express our gratitude to the technical staff (UTM-CSIC) and crew of the R/V García del Cid for supporting our work at sea and to Andrea Cabornero, Esther Capó, Noemí Calafat, Tim Toomey, Juan Gabriel Herná́ndez, David Roque, Albert Miralles, Manuel Rubio, and Marc Torner for their participation in the acquisition of the observations. We thank SOCIB staff for the MIO's glider deployment.

References