A warm jet in a cold ocean

Unprecedented quantities of heat are entering the Pacific sector of the Arctic Ocean through Bering Strait, particularly during summer months. Though some heat is lost to the atmosphere during autumn cooling, a significant fraction of the incoming warm, salty water subducts (dives beneath) below a cooler fresher layer of near-surface water, subsequently extending hundreds of kilometers into the Beaufort Gyre. Upward turbulent mixing of these sub-surface pockets of heat is likely accelerating sea ice melt in the region. This Pacific-origin water brings both heat and unique biogeochemical properties, contributing to a changing Arctic ecosystem. However, our ability to understand or forecast the role of this incoming water mass has been hampered by lack of understanding of the physical processes controlling subduction and evolution of this this warm water. Crucially, the processes seen here occur at small horizontal scales not resolved by regional forecast models or climate simulations; new parameterizations must be developed that accurately represent the physics. Here we present novel high resolution observations showing the detailed process of subduction and initial evolution of warm Pacific-origin water in the southern Beaufort Gyre.

Arctic summer sea ice has been rapidly declining1,2,3,4, with significant implications for global climate through albedo and other feedbacks5. In the Western Arctic, sea ice reductions are largest in areas near and just offshore of the Chukchi shelf, where warmer Pacific water enters through Bering Strait6. Several studies point to the rapidly rising heat content and throughput of this in-flowing Pacific origin water as a primary cause for accelerating sea ice melt in this region7,8. Accurate prediction of sea ice melt rates hence requires a better understanding of the processes and pathways by which warm incoming Pacific water enters the main basin.

Unlike most of the world’s oceans, density in the upper Arctic Ocean is primarily controlled by salinity9,10. The Beaufort Gyre (BG) is home to the freshest surface waters, resulting from a combination of river inflow, precipitation, and sea ice melt which are consolidated in the center of the gyre through the convergence of wind-drive flows11. Strong salinity stratification allows heat to be sequestered sub-surface, in water masses that are warmer but saltier than cool, fresh, surface water.

In the Western Arctic, the largest reservoir of heat in the upper ocean (above ~200 m depth) is Pacific Summer Water (PSW)12,13. PSW originates as Pacific ocean water entering through Bering Strait, and is considerably warmed while transiting across the shallow Chukchi shelf during summer months14. An appreciable amount of the summer warmed water subducts, allowing heat to spread into the basin interior, where it is sequestered from direct interaction with the atmosphere. Once subsurface, PSW is found throughout the Beaufort Gyre15. Here we define PSW to be water with potential density between 23.2 and 25.2 kg m−3, warmer than 0 °C, roughly consistent with other recent characterizations8,13,16. Within the central BG, PSW is typically found 30–100 m below the surface, frequently in the form of patchy eddies and filaments13,17,18.

The growing heat content of the PSW sub-surface layer has a first-order effect on accelerating sea ice melt rates in this region, in several ways6,15,19. In the short term, the fate of subducted heat may impact the timing of sea ice growth in late autumn, if it remains close enough to the surface to be rapidly mixed upwards by strong fall storms20. Longer-term and basin-wide, the observed rate of Arctic sea ice decline thermodynamically requires only a 1 W m−2 imbalance in heat exchange between ocean and atmosphere7. Annually averaged, heat flux rates through Barrow Canyon were 3 TW in 201021, a value that has likely grown in the last decade8. If all the heat entering through Barrow Canyon subducted and spread out evenly sub-surface within the BG, that would lead to a warming rate for that layer of 3 W m−215. Observations show that the heat content of the sub-surface PSW within the BG has nearly doubled over the last 30 years15. If all this heat were turbulently mixed upwards, it could melt more than a meter of sea ice15.

Pacific origin water that enters the Arctic through Bering Strait brings not only heat, but supplies nutrient-rich water to the Chukchi and East Siberian shelves22,23. High productivity combined with enhanced rates of sedimentary denitrification across the shelf24,25 results in primary productivity in the western Arctic Ocean being typically nitrate limited26. The physical processes that control the subduction, stirring and mixing rates of Pacific origin water hence also regulate the supply of nutrients that sustain primary production in the BG. An improved understanding of the PSW subduction process is critical for understanding ecosystem functioning and evolution, especially under changing Arctic conditions.

The primary motivation for this work is thus to observe and identify the processes by which PSW subducts. It may seem intuitive that this incoming warm but salty (denser) water would slide under the fresh (lighter) surface water in the BG. However, water of different densities may stably sit side by side in the ocean when the pressure force associated with density gradients is balanced by the Coriolis force acting on currents that flow perpendicular to density gradients, or along fronts27. This type of geostrophic balance is common especially in large-scale oceanic and atmospheric systems. In order for PSW to subduct, something has to break that geostrophic balance. Some PSW may subduct through wind-driven Ekman pumping along the edge of the Chukchi shelf, analogous to ventilation of the thermocline at mid-latitudes13, though more recent work argues this effect is modest28. However the warmest water that subducts to enter the BG travels through Barrow Canyon in late summer as a surface-intensified current16,29,30. Some of the Barrow Canyon outflow has been observed flowing in currents that hug the continental slope towards both the East and West, depending on complex combinations of wind and buoyancy forcing31,32,33. Those currents may in turn be subject to instabilities that facilitate exchanges with the basin interior33,34,35. Yet the process of warm PSW subduction in the Arctic has not been previously observed in detail. Here we present observations of instabilities and warm PSW subduction occurring directly downstream of the Barrow Canyon outflow, and argue that the process of subduction is directly related to meanders in that outflow.

Accurate prediction of the effect of PSW on Arctic sea ice melt and evolving ecosystems requires a better understanding of the physical processes that mediate subduction, stirring and mixing of that water once it enters the Arctic Ocean. The observations presented here show a subduction process with lateral scales of km or smaller. These so-called submesoscale processes cannot be well represented even with high resolution eddy-resolving models (grid scales of 2–3 km, e.g.,36), let alone the much coarser regional or global forecast models37. Though there have been some initial attempts to parameterize the subduction process38, further development is limited by a lack of clear understanding and observational verification of the physical mechanisms involved. The U.S. Office of Naval Research Stratified Ocean Dynamics of the Arctic (SODA) project was designed to explore and understand the evolving upper ocean stratification of the Western Arctic ocean. The novel measurements reported here from SODA scientists and collaborators represent the most detailed observations to date of the subduction process.

Satellite observations from 15 Sept 2018 show a jet of warm water meandering offshore from Barrow Canyon (Fig. 1b). The warmest surface waters in this image are observed directly over Barrow Canyon. The jet initially turns eastward, but then makes a sharp turn northward/offshore near an indentation in the topographic slope (blue contour). Subsequent meanders in the jet of ~50–100 km scale are consistent with previous observations39,40. The Sea Surface Temperature (SST) anomaly associated with the jet is visible with reduced amplitude over 100 km offshore, impinging upon a remnant patch of multi-year sea ice. Convoluted stirring patterns can be seen between the warmest waters of the jet and cooler Arctic surface water. Upper ocean currents measured during the ship survey (Fig. 1c, arrows) are consistent with the meandering pattern visually apparent in SST, with peak amplitudes of over 1 m s−1 in the core of the jet.

Uniquely high-resolution ship-based profiling (see Methods) reveals the complex sub-surface structure of the meandering jet (Fig. 1c). At the onshore edge of the observed jet, water warmer than 4 °C extends from the surface to more than 100 m depth. On the northern/offshore edge (far right of the lowest portion of Fig. 1), the vertical extent of warm water compresses to less than 30 m, centered around 40 m below the surface. The sloping wedge of PSW is overlain by a thin layer of cooler, fresher surface water. On the north-eastern edge of the survey, within the marginal ice zone, the PSW layer is completely sub-surface and appears to have partially broken up into a series of closely spaced sub-surface warm patches.

A cross-sectional slice reveals the internal structure of the warm-water jet in greater detail (Fig. 2). The core of the jet is bounded by the two isopycnals used for our definition of PSW (Fig. 2a, bold white contours). Between the highlighted density contours and in the main jet (cross-jet distance of less than 17 km), the water has an average temperature of 6 ∘C, and a salinity of 30.9 psu. Velocity in the along-jet direction peaks sub-surface and is concentrated between the highlighted density surfaces (Fig. 2c).

The context for this unusually warm water can be seen by comparing its temperature and salinity values to other regional water masses (Fig. 3). As temperature and salinity are both largely conserved along water pathways in the ocean interior, Temperature–Salinity (T–S) plots are often used to classify water masses of different origins. They are also used to identify evidence of mixing between water masses, which manifests on a T-S plot as data on a straight line between end-point water masses41. The water flowing through Barrow Canyon in late summer (Fig. 3a) shows both typical surface intensified warm PSW and cooler slightly saltier Pacific Winter Water below19. The warm water from Fig. 2 can be identified in T–S space as the incoming PSW, which is mixing with both winter water below and cool-fresh melt-water above (Fig. 3b).

To first order this meandering jet appears to be in geostrophic balance; lateral pressure gradients, set by a combination of sea surface slope and lateral density gradients within the water column, show similar patterns to those of the Coriolis force acting on the along-jet flow. Since we do not have a direct measurement of pressure associated with sea-surface height gradients, the balance is assessed by comparing the vertical gradients of both forces, known as thermal wind balance (Methods, Fig. 2e, f). In this case, the curvature of the meandering warm jet visible in Fig. 1 adds a significant term to the effective Coriolis force (Supplementary Methods).

The strength of these currents puts this feature firmly in the nonlinear realm. Much of the instability theory considered by previous work has been based upon small perturbations to modest amplitude currents. The non-linearity of this current system may be quantified by the strength of both the lateral and vertical shear. Lateral shear and curvature in the along-jet velocity create negative (clockwise) relative vorticity on the southern side of the jet, with positive (counterclockwise) relative vorticity on the northern side (Fig. 2g). This relative vorticity is of roughly the same side as the earth’s rotation rate; their ratio gives a non-dimensional Rossby number, which is often used to measure non-linearity. The component of vertical shear associated with this jet is also strong; normalizing by stratification gives a non-dimensional Richardson number of order one (Fig. 2h). The strength of this laterally and vertically sheared current puts it firmly into the category of submesoscale ocean dynamics42,43. Order one values of both Rossby and Richardson numbers are associated with increased propensity to various instabilities that may enhance stirring and mixing of PSW, as discussed further below.

The sub-surface temperature structure visible in Figs. 1c and 2a showcases several steps of the subduction process. Here we loosely interpret the left to right (South to North) variability in these sections as time evolution connecting the more depth-uniform PSW flow that emerges from Barrow Canyon with the discrete sub-surface PSW features visible in the far lower right of Fig. 1c. First, the warm PSW wedge appears to be squeezing sideways (northward in Fig. 2a) and vertically compressing. Second, it is developing a particular pattern of relative vorticity. Within the main wedge of PSW (southward of ~16 km cross-front distance in Fig. 2g), there are alternating broad regions of negative (green) relative vorticity and narrow semi-vertical stripes of positive (red) relative vorticity. Physically, parts of the PSW wedge are starting to spin at different rates relative to each other. Finally, the wedge of PSW is breaking up into discrete sub-surface eddies, as seen in curved, scalloping anomalies of temperature and salinity that line up with the relative vorticity features.

Though sub-surface PSW eddies and filaments have been observed previously, and discussed in theoretical literature, the subduction process has not been observed in detail. Here we discuss each of these observations in light of relevant dynamical processes and describe how they are linked together. We build upon previous work discussing Arctic eddy formation through baroclinic or frontal instabilities20,44,45,46,47,48,49, though again we comment that the strong magnitude of the meandering jet here puts it firmly into a nonlinear, or submesoscale regime, distinct from most previous work in this region.

The first step in the observational evolution of PSW is a squeezing and vertical compression of warm PSW (Fig. 2a). Such behavior is reminiscent of previous theoretical work on meandering jets50. For a large-scale meandering current, the cross-jet flow is expected to be confluent (laterally squeezing together) upstream of each trough in the meander45,51. The section shown in Fig. 2 was taken just upstream of a trough in the large-scale meandering path of the jet (Fig. 1). In this section, the observed cross-jet velocity (Fig. 2d) is indeed confluent, with northward (pink) velocity on the southern side and southward (blue) velocity on the northern side.

Geometrically, confluent cross-jet flow tends to steepen sloping isopycnals, pushing them away from a state of geostrophic balance. Theoretical work shows that the imbalance leads to a series of adjustments that cumulatively tend to restore the flow to geostrophic balance, which is a preferred, stable state for this type of dynamical system. In particular, theory predicts that a modest, cross-front circulation develops which restores the original isopycnal slope52. The anticipated cross-jet circulation can be calculated using the omega equation (see “Methods”), which represents this multi-step process53.

The inferred cross-front circulation associated with that confluence is shown with contours and arrows in 4e, f. Unlike previous studies of individual sloping fronts, here there are two signs of isopycnal slopes, making a double-front. On the upper side of the PSW wedge, isopycnals slope downward and northward, while below they slope upward and northward. The double-front set-up leads to two calculated overturning cells, a deeper clockwise cell (dotted cyan in Fig. 4e, f) and an upper counterclockwise cell (solid cyan). These cells conspire to push warm water further northward (positive cross-front direction) into the wedge (cyan arrows on Fig. 4e) and spread cool fresh surface water southward atop the subducting jet (pink arrows near surface). The circulation vertically compresses the wedge of warm water from both above and below—the top is pulled down while the bottom is pushed up. The predicted vertical velocities are of order 10 m per day and predicted cross-jet velocities are a few km per day.

While the observations are qualitatively consistent with the calculated double-cell cross-front circulation, the calculated circulation may be an under-estimate given the strength of these currents. The theory underlying the cross-jet circulation shown in Fig. 4e, f formally requires weak currents. For the observed strong flows (order one Rossby and Richardson numbers, Fig. 2) the calculated cross-jet circulation is probably a lower bound54. Furthermore, for strong negative relative vorticity (Ro < −1, “Methods”) a fast-growing, energetic and time-dependent process known as centrifugal instability may develop55,56. There are a few locations on the section shown in Fig. 2 and on the section just upstream where Ro is close to −1 (Fig. 2g). There may be even stronger negative relative vorticity upstream, immediately after the flow exited Barrow canyon (Fig. 1). Evidence for this instability is discussed further in the Supplementary material.

Near the surface, different processes may directly spread cool fresh surface water over the top of PSW. Viewed from above, there is a taffy-like stirring pattern between warmer (coastal origin) and cooler (offshore origin) surface water (Fig. 1). The near-surface water visible in Fig. 2 is one slice through this structure. Near-surface water experiences direct wind and wave forcing and higher levels of turbulence. The combination of that forcing differentially advects fresh near-surface water in a way that is partially decoupled from the stratified ocean below, with a net effect of spreading out freshwater filaments (Supplementary).

The second step in the observed process is the development of the striped patterns of relative vorticity, indicating differential rotation rates, within the PSW wedge (Fig. 2g). These developing patterns are consistent with conservation of a quantity known as Ertel potential vorticity (PV)42,51 (Methods). Ertel PV consists of three primary components: a vertical “squashing” term (Fig. 4a), a relative vorticity “spinning” term (Fig. 4b) and a “tilting” term (Fig. 4c). Here we include curvature of the meandering jet flow in the relative vorticity term57,58,59. As PSW is squeezed sideways into the wedge of sloping isopycnals, it is vertically compressed, increasing the squashing component of PV (Fig. 2a). This leads to enhancement of the spinning component of PV within the PSW wedge (Fig. 2b), in order to conserve total PV (Fig. 2d). Interplay between the different spatially varying components of this conserved quantity (Fig. 4) likely produces the observed spatial patterns in relative vorticity. A further step in the water mass evolution may be seen in the furthest offshore observations, on the far lower right of Fig. 1; here PSW is visible in more distinct features, which may be the precursors of individual intrathermocline eddies18. The horizontal length scales of the developing relative vorticity patterns, and eventual intrathermocline eddies, are thought to be set by geometric and dynamic constraints as those features adjust back towards a state of geostrophic balance48,50.

As warm PSW subducts and vertically compresses, the developing patterns of relative vorticity lead to differentiation and discretization of water mass features in Fig. 2. For example, the noses of two scalloping temperature and salinity anomalies can be seen around 50 m depth and 5 and 7 km in the cross-jet distance coordinate in Fig. 2a, b. They extend both downward and southward, and in more muted form upward and southward. These wispy features involve water that is cooler but fresher than the PSW to either side. A more pronounced scalloping feature is centered around 50 m depth and 10.5 km cross-jet distance (Fig. 2a, b), similarly with water that is cooler but fresher than that on either side. In all cases, the temperature and salinity anomalies have compensating effects on density, as evidenced by the small deviations of the contoured isopycnal lines in both panels. Such temperature and salinity anomalies cannot result from vertical motions, as the water below the jet is cooler but saltier. Instead, we posit that the water mass anomalies result from stirring along isopycnals at the edges of this meandering jet, in and out of the page from this section perspective. These water mass anomalies line up with the relative vorticity anomalies (Fig. 2g), providing strong evidence that the developing relative vorticity features described above are responsible for the lateral stirring.

The relationship between lateral stirring and the temperature or salinity anomalies in Fig. 2 is further elucidated by turning back to the T-S diagrams (Fig. 3). When PSW enters the BG (Fig. 3, panel b) in the form of this warm jet, it meanders through water characterized by typical regional temperature and salinity vertical profiles (panel c). Those background profiles show a straight mixing line between surface melt-water above (left side of Fig. 3c), and Pacific winter water below (right side of panel c). The developing relative vorticity features described above stir and mix PSW with the background water along isopycnals, between the peak of the triangle shown in panel b, and the lower leg of that triangle shown in panel c, especially in the isopycnal range bounded by the two solid black lines identifying PSW density. That stirring fills in the T-S triangle, as can be seen in a section further downstream (Fig. 3d). The scalloping features visible in Fig. 2a, b are the initial stage of this stirring process.

During a subsequent survey of the same region 9 days later, the subduction process appears to be complete, with cool fresh water found everywhere near the surface (Fig. 5a, b). Compared to the original survey, mixing and stirring along and perhaps across isopycnals has now filled in the entire triangle in T-S space (Fig. 3e). The stirring process can be seen in more detail by considering four representative profiles through the second survey (Fig. 5 colors). There is some indication that subducted PSW (green profile) is mixing with background water along the same isopycnal found just outside the eddies (e.g., black profile). Similarly, it appears as though the red and blue profiles have exchanged water properties near the temperature maximum (Fig. 5d). The T-S profiles and the bite taken out of the eddy in the red profile section of the snapshot (Fig. 5a) suggest an intriguing form of eddy cannibalism induced by complex three-dimensional velocity and vorticity fields.

Here we consider both the evolution of oceanic heat within this particular event and the implications for the basin as a whole. Combining observed temperature and velocity in the primary survey gives a lateral heat transport of 35-40 TW for the sections shown in Fig. 1 (“Methods”). This heat flux is probably an underestimate of the true flux, as our measurements do not include the southern edge of the jet or the more complex eddying meanders visible in the satellite image. These heat flux estimates are consistent with but near the high end of previous measurements of peak late summer heat fluxes through Barrow Canyon, consistent with warming PSW inflow16,21.

For the subduction event described here, some incoming heat is lost to the atmosphere. Over the course of the initial survey (14–17 Sept) average ship-measured heat flux from the ocean to the atmosphere was 122 W m−2, similar to previous measurements60 (Supplementary). As a simple scaling, a heat loss of 100 W m−2 cools a 10m surface layer ~0.2 °C/day, which is consistent with what was observed by floats transiting Barrow Canyon over ~10 days. It cannot explain the evolution of the jet in TS space, which must come principally from isopycnal stirring and mixing with background waters of the Beaufort Gyre (Fig. 3). During the early stages of the subduction process, when PSW is near the surface, the measured upward turbulent heat fluxes within the ocean are large, approaching the rate of surface heat loss to the atmosphere. However, once the PSW core has been sequestered further below the surface, measured turbulent heat fluxes are less than 5 W m−2 (Supplementary). By the second survey on 24 Sept, air–sea heat fluxes were only 16 W m−2 from ocean to atmosphere. The net picture is one of warm water steadily being pulled sub-surface, losing a small portion of its heat upwards though mixing and then air-sea heat loss. As a rough upper bound, integrating a heat flux to the atmosphere of 100 W m−2 over the roughly 100 km × 100 km area shown in Fig. 1c gives a net loss of 1 TW, a small portion of the 40+ TW entering the Arctic in this jet. We conclude that the majority of the heat in this event subducts, where it is subsequently sequestered from direct contact with the atmosphere or sea ice. Given that typical maximum winter mixed-layer depths in the BG are only 30-35m61, most of this subducted heat may be insulated for months if not longer.

As subducted PSW spreads into the main gyre, the net effect on sea ice will depend on both the lateral stirring rate and the upward diffusion rate of heat. Some water is entrained into eastward or westward boundary currents. For example, sub-surface PSW with the same T-S characteristics as observed here was observed flowing westward in the Chukchi Slope Current62. Longer term, both PSW subducted directly from the Barrow Canyon outflow, and the portion subsequently caught up in boundary currents, may be drawn into the central Beaufort Gyre and beneath the main ice pack through a combination of advection by the anti-cyclonic circulation and eddy stirring63. Water with similar T-S properties was observed passing by the SODA-A mooring, 120 km to the north-east, in early November 2018 (Fig. 3f). Warm water (>3 °C) was observed at the mooring in a series of pulses from 28 October–13 November, between 45 and 90 m; no measurements are available higher in the water column. The maximum temperature observed at SODA-A is slightly reduced compared to the ship survey(Fig. 3f), perhaps by continued stirring and mixing. By November this region was fully ice covered; these measurements thus showcase persistence of very warm water lurking beneath a frozen surface.

In addition to carrying sub-surface heat into the basin, the warm PSW carries important biogeochemical tracers that reflect the meandering jet’s Pacific origins and transit across the Chukchi shelf in the summer months prior to the field program. Within the warm jet, biogeochemical properties display complex patterns that mirror the intricate subduction, stirring and mixing that controls heat distribution. The profiling, drifting Wirewalker data shows clear signatures in both Chlorophyll concentration and CDOM (Chromophoric Dissolved Organic Matter, Fig. 6b, c) that follow the warmest temperatures (Fig. 6a). CDOM is elevated in a range that maps well onto the 6° C isotherm (contoured cyan), reflecting both the substantial river influence for the warmest incoming PSW64, and potentially also scouring as water travels along the Chukchi shelf and through Barrow Canyon.

Surface waters and the upper layer of the plume are nitrate-depleted (<0.1 mmol m−3), reflective of algal uptake across the shallow shelf, though some nitrate may be restored through turbulent fluxes within Barrow Canyon30. As mentioned above, the stirring of near-surface waters responds to additional wind and wave forcing, so the surface layer may deferentially advect compared to the sub-surface plume. Nitrate concentrations increase within the bottom 20–30 m of the plume and reach ~5 mmol m−3 at its base (the 25.2 kg m−3 isopycnal) (Fig. 7a). The top of the nitracline is thus embedded within the plume.

Chlorophyll-a concentrations are elevated throughout the upper layers of the plume and in the sunlit surface waters, but especially at a depth that maps well onto the top surface of the PSW jet (Figs. 6b, 7b). The depth of the top of the plume and the chlorophyll maximum, deepen in time as the jet flows northward and subducts downward, carrying with it and redistributing shelf-origin waters. Within the upper, nitrate-depleted layers biomass, specifically particulate carbon and particulate organic nitrogen (PC and PON, respectively), is more than twofold higher in the jet (62 ± 29 μg L−1 PC and 12 ± 6 μg L−1 PON, respectively) compared to the surrounding Beaufort Gyre (26 ± 3 μg L−1 PC and 5 ± 0.5 μg L−1 PON, Fig. 7d, e and Supplementary). Evidence of enhanced phytoplankton growth within the jet is also found in the stable carbon and nitrogen isotope signatures (Fig. 7f, g and Supplementary). Carbon isotopes of PC are ~2‰ lower in the jet (−26.2 ± 0.8‰) compared to the surrounding Beaufort Gyre (−28.8 ± 0.1‰), indicating enhanced primary productivity65. Further, nitrogen isotopes of PON are enriched in the jet (6.7‰) compared to the surrounding Beaufort Gyre (4.6 ± 3.2‰), indicating complete assimilation of isotopically-enriched nitrate within the jet, reflective of the Bering–Strait–Chukchi pathway24,25 and incomplete nitrate assimilation in the surrounding waters.

Algal growth and biomass accumulation beyond that which has been exported off the shelf within the plume is likely to be severely nitrate limited. However, lateral stirring induced by the jet’s developing relative vorticity (Fig. 2g) leading to vertical shear and enhanced turbulent kinetic energy dissipation18 may drive fluxes of new nitrate across the nitracline and alleviate nitrate limitation within the euphotic zone. Algal growth within subsurface chlorophyll layers makes an important, yet still poorly constrained contribution to Arctic primary production66,67. Thus the capacity of this shallow warm meandering jet to carry carbon-laden waters from the shelf into the carbon-poor basin68 whilst also providing a mechanism for vertical nutrient fluxes may play an important role in biological production across the Canadian Basin. The net impact is complementary to the nutrient impact carried into the BG interior by deeper, cold-core eddies69.

The novel observations presented here reveal a jet of warm salty water that appears to be (1) subducting beneath a cool fresh surface layer, and (2) vertically compressing and breaking up into smaller, spinning eddies. Though some heat is lost to the atmosphere, the vast majority subducts to be sequestered from direct contact with the atmosphere, and subsequently stirred and spread sub-surface toward the central basin. Both the subduction and development of relative vorticity are qualitatively consistent with cross-jet circulations associated with mesoscale jet meanders, submesoscale instabilities, and potential vorticity conservation.

Meandering warm jets offshore of Barrow Canyon are not unusual33,39,40,60. In previous studies, some of that heat leaving the coast near Barrow Canyon seems to be firmly ejected offshore, while some meanders to rejoin one of the boundary currents, and may subsequently leak into the interior by other stirring processes62. The subduction and intrathermocline eddy development processes observed here may be a regular occurrence, linking the larger-scale influxes of warm, surface-intensified PSW observed within Barrow Canyon21,29 with discrete intrathermocline eddies and filaments observed within the basin17,18.

The net effect on redistributing heat and biogeochemical properties in the main basin will depend on the interplay between lateral stirring, vertical mixing, and the evolving large-scale basin circulation. How long warm intrathermocline eddies last, and thus how far they get under the sea ice, is an open question involving their translation speed, intrinsic stability59, and frictional effects of the ice above36,70. While lateral transport of shelf-origin waters helps maintain the upper Arctic halocline and its nutrient maxima71, upward nutrient fluxes through the halocline depend both on halocline ventilation rates and lateral fluxes within the surface layer; differential changes may alter both biological production and biogeochemical cycling rates72.

As the heat content of PSW is growing8,14,21, the combination of PSW subduction, lateral stirring, and upward vertical mixing should lead to a pattern of accelerating sea ice melt spreading out from the Pacific inflow, as has been observed in recent decades6,73. Related processes have been observed on the other side of the Arctic, associated with warming of subducted Atlantic water74. At the same time, the BG is spinning up, with a growing accumulation of near-surface freshwater that tends to deepen the halocline and nitricline75,76. However, as the PSW warms and lightens, it is subducting along shallower isopycnals, putting it closer to the surface77. The physical insights gained with this work should help develop more accurate ways to represent the subduction and stirring processes in models. Improved understanding and modeling ability of the processes described here will help forecast the detailed geography and time-frame for the changing ecosystem and accelerating Arctic sea ice loss.

Most of the measurements presented here were conducted during September 2018 aboard the R/V Sikuliaq, operated by the University of Alaska. Details of each measurement technique, as well as the analysis calculations presented above, are shown here. Additional views of the intermediates steps of these calculations are presented in Supplementary material.

The hybrid MODIS image in Fig. 1b was created by combining a Level 2 SST product created by the Ocean Biology Processing Group, NASA Goddard Space Flight Center, and a true-color image created at the University of Miami by applying the NASA Blue Marble algorithm to MODIS bands 1, 4, and 3. The two products were interpolated to a common spatial grid, and a threshold was applied to the red band of the interpolated true-color image. Pixels above this threshold, corresponding to sea ice, land, and clouds, were assigned their true color in the hybrid image. Color for the below-threshold pixels, corresponding to open ocean surface, was assigned using the temperature scale shared with the sub-surface measurements in Fig. 1c.

The Fast CTD (conductivity-temperature-depth) system is a rapid profiler built and operated by the Multiscale Ocean Dynamics Group at the Scripps Institution of Oceanography, UC San Diego. The system was operated from a custom direct drive electric winch mounted on the starboard quarter of the ship, and profiled off a 10 m long boom that is deployed outboard and aft, in order to avoid contamination by the ship wake. With the very low drag FCTD fish, vertical profiling speeds up to 5 m/s were conducted while the ship was steaming at up to 5 knots. Temperature, conductivity, and pressure measurements were made with a Seabird SBE49 instrument, sampled at 16 Hz, and binned to 0.5 m vertical resolution. Temperature and conductivity data were adjusted to match phase and resolution before computing salinity, using standard techniques. Profiles to 200 m depth were completed every few minutes, leading to an average horizontal spacing of 160 m between profiles. Salinity and water density were calculated from temperature and conductivity using standard routines.

A Wirewalker drifting, wave-powered profiling system78 conducted one deployment during this survey (Fig. 1c, red track). The vehicle is equipped with CTD (RBR Concerto), velocity (Nortek Signature 1 MHz), optics (WETLabs chlorophyll and CDOM fluorescence, 532 nm backscatter), and a microstructure sensor. The surface buoy is outfitted with a satellite tracker. The profiler was deployed on a 100-m wire. Profiling speed depended on the wave state, with the average profiling frequency 7.6 cycles per hour. The pathway of Wirewalker drift is shown in Fig. 1c. The drift lasted from 07:30 on 15 Sept 2018 to 06:30 on 16 Sept 2018. Wirewalker chlorophyll data have been corrected for non-photochemical quenching79. Although water bottle samples for CDOM were not collected, the vertical gradients shown in Fig. 6 are an order of magnitude larger than the instrument sensitivity (~0.2 ppb) and therefore reflect actual environmental variability.

Ocean currents were measured with a 300 kHz Acoustic Doppler Current Profiler from Teledyne RDI, installed in the centerboard drop keel of the Sikuliaq. Data were sampled with 1 Hz and 2 m vertical resolution. Data presented in Figs. 1 and 2 has been smoothed and sub-sampled to 2 min temporal resolution. Velocity at this resolution has an uncertainty of ±2 cm/s due to instruemnt noise. The shear data shown in Fig. 2 have been further smoothed to 4 min and 4 m resolution, producing an estimated shear uncertainty of 0.003 s−1.

Temperature and salinity profile data from September 2018 used in Fig. 3a. come from a hydrographic database that was compiled and quality-controlled as described by8 with additional data from an array of autonomous ocean profilers deployed jointly by the Woods Hole Oceanographic Institution and the NOAA Arctic Heat experiment80.

Turbulent fluxes over open water were computed from ship meteorological sensors using the COARE 3.5 algorithm81 as refined by82, using the cool-skin option. Additional details are in the Supplementary material.

Upward turbulent heat fluxes are calculated from the Modular Microstructure Profiler (MMP, track shown in Fig. 1c, blue line). MMP is a loosely tethered free-fall instrument ballasted to at 0.7 m s−1. As the instrument falls, air-foil shear probes measure shear with centimeter resolution83,84. The MMP was profiled from the FastCTD winch/boom system, repeatedly cycled to a depth of 100 m. Turbulent dissipation rates (ϵ) were computed from the shear data using a frozen field hypothesis and standard spectral fitting techniques85, with a noise level near 10−10 W kg−1. The MMP is also equipped with a pumped Seabird CTD for measuring temperature, conductivity, and pressure, from which salinity and potential density are calculated using standard techniques. Turbulent heat flux (Jq) was calculated as \({J}_{q}=-\rho {c}_{p}\frac{0.2\epsilon }{{N}^{2}}\frac{dT}{dz},\) where ρ is the ocean density, cp the specific heat of seawater, N2 the buoyancy frequency, and T(z) the measured vertical profile of temperature. See Supplementary for further details.

SWIFT drifters provide estimates of the dissipation rate of turbulent kinetic energy in the upper few meters of the surface ocean86. Profiles of turbulent velocity fluctuations are measured using the center beam of a down-looking Nortek Signature1000 in the hull of the SWIFT, using a pulse-coherent “HR” mode with 4 cm bin size and 8 Hz sampling. These profiles are used to compute a spatial structure-function, from which the dissipation rate is estimated by fitting an r2/3 dependence (r is the distance between bins). These estimates are produced every 12 min, in a surface wave-following reference frame. Conductivity and temperature are measured in situ at the upper most bin (0.3 m below the surface).

For isotopic analysis, seawater was collected using Niskin bottles attached to a CTD rosette frame. Between 1.8 and 3.7 L of seawater was filtered onto pre-combusted glass fiber filters (nominal pore size of 0.7 micron) using vacuum filtration. Filters were stored frozen, then freeze-dried for 12 h and pelletized. Filters were analyzed at the University of Liverpool using a Costech Elemental Analyzer coupled to a Thermo isotope ratio mass spectrometer. Particulate carbon and particulate organic nitrogen concentrations and stable carbon and nitrogen isotope values were corrected using international reference standards (USGS40 and USGS41A) and expressed in δ notation (δ13C-PC (‰ vs VPDB) = (Rsam/Rstd − 1) × 1000 and δ15N-PN (‰ vs AIR) = (Rsam/Rstd − 1) × 1000). Standards were run in triplicate with a reproducibility better than 0.1‰ for carbon and 0.05 ‰ for nitrogen.

Samples for the analysis of dissolved inorganic nutrients (nitrate, nitrite, ammonium, phosphate and silicate) were collected from multiple depths resolving full-depth profiles. The location of the casts is indicated in white dots in Fig. 3g, and presented in Table 1 of the Supplementary material. 100 mL of seawater were collected directly from Niskin bottles into 120 mL HDPE bottles (10% HCl pre-cleaned) using an inline AcroPak (0.45 μm pore size) filter. Samples were then stored frozen (upright) at −20 °C for later analysis on land. Analyses were carried out at the Marine Chemistry Laboratory, University of Washington, using an AA3 Seal Analytical continuous segmented flow nutrient analyzer following standard colorimetric techniques. Overall precision and accuracy were equal or better than 2%, as assessed via measurements of OSIL Scientific certified materials. For the purpose of the current study, only nitrate data are shown.

Lateral heat flux carried by the warm jet was calculated as ρcpU(T − Tf) where ρ is the ocean density, cp the specific heat, U the along-jet velocity, and Tf is the salinity dependent freezing point of seawater8,14,21.

In Fig. 2, the vertical shear (dU/dz) is calculated by first-differencing along-jet velocity (U), and then smoothing with box-car filters of 425 m horizontally and 4 m vertically. Here the shear has been scaled by the quantity \(\left.{f}_{c}=f-2* (U/r)\right)\), i.e., the local inertial frequency f adjusted by twice the curvature vorticity of the jet. This term accounts for the form of thermal wind balance expected for cyclogeostrophic flows57,58,59. We use the best fit for curvature at this point, though the sign and magnitude of the curvature varies throughout the meandering process.

The lateral buoyancy gradient is calculated by first-differencing buoyancy (\(b=-(g/{\rho }_{0})\rho ^{\prime}\) where ρ0 = 1025 kg m−3 is a reference density, and \(\rho ^{\prime}\) is the potential density calculated from temperature and salinity) in the across-jet direction (y), smoothing with the same box-car filters. Lateral gradients in along-jet velocity (dU/dy) are calculated similarly to those of buoyancy. The Rossby number (Ro = [ −dU/dy − U/r]/f) takes into account the curvature of the meandering jet. To estimate r we use the best fit for the radius of curvature of the flow, r = 30 km, at this point.

The Ertel Potential Vorticity (PV, Fig. 4) is a conserved quantity of an inviscid, adiabatic fluid in a rotating reference frame42. Following previous work for currents with a appreciable curvature57,58 we write the three terms comprising the PV as \(PV=f\frac{\partial b}{\partial z}-\left(\frac{\partial U}{\partial y}+\frac{U}{r}\right)\frac{\partial b}{\partial z}+\frac{\partial U}{\partial z}\frac{\partial b}{\partial y}\), which from left to right are the quantities labeled squashing, spinning, and tilting in Fig. 4a–c.

The ageostrophic cross-jet velocity in Fig. 2g is calculated by numerically solving the adiabatic, two-dimensional version of the quasigeostrophic omega equation \({f}^{2}\frac{{\partial }^{2}\psi }{\partial {z}^{2}}+\frac{\partial b}{\partial z}\frac{{\partial }^{2}\psi }{\partial {y}^{2}}=2\frac{\partial V}{\partial y}\frac{\partial b}{\partial y}\) where Ψ is the streamfunction, V the cross-jet velocity. The omega equation was solved over the domain shown in Fig. 2 using boundary conditions of \(\frac{\partial \psi }{\partial x}=0\) at the top and bottom boundaries and \(\frac{\partial \psi }{\partial z}=0\) at the side boundaries, using the Matlab minimal residual method solver with a tolerance of 10−4, following53. From the streamfunction, ageostrophic cross-front and vertical velocities are calculated as \({v}_{a}=\frac{\partial \psi }{\partial z}\), \({w}_{a}=-\frac{\partial \psi }{\partial y}\) respectively. Cross-jet velocity confluence (dV/dy) was extrapolated down to 130 m by assuming it had the same vertical decay as lateral buoyancy gradients, both of which characterize the mesoscale flow. The lower boundary of 130 m is somewhat arbitrary. Imposing a lower boundary at the base of the measured velocity range (~90 m) creates an artificial and unrealistic lower boundary to the streamfunction; density measurements show the sloping front extends deeper. However extrapolating velocity data below the range of measurements is also artificial. Sensitivity tests were performed using lower boundary conditions of 100–200 m, in 20 m increments. For deeper lower boundaries, the magnitude of the cross-jet circulation increases, but all simulations show the same qualitative pattern as in Fig. 4. Given that the sensitivity in magnitude of cross-front circulation strength is likely smaller than the intrinsic caveats associated with using this method for a Ro ~ 1 flow, we emphasize the qualitative pattern.

• Arctic Heat data are available through https://www.pmel.noaa.gov/arctic-heat/data

• All ship-board and mooring data collected as part of the US Office of Naval Research Stratified Ocean Dynamics of the Arctic program is still in the process of being organized and archived. All data will be fully released to the public when the program formally concludes. In the meantime, details about the program and data collected can be seen here (http://www.apl.washington.edu/project/project.php?id=soda). SODA data presented in this study are available from the corresponding author upon reasonable request.

• Data from the ARISE, PEANUTS, and BMBF projects are also available from the corresponding author upon reasonable request.

Perovich, D. K. et al. Arctic sea-ice melt in 2008 and the role of solar heating. Ann. Glaciol. 52, 355–359 (2011).

ADS  Article  Google Scholar 

Serreze, M. C., Crawford, A. D., Stroeve, J. C., Barrett, A. P. & Woodgate, R. A. Variability, trends, and predictability of seasonal sea ice retreat and advance in the Chukchi Sea. J. Geophys. Res. Oceans 121, 7308–7325 (2016).

ADS  Article  Google Scholar 

Onarheim, I. H., Eldevik, T., Smedsrud, L. H. & Stroeve, J. C. Seasonal and regional manifestation of Arctic sea ice loss. J. Clim. 31, 4917–4932 (2018).

ADS  Article  Google Scholar 

Perovich, D. et al. The Arctic sea ice cover in "State of the Climate in 2018”. Bull. Am. Meteorol. Soc. 100, S146–S150 (2019).

Google Scholar 

Pistone, K., Eisenman, I. & Ramanathan, V. Observational determination of albedo decrease caused by vanishing Arctic sea ice. Proc. Natl Acad. Sci. 111, 3322–3326 (2014).

ADS  CAS  PubMed  Article  Google Scholar 

Shimada, K. et al. Pacific Ocean inflow: Influence on catastrophic reduction of sea ice cover in the Arctic Ocean. Geophy. Res. Lett. 33, L08605 https://doi.org/10.1029/2005GL025624 (2006).

Carmack, E. et al. Towards quantifying the increasing role of oceanic heat in sea ice loss in the new Arctic. Bull. Am. Meteorol. Soc. 96, 2079–2105 (2015).

ADS  Article  Google Scholar 

Danielson, S. et al. Manifestation and consequences of warming and altered heat fluxes over the Bering and Chukchi Sea continental shelves. Deep-Sea Res. II 177 (2020).

Carmack, E. C. The alpha/beta ocean distinction: a perspective on freshwater fluxes, convection, nutrients and productivity in high-latitude seas. Deep Sea Res. II 54, 2578–2598 (2007).

ADS  Article  Google Scholar 

Jackson, L., Hallberg, R. & Legg, S. A parameterization of shear-driven turbulence for ocean climate models. J. Phys. Oceanogr. 38, 1033–1053 (2008).

ADS  Article  Google Scholar 

Carmack, E. C. et al. Freshwater and its role in the Arctic marine system: Sources, disposition, storage, export, and physical and biogeochemical consequences in the Arctic and global oceans. J. Geophys. Res.: Biogeosci. 121, 675–717 (2016).

CAS  Article  Google Scholar 

Steele, M. et al. Circulation of summer Pacific halocline water in the Arctic Ocean. J. Geophys. Res. 109, https://doi.org/10.1029/2003JC002009 (2004).

Timmermans, M.-L. et al. Mechanisms of Pacific Summer Water variability in the Arctic’s Central Canada Basin. J. Geophys. Res.: Oceans 119, 7523–7548 (2014).

ADS  Article  Google Scholar 

Woodgate, R. A. Increases in the Pacific inflow to the Arctic from 1990 to 2015, and insights into seasonal trends and driving mechanisms from year-round Bering Strait mooring data. Prog. Oceanogr. 160, 124–154 (2018).

ADS  Article  Google Scholar 

Timmermans, M.-L., Toole, J. & Krishfield, R. Warming of the interior Arctic Ocean linked to sea ice losses at the basin margins. Sci. Adv. 4, eaat6773 https://doi.org/10.1126/sciadv.aat6773 (2018).

Pickart, R. S. et al. Seasonal to mesoscale variability of water masses and atmospheric conditions in Barrow Canyon, Chukchi Sea. Deep Sea Res. II 162, 32–49 (2019).

Article  Google Scholar 

Kawaguchi, Y., Itoh, M. & Nishino, S. Detailed survey of a large baroclinic eddy with extremely high temperatures in the Western Canada Basin. Deep Sea Res. I 66, 90–102 (2012).

Article  Google Scholar 

Fine, E. C., MacKinnon, J. A., Alford, M. H. & Mickett, J. B. Microstructure observations of turbulent heat fluxes in a warm-core Canada Basin eddy. J. Phys. Oceanogr. 48, 2397–2418 (2018).

ADS  Article  Google Scholar 

Timmermans, M.-L. & Marshall, J. Understanding Arctic Ocean circulation: a review of ocean dynamics in a changing climate. J. Geophys. Res.: Oceans 125, e2018JC014378 (2020).

ADS  Article  Google Scholar 

Lu, K. et al. Lateral mixing across ice meltwater fronts of the Chukchi Sea shelf. Geophys. Res. Lett. 42, 6754–6761 (2015).

ADS  Article  Google Scholar 

Itoh, M. et al. Water properties, heat and volume fluxes of Pacific water in Barrow Canyon during summer 2010. Deep-Sea Res. Part I 102, 43–54 (2015).

Article  Google Scholar 

Torres-Valdés, S. et al. Export of nutrients from the Arctic Ocean. J. Geophys. Res.: Oceans 118, 1625–1644 (2013).

ADS  Article  Google Scholar 

Tremblay, J.-É. et al. Global and regional drivers of nutrient supply, primary production and CO2 drawdown in the changing Arctic Ocean. Prog. Oceanogr. 139, 171–196 (2015).

ADS  Article  Google Scholar 

Fripiat, F. et al. Influence of the bordering shelves on nutrient distribution in the Arctic halocline inferred from water column nitrate isotopes. Limnol. Oceanogr. 63, 2154–2170 (2018).

ADS  CAS  Article  Google Scholar 

Granger, J., Sigman, D. M., Gagnon, J., Tremblay, J.-E. & Mucci, A. On the properties of the Arctic halocline and deep water masses of the Canada Basin from nitrate isotope ratios. J. Geophys. Res.: Oceans 123, 5443–5458 (2018).

ADS  Article  Google Scholar 

Mills, M. M. et al. Nitrogen limitation of the summer phytoplankton and heterotrophic prokaryote communities in the Chukchi Sea. Front. Mar. Sci. 5, 362 (2018).

Article  Google Scholar 

Vallis, G. K. Atmospheric and oceanic fluid dynamics (Cambridge University Press, 2017).

Meneghello, G., Marshall, J., Timmermans, M.-L. & Scott, J. Observations of seasonal upwelling and downwelling in the Beaufort Sea mediated by sea ice. J. Phys. Oceanogr. 48, 795–805 (2018).

ADS  Article  Google Scholar 

Brugler, E. T. et al. Seasonal to interannual variability of the Pacific water boundary current in the Beaufort Sea. Prog. Oceanogr. 127, 1–20 (2014).

ADS  Article  Google Scholar 

Beaird, N., Shroyer, E., Juranek, L., Hales, B. & Goñi, M. Nutrient-rich gravity current formed by upwelling in Barrow Canyon: high resolution observations. J.Geophys. Res.: Oceans 125, e2020JC016160 (2020).

Article  Google Scholar 

Nikolopoulos, A. et al. The western arctic boundary current at 152 w: Structure, variability, and transport. Deep Sea Res. Part II: Top. Stud. Oceanogr. 56, 1164–1181 (2009).

ADS  Article  Google Scholar 

Corlett, W. B. & Pickart, R. S. The Chukchi slope current. Prog. Oceanogr. 153, 50–65 (2017).

ADS  Article  Google Scholar 

Spall, M. A. et al. Transport of Pacific water into the Canada Basin and the formation of the Chukchi slope current. J. Geophys. Res.: Oceans 123, 7453–7471 (2018).

ADS  Article  Google Scholar 

Pickart, R. S. Shelfbreak circulation in the alaskan beaufort sea: Mean structure and variability. J. Geophys. Res.: Oceans (1978–2012) 109, C04024 https://doi.org/10.1029/2003JC001912 (2004).

von Appen, W.-J. & Pickart, R. S. Two configurations of the western Arctic shelfbreak current in summer. J. Phys. Oceanogr. 42, 329–351 (2012).

ADS  Article  Google Scholar 

Meneghello, G. et al. Genesis and decay of mesoscale baroclinic eddies in the seasonally ice-covered interior Arctic Ocean. J. Phys. Oceanogr. 51, 115–129 (2021).

ADS  Article  Google Scholar 

Notz, D. & Community, S. Arctic sea ice in cmip6. Geophys. Res. Lett. 47, e2019GL086749 (2020).

Article  Google Scholar 

Lu, K., Danielson, S., Hedstrom, K. & Weingartner, T. Assessing the role of oceanic heat fluxes on ice ablation of the central Chukchi Sea Shelf. Prog. Oceanogr. 184, 102313 (2020).

Watanabe, E. Beaufort shelf break eddies and shelf-basin exchange of Pacific summer water in the western Arctic Ocean detected by satellite and modeling analyses. J. Geophys. Res. 116, C08034 (2011).

ADS  Google Scholar 

Kozlov, I. E., Artamonova, A. V., Manucharyan, G. E. & Kubryakov, A. A. Eddies in the Western Arctic Ocean from spaceborne SAR observations over open ocean and marginal ice zones. J. Geophys. Res.: Oceans 124, 6601–6616 (2019).

ADS  Article  Google Scholar 

Talley, L. D. Descriptive physical oceanography: an introduction (Academic press, 2011).

Thomas, L. N., Tandon, A. & Mahadevan, A. Submesoscale processes and dynamics. Ocean Model. Eddying Regime 177, 17–38 (2008).

Article  Google Scholar 

McWilliams, J. C. Submesoscale currents in the ocean. Proc. R. Soc. A 472, 20160117 (2016).

ADS  PubMed  MATH  Article  Google Scholar 

Manley, T. & Hunkins, K. Mesoscale eddies of the Arctic Ocean. J. Geophys. Res. 90, 4911–4930 (1985).

ADS  Article  Google Scholar 

Spall, M. A., Pickart, R. S., Fratantoni, P. S. & Plueddemann, A. J. Western Arctic shelfbreak eddies: formation and transport. J. Phys. Oceanogr. 38, 1644–1668 (2008).

ADS  Article  Google Scholar 

Pickart, R. S. & Stossmeister, G. Outflow of Pacific water from the Chukchi Sea to the Arctic Ocean. Chin. J. Polar Sci. 19, 135–148 (2008).

Google Scholar 

Timmermans, M.-L., Toole, J., Proshutinsky, A., Krishfield, R. & Plueddemann, A. Eddies in the Canada Basin, Arctic Ocean, observed from ice-tethered profilers. J. Phys. Oceanogr. 38, 133–145 (2008).

ADS  Article  Google Scholar 

Manucharyan, G. E. & Timmermans, M.-L. Generation and separation of mesoscale eddies from surface ocean fronts. J. Phys. Oceanogr. 43, 2545–2562 (2013).

ADS  Article  Google Scholar 

Zhao, M. et al. Characterizing the eddy field in the Arctic Ocean halocline. J. Geophys. Res.: Oceans 119, 8800–8817 (2014).

ADS  Article  Google Scholar 

Spall, M. A. Frontogenesis, subduction, and cross-front exchange at upper ocean fronts. J. Geophys. Res. 100, 2543–2557 (1995).

ADS  Article  Google Scholar 

Thomas, L. N. & Joyce, T. M. Subduction on the northern and southern flanks of the Gulf Stream. J. Phys. Oceanogr. 40, 429–438 (2010).

ADS  Article  Google Scholar 

Hoskins, B. J. & Bretherton, F. P. Atmospheric frontogenesis models: Mathematical formulation and solutions. J. Atmos. Sci. 29, 11–37 (1972).

ADS  Article  Google Scholar 

Thomas, L. N., Lee, C. M. & Yoshikawa, Y. The subpolar front of the Japan/East Sea. Part II: Inverse method for determining the frontal vertical circulation. J. Phys. Oceanogr. 40, 3–25 (2010).

ADS  Article  Google Scholar 

Pallàs-Sanz, E., Johnston, T. M. S. & Rudnick, D. L. Frontal dynamics in a California Current system shallow front: 2. mesoscale vertical velocity. J. Geophys. Res. 115, C12068 https://doi.org/10.1029/2010JC006474 (2010).

Gula, J., Molemaker, M. J. & McWilliams, J. C. Topographic generation of submesoscale centrifugal instability and energy dissipation. Nature Commun. 7, 1–7 (2016).

Article  CAS  Google Scholar 

Wenegrat, J. O. & Thomas, L. N. Centrifugal and symmetric instability during Ekman adjustment of the bottom boundary layer. J. Phys. Oceanogr. 50, 1793–1812 (2020).

ADS  Article  Google Scholar 

Shakespeare, C. J. Curved density fronts: Cyclogeostrophic adjustment and frontogenesis. J. Phys. Oceanogr. 46, 3193–3207 (2016).

ADS  Article  Google Scholar 

Buckingham, C. E., Gula, J. & Carton, X. The role of curvature in modifying frontal instabilities. part I: review of theory and presentation of a nondimensional instability criterion. J. Phys. Oceanogr. 51, 299–315 (2021).

ADS  Article  Google Scholar 

Buckingham, C. E., Gula, J. & Carton, X. The role of curvature in modifying frontal instabilities. part II: application of the criterion to curved density fronts at low richardson numbers. J. Phys. Oceanogr. 51, 317–341 (2021).

ADS  Article  Google Scholar 

Münchow, A. & Carmack, E. C. Synoptic flow and density observations near an Arctic shelf break. J. Phys. Oceanogr. 27, 1402–1419 (1997).

ADS  Article  Google Scholar 

Peralta-Ferriz, C. & Woodgate, R. A. Seasonal and interannual variability of pan-Arctic surface mixed layer properties from 1979 to 2012 from hydrographic data, and the dominance of stratification for multiyear mixed layer depth shoaling. Prog. Oceanogr. 134, 19–53 (2015).

ADS  Article  Google Scholar 

Boury, S. et al. Whither the chukchi slope current? J. Phys. Oceanogr. 50, 1717–1732 (2020).

ADS  Article  Google Scholar 

Meneghello, G., Marshall, J., Cole, S. T. & Timmermans, M.-L. Observational inferences of lateral eddy diffusivity in the halocline of the Beaufort Gyre. Geophys. Res. Lett. 44, 12–331 (2017).

Article  Google Scholar 

Steele, M. et al. Circulation of summer Pacific halocline water in the Arctic Ocean. J. Geophys. Res.: Oceans 109 (2004).

de la Vega, C., Jeffreys, R. M., Tuerena, R., Ganeshram, R. & Mahaffey, C. Temporal and spatial trends in marine carbon isotopes in the Arctic Ocean and implications for food web studies. Glob. Change Biol. 25, 4116–4130 (2019).

ADS  Article  Google Scholar 

Martin, J., Dumont, D. & Tremblay, J.-É. Contribution of subsurface chlorophyll maxima to primary production in the coastal Beaufort Sea (Canadian Arctic): A model assessment. J. Geophys. Res. 118, 5873–5886 (2013).

ADS  Article  Google Scholar 

Brown, Z. W. et al. Characterizing the subsurface chlorophyll a maximum in the Chukchi Sea and Canada Basin. Deep Sea Res. Part II: Topical Stud. Oceanogr. 118, 88–104 (2015).

CAS  Article  Google Scholar 

Bates, N. R., Hansell, D. A., Moran, S. B. & Codispoti, L. A. Seasonal and spatial distribution of particulate organic matter (POM) in the Chukchi and Beaufort Seas. Deep-Sea Res. II 52, 3324–3343 (2005).

ADS  Article  Google Scholar 

Mathis, J. T., Pickart, R. S., Hansell, D. A., Kadko, D. & Bates, N. R. Eddy transport of organic carbon and nutrients from the Chukchi Shelf: impact on the upper halocline of the western Arctic Ocean. J. Geophys. Res. 112, C05011 https://doi.org/10.1029/2006JC003899 (2007).

Ou, H. W. & Gordon, A. L. Spin-down of baroclinic eddies under sea ice. J. Geophys. Res.: Oceans 91, 7623–7630 (1986).

ADS  Article  Google Scholar 

Aagaard, K., Coachman, L. & Carmack, E. On the halocline of the Arctic Ocean. Deep-Sea Res. Part A 28, 529–545 (1981).

ADS  Article  Google Scholar 

Lewis, K., Van Dijken, G. & Arrigo, K. Changes in phytoplankton concentration now drive increased Arctic Ocean primary production. Science 369, 198–202 (2020).

CAS  PubMed  Article  Google Scholar 

Wang, M. & Overland, J. E. Projected future duration of the sea-ice-free season in the Alaskan Arctic. Prog. Oceanogr. 136, 50–59 (2015).

ADS  Article  Google Scholar 

Polyakov, I. V. et al. Weakening of cold halocline layer exposes sea ice to oceanic heat in the eastern Arctic Ocean. J. Clim. 33, 8107–8123 (2020).

ADS  Article  Google Scholar 

Jackson, J., Allen, S., Carmack, E. & McLaughlin, F. Suspended particles in the Canada Basin from optical and bottle data, 2003–2008. Ocean Sci. Discuss. 7, 799–813 (2010).

ADS  Article  Google Scholar 

McLaughlin, F. A. & Carmack, E. C. Deepening of the nutricline and chlorophyll maximum in the Canada Basin interior, 2003–2009. Geophys. Res. Lett. 37, L24602 https://doi.org/10.1029/2010GL045459 (2010).

Timmermans, M.-L. & Jayne, S. R. The Arctic Ocean spices up. J. Phys. Oceanogr. 46, 1277–1284 (2016).

ADS  Article  Google Scholar 

Pinkel, R. et al. The wirewalker: a vertically profiling instrument carrier powered by ocean waves. J. Atmos. Ocean. Technol. 28, 426–435 (2011).

ADS  Article  Google Scholar 

Todd, R. E., Rudnick, D. L. & Davis, R. E. Monitoring the greater San Pedro Bay region using autonomous underwater gliders during fall of 2006. J. Geophys. Res.: Oceans 114, C06001 https://doi.org/10.1029/2008JC005086 (2009).

Wood, K. R. et al. Results of the first Arctic Heat Open Science experiment. Bull. Am. Meteorol. Soc. 99, 513–520 (2018).

ADS  Article  Google Scholar 

Fairall, C. W., Bradley, E. F., Hare, J. E., Grachev, A. A. & Edson, J. B. Bulk parameterization of air-sea fluxes: updates and verification for the COARE algorithm. J. Clim. 16, 571–591 (2003).

ADS  Article  Google Scholar 

Edson, J. B. et al. On the exchange of momentum over the open ocean. J. Phys. Oceanogr. 43, 1589–1610 (2013).

ADS  Article  Google Scholar 

Alford, M. H. & Gregg, M. C. Near-inertial mixing: Modulation of shear, strain and microstructure at low latitude. J. Geophys. Res. 106, 16947–16968 (2001).

ADS  Article  Google Scholar 

MacKinnon, J. A. & Gregg, M. C. Mixing on the late-summer New England Shelf—Solibores, shear, and stratification. J. Phys. Oceanogr. 33, 1476–1492 (2003).

ADS  Article  Google Scholar 

Osborn, T. R. & Crawford, W. R. An airfoil probe for measuring turbulent velocity fluctuations in water. In Dobson, F., Hasse, L. & Davis, R. (eds.) Air–Sea Interactions: Instruments and Methods, pp. 369–386 (Plenum, New York, 1980).

Thomson, J. Wave breaking dissipation observed with SWIFT drifters. J. Atmos. Ocean. Technol. 29, 1866–1882 (2012).

ADS  Article  Google Scholar 

Download references

Support for this work was provided by the US Office of Naval Research Stratified Ocean Dynamics of the Arctic program (grant numbers N000141512903, N000141612378, N000141612377, N000141612379, N0001416123450, N000141612360, N000141612349, N000141812007, N000141812475, and N000141912514). Additional support for biogeochemistry sampling was provided by UK (NERC) and Germany (BMBF) through the Changing Arctic Ocean Program’s ARISE (NE/P006035/1, NE/P006000/2), PEANUTS (NE/R01275X/1, NE/R012547/2, and BMBF 03F0804) projects and the UK-France PhD program DGA/Dstl. Float deployments and hydrographic data compilations were supported in part by North Pacific Research Board grants A91-99a and A91-00a. The deployment of autonomous ocean profilers was supported by ONR, NOAA Research, and the Joint Institute for the Study of the Atmosphere and Ocean (JISAO) under NOAA Cooperative Agreement NA15OAR4320063. We are grateful to the engineers within our research groups and the captain and crew of the R/V Sikuliaq for facilitating this work, particularly Mike Goldin, Ethan Roth, Cris Seaton, and Paul St.Onge.

Scripps Institution of Oceanography, University of California San Diego, San Diego, CA, USA

Jennifer A. MacKinnon, Matthew H. Alford, Nicole Couto & Andrew J. Lucas

University of Alaska Fairbanks, Fairbanks, Alaska, USA

Harper L. Simmons, Seth L. Danielson & Thilo Klenz

Center for Southeastern Tropical Advanced Remote Sensing, University of Miami, Miami, FL, USA

John Hargrove, Hans C. Graber & Björn Lund

Applied Physics Laboratory, University of Washington, Seattle, WA, USA

Jim Thomson, Samuel D. Brenner, John Guthrie, Craig M. Lee, Luc Rainville & Madison M. Smith

Department of Mechanical Engineering, Massachusetts Institute of Technology, Cambridge, MA, USA

Thomas Peacock & Chanhyung Jeon

National Oceanography Centre, Liverpool, UK

Benjamin I. Barton & Joanne E. Hopkins

Courant Institute of Mathematical Sciences, New York University, New York, NY, USA

Samuel Boury

Woods Hole Oceanographic Institution, Woods Hole, MA, USA

Elizabeth C. Fine & Steven R. Jayne

Department of Oceanography, Pusan National University, Busan, South Korea

Chanhyung Jeon

School of Ocean Sciences, Bangor University, Bangor, Wales, UK

Yueng-Djern Lenn

School of Environmental Sciences, University of Liverpool, Liverpool, UK

Claire Mahaffey & Louisa Norman

Department of Earth System Science, Stanford University, Stanford, CA, USA

Leif N. Thomas

Alfred Wegener Institute, Bremerhaven, Germany

Sinhué Torres-Valdés

Joint Institute for the Study of the Atmosphere and Ocean, University of Washington, Seattle, WA, USA

Kevin R. Wood

NOAA Pacific Marine Environmental Laboratory, Seattle, WA, USA

Kevin R. Wood

Experimental design for this survey was led by J.A.M., H.L.S., J.H., J.T., and T.P. Instrument design and development for the primary survey was led by M.H.A., A.J.L., and J.T. Data acquisition on the primary cruise was carried out by J.A.M., H.L.S., J.H., J.T., T.P., B.I.B., S.B., N.C., S.D., E.C.F., J.G., T.K., C.J., and M.M.S. Satellite data acquisition and analysis were carried out by J.H., H.C.G., and B.L. Mooring development, deployment, and analysis were carried out by C.L., L.R., and S.D.B. Glider and drifter deployments and analysis were led by S.L.D., S.R.J., and K.R.W. Biogeochemistry bottle samples were collected, processed and analyzed by B.I.B., Y.D.L., J.E.H., C.M., L.N., and S.T.V. Instability calculations were carried out by J.A.M. and L.N.T. C.L. is the lead PI for the SODA project. All authors contributed to writing and editing the manuscript.

The authors declare no competing interests.

Peer review information Nature Communications thanks Satoshi Kimura and other, anonymous, reviewers for their contributions to the peer review of this work.

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

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 license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license 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 license, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.