# Wind-Forced Submesoscale Symmetric Instability around Deep Convection in the Northwestern Mediterranean Sea

Anthony Bosse, Pierre Testor, Pierre Damien, Claude Estournel, Patrick Marsaleix, Laurent Mortier, Louis Prieur, Vincent Taillandier

► To cite this version:

Anthony Bosse, Pierre Testor, Pierre Damien, Claude Estournel, Patrick Marsaleix, et al., Wind-Forced Submesoscale Symmetric Instability around Deep Convection in the Northwestern Mediterranean Sea. Fluids, 2021, 6 (3), pp.123. 10.3390/fluids6030123. hal-03199242

HAL Id: hal-03199242

https://amu.hal.science/hal-03199242

Submitted on 15 Apr 2021

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
Article

# Wind-Forced Submesoscale Symmetric Instability around Deep Convection in the Northwestern Mediterranean Sea

Anthony Bosse <fnref n="1" />,* , Pierre Testor <fnref n="2" />, Pierre Damien <fnref n="3" />, Claude Estournel <fnref n="4" />, Patrick Marsaleix <fnref n="4" />, Laurent Mortier <fnref n="2" />,<fnref n="5" />, Louis Prieur <fnref n="6" /> and Vincent Taillandier <fnref n="6" />

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

check for
updates

**Citation:** Bosse, A.; Testor, P.; Damien, P.; Estournel, C.; Marsaleix, P.; Mortier, L.; Prieur, L.; Taillandier, V. Wind-Forced Submesoscale Symmetric Instability around Deep Convection in the Northwestern Mediterranean Sea. *Fluids* 2021, *6*, 123. https://doi.org/10.3390/fluids6030123

Academic Editors: Jonathan Gula, James C. McWilliams and Leif N. Thomas

Received: 19 January 2021

Accepted: 11 March 2021

Published: 17 March 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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

**Copyright:** © 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

**Abstract:** During the winter from 2009 to 2013, the mixed layer reached the seafloor at about 2500 m in the northwestern Mediterranean Sea. Intense fronts around the deep convection area were repeatedly sampled by autonomous gliders. Subduction down to 200–300 m, sometimes deeper, below the mixed layer was regularly observed testifying of important frontal vertical movements. Potential Vorticity dynamics was diagnosed using glider observations and a high resolution realistic model at 1-km resolution. During down-front wind events in winter, remarkable layers of negative PV were observed in the upper 100 m on the dense side of fronts surrounding the deep convection area and successfully reproduced by the numerical model. Under such conditions, symmetric instability can grow and overturn water along isopycnals within typically 1–5 km cross-frontal slanted cells. Two important hotspots for the destruction of PV along the topographically-steered Northern Current undergoing frequent down-front winds have been identified in the western part of Gulf of Lion and Ligurian Sea. Fronts were there symmetrically unstable for up to 30 days per winter in the model, whereas localized instability events were found in the open sea, mostly influenced by mesoscale variability. The associated vertical circulations also had an important signature on oxygen and fluorescence, highlighting their under important role for the ventilation of intermediate layers, phytoplankton growth and carbon export.

**Keywords:** deep convection; submesoscale dynamics; frontal instability; subduction

## 1. Introduction

The Gulf of Lion in the northwestern Mediterranean Sea (NW Mediterranean) is one of the few oceanic regions where intense atmospheric forcing and oceanic preconditioning allow vertical mixing to reach large depths (>500 m to the seafloor at 2500 m [1]). This phenomenon, known as "open-ocean deep convection" (see Marshall and Schott [2] for a review), is a key process of the ocean thermohaline circulation. Deep convection is characterized by localized and intense diapycnal mixing events [3] important for the ocean ventilation [4–6], subsequent intense phytoplankton spring blooms [7–9] and carbon sequestration [10]. In accordance with idealized results [11], the net downwelling induced by deep convection however occurs along the boundary circulation where a large amount

<footnote n="1">¹ Mediterranean Institute of Oceanography, Aix Marseille University, Université de Toulon, CNRS, IRD, MIO UM 110, F-13288 Marseille, France</footnote>

<footnote n="2">² Laboratoire d'Océanographie et du Climat: Expérimentations et Approches Numériques, Sorbonne Universités, UPMC Université Paris 06, IPSL, CNRS, IRD, MNHN, LOCEAN UMR 182, F-75252 Paris, France; testor@locean-ipsl.upmc.fr (P.T.); laurent.mortier@ensta-paristech.fr (L.M.)</footnote>

<footnote n="3">³ Department of Atmospheric and Oceanic Sciences, University of California, Los Angeles, CA 90025, USA; pdamien@ucla.edu</footnote>

<footnote n="4">⁴ Laboratoire d'Études en Géophysique et Océanographie Spatiales, CNRS, CNES, IRD, Université de Toulouse, LEGOS UMR 5566, F-31400 Toulouse, France; claude.estournel@legos.obs-mip.fr (C.E.); patrick.marsaleix@legos.obs-mip.fr (P.M.)</footnote>

<footnote n="5">⁵ ENSTA-Paristech, F-91120 Palaiseau, France</footnote>

<footnote n="6">⁶ Laboratoire d'Océanographie de Villefranche, Sorbonne Universités, CNRS, IMEV-LOV UMR 7093, F-06230 Villefranche-sur-Mer, France; prieur@obs-vlfr.fr (L.P.); taillandier@obs-vlfr.fr (V.T.)</footnote>

<footnote>* Correspondence: anthony.bosse@mio.osupytheas.fr</footnote>
of heat is released to the atmosphere [12]. In the western Mediterranean Sea, the coupled
process of deep convection has been variable over the last decades [13] and affected by
abrupt changes in ocean heat and salt contents [14–16]. Numerical models predicted a
decrease in convection activity due to climate change [17,18] with important consequences
on the marine ecosystems [19,20].

Open-ocean deep convection has traditionally been described by three phases [21,22],
however overlapping in space and time: a preconditioning phase characterized by a
basin-scale cyclonic and more localized (sub)mesoscale circulations associated with a
doming of isopycnals toward the surface; a vigorous mixing phase triggered by intense
heat losses and evaporation to the atmosphere, typically 100–1000 W m⁻² [23] and caused
by intense cold/dry wind events; a spreading phase of the newly-formed water mass along
neutral density layers initiated when ocean restratification overcomes the destabilizing
atmospheric forcing.

During the mixing phase, a mixed patch forms in the area of deep mixing separated
from stratified waters by density fronts. Diapycnal mixing increases the available potential
energy of the mixed patch feeding eddy kinetic energy by baroclinic instability of the rim
current [24–27]. In particular, Submesoscale Coherent Vortices [28] are generated. They are
non-linear eddies in cyclogeostrophic balance with O(1) Rossby number, peak velocities at
great depth and a weakly stratified core spreading the newly-formed water mass over long
distances due to their extended lifetime (>1 year) [29–37].

The horizontal circulation in northwestern Mediterranean Sea follows a basin scale cyclonic gyre [38,39]. In particular, the “Northern Current” (NC) establishes in the north-west of Corsica from the merging of the “West Corsica Current” (WCC) with the Corsica Channel throughflow. The NC then flows along the continental slope from the Ligurian Sea to the Catalan Sea (Figure 1). The cyclonic circulation observed as a permanent feature is mostly in geostrophic balance creating a dome of high density waters in the central basin. The central area of the gyre circulation is thus preconditionned to deep convection [2,21]. The frontal zone of the NC separates low salinity surface waters of Atlantic origin (Atlantic Water, AW) from higher salinity waters in the central basin modified through the action of vertical mixing and evaporation. The front has a thickness of more than 200 m throughout the seasons and a 25-km wide geostrophic jet with peak velocities of 0.3–0.5 m s⁻¹ roughly aligned with the 2000-m isobath (see Prieur et al. [39] for a recent review). Under the effects of surface fluxes (heat, evaporation, wind) and topography [40], the NC undergoes baroclinic instability and meanders with a timescale of O (5–10) days [24,41,42]. This instability tends to restore equilibrium by slanting isopycnal and generating mesoscale meanders and eddies.

In addition to baroclinic instability acting at mesoscale, frontal regions are also prone to
an important dynamics at submesoscale generating intense vertical velocities, as suggested
by observations and models [43]. The submesoscale regime classically refers to unbalanced
flows characterized by O(1) Rossby number and scales smaller than the first baroclinic de-
formation radius [44]. Significant vertical water displacements of O (10–100) m day⁻¹ occur
in such regions [45–49]. High levels of turbulent diapycnal mixing in active fronts [50,51]
foster heat transfers from the surface of the ocean to the interior [52,53]. Due to their action
at timescales relevant for biology, the vertical injection of nutrients at submesoscale has
strong major impacts on phytoplankton growth and marine ecosystems [54–57]. Fronts are
also crucial to the formation of low Potential Vorticity (PV) “mode waters” [58,59] important
for the ocean ventilation [60], deep convection [61,62] and carbon sequestration [63].

Potential vorticity is a key variable for the dynamics of frontal regions. By definition, PV is conserved along adiabatic streamlines [64], i.e., in absence of forcing and dissipation. Hence, major sources and sinks of PV are found where isopycnals outcrop at the air–sea interface, or in the bottom boundary layer [65]. Frontal destruction of PV at the surface can be caused by thermodynamic heat losses [66], or mechanical friction due to winds [67]. In particular, Ekman transport by down-front winds (i.e., along the geostrophic flow) advects water from the dense to the light side of a front [68]. This mechanism can cause important
buoyancy loss equivalent to O (1000) W m⁻² and localized at the front. These fluxes can be significantly larger than air/sea net heat fluxes [69–71]. The destruction of PV by winds was also shown to play a role in deep water formation in the northwestern Mediterranean Sea [61].

Alterations of PV may expose frontal regions to instabilities. In particular where potential vorticity has opposite sign of *f* the Coriolis parameter (i.e., negative in the northern hemisphere), symmetric instability (SI) arises [72]. In the ocean, this instability drives O (0.1–1 km) overturning cells aligned with isopycnals, symmetric in the along-front direction and growing in short O (1–10 h) timescale [73,74]. Large-Eddy Simulations have shown an increase in dissipation rate of turbulent kinetic energy due to SI [75], as well as primary production stimulation at high latitude fronts [54]. SI also generates important mixing in the bottom boundary layer of intense currents [76,77]. It is thus an important submesoscale process, yet still poorly observed and characterized.

Favorable conditions for SI were observed along the frontal region of western boundary currents like the Gulf Stream [73,78] and the Kuroshio [79]. However, these observations are difficult to collect since such unstable fronts are caused by intense atmospheric forcing causing rough sea conditions in which research vessels can hardly operate. This challenge has been met by ocean gliders [80] which allow high horizontal resolution autonomous measurements, from 1 km to 2 km between consecutive 1000 m profiles. Gliders have been widely used to characterize frontal dynamics [69,70,78,81–85] and can overcome winter harsh sea conditions, although their displacement through the water of 20–40 km day⁻¹, fast for the mesoscale of the NC [42], can lead to aliasing of rapidly evolving submesoscale structures.

While observational evidence of SI remains scarce, despite growing interest in characterizing and quantifying submesoscale regime from in situ observations, the increase and repeat use of autonomous platforms such as gliders during the last decade has revealed new insights in submesoscale flows. In this study, we use an extensive observational data set of the fronts surrounding the deep convection area in the northwestern Mediterranean Sea, mainly associated with the NC front, see Figure 1, to evaluate and characterize frontal submesoscale dynamics in interaction with strong winds and heat losses. The PV is further diagnosed in both observations and a very high resolution numerical model. We then assessed the role of down-front winds in generated favorable conditions to SI in terms of PV destruction and Ekman buoyancy fluxes. Key regions regularly prone to SI are linked to topographically-steered circulation and dominant winds. Finally, we discussed the role of SI on vertical circulations at fronts and biogeochemical implications.
<figure><img src="image_4.png" /><figcaption>**Figure 1.** Maps of glider observations in the northwestern Mediterranean Sea from November 2008 to October 2013. The four endurance lines (T00, T01, T02 and T03) and mooring sites (LION and DYFAMED) are parts of the MOOSE (Mediterranean Ocean Observing System, https://www.moose-network.fr (accessed on 16 March 2021), see also Coppola et al. [86] for an overview). The color correspond to the number of profiles collected in 10-km boxes in the northwestern Mediterranean Sea. The upper panel further represents the monthly number of profiles collected along each endurance line during the study period. The main features of the circulation are represented by red arrows: Norther Current (NC), North-Balearic Front (NBF, in dashed line because of its variable and unstable pathway at mesoscale), West Corsica Current (WCC). GoL stands for “Gulf of Lion” and Lig. Sea for “Ligurian Sea”.</figcaption></figure>

**2. Materials and Methods**

**2.1. Glider Measurements**

Autonomous ocean gliders are an essential component of Global Ocean Observing System [80]. They sample the ocean along saw-tooth trajectories between the surface and a maximal depth of 1000 m traveling at horizontal speed of 20–40 km per day. Deep profiles are typically separated by 1–3 km and 1–4 h. Gliders can thus sample fronts and eddies characterized by typical horizontal scale of 10–50 km in about 1–3 days. In this study, we used data collected by gliders from November 2008 to October 2013 in the northwestern Mediterranean Sea. This represents a total of about 35,000 profiles from 72 glider missions. Most of the glider data have been collected in the framework of the MOOSE (Mediterranean Ocean Observing System, https://www.moose-network.fr/ (accessed on 16 March 2021), the long-term observatory of the northwestern Mediterranean Sea [86]. Those were collected along repeated transects (see Figure 1: T00, 13061 profiles—33 missions; T01, 3702 profiles—33 missions; T02, 5970 profiles—23 missions; T03, 4685
profiles—21 missions). During the study period, an exceptionally intense deep convection activity was observed in the Gulf of Lion with bottom-reaching mixed layer observed at the LION moored time series during five consecutive winters [1] (Figure 2). The MOOSE endurance glider transects—with experimental transects starting in 2006 and regular missions since 2008—were designed to monitor the convection area, characterize the basin circulation and visit the LION and DYFAMED mooring lines in operation since 2008 [87] and 2009 [88], respectively, (Figure 1).

<figure><img src="image_5.png" /><figcaption>**Figure 2.** (a) Net surface heat loss (left axis) and wind stress amplitude (right axis) from CNRM-ALADIN reanalysis at the LION mooring location. Red and blue bars show monthly net heat flux with mean seasonal cycle represented by the black line. Wind stress amplitude were daily/monthly averaged in orange/gray. (b) Observed potential temperature at LION moored time series with MLD from [1] represented in white. The period of numerical simulation from October 2010 to December 2011 is shown in red. (c) Potential temperature simulated by SYMPHONIE at LION mooring site during the winter 2011. Simulated MLD at LION is shown by solid white line. MLD is here from a 0.1 °C threshold criterion on temperature profile. The dashed line represents the 95th percentile of the MLD in a box 41.5–42.5° N/4–5° E around the mooring position representative of the deep convection area. Observed MLD is shown in red as reference.</figcaption></figure>

Gliders were equipped with Conductivity-Temperature-Depth (CTD) sensor sampling at 2–8 s frequency allowing a vertical resolution of 0.5–2 m. Offset of temperature and salinity were determined for all missions following the approach used in Bosse et al. [34]. In addition, thermal lag effects affecting salinity measurements by unpumped CTD probe in strong thermoclines have been corrected following Garau et al. [89].

Dead reckoning navigation compared to GPS fixes at surface enables gliders to estimate a mean oceanic current averaged during each dive, hereafter referred as depth-
average current (DAC). Compass calibration has been carried out before each deployment allowing DAC to be used as reference for geostrophic velocities with typical precision of 1 cm s⁻¹ [90].

Gliders additionally performed oxygen and/or optical fluorescence/turbidity measurements. The fluorescence data shown here were corrected for non-photochemical quenching [91] and adjusted against satellite measurements [92]. This method has been validated against high-performance liquid chromatography measurements [9].

## 2.2. Numerical Model

SYMPHONIE is an Ocean General Circulation Model (OGCM) based on the Navier–Stokes primitive equations solved on an Arakawa curvilinear C-grid under the hydrostatic and Boussinesq approximations [93]. The regional configuration encompassed a significantly larger region than the deep convection zone from the Balearic Sea to the Tyrrhenian Sea (Figure 1). The vertical grid had 40 generalized σ-levels with varying thickness from 2.5 m near the surface to 150 m to the bottom (and about 15 m at 100 m along the Northern Current). The horizontal grid size was 1 km (effective resolution of about 5 km [94]) in order to resolve mesoscale eddies and fronts as well as the upper bound of the submesoscale (the first baroclinic deformation radius being about 5–10 km in this region [31]). This regional configuration has been validated against numerous in situ observations (Figure 2) and reproduced long-lived 5 km radius SCVs formed by deep convection [36]. It was also coupled to a biogeochemical model in order to study vertical transfers of nutrients and carbon export by deep convection [95,96].

Mixing processes associated with deep convection involves small-scale downwelling plumes in non-hydrostatic balance characterized by <1 km diameter and intense vertical velocities up to 10–15 cm s⁻¹ [97]. These are not resolved by hydrostatic OGCM with kilometric resolution. Here, a non-penetrative adjustment algorithm redistributes surface buoyancy losses throughout the mixed layer [98] while the vertical diffusion coefficient computed by the Gaspar et al. [99] turbulence closure scheme increases up to O (1 m² s⁻²) generating instantaneous and adiabatic parametrized mixing. Radiative, momentum, heat and salt fluxes at the air–sea interface were provided by the dynamical downscaling at 50-km of ERA-Intrerim reanalysis over Europe performed by the CNRM–ARPEGE atmospheric model [100].

The numerical setup is similar to the one described in Damien et al. [36]. Mercator PSY2V4R3 at 1/12° (7 km) resolution [101] were used to prescribed initial and open boundary conditions. The simulation period extended from 1st of October 2010 to 31st of December 2011. Daily average outputs were saved for the whole simulation period. A complete seasonal cycle of winter convection was simulated and a focus was put on the last 12 months (January to December 2011).

## 2.3. Thermodynamic and Wind Effects of Surface Forcing

Surface forcing (air/sea heat fluxes and wind) was also retrieved during the study period from CNRM–ALADIN reanalysis (a downscaling of ERA-Interim atmospheric model at 12-km resolution with 6h outputs) and used to interpolate atmospheric forcing in space and time along glider sampling.

A complex principal component analysis was performed from daily average wind stress from CNRM-ALADIN during the period 2009–2013 in order to evaluate the dominant wind forcing patterns.

During the winter, large amounts of latent and sensible heat are transferred from the ocean to the atmosphere at the basin scale. As mentioned earlier, winds can also have a dominant effect on the upper ocean buoyancy budget through the action of Ekman transport across fronts. The Ekman buoyancy flux can be quantified as an equivalent heat flux according to [102]:

$$Q_{Ek} = \frac{\rho C_p}{\alpha g} \vec{M}_e \cdot \vec{\nabla} b|_{z=0} = \frac{C_p}{\alpha g f} (\vec{\tau} \times \hat{\mathbf{z}}) \cdot \vec{\nabla} b|_{z=0} \quad (1)$$
with $\vec{M}_e$ the Ekman transport, $\rho$ the surface density, $b = -g\rho/\rho_0$ the buoyancy of the fluids, $\vec{\tau}$ the wind stress and $f$ the Coriolis parameter. Non-linear Ekman transport could be considered [103] enhancing (resp. decreasing) fluxes in anticyclonic (resp. cyclonic) areas by 30–50% depending on the local Rossby number. However, as the pattern of fluxes induced by winds friction in frontal regions were not affected, the linear expression was kept for the sake of simplicity.

The equivalent Ekman heat flux ($Q_{Ek}$) was estimated in both model and glider data. Due to the fact that gliders sample along their trajectory and are not able to capture the front orientation, this can lead to a wrong sign or amplitude of this parameter as discussed in Thompson et al. [82] for open-ocean fronts. However in the case of topographically-steered geostrophic fronts and glider sampling perpendicular to isobaths (like the MOOSE endurance lines crossing the NC here), the front orientation can be more confidently defined using glider DAC, as also done by du Plessis et al. [69].

## 2.4. Potential Vorticity

For continuously stratified flows, $q$ (in m$^{-1}$ s$^{-1}$) is defined as [104]:

$$q = (f\hat{z} + \vec{\nabla} \times \vec{u}) \cdot \vec{\nabla}b/g = (\zeta^a N^2 + \partial_z u \partial_y b - \partial_z v \partial_x b)/g \quad (2)$$

where $\vec{u} = (u, v, w)$ is the velocity field, $g$ the gravitational acceleration, $\zeta^a = f + \partial_x v s - \partial_y u$ the absolute vorticity and $N^2 \equiv \partial_z b$ the squared Brunt-Väisälä frequency.

Assuming the flow being mainly in geostrophic balance ($\vec{u} = \vec{u}_g$), the thermal wind balance implies $\partial_z u_g = -f^{-1} \partial_y b$ and $\partial_z v_g = f^{-1} \partial_x b$. The geostrophic PV then simplifies as:

$$q = (f \zeta_g^a N^2 - |\vec{\nabla}_h b|^2) / f g = q_N + q_{bc} \quad (3)$$

with $\zeta_g^a$ the absolute vorticity of the geostrophic flow and $\vec{\nabla}_h = (\partial_x, \partial_y)$ the horizontal gradient operator. PV is here separated in two terms: $q_N$ the stratification term linked to stratification and absolute vorticity of the fluid, and $q_{bc}$ the baroclinic term due to the horizontal buoyancy gradient (i.e., front sharpness).

## 2.5. Flow Instabilities

Potential vorticity can be used to characterize the different types of frontal instability. For instance, symmetric instability arises when $q$ has opposite sign of $f$ (i.e., negative in the northern hemisphere) [72]. This is the case of intense fronts where the vertical stratification term in PV is dominated by horizontal buoyancy gradients: $f(f + \zeta_g)N^2 < |\vec{\nabla}_h b|^2$. This instability generates cross-front perturbations aligned with sloping isopycnals resulting in important vertical exchanges and turbulence [54,73].

The PV expression (3) can then be rewritten in a more compact way using the balanced Richardson ($Ri_B = f^2 N^2 / |\vec{\nabla}_h b|^2$) and Rossby ($Ro = \zeta_g / f$) numbers:

$$q = f N^2 (1 + Ro - Ri_B^{-1}) / g \quad (4)$$

Different regimes of instability can be categorized: slantwise convection or symmetric instability (SI) for $q < 0$ and $N^2 > 0$, vertical convection or gravitational instability (GI) for $N^2 < 0$, inertial instability (II) for $Ro < -1$ and mixed regimes (II/SI and GI/SI). These regimes can more easily be classified introducing the Richardson angle [73]: $\Phi_B = \tan^{-1}(-Ri_B^{-1})$. This angle associates a infinite range of Richardson numbers with finite angles useful to describe the different kind of instability that can result. In particular, the criterion for SI then becomes: $-\pi/2 < \Phi_B < \min(\Phi_C, -\pi/4)$ where $\Phi_C = \tan^{-1}(-1 - Ro)$ (see Figure 1 in Thomas et al. [73] for a sketch summarizing all the different cases). We here described flow instabilities in observations and model following this framework.
## 2.6. Estimating PV with Gliders

In the following, the along-front direction was arbitrary defined as the x axis with y pointing in the glider trajectory (right-handed coordinate system). Computing PV requires the estimation of both vertical and horizontal buoyancy gradients. A first approach neglecting cross-track geostrophic velocities ($v = 0$) and gradients ($\partial_x(\cdot) = 0$) yields: $q_g = (f - \partial_y u)N^2/g - |\partial_y b|^2/f_g$, as done in a number of studies [81–83,105].

Small-scale isopycnal oscillations due to unbalanced flows and waves need to be filtered out in glider data before the geostrophic shear could be estimated [106]. Here, this was done using a Gaussian moving average ($\sigma = 2.5$ km) on the buoyancy field as in Bosse et al. [107]. Absolute geostrophic velocities $u$ were then computed by vertically integrating the thermal wind balance along the glider trajectory and using the DAC component perpendicular to the glider track as absolute reference.

As the glider trajectory will likely not cross fronts at exact right angle, geostrophic shear and vorticity will thus be underestimated by a factor $|\sin\theta|$ where $\theta$ is the angle between the glider DAC and the glider track assuming that DAC are representative of the front orientation. This factor is even squared in both PV terms. To avoid large bias, we here adopted the local streamwise coordinate system orientated by glider DACs (details of the method can be found in Todd et al. [78] and in supplementary information of Bosse and Fer [85]). The correcting factor in PV is lower than 2 when the angle $\theta$ is smaller than 45°. We therefore excluded PV values according to this threshold in order to reduce the uncertainties of glider-based PV estimates.

## 2.7. Validation of Glider-Based PV

In order to validate glider-based PV estimates, we have simulated the glider sampling in outputs of the numerical model. The 1-km horizontal resolution of the model is comparable to the glider sampling and able to reproduce a realistic glider sampling. Temperature, salinity and DAC were interpolated in space and time along an actual glider trajectory with a focus on the NC frontal dynamics during an intense wind event (Figure 3). To do so, we used the glider mission MOOSE T02–02, which lasted for 87 days (from 12 November 2010 to 7 February 2011) and sampled four times the northern part of T02 and T03 sections from the Gulf of Lion’s shelf to the LION mooring line (Figure 1).

The geostrophic assumption was first tested by computing total and geostrophic PV calculated using, respectively, the total horizontal velocities Equation (2), and geostrophic ones Equation (3). Model geostrophic velocities were computed by integrating the thermal-wind balance referenced by surface currents given by gradients of sea surface height. The total and geostrophic PV in the model exhibit differences mainly in the upper layer where ageostrophic motions mostly prevail. In the upper 25 m, their deviated from each other by $0.75 \pm 2.7 \times 10^{-11} \text{ m}^{-1} \text{ s}^{-1}$ and by an order of magnitude lower below between 25 and 1000 m: $0.13 \pm 2.4 \times 10^{-12} \text{ m}^{-1} \text{ s}^{-1}$. This difference is about an order of magnitude lower than the PV magnitude observed in frontal regions of about $10^{-11} \text{ m}^{-1} \text{ s}^{-1}$. Thus, the geostrophic PV can capture the main frontal PV structure and the associated instability regimes (77% of the SI occurrences diagnosed with total PV were also found using geostrophic PV).

The same calculation as the one used with observations was then performed with the virtual glider data and compared with total model geostrophic PV in order to assess the effect of the glider sampling and PV calculation (Figure 3). The stratification term $q_N$ is accurately reproduced without large biases (mean deviation of 1–10%, see Table 1) and the cross-stream correction has no effect. The main bias in the PV calculation is however introduced by the baroclinic term $q_{bc}$, whose intensity is underestimated by 30 to 60% in absence of correction (see Table 1). The front width looks slightly wider as the result of the horizontal averaging. The cross-stream geometric correction seems particularly important to correct it, especially in regions subjected to SI where it results in overestimation by 20% of the baroclinic term on average. In the end, the geometrically-corrected total PV is in
much better agreement (1–5% difference compared to 40% without correction). Stable and symmetrically unstable regions are well captured in 90% to 80% of the occurrences.

<table><caption>**Table 1.** Glider-based to total geostrophic PV ratio in the SYMPHONIE model with and without cross-stream correction considered for the front section shown in Figure 3. Numbers are given for the two most representative regions: stable and symmetrically unstable. Accurate detection of stability regimes (in percentage) are also given with reference to the ones obtained using the total PV.</caption><thead><tr><th>Regime</th><th>q<sub>N</sub><sup>g</sup>/q<sub>N</sub><sup>tot</sup></th><th>q<sub>N</sub><sup>g,cor</sup>/q<sub>N</sub><sup>tot</sup></th><th>q<sub>bc</sub><sup>g</sup>/q<sub>bc</sub><sup>tot</sup></th><th>q<sub>bc</sub><sup>g,cor</sup>/q<sub>bc</sub><sup>tot</sup></th><th>q<sup>g</sup>/q<sup>tot</sup></th><th>q<sup>g,cor</sup>/q<sup>tot</sup></th><th>detec %</th><th>detec %, cor</th></tr></thead><tbody><tr><td>stable</td><td>0.99 ± 0.22</td><td>0.98 ± 0.21</td><td>0.42 ± 1.3</td><td>0.47 ± 2</td><td>1 ± 0.25</td><td>0.99 ± 0.32</td><td>92</td><td>92</td></tr><tr><td>SI</td><td>1.1 ± 2.2</td><td>1.1 ± 1.7</td><td>0.72 ± 0.37</td><td>1.2 ± 0.53</td><td>0.61 ± 0.41</td><td>0.97 ± 0.53</td><td>79</td><td>83</td></tr></tbody></table>

<figure><img src="image_6.png" /><figcaption>**Figure 3.** (a) Stratification term of PV ($q_N$) from model outputs interpolated along a glider trajectory crossing the NC front during a down-front wind event. (b) glider-like reconstruction of $q_N$ from virtual glider in the model. The white line in (a,b) shows the MLD. (c,d) Same as (a,b) but for the baroclinic PV term, $q_{bc}$. (g,h) Same as (a,b) for the total PV, $q = q_N + q_{bc}$. The blue contours represent regions subjective to SI according to the balanced Richardson number (see Section 2.5).</figcaption></figure>
To conclude, glider-based PV seems to accurately characterize PV at fronts and detect SI events. Calculation in cross-stream axis provides a better estimate compared to cross-track ones, especially in intense frontal regions. The main difficulty of observing SI events remains the timing of the front crossing, as these features evolve rapidly.

# 3. Results

## 3.1. Deep Convection and Seasonality of Submesoscale

Dry and cold northerly winds ("Mistral" and "Tramontane") prevail in the north-western Mediterranean Sea in winter resulting in intense heat losses to the atmosphere (monthly average of about 250 W m⁻², Figure 2a). The Gulf of Lion region has been subjected to very intense mixing events from 2009 to 2013. The MLD reached the seafloor ever winter at about 2500 m typically around February, as observed by the LION moored time series (Figure 2b). The SYMPHONIE model realistically reproduced the deep convection of winter 2011. The modeled MLD time series at the nearest grid point from the LION mooring site compared reasonably well with a root mean square error (rmse) of 530 m (white line in Figure 2c). A better representation in timing and amplitude of convection events was however found by considering a larger area around the LION mooring site lowering down the rmse to about 340 m (dashed white line in Figure 2c). This could be due to the shape and position of the mixed patch in the model, which is mostly localized along the NC (blue contours in Figure 4a) and does not extend as much offshore, as ocean color seen by satellite revealed. This highlights the local and basin-scale preconditioning of vertical mixing by (sub)mesoscale eddies [34,40]. During the study period, satellite images showed a maximum deep convection area of 15,800–24,100 km² in winter 2011 [1] compared in the model to a maximum area where the MLD exceeded 500 m of 9100 km² reached on 24 January. This is satisfactory for our study of frontal processes in the presence of deep mixed layer, but illustrates the difficulty to represent deep convection in regional and global numerical models and the importance of convection parameterization (see for instance [108]).

The horizontal pattern of vorticity exhibited (sub)mesoscale features such as meanders, fronts and eddies with notably smaller scales in the northern domain where the vertical mixing was the deepest. The dynamics in the convection area had a marked seasonal cycle in the model (Figure 4). Surface relative vorticity was strongly skewed toward positive values. Vorticity significantly larger than $f$ were observed during the winter, while such high values were almost absent during the summer months. Anticyclonic vorticity stayed bounded at around $-0.6f$ in both seasons. At 500 m, the dynamics was less energetic with vorticity distribution tapering during summer (Figure 4f). The skew toward cyclonic vorticity was still present at subsurface but less marked. The distribution showed a clear seasonal increase toward higher vorticity in winter with values up to $|Ro| \sim 0.5$.

Deep convection sharply modifies the ocean dynamics as vertical mixing diminishes the deformation radius in the mixed layer ($R_d = NH/f$ where N is the buoyancy frequency, H the mixed layer and f the Coriolis frequency) and energize small scale dynamics compared to summer stratified conditions [109]. Winter mixing also feeds the available potential energy that is transferred into eddy kinetic energy during the restratification phase [22].
<figure><img src="image_7.png" /><figcaption>**Figure 4.** Relative vorticity normalized by *f* from SYMPHONIE model at 5 m during (a) the winter (snapshot on 24 January, 2011), (b) the summer (snapshot on 15 July 2011). The blue contour shows the 500 m MLD area. The black dot indicate the LION mooring site. (c) Histogram of relative vorticity during the winter months (December to February) in blue and the summer months (May to July) in red in the box indicated by a dashed line rectangle. (d–f) Same at 500 m depth with different colorbar scale.</figcaption></figure>

## 3.2. Northern Current Frontal Dynamics during a Down-Front Wind Event

Between 15 to 27 January 2011, a glider sampled the NC front back and forth along T03 section (Figure 5a). The two sections showed very different frontal situations following an intense wind event occurring on the 21st of January. For two days, wind stress exceeded 0.4 N m⁻² and heat losses 500 W m⁻². The wind orientation was mostly from the north where the glider sampled the NC, thus being roughly aligned with the NC pathway.

In the first section during calm conditions, the NC front was sampled over a wide area of about 50 km (Figure 5c). This width was larger than the generally observed NC width [39] and could result from the meandering of the NC jet. Isopycnals near the surface were relatively flat except above the continental slope where the NC front separates coastal warm/fresh AW (also influenced by river runoffs [110]) from higher salinity modified AW found offshore. The associated frontal jet was about 10 km wide with geostrophic velocities larger than 0.5 m s⁻¹. An secondary frontal jet was also sampled about 30 km offshore from the slope with a weaker velocity signal. The detachment of a 10-km diameter anticyclonic eddy from the NC can also be observed at about *y* = 55 km.
<figure><img src="image_8.png" /><figcaption>**Figure 5.** (a) Wind and air/sea heat fluxes from ALADIN reanalysis on 21 January 2011. Wind arrows show the direction and intensity in Beaufort wind force scale. Black dots in the west represent glider profiles in the section period with DAC in red. (b) Net surface heat flux (in blue), Ekman equivalent heat flux (in white) and wind stress (in red) interpolated along the glider sampling. The period of intense winds peaking on the 21st of January is highlighted in red, as well as in the glider sections in panel c and d. Glider section of (c) potential temperature and (e) practical salinity with MLD shown by the white line and isopycnals contours in black. The main water masses are labeled (AW: Atlantic Water, LIW: Levantine Intermediate Water, WMDW: Western Mediterranean Deep Water). (g) Glider-based absolute geostrophic currents with isopycnals contours in gray. (d,f,h) Same as (d,e,g) for the second glider crossing of the NC front.</figcaption></figure>

In the second section during and following the wind event (Figure 5d), isopycnals were much steeper. Geostrophic transport increased (wider frontal jet) with no sign of eddies
offshore where the MLD directly dropped below 1000 m. This sharp contrast resulted from the intense surface forcing activating vertical mixing in the convection area. Following the wind event, temperature and salinity interleaving structures were observed in the NC frontal area (y between 20 and 40 km). Below the steepest isopycnals, where Levantine Intermediate Water (LIW) outcrop ($\sigma_{\theta} = 29.05-29.1 \text{ kg m}^{-3}$), along isopycnals up- and downward interleavings of warm/salty and cold/fresh waters can be observed.

The wind event was aligned in the down-front direction with the NC (Figure 5a). As a result, the isopycnals of the front steepened significantly and the stratification generally decreased in the upper 100 m (Figure 6a). In the same area, horizontal buoyancy gradients intensified through the mechanical action of winds (Figure 6b) and the negative baroclinic term in PV became dominant in the vicinity of the NC front. A layer of negative PV (unstable to symmetric instability, SI) was then observed by the glider from the surface down to about 100 m (Figure 6e). This zone of instability coincides with elevated Ekman equivalent heat loss generated by the down-front winds reaching about $1500 \text{ W m}^{-2}$ (estimated from ALADIN reanalysis winds and in situ glider-based buoyancy gradients, see Figure 5b). This flux locally overtook the net surface heat loss, highlighting the important role of wind in the NC front dynamics. The SI layer also corresponded to the area where temperature and salinity interleaving were present (Figure 4d).

<figure><img src="image_9.png" /><figcaption>**Figure 6.** PV decomposition for the glider transect of the right column in Figure 5: (a) stratification term $q_N$ with isopycnal contours in black and (c) baroclinic term $q_{bc}$ term with contours of geostrophic velocities. MLD is shown by the white line. (e) Section of geostrophic PV ($q = q_N + q_{bc}$). The blue contours represent areas prone to symmetric instability. The red lines show where the angle between DAC and glider track exceed 45°. (b,d,f) same as (a,c,d) but evaluated in the numerical model.</figcaption></figure>
The numerical model showed a similar behavior in terms of PV dynamics with the appearance of a SI layer in the upper 100 m. In the model, the thermocline intensity was however slightly less stratified (Figure 6b). This was due to important vertical mixing on the light side of the NC front in the model (slope/shelf region) with MLD reaching almost 300 m compared to about 100 m in observations. Saltier surface waters in the model with practical salinity of 38.3 were observed compared to 38.1 in observations highlighting the importance of coastal processes (e.g, river runoffs). Note also that the front being narrower and shifted offshore by about 20 km, it delayed the front crossing by about a day after the maximum winds. A negative PV layer was still found in the model despite lower Ekman equivalent heat flux at the front of about 500 W m⁻² against 1500 W m⁻² in observations.

### 3.3. *The Basin-Scale PV Response to Typical Winter Mistral Wind Event*

As the numerical model well reproduced the main features of the PV field, further insights in the processes at play and their horizontal distribution can be assessed. The surface forcing at the basin scale at the apex of the wind event on 21 January is shown in Figure 7a,b. The Ekman equivalent heat flux was of similar amplitude as the net surface heat flux, or even locally dominating in frontal areas. The Ekman flux was particularly intense over a large area along the NC front from the Ligurian Sea to the western part of the Gulf of Lion. The wind pattern was indeed a combination of typical northerly Mistral and intense easterly winds corresponding to down-front winds for the NC in the western Gulf of Lion and Ligurian Sea. In the open sea, the fronts orientation mainly influenced by the mesoscale had thus a more variable structure. Consequently, the associated Ekman fluxes were locally intense in negative and positive values over smaller and less organized areas compared to the signal observed along the topographically-steered NC.

<figure><img src="image_10.png" /><figcaption>**Figure 7.** (a) Daily average net surface heat flux and winds from ALADIN re-analysis on 21 January 2011. (b) Ekman equivalent heat flux computed from the modeled surface buoyancy with black arrows proportional to Ekman transport. (c) Modeled geostrophic PV at 50 m depth with black contours corresponding to Ekman equivalent heat flux exceeding 500 W m⁻² (d) Instability regimes at 50 m following Thomas et al. [73]. The dashed line shows the mean NC path, defined as the zero vorticity contour computed from yearly mean currents, used to construct the Hovmöller diagram of Figure 8.</figcaption></figure>
<figure><img src="image_11.png" /><figcaption>**Figure 8.** Hovmöller diagram in the SYMPHONIE model along the NC path (trajectory shown in Figure 7d). The signals are averaged in 10 × 10-km bins, typical of the deformation radius. (a) net surface heat flux, (b) Ekman equivalent heat flux, (c) PV at 50 m depth with black contours where $Q_{Ek} < -500 \text{ W m}^{-2}$, (d) instability regimes where color corresponds to at least 25% of a bin area. (e) Number of days per year along NC with net surface heat flux below −500 W m$^{-2}$ (annual in black, winter DJFM in blue). (f) Same as (e) but for Ekman equivalent heat flux. (g) Number of days per year with negative PV (annual in black, winter DJFM in blue). (h) Number of days per year where different types of instability prevail.</figcaption></figure>

Below the Ekman layer (typically $h_{Ek} = \sqrt{2K_m/|f|} \sim 50 \text{ m}$ with $K_m \sim 0.1 \text{ m}^2 \text{ s}^{-1}$ the turbulent diffusivity), ageostrophic Ekman currents are generally small compared to the geostrophic component. The geostrophic and total PV at 15 m had opposite signs in areas of positive Ekman buoyancy fluxes due to restratifying agesotrophic circulations stabilizing surface fronts. On the other hand, geostrophic PV showed very little difference with the total one deeper at 50 m and is thus be a good proxy to detect symmetrically unstable mesoscale fronts (Figure 7c,d). During the wind event, all areas of negative PV subjected to SI matched strong Ekman heat losses larger than $500 \text{ W m}^{-2}$. SI layers can extend further down to 100 m depth or more (see for instance Figure 6e). At 50 m, they were mostly found in 150–200 km bands along the NC in both the western Gulf of Lion and the Ligurian Sea, but also in the open sea in 10–50 km associated with mesoscale fronts (Figure 7d). Those offshore patches showed a more random distribution and evolved on timescale of a few days typical of the mesoscale (see Video S1 in Supplementary Materials). Strong Ekman heat losses were not always sufficient to trigger SI at 50 m offshore. The destruction of PV at fronts indeed needs to be intense and sustained in time in order to force SI at a mesoscale front. This is then less likely to occur in the open sea compared to the NC aligned with topography.
### 3.4. *Seasonal PV Destruction Along the Northern Current Front*

The numerical model characterized the frontal dynamics of the NC over a complete seasonal cycle. The net surface heat flux followed the seasonal cycle from positive to negative values (−250 to 100 W m⁻² on average) with short periods (1–10 days) of intense heat losses reaching −1000 W m⁻² (Figure 8a). During the simulated year, there were about 10 to 20 days of values below −500 W m⁻² along the NC with more frequent days of strong heat losses in the western Gulf of Lion (Figure 8e). On the contrary, the Ekman equivalent heat flux only depends on the wind intensity and its orientation relative to the NC front. It thus reached intense positive and negative values during all seasons (from −1000 to +1000 W m⁻²) with a slightly calmer period from April to June. During the winter, there were almost 20 days with $Q_{Ek} < -500$ W m⁻² in the western Gulf of Lion and about 10 in the Ligurian Sea. These numbers represent about half of the about 30–40 days found over the year showing the importance of the winter period regarding the wind forcing.

Potential vorticity at 50 m responded to the combined effect of both surface thermodynamic and mechanical forcing. The seasonal restratification promoted positive values by increasing stratification, while PV decreased to low values during the winter at around $10^{-11}$ m s⁻¹. The strong down-front wind events regularly destroyed PV down to at least 50 m causing negative values and generate favorable conditions to SI in ~200-km bands along the NC (Figure 8d). The western part of the Gulf of Lion (around T03) was the main hotspot for SI with almost 40 days of instability per year (Figure 8h) and almost all of them occurring during the winter months (Figure 8g). Farther east in the Ligurian Sea, the NC was also frequently symmetrically unstable during about 20 days (again mostly during the winter). These two regions of the NC were identified as hotspots for SI in the numerical model, in close relationship to the dominant down-front winds prevailing in the area and intense Ekman buoyancy losses.

### 3.5. *Symmetric Instability in the Northern Current Observed by Gliders*

In glider observations, SI events were defined from PV estimates and the balanced Richardson angle (see Section 2.5). Specific SI events were identified where SI was present at 50 m over at least 25 m in the vertical and two successive glider profiles. The characteristics of the corresponding SI fronts in terms of negative PV value, thickness and cross-front width are summarized in Table 2. Bear in mind that profiles with an angle between the DAC and glider track exceeded 45° were excluded (see Section 2.6). On average, the detected SI events had a PV signal of about $-1.6 \times 10^{-11}$ m s⁻¹ over 8.3 km cross-front distance (from 4.7 to 22 km) and 54 m thickness (from 32 to 106 m). The width of SI layers was smaller than the mesoscale of the NC front, even though this submesoscale instability develops within steep fronts whose scale is determined typically by the mesoscale. In line with model detection of SI, it is interesting to note that more events were also detected in the western Gulf of Lion (7/14) than in the Ligurian Sea (1/14), eastern Gulf of Lion (3/14) or the open sea (3/14).

The glider-based estimation of PV reproduced the seasonal cycle (Figure 9a). SI layers at 50 m were only detected during the winters 2011, 2012 and 2013, mostly because of a lack of data during the winters 2009 and 2010. SI events at 50 m were only detected from October to March with a highest occurrence in February (7 events out of 14, Figure 9b). This generally corresponds to the month of deepest MLD and lowest stratification. SI layers at depths shallower than 50 m might still be present during other months, but their detection and characterization by gliders would be less robust due to the sampling so we decided to discard them from the statistics.
<table><caption>**Table 2.** Characteristics of frontal SI events detected by glider: *q*, *L*, *H* correspond to the mean PV, width and thickness; *λ*<sub>SI</sub> refers to the theoretical SI wavelength; *Q*<sub>Ek</sub> is the minimum Ekman equivalent heat flux observed along the glider-track within a 48 h window around each event.</caption><thead><tr><th></th><th>Date</th><th>Location</th><th>q (m<sup>-1</sup> s<sup>-1</sup>)<br>×10<sup>-11</sup></th><th>L (km)</th><th>H (m)</th><th>λ<sub>SI</sub> (km)</th><th>Q<sub>Ek</sub> (W m<sup>-2</sup>)<br>min in 48 h</th></tr></thead><tbody><tr><td>1</td><td>13 December 2010</td><td>NC, GoL west</td><td>-4.7</td><td>8.9</td><td>71</td><td>6.2</td><td>-1810</td></tr><tr><td>2</td><td>22 January 2011</td><td>NC, GoL west</td><td>-2.8</td><td>22.0</td><td>58</td><td>4.0</td><td>-1620</td></tr><tr><td>3</td><td>25 January 2012</td><td>NC, GoL east</td><td>-0.37</td><td>6.2</td><td>32</td><td>0.8</td><td>-30</td></tr><tr><td>4</td><td>2 February 2012</td><td>NC, GoL west</td><td>-0.46</td><td>5.8</td><td>36</td><td>1.0</td><td>-560</td></tr><tr><td>5</td><td>6 February 2012</td><td>NC, GoL west</td><td>-1.8</td><td>12.5</td><td>106</td><td>5.7</td><td>-2930</td></tr><tr><td>6</td><td>13 December 2012</td><td>Open-sea</td><td>-0.79</td><td>4.7</td><td>33</td><td>1.2</td><td>-120</td></tr><tr><td>7</td><td>19 January 2013</td><td>NC, GoL east</td><td>-1.6</td><td>12.1</td><td>75</td><td>3.8</td><td>-110</td></tr><tr><td>8</td><td>19 January 2013</td><td>NC, GoL east</td><td>-1.3</td><td>5.5</td><td>42</td><td>1.9</td><td>-380</td></tr><tr><td>9</td><td>2 February 2013</td><td>NC, GoL west</td><td>-0.66</td><td>4.8</td><td>45</td><td>1.5</td><td>-1870</td></tr><tr><td>10</td><td>7 February 2013</td><td>NC, GoL west</td><td>-0.92</td><td>5.0</td><td>37</td><td>1.5</td><td>-3160</td></tr><tr><td>11</td><td>13 February 2013</td><td>Open-sea</td><td>-5.7</td><td>5.0</td><td>45</td><td>4.7</td><td>-640</td></tr><tr><td>12</td><td>14 February 2013</td><td>NC, GoL west</td><td>-0.8</td><td>10.3</td><td>73</td><td>2.7</td><td>-260</td></tr><tr><td>13</td><td>18 February 2013</td><td>NC, Ligurian</td><td>-3.7</td><td>7.3</td><td>70</td><td>5.3</td><td>-200</td></tr><tr><td>14</td><td>15 March 2013</td><td>Open-sea</td><td>-19</td><td>5.5</td><td>36</td><td>6.6</td><td>-3970</td></tr></tbody><tfoot><tr><td>mean</td><td></td><td></td><td>-1.7 ± 1.1</td><td>8.3 ± 4.8</td><td>54 ± 22</td><td>3.3 ± 2.1</td><td>-1260 ± 1310</td></tr></tfoot></table>

<figure><img src="image_12.png" /><figcaption>**Figure 9.** (a) Time series of PV observed by gliders in the study area. The colored dots correspond to SI events. (b) Monthly histogram of the number of observed SI events. (c) Maps of SI events with gray dots showing all glider profiles. SI events are represented by a dot color-coded by the observed mean negative PV at 50 m and whose size is proportional to the depth of the SI layer. The blue (resp. red) contour shows where the number of days prone to SI in the model exceeds 10 (resp. 20 days). Light gray contours represent isobaths every 500-m.</figcaption></figure>
Regarding the characteristics of symmetric instability in unstable fronts, Thomas and Lee [102] inferred a wavelength in relationship with the negative PV signal $q$ and thickness $H$ of the SI layers: $\lambda_{SI} = 4H\sqrt{-fq}/f^2$. On average, $\lambda_{SI}$ was about 3.3 km. Note that length is smaller than the cross-front scale of SI layers of 8.3 km on average. Even if this number might be underestimated given the short scales in space and time of negative PV signals, it corresponds on average to about 2 wavelengths. This is more or less what can be observed in Figure 5 corresponding to event 3 in Table 2 showing SI extending over about 20 km and 60 m deep with two unstable wavelengths revealed by temperature and salinity interleavings.

### 3.6. Wind Forcing of Symmetric Instability

Northerly dry Mistral events seem to dominate the wind forcing of the northwestern Mediterranean Sea, as exemplified by the first empirical orthogonal function (EOF) explaining about 66.8% of the wind stress variability (Figure 10). The second EOF corresponding to 13.5% of the variance represents easterly winds in the Ligurian Sea combined with westerly and northerly winds in the western GoL. The first two EOFs are both situations of down-front winds along the NC in the two hotspots of SI identified in the model (western Gulf of Lion and Ligurian Sea).

<figure><img src="image_13.png" /><figcaption>**Figure 10.** Principal component analysis of daily wind stress from CNRM-ALADIN reanalysis (2009–2013). (a) First EOF representing 66.8% of the signal variance. Dots represent the spatial distribution of SI events with blue dots during period when pc exceeds one standard deviation. The red arrows shows the yearly average surface circulation from the model. (b) Time series of the eigenvalues for the first EOF. The dashed blue lines are drawn at $\pm$ one std. Dots show pc values during SI events (in blue for those exceeding one std). (c,d) same as (a,b) but for the second EOF (13.5% of the variance).</figcaption></figure>

In terms of SI events, the majority of them occurred during strong Mistral or easterly winds: 10/14 (reps. 8/14) of the detected SI events have a principal component for the first (resp. second) EOF larger than one standard deviation (std). In a general manner, SI events were closely linked to the dominant wind patterns in the area and presence of topographically-steered front of the Northern Current. Over the study period, 25% of the year the wind were favorable to PV destruction in the two hotspots for SI along the NC
(western Gulf of Lion and Ligurian Sea) with one of the first two EOFs exceeding one std. This was also confirmed by the strongly negative Ekman equivalent heat flux reported by gliders during the sampling SI events. The minimum Ekman buoyancy flux found within a 48 h window was on average −1260 W m⁻² and 8/14 events had values dropping below −500 W m⁻². At least half of the SI events were thus associated with intense down-front wind events. Observing strong negative Ekman equivalent heat flux requires the front sampling to be synchronous with the down-front winds, which might not be the case for every SI events.

### 3.7. Vertical Circulations and Biogeochemical Implications

Symmetric instability forced by down-front winds can generate important vertical velocities. This is also suggested by glider temperature and salinity sections (see Figures 5 and 11 showing another SI layer triggered by a strong Mistral event in February 2013). A glider crossed the NC front along the T03 section at the time of maximum wind. The wind stress reached 1.5 N m⁻², while net surface heat losses peaked at about 1200 W m⁻² and Ekman equivalent buoyancy fluxes reached almost −2000 W m⁻² (Figure 11b). A SI layer of about 50 m in thickness and 5 km in cross-front direction was captured at the NC front by the glider. The temperature and salinity sections exhibit a clear pattern of alternated vertical displacements below the mixed layer depth (Figure 11g,h).

Cold/fresh water was subducted down to 200 to 400 m below the NC front in two cells transporting high oxygen and chlorophyll waters (Figure 11e,f). In the mixed layer, especially when the winds were strong and the mixing active, phytoplankton was diluted on the vertical by convective convective plumes whose amplitude can reach 0 (1000 m days⁻¹) in the deep convection area [97]. Following Niewiadomska et al. [111] who described the optical properties of subduction along the one of the first glider deployment along the MOOSE T00 endurance line, the sign of frontal ageostrophic secondary circulation can easily be identified. At short timescale, vertical circulations modulate the chlorophyll distribution as a passive tracer with patches transferred at depth and low-chlorophyll water lifted toward the surface. Upward movements could also inject nutrients into the sunlit layer and fuel phytoplankton growth. Higher chlorophyll concentrations could well be observed at the surface near the unstable front. However, it is hard to disentangle the dynamical processes that can affect chlorophyll concentration (dilution, restratification, lateral advection, ...) from local production.

Phytoplankton concentration should decay below the sunlit layer in the absence of fluxes that can sustain the phytoplankton population. For a surface chlorophyll concentration of 0.5–1.5 mg m⁻³ observed at the front, the euphotic layer, where only 1% of the surface photosynthetic available radiation remains, can be estimated to 30–45 m following Morel and Berthon [112]. In the deepest patch of chlorophyll observed at around 200–400 m, phytoplankton cells were thus well below the euphotic layer. This signal could not be caused by local production but should result from important vertical displacements. Considering a survival timescale of phytoplankton in absence of light of about 1–10 days, this would indicate intense vertical circulations of O (10–100 m day⁻¹) along the NC front below the unstratified mixed layer.

Numerical model suggested that frontal regions unstable to SI are organized in O (100 km) bands along the NC front. Considering vertical displacements *w*′ of O (10–100 m day⁻¹) roughly estimated from the subduction of phytoplankton below the mixed layer in O (1 km) cross-front patches, this would represent a subduction/obduction rate of 0.01–0.1 Sv along the whole unstable frontal area. This would apply only when SI is active (i.e., 20–40 days in the two hotspots of the western Gulf of Lion and Ligurian Sea). Along isopycnal temperature and salinity anomalies *T*′/*S*′ associated with interleaving could also be observed with the order of magnitude of 0.2 °C/0.05 (see for instance Figure 5 or Figure 11). The sign of the anomalies was correlated to the vertical velocities, so all contribute to an upward heat flux of typically *F*<sub>Q</sub> = *ρC*<sub>p0</sub> *T*<sub>0</sub> *T*′ > ~100–1000 W m⁻². This rough scaling shows that vertical heat fluxes generated by SI in the NC could be
comparable to fluxes due to surface heat flux or mechanical effect of winds. This process thus needs a particular attention in order to close the heat/salt budgets of the NC. Similarly, the oxygen budget of intermediate layers might also be importantly affected by SI, as well as vertical fluxes of nutrients and carbon. Bear in mind that other (sub)mesoscale processes, for instance slower but more persistent baroclinic mixed layer instabilities [113], might also played an important role in the front heat, salt and momentum balance.

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

**Figure 11.** Glider sampling of SI event 9 in Table 2. (a) Net surface heat flux and winds on February 3rd from CNRM-ALADIN reanalysis. (b) Net surface heat flux (i blue), Ekman equivalent heat flux (in white) and wind stress (in red) along the glider track. The red bar shows the peak of the wind event on February 3. Glider sections of (c) potential temperature, (d) practical salinity, (e) dissolved oxygen (uncalibrated), (f) chlorophyll-a fluorescence (quality-controlled with satellite, see Section 2.1), (g) potential vorticity with contours of SI layer in blue, (h) instability categories. On panels (c–h), the gray contours are isopycnals and the white line shows the mixed layer depth.
## 4. Discussion

The topographically-steered path of the NC in the northwestern Mediterranean Sea enables sustained destruction of PV over large areas regularly subjected to down-front winds. This current system also interact with deep convection where vertical mixing reaching depths from 500 to 2500 m intensifies intensity and depth of horizontal density gradients between the NC, an area subjected to lateral advection, and the preconditioned convection area offshore. This interaction contributes to important PV destruction by weakening the water column stratification. The main winds (Mistral and Tramontane in the Gulf of Lion, and easterlies in the Ligurian Sea) foster down-front winds in 100–200 km bands along the NC and symmetric instability especially in two hotspots: the western Gulf of Lion and the Ligurian Sea. Whether this instability remains important in the open sea, and especially along the North Balearic Front, needs to be further investigated. Previous studies reported important effect of down-front winds in the Atlantic and the Southern Ocean [69,70,82]. However, these regions are characterized by larger deformation radius (50–100 km). In the northwestern Mediterranean Sea, the mesoscale is generally smaller (10–30 km) and evolves quickly (1–5 days). This might be a limiting factor to sustain destruction of PV outside of the topographically-steered circulation. A Lagrangian water parcel caught in a front meandering at the mesoscale would indeed stay aligned in the down-front direction during period of time possibly shorter than the timescale of the wind event. The quantification of the relative importance of other sources of vertical circulations (e.g., frontogenesis [43,46]) and their seasonal dependence and forcing deserves to be studied specifically.

The horizontal and vertical resolution and viscosity are key ingredients for a numerical model to be able to resolve SI, as shown by Bachman and Taylor [114] in an idealized front. In a realistic model, the parameterized vertical mixing scheme makes it however harder to disentangle the relative contribution of restratification by SI from the increasing vertical diffusivity. Following on the example of the NC front after the intense down-front event on 21 January 2011 showing deep and extended negative PV layers, a closer inspection at hourly outputs showed how unbalanced motions rapidly ceased after restoring the front in geostrophic balance within about a day after the peak winds (not shown). Within the SI layer, important positive and negative vertical velocities also emerged in the model (up to ~100 m per day) suggesting that the SYMPHONIE numerical model at 1 km resolution could at least marginally resolve some submesoscale features of the front. However, a horizontal resolution of 0 (100 m) seems required to better resolve SI because of the isopycnal slope in the Northern Current can reach value close to 0.1. The intensity and thickness of the negative PV layers can indeed lead to larger wavelengths of the instability [102], easier to be resolved numerically.

Gliders were able to characterize symmetrically unstable fronts in the northwestern Mediterranean Sea, adding to regions such as the Equatorial Pacific [81], the Southern Ocean [69,70,83], the Atlantic Ocean [82,84] and the Gulf Stream [78]. The key to capture negative PV was to sample the front during or shortly after a down-front wind event. Glider track orthogonal to the frontal jet is also important, although cross-stream coordinate system can correct the geometrical bias in horizontal gradients observed by gliders [78,85]. The constrained path of topographically-steered currents (like the Northern Current) helped gliders sample perpendicular to the flow. A sustained effort in deployment during several year (as part of the MOOSE observatory [86]) was also important in order to cover the seasonality of the instabilities and to draw statistics from a significant number of events.

Observations further showed the impacts of the vertical exchanges generated by SI on physical (temperature, salinity) and biogeochemical (oxygen, chlorophyll-a fluorescence) properties. The impact was important below the mixed layer at intermediate depths (200–600 m) where one finds the warm, salty and relatively poor in oxygen Levantine Intermediate Water. The subduction of cold, fresh and oxygenated patches below unstable fronts could thus have an important impact on the heat, salt and oxygen budgets of the Northern Current. In addition SI being a recurrent feature in the NC during the winter, it
could importantly impact the vertical fluxes of nutrients, boost phytoplankton growth and drive the export of carbon. The present work could be extended pending a specific and important work on the harmonization of oxygen and optical measurements (fluorescence and backscatter) from the MOOSE gliders. Finally, among the processes impacting vertical fluxes, the enhanced turbulence during SI events numerically predicted by Taylor and Ferrari [75] could also be resolved using gliders equipped with high frequency shear probes [115].

High resolution observations and simulations were shown to be complementary. While glider observed episodic signatures of SI at fronts, the model reproducing those signals could be used to map the processes spatially and extend the analysis over a full seasonal cycle. Based on these results, a similar synergy could be applied with a coupled physics-biogeochemistry numerical model to tackle the open questions about the impact of SI along the NC for the ventilation of intermediate layers and carbon sequestration.

More generally, the regular instances of symmetric instability during the winter around the deep convection region of the northwestern Mediterranean Sea raises the question of whether it is also the case in other deep convection sites (e.g., the Labrador, Irminger, Greenland and Weddel Seas). The role of submesoscale dynamics on deep convection has been shown to be important in the Labrador Sea highlighting the need to account for their effects in a large-scale earth system model [62]. All those sites are regions of intense dry and cold winds occurring in bursts. The interaction of those winds with the topographically-steered circulation could reveal hotspots of SI. The large buoyancy fluxes generated at the front by mechanical effect near deep convection regions further promote restratification or instabilities (SI, gravitational or mixed). Furthermore, given the scale of the associated heat fluxes, wind-driven submesoscale likely influences the intensity and position of deep convection at a larger scale. Numerical experiment including parametrization of SI (e.g., Bachman et al. [74]) with regional ocean model could also be used in order to evaluate the future evolution of deep convection in a more realistic approach.

**Supplementary Materials:** The following are available online at www.mdpi.com/xxx/s1, Video S1: Seasonal evolution of Surface buoyancy fluxes, Potential Vorticity and instability regimes at 50 m.

**Author Contributions:** A.B., P.T., L.P., V.T. conceptualized the analysis; A.B., P.T. developed the methodology; A.B. did the formal analysis; P.D., C.E., P.M. developed the numerical simulation; A.B. did the writing and original draft preparation; all authors contributed to reviewing and editing the manuscript; P.T., L.M. supervised A.B during their PhD and acquired funding of this research. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by the French "Chantier Méditranée" MISTRALS (HyMeX and MERMeX programs) and the French long-term Mediterranean Observatory (MOOSE), apart of Research Infrastructures (SOERE/AllEnvi-SNO/INSU, IR-ILICO-EMSO-ERIC), as well as by the EU projects FP7 GROOM (Grant Agreement No. 284321), FP7 PERSEUS (Grant Agreement No. 287600) and FP7 JERICO (Grant Agreement No. 262584) and the COST Action ES0904 "EGO" (Everyone's Gliding Observatories). The authors acknowledge the support of the HPC resources from CALMIP (Grants P1325, P09115) and from GENCI and CINES (Grand Equipement National de Calcul Intensif, project A0040110088).

**Data Availability Statement:** The data collected by the MOOSE network are freely accessible: in particular, the data used here collected by the gliders (https://doi.org/10.17882/52027 (accessed on 16 March 2021), [116]) and the LION mooring line (https://doi.org/10.17882/44411 (accessed on 16 March 2021), [87]). The model outputs are available upon request to the corresponding author.

**Acknowledgments:** We thank the staff of the French National Gliders Facility at DT-INSU (CNRS/Ifremer) for successful gliders missions carried out over the years in the framework of MOOSE, as well as the intensive deployments during this 2012–2013 DEWEX experiment. Captains and crew of R/V *Le Tethys II*, *Le Suroit*, as well as all scientists and technicians who participated to MOOSE-GE cruises maintaining the LION and DYFAMED moored time series are also warmly thanked.
**Conflicts of Interest:** The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

## References

1. Houpert, L.; Durrieu de Madron, X.; Testor, P.; Bosse, A.; D'Ortenzio, F.; Bouin, M.N.; Dausse, D.; Le Goff, H.; Kunesch, S.; Labaste, M.; et al. Observations of open-ocean deep convection in the northwestern Mediterranean Sea: Seasonal and interannual variability of mixing and deep water masses for the 2007-2013 Period. *J. Geophys. Res. Ocean.* **2016**, *121*, 8139–8171. [CrossRef]

2. Marshall, J.; Schott, F. Open-ocean convection: Observations, theory, and models. *Rev. Geophys.* **1999**, *37*, 1–64. [CrossRef]

3. Send, U.; Marshall, J. Integral Effects of Deep Convection. *J. Phys. Oceanogr.* **1995**, *25*, 855–872. [CrossRef]

4. Koelling, J.; Wallace, D.W.R.; Send, U.; Karstensen, J. Intense oceanic uptake of oxygen during 2014–2015 winter convection in the Labrador Sea. *Geophys. Res. Lett.* **2017**, *44*, 7855–7864. [CrossRef]

5. Coppola, L.; Prieur, L.; Taupier-Letage, I.; Estournel, C.; Testor, P.; Lefèvre, D.; Belamari, S.; LeReste, S.; Taillandier, V. Observation of oxygen ventilation into deep waters through targeted deployment of multiple Argo-O2 floats in the north-western Mediterranean Sea in 2013. *J. Geophys. Res. Ocean.* **2017**, *122*, 6325–6341. [CrossRef]

6. Ulses, C.; Estournel, C.; Fourrier, M.; Coppola, L.; Kessouri, F.; Lefèvre, D.; Marsaleix, P. Oxygen budget of the north-western Mediterranean deep-convection region. *Biogeosciences* **2021**, *18*, 937–960. [CrossRef]

7. D'Ortenzio, F.; Lavigne, H.; Besson, F.; Claustre, H.; Coppola, L.; Garcia, N.; Laës-Huon, A.; Reste, S.L.; Malardé, D.; Migon, C.; et al. Observing mixed layer depth, nitrate and chlorophyll concentrations in the northwestern Mediterranean: A combined satellite and NO3 profiling floats experiment. *Geophys. Res. Lett.* **2014**, *41*, 6443–6451. [CrossRef]

8. Lacour, L.; Claustre, H.; Prieur, L.; D'Ortenzio, F. Phytoplankton biomass cycles in the North Atlantic subpolar gyre: A similar mechanism for two different blooms in the Labrador Sea. *Geophys. Res. Lett.* **2015**, *42*, 5403–5410. [CrossRef]

9. Mayot, N.; D'Ortenzio, F.; Taillandier, V.; Prieur, L.; de Fommervault, O.P.; Claustre, H.; Bosse, A.; Testor, P.; Conan, P. Physical and Biogeochemical Controls of the Phytoplankton Blooms in North Western Mediterranean Sea: A Multiplatform Approach Over a Complete Annual Cycle (2012–2013 DEWEX Experiment). *J. Geophys. Res. Ocean.* **2017**, *122*, 9999–10019. [CrossRef]

10. Fröb, F.; Olsen, A.; Våge, K.; Moore, G.W.K.; Yashayaev, I.; Jeansson, E.; Rajasakaren, B. Irminger Sea deep convection injects oxygen and anthropogenic carbon to the ocean interior. *Nat. Commun.* **2016**, *7*, 13244. [CrossRef]

11. Spall, M.A. Dynamics of Downwelling in an Eddy-Resolving Convective Basin. *J. Phys. Oceanogr.* **2010**, *40*, 2341–2347. [CrossRef]

12. Waldman, R.; Brüggemann, N.; Bosse, A.; Spall, M.; Somot, S.; Sevault, F. Overturning the Mediterranean Thermohaline Circulation. *Geophys. Res. Lett.* **2018**, *45*, 8407–8415. [CrossRef]

13. Somot, S.; Houpert, L.; Sevault, F.; Testor, P.; Bosse, A.; Taupier-Letage, I.; Bouin, M.N.; Waldman, R.; Cassou, C.; Sanchez-Gomez, E.; et al. Characterizing, modelling and understanding the climate variability of the deep water formation in the North-Western Mediterranean Sea. *Clim. Dyn.* **2018**, *51*, 1179–1210. [CrossRef]

14. Schroeder, K.; Ribotti, A.; Borghini, M.; Sorgente, R.; Perilli, A.; Gasparini, G.P. An extensive western Mediterranean deep water renewal between 2004 and 2006. *Geophys. Res. Lett.* **2008**, *35*, 1–7. [CrossRef]

15. Schroeder, K.; Chiggiato, J.; Bryden, H.L.; Borghini, M.; Ben Ismail, S. Abrupt climate shift in the Western Mediterranean Sea. *Sci. Rep.* **2016**, *6*, 23009. [CrossRef]

16. Margirier, F.; Testor, P.; Heslop, E.; Mallil, K.; Bosse, A.; Houpert, L.; Mortier, L.; Bouin, M.N.; Coppola, L.; D'ortenzio, F.; et al. Abrupt warming and salinification of intermediate waters interplays with decline of deep convection in the Northwestern Mediterranean Sea. *Sci. Rep.* **2020**, *10*, 1–11. [CrossRef]

17. Somot, S.; Sevault, F.; Déqué, M. Transient climate change scenario simulation of the Mediterranean Sea for the twenty-first century using a high-resolution ocean circulation model. *Clim. Dyn.* **2006**, *27*, 851–879. [CrossRef]

18. Somot, S.; Sevault, F.; Déqué, M.; Crépon, M. 21st century climate change scenario for the Mediterranean using a coupled atmosphere-ocean regional climate model. *Glob. Planet. Chang.* **2008**, *63*, 112–126. [CrossRef]

19. Herrmann, M.; Estournel, C.; Adloff, F.; Diaz, F. Impact of climate change on the northwestern Mediterranean Sea pelagic planktonic ecosystem and associated carbon cycle. *J. Geophys. Res. Ocean.* **2014**, *119*, 5815–5836. [CrossRef]

20. Macias, D.; Garcia-Gorriz, E.; Stips, A. Deep winter convection and phytoplankton dynamics in the NW Mediterranean Sea under present climate and future (horizon 2030) scenarios. *Sci. Rep.* **2018**, *8*, 6626. [CrossRef]

21. MEDOC Group. Observation of Formation of Deep Water in the Mediterranean Sea, 1969. *Nature* **1970**, *225*, 1037–1040. [CrossRef]

22. Testor, P.; Bosse, A.; Houpert, L.; Margirier, F.; Mortier, L.; Legoff, H.; Dausse, D.; Labaste, M.; Karstensen, J.; Hayes, D.; et al. Multiscale Observations of Deep Convection in the Northwestern Mediterranean Sea During Winter 2012–2013 Using Multiple Platforms. *J. Geophys. Res. Ocean.* **2018**, *123*, 1745–1776. [CrossRef]

23. Leaman, K.D.; Schott, F.A. Hydrographic Structure of the Convection Regime in the Gulf of Lions: Winter 1987. *J. Phys. Oceanogr.* **1991**, *21*, 575–598. [CrossRef]

24. Gascard, J.C. Mediterranean deep water formation baroclinic instability and oceanic eddies. *Oceanol. Acta* **1978**, *1*, 315–330.

25. Legg, S.; Marshall, J.C. A Heton Model of the Spreading Phase of Open-Ocean Deep Convection. *J. Phys. Oceanogr.* **1993**, *23*, 1040–1056. [CrossRef]

26. Visbeck, M.; Marshall, J.; Jones, H. Dynamics of Isolated Convective Regions in the Ocean. *J. Phys. Oceanogr.* **1996**. [CrossRef]

27. Jones, H.; Marshall, J. Restratification after Deep Convection. *J. Phys. Oceanogr.* **1997**, *27*, 2276–2287. [CrossRef]
28. McWilliams, J.C. Submesoscale, coherent vortices in the ocean. *Rev. Geophys.* **1985**, *23*, 165. [CrossRef]

29. Lilly, J.M.; Rhines, P.B. Coherent Eddies in the Labrador Sea Observed from a Mooring. *J. Phys. Oceanogr.* **2002**, *32*, 585–598.
[CrossRef]

30. Testor, P.; Gascard, J.C. Large-Scale Spreading of Deep Waters in the Western Mediterranean Sea by Submesoscale Coherent Eddies. *J. Phys. Oceanogr.* **2003**, *33*, 75–87. [CrossRef]

31. Testor, P.; Gascard, J.C.C. Post-convection spreading phase in the Northwestern Mediterranean Sea. *Deep. Sea Res. Part I* **2006**, *53*, 869–893. [CrossRef]

32. Bower, A.S.; Hendry, R.M.; Amrhein, D.E.; Lilly, J.M. Direct observations of formation and propagation of subpolar eddies into the subtropical North Atlantic. *Deep Sea Res. II* **2013**, *85*, 15–41. [CrossRef]

33. Zhang, Z.; Li, P.; Xu, L.; Li, C.; Zhao, W.; Tian, J.; Qu, T. Subthermocline eddies observed by rapid-sampling Argo floats in the subtropical northwestern Pacific Ocean in Spring 2014. *Geophys. Res. Lett.* **2015**, *42*, 6438–6445. [CrossRef]

34. Bosse, A.; Testor, P.; Houpert, L.; Damien, P.; Prieur, L.; Hayes, D.; Taillandier, V.; Durrieu de Madron, X.; d'Ortenzio, F.; Coppola, L.; et al. Scales and dynamics of Submesoscale Coherent Vortices formed by deep convection in the northwestern Mediterranean Sea. *J. Geophys. Res. Ocean.* **2016**, *121*, 7716–7742. [CrossRef]

35. Bosse, A.; Testor, P.; Mayot, N.; Prieur, L.; D'Ortenzio, F.; Mortier, L.; Le Goff, H.; Gourcuff, C.; Coppola, L.; Lavigne, H.; et al. A submesoscale coherent vortex in the Ligurian Sea: From dynamical barriers to biological implications. *J. Geophys. Res. Ocean.* **2017**, *122*, 6196–6217. [CrossRef]

36. Damien, P.; Bosse, A.; Testor, P.; Marsaleix, P.; Estournel, C. Modeling Postconvective Submesoscale Coherent Vortices in the Northwestern Mediterranean Sea. *J. Geophys. Res. Ocean.* **2017**, *122*, 9937–9961. [CrossRef]

37. Li, C.; Zhang, Z.; Zhao, W.; Tian, J. A statistical study on the subthermocline submesoscale eddies in the northwestern Pacific Ocean based on Argo data. *J. Geophys. Res. Ocean.* **2017**, *122*, 3586–3598. [CrossRef]

38. Millot, C. Circulation in the Western Mediterranean Sea. *J. Mar. Syst.* **1999**, *20*, 423–442. [CrossRef]

39. Prieur, L.; d'ortenzio, F.; Taillandier, V.; Testor, P. Physical Oceanography of the Ligurian Sea. In *The Mediterranean Sea in the Era of Global Change 1: 30 Years of Multidisciplinary Study of the Ligurian Sea*; ISTE Ltd: London, UK, 2020; pp. 49–78. [CrossRef]

40. Waldman, R.; Somot, S.; Herrmann, M.; Sevault, F.; Isachsen, P.E. On the chaotic variability of deep convection in the Mediterranean Sea. *Geophys. Res. Lett.* **2018**, *45*, 2433–2443. [CrossRef]

41. Crépon, M.; Wald, L.; Monget, J.M. Low-Frequency Waves in the Ligurian Sea During December 1977. *J. Geophys. Res.* **1982**, *87*, 595–600. [CrossRef]

42. Piterbarg, L.; Taillandier, V.; Griffa, A. Investigating frontal variability from repeated glider transects in the Ligurian Current (North West Mediterranean Sea). *J. Mar. Syst.* **2013**. [CrossRef]

43. Capet, X.J.; McWilliams, J.C.; Molemaker, M.J.; Shchepetkin, A.F. Mesoscale to Submesoscale Transition in the California Current System. Part II: Frontal Processes. *J. Phys. Oceanogr.* **2008**, *38*, 44–64. [CrossRef]

44. McWilliams, J.C. Submesoscale currents in the ocean. *Proc. R. Soc. Lond. Math. Phys. Eng. Sci.* **2016**, *472*, 20160117. [CrossRef]
[PubMed]

45. Pollard, R.T.; Regier, L.A. Vorticity and Vertical Circulation at an Ocean Front. *J. Phys. Oceanogr.* **1992**, *22*, 609–625. [CrossRef]

46. Spall, M.A. Frontogenesis, subduction, and cross-front exchange at upper ocean fronts. *J. Geophys. Res.* **1995**, *100*, 2543. [CrossRef]

47. Sanz, E.P.; Viúdez, A. Diagnosing Mesoscale Vertical Motion from Horizontal Velocity and Density Data. *J. Phys. Oceanogr.* **2005**, *35*, 1744–1762. [CrossRef]

48. Ruiz, S.; Pascual, A.; Garau, B.; Pujol, I.; Tintoré, J. Vertical motion in the upper ocean from glider and altimetry data. *Geophys. Res. Lett.* **2009**, *36*, 1–6. [CrossRef]

49. Zhang, Z.; Zhang, X.; Qiu, B.; Zhao, W.; Zhou, C.; Huang, X.; Tian, J. Submesoscale currents in the subtropical upper ocean observed by long-term high-resolution mooring arrays. *J. Phys. Oceanogr.* **2021**, *51*, 187–206. [CrossRef]

50. D'Asaro, E.; Lee, C.M.; Rainville, L.; Harcourt, R.R.; Thomas, L.N. Enhanced Turbulence and Energy Dissipation at Ocean Front. *Science* **2010**, *332*, 318–322. [CrossRef]

51. Johnston, T.M.S.; Rudnick, D.L.; Pallàs-Sanz, E. Elevated mixing at a front. *J. Geophys. Res.* **2011**, *116*, 1–14. [CrossRef]

52. Su, Z.; Wang, J.; Klein, P.; Thompson, A.F.; Menemenlis, D. Ocean submesoscales as a key component of the global heat budget. *Nat. Commun.* **2018**, *9*, 775. [CrossRef]

53. Wenegrat, J.O.; Thomas, L.N.; Sundermeyer, M.A.; Taylor, J.R.; D'Asaro, E.A.; Klymak, J.M.; Shearman, R.K.; Lee, C.M. Enhanced mixing across the gyre boundary at the Gulf Stream front. *Proc. Natl. Acad. Sci. USA* **2020**, *117*, 17607–17614. [CrossRef] [PubMed]

54. Taylor, J.R.; Ferrari, R. Ocean fronts trigger high latitude phytoplankton blooms. *Geophys. Res. Lett.* **2011**, *38*, 1–5. [CrossRef]

55. Lévy, M.; Ferrari, R.; Franks, P.J.S.; Martin, A.P.; Rivière, P. Bringing physics to life at the submesoscale. *Geophys. Res. Lett.* **2012**, *39*, 1–13. [CrossRef]

56. Lévy, M.; Franks, P.J.S.; Smith, K.S. The role of submesoscale currents in structuring marine ecosystems. *Nat. Commun.* **2018**, *9*, 4758. [CrossRef] [PubMed]

57. Ruiz, S.; Claret, M.; Pascual, A.; Olita, A.; Troupin, C.; Capet, A.; Tovar-Sánchez, A.; Allen, J.; Poulain, P.M.; Tintoré, J.; et al. Effects of oceanic mesoscale and submesoscale frontal processes on the vertical transport of phytoplankton. *J. Geophys. Res. Ocean.* **2019**, *124*, 5999–6014. [CrossRef]

58. Luyten, J. The ventilated thermocline. *J. Phys. Oceanogr.* **1983**, *13*, 292–309. [CrossRef]
59. Maze, G.; Deshayes, J.; Marshall, J.C.; Tréguier, A.M.; Chronis, A.; Vollmer, L. Surface vertical PV fluxes and subtropical mode water formation in an eddy-resolving numerical simulation. *Deep. Sea Res. Part II Top. Stud. Oceanogr.* **2013**, *91*, 128–138. [CrossRef]

60. McCartney, M.S. The subtropical recirculation of mode waters. *J. Mar. Res.* **1982**, *40*, 427–464.

61. Giordani, H.; Lebeaupin-Brossier, C.; Léger, F.; Caniaux, G. A PV-approach for dense water formation along fronts: Application to the Northwestern Mediterranean. *J. Geophys. Res. Ocean*. **2017**, *122*, 995–1015. [CrossRef]

62. Tagklis, F.; Bracco, A.; Ito, T.; Castelao, R. Submesoscale modulation of deep water formation in the Labrador Sea. *Sci. Rep.* **2020**, *10*, 1–13.

63. Bopp, L.; Lévy, M.; Resplandy, L.; Sallée, J.B. Pathways of anthropogenic carbon subduction in the global ocean. *Geophys. Res. Lett.* **2015**, *42*, 6416–6423. [CrossRef]

64. Ertel, H. Ein neuer hydrodynamischer Wirbelsatz. *Meteorol. Z.* **1942**, *59*, 277–281.

65. Haynes, P.H.; McIntyre, M.E. On the Evolution of Vorticity and Potential Vorticity in the Presence of Diabatic Heating and Frictional or Other Forces. *J. Atmos. Sci.* **1987**, *44*, 828–841. [CrossRef]

66. Yoshikawa, Y.; Akitomo, K.; Awaji, T. Formation process of intermediate water in baroclinic current under cooling. *J. Geophys. Res.* **2001**, *106*, 1033. [CrossRef]

67. Thomas, L.N. Intensification of Ocean Fronts by Down-Front Winds. *J. Phys. Oceanogr.* **2005**, *35*, 1086–1102. [CrossRef]

68. Thomas, L.N. Destruction of Potential Vorticity by Winds. *J. Phys. Oceanogr.* **2005**, *35*, 2457–2466. [CrossRef]

69. du Plessis, M.; Swart, S.; Ansorge, I.J.; Mahadevan, A. Submesoscale processes promote seasonal restratification in the Subantarctic Ocean. *J. Geophys. Res. Ocean*. **2017**, *122*, 2960–2975. [CrossRef]

70. Viglione, G.A.; Thompson, A.F.; Flexas, M.M.; Sprintall, J.; Swart, S. Abrupt Transitions in Submesoscale Structure in Southern Drake Passage: Glider Observations and Model Results. *J. Phys. Oceanogr.* **2018**, *48*, 2011–2027. [CrossRef]

71. Seyfried, L.; Estournel, C.; Marsaleix, P.; Richard, E. Dynamics of the North Balearic Front during an autumn tramontane and mistral storm: Air–sea coupling processes and stratification budget diagnostic. *Ocean. Sci.* **2019**, *15*, 179–198. [CrossRef]

72. Hoskins, B.J. The role of potential vorticity in symmetric stability and instability. *Q. J. R. Meteorol. Soc.* **1974**, *100*, 480–482. [CrossRef]

73. Thomas, L.N.; Taylor, J.R.; Ferrari, R.; Joyce, T.M. Deep-Sea Research II Symmetric instability in the Gulf Stream. *Deep. Sea Res. Part II* **2013**, 1–15. [CrossRef]

74. Bachman, S.D.; Fox-Kemper, B.; Taylor, J.R.; Thomas, L.N. Parameterization of frontal symmetric instabilities. I: Theory for resolved fronts. *Ocean. Model.* **2017**, *109*, 72–95. [CrossRef]

75. Taylor, J.R.; Ferrari, R. On the equilibration of a symmetrically unstable front via a secondary shear instability. *J. Fluid Mech.* **2009**, *622*, 103–113. [CrossRef]

76. Molemaker, M.J.; McWilliams, J.C.; Dewar, W.K. Submesoscale Instability and Generation of Mesoscale Anticyclones near a Separation of the California Undercurrent. *J. Phys. Oceanogr.* **2015**, *45*, 613–629. [CrossRef]

77. Gula, J.; Molemaker, M.J.; McWilliams, J.C. Topographic generation of submesoscale centrifugal instability and energy dissipation. *Nat. Commun.* **2016**, *7*, 1–7. [CrossRef] [PubMed]

78. Todd, R.E.; Owens, W.B.; Rudnick, D.L. Potential Vorticity Structure in the North Atlantic Western Boundary Current from Underwater Glider Observations. *J. Phys. Oceanogr.* **2016**, *46*, 327–348. [CrossRef]

79. Yoshikawa, Y.; Lee, C.M.; Thomas, L.N. The Subpolar Front of the Japan/East Sea. Part III: Competing Roles of Frontal Dynamics and Atmospheric Forcing in Driving Ageostrophic Vertical Circulation and Subduction. *J. Phys. Oceanogr.* **2012**, *42*, 991–1011. [CrossRef]

80. Testor, P.; de Young, B.; Rudnick, D.; Glenn, S.; Hayes, D.; Lee, C.M.; Pattiaratchi, C.; Hill, K.; Heslop, E.; Turpin, V.; et al. OceanGliders: A component of the integrated GOOS. *Front. Mar. Sci.* **2019**, *6*, 422. [CrossRef]

81. Pietri, A.; Testor, P.; Echevin, V.; Chaigneau, A.; Mortier, L.; Eldin, G.; Grados, C. Finescale Vertical Structure of the Upwelling System off Southern Peru as Observed from Glider Data. *J. Phys. Oceanogr.* **2013**, *43*, 631–646. [CrossRef]

82. Thompson, A.F.; Lazar, A.; Buckingham, C.; Naveira Garabato, A.C.; Damerell, G.M.; Heywood, K.J. Open-ocean submesoscale motions: A full seasonal cycle of mixed layer instabilities from gliders. *J. Phys. Oceanogr.* **2016**, *46*, 1285–1307. [CrossRef]

83. Ruan, X.; Thompson, A.F.; Flexas, M.M.; Sprintall, J. Contribution of topographically generated submesoscale turbulence to Southern Ocean overturning. *Nat. Geosci.* **2017**, *10*, 840–845. [CrossRef]

84. Yu, X.; Naveira Garabato, A.C.; Martin, A.P.; Gwyn Evans, D.; Su, Z. Wind-forced symmetric instability at a transient mid-ocean front. *Geophys. Res. Lett.* **2019**, *46*, 11281–11291. [CrossRef]

85. Bosse, A.; Fer, I. Mean Structure and Seasonality of the Norwegian Atlantic Front Current Along the Mohn Ridge From Repeated Glider Transects. *Geophys. Res. Lett.* **2019**, *46*, 13170–13179. [CrossRef]

86. Coppola, L.; Raimbault, P.; Mortier, L.; Testor, P. Monitoring the environment in the northwestern Mediterranean Sea. *Eos Trans. Am. Geophys. Union* **2019**, *100*. [CrossRef]

87. Testor, P.; Durrieu de Madron, X.; Mortier, L.; d'Ortenzio, F.; Legoff, H.; Dausse, D.; Labaste, M.; Houpert, L. *LION Observatory Data*; Technical Report; SEANOE (SEA scieNtific Open Data Edition), IFREMER: Issy-les-Moulineaux, France, 2020; [CrossRef]

88. Coppola, L.; Diamond Riquier, E.; Carval, T. *Dyfamed Observatory Data*; Technical Report; SEANOE (SEA scieNtific Open Data Edition), IFREMER: Issy-les-Moulineaux, France, 2020; [CrossRef]
89. Garau, B.; Ruiz, S.; Zhang, W.G.; Pascual, A.; Heslop, E.; Kerfoot, J.; Tintoré, J. Thermal lag correction on Slocum CTD glider data. *J. Atmos. Ocean. Technol.* **2011**, *28*, 1065–1071. [CrossRef]

90. Rudnick, D.L.; Sherman, J.T.; Wu, A.P. Depth-Average Velocity from Spray Underwater Gliders. *J. Atmos. Ocean. Technol.* **2018**, *35*, 1665–1673. [CrossRef]

91. Xing, X.; Claustre, H.; Blain, S.; d'Ortenzio, F.; Antoine, D.; Ras, J.; Guinet, C. For chlorophyll fluorescence acquired by autonomous platforms: A case study with instrumented elephant seals in the Kerguelen region (Southern Ocean. *Limnol. Oceanogr.* **2012**, *10*, 483–495. [CrossRef]

92. Lavigne, H.; D'Ortenzio, F.; Claustre, H.; Poteau, A. Towards a merged satellite and in situ fluorescence ocean chlorophyll product. *Biogeosciences* **2012**, *9*, 2111–2125. [CrossRef]

93. Marsaleix, P.; Auclair, F.; Floor, J.W.; Herrmann, M.J.; Estournel, C.; Pairaud, I.; Ulses, C. Energy conservation issues in sigma-coordinate free-surface ocean models. *Ocean. Model.* **2008**, *20*, 61–89. [CrossRef]

94. Lévy, M.; Resplandy, L.; Klein, P.; Capet, X.; Iovino, D.; Éthé, C. Grid degradation of submesoscale resolving ocean models: Benefits for offline passive tracer transport. *Ocean. Model.* **2012**, *48*, 1–9. [CrossRef]

95. Kessouri, F.; Ulses, C.; Estournel, C.; Marsaleix, P.; Séverin, T.; Pujo-Pay, M.; Caparros, J.; Raimbault, P.; de Fommervault, O.P.; d'Ortenzio, F.; et al. Nitrogen and phosphorus budgets in the northwestern Mediterranean deep convection region. *J. Geophys. Res. Ocean.* **2017**, *122*, 9429–9454. [CrossRef]

96. Kessouri, F.; Ulses, C.; Estournel, C.; Marsaleix, P.; d'Ortenzio, F.; Severin, T.; Taillandier, V.; Conan, P. Vertical mixing effects on phytoplankton dynamics and organic carbon export in the western Mediterranean Sea. *J. Geophys. Res. Ocean.* **2018**, *123*, 1647–1669. [CrossRef]

97. Margirier, F.; Bosse, A.; Testor, P.; L'Hévéder, B.; Mortier, L.; Smeed, D. Characterization of Convective Plumes Associated With Oceanic Deep Convection in the Northwestern Mediterranean From High-Resolution In Situ Data Collected by Gliders. *J. Geophys. Res. Ocean.* **2017**, *122*, 9814–9826. [CrossRef]

98. Estournel, C.; Testor, P.; Damien, P.; D'Ortenzio, F.; Marsaleix, P.; Conan, P.; Kessouri, F.; Durrieu de Madron, X.; Coppola, L.; Lellouche, J.M.; et al. High resolution modeling of dense water formation in the north-western Mediterranean during winter 2012–2013: Processes and budget. *J. Geophys. Res. Ocean.* **2016**, *121*, 5367–5392. [CrossRef]

99. Gaspar, P.; Grégoris, Y.; Lefevre, J.M. A simple eddy kinetic energy model for simulations of the oceanic vertical mixing: Tests at station Papa and Long-Term Upper Ocean Study site. *J. Geophys. Res. Ocean.* **1990**, *95*, 16179–16193. [CrossRef]

100. Herrmann, M.J.; Estournel, C.; Déqué, M.; Marsaleix, P.; Sevault, F.; Somot, S. Dense water formation in the Gulf of Lions shelf: Impact of atmospheric interannual variability and climate change. *Cont. Shelf Res.* **2008**, *28*, 2092–2112. [CrossRef]

101. Lellouche, J.M.; Le Galloudec, O.; Drévillon, M.; Régnier, C.; Greiner, E.; Garric, G.; Ferry, N.; Desportes, C.; Testut, C.E.; Bricaud, C.; et al. Evaluation of global monitoring and forecasting systems at Mercator Océan. *Ocean. Sci.* **2013**, *9*, 57–81. [CrossRef]

102. Thomas, L.N.; Lee, C.M. Intensification of ocean fronts. *J. Phys. Oceanogr.* **2005**, *86*, 783. [CrossRef]

103. Stern, M.E. Interaction of a uniform wind stress with a geostrophic vortex. *Deep. Sea Res.* **1965**, *12*, 355–367. [CrossRef]

104. Pedlosky, J. *Geophysical Fluid Dynamics*; Springer: Berlin/Heidelberg, Germany, 1979.

105. Thomsen, S.; Kanzow, T.; Krahmann, G.; Greatbatch, R.J.; Dengler, M.; Lavik, G. The formation of a subsurface anticyclonic eddy in the Peru-Chile Undercurrent and its impact on the near-coastal salinity, oxygen, and nutrient distributions. *J. Geophys. Res. Ocean.* **2016**, *121*, 476–501. [CrossRef]

106. Rudnick, D.L.; Cole, S.T. On sampling the ocean using underwater gliders. *J. Geophys. Res.* **2011**, *116*, 1–12. [CrossRef]

107. Bosse, A.; Testor, P.; Mortier, L.; Prieur, L.; Taillandier, V.; D'Ortenzio, F.; Coppola, L. Spreading of Levantine Intermediate Waters by submesoscale coherent vortices in the northwestern Mediterranean Sea as observed with gliders. *J. Geophys. Res. Ocean.* **2015**, *120*, 1599–1622. [CrossRef]

108. Giordani, H.; Bourdallé-Badie, R.; Madec, G. An Eddy-Diffusivity Mass-Flux Parameterization for Modeling Oceanic Convection. *J. Adv. Model. Earth Syst.* **2020**, *12*, e2020MS002078. [CrossRef]

109. Brannigan, L.; Marshall, D.P.; Naveira-Garabato, A.; Nurser, A.G. The seasonal cycle of submesoscale flows. *Ocean. Model.* **2015**, *92*, 69–84. [CrossRef]

110. Béthoux, J.; Prieur, L. Hydrologie et circulation en Méditerranée nord-occidentale. *Pétroles Tech.* **1983**, *299*, 25–34.

111. Niewiadomska, K.; Claustre, H.; Prieur, L.; D'Ortenzio, F. Submesoscale Physical-Biogeochemical Coupling across the Ligurian Current (Northwestern Mediterranean) Using a Bio-Optical Glider. *Limnol. Oceanogr.* **2008**, *53*, 2210–2225. [CrossRef]

112. Morel, A.; Berthon, J.F. Surface pigments, algal biomass profiles, and potential production of the euphotic layer: Relationships reinvestigated in view of remote-sensing applications. *Limnol. Oceanogr.* **1989**, *34*, 1545–1562. [CrossRef]

113. Boccaletti, G.; Ferrari, R.; Fox-Kemper, B. Mixed Layer Instabilities and Restratification. *J. Phys. Oceanogr.* **2007**, *37*, 2228–2250. [CrossRef]

114. Bachman, S.; Taylor, J. Modelling of partially-resolved oceanic symmetric instability. *Ocean. Model.* **2014**, *82*, 15–27. [CrossRef]

115. Fer, I.; Peterson, A.K.; Ullgren, J.E. Microstructure measurements from an underwater glider in the turbulent Faroe Bank Channel Overflow. *J. Atmos. Ocean. Technol.* **2014**, *31*, 1128–1150. [CrossRef]

116. Testor, P.; Mortier, L.; Coppola, L.; Claustre, H.; D'Ortenzio, F.; Bourrin, F.; Durrieu De Madron, X.; Raimbault, P. *Glider MOOSE Sections*; Technical Report; SEANOE (SEA scieNtific Open Data Edition), IFREMER: Issy-les-Moulineaux, France, 2017; [CrossRef]