The Sub-Mesoscale Ocean Dynamics Experiment (S-MODE) is a NASA Earth Ventures Suborbital investigation designed to test the hypothesis that oceanic frontogenesis and the kilometer-scale ("submesoscale") instabilities that accompany it make important contributions to vertical exchange of climate and biological variables in the upper ocean. These processes have been difficult to resolve in observations, making model validation challenging. A necessary step toward testing the hypothesis was to make accurate measurements of upper-ocean velocity fields over a broad range of scales and to relate them to the observed variability of vertical transport and surface forcing. A further goal was to examine the relationship between surface velocity, temperature, and chlorophyll measured by remote sensing and their depth-dependent distributions, within and beneath the surface boundary layer. To achieve these goals, we used aircraft-based remote sensing, satellite remote sensing, ships, drifter deployments, and a fleet of autonomous vehicles. The observational component of S-MODE consisted of three campaigns, all conducted in the Pacific Ocean approximately 100-km west of San Francisco during 2021–23 fall and spring. S-MODE was enabled by recent developments in remote sensing technology that allowed operational airborne observation of ocean surface velocity fields and by advances in autonomous instrumentation that allowed coordinated sampling with dozens of uncrewed vehicles at sea. The coordinated use of remote sensing measurements from three aircraft with arrays of remotely operated vehicles and other in situ measurements is a major novelty of S-MODE. All S-MODE data are freely available, and their use is encouraged.

Surface drifters (SWIFTs) equipped with down-looking high-resolution acoustic doppler current profilers (ADCPs) were used to estimate the turbulent kinetic energy (TKE) dissipation rate
(ε) within highly stratified diurnal warm layers (DWLs) in the Southern California Bight. Over a 10-day period, five instances of DWLs were observed with strong surface temperature anomalies up to 3°C and velocity anomalies up to 0.3 m s^-1. Profiles of
ε in the upper 5 m suggest turbulence is strongly modulated by the DWL stratification. Burst-averaged (8.5 min) ε is stronger than predicted by law-of-the-wall boundary layer scaling within the DWLs and suppressed below. Predictions for ε within the DWLs are improved by a shear-production scaling using observed shear and linearly decaying turbulent stress. However, ε is still under-predicted. Examination of the un-averaged acoustic backscatter data suggests elevated ε is related to the presence of turbulent structures in the DWLs which span the layer height and strongly modulate TKE. Evolution in the bulk Richardson number each day suggests the DWLs become unstable to layer-scale overturning and entrainment each afternoon, thus the turbulent structures may result from shear-driven instability. This interpretation is supported by a conditional average of the data during a burst characterized by strongly periodic structures. The structures resemble high-frequency internal waves with strong asymmetry in the along-flow direction (steepening) which suggests they are unstable. Coincident asymmetric patterns in upwelling/downwelling and corresponding regions of strong vertical convergence/divergence suggest that both vertical transport and local TKE generation are plausible sources of elevated ε in the DWLs.

Surface drifters (SWIFTs) equipped with down-looking high-resolution acoustic doppler current profilers (ADCPs) were used to estimate the turbulent kinetic energy (TKE) dissipation rate (ε) within highly stratified diurnal warm layers (DWLs) in the Southern California Bight. Over a 10-day period, five instances of DWLs were observed with strong surface temperature anomalies up to 3°C and velocity anomalies up to 0.3 m s^-1. Profiles of ε in the upper 5 m suggest turbulence is strongly modulated by the DWL stratification. Burst-averaged (8.5 min) ε is stronger than predicted by law-of-the-wall boundary layer scaling within the DWLs and suppressed below. Predictions for ε within the DWLs are improved by a shear-production scaling using observed shear and linearly decaying turbulent stress. However, ε is still under-predicted. Examination of the un-averaged acoustic backscatter data suggests elevated ε is related to the presence of turbulent structures in the DWLs which span the layer height and strongly modulate TKE. Evolution in the bulk Richardson number each day suggests the DWLs become unstable to layer-scale overturning and entrainment each afternoon, thus the turbulent structures may result from shear-driven instability. This interpretation is supported by a conditional average of the data during a burst characterized by strongly periodic structures. The structures resemble high-frequency internal waves with strong asymmetry in the along-flow direction (steepening) which suggests they are unstable. Coincident asymmetric patterns in upwelling/downwelling and corresponding regions of strong vertical convergence/divergence suggest that both vertical transport and local TKE generation are plausible sources of elevated ε in the DWLs.

A neutrally buoyant float deployed in an atmospherically driven turbulent ocean boundary layer on the dense side of a submesoscale front was repeatedly carried across the boundary layer by the turbulence and then trapped beneath the slumping front. Lagrangian particles in a large-eddy simulation of a similar baroclinically unstable front forced by surface cooling move along convergent surface filaments toward filament junctions. They are also caught by convective plumes that downwell them at speeds similar to those of the float. Subsequently, some are trapped in the pycnocline by frontal slumping due to ageostrophic secondary frontal circulations. In both observations and simulations, boundary layer turbulence and frontal circulations work together to trap and subduct particles from the mixed layer. The small-scale boundary layer motions move them vertically within the boundary layer and larger, submesoscale frontal circulations move them laterally out of the boundary layer and under the slumping fronts.

As our planet warms, Arctic sea ice coverage continues to decline, resulting in complex feedbacks with the climate system. The core objective of NASA's Salinity and Stratification at the Sea Ice Edge (SASSIE) mission is to understand how ocean salinity and near-surface stratification affect upper-ocean heat content and thus sea ice freeze and melt. SASSIE specifically focuses on the formation of Arctic Sea ice in autumn. The SASSIE field campaign in 2022 collected detailed observations of upper-ocean properties and meteorology near the sea ice edge in the Beaufort Sea using ship-based and piloted and drifting assets. The observations collected during SASSIE include vertical profiles of stratification up to the sea surface, air–sea fluxes, and ancillary measurements that are being used to better understand the role of salinity in coupled Arctic air–sea–ice processes. This publication provides a detailed overview of the activities during the 2022 SASSIE campaign and presents the publicly available datasets generated by this mission (available at https://podaac.jpl.nasa.gov/SASSIE, last access: 29 May 2024; DOIs for individual datasets in the "Data availability" section), introducing an accompanying repository that highlights the numerical routines used to generate the figures shown in this work.

We perform a statistical characterization of the 2016–2017 SPURS-2 field campaign in situ data and coincident satellite data spanning 8°–12°N, 120°–130°W to quantify the spatial and temporal scales of variability of rain and near-surface salinity in the Eastern Pacific Fresh Pool. Observations of rain rate and near-surface to surface salinity are obtained from ships, moorings, autonomous platforms, and satellite remote sensing: Integrated Multi-satellitE Retrievals for GPM (IMERG); and Soil Moisture Active Passive (NASA SMAP L3 V5). The integral length and time scales of rain and near-surface salinity vary seasonally. In the rainy season (August–October) when the Intertropical Convergence Zone (ITCZ) migrates over the SPURS-2 study site, the integral time scales of rain were about 30–60 min and those of near-surface salinity were closer to that of the rain, 1–2 days, indicating forcing by rain. Meanwhile, the zonal integral length scale of in situ near-surface salinity was twice as large as the meridional scale (50 vs. 20 km), consistent with the ITCZ's zonally-propagating and -organized rain features. The magnitude and seasonal variation of the sea surface salinity integral time scale were not captured by SMAP since the rainy ITCZ-period scales were smaller than SMAP resolution (70 km, 8-day running mean). In the dry season (February–May), the in situ rain integral time scale reduced to less than 30 min while that of the near-surface salinity increased to 1–5 days, the ocean mesoscale. IMERG overestimated the rain integral time scale by a factor of two to ten in both seasons.

We show that the high seas are colonized by a diverse array of coastal species, which survive and reproduce in the open ocean, contributing strongly to its floating community composition. Analysis of rafting plastic debris in the eastern North Pacific Subtropical Gyre revealed 37 coastal invertebrate taxa, largely of Western Pacific origin, exceeding pelagic taxa richness by threefold. Coastal taxa, including diverse taxonomic groups and life history traits, occurred on 70.5% of debris items. Most coastal taxa possessed either direct development or asexual reproduction, possibly facilitating long-term persistence on rafts. Our results suggest that the historical lack of available substrate limited the colonization of the open ocean by coastal species, rather than physiological or ecological constraints as previously assumed. It appears that coastal species persist now in the open ocean as a substantial component of a neopelagic community sustained by the vast and expanding sea of plastic debris.

Vertical transport pathways in the ocean are still only partially understood despite their importance for biogeochemical, pollutant, and climate applications. Detailed measurements of a submesoscale frontal jet in the Alboran Sea (Mediterranean Sea) during a period of highly variable winds were made using cross-frontal velocity, density sections and dense arrays of surface drifters deployed across the front. The measurements show divergences as large as ±f implying vertical velocities of order 100 m/day for a ≈ 20 m thick surface layer. Over the 20 hr of measurement, the divergences made nearly one complete oscillation, suggesting an important role for near-inertial oscillations. A wind-forced slab model modified by the observed background frontal structure and with initial conditions matched to the data produces divergence oscillations and pattern compatible with that observed. Significant differences, though, are found in terms of mean divergence, with the data showing a prevalence of negative, convergent values. Despite the limitations in data sampling and model uncertainties, this suggests the contribution of other dynamical processes. Turbulent boundary layer processes are discussed, as a contributor to enhance the observed convergent phase. Water mass properties suggest that symmetric instabilities might also be present but do not play a crucial role, while downward stirring along displaced isopycnals is observed.

Intensive hydrographic and horizontal velocity measurements collected in the Alboran Sea enabled us to diagnose the three-dimensional dynamics of a frontal system. The sampled domain was characterized by a 40 km diameter anticyclonic eddy, with an intense front on its eastern side, separating the Atlantic and Mediterranean waters. Here, we implemented a multi-variate variational analysis (VA) to reconstruct the hydrographic fields, combining the 1-km horizontal resolution of the Underway Conductivity-Temperature-Depth (CTD) system with information on the flow shape from the Acoustic Doppler Current Profiler velocities. One advantage of the VA is given by the physical constraint, which preserves fine-scale gradients better than the classical optimal interpolation (OI). A comparison between real drifter trajectories and virtual particles advected in the mapping quantified the improvements in the VA over the OI, with a 15% larger skill score. Quasi-geostrophic (QG) and semi-geostrophic (SG) omega equations enabled us to estimate the vertical velocity (w) which reached 40 m/day on the dense side of the front. How nutrients and other passive tracers leave the mixed-layer and subduct is estimated with 3D advection from the VA, which agreed with biological sampling from traditional CTD casts at two eddy locations. Downwelling warm filaments are further evidence of subduction, in line with the w from SG, but not with QG. SG better accounted for the along-isopycnal component of w in agreement with another analysis made on isopycnal coordinates. The multi-platform approach of this work and the use of variational methods improved the characterization and understanding of (sub)-mesoscale frontal dynamics.

Discoveries of persistent coastal species in the open ocean shift our understanding of biogeographic barriers. Floating plastic debris from pollution now supports a novel sea surface community composed of coastal and oceanic species at sea that might portend significant ecological shifts in the marine environment.

A crucial region of the ocean surface boundary layer (OSBL) is the strongly-sheared and -stratified transition layer (TL) separating the mixed layer from the upper pycnocline, where a diverse range of waves and instabilities are possible. Previous work suggests that these different waves and instabilities will lead to different OSBL behaviours. Therefore, understanding which physical processes occur is key for modelling the TL. Here we present observations of the TL from a Lagrangian float deployed for 73 days near Ocean Weather Station Papa (50°N, 145°W) during Fall 2018. The float followed the vertical motion of the TL, continuously measuring profiles across it using an ADCP, temperature chain and salinity sensors. The temperature chain made depth/time images of TL structures with a resolution of 6 cm and 3 seconds. These showed the frequent occurrence of very sharp interfaces, dominated by temperature jumps of O(1)°C over 6 cm or less. Temperature inversions were typically small (less than about 10 cm), frequent, and strongly-stratified; very few large overturns were observed. The corresponding velocity profiles varied over larger length scales than the temperature profiles. These structures are consistent with scouring behaviour rather than Kelvin–Helmholtz-type overturning. Their net effect, estimated via a Thorpe-scale analysis, suggests that these frequent small temperature inversions can account for the observed mixed layer deepening and entrainment flux. Corresponding estimates of dissipation, diffusivity, and heat fluxes also agree with previous TL studies, suggesting that the TL dynamics is dominated by these nearly continuous 10-cm scale mixing structures, rather than by less frequent larger overturns.

Vertical transport generated by mesoscale and submesoscale flows plays a key role in the exchange of physical and biogeochemical properties between the surface and the ocean interior. Using multiple simultaneous drifter observations, we compute spatial gradients of velocity to obtain estimates of the divergence field. Thanks to the fact that drifters were deployed at two different depths, we can observe the vertical dependence of divergence in the upper 15 m and estimate the associated vertical velocity. In this study, we estimate divergence and vertical velocity in a ~1-m/s semipermanent frontal jet in the Alboran Sea by making use of a multiplatform data set including 82 drifters, a Lagrangian float, and along-shiptrack profile timeseries of temperature and salinity.

A coordinated survey between a subsurface Lagrangian float and a ship-towed Triaxus profiler obtained detailed measurements of a restratifying surface intensified front (above 30 m) within the California Current System. The survey began as downfront winds incited mixing in the boundary layer. As winds relaxed and mixing subsided, the system entered a different dynamical regime as the front developed an overturning circulation with large vertical velocities that tilted isopycnals and stratified the upper ocean within a day. The horizontal buoyancy gradient was 1.5 x 10^-6 s^-2 and associated with vorticity, divergence, and strain that approached the Coriolis frequency. Estimates of vertical velocity from the Lagrangian float reached 1.2 x 10^-3 m s^-1. These horizontal gradients and vertical velocities were consistent with submesoscale dynamics that are distinct from the classic quasigeostrophic framework used to describe larger-scale flows. Vertical and horizontal gradients of velocity and buoyancy in the vicinity of the float revealed that sheared currents differentially advected the horizontal buoyancy gradient to increase vertical stratification. This was supported by analyses of temperature and salinity gradients that composed the horizontal and vertical stratification. Potential vorticity was conserved during restratification at 16 m, consistent with adiabatic processes. Conversely, potential vorticity near the surface (8 m) increased, highlighting the role of friction in modulating near-surface stratification. The observed increase in stratification due to these submesoscale processes was equivalent to a heat flux of 2000 W m^-2, which is an order-of-magnitude larger than the average observed surface heat flux of 100 W m^-2.

SPURS-2 (Salinity Processes in the Upper-ocean Regional Study 2) used the schooner Lady Amber, a small sailing research vessel, to deploy, service, maintain, and recover a variety of oceanographic and meteorological instruments in the eastern Pacific Ocean. Low operational costs allowed us to frequently deploy floats and drifters to collect data necessary for resolving the regional circulation of the eastern tropical Pacific. The small charter gave us the opportunity to deploy drifters in locations chosen according to current conditions, to recover and deploy various autonomous instruments in a targeted and adaptive manner, and to collect additional near-surface and atmospheric measurements in the remote SPURS-2 region.

An autonomous Lagrangian float equipped with a high-resolution acoustic Doppler current profiler observed the evolution of upper-ocean stratification and velocity in the Eastern Pacific Fresh Pool for over 100 days in August–November 2016. Although convective mixing homogenized the water column to 40 m depth almost every night, the combination of diurnal warming on clear days and rainfall on cloudy days routinely produced strong stratification in the upper 10 m. Whether due to thermal or freshwater effects, the initial strong stratification was mixed downward and incorporated in the bulk of the mixed layer within a few hours. Stratification cycling was associated with pronounced variability of ocean surface boundary layer turbulence and vertical shear of wind-driven (Ekman) currents. Decoupled from the bulk of the mixed layer by strong stratification, warm and fresh near-surface waters were rapidly accelerated by wind, producing the well-known "slippery layer" effect, and leading to a strong downwind near-surface distortion of the Ekman profile. A case study illustrates the ability of the new generation of Lagrangian floats to measure rapidly evolving temperature, salinity, and velocity, including turbulent and internal wave components. Quantitative interpretation of the results remains a challenge, which can be addressed with high-resolution numerical modeling, given sufficiently accurate air-sea fluxes.

The Lagrangian Submesoscale Experiment (LASER) was designed to study surface flows during winter conditions in the northern Gulf of Mexico. More than 1000 mostly biodegradable drifters were launched. The drifters consisted of a surface floater extending 5 cm below the surface, containing the satellite tracking system, and a drogue extending 60 cm below the surface, hanging beneath the floater on a flexible tether. On some floats, the drogue separated from the floater during storms. This paper describes methods to detect drogue loss based on two properties that distinguish drogued from undrogued drifters. First, undrogued drifters often flip over, pointing their satellite antenna downward and thus intermittently reducing the frequency of GPS fixes. Second, undrogued drifters respond to wind forcing more than drogued drifters. A multistage analysis is used: first, two properties are used to create a preliminary drifter classification; then, the motion of each unclassified drifter is compared to that of its classified neighbors in an iterative process for nearly all of the drifters. The algorithm classified drifters with a known drogue status with an accuracy of virtually 100%. Drogue loss times were estimated with a precision of less than 0.5 and 3 h for 60% and 85% of the drifters, respectively. An estimated 40% of the drifters lost their drogues in the first 7 weeks, with drogue loss coinciding with storm events, particularly those with steep waves. Once the drogued and undrogued drifters are classified, they can be used to quantify the differences in material dispersion at different depths.

This paper describes the instrumentation and techniques for long-term targeted observation of the centimeter-scale velocity structure within the oceanic surface boundary layer, made possible by the recent developments in capabilities of autonomous platforms and self-contained pulse-coherent acoustic Doppler current profilers (ADCPs). Particular attention is paid to the algorithms of ambiguity resolution ("unwrapping") of pulse-coherent Doppler velocity measurements. The techniques are demonstrated using the new Nortek Signature1000 ADCP mounted on a Lagrangian float, a combination shown to be capable of observing ocean turbulence in a number of recent studies. Statistical uncertainty of the measured velocities in relation to the ADCP setup is also evaluated. Described techniques and analyses should be broadly applicable to other autonomous and towed applications of pulse-coherent ADCPs.

Plastics and spilled oil pose a critical threat to marine life and human health. As a result of wind forcing and wave motions, theoretical and laboratory studies predict very strong velocity variation with depth over the upper few centimeters of the water column, an observational blind spot in the real ocean. Here we present the first-ever ocean measurements of the current vector profile defined to within 1 cm of the free surface. In our illustrative example, the current magnitude averaged over the upper 1 cm of the ocean is shown to be nearly four times the average over the upper 10 m, even for mild forcing. Our findings indicate that this shear will rapidly separate pieces of marine debris which vary in size or buoyancy, making consideration of these dynamics essential to an improved understanding of the pathways along which marine plastics and oil are transported.

Ocean currents move material released on the ocean surface away from the release point and, over time, spread it over an increasingly large area. However, observations also show high concentrations of the material even after significant spreading. This work examines a mechanism for creating such concentrations: downwelling of water at the boundaries of different water masses concentrates floating material at this boundary. Hundreds of satellite-tracked drifters were released near the site of the 2010 Deepwater Horizon oil spill. Surprisingly, most of these gathered into a single cluster less than 100 m in size, dramatically demonstrating the strength of this mechanism.

The Salinity Processes in the Upper-ocean Regional Study (SPURS) aims to understand the patterns and variability of sea surface salinity. In order to capture the wide range of spatial and temporal scales associated with processes controlling salinity in the upper ocean, research vessels delivered autonomous instruments to remote sites, one in the North Atlantic and one in the Eastern Pacific. Instruments sampled for one complete annual cycle at each of these two sites, which are subject to contrasting atmospheric forcing. The SPURS field programs coordinated sampling from many different platforms, using a mix of Lagrangian and Eulerian approaches. This article discusses the motivations, implementation, and first results of the SPURS-1 and SPURS-2 programs.

The Arabian Sea circulation is forced by strong monsoonal winds and is characterized by vigorous seasonally reversing currents, extreme differences in sea surface salinity, localized substantial upwelling, and widespread submesoscale thermohaline structures. Its complicated sea surface temperature patterns are important for the onset and evolution of the Asian monsoon. This article describes a program that aims to elucidate the role of upper-ocean processes and atmospheric feedbacks in setting the sea surface temperature properties of the region. The wide range of spatial and temporal scales and the difficulty of accessing much of the region with ships due to piracy motivated a novel approach based on state-of-the-art autonomous ocean sensors and platforms. The extensive data set that is being collected, combined with numerical models and remote sensing data, confirms the role of planetary waves in the reversal of the Somali Current system. These data also document the fast response of the upper equatorial ocean to monsoon winds through changes in temperature and salinity and the connectivity of the surface currents across the northern Indian Ocean. New observations of thermohaline interleaving structures and mixing in setting the surface temperature properties of the northern Arabian Sea are also discussed.

Where, when, and what to sample, and how to optimally reach the sampling locations, are critical questions to be answered by autonomous and Lagrangian platforms and sensors. For a reproducible scientific sampling approach, answers should be quantitative and provided using fundamental principles. This article reviews concepts and recent progress toward this principled approach, focusing on reachability, path planning, and adaptive sampling, and presents results of a real-time forecasting and planning experiment completed during February–April 2017 for the Northern Arabian Sea Circulation-autonomous research program. The predictive skill, layered fields, and uncertainty estimates obtained using the MIT MSEAS multi-resolution ensemble ocean modeling system are first studied. With such inputs, deterministic and probabilistic three-dimensional reachability forecasts issued daily for gliders and floats are then showcased and validated. Finally, a Bayesian adaptive sampling framework is shown to forecast in real time the observations that are most informative for estimating classic ocean fields and also secondary variables such as Lagrangian coherent structures.

We propose a novel method for constructing a pseudo-synoptic view of estuarine features from non-synoptic observations captured by mobile platforms. The model-aided Lagrangian interpretation (MALI) method is based on relocating observations to a common reference moment in time along three-dimensional Lagrangian trajectories derived from a numerical model of estuarine circulation. The method relies on the model skill to capture large-scale circulation features, and on high-resolution in situ observations to characterize small-scale hydrographic structure. We demonstrate our technique by applying MALI to autonomous underwater vehicle observations in the Columbia River estuary, with the aid of a validated unstructured-grid finite-element numerical simulation. The method can be readily adapted to a broader range of environments, observational platforms, and model-data combinations.

The Gulf Stream is a major conduit of warm surface water from the tropics to the subpolar North Atlantic. Here we observe and simulate a submesoscale (<20 km) mechanism by which the Gulf Stream exchanges water with subpolar water to the north. Along isopycnals, the front has a sharp compensated temperature-salinity contrast, with distinct mixed water between the two water masses 2 and 4 km wide. This mixed water does not increase downstream despite substantial energy available for mixing. A series of streamers detrain this water at the crest of meanders. Subpolar water replaces the mixed water and resharpens the front. The water mass exchange accounts for a northward flux of salt of 0.5–2.5 psu m² s^-1, (large-scale diffusivity O (100 m² s^-1)). This is similar to bulk-scale flux estimates of 1.2 psu m² s^-1 and supplies fresher water to the Gulf Stream required for the production of 18° subtropical mode water.

The passage of a winter storm over the Gulf Stream observed with a Lagrangian float and hydrographic and velocity surveys provided a unique opportunity to study how the interaction of inertial oscillations, the front, and symmetric instability (SI) shapes the stratification, shear, and turbulence in the upper ocean under unsteady forcing. During the storm, the rapid rise and rotation of the winds excited inertial motions. Acting on the front, these sheared motions modulate the stratification in the surface boundary layer. At the same time, cooling and downfront winds generated a symmetrically unstable flow. The observed turbulent kinetic energy dissipation exceeded what could be attributed to atmospheric forcing, implying SI drew energy from the front. The peak excess dissipation, which occurred just prior to a minimum in stratification, surpassed that predicted for steady SI turbulence, suggesting the importance of unsteady dynamics. The measurements are interpreted using a large-eddy simulation (LES) and a stability analysis configured with parameters taken from the observations. The stability analysis illustrates how SI more efficiently extracts energy from a front via shear production during periods when inertial motions reduce stratification. Diagnostics of the energetics of SI from the LES highlight the temporal variability in shear production but also demonstrate that the time-averaged energy balance is consistent with a theoretical scaling that has previously been tested only for steady forcing. As the storm passed and the winds and cooling subsided, the boundary layer restratified and the thermal wind balance was reestablished in a manner reminiscent of geostrophic adjustment.

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) m² 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.

The North Atlantic subtropical salinity maximum harbors the saltiest surface waters of the open world ocean. Subduction of these waters gives rise to Subtropical Underwater, spreading the high-salinity signature over the entire basin. The Salinity Processes in the Upper-ocean Regional Study (SPURS) is aimed at understanding the physics controlling the thermohaline structure in the salinity maximum region. A combination of moored and autonomous float observations is used here to describe the vertical water mass interleaving in the area. Seasonal intensification of interleaving in late spring and the abundance of small-scale thermohaline intrusions point to an important role for submesoscale processes in the initial subduction and subsequent evolution of Subtropical Underwater.

We examine the variability of near-surface salinity in a 10° x 10° region of the eastern North Atlantic, the location of the first part of the Salinity Processes in the Upper-ocean Regional Study (SPURS-1). The data used were collected over a two-year period, largely by a group of two types of profiling floats equipped with sensors that record high-resolution temperature and salinity measurements in the upper few meters of the water column. By comparing the SPURS-1 measurements to observations in the area from previous decades, we examine variability at time scales ranging from a few hours (mostly consisting of rainfall-driven decreases in salinity) to diurnal cycles in temperature and salinity, seasonal variability and the annual cycle, and finally to decadal-scale changes. The relationship of near-surface salinity to the hydrological cycle suggests a continuous spectrum of variability in this cycle from hours to decades.

Nearly all operational ocean models use air-sea fluxes and the ocean shear and stratification to estimate upper ocean boundary layer mixing rates. This approach implicitly parameterizes surface wave effects in terms of these inputs. Here, we test this assumption using parallel experiments in a lake with small waves and in the open ocean with much bigger waves. Under the same wind stress and adjusting for buoyancy flux, we find the mixed layer average turbulent vertical kinetic energy in the open ocean typically twice that in the lake. The increase is consistent with models of Langmuir turbulence, in which the wave Stokes drift, and not wave breaking, is the dominant mechanism by which waves energize turbulence in the mixed layer. Applying these same theories globally, we find enhanced mixing and deeper mixed layers resulting from the inclusion of Langmuir turbulence in the boundary layer parameterization, especially in the Southern Ocean.

Wave and wind measurements at Ocean Weather Station P (OWS-P, 50°N 145°W) are used to evaluate the equilibrium range of surface wave energy spectra. Observations are consistent with a local balance between wind input and breaking dissipation, as described by Philips (1985). The measurements include direct covariance wind stress estimates and wave breaking dissipation rate estimates during a 3 week research cruise to OWS-P. The analysis is extended to a wider range of conditions using observations of wave energy spectra and wind speed during a 2 year mooring deployment at OWS-P. At moderate wind speeds (5–15 m/s), mooring wave spectra are in agreement, within 5% uncertainty, with the forcing implied by standard drag laws and mooring wind measurements. At high wind speeds (>15 m/s), mooring wave spectra are biased low, by 13%, relative to the forcing implied by standard drag laws and mooring wind measurements. Deviations from equilibrium are associated with directionality and variations at the swell frequencies. A spectral wave hindcast accurately reproduces the mooring observations, and is used to examine the wind input.

A detailed view of upper ocean vorticity, divergence, and strain statistics was obtained by a two-vessel survey in the North Atlantic Mode Water region in winter 2012. Synchronous Acoustic Doppler Current Profiler sampling provided the first in situ estimates of the full velocity gradient tensor at O(1 km) scale without the usual mix of spatial and temporal aliasing. The observed vorticity distribution in the mixed layer was markedly asymmetric (skewness 2.5), with sparse strands of strong cyclonic vorticity embedded in a weak, predominantly anticyclonic background. Skewness of the vorticity distribution decreased linearly with depth, disappearing completely in the pycnocline. Statistics of divergence and strain rate generally followed the normal and χ distributions, respectively. These observations confirm a high-resolution numerical model prediction for the structure of the active submesoscale turbulence field in this area.

Shipboard ADCP and towed CTD measurements are presented of a near-inertial internal gravity wave radiating away from a zonal jet associated with the Subtropical Front in the North Pacific. Three-dimensional spatial surveys indicate persistent alternating shear layers sloping downward and equatorward from the front. As a result, depth-integrated ageostrophic shear increases sharply equatorward of the front. The layers have a vertical wavelength of about 250 m and a slope consistent with a wave of frequency 1.01 f. They extend at least 100 km south of the front. Time series confirm that the shear is associated with a downward-propagating near-inertial wave with frequency within 20% of f. A slab mixed layer model forced with shipboard and NCEP reanalysis winds suggests that wind forcing was too weak to generate the wave. Likewise, trapping of the near-inertial motions at the low-vorticity edge of the front can be ruled out because of the extension of the features well south of it. Instead, the authors suggest that the wave arises from an adjustment process of the frontal flow, which has a Rossby number about 0.2–0.3.

Interactions between vortices and a shelfbreak current are investigated, with particular attention to the exchange of waters between the continental shelf and slope. The nonlinear, three-dimensional interaction between an anticyclonic vortex and the shelfbreak current is studied in the laboratory while varying the ratio ε of the maximum azimuthal velocity in the vortex to the maximum alongshelf velocity in the shelfbreak current. Strong interactions between the shelfbreak current and the vortex are observed when ε > 1; weak interactions are found when ε < 1. When the anticyclonic vortex comes in contact with the shelfbreak front during a strong interaction, a streamer of shelf water is drawn offshore and wraps anticyclonically around the vortex. Measurements of the offshore transport and identification of the particle trajectories in the shelfbreak current drawn offshore from the vortex allow quantification of the fraction of the shelfbreak current that is deflected onto the slope; this fraction increases for increasing values of ε. Experimental results in the laboratory are strikingly similar to results obtained from observations in the Middle Atlantic Bight (MAB); after proper scaling, measurements of offshore transport and offshore displacement of shelf water for vortices in the MAB that span a range of values of ε agree well with laboratory predictions.

In late summer, satellite ocean color data consistently show localized chlorophyll blooms in the oligotrophic NE Pacific. Based on historical data and the results from recent cruises, these blooms are associated with elevated diatom abundance. However, the physical dynamics that stimulate the blooms remain unknown. Mechanisms suggested to be driving the blooms include mixing at the subtropical front, breaking of internal waves at the critical latitude, shoaling of the mixed layer depth, eddy interactions, and winter mixing of nutrients. To examine these hypotheses, we use data from four summer cruises (2002, 2007, 2008, and 2009) in this region that sampled near a bloom temporally and/or spatially. Conditions associated with five blooms (two blooms were sampled in 2009) are examined. Each area was sampled at a different stage in bloom development, including prebloom, initiation, full bloom, decline, and postbloom conditions. No one variable is found which can explain unequivocally the development of a chlorophyll bloom at a certain location. We describe a set of conditions that could result in the injection of nutrients into the surface water to stimulate a bloom. This "perfect storm" of conditions requires a subsurface stratification minimum layer that intersects the nutricline and that this minimum is close to the base of the mixed layer. These conditions are not predictable in the sense of an annual climatology; however, they do occur often enough to create a reasonably certain, if spatially variable, summer NE Pacific bloom.

Four instances of persistent intrusive deformation of the North Pacific Subtropical Front were tagged individually by a Lagrangian float and tracked for several days. Each feature was mapped in three dimensions using repeat towed observations referenced to the float. Isohaline surface deformations in the frontal zone included sheetlike folds elongated in the alongfront direction and narrow tongues extending across the front. All deformations appeared as protrusions of relatively cold, and fresh, water across the front. No corresponding features of the opposite sign or isolated lenslike structures were observed. The sheets were O(10 m) thick, protruded about 10 km into the warm saline side of the front, and were coherent for 10–30 km along the front. Having about the same thickness and cross-frontal extent as the sheets, tongues extended less than 5 km along the front.

All of the intrusions persisted as long as they were followed, several days to one week. Their structures evolved on both inertial (23 h) and subinertial (10 days) time scales in response to differential lateral advection. The water mass surrounding the intrusions participated in gradual anticyclonic rotation as a part of a mesoscale meander of the subtropical front. The intrusions may be interpreted as a manifestation of three-dimensional submesoscale turbulence of the frontal zone, driven by the mesoscale. Absence of large features of the opposite sign may be indicative of the asymmetry of the underlying dynamics.

A monthlong field survey in July 2007, focused on the North Pacific subtropical frontal zone (STFZ) near 30°N, 158°W, combined towed depth-cycling conductivity-temperature-depth (CTD) profiling with shipboard current observations. Measurements were used to investigate the distribution and structure of thermohaline intrusions. The study revealed that local extrema of vertical salinity profiles, often used as intrusion indicators, were only a subset of a wider class of distortions in thermohaline fields due to interleaving processes. A new method to investigate interleaving based on diapycnal spiciness curvature was used to describe an expanded class of laterally coherent intrusions. STFZ intrusions were characterized by their overall statistics and by a number of case studies. Thermohaline interleaving was particularly intense within 5 km of two partially compensated fronts, where intrusions with both positive and negative salinity anomalies were widespread. The vertical and cross-frontal scales of the intrusions were on the order of 10 m and 5 km, respectively. Though highly variable, the slopes of these features were typically intermediate between those of isopycnals and isohalines. Although the influence of double-diffusive processes sometime during the evolution of intrusions could not be excluded, the broad spectrum of the observed features suggests that any role of double diffusion was secondary.