Enhanced mixing across the gyre boundary at the Gulf Stream front
- aDepartment of Atmospheric and Oceanic Science, University of Maryland, College Park, MD 20742;
- bEarth System Science, Stanford University, Stanford, CA 94305;
- cDepartment of Estuarine and Ocean Sciences, University of Massachusetts Dartmouth, North Dartmouth, MA 02747;
- dDepartment of Applied Mathematics and Theoretical Physics, University of Cambridge, Cambridge CB2 1TN, United Kingdom;
- eApplied Physics Laboratory, University of Washington, Seattle, WA 98105;
- fSchool of Earth and Ocean Sciences, University of Victoria, Victoria, BC V8P 5C2, Canada;
- gCollege of Earth, Ocean, and Atmospheric Sciences, Oregon State University, Corvallis, OR 97331
See allHide authors and affiliations
Edited by Carl Wunsch, Harvard University, Cambridge, MA, and approved June 11, 2020 (received for review March 24, 2020)

Significance
The Gulf Stream front separates the North Atlantic subtropical and subpolar ocean gyres, which have very distinct physical, biological, and chemical properties. Exchange across the front is believed to be important for both ocean circulation and biological productivity; however, the physical mechanisms responsible are not fully understood. Here, we present direct observational evidence of enhanced mixing across the north wall of the Gulf Stream. High-resolution simulations suggest that the presence of the sharp density front, not typically resolved in ocean-circulation models, is central to generating the observed mixing. Interpreting the observations in terms of the flux of freshwater and nutrients across the Gulf Stream north wall suggests an important mechanism of exchange that has not been fully explored.
Abstract
The Gulf Stream front separates the North Atlantic subtropical and subpolar ocean gyres, water masses with distinct physical and biogeochemical properties. Exchange across the front is believed to be necessary to balance the freshwater budget of the subtropical gyre and to support the biological productivity of the region; however, the physical mechanisms responsible have been the subject of long-standing debate. Here, the evolution of a passive dye released within the north wall of the Gulf Stream provides direct observational evidence of enhanced mixing across the Gulf Stream front. Numerical simulations indicate that the observed rapid cross-frontal mixing occurs via shear dispersion, generated by frontal instabilities and episodic vertical mixing. This provides unique direct evidence for the role of submesoscale fronts in generating lateral mixing, a mechanism which has been hypothesized to be of general importance for setting the horizontal structure of the ocean mixed layer. Along the Gulf Stream front in the North Atlantic, these observations further suggest that shear dispersion at sharp fronts may provide a source of freshwater flux large enough to explain much of the freshwater deficit in the subtropical-mode water budget and a flux of nutrients comparable to other mechanisms believed to control primary productivity in the subtropical gyre.
The boundary between the subtropical and subpolar gyres of the North Atlantic is marked by the Gulf Stream front, where lateral gradients in temperature, salinity, density, and nutrients are extremely large (Fig. 1). It is thought that lateral fluxes of nutrients across the front play a critical role in closing nutrient budgets in the subtropical gyre and in providing the phosphate needed to maintain nitrogen fixation there (1⇓⇓⇓–5). Analogously, the flux of cool, low-salinity water across the front has been implicated in the heat and salinity budgets of Eighteen Degree Water—a weakly stratified, voluminous water mass in the subtropical gyre thought to contribute to decadal climate variability (6). However, the physical processes that are responsible for mixing these tracers across the front are not well understood. Fronts by definition are characterized by enhanced tracer gradients—which would tend to be eroded in the presence of significant cross-frontal mixing. Strain associated with large-scale confluent flow provides a mechanism for generating and maintaining the strong gradients; however, these flows are often approximately two-dimensional in the horizontal, suggesting that the fronts will act as dynamical barriers to cross-frontal transport and mixing (7). On the other hand, their strong currents provide a source of energy for turbulent motions such as eddies that can blend the disparate water masses across the front. The debate as to whether the Gulf Stream front in particular acts as a barrier or blender has been ongoing for decades (8, 9).
Overview of the LatMix dye-release experiment along the Gulf Stream front. (A) Tracks of the Lagrangian float (black), R/V Atlantis (green), and R/V Knorr (blue) from yeardays 64 to 68, superimposed on sea-surface temperature from satellite observations on yearday 71 (color). The location of the dye release is indicated, and the R/V Knorr sections shown in Fig. 2 are indicated with dashed white lines. (B) Cross-frontal section of salinity (color) and density (black contours, with values indicated as σ = ρ −1,000 kg
In this article, we describe direct measurements of lateral mixing in the Gulf Stream front, made by using the release and tracking of fluorescein dye. The observations were made as part of the Scalable Lateral Mixing and Coherent Turbulence (LatMix) field campaign, which was devised to study the physics of currents at the submesoscale (0.1 to 10 km)—a range of motions between mesoscale eddies and small-scale, fully three-dimensional turbulence—and their effect on mixing (10, 11). At these scales, the circulation is characterized by fast-growing instabilities, strong nonlinearities, and vigorous overturning circulations in the surface boundary layer. These processes are thought to play an important role in setting the depth of the ocean mixed layer, dissipating the kinetic energy of the ocean circulation, facilitating the growth of phytoplankton, and sequestering carbon in the ocean (12⇓⇓⇓⇓–17). The submesoscale is well below the resolution of most prior work on lateral mixing in the Gulf Stream region—but it had been posited that submesoscale currents could also be crucial to the horizontal dispersal of tracers in the ocean, and one of the aims of the LatMix field campaign was to test this hypothesis. Here, we present evidence from a dye-release experiment in the Gulf Stream, which shows that near the surface, there can be significant exchange across the front, and we use high-resolution numerical simulations to understand the underlying submesoscale processes responsible for the observed dispersion. The essential mechanism we identify—shear dispersion generated by the interaction between submesoscale instabilities and time-varying boundary-layer turbulence—is likely ubiquitous across many regions of the world’s oceans and adds another dimension to the long-standing question of whether the Gulf Stream acts as a barrier or a blender for intergyre exchange in the North Atlantic.
Dye-Release Experiment.
The LatMix observational campaign was conducted along the north wall of the Gulf Stream front (Fig. 1) from 19 February to 17 March 2012 (yeardays 50 to 77). This work involved two global-class research vessels, the R/V Atlantis and the R/V Knorr, and an autonomously profiling Lagrangian (water-following) float (18). A neutrally buoyant mix of approximately 200 gallons of water and 100 kg of fluorescein dye was released at 25-m depth around the Lagrangian float. The evolution of the dye patch was then tracked by following the acoustically tracked position of the Lagrangian float. Measurements of dye concentration, water temperature, salinity, and velocity were taken from a ship-towed profiler in a series of cross-frontal surveys following the float location (Fig. 1) (19).
The dye was released on yearday 64.87 and immediately began spreading across the front (Fig. 2 and Movie S1). The dye center-of-mass remained on a density surface close to the initial injection density for the first 12 h following release, indicating that the initial cross-frontal spreading occurred predominantly due to mixing along isopycnals, which were flattening out during the period of weak surface forcing (Fig. 1). On yearday 65.2, the surface winds intensified and rotated in the clockwise direction, exciting strong currents in the surface mixed layer, and enabling the explosive growth of symmetric instability across the ocean-surface boundary layer (11), a form of small-scale instability active at submesoscale ocean fronts (20). Symmetric instability is associated with an overturning circulation in the across-front direction, which can mix tracers along isopycnals (21) and enhance boundary-layer turbulence by generating secondary Kelvin–Helmholtz instabilities (13, 22). As the wind stress increased and the flow became more unstable to symmetric instability, the dye spread rapidly in the cross-frontal direction, and the center-of-mass moved toward colder and denser water (Figs. 2 and 3). The observations from the Lagrangian float evolved similarly; however, the changes in water-mass properties were more pronounced (Fig. 3). These differences are likely a consequence of the dye reflecting the integrated contributions of several pathways of mixing (SI Appendix, Fig. S1), including vertical mixing across the boundary layer and horizontal mixing between the warm and salty Gulf Stream core and the cold and fresh subpolar gyre (Fig. 2).
Across Gulf Stream sections of observed dye concentration (color) and density (white contours, with values indicated as σ = ρ −1,000 kg
Evolution of the dye from the R/V Knorr observations, showing the spread of the dye in density (Top), temperature (Middle), and salinity (Bottom) space as a function of time. The mean dye concentration of the binned observations is shown (color scale), normalized by the maximum observed concentration across all surveys. The values at the dye center-of-mass (dashed line) and from the Lagrangian float (dashed-dotted line) are also shown, both smoothed by using a cubic-spline over an approximately 8-h window. Gray areas indicate bins where data were not available.
The observed increase in horizontal variance of the dye concentration can be used to estimate the effective horizontal diffusivity of the flow (SI Appendix) (23). The rapid dispersion of the dye in the first 16 h after release suggests an average effective horizontal diffusivity in the across-front direction of
Numerical Simulations.
High-resolution large-eddy simulations (LESs) were used to isolate the physical processes responsible for the observed dye dispersion. The model was run in an idealized, horizontally periodic configuration with an imposed mean background horizontal density gradient that did not evolve in time (described in more detail in SI Appendix). Simulations were initialized on yearday 64.5 with a weakly stratified surface layer overlying a uniformly stratified interior, with parameters chosen to approximate the observations (11) and surface forcing provided by the surface heat fluxes and wind stress calculated from the shipboard observations (Fig. 1C). The most physically realistic “FRONT” run was configured with a mean background-buoyancy gradient consistent with the observations (
Evolution of the passive numerical tracer in the LES run with a horizontal buoyancy gradient (FRONT; Left) and without a horizontal buoyancy gradient (NO-FRONT; Right). Snapshots of along-front averaged dye concentration (colorscale) and buoyancy (white contours, with interval equivalent to a change in density of 0.1 kg
The FRONT run shows a dye evolution consistent with the observations (Fig. 4 and Movie S2). After release, the dye was rapidly mixed vertically throughout the weakly stratified upper 75 m. Vertically sheared cross-frontal flow then tilted over isopycnals and sheared the dye patch, converting horizontal dye gradients into vertical gradients, similar to the observed dye evolution (Fig. 2). Turbulence forced by strong surface winds, and secondary instabilities of symmetric instability, beginning around yearday 65.2 then homogenized the dye concentration in the vertical, causing an irreversible horizontal dispersion of the dye patch. The effective horizontal diffusivity reached peak values of
In contrast to the FRONT run, in the NO-FRONT case, the dye patch was not tilted as strongly in the horizontal and, hence, exhibited very little horizontal dispersion through the wind-mixing event (Fig. 5). The maximum estimated effective horizontal diffusivity was
Across-frontal dye variance,
The sequence of dye evolution in the FRONT run was consistent with the classic Taylor shear-dispersion mechanism (27), whereby advection of a tracer by a sheared velocity field, plus turbulent mixing in the direction of the velocity gradient, generates an effective turbulent dispersion in the direction of the flow (Fig. 6). Here, the generation mechanism for the cross-frontal shear was vertical mixing of the geostrophically balanced along-front flow, which through the thermal wind balance had a uniform vertical shear of magnitude
Schematic of the shear-dispersion mechanism active at a submesoscale front. (A) Initially, the flow is in geostrophic balance, with Coriolis (COR) and across-front pressure gradient forces (PGF) canceling one another. A patch of passive tracer is shown in red. (B) Turbulence, and symmetric instability overturning cells (dashed lines), in the surface boundary layer mix away the thermal wind shear, disrupting the geostrophic balance by reducing the Coriolis force and allowing flow down the pressure gradient. Cross-frontal flow tilts over isopycnals and shears out the tracer patch. (C) Surface heat fluxes (orange arrows) or destablizing winds (gray arrows) generate near-surface turbulence. When this becomes sufficiently strong to overcome the restratifying tendency of the cross-frontal advection of density, the tracer is mixed vertically across isopcynals. In each image, the profile of along-front velocity in the center of the front is shown projected on the right plane.
In simulation NO-GEOMIX, this disruption of the geostrophic balance was artificially suppressed by removing the vertical advection of geostrophic momentum, which stops mixing due to resolved turbulence. As a result, the cross-frontal flow in NO-GEOMIX developed little vertical shear, and, hence, the dye evolved similarly to the NO-FRONT run with minimal tilting of the dye patch through day 65.2 (Fig. 7). The subsequent wind forcing, in turn, generated limited across-isopcynal mixing of the dye, leading to effective diffusivities in NO-GEOMIX that were similar to the NO-FRONT case (Fig. 5). The mixing of geostrophic momentum is therefore central to generating the across-frontal shear flow and the concomitant shear dispersion, evident in the observations and FRONT run. More generally, however, the essential mechanism we have outlined—shear dispersion due to a combination of cross-frontal shear flow and vertical mixing—can arise through a variety of dynamical mechanisms at sharp fronts and is therefore likely a generic feature of regions with active submesoscales. We comment more on this below.
Evolution of the passive tracer in the NO-GEOMIX simulation. Snapshots are shown at three times (indicated to the left of each row) for dye concentration (colorscale) and buoyancy (white contours). Profiles of the horizontally averaged velocity are also shown (blue, across-front velocity; orange, along-front velocity; and dashed gray in top row, background geostrophic velocity). Note that the horizontal spatial coordinate (y), which is periodic, shifts between plots to center the dye patch.
Discussion
The observations presented here provide unique direct wintertime observations of cross-frontal dispersion across the Gulf Stream front, the boundary between the North Atlantic subtropical and subpolar gyres. These water masses have distinct physical and biogeochemical properties (Fig. 1 and SI Appendix, Fig. S4); hence, the observed dye dispersion suggests that submesoscale fronts are acting to enhance fluxes of freshwater and nutrients between the subpolar gyre and the Gulf Stream core. In the North Atlantic, these cross-frontal fluxes have the potential to significantly affect the formation of subtropical-mode water, which provides a pathway for the export of heat and carbon from the surface ocean to the interior, and to affect carbon sequestration through the ocean biological pump (3⇓⇓–6).
As an example, it is informative to evaluate the potential contribution of shear dispersion at submesoscale fronts to the salinity budget of the subtropical mode water formed in the Gulf Stream, which is not currently well understood, despite the implications for decadal climate variability (6). The horizontal diffusive flux of salinity in the mixed layer will scale as
Similar estimates can also be formed for the flux of nutrients such as phosphate and nitrate, which indicates that submesoscale shear dispersion may likewise provide a contribution to exchange across the north wall front that is comparable to prior estimates of lateral advective nutrient fluxes (SI Appendix). Current uncertainties in the spatial, and temporal, distribution of these frontal processes make it difficult to infer the ultimate fate of these tracers, which will have trajectories determined by a combination of lateral mixing processes and advection by the fast-flowing Gulf Stream jet. However, the observed rapid dispersion, and irreversible mixing across isopycnals, suggest that submesoscale shear dispersion should be considered alongside the more well-studied advective mechanisms, such as mesoscale eddy stirring and wind-driven transport (3, 4, 8), for intergyre mixing at the Gulf Stream front.
The physical mechanism identified here as responsible for the observed rapid tracer dispersion, shear dispersion generated by the mixing of geostrophic momentum, will be enhanced at sharp submesoscale fronts where strong buoyancy gradients imply large thermal wind shears. These fronts are also often susceptible to a variety of fast-growing instabilities that contribute to the mixing of geostrophic momentum and the generation of turbulence (20, 22, 32), further enhancing the cross-frontal shear dispersion. However, beyond the particular dynamical progression outlined here, the basic mechanism of shear dispersion due to a combination of sheared cross-frontal flows and time-varying turbulent mixing of the boundary layer is likely a generic feature of submesoscale fronts and is robust to both the details of the turbulent mixing and the origin of the shear flow (28, 29, 33). For example, strong cross-frontal shear flows can be generated at sharp fronts through symmetric instability (as documented here), baroclinic mixed-layer instabilities (34), turbulent mixing (35), or simply lateral slumping due to gravity (36). The observed rapid dye dispersion presented here therefore also provides direct observational confirmation of the hypothesized role of shear dispersion at the submesoscale in reducing tracer gradients in the turbulent surface mixed layer (e.g., compare the salinity gradient above and below 75-m depth in Fig. 1B; and see ref. 37). Along the Gulf Stream, this mechanism has the potential to play an important role in driving exchange across the north wall front, motivating the need for the development of a more complete understanding of submesoscale-turbulence interactions and the cumulative effects of shear dispersion at sharp ocean fronts.
Acknowledgments
We thank the crews of the R/V Knorr and R/V Atlantis for making these observations possible. This work was supported by the Scalable Lateral Mixing and Coherent Turbulence Departmental Research Initiative and the Physical Oceanography Program of the Office of Naval Research, program officers Terri Paluszkiewicz and Scott Harper. Comments and suggestions from Baylor Fox-Kemper and an anonymous reviewer improved this manuscript.
Footnotes
- ↵1To whom correspondence may be addressed. Email: wenegrat{at}umd.edu.
Author contributions: J.O.W., L.N.T., and J.R.T. designed the numerical simulations; J.R.T. performed numerical simulations; M.A.S. was responsible for the collection and analysis of the in situ dye observations; L.N.T., M.A.S., E.A.D., J.M.K., R.K.S., and C.M.L. designed and implemented the LatMix 2012 field campaign; J.O.W. led the analysis of results; J.O.W., L.N.T., M.A.S., J.R.T., E.A.D., J.M.K., R.K.S., and C.M.L. contributed to interpreting results; and J.O.W. wrote the paper with contributions from all authors.
The authors declare no competing interest.
This article is a PNAS Direct Submission.
Data deposition: Data related to this paper are available in Figshare (https://doi.org/10.6084/m9.figshare.12567740). Analysis code is publicly available at GitHub (http://github.com/wenegrat/LATMIXDYE/releases/latest).
This article contains supporting information online at https://www.pnas.org/lookup/suppl/doi:10.1073/pnas.2005558117/-/DCSupplemental.
Published under the PNAS license.
References
- ↵
- R. G. Williams,
- M. J. Follows
- ↵
- A. Oschlies
- ↵
- J. B. Palter,
- M. S. Lozier,
- J. L. Sarmiento,
- R. G. Williams
- ↵
- R. T. Letscher,
- F. Primeau,
- J. K. Moore
- ↵
- A. Yamamoto et al.
- ↵
- T. M. Joyce,
- L. N. Thomas,
- W. K. Dewar,
- J. B. Girton
- ↵
- G. Haller,
- D. Karrasch,
- F. Kogelbauer
- ↵
- ↵
- J. M. Klymak et al.
- ↵
- A. Y. Shcherbina et al.
- ↵
- L. N. Thomas et al.
- ↵
- M. W. Hecht,
- H. Hasumi
- L. N. Thomas,
- A. Tandon,
- A. Mahadevan
- ↵
- E. D’Asaro,
- C. Lee,
- L. Rainville,
- R. Harcourt,
- L. Thomas
- ↵
- M. Lévy,
- R. Ferrari,
- P. J. S. Franks,
- A. P. Martin,
- P. Rivière
- ↵
- M. M. Omand et al.
- ↵
- A. Mahadevan
- ↵
- ↵
- E. A. D’Asaro
- ↵
- J. Wenegrat et al.
- ↵
- ↵
- S. Bachman,
- B. Fox-Kemper,
- J. Taylor,
- L. Thomas
- ↵
- J. R. Taylor,
- R. Ferrari
- ↵
- M. A. Sundermeyer,
- J. R. Ledwell
- ↵
- R. Lumpkin,
- S. Elipot
- ↵
- I. I. Rypina,
- I. Kamenkovich,
- P. Berloff,
- L. J. Pratt
- ↵
- W. G. Large,
- J. C. McWilliams,
- S. C. Doney
- ↵
- G. I. Taylor
- ↵
- W. R. Young,
- P. B. Rhines,
- C. J. R. Garrett
- ↵
- W. R. Young
- ↵
- L. Chen,
- W. R. Young
- ↵
- M. N. Crowe,
- J. R. Taylor
- ↵
- A. F. Thompson et al.
- ↵
- ↵
- ↵
- J. O. Wenegrat,
- M. J. McPhaden
- ↵
- ↵
- D. L. Rudnick,
- R. Ferrari
Citation Manager Formats
Article Classifications
- Physical Sciences
- Earth, Atmospheric, and Planetary Sciences