Search

Moored observations and a realistic, tidally forced 3D model are presented of flow and internal-tide-driven turbulence over a supercritical 3D fan in southeastern Luzon Strait. Two stacked moored profilers, an acoustic Doppler current profiler, and a thermistor string measured horizontal velocity, density, and salinity over nearly the entire water column every 1.5 h for 50 days. Observed dissipation rate computed from Thorpe scales decays away from the bottom and shows a strong spring–neap cycle; observed depth-integrated dissipation rate scales as where UBT is the barotropic velocity. Vertical velocities are strong enough to be comparable at times to the vertical profiling speed of the moored profilers, requiring careful treatment to quantify bias in dissipation rate estimates. Observations and the model are in reasonable agreement for velocity, internal wave displacement and depth-integrated dissipation rate, allowing the model to be used to understand the 3D flow. Turbulence is maximum following the transition from up-fan to down-fan flow, consistent with breaking lee waves advected past the mooring as seen previously at the Hawaiian Ridge, but asymmetric flow arises because of the 3D topography. Observed turbulence varies by a factor of 2 over the four observed spring tides as low-frequency near-bottom flow changes, but the exact means for inclusion of such low-frequency effects is not clear. Our results suggest that for the extremely energetic turbulence associated with breaking lee waves, dissipation rates may be quantitatively predicted to within a factor of 2 or so using numerical models and simple scalings. Significance Statement This paper describes deep ocean turbulence caused by strong tidal and low-frequency meandering flows over and around a three-dimensional bump, using moored observations and a computer simulation. Such information is important for accurately including these effects in climate simulations. The observations and model agree well enough to be able to use both to synthesize a coherent picture. The observed and modeled turbulence scale as the cube of the tidal speed as expected from theory, but low-frequency flows complicate the picture. We also demonstrate the underestimation of the turbulence that can result when vertical profiling rates are comparable to the internal wave velocities.

Turbulent mixing in the ocean is often parameterized in terms of the downscale energy transfer by internal waves. Expressed in terms of the vertical wavenumber spectrum of velocity shear and isopycnal strain, the “finescale parameterization” relies on several parameters that set the shape of the spectrum. Here, we use an unsupervised learning model to identify dominant spectral shapes and their geographical dependencies in data from 1875 global hydrographic profiles. The model reveals clustered regions with underlying spectra deviating significantly from the canonical Garrett-Munk (GM) spectrum, suggesting potential departures from the underlying assumption about the downscale energy cascade. These clusters are further found to be associated with regions where parameterized estimates of mixing differ significantly from temperature microstructure observations. Extending the methodology to other hydrographic datasets has the potential to reveal reasons for parameterization bias and to identify the dynamical underpinnings leading to more robust parameterizations of oceanic turbulent mixing.

Mode-1 internal tides can propagate far away from their generation sites, but how and where their energy is dissipated is not well understood. One example is the semidiurnal internal tide generated south of New Zealand, which propagates over a thousand kilometers before impinging on the continental slope of Tasmania. In situ observations and model results from a recent process-study experiment are used to characterize the spatial and temporal variability of the internal tide on the southeastern Tasman slope, where previous studies have quantified large reflectivity. As expected, a standing wave pattern broadly explains the cross-slope and vertical structure of the observed internal tide. However, model and observations highlight several additional features of the internal tide on the continental slope. The standing wave pattern on the sloping bottom as well as small-scale bathymetric corrugations lead to bottom-enhanced tidal energy. Over the corrugations, larger tidal currents and isopycnal displacements are observed along the trough as opposed to the crest. Despite the long-range propagation of the internal tide, most of the variability in energy density on the slope is accounted by the spring–neap cycle. However, the timing of the semidiurnal spring tides is not consistent with a single remote wave and is instead explained by the complex interference between remote and local tides on the Tasman slope. These observations suggest that identifying the multiple waves in an interference pattern and their interaction with small-scale topography is an important step in modeling internal energy and dissipation.

The scale-dependent variance of tracer properties in the ocean bears the imprint of the oceanic eddy field. Anomalies in spice (which combines anomalies in temperature T and salinity S on isopycnal surfaces) act as passive tracers beneath the surface mixed layer (ML). We present an analysis of spice distributions along isopycnals in the upper 200 m of the ocean, calculated with over 9000 vertical profiles of T and S measured along ~4800 km of ship tracks in the Bay of Bengal. The data are from three separate research cruises—in the winter monsoon season of 2013 and in the late and early summer monsoon seasons of 2015 and 2018. We present a spectral analysis of horizontal tracer variance statistics on scales ranging from the submesoscale (~1 km) to the mesoscale (~100 km). Isopycnal layers that are closer to the ML-base exhibit redder spectra of tracer variance at scales km than is predicted by theories of quasigeostrophic turbulence or frontogenesis. Two plausible explanations are postulated. The first is that stirring by submesoscale motions and shear dispersion by near-inertial waves enhance effective horizontal mixing and deplete tracer variance at horizontal scales km in this region. The second is that the spice anomalies are coherent with dynamical properties such as potential vorticity, and not interpretable as passively stirred.

An along-channel momentum budget is quantified in the near-field plume region of a controlled river flow entering Doubtful Sound, New Zealand. Observations include highly resolved density, velocity and turbulence, enabling a momentum budget to be constructed over a control volume. Estimates of internal stress (τ) were made from direct measurements of turbulence dissipation rates (ϵ) using vertical microstructure profiles. High flow speeds of the surface plume over 2 m s−1 and strong stratification (N2∼10-1 s−2) resulted in enhanced turbulence dissipation rates (ϵ>10-3 W kg−1) and internal stress (τ>10-2 m2 s−2) at the base of the surface layer. Internal waves were observed propagating along the base of the plume, potentially released subsequent to a hydraulic jump in the initial 1 km downstream of the plume discharge point. The momentum flux divergence of these internal waves suggests that almost 15 % of the total plume momentum can be transported out of the system by wave radiation, therefore playing a crucial role in the redistribution of momentum within the near-field plume. Observations illustrate that the evolution of the momentum budget components vary between the distinct surface plume layer and the turbulent, shear-stratified interfacial layer. Within the surface plume, a momentum balance was achieved. The dynamical balance demonstrates that the deceleration of the plume, driven by along-channel advection, is controlled by turbulence stress from the plume discharge point to as far as 3 km downstream. In the interfacial layer, however, the momentum equation was dominated by the turbulence stress term and the balance was not closed. The redistribution of momentum within the shear-stratified layer by internal wave radiation and other hydraulic features could account for the discrepancy in the budget.

Ice mélange has been postulated to impact glacier and fjord dynamics through a variety of mechanical and thermodynamic couplings. However, observations of these interactions are very limited. Here, we report on glaciological and oceanographic data that were collected from 2016 to 2017 at LeConte Glacier and Bay, Alaska, and serendipitously captured the formation, flow and break-up of ephemeral ice mélange. Sea ice formed overnight in mid-February. Over the subsequent week, the sea ice and icebergs were compacted by the advancing glacier terminus, after which the ice mélange flowed quasi-statically. The presence of ice mélange coincided with the lowest glacier velocities and frontal ablation rates in our record. In early April, increasing glacier runoff and the formation of a sub-ice-mélange plume began to melt and pull apart the ice mélange. The plume, outgoing tides and large calving events contributed to its break-up, which took place over a week and occurred in pulses. Unlike observations from elsewhere, the loss of ice mélange integrity did not coincide with the onset of seasonal glacier retreat. Our observations provide a challenge to ice mélange models aimed at quantifying the mechanical and thermodynamic couplings between ice mélange, glaciers and fjords.

Towed shipboard and moored observations show internal gravity waves over a tall, supercritical submarine ridge that reaches to 1000 m below the ocean surface in the tropical western Pacific north of Palau. The lee-wave or topographic Froude number, Nh0/U0 (where N is the buoyancy frequency, h0 the ridge height, and U0 the farfield velocity), ranged between 25 and 140. The waves were generated by a superposition of tidal and low-frequency flows and thus had two distinct energy sources with combined amplitudes of up to 0.2 m s−1. Local breaking of the waves led to enhanced rates of dissipation of turbulent kinetic energy reaching above 10−6 W kg−1 in the lee of the ridge near topography. Turbulence observations showed a stark contrast between conditions at spring and neap tide. During spring tide, when the tidal flow dominated, turbulence was approximately equally distributed around both sides of the ridge. During neap tide, when the mean flow dominated over tidal oscillations, turbulence was mostly observed on the downstream side of the ridge relative to the mean flow. The drag exerted by the ridge on the flow, estimated to for individual ridge crossings, and the associated power loss, thus provide an energy sink both for the low-frequency ocean circulation and the tidal flow.

Despite sufficient wind forcing, internal waves in the South China Sea do not exhibit the strong near-inertial wave (NIW) peak that is typical in most of the world oceans. Using data from 10 contemporaneous moorings deployed in summer 2011, we show that strong isopycnal vertical tidal displacements transfer most of the near-inertial (NI) kinetic energy (KE) to frequencies higher than the inertial frequency in an Eulerian reference frame. Transforming to an isopycnal-following reference frame increases the KE at NI frequencies, suggesting the presence of NIWs. However, the projection onto a semi-Lagrangian coordinate system still underestimates the expected NI peak. To fully resolve NIWs requires the use of time-dependent vertical wavenumber–frequency spectra because the intrinsic frequency of the NIWs varies substantially, owing to Doppler shifting by lateral mesoscale flows. Here, we show NIW intrinsic frequency variations of ±0.2 cpd within few days, of similar magnitude as the observed variations of relative vorticity associated with the meandering Kuroshio.

This chapter gives an overview of recent advances in in situ observations of ocean mixing. It starts by providing a brief history of measuring ocean mixing. It then describes adaptations of traditional measurements and methods to autonomous platforms, including profiling floats and autonomous underwater vehicles. For each approach, the method is described as well as its limitations. Evidence for successful application of each method is provided.

The abyssal Southwest Pacific Basin has warmed significantly between 1992-2017, consistent with warming along the bottom limb of the meridional overturning circulation seen throughout the global oceans. Here we present a framework for assessing the abyssal heat budget that includes the time-dependent unsteady effects of decadal warming and direct and indirect estimates of diapycnal mixing from microscale temperature measurements and finescale parameterizations. The unsteady terms estimated from the decadalwarming rate are shown to be within a factor of 3 of the steady state terms in the abyssal heat budget for the coldest portion of the water column and therefore, cannot be ignored. We show that a reduction in the lateral heat flux for the coldest temperature classes compensated by an increase in warmer waters advected into the basin has important implications for the heat balance and diffusive heat fluxes in the basin. Finally, vertical diffusive heat fluxes are estimated in different ways: using the newly available CTD-mounted microscale temperature measurements, a finescale strain parameterization, and a vertical kinetic energy parameterization from data along the P06 transect along 32.5°S. The unsteady-state abyssal heat budget for the basin shows closure within error estimates, demonstrating that (i) unsteady terms have become consequential for the heat balance in the isotherms closest to the ocean bottom and (ii) direct and indirect estimates from full depth GO-SHIP hydrographic transects averaged over similarly large spatial and temporal scales can capture the basin-averaged abyssal mixing needed to close the deep overturning circulation.

An edited version of this paper was published by AGU. Copyright 2019 American Geophysical Union.

We present a workflow to track icebergs in proglacial fjords using oblique time-lapse photos and the Lucas-Kanade optical flow algorithm. We employ the workflow at LeConte Bay, Alaska, where we ran five time-lapse cameras between April 2016 and September 2017, capturing more than 400 000 photos at frame rates of 0.5–4.0 min−1. Hourly to daily average velocity fields in map coordinates illustrate dynamic currents in the bay, with dominant downfjord velocities (exceeding 0.5 m s−1 intermittently) and several eddies. Comparisons with simultaneous Acoustic Doppler Current Profiler (ADCP) measurements yield best agreement for the uppermost ADCP levels (~ 12 m and above), in line with prevalent small icebergs that trace near-surface currents. Tracking results from multiple cameras compare favorably, although cameras with lower frame rates (0.5 min−1) tend to underestimate high flow speeds. Tests to determine requisite temporal and spatial image resolution confirm the importance of high image frame rates, while spatial resolution is of secondary importance. Application of our procedure to other fjords will be successful if iceberg concentrations are high enough and if the camera frame rates are sufficiently rapid (at least 1 min−1 for conditions similar to LeConte Bay).

Low-mode internal tides, a dominant part of the internal wave spectrum, carry energy over large distances, yet the ultimate fate of this energy is unknown. Internal tides in the Tasman Sea are generated at Macquarie Ridge, south of New Zealand, and propagate northwest as a focused beam before impinging on the Tasmanian continental slope. In situ observations from the Tasman Sea capture synoptic measurements of the incident semidiurnal mode-1 internal-tide, which has an observed wavelength of 183 km and surface displacement of approximately 1 cm. Plane-wave fits to in situ and altimetric estimates of surface displacement agree to within a measurement uncertainty of 0.3 cm, which is the same order of magnitude as the nonstationary (not phase locked) mode-1 tide observed over a 40-day mooring deployment. Stationary energy flux, estimated from a plane-wave fit to the in situ observations, is directed toward Tasmania with a magnitude of 3.4 ± 1.4 kW m−1, consistent with a satellite estimate of 3.9 ± 2.2 kW m−1. Approximately 90% of the time-mean energy flux is due to the stationary tide. However, nonstationary velocity and pressure, which are typically 1/4 the amplitude of the stationary components, sometimes lead to instantaneous energy fluxes that are double or half of the stationary energy flux, overwhelming any spring–neap variability. Despite strong winds and intermittent near-inertial currents, the parameterized turbulent-kinetic-energy dissipation rate is small (i.e., 10−10 W kg−1) below the near surface and observations of mode-1 internal tide energy-flux convergence are indistinguishable from zero (i.e., the confidence intervals include zero), indicating little decay of the mode-1 internal tide within the Tasman Sea.

Lateral submesoscale processes and their influence on vertical stratification at shallow salinity fronts in the central Bay of Bengal during the winter monsoon are explored using high-resolution data from a cruise in November 2013. The observations are from a radiator survey centered at a salinity-controlled density front, embedded in a zone of moderate mesoscale strain (0.15 times the Coriolis parameter) and forced by winds with a downfront orientation. Below a thin mixed layer, often ≤10 m, the analysis shows several dynamical signatures indicative of submesoscale processes: (i) negative Ertel potential vorticity (PV); (ii) low-PV anomalies with O(1–10) km lateral extent, where the vorticity estimated on isopycnals and the isopycnal thickness are tightly coupled, varying in lockstep to yield low PV; (iii) flow conditions susceptible to forced symmetric instability (FSI) or bearing the imprint of earlier FSI events; (iv) negative lateral gradients in the absolute momentum field (inertial instability); and (v) strong contribution from differential sheared advection at O(1) km scales to the growth rate of the depth-averaged stratification. The findings here show one-dimensional vertical processes alone cannot explain the vertical stratification and its lateral variability over O(1–10) km scales at the radiator survey.

The classic model of fjord renewal is complicated by tidewater glacier fjords, where submarine melt and subglacial discharge provide substantial buoyancy forcing at depth. Here we use a suite of idealized, high-resolution numerical ocean simulations to investigate how fjord circulation driven by subglacial plumes, tides, and wind stress depends on fjord width, grounding line depth, and sill height. We find that the depth of the grounding line compared to the sill is a primary control on plume-driven renewal of basin waters. In wide fjords the plume exhibits strong lateral recirculation, increasing the dilution and residence time of glacially-modified waters. Rapid drawdown of basin waters by the subglacial plume in narrow fjords allows for shelf waters to cascade deep into the basin; wide fjords result in a thin, boundary current of shelf waters that flow toward the terminus slightly below sill depth. Wind forcing amplifies the plume-driven exchange flow; however, wind-induced vertical mixing is limited to near-surface waters. Tidal mixing over the sill increases in-fjord transport of deep shelf waters and erodes basin stratification above the sill depth. These results underscore the first-order importances of fjord-glacier geometry in controlling circulation in tidewater glacier fjords and, thus, ocean heat transport to the ice.

Shipboard measurements of velocity and density were obtained in the vicinity of a small channel in the Mendocino Ridge, where flows were predominantly tidal. Measured daily inequalities in transport are much greater than those predicted by a barotropic tide model, with the strongest transport associated with full depth flows and the weakest with shallow, surface-confined flows. A regional numerical model of the area finds that the subinertial K1 (diurnal) tidal constituent generates topographically trapped waves that propagate anticyclonically around the ridge and are associated with enhanced near-topographic K1 transports. The interaction of the baroclinic trapped waves with the surface tide produces a tidal flow whose northward transports alternate between being surface confined and full depth. Full depth flows are associated with the generation of a large-amplitude tidal lee wave on the northward face of the ridge, while surface-confined flows are largely nonturbulent. The regional model demonstrates that, consistent with field observations, near-topographic dissipation over the entire ridge is diurnally modulated, despite the semidiurnal tidal constituent having larger barotropic velocities. It is concluded that at this location it is the bottom-trapped subinertial internal tide that governs near-topographic dissipation and mixing. The effect of the trapped wave on regional energetics is to increase the fraction of converted barotropic–baroclinic tidal energy that dissipates locally.

At tidewater glaciers, plume dynamics affect submarine melting, fjord circulation, and the mixing of meltwater. Models often rely on buoyant plume theory to parameterize plumes and submarine melting; however, these parameterizations are largely untested due to a dearth of near-glacier measurements. Here, we present a high-resolution ocean survey by ship and remotely-operated boat near the terminus of Kangerlussuup Sermia in west Greenland. These novel observations reveal the 3D structure and transport of a near-surface plume, originating at a large undercut conduit in the glacier terminus, that is inconsistent with axisymmetric plume theory, the most common representation of plumes in ocean-glacier models. Instead, the observations suggest a wider upwelling plume – a ‘truncated’ line plume of ∼200 m width – with higher entrainment and plume-driven melt compared to the typical axisymmetric representation. Our results highlight the importance of a subglacial outlet's geometry in controlling plume dynamics, with implications for parameterizing the exchange flow and submarine melt in glacial fjord models.

Observations from Eel Canyon, located on the north coast of California, show that elevated turbulence in the full water column arises from the convergence of remotely generated internal wave energy. The incoming semidiurnal and bottom-trapped diurnal internal tides generate complex interference patterns. The semidiurnal internal tide sets up a partly standing wave within the canyon due to reflection at the canyon head, dissipating all of its energy within the canyon. Dissipation in the near bottom is associated with the diurnal trapped tide, while midwater isopycnal shear and strain is associated with the semidiurnal tide. Dissipation is elevated up to 600 m off the bottom, in contrast to observations over the flat continental shelf where dissipation occurs closer to the topography. Slope canyons are sinks for internal wave energy and may have important influences on the global distribution of tidally driven mixing.

An along-channel momentum budget is quantified in the near-field plume region of a controlled river flow entering Doubtful Sound, New Zealand. Observations include highly resolved density, velocity and turbulence, enabling a momentum budget to be constructed over a control volume. Estimates of internal stress (τ) were made from direct measurements of turbulence dissipation rates (ε) using vertical microstructure profiles. High flow speeds of the surface plume over 2 m s−1 and strong stratification (N2 ~ 10−1 s−2) resulted in enhanced turbulence dissipation rates (ε > 10−3 W kg−1) and internal stress (τ > 10−2 m2 s−2) at the base of surface layer. An observed transition from a supercritical to sub-critical flow regime in the initial 1 km indicates the presence of an internal hydraulic jump and the subsequent release of internal gravity waves. The momentum flux divergence of these internal waves suggests that almost 15 % of the total plume momentum can be transported out of the system by wave radiation, therefore playing a crucial role in the redistribution of momentum within the near-field plume. Observations illustrate that the evolution of the momentum budget components vary between the distinct surface plume layer and the turbulent, shear-stratified interfacial layer. Within the surface plume, a momentum balance was achieved. The dynamical balance demonstrates that the deceleration of the plume, driven by along-channel advection, is controlled by turbulence stress from the plume discharge point to as far as 3 km downstream. In the interfacial layer however, the momentum equation was dominated by the turbulence stress term and the balance was not closed. The redistribution of momentum within the shear-stratified layer by the observed hydraulic jump and internal wave radiation could account for the discrepancy in the budget.

The Argo Program has been implemented and sustained for almost two decades, as a global array of about 4000 profiling floats. Argo provides continuous observations of ocean temperature and salinity versus pressure, from the sea surface to 2000 dbar. The successful installation of the Argo array and its innovative data management system arose opportunistically from the combination of great scientific need and technological innovation. Through the data system, Argo provides fundamental physical observations with broad societally-valuable applications, built on the cost-efficient and robust technologies of autonomous profiling floats. Following recent advances in platform and sensor technologies, even greater opportunity exists now than 20 years ago to (i) improve Argo’s global coverage and value beyond the original design, (ii) extend Argo to span the full ocean depth, (iii) add biogeochemical sensors for improved understanding of oceanic cycles of carbon, nutrients, and ecosystems, and (iv) consider experimental sensors that might be included in the future, for example to document the spatial and temporal patterns of ocean mixing. For Core Argo and each of these enhancements, the past, present, and future progression along a path from experimental deployments to regional pilot arrays to global implementation is described. The objective is to create a fully global, top-to-bottom, dynamically complete, and multidisciplinary Argo Program that will integrate seamlessly with satellite and with other in situ elements of the Global Ocean Observing System (Legler et al., 2015). The integrated system will deliver operational reanalysis and forecasting capability, and assessment of the state and variability of the climate system with respect to physical, biogeochemical, and ecosystems parameters. It will enable basic research of unprecedented breadth and magnitude, and a wealth of ocean-education and outreach opportunities., Sí

Underwater melting How fast does warm ocean water melt glaciers that terminate in the sea? That question is central to understanding how fast ice sheets may lose mass, and thus how fast sea level will rise, in response to global warming, but there are few data about the process. Sutherland et al. used repeat multibeam sonar surveys to observe an Alaskan subsurface tidewater glacier face to create a time series of its melting and calving patterns. They observed melt rates up to a hundred times larger than those predicted by theory, observations that compel us to reevaluate predictions of such ice loss. Science , this issue p. 369

Turbulence in geophysical flows tends to organize itself so that the mean flow remains close to a stability boundary in parameter space. That characteristic suggests self-organized criticality (SOC), a statistical property that has been identified in a range of complex phenomena including earthquakes, forest fires and solar flares. This note explores the relationship between the properties of forced, sheared, stratified turbulence (as found in oceans, atmospheres and other geophysical fluids) and those of SOC. Self-organization to the critical state is demonstrated in a wide range of cases drawn mostly (but not entirely) from in situ observations of ocean turbulence. Turbulent events in the ocean also exhibit a second characteristic associated with SOC: their sizes follow a power-law distribution indicating self-similarity. These results suggest SOC as a new conceptual foundation for the study of geophysical turbulence, an explanation for the mixing efficiency of ocean turbulence and a potential for cross-fertilization with other areas of geophysics.

The reflection of a low-mode internal tide on the Tasman continental slope is investigated using simulations of realistic and simplified topographies. The slope is supercritical to the internal tide, which should predict a large fraction of the energy reflected. However, the response to the slope is complicated by a number of factors: the incoming beam is confined laterally, it impacts the slope at an angle, there is a roughly cylindrical rise directly offshore of the slope, and a leaky slope-mode wave is excited. These effects are isolated in simulations that simplify the topography. To separate the incident from the reflected signal, a response without the reflector is subtracted from the total response to arrive at a reflected signal. The real slope reflects approximately 65% of the mode-1 internal tide as mode 1, less than two-dimensional linear calculations predict, because of the three-dimensional concavity of the topography. It is also less than recent glider estimates, likely as a result of along-slope inhomogeneity. The inhomogeneity of the response comes from the Tasman Rise that diffracts the incoming tidal beam into two beams: one focused along beam and one diffracted to the north. Along-slope inhomogeneity is enhanced by a partially trapped, superinertial slope wave that propagates along the continental slope, locally removing energy from the deep-water internal tide and reradiating it into the deep water farther north. This wave is present even in a simplified, straight slope topography; its character can be predicted from linear resonance theory, and it represents up to 30% of the local energy budget.