Air–sea and air–sea–ice fluxes in the Southern Ocean play a critical role in global climate through their impact on the overturning circulation and oceanic heat and carbon uptake. The challenging conditions in the Southern Ocean have led to sparse spatial and temporal coverage of observations. This has led to a 'knowledge gap' that increases uncertainty in atmosphere and ocean dynamics and boundary-layer thermodynamic processes, impeding improvements in weather and climate models. Improvements will require both process-based research to understand the mechanisms governing air-sea exchange and a significant expansion of the observing system. This will improve flux parameterizations and reduce uncertainty associated with bulk formulae and satellite observations. Improved estimates spanning the full Southern Ocean will need to take advantage of ships, surface moorings, and the growing capabilities of autonomous platforms with robust and miniaturized sensors. A key challenge is to identify observing system sampling requirements. This requires models, Observing System Simulation Experiments (OSSEs), and assessments of the specific spatial-temporal accuracy and resolution required for priority science and assessment of observational uncertainties of the mean state and direct flux measurements. Year-round, high-quality, quasi-continuous in situ flux measurements and observations of extreme events are needed to validate, improve and characterize uncertainties in blended reanalysis products and satellite data as well as to improve parameterizations. Building a robust observing system will require community consensus on observational methodologies, observational priorities, and effective strategies for data management and discovery.

We describe a process called "squeeze dispersion" in which the squeezing of oceanic tracer gradients by waves, eddies, and bathymetric flow modulates diapycnal diffusion by centimeter to meter‐scale turbulence. Due to squeeze dispersion, the effective diapycnal diffusivity of oceanic tracers is different and typically greater than the average "local" diffusivity, especially when local diffusivity correlates with squeezing. We develop a theory to quantify the effects of squeeze dispersion on diapycnal oceanic transport, finding formulas that connect density‐averaged tracer flux, locally measured diffusivity, large‐scale oceanic strain, the thickness‐weighted average buoyancy gradient, and the effective diffusivity of oceanic tracers. We use this effective diffusivity to interpret observations of abyssal flow through the Samoan Passage reported by Alford et al. (2013) and find that squeezing modulates diapycnal tracer dispersion by factors between 0.5 and 3.

A methodology is described for control of vertically profiling floats that uses an imperfect predictive model of ocean currents. In this approach, the floats have control only over their depth. This control authority is combined with an imperfect model of ocean currents to attempt to force the floats to maintain position. First, the impact of model accuracy on the ability to station keep (e.g. maintain X–Y position) using simulated planning and nature (ground-truth in simulation) models is studied. In this study, the impact of batch versus continuous planning is examined. In batch planning the float depth plan is derived for an extended period of time and then executed open loop. In continuous planning the depth plan is updated with the actual position and the remainder of the plan re-planned based on the new information. In these simulation results are shown that (a) active control can significantly improve station keeping with even an imperfect predictive model and (b) continuous planning can mitigate the impact of model inaccuracy. Second, the effect of using heuristic path completion estimators in search are studied. In general, using a more conservative estimator increases search quality but commensurately increases the amount of search and therefore computation time. Third are presented results from an April 2015 deployment in the Pacific Ocean that show that even with an imperfect model of ocean currents, model-based control can enhance float control performance.

Mixed layers are defined to have homogeneous density, temperature, and salinity. However, bio‐optical profiles may not always be fully homogenized within the mixed layer. The relative timescales of mixing and biological processes determine whether bio‐optical gradients can form within a uniform density mixed layer. Vertical profiles of bio‐optical measurements from biogeochemical Argo floats and elephant seal tags in the Southern Ocean are used to assess biological structure in the upper ocean. Within the hydrographically defined mixed layer, the profiles show significant vertical variance in chlorophyll‐a (Chl‐a) fluorescence and particle optical backscatter. Biological structure is assessed by fitting Chl‐a fluorescence and particle backscatter profiles to functional forms (i.e., Gaussian, sigmoid, exponential, and their combinations). In the Southern Ocean, which characteristically has deep mixed layers, only 40% of nighttime bio‐optical profiles were characterized by a sigmoid, indicating a well‐mixed surface layer. Of the remaining 60% that showed structure, ∼40% had a deep fluorescence maximum below 20‐m depth that correlated with particle backscatter. Furthermore, a significant fraction of these deep fluorescence maxima were found within the mixed layer (20–80%, depending on mixed‐layer depth definition and season). Results suggest that the timescale between mixing events that homogenize the surface layer is often longer than biological timescales of restratification. We hypothesize that periods of quiescence between synoptic storms, which we estimate to be ∼3–5 days (depending on season), allow bio‐optical gradients to develop within mixed layers that remain homogeneous in density.

A three-dimensional, near real-time data-assimilative modeling system for the California coastal ocean is presented. The system consists of a Regional Ocean Modeling System (ROMS) forced by the North American Mesoscale Forecast System (NAM). The ocean model has a horizontal resolution of approximately three kilometers and utilizes a multi-scale three-dimensional variational (3DVAR) data assimilation methodology. The system is run in near real-time to produce a nowcast every six hours and a 72-hour forecast every day. The performance of this nowcast system is presented using results from a six-year period of 2009–2015.

Methods for measuring waves and winds from a Wave Glider Autonomous Surface Vehicle (ASV) are described and evaluated. The wave method utilizes the frequency spectra of orbital velocities measured by GPS, and the wind stress method utilizes the frequency spectra of turbulent wind fluctuations measured by ultrasonic anemometer. Both methods evaluate contaminations from vehicle motion. The methods were evaluated with 68 days of data over a full range of open ocean conditions, in which wave heights varied from 1 to 8 m and wind speeds varied from 1 to 17 m/s. Reference data were collected using additional sensors onboard the vehicle. For the waves method, several additional datasets are included which use independently moored Datawell waverider buoys as reference data. Bulk wave parameters are determined waverider buoys as reference data. Bulk wave parameters are determined within 5% error, with biases of less than 5%. Wind stress is determined within 4% error, with 1% bias. Wave directional spectra also compare well, although the Wave Glider results have more spread at low frequencies.

A model devised by Thorpe & Li (2014) that predicts the conditions in which stationary turbulent hydraulic jumps can occur in the flow of a continuously stratified layer over a horizontal rigid bottom is applied to, and its results compared with, observations made at several locations in the ocean. The model identifies two positions in the Samoan Passage at which hydraulic jumps should occur and where changes in the structure of the flow are indeed observed. The model predicts the amplitude of changes and the observed mode 2 form of the transitions. The predicted dissipation of turbulent kinetic energy is also consistent with observations. One location provides a particularly well-defined example of a persistent hydraulic jump. It takes the form of a 390 m thick and 3.7 km long mixing layer with frequent density inversions separated from the seabed by some 200 m of relatively rapidly moving dense water, thus revealing the previously unknown structure of an internal hydraulic jump in the deep ocean. Predictions in the Red Sea Outflow in the Gulf of Aden are relatively uncertain. Available data, and the model predictions, do not provide strong support for the existence of hydraulic jumps. In the Mediterranean Outflow, however, both model and data indicate the presence of a hydraulic jump.

The four-month mission of a Wave Glider in the Southern Ocean has demonstrated the capability for an autonomous surface vehicle to make sustained measurements of air-sea interactions in remote regions. Several new sensor payloads were integrated for this mission, including a three-axis sonic anemometer for turbulent wind stress estimation and a high-resolution atmospheric pressure gage. The mission focused on Drake Passage, where strong gradients are common along the Antarctic Circumpolar Current (ACC) fronts. Using satellite data products, pilots ashore were able to remotely navigate the Wave Glider across the ACC Polar Front and measure changes in air-sea coupling. The resulting data set combines the persistence of a mooring with the adaptability of a ship-based survey.

Lee waves are thought to play a prominent role in Southern Ocean dynamics, facilitating a transfer of energy from the jets of the Antarctic Circumpolar Current to microscale, turbulent motions important in water mass transformations. Two EM-APEX profiling floats deployed in the Drake Passage during the Diapycnal and Isopycnal Mixing Experiment (DIMES) independently measured a 120 ± 20-m vertical amplitude lee wave over the Shackleton Fracture Zone. A model for steady EM-APEX motion is developed to calculate absolute vertical water velocity, augmenting the horizontal velocity measurements made by the floats. The wave exhibits fluctuations in all three velocity components of over 15 cm s^-1 and an intrinsic frequency close to the local buoyancy frequency. The wave is observed to transport energy and horizontal momentum vertically at respective peak rates of 1.3 ± 0.2 W m^-2 and 8 ± 1 N m^-2. The rate of turbulent kinetic energy dissipation is estimated using both Thorpe scales and a method that isolates high-frequency vertical kinetic energy and is found to be enhanced within the wave to values of order 10^-7 W kg^-1. The observed vertical flux of energy is significantly larger than expected from idealized numerical simulations and also larger than observed depth-integrated dissipation rates. These results provide the first unambiguous observation of a lee wave in the Southern Ocean with simultaneous measurements of its energetics and dynamics.

Using 111 shipboard hydrographic sections across Denmark Strait occupied between 1990 and 2012, we characterize the mean conditions at the sill, quantify the water mass constituents, and describe the dominant features of the Denmark Strait Overflow Water (DSOW). The mean vertical sections of temperature, salinity, and density reveal the presence of circulation components found upstream of the sill, in particular the shelfbreak East Greenland Current (EGC) and the separated EGC. These correspond to hydrographic fronts consistent with surface-intensified southward flow. Deeper in the water column the isopycnals slope oppositely, indicative of bottom-intensified flow of DSOW. An end-member analysis indicates that the deepest part of Denmark Strait is dominated by Arctic-Origin Water with only small amounts of Atlantic-Origin Water. On the western side of the strait, the overflow water is a mixture of both constituents, with a contribution from Polar Surface Water. Weakly stratified "boluses" of dense water are present in 41% of the occupations, revealing that this is a common configuration of DSOW. The bolus water is primarily Arctic-Origin Water and constitutes the densest portion of the overflow. The boluses have become warmer and saltier over the 22 year record, which can be explained by changes in end-member properties and their relative contributions to bolus composition.

The abyssal flow of water through the Samoan Passage accounts for the majority of the bottom water renewal in the North Pacific, thereby making it an important element of the meridional overturning circulation. Here the authors report recent measurements of the flow of dense waters of Antarctic and North Atlantic origin through the Samoan Passage. A 15-month long moored time series of velocity and temperature of the abyssal flow was recorded between 2012 and 2013. This allows for an update of the only prior volume transport time series from the Samoan Passage from WOCE moored measurements between 1992 and 1994. While highly variable on multiple time scales, the overall pattern of the abyssal flow through the Samoan Passage was remarkably steady. The time-mean northward volume transport of about 5.4 Sv (1 Sv = 10⁶ m³ s⁻¹) in 2012/13 was reduced compared to 6.0 Sv measured between 1992 and 1994. This volume transport reduction is significant within 68% confidence limits (±0.4 Sv) but not at 95% confidence limits (±0.6 Sv). In agreement with recent studies of the abyssal Pacific, the bottom flow through the Samoan Passage warmed significantly on average by 1 x 10⁻³°C yr⁻¹ over the past two decades, as observed both in moored and shipboard hydrographic observations. While the warming reflects the recently observed increasing role of the deep oceans for heat uptake, decreasing flow through Samoan Passage may indicate a future weakening of this trend for the abyssal North Pacific.

A global map of open-ocean mode-1 M₂ internal tides is constructed using sea-surface height (SSH) measurements from multiple satellite altimeters during 1992–2012, representing a 20-year coherent internal tide field. A two-dimensional plane wave fit method is employed to (1) suppress mesoscale contamination by extracting internal tides with both spatial and temporal coherence, and (2) separately resolve multiple internal tidal waves. Global maps of amplitude, phase, energy and flux of mode-1 M₂ internal tides are presented. M₂ internal tides are mainly generated over topographic features including continental slopes, mid-ocean ridges and seamounts. Internal tidal beams of 100–300 km width are observed to propagate hundreds to thousands of km. Multi-wave interference of some degree is widespread, due to the M₂ internal tide's numerous generation sites and long-range propagation. The M₂ internal tide propagates across the critical latitudes for parametric subharmonic instability (28.8°S/N) with little energy loss, consistent with field measurements by MacKinnon et al. (2013). In the eastern Pacific Ocean, the M₂ internal tide loses significant energy in propagating across the Equator; in contrast, little energy loss is observed in the equatorial zones of the Atlantic, Indian, and western Pacific oceans. Global integration of the satellite observations yields a total energy of 36 PJ (1 PJ = 10¹⁵ J) for the coherent mode-1 M₂ internal tide. The satellite observed M₂ internal tides compare favorably with field mooring measurements and a global eddy-resolving numerical model.

The flow of dense water through the Samoan Passage accounts for the major part of the bottom water renewal in the North Pacific and is thus an important element of the Pacific meridional overturning circulation. A recent set of highly resolved measurements used CTD/LADCP, a microstructure profiler, and moorings to constrain the complex pathways and variability of the abyssal flow. Volume transport estimates for the dense northward current at several sections across the passage, calculated using direct velocity measurements from LADCPs, range from 3.9 x 10⁶ to 6.0 x 10⁶ ± 1 x 10⁶ m³ s^-1. The deep channel to the east and shallower pathways to the west carried about equal amounts of this volume transport, with the densest water flowing along the main eastern channel. Turbulent dissipation rates estimated from Thorpe scales and direct microstructure agree to within a factor of 2 and provide a region-averaged value of O(10^-8) W kg^-1 for layers colder than 0.8°C. Associated diapycnal diffusivities and downward turbulent heat fluxes are about 5 x 10^-3 m² s^-1 and O(10) W m^-2, respectively. However, heat budgets suggest heat fluxes 2–6 times greater. In the vicinity of one of the major sills of the passage, highly resolved Thorpe-inferred diffusivity and heat flux were over 10 times larger than the region-averaged values, suggesting the mismatch is likely due to undersampled mixing hotspots.

Wind-forced internal waves close to the inertial frequency are ubiquitous throughout the world's oceans, but observational constraints on their global energetics and impact on subsurface mixing remain scarce. This study reports on velocity measurements from three Electromagnetic Autonomous Profiling Explorers (EM-APEX) deployed in February 2009. These floats observed downward-propagating near-inertial internal waves near the Subantarctic and Polar Fronts of the Antarctic Circumpolar Current. These waves were episodic and enhanced at middepth between 500 and 1000 m. Depth-integrated kinetic energy varied between 1 and 7 kJ m^-2 and averaged 1.6 kJ m^-2 with typical group velocities of 40 m day^-1, implying an average energy flux of 3 mW m^-2 at the mixed layer base decreasing to approximately 25% of that value at 1500 m. Modeled currents forced by reanalysis winds along each float track agree with observed surface currents from EM-APEX, provided that mixed layer depth is restricted to the layer of weakest observable stratification (interpreted as the maximum depth that can remain mixed over an inertial period given the continual balance between mixing and restratification). This model estimates an average wind power of 3 mW m^-2. Shipboard wind and current observations during a strong storm show an integrated wind work of 3.5 kJ m^-2, comparable to the vertically integrated kinetic energy over the following month. Model wind work estimates are considerably less, likely because of the mixed layer depth used. A model with varying stratification in response to the wind provides a better match to the observations, emphasizing the importance of stratification within the mixed layer in amplifying wind energy input.

Observations of the semidiurnal internal tide on the California continental margin between Monterey Bay and Point Sur confirm the existence of northward energy flux predicted by numerical models of the region. Both a short-duration tide-resolving survey with expendable profilers and a multi-week timeseries from FLIP measured northward flux in the mean, supporting the hypothesis that topographic features off Point Sur are the source of the strong internal tides observed in Monterey Canyon. However, the observed depth-integrated semidiurnal flux of 450±200 W m^-1 is approximately twice as large as the most directly-comparable model and FLIP results. Though dominated by low modes with O(100 km) horizontal wavelengths, a number of properties of the semidiurnal internal tide, including kinetic and potential energy, as well as energy flux, show lateral variability on O(5 km) scales. Potential causes of this spatial variability include interference of waves from multiple sources, the sharp delineation of beams generated by abrupt topography due to limited azimuthal extent, and local generation and scattering of the internal tide into higher modes by small-scale topography. A simple two-source model of a first-mode interference pattern reproduces some of the most striking aspects of the observations.

We report the first direct turbulence observations in the Samoan Passage (SP), a 40-km wide notch in the South Pacific bathymetry through which flows most of the water supplying the North Pacific abyssal circulation. The observed turbulence is 1000 to 10,000 times typical abyssal levels — strong enough to completely mix away the densest water entering the passage — confirming inferences from previous coarser temperature and salinity sections. Accompanying towed measurements of velocity and temperature with horizontal resolution of about 250 m indicate the dominant processes responsible for the turbulence. Specifically, the flow accelerates substantially at the primary sill within the passage, reaching speeds as great as 0.55 m s^-1. A strong hydraulic response is seen, with layers first rising to clear the sill and then plunging hundreds of meters downward. Turbulence results from high shear at the interface above the densest fluid as it descends and from hydraulic jumps that form downstream of the sill. In addition to the primary sill, other locations along the multiple interconnected channels through the Samoan Passage also have an effect on the mixing of the dense water. In fact, quite different hydraulic responses and turbulence levels are observed at seafloor features separated laterally by a few kilometers, suggesting that abyssal mixing depends sensitively on bathymetric details on small scales.

Electromagnetic current instrumentation has been added to the Bathy Systems, Inc. POGO transport sondes to produce a free-falling absolute velocity profiler called EM-POGO. The POGO is a free-fall profiler that measures a depth-averaged velocity using GPS fixes at the beginning and end of a round trip to the ocean floor (or a pre-set depth). The EM-POGO adds a velocity profile determined from measurements of motionally-induced electric fields generated by the ocean current moving through the vertical component of the Earth's magnetic field. In addition to providing information about the vertical structure of the velocity, the depth-dependent measurements improve transport measurements by correcting for the non-constant fall-rate. Neglecting the variable fall rate results in errors O(1 cm s^-1). The transition from POGO to EM-POGO included electrically isolating the POGO and electric-field-measuring circuits, installing a functional GPS receiver, finding a pressure case that provided an optimal balance among crush-depth, price and size, and incorporating the electrodes, electrode collar, and the circuitry required for the electric field measurement. The first EM-POGO sea-trial was in July 1999. In August 2006 a refurbished EM-POGO collected 15 absolute velocity profiles; relative and absolute velocity uncertainty was ~ 1 cm s^-1 and 0.5–5 cm s^-1, respectively, at a vertical resolution of 25 m. Absolute velocity from the EM-POGO compared to shipboard ADCP measurements differed by ~ 1–2 cm s^-1, comparable to the uncertainty in absolute velocity from the ADCP. The EM-POGO is thus a low-cost, easy to deploy and recover, and accurate velocity profiler.

Low-mode internal tides propagate over thousands of kilometers from their generation sites, distributing tidal energy across the ocean basins. Though internal tides can have large vertical displacements (often tens of meters or more) in the ocean interior, they deflect the sea surface only by several centimeters. Because of the regularity of the tidal forcing, this small signal can be detected by state-of-the-art, repeat-track, high-precision satellite altimetry over nearly the entire world ocean. Making use of combined sea surface height measurements from multiple satellites (which together have denser ground tracks than any single mission), it is now possible to resolve the complex interference patterns created by multiple internal tides using an improved plane-wave fit technique. As examples, we present regional M2 internal tide fields around the Mariana Arc and the Hawaiian Ridge and in the North Pacific Ocean. The limitations and some perspective on the multisatellite altimetric methods are discussed.

Microstructure measurements along the axes of Monterey and Soquel Submarine Canyons reveal 200%u2013300-m-thick well-stratified turbulent near-bottom layers with average turbulent kinetic energy dissipation rates 4 x 10^-8 W kg^-1 and eddy diffusivities 16 x 10^-4 m² s^-1 (assuming mixing efficiency γ = 0.2) to at least thalweg depths of 1200 m. Turbulent dissipation rates are an order of magnitude lower in overlying waters, whereas buoyancy frequencies are only 25% higher. Well-mixed bottom boundary layer thicknesses h_N are an order of magnitude thinner than the stratified turbulent layer. Because well-stratified turbulent layers are commonly observed above slopes, arguments that mixing efficiency should be reduced on sloping boundaries do not hold in cases of energetic internal-wave generation or interaction with topography. An advective–diffusive balance is used to infer velocities and transports, predicting horizontal upslope flows of 10–50 m day^-1. Extrapolating this estimate globally suggests that canyon turbulence may contribute 2–3 times as much diapycnal transport to the World Ocean as interior mixing. The upcanyon turbulence-driven transports are not uniform, and the resulting upslope convergences will drive exchange between the turbulent layer and more quiescent interior. Predicted density surfaces of detrainment and entrainment are consistent with observed isopycnals of intermediate nepheloid and clear layers. These data demonstrate that turbulent mixing dynamics on sloping topography are fundamentally 2D or 3D in the ocean, so they cannot be accurately described by 1D models.

Internal waves contain a significant fraction of the kinetic energy in the ocean and are important intermediaries between the forcing (by wind and tide) and interior diapycnal mixing. We report here on measurements from Mindoro Strait in the Philippines (connecting the South China Sea to the Sulu Sea) of an internal wave field with a number of surprising properties that point to previously-unrecognized processes at work in the region. Continuum spectral levels are very close to typical "background" values found in the open ocean, but internal tide energy in both the diurnal and semidiurnal frequency bands is significantly elevated—and higher at the northern mooring (MP1) than the southern (MP2). Two particularly energetic depth ranges stand out at MP1: an upper layer centered near 300 m, and one at the bottom of the water column, near 1800 m. The upper layer contains both internal tides and a near-inertial band with upward and downward propagating waves and an apparent spring-neap cycle. The combination is suggestive of Parametric Subharmonic Instability as the forcing for the near-inertial band—a conclusion supported by bicoherence estimates. Mixing, estimated from density overturns, is weak over much of the water column but enhanced by about an order of magnitude in the deep layer and closely tied to the internal tide—both diurnal and semidiurnal. Near-inertial currents in this deep layer are dominantly rectilinear and not well-correlated with the mixing. Bulk mixing rates at the two sites are less than required to produce property changes seen in hydrography, suggesting additional enhancement elsewhere in the archipelago.

Satellite altimetric sea surface height anomaly (SSHA) data from Geosat Follow-on (GFO) and European Remote Sensing (ERS), as well as TOPEX/Poseidon (T/P), are merged to estimate M₂ internal tides around the Hawaiian Ridge, with higher spatial resolution than possible with single-satellite altimetry. The new estimates are compared with numerical model runs. Along-track analyses show that M₂ internal tides can be resolved from both 8 years of GFO and 15.5 years of ERS SSHA data. Comparisons at crossover points reveal that the M₂ estimates from T/P, GFO, and ERS agree well. Multisatellite altimetry improves spatial resolution due to its denser ground tracks. Thus M₂ internal tides can be plane wave fitted in 120 km x 120 km regions, compared to previous single-satellite estimates in 4° lon x 3° lat or 250 km x 250 km regions. In such small fitting regions the weaker and smaller-scale mode 2 M₂ internal tides can also be estimated.

The higher spatial resolution leads to a clearer view of the M₂ internal tide field around the Hawaiian Ridge. Discrete generation sites and internal tidal beams are clearly distinguishable, and consistent with the numerical model runs. More importantly, multisatellite altimetry produces larger M₂ internal tidal energy fluxes, which agree better with model results, than previous single-satellite estimates. This study confirms that previous altimetric underestimates are partly due to the more widely spaced ground tracks and consequently larger fitting region. Multisatellite altimetry largely overcomes this limitation.

Three autonomous profiling Electromagnetic Autonomous Profiling Explorer (EM-APEX) floats were air deployed one day in advance of the passage of Hurricane Frances (2004) as part of the Coupled Boundary Layer Air-Sea Transfer (CBLAST)-High field experiment. The floats were deliberately deployed at locations on the hurricane track, 55 km to the right of the track, and 110 km to the right of the track. These floats provided profile measurements between 30 and 200 m of in situ temperature, salinity, and horizontal velocity every half hour during the hurricane passage and for several weeks afterward. Some aspects of the observed response were similar at the three locations - the dominance of near-inertial horizontal currents and the phase of these currents - whereas other aspects were different. The largest-amplitude inertial currents were observed at the 55-km site, where SST cooled the most, by about 2.2C, as the surface mixed layer deepened by about 80 m. Based on the time-depth evolution of the Richardson number and comparisons with a numerical ocean model, it is concluded that SST cooled primarily because of shear-induced vertical mixing that served to bring deeper, cooler water into the surface layer. Surface gravity waves, estimated from the observed high-frequency velocity, reached an estimated 12-m significant wave height at the 55-km site. Along the track, there was lesser amplitude inertial motion and SST cooling, only about 1.2C, though there was greater upwelling, about 25-m amplitude, and inertial pumping, also about 25-m amplitude. Previously reported numerical simulations of the upper-ocean response are in reasonable agreement with these EM-APEX observations provided that a high wind speed-saturated drag coefficient is used to estimate the wind stress. A direct inference of the drag coefficient CD is drawn from the momentum budget. For wind speeds of 32-47 m s^-1, CD ~ 1.4 x 10^-3.

This article discusses the challenges of developing a regional ocean prediction model for the Philippine Archipelago, a complex area in terms of geometry, bathymetry-dominated dynamics and variability, and strong local and remote wind forcing, where there are limited temporal and spatial ocean measurements. We used the Regional Ocean Modeling System (ROMS) for real-time forecasting during the Philippine Straits Dynamics Experiment (2007-2009) observational program. The article focuses on the prediction experiments before and during the exploratory cruise period, June 6 - July 3, 2007. The gathered observations were not available in real time, so the 4-Dimensional Variational (4D-Var) data assimilation experiments were carried out in hindcast mode. The best estimate of ocean state (nowcast) is determined by combining satellite-derived products for sea surface temperature and height, and subsurface temperature and salinity measurements from several hydrographic assets over a sequential five-day data assimilation window. The largest source of forecast uncertainty is from the prescribed lateral boundary conditions in the nearby Pacific Ocean, especially excessive salt flux. This result suggests that remote forcing and inflows from the Pacific are crucial for predicting ocean circulation in the Philippine Archipelago region. The lateral boundary conditions are derived from 1/12 degree global HYbrid Coordinate Ocean Model (HYCOM) daily snapshots. The incremental, strong-constraint 4D-Var data assimilation successfully decreased temperature and salinity errors of the real-time, nonassimilative control forecast by 38% and 49%, respectively.

Internal wave measurements from moorings and profiling floats throughout the Philippine Archipelago, collected as part of the Office of Naval Research Philippine Straits Dynamics Experiment, reveal a wealth of subsurface processes, some of which have not been observed previously (in the Philippines or elsewhere). Complex bathymetry and spatially varying tide and wind forcing produce distinct internal wave environments within the network of seas and channels, ranging from quiescent interior basins to remotely forced straits. Internal tides in both the diurnal and semidiurnal bands dominate much of the velocity structure and are likely the dominant source of energy for mixing in the region. In addition, the transfer of energy from the internal tide directly to near-inertial motions through parametric subharmonic instability appears to be active and, rather than wind forcing, is the dominant source of near-inertial band energy.

The vertical dispersion of a tracer released on a density surface near 1500-m depth in the Antarctic Circumpolar Current west of Drake Passage indicates that the diapycnal diffusivity, averaged over 1 yr and over tens of thousands of square kilometers, is (1.3 ± 0.2) x 10^-5 m² s^-1. Diapycnal diffusivity estimated from turbulent kinetic energy dissipation measurements about the area occupied by the tracer in austral summer 2010 was somewhat less, but still within a factor of 2, at (0.75 ± 0.07) x 10^-5 m² s^-1. Turbulent diapycnal mixing of this intensity is characteristic of the midlatitude ocean interior, where the energy for mixing is believed to derive from internal wave breaking. Indeed, despite the frequent and intense atmospheric forcing experienced by the Southern Ocean, the amplitude of finescale velocity shear sampled about the tracer was similar to background amplitudes in the midlatitude ocean, with levels elevated to only 20%–50% above the Garrett–Munk reference spectrum. These results add to a long line of evidence that diapycnal mixing in the interior middepth ocean is weak and is likely too small to dictate the middepth meridional overturning circulation of the ocean.

An autonomous, profiling float called EM-APEX was developed to provide a quantitative and comprehensive description of the ocean side of hurricane-ocean interaction. EM-APEX measures temperature, salinity and pressure to CTD quality and relative horizontal velocity with an electric field sensor. Three prototype floats were air-deployed into the upper ocean ahead of Hurricane Frances (2004). All worked properly and returned a highly resolved description of the upper ocean response to a category 4 hurricane. At a float launched 55 km to the right of the track, the hurricane generated large amplitude, inertially rotating velocity in the upper 120 m of the water column. Coincident with the hurricane passage there was intense vertical mixing that cooled the near surface layer by about 2.2°C. We find consistent model simulations of this event provided the wind stress is computed from the observed winds using a high wind-speed saturated drag coefficient.

The overflow of dense water from the Nordic Seas through the Faroe Bank Channel (FBC) has attributes suggesting hydraulic control—primarily an asymmetry across the sill reminiscent of flow over a dam. However, this aspect has never been confirmed by any quantitative measure, nor is the position of the control section known. This paper presents a comparison of several different techniques for assessing the hydraulic criticality of oceanic overflows applied to data from a set of velocity and hydrographic sections across the FBC. These include 1) the cross-stream variation in the local Froude number, including a modified form that accounts for stratification and vertical shear, 2) rotating hydraulic solutions using a constant potential vorticity layer in a channel of parabolic cross section, and 3) direct computation of shallow water wave speeds from the observed overflow structure. Though differences exist, the three methods give similar answers, suggesting that the FBC is indeed controlled, with a critical section located 20–90 km downstream of the sill crest. Evidence of an upstream control with respect to a potential vorticity wave is also presented. The implications of these results for hydraulic predictions of overflow transport and variability are discussed.

An integrated analysis of turbulence observations from four unique instrument platforms obtained over the Hawaiian Ridge leads to an assessment of the vertical, cross-ridge, and along-ridge structure of turbulence dissipation rate and diffusivity. The diffusivity near the seafloor was, on average, 15 times that in the midwater column. At 1000-m depth, the diffusivity atop the ridge was 30 times that 10 km off the ridge, decreasing to background oceanic values by 60 km. A weak (factor of 2) spring–neap variation in dissipation was observed. The observations also suggest a kinematic relationship between the energy in the semidiurnal internal tide (E) and the depth-integrated dissipation (D), such that D ~ E^1±0.5 at sites along the ridge. This kinematic relationship is supported by combining a simple knife-edge model to estimate internal tide generation, with wave–wave interaction time scales to estimate dissipation. The along-ridge kinematic relationship and the observed vertical and cross-ridge structures are used to extrapolate the relatively sparse observations along the length of the ridge, giving an estimate of 3 ± 1.5 GW of tidal energy lost to turbulence dissipation within 60 km of the ridge. This is roughly 15% of the energy estimated to be lost from the barotropic tide.

This study examines water property distributions in the deep South China Sea and adjoining Pacific Ocean using all available hydrographic data. Our analysis reveals that below about 1500 m there is a persistent baroclinic pressure gradient driving flow from the Pacific into the South China Sea through Luzon Strait. Applying hydraulic theory with assumptions of zero potential vorticity and flat bottom to the Luzon Strait yields a transport estimate of 2.5 Sv (1 Sv = 10⁶ m³ s^-1). Some implications of this result include: (1) a residence time of less than 30 years in the deep South China Sea, (2) a mean diapycnal diffusivity as large as 10^-3 m² s^-1, and (3) an abyssal upwelling rate of about 3 x 10^-6 m s^-1. These quantities are consistent with residence times based on oxygen consumption rates. The fact that all of the inflowing water must warm up before leaving the basin implies that this marginal sea contributes to the water mass transformations that drive the meridional overturning circulation in the North Pacific. Density distributions within the South China Sea basin suggest a cyclonic deep boundary current system, as might be expected for an overflow-driven abyssal circulation.

A series of idealised numerical simulations of dense water flowing down a broad uniform slope are presented, employing both a z-coordinate model (the MIT general circulation model) and an isopycnal coordinate model (the Hallberg Isopycnal Model). Calculations are carried out at several different horizontal and vertical resolutions, and for a range of physical parameters. A subset of calculations are carried out at very high resolution using the non-hydrostatic variant of the MITgcm. In all calculations dense water descends the slope while entraining and mixing with ambient fluid. The dependence of entrainment, mixing and down-slope descent on resolution and vertical coordinate are assessed. At very coarse resolutions the z-coordinate model generates excessive spurious mixing, and dense water has difficulty descending the slope. However, at intermediate resolutions the mixing in the z-coordinate model is less than found in the high-resolution non-hydrostatic simulations, and dense water descends further down the slope. Isopycnal calculations show less resolution dependence, although entrainment and mixing are both reduced slightly at coarser resolution. At intermediate resolutions the z-coordinate and isopycnal models produce similar levels of mixing and entrainment. These results provide a benchmark against which future developments in overflow entrainment parameterizations in both z-coordinate and isopycnal models may be compared.

We developed an autonomous ocean profiling velocity and density float that provides exceptional vertical coverage and temporal resolution to depths of 2000 m for deployments of many years. Electrodes were added to the exterior of standard WRC APEX floats, and electronics were added inside. The electrode voltages result from the motion of seawater and the instrument through the Earth's magnetic field. Other systems included magnetic compass, tilt, CTD, GPS, and Iridium (providing sampling/mission changes). Three EM-APEX floats were deployed from a C-130 aircraft ahead of Hurricane Frances. The floats profiled for 10 hr from the surface to 200 m, then continued profiling between 30 m and 200 m with excursions to 500 m every half inertial period. The velocity computations were performed onboard and saved for later transmission. After five days, the floats surfaced and transmitted the accumulated processed observations, then the floats profiled from 500 m every half inertial period until recovered early in October located by GPS and Iridium.

Bulk properties of the Denmark Strait overflow (DSO) plume observed in velocity and hydrography surveys undertaken in 1997 and 1998 are described. Despite the presence of considerable short-term variability, it is found that the pathway and evolution of the plume density anomaly are remarkably steady. Bottom stress measurements show that the pathway of the plume core matches well with a rate of descent controlled by friction. The estimated entrainment rate diagnosed from the rate of plume dilution with distance shows a marked increase in entrainment at approximately 125 km from the sill, leading to a net dilution consistent with previous reports of a doubling of overflow transport measured by current meter arrays. The entrainment rate increase is likely related to the increased topographic slopes in the region, compounded by a decrease in interface stratification as the plume is diluted and enters a denser background.

We report on a combined modeling and observational effort to understand the Denmark Strait Overflow (DSO). Four cruises over the course of 3 years mapped hydrographic properties and velocity fields with high spatial resolution. The observations reveal the mean path of the dense water, as well as the presence of strong barotropic flows, energetic variability, and strong bottom friction and entrainment. A regional sigma coordinate numerical model of interbasin exchange using realistic bottom topography and an overflow forced only by an upstream reservoir of dense fluid is compared with the observations and used to further investigate these processes. The model successfully reproduces the volume transport of dense water at the sill, as well as the 1000-m descent of the dense water in the first 200 km from the sill and the intense eddies generated at 1–3 day intervals. Hydraulic control of the mean flow is indicated by a region supercritical to long gravity waves in the dense layer located approximately 100 km downstream of the sill in both model and observations. In addition, despite the differences in surface forcing, both model and observations exhibit similar transitions from mostly barotropic flow at the sill to a bottom-trapped baroclinic flow downstream, indicating the dominant role of the overflow in determining the full water column dynamics.

We report on a rapid high-resolution survey of the Denmark Strait overflow (DSO) as it crosses the sill, the first such program to incorporate full-water-column velocity profiles in addition to conventional hydrographic measurements. Seven transects with expendable profilers over the course of one week are used to estimate volume transport as a function of density. Our observations reveal the presence of a strongly barotropic flow associated with the nearly-vertical front dividing the Arctic and Atlantic waters. The seven-section mean transport of water denser than σ_Θ=27.8 is 2.7±0.6Sv, while the mean transport of water colder than 2.0°C is 3.8±0.8 Sv. Although this is larger than the 2.9 Sv of Θ < 2°C water measured by a 1973 current meter array, we find that a sampling of our sections equivalent to the extent of that array also measures 2.9Sv of cold water. Both the structure and magnitude of the measured flow are reproduced well by a high-resolution numerical model of buoyancy-driven exchange with realistic topography.