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.