The role of submesoscale currents in structuring marine ecosystems

From microbes to large predators, there is increasing evidence that marine life is shaped by short-lived submesoscales currents that are difficult to observe, model, and explain theoretically. Whether and how these intense three-dimensional currents structure the productivity and diversity of marine ecosystems is a subject of active debate. Our synthesis of observations and models suggests that the shallow penetration of submesoscale vertical currents might limit their impact on productivity, though ecological interactions at the submesoscale may be important in structuring oceanic biodiversity. Short-lived three-dimensional submesoscale currents, responsible for swirling ocean color chlorophyll filaments, have long been thought to affect productivity. Current research suggests they may not be effective in enhancing phytoplankton growth, but may have important contributions to biodiversity.

creates streaks much thinner than the spoon. Importantly, stirring alone only implies the redistribution of the biological constituents, which are transported by the currents, with no change in their concentrations. Thus in principle stirring should have no consequences on the productivity or diversity of marine ecosystems.
During the last decade or so, observations and models have revealed evidence for direct forcing at scales significantly smaller than mesoscale eddies [5][6][7] . These submesoscale motions (Box 2) are characterized by small vortices and a plethora of rapidly changing small-scale density filaments and fronts. The lateral spatial scales of these features are in the range of a few hundred meters up to a few kilometers 8 , about an order of magnitude smaller than the latitude-dependent deformation scale. Upperocean submesoscale currents are three-dimensional 9 : intense along front lateral currents, combined with a secondary vertical circulation. These currents may drive nutrients from the deep pool into the euphotic zone-and drive phytoplankton into the dark. Submesoscale processes generally also reduce mixed-layer depth, increase vertical stratification, and decrease vertical mixing, with consequences for the residence time of phytoplankton in the euphotic zone. Thus, in contrast to mesoscale stirring, submesoscale forcing can affect growth rates, biomasses, biogeochemical fluxes, and community structure.
Swirling phytoplankton features are often characterized by the presence of a local maximum (of species, taxa, or types), bounded by strong horizontal concentration gradients (Figs. 1 and 3). Because at least one of the two horizontal dimensions of such features falls in the submesoscale range, they are often referred to as submesoscale and we will use this terminology throughout the paper. We make a distinction here between submesoscale patchiness and submesoscale variability. Submesoscale features, or patches, will appear as coherent blobs, streaks, and whorls of a property as shown in Fig. 1. Variability, on the other hand, is the magnitude of variation of a property at a given spatial scale.
Disentangling whether submesoscale phytoplankton patches are generated by stirring of existing biological gradients, or are an active response to submesoscale physics, is therefore key to quantifying their global impact on marine ecosystems: are they nothing more than astonishing rearrangements, or do they reflect a complex response of the ecosystem to forcings that are difficult to observe, model and quantify? Submesoscale currents are continuously forming, moving, and dissipating over time scales ranging from days to weeks, making them particularly difficult to sample and model. Because biological response time scales are long enough to allow displacement of the biological response from the proximate physical forcing, submesocale forcings are even more difficult to relate to ecosystem dynamics. Satellite altimeters, used to provide global maps of upper-ocean currents, do not presently resolve velocity features smaller than about 100 km, and so entirely miss submesoscale currents. Numerical models require extremely fine grid resolution to explicitly resolve these dynamics, which is computationally challenging. And the logistical issues of sampling and modeling short time-scale, small spatial-scale dynamics are compounded by a lack of tools to quantify biological properties-particularly rates, such as growth and grazing-in a fast-moving frame of reference.
Submesoscale phytoplankton features form through an assortment of often-concurrent, highly interconnected mechanisms that drive planktonic ecosystem reactions, affecting diversity, competition, and marine food-web structure. We propose here a conceptual framework to sort the mechanisms responsible for the in a. c A false-color rendering of the data obtained from a two-dimensional inverse Fourier transform of the spatial spectrum from panel b, but with a randomized phase. The median radial spatial spectra of b and c are shown in d, where the spectrum from d has been divided by 2 to make it visible: the two spectra are statistically identical. It is clear that the coherent phytoplankton patches in a and b are a consequence of the phase relationships among the Fourier components making up the spectrum, rather than the relative magnitudes of the Fourier components (the spectral slope). Panels b and c have identical variability-their spectra are the same, but the patchiness in b forms coherent structures while it is random in c formation of submesoscale patchiness, guide our synthesis of observations and high-resolution models, and focus our discussion of their biogeochemical and ecological implications. We begin by analyzing the ways in which submesoscale dynamics drive active responses in phytoplankton growth rates-the central topic of this review. For example, strong submesoscale vertical velocities may bring limiting nutrients into the euphotic zone, causing a rapid, local increase in phytoplankton growth rates (Fig. 4b). We then explore how submesoscale biological patchiness can be generated by passive processes: the stirring and mixing of existing biological features by ocean currents (Fig. 4a). This occurs with phytoplankton blooms, for example, when stirring is faster than the local biologically driven rates of change. Finally, we discuss biological processes driven by a reaction to active and passive forcings (Fig. 4d). These biological reactions can be behavioral (i.e., swimming) or ecosystem processes (i.e., changes in community structure), and may themselves lead to biological features with distinct spatial scales. This analysis leads us to discuss why submesoscale dynamics might be significant for ecosystem diversity, but less important for marine productivity.
Active processes alter phytoplankton growth rates The temporal and spatial scales of submesoscale currents are often similar to those of phytoplankton growth time scales and patch scales, suggesting the possibility of close coupling of phytoplankton growth with submesoscale forcings. One feature of submesoscale dynamics that makes them particularly relevant to planktonic ecosystems is that they drive strong local vertical velocities at fronts 9-12 (Box 2). These vertical velocities may drive locally enhanced nutrient fluxes up into the euphotic zone, and pull phytoplankton down into the dark ocean interior.
Frontal dynamics at submesoscales is a topic of current research, and categorizing the various types of fronts is a necessary precursor to discussing the role of submesoscales in biology. Broadly, there are a range of fronts, with persistent fronts and ephemeral fronts representing the two extremes (Box 2). Persistent fronts include, for example, the intense horizontal density gradient associated with the Gulf Stream and Kuroshio (and many others), which have lateral scales on order 10-50 km, and which are relatively steady on timescales of a month or so. They are held in place by coastal boundaries and atmospheric forcing, have strong along-front currents, and large vertical velocities associated with a cross-frontal ageostrophic circulation that reaches well into the ocean's interior. Fast-changing, smaller, ephemeral (or, simply, submesoscale) fronts and filaments form through both secondary instabilities associated with forced, persistent fronts, and in the open ocean, through the generic process of straining by mesoscale currents and open-ocean submesoscale instabilities. These fronts evolve on timescales of days to weeks, and like their larger, forced cousins, have strong cross-frontal circulations that are out of geostrophic balance. However, the associated vertical velocities are typically trapped within the ocean's well-mixed surface layer. Upward nutrient supplies by submesoscale transport. In midlatitude and high-latitude, wintertime convective mixing brings nutrients into the euphotic zone, and these are consumed when light conditions become favorable, leading to strong seasonality in phytoplankton abundance with a pronounced springtime maximum; it is in only in summer that the environment is oligotrophic, i.e. that the euphotic layer is deprived of nutrients. In contrast, the subtropics are oligotrophic "ocean deserts"-despite having plenty of light year-round, because large-scale downwelling and shallow mixed layers limit nutrient supplies from the deep. In such nutrient-depleted oligotrophic environments, phytoplankton populations respond to upward nutrient inputs into the euphotic zone by submesoscale vertical velocities on time scales of hours to days, exhibiting locally increased growth rates and subsequent increased biomass 13 . Though sampling such transient and small-scale events is difficult, new technologies such as profiling floats are providing direct evidence of short-lived, submesoscale nutrient patches, presumably indicative of vertical nutrient fluxes at the same lateral scale [14][15][16][17][18] . Associations between frontal vertical velocities and chlorophyll patches have been observed in high-resolution surveys conducted in frontal regions 6,19 . In these studies, vertical velocities were derived from the quasigeostrophic omega-equation, and mapped onto chlorophyll distributions at horizontal resolution~20 km. Other   Fig. 4 Schematic representation of passive, active, and reactive mechanisms. a Passive deformations of existing phytoplankton patches by mesoscale currents can form submesoscale filaments. b Changes in phytoplankton growth rates driven by submesoscale nutrient fluxes into the euphotic zone can lead to the formation of phytoplankton patches at deep fronts. If the submesoscale circulations do not penetrate into the nutricline below the euphotic zone at shallow fronts, then it is unlikely that phytoplankton patches will form c. d Phytoplankton reactions (e.g., growth, competition, swimming, loss by grazing, etc.) and behavioral responses of higher trophic-level organisms (including fish, birds, and mammals) in response to passive and active processes studies have shown that specific submesoscale chlorophyll patches present in ocean color images were located over regions of strong stretching 20 or divergence (and thus of upwelling) 21 , identified through Lagrangian analyses of mesoscale surface currents derived from satellite altimetry. More recently, intense submesoscale chlorophyll features in ocean color images have been statistically connected to the presence of temperature fronts in the North Pacific subtropical gyre 22 , another piece of indirect evidence for a submesoscale phytoplanktonic growth response to nutrient upwelling at submesoscale fronts.
A key factor determining the magnitude of vertical nutrient fluxes is the depth of penetration of vertical velocities into the nutricline: vertical velocities confined to a mixed layer above the nutricline (such as those often associated with ephemeral fronts) do little besides mix the nutrient-poor surface layer 23 (Fig. 4c). On the other hand, deep, dynamic, persistent fronts that extend into the nutricline below the mixed layer can provide a nutrientflux pathway to the euphotic zone from the interior 24,25 , and so may be effective in altering budgets of biogeochemical tracers (Fig. 4b). This upward nutrient flux may reach the surface mixed layer to drive surface blooms, or, may only reach the base of the euphotic zone to fuel subsurface blooms 5,17 . Examples of deep and shallow fronts were examined in one submesoscale-resolving model of a generic oligotrophic subtropical gyre 26 . In this model, the vertical velocities associated with submesoscale dynamical fronts were largely confined to the surface mixed layer, contributing little to the vertical nutrient flux. In contrast, the vertical velocities associated with the quasi-permanent Gulf Stream-like front reached well below the mixed-layer and into the nutricline, driving strong vertical nutrient fluxes (Fig. 5).
Whether submesoscale-driven nutrient fluxes have a broader spatial relevance is a topic of current investigation. Biogeochemical models run with horizontal grid-cell sizes differing by twoorders of magnitude (typically from 1 to 100 km) have consistently shown that productivity increased with increased model resolution over the first months to years of model integration 5,[26][27][28][29][30] . This could be interpreted as support for the global importance of active responses to the submesoscale dynamics resolved by the high-resolution models. However, while improved model resolution allows stronger submesoscale fronts, it also leads to stronger persistent fronts that may account for much of the enhanced nutrient flux and increase in basinintegrated productivity (Box 3). We must therefore be careful to distinguish whether model simulations with different resolutions drive changes in productivity through the emergence of submesoscale fronts, or through the strengthening of persistent fronts.
Though they act at small scales, submesoscale dynamics appear to play a role in models of the basin-scale ocean circulation and tracer distribution. For example, comparing solutions from otherwise-identical high-resolution (1/54°) and coarseresolution (1/9°) models 26,31 showed that after more than 50 years the models diverged, with the high-resolution solution showing a deeper nutricline and shallower winter mixed layer in the subtropical gyre. This resulted, counter-intuitively, in less productivity in the subtropics, despite more submesoscale activity. This result illustrates a basin-scale ecosystem response to submesoscale forcing (productivity decrease) in opposition to the submesoscale response (locally increased productivity).
There is also a strong seasonal component to the structure and intensity of submesoscale vertical velocities at higher latitudes: the horizontal buoyancy gradients that underlie vertical velocities at submesoscale fronts are strongest in deep wintertime mixed layers 32 . However, over vast regions of the ocean, including the subtropical gyres, temperate zones, and high latitudes, wintertime mixed layers are already replete with nutrients through convective mixing. Thus in wintertime, submesoscale vertical velocities may not be effective in driving nutrient transport to the euphotic zone. In summer, shallow mixed layers are nutrient-depleted and phytoplankton growth rates can be enhanced by vertical nutrient fluxes. However, thin surface mixed layers create weaker summer submesoscale vertical velocities that may not penetrate down to the nutricline 8,32,33 . Satellite observations of chlorophyll in the North Pacific subtropical gyre tend to confirm the strong seasonality of the impact of submesoscale dynamics on phytoplankton biomass, showing a negligible effect in summer 22 . This annual cycle suggests that over large regions of the world's ocean, submesoscale-driven vertical transport could be out of phase with the seasonal nutrient requirements of phytoplankton.
The seasonality and shallow penetration of vertical velocities at submesoscale fronts bring into question their efficacy for upward transport of nutrients, particularly in comparison to persistent fronts. Nevertheless, modeling studies suggest that vertical transport at submesoscale fronts could be significantly enhanced in the presence of variable winds 34,35 , and that the presence or absence of high-frequency winds might be a critical factor influencing the depth reached by nutrient fluxes at fronts 36,37 . More systematic analysis is required to disentangle the relative contributions and scales of the fluxes associated with different varieties of front 38 , in presence or absence of winds.
Growth limitation by downward submesoscale transport. Submesoscale vertical velocities may also limit growth by moving phytoplankton and nutrients out of nutrient-rich environments 26,28,39 . The downward branch of frontal circulation cells that reach below the base of the euphotic zone may be especially effective in this regard 26,28,39 . In eastern boundary upwelling systems and equatorial regions, for example, nutrients are abundantly supplied to the euphotic zone by large-scale upwelling, sustaining elevated levels of phytoplankton growth. In such systems, model simulations support the suggestion that the observed reduction of biological production in regions of elevated eddy kinetic energy might be due to the downward leakage of both phytoplankton and nutrients from the nearshore to the open ocean 40,41 . In the North Atlantic, where seasonal productivity is characterized by the occurrence of a strong spring bloom, the observational and high-resolution modeling study of Omand et al. 42 suggests that the most favorable conditions for efficient export are during the transition from winter to spring.
Light exposure enhanced by submesoscale stratification. Several field and modeling studies have shown that submesoscale dynamics may also increase phytoplankton residence times in the euphotic zone of deep mixing layers 43 by tilting existing horizontal density gradients to increase the vertical stratification of surface mixed layers 7,44,45 . This vertical-shear-driven increase in vertical stratification reduces the local turbulence; in mid and high latitudes (with most observations from the North Atlantic) the increased residence time of phytoplankton in the euphotic zone leads to early local phytoplankton blooms compared to surrounding, actively mixing waters 7,46 . The horizontal scale of such changes in vertical turbulence is set by the submesoscale circulation; this leads to submesoscale patches of enhanced phytoplankton growth rates, which can ultimately form submesoscale patches of phytoplankton earlier than the regional average. Such early, local blooms may not change the seasonal productivity budget 47 but presumably affect the seasonal species succession and local community structure.
Finally, submesoscale vertical velocities also have the potential to penetrate into the deep chlorophyll maximum, displacing lowlight-adapted cells upward and high-light-adapted cells downward, altering the community structure at the base of the euphotic zone, and potentially altering the net growth rate there 5 .

Passive processes reorganize phytoplankton without affecting biomass
We now consider the complementary issue of the generation of submesoscale biological features through passive stirring. Stirring by larger-scale currents can readily generate horizontal concentration gradients at scales much smaller than the flow itself ( Fig. 4a), complicating the attribution of submesoscale biological patches to their forcing.
Surface stirring stretches phytoplankton into submesoscale patches. Consider an example of a few smooth mesoscale vortices a few hundred kilometers in diameter, with no dynamic features at smaller scales. A tracer patch in the surface mixed layer stirred by this current will be stretched out in one direction, and laterally compressed in the other by the strain between eddies (Fig. 4a). , with nutrient contours (black lines) and mixed-layer depth (heavy blue lines) overlaid. Vertical velocities associated with the persistent front (meandering Gulf Stream, large black arrows) extend from the surface and penetrate below the nutricline (delinated as the first nutrient contour). The near-surface vertical velocity field also exhibits thin, elongated submesoscale structures, located at submesoscale temperature fronts (small black arrows). Unlike the vertical velocities associated with the persistent front, the vertical velocities at submesoscale fronts do not extend much below the surface mixed layer, and do not always reach the nutricline. This results in a response of surface phytoplankton at the deep persistent front (yellow color show maximum phytoplankton concentrations at the front, large black arrows) but not at the shallow submesoscale front (small black arrows). Adapted from Lévy et al. 26 The thinnest width of a tracer filament will be set by a competition between the rate of compression/stretching, which thins the tracer, and the rate of horizontal diffusion, which broadens it. The resulting spatial patterns of chlorophyll will have smaller-scale variability than the flow itself. d'Ovidio et al. 48 (Fig. 2) used satellite altimetry to generate mesoscale surface velocity fields in an eddy-active region of the ocean. These smooth fields were used to advect virtual particles, representing the four broad groups of phytoplankton types that were identified from ocean color data in this region. The stirring stretched the initial large-scale surface patches of virtual phytoplankton types into submesoscale patches whose shapes were visually similar to those seen in satellite color and SST data, at least at the mesoscale. At submesoscales (below roughly 1/2°in Fig. 2), the correlation is more questionable. This may be due to the limited (~100 km) resolution of the underlying velocity field, or to inhomogenous biological reactions, which can alter the structure of the variance 49,50 . Similarities between the satellite data and the modeled surface phytoplankton distributions, however, suggest that stirring by smooth mesoscale currents may be a key player here, distorting the large-scale phytoplankton landscape into submesoscale patches delineated by sharp gradients.
The rearrangement of spatial gradients by stirring can sometimes be detected in the different slopes of the variance spectra (a Fourier transform of the squared spatial property anomaly) of the underlying flows and the property being stirred. Oceanic velocity power spectra often exhibit power-law behavior, with spatial variance proportional to k −a , where k is the wavenumber, and −a the spectral slope on a log-log plot. Flatter slopes (less negative, e.g., k −1 ) indicate more variance at small scales than steeper Box 1 A brief history of the mesoscale A discussion of submesoscale dynamics requires some historical and conceptual understanding of the mesoscales. Up through the 1950s, basic theory considered the ocean to be characterized by gyre-scale currents and meters-scale turbulence, with little in between. Oceanic "eddies" were first inferred from the erratic patterns of John Swallow's trackable floats 129 , and from time series of moored current arrays in the early 1960s 130 . These indicated the presence of time-variable currents with amplitudes relatively large compared to their mean values. It took another decade for a shadow of the spatial structure of these disturbances to be revealed by the USSR POLYGON experiment and the US Mid-Ocean Dynamics Experiment MODE 131 . Both of these used tightly coordinated campaigns of ship-based observations and current-meter arrays to piece together a picture of "mesoscale" dynamics; both programs focused on relatively small regions in the North Atlantic. They revealed elliptical eddies with lateral scales of about 100-200 km and time scales of months. As first explained by Gill et al. 132 , the spatial and temporal scales of mesoscale eddies are consistent with baroclinic instability (also the dominant source of midlatitude weather systems) which converts the massive potential energy stored in the ocean's sloping density surfaces-an energy source 1000 times as large as the kinetic energy of oceanic currents-into turbulent eddies. The dominant lateral scale of baroclinic instability is near the Rossby deformation scale, L D = NH/f, where N is the buoyancy frequency, f is the Coriolis frequency (which vanishes at the equator and is maximum at the poles), and H is roughly the depth of the thermocline. This is the scale at which rotational and buoyancy effects play equally important roles in the dynamics. Typical values of L D range from tens of kilometers in the polar oceans to hundreds of kilometers in the tropics. Satellite altimetry, beginning with the short-lived SEASAT program in 1978, revolutionized oceanography, enabling synoptic views of entire ocean basins 133 . Altimetric measurements show the spatial patterns of sea-surface height, from which currents can be estimated via geostrophic balance: on a rotating planet horizontal pressure gradients are approximately balanced by the Coriolis force, causing currents to flow around elevated regions (clockwise around elevated regions in the Northern Hemisphere, and so forth). As technology improved, a global synoptic view of the mesoscale ocean at last emerged with the TOPEX/Poseidon program in the 1990s 134 . Observations from this satellite revealed eddies nearly everywhere (as seen, for example, in these two plots of AVISO sea-level anomaly from January 31, 2013 over the Atlantic and Pacific oceans, image credit: AVISO/CNES/CLS), with scales of order 100 km, decreasing with increasing latitude, roughly proportional to but larger than the latitudinal scaling of L D . However, ocean variability on smaller scales is difficult to infer from current satellite-based altimetric instruments, which have an effective resolution near the deformation scale-about 100 km. Features seen in satellite altimetry are intrinsically dynamic, thus the absence of such global measurements at smaller scales is a severe limitation on our ability to separate dynamic from kinematic effects in submesoscale structures. Submesoscale motions, with lateral scales a decade or so below the mesoscale, are best defined in terms of dynamical balances 11 : unlike mesoscale dynamics, submesoscale dynamics are not well-approximated by geostrophic balance, and so all the terms in the equations of motion are potentially important. Submesoscale dynamics are often closely linked with fronts. The purpose of this box is to give a brief overview of the physical forcings that lead to fronts and submesoscales, and the phenomenology of the flows themselves. Early evidence of dynamical features much smaller than the mesoscale came from sun-glitter in photographs taken during the 1970s Apollo Mission, which revealed predominantly cyclonic spiral eddies in the ocean, with scales of around 10-25 km 135 . All such submesoscale motions derive their energy from lateral density gradients in the surface ocean (see schematic). These density gradients are driven by atmospheric forcing: heat and freshwater exchanges, and winds, which inhomogeneously mix the upper ocean. Horizontal straining by mesoscale eddies (red arrows) can squeeze these density contrasts into nearly vertical planes (gray surfaces), intensifying their geostrophic along-front, sheared currents (green arrows)-this is an example of frontogenesis. Fluid instabilities feed off these localized sources of potential and kinetic energy, generating submesoscale vortices and filaments that carry of some of this energy away. A sufficiently strong front will generate a cross-frontal ageostrophic secondary circulation-an overturning circulation directed in the sense of trying to flatten the density surfaces in the front (yellow arrows). Next we describe persistent fronts, submesoscale fronts, and a more complicated scenario that inextricably links mesoscale and submesoscale motions. Persistent fronts like those associated with the Gulf Stream and Kuroshio are locked in place by the coastal boundary and large-scale atmospheric forcing. Strong winter-time winds are directed nearly parallel to the axis of the jet, driving cold dense water southward across the front through Ekman forcing-a dynamical balance between surface friction, pressure, and Coriolis acceleration that moves upper ocean waters to the right of the wind in the Northern Hemisphere-keeping the density surfaces steep, the stratification low, and strengthening the jet 24,136 . This forcing is directly balanced by submesoscale symmetric instability, which takes energy mostly from the kinetic energy of the jet, and baroclinic instability, which converts the potential energy of the sloping density surfaces into large meanders and eddies. Associated with these highly energetic processes is a deep vertical velocity structure 137 that reaches into the thermocline. The ubiquitous submesoscale density fronts and filaments that are continuously created at the ocean's surface by mesoscale strain differ from persistent fronts like the Kuroshio and Gulf Stream. Unlike their more powerful cousins, submesoscale fronts are continuously forming, moving, and dissipatingthey are not in a steady-state balance between forcing and dissipation at a fixed location. These fronts have associated cross-frontal secondary circulations that are generally confined to the vertically well-mixed upper layer of the ocean 138 . In the presence of shallow summer-time mixed layers, these features do not carry much energy and thus have little impact on the local circulation. In the winter, on the other hand, the nearly vertical density surfaces associated with these fronts and filaments extend through the deep mixed layer, and so contain significant potential energy in addition to the intense kinetic energy of their cross-front currents. In the last few years, ample evidence 8,[139][140][141] indicates that this energy is converted to kinetic energy by mixed layer instability 142, 143 : a submesoscale baroclinic instability with maximal growth at the mixed-layer deformation scale, L ML = N ML h ML /f, where N ML and h ML are, respectively, the mixed-layer buoyancy frequency and mixed-layer depth. In mid-latitudes in winter, with N ML~1 0 −3 s −1 , h ML~2 00 m and f~2 × 10 −4 s −1 , this gives L ML~1 km, and like mesoscale baroclinic instability, the result is the formation of submesoscale vortices with sizes close to this scale 143 . Sufficiently steep density contrasts confined to the mixed layer can generate significant submesoscale energy, as above. But less-steep density surfaces that intersect the surface and reach deep into the interior can also generate significant energy at submesoscales, through a process called Charney-type baroclinic instability, which is likely widespread in the ocean 144 . This mesoscale-driven process provides an alternate but potentially important pathway between the surface and interior 23,25 . In all cases, submesoscale processes act against fronts, reducing mixed-layer depth, increasing stratification, and decreasing vertical mixing 145 . The eddy kinetic energy produced in these processes likely enhances lateral turbulent mixing at the submesoscale 146  slopes (more negative, e.g., k −3 ). The spectrum of passive tracer variance, with no injection or removal and initially spread over a region large relative to the eddy scale, will depend on the spectrum of the velocity field that is stirring it. Mesoscale flows with little variance at submesoscales (thus a steep spectrum) will tend to produce tracer variance spectrum with a slope near k −1this is the well-known Batchelor cascade 51,52 . Counter-intuitively, a velocity field with a great deal of submesoscale variability (Box 2) will produce a tracer field with relatively less submesoscale variance than a mesoscale velocity field would, resulting in a steeper tracer variance spectrum with slopes closer to k −2 (e.g., refs. 53,54 ). Thus, under the assumption of weak biological reactions, spectral slopes of the variance of biological tracers around k −1 could be indicative of stirring by mesoscale flows, while slopes around k −2 might indicate local stirring by submesoscale flows. Spatial spectra of phytoplankton have been generated from data gathered through 1D in situ transects, 2D satellite ocean Box 3 Capturing submesoscale dynamics with Ocean General Circulation Models (OGCMs) OGCMs are computer codes that estimate solutions to the oceanic primitive differential equations of fluid motion and thermodynamics discretized on a spatial grid that covers the word's ocean. These models form the ocean core of Earth System Models (ESMs). Their computational cost depends on the number of grid points. A decrease of the horizontal grid size by a factor of N, equivalent to an increase in the number of grid points by a factor of N 2 , increases the computational cost by N 3 , because higher resolution also requires smaller time steps 31 . Since the 1970s when OGCMs were first developed, their horizontal grid resolution has always been constrained by the capacity of supercomputers. This capacity has increased tremendously over the last decade. While ocean grids used for the 5th Coupled Model Intercomparison Project (CMIP5) of the World Climate Research Program (WCRP) ranged between 1°and 2°, ocean grids with~1/10°resolution are emerging in CMIP6. Mesoscales and submesoscales are thus smaller than the grid size of current ESMs and parameterizations are used to mimic their large-scale effects 143,148 . Reducing uncertainties linked to such subgridscale processes is one of the biggest challenges in ocean modeling 149 . Nevertheless, explicit resolution of submesoscale processes is achievable with regional or idealized models, which use the same computer codes as OGCMs but with reduced domains, allowing for finer grid cells. One example taken from Lévy et al. 150 is reproduced here: the model equations are solved over a rectangular ocean basin of~20°× 30°with closed boundaries. The circulation and thermodynamics are forced by seasonally varying winds and heat fluxes at the sea surface. The model equations are integrated on three different horizontal grids of resolution 1°, 1/9°, and 1/54°, and with the same vertical resolution of 30 vertical layers, 10 of which are in the upper 100 m. Common features emerged from the three model solutions including a strong persistent front (marked by the white SST contour) extending toward the northeast, separating a warm subtropical gyre from a cold subpolar gyre. Strong horizontal turbulence emerged as the resolution was increased. At 1/9°, the turbulence appeared as the meandering of the persistent front and in the formation of a few mesoscale eddies. At 1/54°, sharp submesoscale temperature fronts were present and the meandering of the persistent front was more pronounced. In addition, the persistent front was progressively displaced southward when the resolution increased. Thus, resolving small-scale ocean processes can lead to substantial changes in large-scale patterns. Importantly, the vertical velocity field close to the ocean surface was strongly modified both in its structure and intensity as the resolution increased. While at coarse resolution the main patterns in vertical velocity were coastal upwellings with subsidence in the center of the subtropical gyre, energetic fine structures in vertical velocity became more prominent as the resolution increased. These changes in SST and vertical velocity were also reflected in their spectral slopes (red and blue curves show the spectra at 1/54°and 1/9°, respectively), with much flatter spectra at 1/54°than at 1/9°resolution. One might question to what extent submesoscale processes are adequately represented at 1/54°resolution. One indication is provided by the power spectra, which show a sharp decrease of variance at scales smaller than 10 km (gray shading). This is the dissipative range for this model resolution: the effective model resolution is closer to 1/9°when the grid resolution is equal to 1/54°. Thus it appears that a submesoscale-resolving grid enables much better resolution of the mesoscale dynamics, providing the straining and stirring environment favorable for the emergence of submesoscale dynamics. This is clearly seen in the spatial spectra in the 10-100 km range: there is more energy with a 1/54°resolution than with a 1/9°resolution. The submesoscale dynamics, however, are still strongly damped at 1/54°resolution, which explains why this resolution range is often referred to as submesoscale permitting rather than submesoscale resolving. Recent studies suggest that an order of magnitude increase in resolution is required to adequately capture the full strength of the submesoscale dynamics, both in the vertical 8 and horizontal 151 dimensions. This resolution is not yet achievable with present-day computer capacities for regional biogeochemical studies; improved understanding will continue to depend on technological progress in high-performance computing. color data, or model results, over length scales ranging from meters to 100s of kilometers 1,55,56 . Phytoplankton spectra are often found to exhibit power-law scaling, with spectral slopes between k −2 and k −1 over mesoscale to submesoscale length scales (e.g., refs. 1,56-60 ), suggesting that passive stirring can create submesoscale phytoplankton variability by deforming large-scale regional gradients 3,61 . This is possible because the persistence time scale of phytoplankton surface blooms, 3-6 months on average 61 , is similar to the physical time over which stirring transfers variance to small scales 62 . A complexity, however, is that phytoplankton variability can also be created at small scales through active and reactive processes, yet still exhibit power-law scaling of spatial variance 49,50,63 . This complexity is reflected, for instance, in the strong temporal and regional variability of spectral slopes of ocean color variance 64 ; assigning particular mechanisms to specific values of plankton spectral slopes is therefore often ill-advised. Moreover, an important caveat is that the spectrum alone does not contain information about the coherence of the underlying patchiness (Fig. 1). It is the phase relationships among the Fourier components of the spectrum that determine the presence of coherent structures-patches-of properties, such as chlorophyll concentration in the ocean 65,66 .
Identical spectra can give a very patchy distribution, or a rather diffuse distribution of the property (Fig. 1), depending on the relative phases of the Fourier components.
Subsurface stirring forms subsurface phytoplankton layers.
The introduction of new instruments and sampling platforms, such as towed undulating vehicles (e.g., SeaSoar), gliders, and rapid CTD profilers (e.g., Moving Vessel Profiler) has begun to reveal the richness of the three-dimensional structure of biological properties in situ 6,67,68 . For instance, undulating glider transects in the vicinity of the intense upwelling located off the southern coast of Peru 69 revealed several submesoscale tongues of high-salinity, high-fluorescence water extending downward from the euphotic zone (Fig. 3). The coincidence of the salinity and fluorescence supports the notion of a passive stirring mechanism in forming the layers. When plotted using density as the vertical axis, it becomes clear that these layers extend vertically across isopycnals, and tilt horizontally, consistent with the tilting and stretching of existing patches by a cross-frontal vertical shear 70 . These phytoplankton layers formed patches whose edges were well defined by their spatial gradients (Fig. 3); the patches had 2-10 km horizontal scales across the front, and persisted for weeks. The presence of high fluorescence in these layers well below the euphotic zone suggests that they originated near the surface, and became subducted at the front. Though the orientation and vertical gradients of chlorophyll in the submesoscale patches observed by Pietri et al. 69 suggested an origin in the euphotic zone, submesoscale subsurface patches have also been observed originating well below the surface. The formation of a subsurface thin layer, far example, was inferred from an initial phytoplankton/salinity patch at~125 m depth, with cross-patch and along-patch horizontal scales of~100 m and 75 km, respectively 71 . Vertical shearing of the horizontal patch gradients by low-frequency internal waves displaced the top and bottom of the patch relative to each other, forming a layer that tilted and thinned from 25 to 12 m over the course of~20 h. Though the origin of this low-salinity/high-fluorescence layer is unknown, it is clear that stirring and shearing by the ambient flow created submesoscale horizontal and vertical structure from larger-scale gradients of the original phytoplankton patch 72 .
It is generally difficult to tease apart the relative contributions of along isopycnal stirring vs. purely horizontal stirring in forming subsurface patches. As with horizontal stirring, by deforming larger-scale gradients, the downward branch of the circulation at submesoscale dynamic fronts (Box 2) will create smaller-scale vertical structure in the tracer 40 . These downward motions will deform tracer gradients-temperature, salinity, nitrate, and phytoplankton-generating interleaved layers. Such interleaved layers in cross-frontal sections, with slopes often steeper than those of isopycnals [73][74][75][76] , can also form through horizontal deformations of existing plankton patches that extend subsurface, with no need to invoke vertical velocities.
Reactive processes drive changes in community structure While submesoscale phytoplankton patches can be formed by active increases in growth rate resulting from locally favorable environmental conditions, or passive stirring of existing gradients, not all members of the community respond identically to the underlying physical forcing. We would thus expect to observe altered community structure. These reactions could result from differential growth rates or different behaviors (swimming, sinking, or floating) within the phytoplankton community, or from gradients in mortality rates as predators respond to the phytoplankton distribution, and could thus propagate up the trophic chain to zooplankton and larger predators.
Phytoplankton diversity is affected by active and passive processes. One useful measure of the reaction of planktonic communities to physical and biological forcings is their diversity. In terrestrial ecology, diversity is usefully divided into α-diversity and γ-diversity: α-diversity refers to the number of species within a local habitat, while γ-diversity describes the total number of species across multiple habitats in the region. In the ocean it can be difficult to define habitat and region. Here we will consider habitats to be waters that contain similar communities-typically determined by hydrographic properties (temperature, salinity), light availability, nutrient concentrations, and nutrient fluxes. For example, the colored patches in Fig. 3 delineate distinct habitats. In our context a region will usually refer to a frontal region, where different habitats can be brought together through stirring, or created locally by physical processes. For example, the ensemble of the colored patches in Fig. 3 represents a region. Thus, with the above definition, α-diversity will usually be found in submesoscale patches, while γ-diversity occurs at the mesoscale and larger.
Some studies have found that high planktonic γ-diversity in phytoplankton patches in frontal regions simply represents the abutment of the communities found in the water masses making up the front 48,77-79 . However, submesoscale nutrient pulses might be expected to give rise to submesoscale community patchiness (Fig. 2b, d), enhancing α-diversity and γ-diversity over that expected from purely passive stirring of regional communities, with consequences for trophic dynamics, and export fluxes. Short-term (days) responses of the nutrient-limited planktonic community to these nutrient pulses often include elevated phytoplankton growth and grazing rates, and a shift of the community toward dominance by larger phytoplankton-usually diatoms. Many studies have indeed found unique planktonic communities in submesoscale patches at fronts 6,77-82 . Distinct coastal and oceanic cyanobacterial communities, for example, were separated by a front in the Southern California Current System 81 , with the cyanobacterium Synechococcus dominating in the colder mesotrophic waters, and Prochlorococcus dominating the oligotrophic offshore waters. However, a patch with a unique phytoplankton community dominated by diatoms was found in the frontal waters. This community of large phytoplankton was coincident with increased local nitrate fluxes and decreased microzooplankton grazing 83 -all indicative of physiological and trophic responses to the frontal physical dynamics that enhanced the γ-diversity through the creation of a new, submesoscale frontal habitat.
Exploring a large number of fronts, a model including 100 phytoplankton types embedded into a submesoscale-resolving flow suggested more complex changes 84 . Overall, the model generated statistically larger α-diversity at fronts than elsewhere (Fig. 6a). However, this was not systematic, with α-diversity varying among fronts, and some fronts being even less diverse than their surrounding waters. These results showed that not all fronts were efficient at driving nutrient fluxes or allowing the exploitation of new habitats. Interestingly, while in this model the fronts emerged as α-diversity and γ-diversity hot spots, stirring could in some instances decrease γ-diversity 85 by causing the local extinction of rare species with very specific environmental niches (Fig. 6b and c).
Plankton behavioral response interacts with active and passive processes. Phytoplankton exhibit a variety of behaviors that allow them to move relative to the water, from vertical migration to avoid predators 86 or to gain access to light or nutrients 87 , to floating or sinking in response to nutrient stress 88 . Phytoplankton swimming behaviors can interact with the 3D flows at submesoscale fronts to generate accumulations of plankton that form submesoscale patches 89,90 . Such patches form when horizontal velocity gradients sweep the organisms toward convergence zones, where their vertical swimming (and floating) behaviors cause them to accumulate 91 . For floating organisms, these patches form within submesoscale downwelling regions; the surface patches have subsurface extensions that form stronger horizontal gradients with stronger floating speeds 91 .
Predator-prey interactions are affected by active and passive processes. Stirring 4 , plankton motility (e.g., diel vertical migration), diffusion 63 , and active nutrient fluxes 92 can all decorrelate phytoplankton and zooplankton distributions and lead to nonlinear changes in the grazing of phytoplankton by zooplankton [93][94][95] . Similar to phytoplankton, gradients of zooplankton are often (but not always 96,97 ) associated with fronts. Composite analyses of hundreds of glider transects across fronts in the California Current System, for example, showed persistently elevated acoustic backscatter (proportional to zooplankton abundance/biomass) on the cold sides of the fronts 98 . The abrupt change in acoustic backscatter would seem to be indicative of abutment of two dissimilar zooplanktonic communities at the fronts, suggestive of increased γdiversity simply through passive stirring of distinct habitats. However, taxon-resolving sampling has shown distinct that submesoscale zooplankton communities form at fronts 96,99-101 , presumably in reaction to the frontal dynamics, indicating frontenhanced zooplankton α-diversity.
In a submesoscale patch in the North Pacific subtropical gyre, lateral stirring led to a large-scale diffusion which decreased the encounter rates of phytoplankton and zooplankton, thereby reducing the grazing pressure on the phytoplankton 102 . In general, however, the magnitude of these non-linear interaction effects remains poorly quantified; they are particularly difficult to assess from observations, requiring a synoptic description of phytoplankton and zooplankton distributions, together with spatially and temporally resolved measurements of their growth, grazing, and mortality rates 103 .
Passive and active forcings propagate up the trophic web. The enhanced abundances of eggs, nauplii, and copepodids sometimes observed at fronts 99,100 are indicative of locally enhanced zooplankton production, presumably in response to enhanced phytoplankton growth fueled by increased nitrate fluxes at the front 83 . Though identifying the proximal cause of the enhanced nitrate fluxes is difficult, it is clear that such forcings propagate through the ecosystem, creating distinct local communities with enhanced α-diversity and γ-diversity in the frontal region.  6 α-plankton and γ-plankton diversity. a The distribution of phytoplankton α-diversity over fronts (green) and away from fronts (black). Higher number of types demonstrates the increased coexistence of phytoplankton at fronts. Shown are outputs from the Darwin model comprising 100 phytoplankton types embedded in the high-resolution submesoscale flow shown in Box 3. In this model, the community structure is not imposed; rather it "self-assembles" according to the relative fitness of the phytoplankton types. The consequence is that at each time step, each model grid point is characterized by a specific community of plankton and α-diversity. b The abundances of the phytoplankton types that emerged over the entire region with the same simulation (with stirring) and c in a simulation without stirring. Stirring was suppressed by embedding the phytoplankton model in the annual time-mean flow instead of the submesoscale flow. Without stirring, the community contained more different types with more even abundances, thus the γ-diversity was larger.
Thus at a local level there were more types over fronts (more α-diversity); at the regional level the stirring decreased the number of coexisting types (less γ-diversity). Adapted from Lévy et al. 84,85 However, while plankton mainly drift with the ocean currents, larger organisms are motile and might be predicted to have a weaker relationship with local physical forcing. Nevertheless, thanks to the development of miniaturized animal tracking devices, associations of top predators with fronts have been reported for many taxa (for example basking sharks 104 , albacore tuna 105 , whales 106 , elephant seals 107 , frigatebirds 108,109 , gannets 110 ). The reasons why predators target frontal structures remain unclear and could have a variety of origins, including foraging and entrainment effects 62 , and bioenergetics 105 . One impediment to a better understanding of how submesoscale patchiness extends up the trophic web is the lack of information on the submesoscale distribution of intermediate-level fish that prey upon zooplankton and larval fish, and are foraged by top predators. Frequency-differencing echograms have revealed spatial coherence between top predators and forage fish in the submesoscale range 111 , and open the way to a better understanding of the reactive responses of ecosystems at fronts.
Given that it would likely take days for larger organisms to find and exploit physical forcings, persistent fronts lasting weeks and months present a biologically predictable ecotone. However, given the relative transience of submesoscale events, it remains unclear whether such features would be predictable enough (in an evolutionary sense) to contribute significantly to global ocean biological dynamics.

Discussion
Field, laboratory, and modeling studies have provided extensive evidence that submesoscale oceanic currents affect not only phytoplankton concentrations, but also plankton diversity and the entire marine trophic web. Still, the spatial scales of submesoscale phenomena pose an ongoing, observational challenge, particularly when it comes to sampling biodiversity or biogeochemical fluxes. In addition, due to computational limitations, the grid resolution of global models is too coarse to assess the contribution of submesoscale dynamics to biogeochemical budgets. Today, we can only speculate on how much these dynamics affect global biogeochemical cycles, marine productivity, and ocean biodiversity. The categorization of the processes into active, passive, and reactive used in this review, although imperfect, is intended to provide a conceptual framework to guide our speculation and will hopefully open the path for more quantitative approaches in the future. Based on our synthesis, we conclude that submesoscale dynamics might not be as important as previously thought for marine productivity, but might be more important than anticipated for plankton diversity and trophic web dynamics. We want to stress that our arguments, discussed in the following, are more qualitative than quantitative.
Submesoscale vertical nutrient transport may not be an important term in nutrient budgets. While passive stirring mainly redistributes constituents, active dynamics have the potential to strongly affect biogeochemical fluxes. This raises the question of how effective submesoscale vertical transport processes are in mediating net fluxes of nutrients within the euphotic zone. Early modeling studies of vertical nutrient transport at fronts explored persistent fronts where vertical velocities reached hundreds of meters below the surface, bringing them well into the deep nutrient pool. Such studies suggested that fronts could make a large (20-50%) contribution to total biogeochemical fluxes 5,27,112 . In addition to this modeling bias toward deepreaching persistent fronts, there may have been a bias toward in situ observations at such fronts because they are more predictable than submesoscale fronts.
A prevalent feature of submesoscale dynamics is their perturbation through mixed-layer instability (Box 2). These instabilities are strongest in deep mixed layers, found during winter at mid-latitude and high-latitude 32 , and their associated vertical velocities tend to confined to the mixed layer 25 . Such motions are thus strongest when phytoplankton growth could be limited by the lack of light through short residence times in the euphotic zone, and weak later in the season when growth is nutrient-limited; their contribution may be significant during the transition period at the onset of the spring bloom. While this is just one of many possible dynamical processes that characterize the submesoscale, it suggests that the combination of shallow submesoscale vertical fluxes and their phasing relative to phytoplankton growth cycles over vast open ocean oceanic regions may make them inefficient at driving reactive biological responses, and thus less important in global budgets than persistent fronts.
Quantifying the importance of vertical nutrient transport at ephemeral and persistent fronts will require better global sampling and statistics. Unfortunately, the surface detection of fronts -which is feasible at the global scale with sea-surface temperature and altimetry satellite data 113 -is not sufficient to quantitatively evaluate their potential for enhanced nutrient fluxes. Regional surveys using moored arrays, gliders, and highresolution towed instruments, such as the OSMOSIS project 114 have greatly improved our understanding of submesoscale processes; however, one cannot confidently extrapolate the results of a North Atlantic study to the global ocean. The large array (3800 as of today) of profiling Argo floats that populates the word's ocean providing continuous monitoring of the surface ocean vertical structure is a first step; however, this array undersamples both mesoscale and submesoscale features.
Passive and active processes shape planktonic biodiversity. Passive processes such as mesoscale stirring do not, in principle, change planktonic α-diversity, as habitats are simply advected and distorted by the ambient flow 48 . However, mesoscale stirring can enhance γ-diversity by bringing formerly geographically distinct habitats and communities close together at frontal regions, increasing the local number of habitats. This simple consideration has started to allow us to detect γ-diversity hotspots using ocean color satellites 115,116 . In addition, the presence of unique communities at fronts-which are sites of enhanced submesoscale physical forcing-may also reflect the effects of submesoscale dynamics on phytoplankton biomass and community structure at fronts. Such active responses of the planktonic community to submesoscale transports will affect both αdiversity and γ-diversity by creating new habitats at the fronts. Again, because of the logistical issues mentioned earlier, most in situ evidence of patterns in α-diversity at fronts comes from measurements at persistent fronts with strong nutrient fluxes. The deployment of continuously sampling underway systems offers the potential to better explore subsurface submesoscale patchiness in planktonic α-diversity.
Investigations of the relative roles of the different processesstirring, mixing, nutrient supply, stratification, swimming-in shaping planktonic biodiversity at the submesoscale are very preliminary, and are areas of active research. Nevertheless, two characteristics lead us to speculate that these processes play crucial roles in shaping biodiversity over much broader spatial and temporal scales than the immediate frontal region. The first characteristic is the constant submesoscale perturbation of marine ecosystems: organisms that succeed in this environment must be adapted to these rapid fluctuations. Thus the environmental variability is a key characteristic, and determinant of the community structure. Geographical variations in eddy stirring 117 and associated horizontal diffusivities 118 show strong contrasts over broad (>1000 km) regions, suggesting a similar scale of variability for ecosystem diversity. The second characteristic is that stirring and mixing of communities that were initially apart sets them in competition with one another. When the mixing time scale is faster than the competitive exclusion time scale, such stirring may be expected to increase α-diversity and γ-diversity. The degree to which different habitats are being constantly brought together, mixed and modified by the rapidly evolving submesoscale flow will determine the outcome of species competition over longer time and larger space scales.
Future directions. To distinguish among active, passive, and reactive processes requires high-frequency, high-resolution, 3D sampling of physical, chemical, and biological variables. Measuring down to submesoscale resolution in the field requires significant infrastructure, such as multiple ships, gliders, or large swarms of drifting floats 54,119 .
The greatest physical oceanographic challenge is to estimate submesoscale vertical velocities. Direct measurements are difficult due to the presence of inertial motions and internal waves at similar scales. Indirect estimates rely on approximations; these provide limited insight into the vertical extent of submesoscaledriven fluxes that determines their biogeochemical effectiveness 120 . The greatest biological oceanographic challenge is to quantify community composition, physiological rates, and trophic interactions with 3D submesoscale resolution.
Because fronts are regions of strong horizontal shear, vertical transport and horizontal stirring are strongly linked in the upper ocean. Thus, through this phase relationship, passive, active, and reactive processes are highly interconnected. It follows that it is difficult to determine whether a submesoscale phytoplankton patch formed as an active response to submesoscale flow or is simply the result of stirring of a larger patch that formed in response to larger-scale processes. An improved understanding of submesoscale processes will thus rely on our ability to reconstruct the Lagrangian histories of submesoscale patches from observations that are often Eulerian 62,119,121-124 . Real-time estimates of stirring from satellite altimetry can be useful for guiding adaptive sampling strategies during ship-based field campaigns to observe the evolution of tracer fields 125 . Combined Eulerian/Lagrangian sampling schemes in the field, such as making measurements while following drifters, along with supplementary transect information 126 , can help in tracking the evolution of biological communities as they are advected away from the submesoscale physical features that formed them. The anticipated 2021 launch of the new, high-resolution Surface Water and Ocean Topography (SWOT, https://swot.cnes.fr) satellite altimeter mission may provide a significant new remote-sensing tool for resolving Lagrangian trajectories with a better resolution than current altimeters 127,128 .
The contributions of submesoscale physical forcing to marine ecosystem dynamics remain unconstrained and uncertain. We are confronted with a problem at the interface of two fields-ocean physics and ecology-that lies at the frontier of knowledge in both fields, with new interaction mechanisms still being identified and their relative importance being evaluated. The strong vertical velocities associated with submesoscale fronts in the upper ocean would appear to position them to have profound biological effects. However, the vertical penetration and seasonal phasing of such dynamics may limit their efficacy in stimulating local phytoplankton growth. The considerable developments in modeling and instrumentation that have occurred over the last decades have helped in quantifying the occurrence of submesoscale biological variability. However, we will require coordinated, well-designed interdisciplinary sampling programs to disentangle the physical and biological dynamics, and in particular, the relative contributions of mesoscale and submesoscale processes. These programs must investigate both the immediate, local biological responses to submesoscale physical forcing, as well as how such responses emerge to affect global processes, such as biogeochemical cycles and biological diversity.