Description: Cruise data collected across three separate cruises MeRMEED-1 (WS16336; 1-7 December 2016), MeRMEED-2 (WS17305; 31 October - 10 November 2017) and MeRMEED-3 (WS18066; 4-16 March 2018). All cruises were aboard the R/V F. G. Walton Smith. Cruise region: 76.5W-77.2W, 26.15N-26.8N. Cruises consisted of zonal tethered vertical microstructure profiler (VMP) and vessel mounted acoustic Doppler current profiler (ADCP) surveys, yielding sections of the turbulent dissipation rate (units: W/kg) and zonal and meridional velocity (units: m/s). The VMP used was the tethered Rockland Scientific International VMP-2000, and the ADCP was a vessel-mounted RDI 75 kHz Ocean Surveyor ADCP configured and run through the University of Hawai'i Data Acquisition System (UHDAS) software suite. The VMP was also mounted with a CTD (conductivity-temperature-depth) sensor package, including SBE-3 and SBE-4 sensors. The CTD data was processed according to the manufacturer specifications, using recommended values for the cell thermal mass coefficients (alpha=0.03 and beta=7.0). Files are named as follows: WSxxxxx_CTD.mat, WSxxxxx_vmp_all.mat and WSxxxxx_os75nb.nc for CTD, VMP and ADCP data respectively. Mooring-based data is from two 75 kHz ADCPs, one upward facing and the other downward facing, mounted at approximately 700 m on the RAPID/MOCHA (Rapid Climate Change / Meridional Overturning Circulation and Heat flux Array) mooring WB1 at 76.8W and 26.5N. This configuration gave full depth meridional and zonal velocity profiles (units: m/s) in a water depth of approximately 1300 m. The data spanned December 2015 to March 2017. The data are binned in 16 m depth bins and are bin averaged into 1 hour time bins. File name: WB1_adcp_UV_1hr.mat.

Description: Author Posting. © American Meteorological Society, 2010. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 40 (2010): 789-801, doi:10.1175/2009JPO4039.1.

Description: The issue of internal wave–mesoscale eddy interactions is revisited. Previous observational work identified the mesoscale eddy field as a possible source of internal wave energy. Characterization of the coupling as a viscous process provides a smaller horizontal transfer coefficient than previously obtained, with vh 50 m2 s−1 in contrast to νh 200–400 m2 s−1, and a vertical transfer coefficient bounded away from zero, with νυ + (f2/N2)Kh 2.5 ± 0.3 × 10−3 m2 s−1 in contrast to νυ + (f2/N2)Kh = 0 ± 2 × 10−2 m2 s−1. Current meter data from the Local Dynamics Experiment of the PolyMode field program indicate mesoscale eddy–internal wave coupling through horizontal interactions (i) is a significant sink of eddy energy and (ii) plays an O(1) role in the energy budget of the internal wave field.

Description: Author Posting. © American Meteorological Society, 2008. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 38 (2008): 2556-2574, doi:10.1175/2008JPO3666.1.

Description: Vertical profiles of horizontal velocity obtained during the Mid-Ocean Dynamics Experiment (MODE) provided the first published estimates of the high vertical wavenumber structure of horizontal velocity. The data were interpreted as being representative of the background internal wave field, and thus, despite some evidence of excess downward energy propagation associated with coherent near-inertial features that was interpreted in terms of atmospheric generation, these data provided the basis for a revision to the Garrett and Munk spectral model. These data are reinterpreted through the lens of 30 years of research. Rather than representing the background wave field, atmospheric generation, or even near-inertial wave trapping, the coherent high wavenumber features are characteristic of internal wave capture in a mesoscale strain field. Wave capture represents a generalization of critical layer events for flows lacking the spatial symmetry inherent in a parallel shear flow or isolated vortex.

Description: Salary support for this analysis was provided by Woods Hole Oceanographic Institution bridge support funds.

Description: Author Posting. © American Meteorological Society, 2010. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 40 (2010): 2605–2623, doi:10.1175/2010JPO4132.1.

Description: Steady scale-invariant solutions of a kinetic equation describing the statistics of oceanic internal gravity waves based on wave turbulence theory are investigated. It is shown in the nonrotating scale-invariant limit that the collision integral in the kinetic equation diverges for almost all spectral power-law exponents. These divergences come from resonant interactions with the smallest horizontal wavenumbers and/or the largest horizontal wavenumbers with extreme scale separations. A small domain is identified in which the scale-invariant collision integral converges and numerically find a convergent power-law solution. This numerical solution is close to the Garrett–Munk spectrum. Power-law exponents that potentially permit a balance between the infrared and ultraviolet divergences are investigated. The balanced exponents are generalizations of an exact solution of the scale-invariant kinetic equation, the Pelinovsky–Raevsky spectrum. A small but finite Coriolis parameter representing the effects of rotation is introduced into the kinetic equation to determine solutions over the divergent part of the domain using rigorous asymptotic arguments. This gives rise to the induced diffusion regime. The derivation of the kinetic equation is based on an assumption of weak nonlinearity. Dominance of the nonlocal interactions puts the self-consistency of the kinetic equation at risk. However, these weakly nonlinear stationary states are consistent with much of the observational evidence.

Description: This research is supported by NSF CMG Grants 0417724, 0417732 and 0417466. YL is also supported by NSF DMS Grant 0807871 and ONR Award N00014-09-1-0515.

Description: Author Posting. © American Meteorological Society, 2013. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 43 (2013): 283–300, doi:10.1175/JPO-D-11-0240.1.

Description: Motivated by the recent interest in ocean energetics, the widespread use of horizontal eddy viscosity in models, and the promise of high horizontal resolution data from the planned wide-swath satellite altimeter, this paper explores the impacts of horizontal eddy viscosity and horizontal grid resolution on geostrophic turbulence, with a particular focus on spectral kinetic energy fluxes Π(K) computed in the isotropic wavenumber (K) domain. The paper utilizes idealized two-layer quasigeostrophic (QG) models, realistic high-resolution ocean general circulation models, and present-generation gridded satellite altimeter data. Adding horizontal eddy viscosity to the QG model results in a forward cascade at smaller scales, in apparent agreement with results from present-generation altimetry. Eddy viscosity is taken to roughly represent coupling of mesoscale eddies to internal waves or to submesoscale eddies. Filtering the output of either the QG or realistic models before computing Π(K) also greatly increases the forward cascade. Such filtering mimics the smoothing inherent in the construction of present-generation gridded altimeter data. It is therefore difficult to say whether the forward cascades seen in present-generation altimeter data are due to real physics (represented here by eddy viscosity) or to insufficient horizontal resolution. The inverse cascade at larger scales remains in the models even after filtering, suggesting that its existence in the models and in altimeter data is robust. However, the magnitude of the inverse cascade is affected by filtering, suggesting that the wide-swath altimeter will allow a more accurate determination of the inverse cascade at larger scales as well as providing important constraints on smaller-scale dynamics.

Description: BKA received support from Office of Naval Research Grant N00014-11-1-0487, National Science Foundation (NSF) Grants OCE-0924481 and OCE- 09607820, and University of Michigan startup funds. KLP acknowledges support from Woods Hole Oceanographic Institution bridge support funds. RBS acknowledges support from NSF grants OCE-0960834 and OCE-0851457, a contract with the National Oceanography Centre, Southampton, and a NASA subcontract to Boston University. JFS and JGR were supported by the projects ‘‘Global and remote littoral forcing in global ocean models’’ and ‘‘Agesotrophic vorticity dynamics of the ocean,’’ respectively, both sponsored by the Office of Naval Research under program element 601153N.

Description: 2013-08-01

Description: Author Posting. © American Geophysical Union, 2014. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Oceans 119 (2014): 8627–8645, doi:10.1002/2014JC010099.

Description: One of the most remarkable features of contemporary oceanic climate change is the warming and contraction of Antarctic Bottom Water over much of global ocean abyss. These signatures represent changes in ventilation mediated by mixing and entrainment processes that may be location-specific. Here we use available data to document, as best possible, those mixing processes as Weddell Sea Deep and Bottom Waters flow along the South Orkney Plateau, exit the Weddell Sea via Orkney Passage and fill the abyssal Scotia Sea. First, we find that an abrupt transition in topography upstream of Orkney Passage delimits the extent of the coldest waters along the Plateau's flanks and may indicate a region of especially intense mixing. Second, we revisit a control volume budget by Heywood et al. (2002) for waters trapped within the Scotia Sea after entering through Orkney Passage. This budget requires extremely vigorous water mass transformations with a diapycnal transfer coefficient of inline image m2 s−1. Evidence for such intense diapycnal mixing is not found in the abyssal Scotia Sea interior and, while we do find large rates of diapycnal mixing in conjunction with a downwelling Ekman layer on the western side of Orkney Passage, it is insufficient to close the budget. This leads us to hypothesize that the Heywood budget is closed by a boundary mixing process in which the Ekman layer associated with the Weddell Sea Deep Water boundary current experiences relatively large vertical scale overturning associated with tidal forcing along the southern boundary of the Scotia Sea.

Description: KLP gratefully acknowledges salary support from Woods Hole Oceanographic Institution bridge support funds. ACNG acknowledges the support of a Philip Leverhulme Prize. LJ and MPM were supported by the ANDREX project, funded by the U.K. National Environment Research Council (NE/E01366X/1).

Description: 2015-06-16

Description: Submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy at the Massachusetts Institute of Technology and the Woods Hole Oceanographic Institution September 1992

Description: Oceanic profiles of temperature, salinity, horizontal velocity, rate of dissipation of turbulent kinetic energy (ε) and rate of dissipation of thermal variance (χ) are used to examine the parameterization of turbulent mixing in the ocean due to internal waves. Turbulent mixing is quantified through eddy diffusivity parameterizations of the mass (Kρ; Osborn, 1980) and heat fluxes (Kτ; Osborn and Cox, 1972) in turbulent production/dissipation balances. Turbulence in the ocean is generally held to result from the occurrence of shear instability in regions where the Richardson number is locally supercritical (i.e. Ri ≤ 1/4), permitting the growth of small-scale waves which break and result in turbulent mixing. The occurrence of shear instability results from the local intensification of the shear in the internal wave field. The energy dissipated in such events is provided by the energy flux to higher wavenumber due to nonlinear wave/wave interactions on scales of 10's to 100's of meters. In turn, the strength of the wave/wave interactions depends generally on the energy content of the internal wave field, which can vary considerably over even larger scales due to the presence of topography or background flows. The magnitude of turbulent mixing is linked to internal wave dynamics by equating the turbulent dissipation with the energy flux through the vertical wavenumber spectrum under the priviso that the model spectrum which forms the basis for the analysis is statistically stationary with respect to the nonlinear interactions. Dynamical models (McComas and Muller, 1981; Henyey et al., 1986) indicate that the Garrett and Munk (GM; Munk, 1981) spectrum is stationary. Observations from the far field of a seamount in a region of negligible large-scale flow were examined to address the issue of the buoyancy scaling of ε. These data exhibited large variations in background stratification with depth, but the internal wave characteristics were not substantially differentiable from the GM prescription. The magnitude of ε and its functional dependence upon internal wave energy levels (E) and buoyancy frequency (N) was best described by the dynamical model ofHenyey et al. (1986) (ε ~ E2N2). The Richardson number scaling model of Kunze et al. (1990) produced consistent estimates. A second dynamical model, McComas and Muller (1981), predicted an appropriate (E,N) scaling, but overestimated the observed dissipation rates by a factor of five. Two kinematical dissipation parameterizations (Garmett and Holloway (1984) and Munk (1981)) predicted buoyancy scalings of N3/2 which were inconsistent with the observed scaling. Data from an upper-ocean front, a warm core ring and a region of steep topography were analyzed in order to examine the parameter dependence of E in internal wave fields which exhibited potentially nonstationary characteristics. Evidence was provided which implied the internal wave field in an upper ocean front was interacting with and modified by the background flow. Inhomogeneity and anisotropy of the internal wave field were noted in that data set. The model of Gregg (1989), which in turn was based upon the model of Henyey et al., effectively collapsed the observed diffusivity estimates from the front. The warm core ring profiles were noted to be anisotropic, dominated by near-inertial frequencies and to have a peaked vertical wavenumber shear spectrum. The data from a region of steep topography were noted to have a peaked vertical wavenumber spectrum and were characterized by higher than GM frequency motions. For the latter two data sets, application of a frequency based correction to the Henyey et al. model (Henyey, 1991) reduced more than an order of magnitude scatter in the parameterized estimates of E to less than a factor of four. Of the possible non-equilibrium conditions in the internal wave field, the (E,N) scaled dissipation rates were most sensitive to deviations in wave field frequency content. On the basis of a number of theoretical Richardson number probability distributions (Ri = N2/S2, where S2 is the sum of the squared vertical derivatives of horizontal velocity), the nominal dissipation scaling of the Kunze et al. model was determined to be E2N3. This scaling is altered to the observed ε ~ E2N2 scaling by a statistical dependence between N2 and S2 which reduces the occurrence of supercritical Ri values. This statistical dependence is hypothesized to be an effect of the turbulent momentum and buoyancy fluxes on the internal wave shear and strain profiles caused by shear instability. The statistical dependence between N2 and S2 exhibited a buoyancy scaling which was interpreted as resulting from the decreasing ratio between the time scale of the shear instability mechanism [T- 2π/N] and the adiabatic time scale [T - 2π/(Nf)1/2] of the internal wave field (f is the Coriolis parameter). This phenomenology is interpreted in light of saturated spectral theories which suggest that the magnitude and shape of the vertical wavenumber spectrum is controlled by instability mechanisms at large wavenumber ( ≥ .1 cpm). We argue that saturated spectral theories are valid only in the limit where a separation exists between the two time scales, i.e. for large N, low internal wave frequency content, and small f. These results have immediate implications for oceanic mixing driven by internal wave motions. First, background diffusivities are small: at GM energy levels, Kρ - .03x10-4 m2/s (Kρ = .25ε/N2). Secondly, since Kρ is independent of N at constant E, some process or collection of processes must be responsible for heightened E values in the abyss if internal waves cause the 0(1-10x10-4 m2/s) diffusivities generally inferred from deep ocean hydrographic data. We view internal wave reflection and/or internal wave generation associated with topographic features to be likely candidates.

Description: Author Posting. © American Meteorological Society, 2015. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 45 (2015): 966–987, doi:10.1175/JPO-D-14-0110.1.

Description: A key remaining challenge in oceanography is the understanding and parameterization of small-scale mixing. Evidence suggests that topographic features play a significant role in enhancing mixing in the Southern Ocean. This study uses 914 high-resolution hydrographic profiles from novel EM-APEX profiling floats to investigate turbulent mixing north of the Kerguelen Plateau, a major topographic feature in the Southern Ocean. A shear–strain finescale parameterization is applied to estimate diapycnal diffusivity in the upper 1600 m of the ocean. The indirect estimates of mixing match direct microstructure profiler observations made simultaneously. It is found that mixing intensities have strong spatial and temporal variability, ranging from O(10−6) to O(10−3) m2 s−1. This study identifies topographic roughness, current speed, and wind speed as the main factors controlling mixing intensity. Additionally, the authors find strong regional variability in mixing dynamics and enhanced mixing in the Antarctic Circumpolar Current frontal region. This enhanced mixing is attributed to dissipating internal waves generated by the interaction of the Antarctic Circumpolar Current and the topography of the Kerguelen Plateau. Extending the mixing observations from the Kerguelen region to the entire Southern Ocean, this study infers a large water mass transformation rate of 17 Sverdrups (Sv; 1 Sv ≡ 106 m3 s−1) across the boundary of Antarctic Intermediate Water and Upper Circumpolar Deep Water in the Antarctic Circumpolar Current. This work suggests that the contribution of mixing to the Southern Ocean overturning circulation budget is particularly significant in fronts.

Description: AM was supported by the joint CSIRO–University of Tasmania Quantitative Marine Science (QMS) program and the 2009 CSIRO Wealth from Ocean Flagship Collaborative Fund. BMS was supported by the Australian Climate Change Science Program, jointly funded by the Department of the Environment and CSIRO. KLPs salary support was provided by Woods Hole Oceanographic Institution bridge support funds.

Description: 2015-10-01

Description: © The Author(s), 2015. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Journal of Geophysical Research: Oceans 120 (2015): 2096–2111, doi:10.1002/2014JC010573.

Description: In this study, we address the question whether eddy-driven transports in the Argentine Basin can be held responsible for enhanced sediment accumulation over the Zapiola Rise, hence accounting for the existence and growth of this sediment drift. To address this question, we perform a 6 year simulation with a strongly eddying ocean model. We release two passive tracers, with settling velocities that are consistent with silt and clay size particles. Our experiments show contrasting behavior between the silt fraction and the lighter clay. Due to its larger settling velocity, the silt fraction reaches a quasisteady state within a few years, with abyssal sedimentation rates that match net input. In contrast, clay settles only slowly, and its distribution is heavily stratified, being transported mainly along isopycnals. Yet, both size classes display a significant and persistent concentration minimum over the Zapiola Rise. We show that the Zapiola Anticyclone, a strong eddy-driven vortex that circulates around the Zapiola Rise, is a barrier to sediment transport, and hence prevents significant accumulation of sediments on the Rise. We conclude that sediment transport by the turbulent circulation in the Argentine Basin alone cannot account for the preferred sediment accumulation over the Rise. We speculate that resuspension is a critical process in the formation and maintenance of the Zapiola Rise.

Description: This research was supported by the Regional and Global Climate Modeling Program of the US Department of Energy Office of Science (WW). Los Alamos National Laboratory is operated by the Los Alamos National Security, LLC for the National Nuclear Security Administration of the U.S. Department of Energy under contract DE-AC52-06NA25396.

Description: Author Posting. © American Geophysical Union, 2015. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Geophysical Research Letters 42 (2015): 3431–3441, doi:10.1002/2015GL063494.

Description: The impact on the upper ocean of the passage of a short, intense storm over a Southern Ocean site, in proximity to an Antarctic Circumpolar Current front, is characterized. The storm causes a wind-induced deepening of the mixed layer and generates an inertial current. Immediate poststorm observations indicate a mixed layer extending to approximately 50 m depth. Subsequent measurements show the upper ocean to have restratified, injecting near-inertial shear in stratified waters within 1 day of the storm's passage. This time scale for the development of near-inertial shear is 1 order of magnitude shorter than that predicted by the β dispersion paradigm. The observed rapid changes in upper ocean stratification point to the existence of an as yet undocumented, efficient mechanism for injection of near-inertial shear into the stratified ocean that is in turn associated with enhanced turbulence and mixing.

Description: The SOFine project is funded by the UK Natural Environmental Research Council (NERC) (grant NE/G001510/1).

Description: 2015-11-07