Global observations show that the ocean lost approximately 2% of its oxygen inventory over the past five decades^1,2,3, with important implications for marine ecosystems^4,5. The rate of change varies regionally, with northwest Atlantic coastal waters showing a long-term drop^6,7 that vastly outpaces the global and North Atlantic basin mean deoxygenation rates^5,8. However, past work has been unable to differentiate the role of large-scale climate forcing from that of local processes. Here, we use hydrographic evidence to show that a Labrador Current retreat is playing a key role in the deoxygenation on the northwest Atlantic shelf. A high-resolution global coupled climate–biogeochemistry model⁹ reproduces the observed decline of saturation oxygen concentrations in the region, driven by a retreat of the equatorward-flowing Labrador Current and an associated shift towards more oxygen-poor subtropical waters on the shelf. The dynamical changes underlying the shift in shelf water properties are correlated with a slowdown in the simulated Atlantic Meridional Overturning Circulation (AMOC)¹⁰. Our results provide strong evidence that a major, centennial-scale change of the Labrador Current is underway, and highlight the potential for ocean dynamics to impact coastal deoxygenation over the coming century.

Lateral stirring is a basic oceanographic phenomenon affecting the distribution of physical, chemical, and biological fields. Eddy stirring at scales on the order of 100 km (the mesoscale) is fairly well understood and explicitly represented in modern eddy-resolving numerical models of global ocean circulation. The same cannot be said for smaller-scale stirring processes. Here, the authors describe a major oceanographic field experiment aimed at observing and understanding the processes responsible for stirring at scales of 0.1–10 km. Stirring processes of varying intensity were studied in the Sargasso Sea eddy field approximately 250 km southeast of Cape Hatteras. Lateral variability of water-mass properties, the distribution of microscale turbulence, and the evolution of several patches of inert dye were studied with an array of shipboard, autonomous, and airborne instruments. Observations were made at two sites, characterized by weak and moderate background mesoscale straining, to contrast different regimes of lateral stirring. Analyses to date suggest that, in both cases, the lateral dispersion of natural and deliberately released tracers was O(1) m2 s–1 as found elsewhere, which is faster than might be expected from traditional shear dispersion by persistent mesoscale flow and linear internal waves. These findings point to the possible importance of kilometer-scale stirring by submesoscale eddies and nonlinear internal-wave processes or the need to modify the traditional shear-dispersion paradigm to include higher-order effects. A unique aspect of the Scalable Lateral Mixing and Coherent Turbulence (LatMix) field experiment is the combination of direct measurements of dye dispersion with the concurrent multiscale hydrographic and turbulence observations, enabling evaluation of the underlying mechanisms responsible for the observed dispersion at a new level.

Oceanic fronts are dynamically active regions of the global ocean that support upwelling and downwelling with significant implications for phytoplankton production and export. However (on time scales the inertial time scale), the vertical velocity is 10³–10⁴ times weaker than the horizontal velocity and is difficult to observe directly. Using intensive field observations in conjunction with a process study ocean model, we examine vertical motion and its effect on phytoplankton fluxes at multiple spatial horizontal scales in an oligotrophic region in the Western Mediterranean Sea. The mesoscale ageostrophic vertical velocity (∼10 m/day) inferred from our observations shapes the large‐scale phytoplankton distribution but does not explain the narrow (1–10 km wide) features of high chlorophyll content extending 40–60 m downward from the deep chlorophyll maximum. Using modeling, we show that downwelling submesoscale features concentrate 80% of the downward vertical flux of phytoplankton within just 15% of the horizontal area. These submesoscale spatial structures serve as conduits between the surface mixed layer and pycnocline and can contribute to exporting carbon from the sunlit surface layers to the ocean interior.

The challenges associated with meso- and submesoscale variability (between 1-100 km) require high-resolution observations and integrated approaches. Here we describe a major oceanographic experiment designed to capture the intense but transient vertical motions associated with mesoscale and submesoscale features in an area characterized by strong fronts. Finescale processes were studied in the eastern Alboran Sea (Western Mediterranean) about 400 km east of the Strait of Gibraltar, a relatively sparsely sampled area. In-situ systems were coordinated with satellite data to provide a full description of the physical and biogeochemical variability. Hydrographic data confirmed the presence of an intense salinity front formed by the confluence of Atlantic Waters, entering from Gibraltar, with the local Mediterranean waters. The drifters coherently followed the northeastern limb of an anticyclonic gyre. Near real time data from acoustic current meter data profiler showed consistent patterns with currents of up to 1m/s in the southern part of the sampled domain. High-resolution glider data revealed submesoscale structures with tongues of chlorophyll-a and oxygen associated with the frontal zone. Numerical results show large vertical excursions of tracers that could explain the subducted tongues and filaments captured by ocean gliders. A unique aspect of AlborEx is the combination of high-resolution synoptic measurements of vessel-based measurements, autonomous sampling, remote sensing and modeling, enabling the evaluation of the underlying mechanisms responsible for the observed distributions and biogeochemical patchiness. The main findings point to the importance of fine-scale processes enhancing the vertical exchanges between the upper ocean and the ocean interior.

Earth System Models increasingly include ocean biogeochemistry models in order to predict changes in ocean carbon storage, hypoxia, and biological productivity under climate change. However, state‐of‐the‐art ocean biogeochemical models include many advected tracers, that significantly increase the computational resources required, forcing a trade‐off with spatial resolution. Here, we compare a state‐of‐the art model with 30 prognostic tracers (TOPAZ) with two reduced‐tracer models, one with 6 tracers (BLING), and the other with 3 tracers (miniBLING). The reduced‐tracer models employ parameterized, implicit biological functions, which nonetheless capture many of the most important processes resolved by TOPAZ. All three are embedded in the same coupled climate model. Despite the large difference in tracer number, the absence of tracers for living organic matter is shown to have a minimal impact on the transport of nutrient elements, and the three models produce similar mean annual preindustrial distributions of macronutrients, oxygen, and carbon. Significant differences do exist among the models, in particular the seasonal cycle of biomass and export production, but it does not appear that these are necessary consequences of the reduced tracer number. With increasing CO₂, changes in dissolved oxygen and anthropogenic carbon uptake are very similar across the different models. Thus, while the reduced‐tracer models do not explicitly resolve the diversity and internal dynamics of marine ecosystems, we demonstrate that such models are applicable to a broad suite of major biogeochemical concerns, including anthropogenic change. These results are very promising for the further development and application of reduced‐tracer biogeochemical models that incorporate “sub‐ecosystem‐scale” parameterizations.

The long‐term trend of increasing phytoplankton net primary production (NPP) in the Arctic correlates with increasing light penetration due to sea ice loss. However, recent studies suggest that enhanced stormy wind mixing may also play a significant role enhancing NPP. Here, we isolate the role of sea ice and stormy winds (hereafter high‐frequency winds) using an eddy‐permitting ice‐ocean‐biogeochemical model configured for the North Atlantic and the Arctic. In the model, the presence of high‐frequency winds stimulates nutrient upwelling by producing an earlier and longer autumn‐winter mixing period with deeper mixing layer. The early onset of autumn mixing results in nutrients being brought‐up to near‐surface waters before the light becomes the dominant limiting factor, which leads to the autumn bloom. The enhanced mixing results in higher nutrient concentrations in spring and thus a large spring bloom. The model also shows significant iron limitation in the Labrador Sea, which is intensified by high‐frequency winds. The effect of sea ice loss on NPP was found to be regionally dependent on the presence of high‐frequency winds. This numerical study suggests high‐frequency winds play significant role increasing NPP in the Arctic and sub‐Arctic by alleviating phytoplankton nutrient limitation and that the isolated effect of sea ice loss on light plays a comparatively minor role.