# Fine-scale structures as spots of increased fish concentration in the open ocean

Alberto Baudena, Enrico Ser-Giacomi, Donatella d'Onofrio, Xavier Capet, Cédric Cotté, Yves Cherel, Francesco d'Ovidio

► To cite this version:

Alberto Baudena, Enrico Ser-Giacomi, Donatella d'Onofrio, Xavier Capet, Cédric Cotté, et al. Fine-scale structures as spots of increased fish concentration in the open ocean. Scientific Reports, 2021, 11, pp.15805. 10.1038/s41598-021-94368-1. hal-02431907

HAL Id: hal-02431907

https://hal.science/hal-02431907

Submitted on 7 Sep 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
<figure><img src="image_2.png" /></figure>

Check for updates

OPEN

# Fine-scale structures as spots of increased fish concentration in the open ocean

Alberto Baudena<fnref n="1" />,²✉, Enrico Ser-Giacomi<fnref n="1" />,⁶, Donatella D'Onofrio³,⁴, Xavier Capet<fnref n="1" />, Cedric Cotté<fnref n="1" />, Yves Cherel⁵ & Francesco D'Ovidio<fnref n="1" />

Oceanic frontal zones have been shown to deeply influence the distribution of primary producers and, at the other extreme of the trophic web, top predators. However, the relationship between these structures and intermediate trophic levels is much more obscure. In this paper we address this knowledge gap by comparing acoustic measurements of mesopelagic fish concentrations to satellite-derived fine-scale Lagrangian Coherent Structures in the Indian sector of the Southern Ocean. First, we demonstrate that higher fish concentrations occur more frequently in correspondence with strong Lagrangian Coherent Structures. Secondly, we illustrate that, while increased fish densities are more likely to be observed over these structures, the presence of a fine-scale feature does not imply a concomitant fish accumulation, as other factors affect fish distribution. Thirdly, we show that, when only chlorophyll-rich waters are considered, front intensity modulates significantly more the local fish concentration. Finally, we discuss a model representing fish movement along Lagrangian features, specifically built for mid-trophic levels. Its results, obtained with realistic parameters, are qualitatively consistent with the observations and the spatio-temporal scales analysed. Overall, these findings may help to integrate intermediate trophic levels in trophic models, which can ultimately support management and conservation policies.

Marine biomass distribution is highly patchy and variable in time across the entire trophic web<fnref n="1" />,². Discerning the factors underpinning ocean patchiness is fundamental to understand how they influence biogeochemical reactions and ecosystem stability³,⁴. These issues are pivotal for conservation purposes⁵ and for assessing the impact of climate change on the marine environment⁶.

One of the origins of the heterogeneity of biotic fields is the dynamic nature of ocean environments, which transports water masses and affects their properties on a large range of temporal scales, including those of ecological relevance. In this regard, the mesoscale and submesoscale processes⁷,⁸, now commonly referred to together as “fine-scales” (from a few to hundreds of kilometers) influence importantly the ecological landscape. One fruitful approach for capturing the structuring effect of fine-scale dynamics is the extraction of so-called Lagrangian Coherent Structures, (LCS⁹,¹⁰). LCSs provide several types of information regarding flow properties, such as the location of barriers to transport¹¹, or retentive and coherent regions¹². In particular, LCSs permit to identify frontal features (hereafter, Lagrangian fronts), usually associated with the presence of strong environmental gradients (see¹³ for a strict definition). One of the most common Lagrangian diagnostics used to determine LCSs is the Finite-size Lyapunov Exponent (FSLE¹⁴). This quantity measures the exponential rate of water parcel deformation and has maximal values (ridges) over frontal regions.

By shaping and elongating water patches, Lagrangian Coherent Structures have been demonstrated to set the frontiers of phytoplanktonic patches in terms of chlorophyll concentration¹⁰, and even functional type¹⁵. This in turn enhances contacts between different communities, regulating plankton diversity¹⁶,¹⁷.

More recently, advances in biologging programs provided evidence on the impact of fine-scale structures on top predators behavior. The concentration of predators foraging efforts has been observed in association with

<footnote n="1">¹Sorbonne Université, CNRS, IRD, MNHN, Laboratoire d'Océanographie et du Climat: Expérimentations et Approches Numériques (LOCEAN-IPSL), Paris, France. ²Sorbonne Université, CNRS, Laboratoire d'Océanographie de Villefranche, UMR 7093 LOV, Villefranche-sur-Mer, France. ³Institute of Atmospheric Sciences and Climate, National Research Council (CNR-ISAC), Torino, Italy. ⁴Copernicus Institute of Sustainable Development, Utrecht University, Utrecht, The Netherlands. ⁵Centre d'Etudes Biologiques de Chizé (CEBC), UMR 7372 du CNRS-La Rochelle Université, 79360 Villiers-en-Bois, France. ⁶Department of Earth, Atmospheric and Planetary Sciences, Massachusetts Institute of Technology, 02139 Cambridge, MA, USA. ✉email: alberto.baudena@gmail.com</footnote>
<table><caption>**Table 1.** Details of the acoustic transects analyzed.</caption><thead><tr><th colspan="5">Acoustic campaigns details</th></tr><tr><th>Cruise</th><th>Season</th><th>St. Date</th><th>End Date</th><th>Distance (km)</th></tr></thead><tbody><tr><td>LOGIPEV193_RUNKER</td><td>Summer</td><td>09/02/2013</td><td>17/02/2013</td><td>2752</td></tr><tr><td>LOGIPEV193_KERMAU</td><td>Summer</td><td>04/03/2013</td><td>10/03/2013</td><td>3781</td></tr><tr><td>OP2013-2_RUNKER</td><td>Winter</td><td>30/08/2013</td><td>10/09/2013</td><td>3310</td></tr><tr><td>LOGIPEV197_RUNKER</td><td>Summer</td><td>06/01/2014</td><td>13/01/2014</td><td>2800</td></tr><tr><td>LOGIPEV197_KERMAU</td><td>Summer</td><td>06/02/2014</td><td>18/02/2014</td><td>2045</td></tr><tr><td>OP2014-2_RUNKER</td><td>Winter</td><td>24/08/2014</td><td>04/09/2014</td><td>3677</td></tr></tbody></table>

Lagrangian fronts<sup>18,19</sup>. Furthermore, fronts detected by Lagrangian Coherent Structures have been observed to influence predators movements<sup>20</sup>. This could enhance energy transfer and gain<sup>21</sup>.

However, while the influence of Lagrangian fronts has been observed on both extremes of the trophic web, much less is known about mid-trophic levels. With this term we will refer to micronektonic organisms such as mesopelagic fish, which are able to actively move along the horizontal. Prants et al. (see in particular<sup>22</sup>) demonstrated a correlation between Pacific saury catches and Lyapunov exponents, and<sup>23</sup> found that several fishery vessels track LCSs when targeting fishery spots. However, these results leave some concerns about possible biases because commercial fisheries provides only isolated observations. In addition, fishing vessels may intentionally target frontal systems through the satellite images they are known to use. Unbiased fish measurements have been instead recently used by<sup>24</sup> to analyse the relationship between a frontal system and acoustic measurements in a coastal upwelling system. This allowed the authors to highlight the different role played by in-shore and off-shore waters. In terms of the mechanisms which can explain how fine-scale structures influence mid-trophic biomass distribution, even less is known. Classical explanations are based primarily on bottom-up mechanisms along fronts with intense upwelling<sup>8,25,26</sup>. However, these hypotheses do not take into account the necessity of a maturation time, which in the case of fish is consistently longer than both the growth response of lower trophic levels (such as phyto- and zooplankton) and the lifetime of the front<sup>27,28</sup>. Nor is the fish behavior considered, despite the fact that fish possess efficient sensorial and horizontal swimming capacities (<sup>29</sup> and Supporting Information SI.2).

Note that fine-scale dynamics include mixed-layer frontal activity generated by kilometer-scale instability processes<sup>30</sup>. These instabilities are a priori unlikely to influence mesopelagic fish because by definition they are restricted to the mixed layer<sup>8</sup>. Therefore, we consider here only the deep-reaching fronts associated with mesoscale stirring captured by altimetry-derived Lagrangian methods.

This study was conducted in the subantarctic area of the Southern Indian Ocean. The functioning of this region is mainly regulated by the Kerguelen plateau, a major topographic barrier for the Antarctic Circumpolar Current (ACC). The plateau fertilizes in iron, a limiting nutrient, the so-called "high-nutrients-low-chlorophyll" waters advected by the ACC<sup>31</sup>. Depending on seasonal light conditions and stratification of the water column, this provokes a large annual phytoplanktonic bloom, which supports a rich trophic web. This is one of the reasons for which the Kerguelen archipelago and its surrounding waters are part of one of the ten largest marine protected areas in the world (http://www.mpatlas.org/).

In this region, myctophids, also known as lantern fish, are one of the most abundant groups of mesopelagic fish. They are also present in other oceans worldwide and are thought to constitute one the largest portions of world fish biomass<sup>32</sup>. They also represent important prey items for numerous predators<sup>33</sup>. Myctophids are reported to play a central role in the carbon export to deep sea depths, and are suspected to affect the climate<sup>34</sup>. Constituting a potentially massive harvestable resource, they are threatened to be exploited in the near future<sup>35</sup>. In addition, their description in trophic models should be improved, and a better understanding of the mechanisms regulating their biomass is pivotal for fisheries management and sound marine spatial planning<sup>36</sup>.

The objective of the present study is to analyze the relationship between fine-scale structures in the open ocean and mid-trophic levels with unbiased, direct observations of fish concentrations. More specifically, the aim is to explore the degree to which fish distribution is shaped by fine-scale features. In order to achieve this, we firstly analyze the difference in fish concentrations observed in presence or absence of fronts. Secondly, we investigate the relationship between the intensity of the frontal dynamics and the fish distribution. Third, we examine this relationship on chlorophyll-rich waters only. Finally, we propose a mechanism by which fine-scale processes could potentially organize mid-trophic level biomass. This model considers explicitly active fish movement, and is calibrated with realistic, non-fitted parameters, derived from observations. The objective of the model is to assess (i) whether the modeled fish concentration varies significantly in presence of frontal features and (ii) if the spatio-temporal scales of the modeled process are comparable with those observed in the present study.

## Materials and methods

**Acoustic measurements.** A set of acoustic echo sounder data was used to analyze fish density. This was collected within the Mycto-3D-MAP program using split-beam echo sounders at 38 and 120 kHz. The Mycto-3D-MAP program included multiple large-scale oceanographic surveys during 2 years and a dedicated cruise in the Kerguelen area. The dataset was collected during 4 large-scale surveys in 2013 and 2014, both in summer (including both northward and southward transects) and in winter, corresponding to 6 acoustic transects of 2860 linear kilometers (see Table 1 for more details). Note that all legs except summer 2014 (MYCTO-3D-Map
cruise) were logistic operations, during which environmental in situ data (such as temperature or salinity profiles) could not be collected. The data were then treated with a bi-frequency algorithm, applied to the 38 and 120 kHz frequencies (details of data collection and processing are provided in³⁷). This provides a quantitative estimation of the concentration of gas-bearing organisms, mostly attributed to fish containing a gas-filled swimbladder in the water column³⁸. Most mesopelagic fish present swimbladders and several works pointed out that myctophids are the dominant mesopelagic fish in the region³⁹. Therefore, we considered the acoustic signal as mainly representative of myctophids concentration. Data were organized in acoustic units, averaging acoustic data over 1.1 km along the ship trajectory on average. Myctophid school length is in the order of tens of meters⁴⁰. For this reason, acoustic units were considered as not autocorrelated. Every acoustic unit is composed of 30 layers, ranging from 10 to 300 meters (30 layers in total).

The dataset was used to infer the Acoustic Fish Concentration (AFC) in the water column. We considered as AFC of the point $(x_i, y_i)$ of the ship trajectory the average of the bifrequency acoustic backscattering of 29 out of 30 layers (the first layer, 0-10 m, was excluded due to surface noise) AFC quantity is dimensionless.

As the previous measurements were performed through acoustic measurements, a non-invasive technique, fishes were not handled for this study.

**Regional data selection.** The geographic area of interest of the present study is the Southern Ocean. To select the ship transects belonging to this region, we used the ecopartition of⁴¹. Only points falling in the *Antarctic Southern Ocean* region were considered. We highlight that this choice is consistent with the ecopartition of⁴² (group 5), which is specifically designed for myctophids, the reference fish family (Myctophidae) of this study. Furthermore, this choice allowed us to exclude large scale fronts (i.e., fronts that are visible on monthly or yearly averaged maps) which have been the subject of past research works⁴³,⁴⁴. In this way our analysis is focused specifically on fine-scale fronts.

**Day-night data separation.** Several species of myctophids present a diel vertical migration. They live at great depths during the day (between 500 and 1000 m), ascending at dusk in the upper euphotic layer, where they spend the night. Since the maximal depth reached by the echo sounder we used is 300 m, in the analysis reported in Figs. 2 and 3 we excluded data collected during the day. However, their analysis is reported in SI.1. A restriction of our acoustic analyses to the ocean upper layer is also consistent with the fact that the fine-scale features we computed are derived in this work by satellite altimetry, thus representative of the upper part of the water column (~ 50 m). Finally, we note that the choice of considering the echo sounder data in the first 300 m of the water columns is coherent with the fact that LCS may extend almost vertically in depth even at 600 m depth⁴⁵,⁴⁶ and with the fact that altimetry-derived velocity fields are consistent with subsurface currents in our region of interest down to 500 m²⁰.

**Satellite data.** Velocity current data and Finite-Size Lyapunov Exponent (FSLE) processing. Velocity currents are obtained from Sea Surface Height (SSH), which is measured by satellite altimetry, through geostrophic approximation. Data, which were downloaded from E.U. Copernicus Marine Environment Monitoring Service (CMEMS, http://marine.copernicus.eu/), were obtained from DUACS (Data Unification and Altimeter Combin nation System) delayed-time multi-mission altimeter, and displaced over a regular grid with spatial resolution of 1/4 × 1/4° and a temporal resolution of 1 day.

Trajectories were computed with a Runge-Kutta scheme of the 4th order with an integration time of 6 hours. Finite-size Lyapunov Exponents (FSLE) were computed following the methods in¹⁴, with initial and final separa tion of 0.04° and 0.4° respectively. This Lagrangian diagnostic is commonly used to identify Lagrangian Coherent Structures. This method determines the location of barriers to transport, and it is usually associated with oceanic fronts⁹. Details on the Lagrangian techniques applied to altimetry data, including a discussion of its limitation, can be found in¹⁰.

*Temperature field and gradient computation.* The Sea Surface Temperature (SST) field was produced from the OSTIA global foundation Sea Surface Temperature (product id: SST_GLO_SST_L4_NRT_ OBSERVATIONS_010_001) from both infrared and microwave radiometers, and downloaded from CMEMS website. The data are represented over a regular grid with spatial resolution of 0.05 × 0.05° and daily-mean maps. The SST gradient was obtained from:

$$ GradSST = \sqrt{g_x^2 + g_y^2} $$

where $g_x$ and $g_y$ are the gradients along the west-east and the north-south direction, respectively. To compute $g_x$, the following expression was used:

$$ g_x = \frac{1}{2dx} \cdot (SST_{i+1} - SST_{i-1}) $$

where the SST values of the adjacent grid cells (along the west-east direction: cells $i + 1$ and $i - 1$) were employed. $dx$ identifies the kilometric distance between two grid points along the longitude (which varies with latitude). The definition is analog for $g_y$, considering this time the north-south direction and $dy \simeq 5$ km (0.05°).

*Chlorophyll field.* Chlorophyll estimations were obtained from the Global Ocean Color product (OCEAN-COLOUR_GLO_CHL_L4_REP_ OBSERVATIONS_009_082-TDS) produced by ACRI-ST. This was down-
loaded from CMEMS website. Daily observations were used, in order to match the temporal resolution of the acoustic and satellite observations. The spatial resolution of the product is 1/24°.

**Estimation of satellite data along ship trajectory.** For each point $(x_i, y_i)$ of the ship trajectory, we extracted a local value of FSLE, SST gradient, and chlorophyll concentration. These were obtained by considering the respective average value in a circular around of radius $\sigma$ of each point $(x_i, y_i)$. Different $\sigma$ were tested (ranging from 0.1° to 0.6°), and the best results were obtained with $\sigma = 0.2°$, reference value reported in the present work. This value is consistent with the resolution of the altimetry data.

**Statistical processing. Front identification.** FSLE and SST gradient were used as diagnostics to detect frontal features, since they are typically associated with front intensity and Lagrangian Coherent Structures¹⁰. Note that the two diagnostics provide similar but not identical information. FSLEs are used to analyze the horizontal transport and to identify material lines along which a confluence of waters with different origins occur. These lines typically display an enhanced SST gradient because water masses of different origin have often contrasted SST signature. However, this is not a necessary condition. FSLE ridges may be invisible in SST maps if transport occurs in a region of homogeneous SST. Conversely, SST gradient unveils structures separating waters of different temperatures, whose contrast is often – but not always – associated with horizontal transport. Therefore, even if they usually detect the same structures, these two metrics are complementary. Frontal features were identified by considering a local FSLE (or SST gradient, respectively) value larger than a given threshold. The threshold values have been chosen heuristically but consistently with the ones found in previous works. For the FSLEs, we used 0.08 days⁻¹, a threshold value in the range of the ones chosen in¹⁸ and⁴⁷. For the SST gradient, we considered representative of thermal front values greater than 0.009°C/km, which is about half the value chosen in⁴⁷. However, in that work, the SST gradient was obtained from the advection of the SST field with satellite-derived currents for the previous 3 days, a procedure which is known to enhance tracer gradients.

**Bootstrap method.** The threshold value splits the AFC into two groups: AFC co-located with FSLE or SST gradient values over the threshold are considered as measured in proximity of a front (i.e., statistically associated with a front), while AFC values below the threshold are considered as not associated with a frontal structure. The statistical independence of the two groups was tested using a Mann-Whitney or U test. Finally, bootstrap analysis is applied following the methodologies used in⁴⁷. This allows estimating the probability that the difference in the mean AFC values, over and under the threshold, is significant, and not the result of statistical fluctuations. Other diagnostics tested are reported in SI.1.

**Linear quantile regression.** Linear quantile regression method⁴⁸ was employed as no significant correlation was found between the explanatory and the response variables. This can be due to the fact that multiple factors (such as prey or predator distributions) influence the fish distribution other than the frontal activity considered. The presence of these other factors can shadow the relationship of the explanatory variables (in this case, the FSLE and the SST gradient) with the mean value of the response variable (the AFC). A common method to address this problem is the use of the quantile regression⁴⁸,⁴⁹, that explores the influence of the explanatory variables on other parts of the response variable distribution. Previous studies, adopting this methodology, revealed the limiting role played by the explanatory variables in the processes considered⁵⁰. The percentiles values used are 75th, 90th, 95th, and 99th. The analysis is performed using the statistical package QUANTREG in R (https://CRAN.R-project.org/package=quantreg, v.5.38⁴⁸,⁵¹), using the default settings.

**Chlorophyll-rich waters selection.** The AFC observations were considered in chlorophyll-rich waters if the corresponding chlorophyll concentration was higher than a given threshold. All the other AFC measurements were excluded, and a linear regression performed between the remaining AFC and FSLE (or SST gradient) values. The corresponding thresholds (one for FSLE and one for SST gradient case) were selected though a sensitivity test reported in SI.1. These resulted in 0.22 and 0.17 mg/m³ for FSLE and for SST gradient, respectively. These values are consistent among them and, in addition, they are in coherence with previous estimates of chlorophyll concentration used to characterise productive waters in the Southern Ocean (0.26mg/m³⁵²).

**Gradient climbing model.** An individual-based mechanistic model is built to describe how fish could move along frontal features. We assume that the direction of fish movement along a frontal gradient is influenced by the noise of the prey field (SI. 2). Specifically, we consider a Markovian process along the (one dimensional) cross-front direction. Time is discretized in timesteps of length $\Delta\tau$. We presuppose that, at each timestep, the fish chooses between swimming in one of the two opposite cross-front directions (“left” and “right”). This decision depends on the comparison between the quantity of a tracer (a cue) present at its actual position and the one perceived at a distance $p_R$ from it, where $p_R$ is the perceptual range of the fish. We consider the fish immersed in a tracer whose concentration is described by the function $T(x)$.

An expression for the average velocity of the fish, $U_F(x)$, can now be derived. We assume that the fish is able to observe simultaneously the tracer to its right and its left. Fish sensorial capacities are discussed in SI.2. The tracer observed is affected by noise. Noise distribution is considered uniform, with $-\xi_{MAX} < \xi < \xi_{MAX}$, $\xi_{MAX} > 0$.

The effective values perceived by the fish, at its left and its right, will be, respectively:
$$ 
\begin{aligned}
\tilde{T}(x_0 - \Delta x) &= T(x_0 - \Delta x) + \xi_1 \\
\tilde{T}(x_0 + \Delta x) &= T(x_0 + \Delta x) + \xi_2.
\end{aligned}
$$

We assume that, if $\tilde{T}(x_0 + \Delta x) > \tilde{T}(x_0 - \Delta x)$, the fish will move to the right, and, vice versa, to the left. We hypothesize that the observational range is small enough to consider the tracer variation as linear, as illustrated in Fig. S7 (SI. 3). In this way:

$$ 
\begin{aligned}
\tilde{T}(x_0 + \Delta x) &= T(x_0) + p_R \frac{\partial T}{\partial x} + \xi_1 \\
\tilde{T}(x_0 - \Delta x) &= T(x_0) - p_R \frac{\partial T}{\partial x} + \xi_2.
\end{aligned}
$$

In case of absence of noise, or with $\xi_{MAX} < p_R \frac{\partial T}{\partial x}$, the fish will always move in the correct direction, in that it will climb the gradient. Assuming V as the cruising swimming velocity of the fish, this means $U_F(x) = V$.

Let's now assume $\xi_{MAX} > p_R \frac{\partial T}{\partial x}$. If $T(x_0 + \Delta x) > T(x_0 - \Delta x)$ (as in Fig. S7), and the fish will swim leftward if

$$ \xi_1 - \xi_2 > 2p_R \frac{\partial T}{\partial x}. $$

Since $\xi_1$ and $\xi_2$ range both between $-\xi_{MAX}$ and $\xi_{MAX}$, we can obtain the probability of leftward moving P(L). This will be the probability that the difference between $\xi_1$ and $\xi_2$ is greater than $2p_R \frac{\partial T}{\partial x}$.

$$ 
\begin{aligned}
P(L) &= \frac{1}{8\xi_{MAX}^2} \left( 2\xi_{MAX} - 2p_R \frac{\partial T}{\partial x} \right)^2 \\
&= \frac{1}{2} \left( 1 - \frac{p_R}{\xi_{MAX}} \frac{\partial T}{\partial x} \right)^2
\end{aligned}
$$

The probability of moving right will be

$$ P(R) = 1 - P(L) $$

and their difference gives the frequency of rightward moving

$$ 
\begin{aligned}
P(R) - P(L) &= 1 - 2P(L) = 1 - \left( 1 - \frac{p_R}{\xi_{MAX}} \frac{\partial T}{\partial x} \right)^2 \\
&= \frac{p_R}{\xi_{MAX}} \frac{\partial T}{\partial x} \left( 2 - \frac{p_R}{\xi_{MAX}} \left| \frac{\partial T}{\partial x} \right| \right),
\end{aligned}
$$

where the absolute value of $\frac{\partial T}{\partial x}$ has been added to preserve the correct climbing direction in case of negative gradient. The above expression leads to:

$$ U_F(x) = \frac{V p_R}{\xi_{MAX}} \frac{\partial T}{\partial x} \left( 2 - \frac{p_R}{\xi_{MAX}} \left| \frac{\partial T}{\partial x} \right| \right). \quad (1) $$

We then assume that, over a certain value of tracer gradient $\frac{\partial T}{\partial x}$, the fish are able to climb the gradient without being affected by the noise. This assumption, from a biological perspective, means that the fish does not live in a completely noisy environment, but that under favorable circumstances it is able to correctly identify the swimming direction that leads to higher prey availability. This means that

$$ p_R * \frac{\partial T}{\partial x} = \xi_{MAX}. \quad (2) $$

Substituting then (2) into (1) gives:

$$ U_F(x) = V \frac{\frac{\partial T}{\partial x}}{\frac{\partial T}{\partial x} \max} \left( 2 - \frac{\left| \frac{\partial T}{\partial x} \right|}{\frac{\partial T}{\partial x} \max} \right). \quad (3) $$

Finally, we can include an eventual effect of transport by the ocean currents, considering that the tracer is advected passively by them, simply adding the current speed $U_C$ to Expr. (3).

The representations provided are individual based, with each individual representing a single fish. However, we note that all the considerations done are also valid if we consider an individual representing a fish school. $U_F$ will then simply represent the average velocity of the fish schools. This aspect should be stressed since many fish species live and feed in groups, especially myctophids (see SI.2 for further details).

**Continuity equation in one dimension.** The solution of this model will now be discussed. The continuity equation is first considered in one dimension. The starting scenario is simply an initially homogeneous distribution of fish, that are moving in a one dimensional space with a velocity given by $U_F(x)$.
We assume that in the time scales considered (few days to some weeks), the fish biomass is conserved, so for instance fishing mortality or growing rates are neglected. In that case, we can express the evolution of the concentration of the fish ρ by the continuity equation

$$ \frac{\partial \rho}{\partial t} + \nabla \cdot \mathbf{j} = 0 \qquad (4) $$

in which $\mathbf{j} = \rho U_F(x)$, so that Eq. (4) becomes

$$ \frac{\partial \rho}{\partial t} + \nabla \cdot (\rho U_F(x)) = 0. \qquad (5) $$

In one dimension, the divergence is simply the partial derivate along the x-axis. Eq. (5) becomes

$$ \frac{\partial \rho}{\partial t} = -\frac{\partial}{\partial x} \left( \rho U_F \right) \qquad (6) $$

Now, we decompose the fish concentration ρ in two parts, a constant one and a variable one $\rho = \rho_0 + \tilde{\rho}$. Eq. (6) will then become

$$ \frac{\partial \rho}{\partial t} = -U_F \frac{\partial \tilde{\rho}}{\partial x} - \rho \frac{\partial U_F}{\partial x}. \qquad (7) $$

Using Expr. (3), Eq. (7) is numerically simulated with the Lax method. In Expr. (3) we impose that $U_F(x) = V$ when $U_F(x) > V$.

This biological assumption means that fish maximal velocity is limited by a physiological constraint rather than by gradient steepness. Details of the physical and biological parameters are provided in SI.6.

## Results

**Relationships between acoustic fish concentration and satellite-derived diagnostics.** A good qualitative agreement is found when comparing the AFC along some portions of the ship transect with the local FSLE and SST gradient fields (Fig. 1). AFCs seem to increase in correspondence with the frontal features identified by high values of the Lyapunov exponents and temperature gradient. A visual exploration of the relation between AFC and ridges of in the FSLE field is difficult to interpret, due to the intermittent nature of the AFC and to the fact that the ship transects intersect FSLE ridges with various angles.

**Bootstrap analysis.** In order to quantify whether the AFC values present significant differences in proximity of fine-scale features, a bootstrap analysis was conducted. Therefore, AFCs were separated in two groups: those falling over a front, identified by FSLE (or SST gradient) values over a threshold, and those falling outside it (further details are provided in Materials and Methods). Bootstrap analysis indicates that significantly higher AFC values are detected in presence of fronts (p-value < 0.001) for both criteria of front detection (FSLE or SST gradient, Fig. 2). This result demonstrates that high fish concentrations occur preferably over fronts detected with the diagnostics employed.

**Quantile regression analysis.** Subsequently, we analyzed the relationship between the intensity of the fronts (provided by the FSLE and SST gradient) and the AFC. When interpolating the data with a linear fit, no significant correlation was found between AFC values and FSLE or SST gradient. Therefore, linear quantile regression method was applied using the 75th, 90th, 95th, and 99th percentiles (Fig. 3; refer to Materials and Methods for further details). All the quantile regression slopes are significantly different from zero (Table S1 and sensitivity test in SI.1), meaning that FSLE and SST gradient intensities play a role in limiting fish concentration. This means that large fish concentrations are associated with the nearby presence of a strong FSLE or SST gradient value. Conversely, not all of the strong fronts are always associated with large fish densities.

**Chlorophyll-rich waters analysis.** The previous results evidence that, together with front intensity, other unaccounted factors play a role in shaping fish distribution. In the following, we investigate the importance of one of these elements, notably the presence of prey. We use the local chlorophyll concentration as a proxy to determine whether a front is found in a region potentially rich in prey for fish or not. This choice is driven by the fact that areas with enhanced primary production are also known to host higher trophic organisms²⁷,²⁸ such as zooplankton (a prey for myctophids). Therefore, we consider only the points of the ship trajectory whose local value of chlorophyll is larger than a threshold. The threshold value is identified through a sensitivity test in SI.1, and it is found consistent for both FSLE and SST gradient. When this selection is applied, a significant linear relationship emerges between the AFC and the intensity of the frontal dynamics (Fig. 4; FSLE: R = 0.66, p < 0.001; SST gradient: R = 0.54, p < 0.001).

**A fine-scale mechanism of fish aggregation.** Why do fish aggregate along frontal features? To try to address this question, we propose a simple mathematical model (see “Gradient climbing model” section for details). The model assumes that fish have a gradient climbing capacity, which is one of the most widespread movement mechanisms used in other biological contexts (e.g., chemotaxis⁵³). This gradient climbing capacity is specifically tuned for mid-trophic levels and myctophids and is based on a cue-pursuing dynamic. Fish try to climb a gradient of tracer. We considered this tracer a proxy of zooplankton concentration, main prey of several
<figure><img src="image_3.png" /><figcaption>**Figure 1.** (A) Overview of the study region. Blue dots show the position of the survey transects conducted during the Mycto-3D-MAP program. They are superimposed on the local bathymetry (gray scale at right). (B) Illustrative example of a transect of the ship trajectory on August, 29th, 2014, superimposed on a simultaneous field of Finite-size Lyapunov Exponents (i.e., frontal structures; upper gray scale). The color of each dot is proportional to the local Acoustic Fish Concentration (adimensional, scale at right). (C) Same as for panel (B), with the difference that the transect is superimposed on a simultaneous field of Sea Surface Temperature gradient (upper gray scale). Note that panels (B) and (C) show a region at the border of the domain analysed.</figcaption></figure>

fish species, including myctophids⁵⁴. At the scales considered in this study (10s of kilometers) zooplankton swimming capacities are restricted to the vertical axis⁵⁵. They can thus be considered as passive tracers along the horizontal. Along this spatial dimension, zooplankton aggregation and growth is usually driven by a relatively fast response to nutrients presence, in the order of days to weeks²⁷,²⁸,⁵⁶. In particular, this is valid also for zooplankton species present in our study region⁵⁷,⁵⁸. Conversely, fish have growth rates consistently slower: in particular, pelagic fishes and myctophids are considered as "slow-growing fish", with lifetimes spanning a few years⁵⁹. Therefore, aggregation can not be explained by a quick increase of their biomass. Alongside this, fish have extremely developed sensorial capacities and, in contrast to zooplankton, they can actively swim, with both capacities involved in many functional activities, including feeding²⁹. These arguments support our approach of modelling zooplankton as a passive tracer and fish as active swimmers (we invite the reader to refer to SI.2 for further details).

To account for ocean patchiness and small scale turbulence, a noise term perturbed the ability of the fish to properly identify the elevated zooplankton regions. However, we assumed that fish are able to orientate without problems over a given threshold of the zooplankton gradient. This threshold was estimated from the zooplankton concentrations (SI.5).

The distribution of the tracer used as a proxy of zooplankton concentration follows a sigmoid function (red line in Figs. 5 and 6), which models a generic local gradient. The transition zones at the edges of the tracer plateau are ~5 km wide, as typically found for fine-scale features such as eddies or upwelling structures. Conversely, the plateau presents a local maximum at its center. We consider the modeled gradient as representative of a frontal feature.
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Bootstrap method results. AFC values were separated according to the respective FSLE and SST gradient values, using a threshold as a limit to define a front (0.08 days⁻¹ and 0.009°C/km respectively;¹⁸,⁴⁷) Light gray columns represent the mean AFC under the respective threshold, considered outside a front. Dark gray columns represent the mean AFC over the threshold, thus in proximity of a front. Left columns refer to FSLE analysis, right columns to SST gradient. Error bars indicate the standard deviation, while black stars indicate the significance of the bootstrap test (***: *p*<0.001).</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>**Figure 3.** Scatter plot of AFC against FSLE (left panel) and SST gradient (right panel) for all the points sampled (thus not filtered based on chlorophyll). The lines, from the bottom to the top, indicate the linear quantile regressions at 75th, 90th, 95th and 99th percentiles. The analysis is used to investigate the relationship between the front intensity and just the higher values of fish concentration. Both axes are in logarithmic scale (values equal to zero are therefore not depicted). Values of the quantile regression coefficients are reported in Table S1.</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>**Figure 4.** Left panel: scatter plot of AFC against FSLE (left panel) for the sub-set of points whose local chlorophyll concentration is higher than 0.22 mg/m³. Right panel: scatter plot of AFC against SST gradient (right panel) for the sub-set of points whose local chlorophyll concentration is higher than 0.17 mg/m³. For both panels, the red line is obtained from the respective linear interpolation of the selected points (note that axes are reported in logarithmic scale). Left panel: R = 0.66, R² = 0.43, *p* < 0.001. Right panel: R = 0.54, R² = 0.29, *p* < 0.001. Sensitivity of the results with respect to the chlorophyll threshold chosen is reported in SI.1.</figcaption></figure>
<figure><img src="image_7.png" /><figcaption>**Figure 5.** Time evolution of the fish concentration (blue line, adimensional) according to the continuity equation. The tracer (red line, adimensional) describes a plateau of 8 km in width. At its limits, its values range from 1 to 0 in about 5 km. Each panel represents a different snapshot at 0, 6 hours, 1 day, and 4 days.</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>**Figure 6.** Time evolution of the fish concentration as reported in Fig. 5. This time, the plateau width has been set to 70 km.</figcaption></figure>

Four different snapshots (at 0, 6 hours, 1 day and 4 days) of the fish modelled concentration are illustrated in Figs. 5 and 6. By construction, fish displacements are most conspicuous in regions of elevated tracer gradients, i.e., at the plateau edges. Indeed, after only 6 hours, two peaks of doubled concentration are present in correspondence with the margins of the tracer plateau. In the following days, the two peaks decrease their growth rate, while the concentration between them increases until they merge together. The fish concentration is thus homogeneous over the plateau, presenting values between 2.5 and 3.5 times higher than the initial concentration. A larger plateau (Fig. 6) produces two peaks of fish concentration in correspondence with its edges, while the merging between the peaks occurs over longer timescales (not reported). Changing the type of fish behavior leads to similar results (see SI.3).

At the same time during which this mechanism occurs, the tracer can evolve. As a sensitivity test, we numerically analyzed a scenario in which the tracer, subjected to a typical frontal dynamics, is stretched in a filament and eroded by diffusion (SI.4). Using realistic bio-physical parameters representative of the study area, we found that these tracer dynamics do not compromise the aggregation mechanism presented above, but may even facilitate the aggregation of fish along frontal zones. Indeed, the model developed predicts several quantitative information, such as an estimate of the fish aggregation over time. The latter highlighted a final concentration, on average, an order of magnitude stronger than the initial one. In addition, it was possible to estimate a "fish concentration life time". This quantity indicates the amount of time during which a group of fish is able to follow a patch of interest before it vanishes due to the frontal dynamics. The school life times we obtained ranged between 7 to 25 days, with an average value of around 2 weeks. Remarkably, this amount of time matches with that of fine-scale processes (one day to several weeks).

The observations from the acoustic echo sounder have a coverage in space and time that is too limited to be used to test the validity of the model outputs. However, we highlight that the model predictions, obtained from realistic parameters without optimisation or fitting, are coherent with the spatio temporal scale analysed. Finally, we note that the modeled aggregation can occur only within prescribed conditions, notably the presence of prey. This is consistent with the results of both bootstrap, quantile, and chlorophyll-rich waters regression analyses, which pointed to the fact that fronts are not the only factor structuring fish abundance.

## Discussion

In the present work, we compare in situ measurements of Acoustic Fish Concentration with satellite-derived frontal structures, which were associated with the presence of Lagrangian fronts (detected through FSLE ridges) and strong Sea Surface Temperature gradients. When comparing the AFC with the presence (or absence) of a front, our findings illustrate that high AFCs are mostly associated with important frontal structures (Fig. 2, p< 0.001 both for FSLE and SST gradient). Subsequently, we analyze the relationship between AFC and front intensity. In this regard, we applied a quantile regression to investigate the relationship of FSLE and SST gradient with the upper part of the AFC distribution. This approach is consistent with the assumption that AFC is affected by several other processes (such as predators or prey presence), other than the presence of fronts. All the percentiles used (75th, 90th, 95th, and 99th) resulted in a quantile slope significantly different from zero (Fig. 3 and Table S1). This allowed us to unveil the important role played by front intensity on AFC, which could not have
been detected through a standard linear correlation. In this regard, the presence of fronts can be seen as a limiting condition for high fish concentrations. This means that strong fish aggregation is preconditioned by the intensity of a frontal feature. Conversely, not all of the strong fronts detected indicate high AFCs, since other processes are involved in shaping fish abundance. Therefore, we examined the effect of one of these factors, notably the presence of prey. This was identified by analysing chlorophyll fields, which were considered as proxy of potential fish preys (such as zooplankton). Indeed, fronts sustaining a strong primary production can support also large zooplankton densities²⁷,²⁸. When considering only the fronts which are co-located with the presence of potential fish prey (i.e., large chlorophyll concentrations), our results show that such fronts modulate significantly the concomitant fish distribution (*p*<0.001 for both FSLE and SST gradient, R=0.66 and R=0.54 respectively, Fig. 4).

The AFC response to frontal features cannot be explained with traditional mechanisms usually prescribed to lower trophic levels, such as rapid growth associated to the presence of nutrients, because fish have slower growth rates. To address this question, we proposed a gradient-climbing model specifically calibrated for the study of mid-trophic levels. Importantly, this model considers explicitly active fish movement which, to our knowledge, has not been included directly in the study of fish aggregation to date (e.g.²⁶,⁶⁰). Once parameterized with values typical for the study region, the model can produce fish behavior characterized by a quick response to gradient structures. Simulated peaks of doubled fish concentration appear after just 6 hours in response to a front, and grow afterwards with lifetimes of several days. Interestingly, fish do not converge at the center of the patch of interest, where the maximal concentration of the patch is found. When the width of the plateau is in the order of hundred kilometers (Fig. 6), two peaks of fish concentration are produced at its edges, where the tracer gradient is maximum. Instead, when the width of the plateau is in the order of ten kilometers (Fig. 5), after a few days fish tend to be homogeneously distributed over it. On-going studies on myctophids' response to food concentration presuppose that myctophids, over a certain prey density, ingest always the same quantity of food (Holling type III functional response; A. Hulley, personal communication). Therefore, above this threshold, they are expected to be homogeneously distributed. This is in qualitative agreement with the fish distribution predicted by the model.

The proposed mechanism of aggregation needs two obvious initial conditions: the presence of fish, and the presence of a zooplankton patch. It is presumably that fine-scale fronts that lack one or both these conditions can- not act as aggregating spots. Furthermore, while we assumed that the aggregation occurs after a certain amount of time, the environment studied is dynamic. Thus, it is likely that the aggregation mechanism was observed during different stages. Therefore, our model suggests that not all of the strong fronts should present high fish concen- trations, but that, conversely strong fish aggregation is preconditioned by the intensity of a frontal feature. This is in accordance with the results obtained from the quantile regression and the chlorophyll-rich water analyses.

Tracer patchiness is known to be associated with frontal features¹⁰,²⁵,⁶¹. However, patchiness evolves, and under the double effect of stretching and diffusion, local gradients can be eroded. Thus, we tested the robustness of our model to this feature in SI.4. As in the former case, none of the parameters employed has been optimized nor fitted, but they all represent rigorous estimations of Southern Ocean physical and biological conditions, such as stretching and mixing rates or fish cruising speed. Results allowed us to estimate a typical lifetime for a fine- scale patch of around two weeks, much longer than the peak doubling time (~6 hours). This robustness of the lifetime of the fine-scale patch is consistent with the hypothesis of a fixed tracer assumed in the gradient climbing scenario. Furthermore, we demonstrated that the stretching and diffusion dynamics can potentially enhance fish aggregations. In addition, the model provided estimations of other aggregation dynamics, such as the dimen- sions of the aggregation patterns or its intensity, with timescales comparable with those of fine-scale processes.

The modeled patch of interest is considered as a proxy of zooplankton concentration. However, being param- eterized as a passive tracer, it can be considered, more simply, as a generic water property detected by fish. Indeed, physical characteristics, such as temperature, have been proven to be used by predators to find favorable condi- tions where concentrated food occurs⁶².

Dedicated missions are required to quantitatively validate the model results. This is currently not possible due to the sparsity, in time and space, of our acoustic dataset. First, continuous measurement are necessary to follow the potential aggregating mechanism along its different phases. Secondly, the displacement of the observed water mass, during the sampling period, must be taken into account.

Within the ACC, the information available on mid-trophic levels reported that the large circumpolar fronts are known to host (i) large densities of zooplankton and myctophids and (ii) that these organisms are patchily distributed⁶³. In this study we suggest the potential mechanisms driving the patchiness observed at fine-scale. Assessing the preconditions and the other dynamics necessary for front aggregation is a new, open challenge emerging from the present work, whose findings look promising.

Note that in our work we focus only on open ocean Lagrangian fronts, that is, fine-scale frontal features induced by the mesoscale open ocean activity. In particular, by considering only the *Antarctic Southern Ocean* region defined by⁴¹, we intentionally exclude coastal fronts as well as large scale fronts. These can represent important productivity spots, and their dynamics and ecological role may be different (e.g.,²⁴,⁴³,⁴⁴). Finally, the limitations of our analysis must be discussed. No vertical dynamics has been included in the model presented. However, these play an important role in the organization of marine biota⁸, and, typically, stronger gradients are present along the vertical scale⁶⁴. This assumption is due to the fact that a 3D Lagrangian analysis, while appealing⁶⁵, would add a level of complexity and a number of further parameters difficult to deal with, due to a lack of observational information on the vertical velocities of the ocean. Future satellite missions (such as SWOT: https://swot.cnes.fr/en) will possibly help to mitigate this problem, by helping the assimilation scheme to better reconstruct the three dimensional dynamics⁶⁶. At the same time, they will also improve satellite resolution, providing a more precise location of fine-scale fronts.

The model studied does not consider the diel vertical migration of myctophids either. This choice is driven by the difficulty in parameterizing such a non linear behavior and in assimilating different migratory diel patterns. However, we note that many zooplankton species exhibit a diel cycle as well. Zooplankton consumption due to
fish foraging is considered to be negligible, or compensated by source terms (such as rapid growth subsequent to a phytoplanktonic bloom). Finally, the density of fish predators in space and time, as well as their feeding frequency should be considered in future studies. The limitations presented can open the way for future investigations.

The results presented here may be useful for improving the representation of intermediate trophic levels in coupled ecological and physical models⁶⁰, habitat models⁶⁷, targeting the mesopelagic compartment in particular. At the same time, the possibility of using Lagrangian Coherent Structures combined with ocean colour information as a proxy of higher fish concentrations may further improve the integration of satellite-derived Lagrangian tools in conservation planning⁶⁸.

## Data availability

All data used in the present study are available upon request.

Received: 22 December 2020; Accepted: 24 June 2021

Published online: 04 August 2021

## References

1. Bertrand, A. et al. Broad impacts of fine-scale dynamics on seascape structure from zooplankton to seabirds. *Nat. Commun.* **5**, 5239 (2014).

2. McManus, M. A. & Woodson, C. B. Plankton distribution and ocean dispersal. *J. Exp. Biol.* **215**(6), 1008 (2012).

3. Martin, A. Phytoplankton patchiness: the role of lateral stirring and mixing. *Prog. Oceanogr.* **57**(2), 125 (2003).

4. Lévy, M. & Martin, A. P. The influence of mesoscale and submesoscale heterogeneity on ocean biogeochemical reactions. *Global Biogeochem. Cycles* **27**(4), 1139 (2013).

5. Gaines, S. D., White, C., Carr, M. H. & Palumbi, S. R. Designing marine reserve networks for both conservation and fisheries management. *Proceedings of the National Academy of Sciences* (2010).

6. Hoegh-Guldberg, O. & Bruno, J. F. The impact of climate change on the world's marine ecosystems. *Science* **328**(5985), 1523 (2010).

7. McGillicuddy, D. J. Mechanisms of physical-biological-biogeochemical interaction at the oceanic mesoscale. *Ann. Rev. Marine Sci.* **8**(1), 125 (2016) (PMID: 26359818).

8. Lévy, M., Franks, P. J. S. & Smith, K. S. The role of submesoscale currents in structuring marine ecosystems. *Nat. Commun.* **9**(1), 1 (2018).

9. Haller, G. Lagrangian coherent structures. *Ann. Rev. Fluid Mech.* **47**(1), 137 (2015).

10. Lehahn, Y., d'Ovidio, F. & Koren, I. A satellite-based lagrangian view on phytoplankton dynamics. *Ann. Rev. Marine Sci.* **10**(1), 99 (2018).

11. Boffetta, G., Lacorata, G., Redaelli, G. & Vulpiani, A. Detecting barriers to transport: a review of different techniques. *Physica D Nonlinear Phenomena* **159**(1), 58 (2001).

12. d'Ovidio, F., De Monte, S., Della Penna, A., Cotté, C. & Guinet, C. Ecological implications of eddy retention in the open ocean: a Lagrangian approach. *J. Phys. A Math. Theor.* **46**(25), 254023 (2013).

13. Prants, S. V. Dynamical systems theory methods to study mixing and transport in the ocean. *Physica Scripta* **87**(3), 038115 (2013).

14. d'Ovidio, F., Fernández, V., Hernández-García, E. & López, C. Mixing structures in the Mediterranean Sea from finite-size Lyapunov exponents. *Geophys. Res. Lett.* **31**(17), (2004).

15. d'Ovidio, F., Monte, S. D., Alvain, S., Dandonneau, Y. & Lévy, M. Fluid dynamical niches of phytoplankton types. *Proc. Natl. Acad. Sci.* **107**(43), 18366 (2010).

16. De Monte, S., Soccodato, A., Alvain, S. & d'Ovidio, F. Can we detect oceanic biodiversity hotspots from space?. *ISME J.* **7**(10), 2054 (2013).

17. Lévy, M., Jahn, O., Dutkiewicz, S., Follows, M. J. & d'Ovidio, F. The dynamical landscape of marine phytoplankton diversity. *J. Royal Soc. Interface* **12**(111), (2015).

18. Kai, E. T. et al. Top marine predators track Lagrangian coherent structures. *Proc. Natl. Acad. Sci. USA* **106**(20), 8245 (2009).

19. Scales, K. L. et al. Fisheries bycatch risk to marine megafauna is intensified in Lagrangian coherent structures. *Proceedings of the National Academy of Sciences of the United States of America* **6**, (2018).

20. Della Penna, A., De Monte, S., Kestenare, E., Guinet, C., & d'Ovidio, F. Quasi-planktonic behavior of foraging top marine predators. *Sci. Rep.* **5** (2015).

21. Abrahms, B. et al. Mesoscale activity facilitates energy gain in a top predator. *Proc. R. Soc. B* **285**(1885), 20181101 (2018).

22. Prants, S., Budyansky, M. & Uleysky, M. Identifying Lagrangian fronts with favourable fishery conditions. *Deep Sea Research Part I Oceanogr. Res. Papers* **90**, 27 (2014).

23. Watson, J. R., Fuller, E. C., Castruccio, F. S. & Samhouri, J. F. Fishermen follow fine-scale physical ocean features for finance. *Front. Marine Sci.* **5**, 1 (2018).

24. Sato, M. et al. Coastal upwelling fronts as a boundary for planktivorous fish distributions. *Marine Ecol. Progress Ser.* **595**, 171 (2018).

25. Yoder, J. A., Ackleson, S. G., Barber, R. T., Flament, P. & Balch, W. M. A line in the sea. *Nature* **371**(6499), 689 (1994).

26. Woodson, C. B. & Litvin, S. Y. Ocean fronts drive marine fishery production and biogeochemical cycling. *Proc. Natl. Acad. Sci.* **112**(6), 1710 (2015).

27. Jansen, S., Klaas, C., Krägefsky, S., von Harbou, L. & Bathmann, U. Reproductive response of the copepod Rhincalanus gigas to an iron-induced phytoplankton bloom in the Southern Ocean. *Polar Biol.* **29**(12), 1039 (2006).

28. Druon, J.-N. et al. Satellite-based indicator of zooplankton distribution for global monitoring. *Sci. Rep.* **9**(1), 4732 (2019).

29. Kasumyan, A. The olfactory system in fish: structure, function, and role in behavior. *J. Ichthyol.* **44**(2), S180 (2004).

30. Boccaletti, G., Ferrari, R. & Fox-Kemper, B. Mixed layer instabilities and restratification. *J. Physical Oceanogr.* **37**(9), 2228 (2007).

31. Blain, S. et al. Effect of natural iron fertilization on carbon sequestration in the Southern Ocean. *Nature* **446**(7139), 1070 (2007).

32. Irigoien, X. et al. Large mesopelagic fishes biomass and trophic efficiency in the open ocean. *Nat. Commun.* **5**, (2014).

33. Cherel, Y., Fontaine, C., Richard, P. & Labatc, J.-P. Isotopic niches and trophic levels of myctophid fishes and their predators in the Southern Ocean. *Limnol. Oceanogr.* **55**(1), 324 (2010).

34. Kwon, E. Y., Primeau, F. & Sarmiento, J. L. The impact of remineralization depth on the air-sea carbon balance. *Nat. Geosci.* **2**(9), 630 (2009).

35. M. A. St. John *et al.*, A dark hole in our understanding of marine ecosystems and their services: perspectives from the mesopelagic community. *Front. Marine Sci.* **3**, 31 (2016).

36. Sala, E. et al. *Nature* **592**(7854), 397 (2021).

37. Béhagle, N. et al. Acoustic distribution of discriminated micronektonic organisms from a bi-frequency processing: The case study of eastern Kerguelen oceanic waters. *Prog. Oceanogr.* **156**, 276 (2017).

38. Kloser, R. J., Ryan, T. E., Young, J. W. & Lewis, M. E. Acoustic observations of micronekton fish on the scale of an ocean basin: potential and challenges. *ICES J. Marine Sci.* **66**(6), 99804 (2009).
39. G. Duhamel et al., Biogeographic Patterns Of Fish (Biogeographic Atlas of the Southern Ocean) chapter Biogeographic Patterns Of Fish, pages 328–362, Scientific Committee on Antarctic Research Cambridge UK (2014).

40. Saunders, R. A., Fielding, S., Thorpe, S. E. & Tarling, G. A. School characteristics of mesopelagic fish at South Georgia. *Deep Sea Research Part I Oceanogr. Res. Papers* **81**, 62 (2013).

41. Sutton, T. T. et al. A global biogeographic classification of the mesopelagic zone. *Deep Sea Research Part I Oceanogr. Res. Papers* **126**, 85 (2017).

42. P. Koubbi et al., Ecoregionalization of myctophid fish in the Indian sector of the Southern Ocean: Results from generalized similarity models. *Deep Sea Research Part II: Topical Studies in Oceanography* **58**(1), 170 (2011), Census of Antarctic Marine Life: Diversity and Change in the Southern Ocean Ecosystems.

43. Lara-Lopez, A. L., Davison, P. & Koslow, J. A. Abundance and community composition of micronekton across a front off Southern California. *J. Plankton Res.* **34**(9), 82803 (2012).

44. Netburn, A. N. & Koslow, J. A. Mesopelagic fish assemblages across oceanic fronts: a comparison of three frontal systems in the southern California Current Ecosystem. *Deep Sea Research Part I Oceanogr. Res. Papers* **134**, 80 (2018).

45. Bettencourt, J. H., López, C. & Hernández-García, E. Characterization of coherent structures in three-dimensional turbulent flows using the finite-size Lyapunov exponent. *J. Phys. A Math. Theor.* **46**(25), 254022 (2013).

46. Bettencourt, J. H., Rossi, V., Hernández-García, E., Marta-Almeida, M. & López, C. Characterization of the structure and cross-shore transport properties of a coastal upwelling filament using three-dimensional finite-size Lyapunov exponents. *J. Geophys. Res. Oceans* **122**(9), 7433 (2017).

47. De Monte, S. et al. Frigatebird behaviour at the ocean-atmosphere interface: integrating animal behaviour with multi-satellite data. *J. Royal Soc. Interface* **9**(77), 3351 (2012).

48. Koenker, R. Quantile Regression Econometric Society Monographs (Cambridge University Press, 2005).

49. Cade, B. S. & Noon, B. R. A gentle introduction to quantile regression for ecologists. *Front. Ecol. Environ.* **1**(8), 412 (2003).

50. Sankaran, M. et al. Determinants of woody cover in African savannas. *Nature* **438**(7069), 846 (2005).

51. R Core Team. *R: A Language and Environment for Statistical Computing*. R Foundation for Statistical Computing Vienna, Austria (2015).

52. Sergi, S. et al. Interaction of the antarctic circumpolar current with seamounts fuels moderate blooms but vast foraging grounds for multiple marine predators. *Front. Marine Sci.* **7**, 416 (2020).

53. Adler, J. Chemotaxis in bacteria. *Ann. Rev. Biochem.* **44**(1), 341 (1975).

54. Pakhomov, E., Perissinotto, R. & McQuaid, C. Prey composition and daily rations of myctophid fishes in the Southern Ocean. *Marine Ecol. Prog. Ser.* **1**, (1996).

55. Genin, A., Jaffe, J. S., Reef, R., Richter, C. & Franks, P. J. S. Swimming against the flow: a mechanism of zooplankton aggregation. *Science* **308**(5723), 860 (2005).

56. Vidal, J. Physioecology of zooplankton. I. Effects of phytoplankton concentration, temperature, and body size on the growth rate of Calanus pacificus and Pseudocalanus sp. *Marine Biol.* **56**(2), 111 (1980).

57. Alonzo, F., Mayzaud, P., Razouls, S., Bocher, P. & Cherel, Y. Seasonal changes in biomass, growth rates and production of subantarctic calanoid copepods in the Bay of Morbihan. *Kerguelen Islands Marine Biol.* **142**(3), 525 (2003).

58. Labat, J.-P., Mayzaud, P. & Sabini, S. Population dynamics of Themisto gaudichaudii in Kerguelen Islands waters. *Southern Indian Ocean Polar Biol.* **28**(10), 776 (2005).

59. Greely, T. M., Gartner, J. V. Jr. & Torres, J. J. Marine Biol. **133**(1), 145 (1999).

60. Robinson, L. M. et al. Pushing the limits in marine species distribution modelling: lessons from the land present challenges and opportunities. *Global Ecol. Biogeogr.* **20**(6), 789 (2011).

61. Abraham, E. R. The generation of plankton patchiness by turbulent stirring. *Nature* **391**(6667), 577 (1998).

62. Snyder, S., Franks, P. J. S., Talley, L. D., Xu, Y. & Kohin, S. Crossing the line: Tunas actively exploit submesoscale fronts to enhance foraging success. *Limnol. Oceanogr. Lett.* **2**(5), 187 (2017).

63. Pakhomov, E. A. & Froneman, P. W. Composition and spatial variability of macroplankton and micronekton within the Antarctic Polar Frontal Zone of the Indian Ocean during austral autumn 1997. *Polar Biol.* **23**(6), 410 (2000).

64. Ohman, M. D. Sources of variability in measurements of copepod lipids and gut fluorescence in the California Current coastal zone. *Marine Ecol. Progress Ser.* **42**(2), 143 (1988).

65. Sulman, M. H., Huntley, H. S., Lipphardt, B. & Kirwan, A. Leaving flatland: diagnostics for Lagrangian coherent structures in three-dimensional flows. *Physica D Nonlinear Phenomena* **258**, 77 (2013).

66. Morrow, R. et al. Global observations of fine-scale ocean surface topography with the Surface Water and Ocean Topography (SWOT) Mission. *Front. Marine Sci.* **6**, 232 (2019).

67. Escobar-Flores, P., O'Driscoll, R. & Montgomery, J. Predicting distribution and relative abundance of mid-trophic level organisms using oceanographic parameters and acoustic backscatter. *Marine Ecol. Prog. Ser.* **592**, 37 (2018).

68. A. Della Penna et al., Lagrangian analysis of multi-satellite data in support of open ocean Marine Protected Area design. *Deep Sea Research Part II: Topical Studies in Oceanography* **140**, 212 (2017), Future of oceanic animals in a changing ocean.

## Acknowledgements

This work is a contribution to the CNES/TOSCA and LEFE project KERTREND and BIOSWOT, and was partly funded by the Copernicus Marine Environment Monitoring Service (CMEMS) Sea Level Thematic Assembly Centre (SL-TAC). This work was supported financially and logistically by the Agence Nationale de la Recherche (ANR MyctO-3D-MAP, Programme Blanc SVSE 7 2011, Y. Cherel), the Institut Polaire Français Paul Emile Victor, and the Terres Australies et Antarctiques Françaises and benefited from ANR-11-IDEX-0004-17-EURE-0006. The authors would also like to thank Mara Baudena, Bettina Fach, Philippe Koubbi, Mark Ohman, Sara Sergi, Lars Stemmann and Jost von Hardenberg for their helpful advice, and the editor and the anonymous referees for reviewing the manuscript, which helped us to improve considerably the quality of our work.

## Author contributions

A.B., E.S.G. and F.d.O designed the research with assistance from X.C. A.B. performed the research. A.B. and D.D. analyzed the data. C.C. and Y.C. provided the data. A.B. wrote the paper, with assistance from all authors.

## Competing interests

The authors declare no competing interests.

## Additional information

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1038/s41598-021-94368-1.
**Correspondence** and requests for materials should be addressed to A.B.

**Reprints and permissions information** is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2021