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.