Printer Friendly

CLIMATE PROCESS TEAM ON INTERNAL WAVE-DRIVEN OCEAN MIXING: The study summarizes recent advances in our understanding of internal wave-driven turbulent mixing in the ocean interior and introduces new parameterizations for global climate ocean models and their climate impacts.

Diapycnal mixing plays a primary role in the thermodynamic balance of the ocean and, consequently, in oceanic heat and carbon uptake and storage. Though observed mixing rates are on average consistent with values required by inverse models, recent attention has focused on the dramatic spatial variability, spanning several orders of magnitude, of mixing rates in both the upper and deep ocean. Away from ocean boundaries, the spatiotemporal patterns of mixing are largely driven by the geography of generation, propagation, and dissipation of internal waves, which supply much of the power for turbulent mixing. Over the last 5 years and under the auspices of U.S. Climate Variability and Predictability Program (CLIVAR), a National Science Foundation (NSF)- and National Oceanic and Atmospheric Administration (NOAA)-supported Climate Process Team has been engaged in developing, implementing, and testing dynamics-based parameterizations for internal wave-driven turbulent mixing in global ocean models. The work has primarily focused on turbulence 1) near sites of internal tide generation, 2) in the upper ocean related to wind-generated near inertial motions, 3) due to internal lee waves generated by low-frequency mesoscale flows over topography, and 4) at ocean margins. Here, we review recent progress, describe the tools developed, and discuss future directions.

**********

Ocean turbulence influences the transport of heat, freshwater, dissolved gases such as C[O.sub.2], pollutants, and other tracers. It is central to understanding ocean energetics and reducing uncertainties in global circulation and simulations from climate models. The dissipation of turbulent energy in stratified water results in irreversible diapycnal (across density surfaces) mixing. Recent work has shown that the spatial and temporal inhomogeneity in diapycnal mixing may play a critical role in a variety of climate phenomena. Hence, a quantitative understanding of the physics that drive the distribution of diapycnal mixing in the ocean interior is fundamental to understanding the oceans role in climate.

Diapycnal mixing is very difficult to accurately parameterize in numerical ocean models for two reasons. The first one is due to the discrete representation of tracer advection in directions that are not perfectly aligned with isopycnals, which can result in numerically induced mixing from truncation errors that is larger than observed diapycnal mixing (Griffies et al. 2000; Ilicak et al. 2012). The second reason is related to the intermittency of turbulence, which is generated by complex and chaotic motions that span a large space-time range. Furthermore, this mixing is driven by a wide range of processes with distinct governing physics that create a rich global geography [see MacKinnon et al. (2013c) for a review]. The difficulty is also related to the relatively sparse direct sampling of ocean mixing, whereby sophisticated ship-based measurements are generally required to accurately characterize ocean mixing processes. Nonetheless, we have sufficient evidence from theory, process models, laboratory experiments, and field measurements to conclude that away from ocean boundaries (atmosphere, ice, or the solid ocean bottom), diapycnal mixing is largely related to the breaking of internal gravity waves, which have a complex dynamical underpinning and associated geography. Consequently, in 2010, a Climate Process Team (CPT), funded by the National Science Foundation (NSF) and the National Atmospheric and Oceanic Administration (NOAA), was convened to consolidate knowledge on internal wave-driven turbulent mixing in the ocean, develop new and more accurate parameterizations suitable for global ocean models, and consider the consequences for global circulation and climate. Here, we report on the major findings and products from this CPT.

Ocean internal gravity waves propagate through the stratified interior of the ocean. They are generated by a variety of mechanisms, with the most important being tidal flow over topography, wind variations at the sea surface, and flow of ocean currents and eddies over topography leading to lee waves (see schematic in Fig. 1). As waves propagate horizontally and vertically away from their generation sites, they interact with each other, producing an internal gravity wave continuum consisting of energy in many frequencies and wavenumbers. The waves with high vertical wavenumbers (small vertical scales) are more likely to break, leading to turbulent mixing. The distribution of diapycnal mixing therefore depends on the entire chain of processes shown in Fig. 1.

A brief history of vertical mixing parameterizations used by ocean models. Ocean models often approximate diapycnal mixing processes through vertical Fickian diffusion, which takes the mathematical form

Fickian diffusion = [partial derivative]/[partial derivative]z ([kappa] [partial derivative][psi]/[partial derivative]z), (1)

where [psi] is the tracer concentration, z is the geopotential vertical coordinate, and k is the diapycnal diffusivity (dimensions of [L.sup.2][T.sup.-1], where I is length and Tis time). Through the 1990s, global models routinely used space-time constant vertical diffusivities. A notable exception was Bryan and Lewis (1979), who prescribed a horizontally uniform diffusivity that increased with depth, reflecting the observed larger vertical mixing in the deep ocean and reduced mixing in the pycnocline. By the mid-1990s, ocean climate models began to separate diapycnal mixing into surface boundary layer and interior processes. In and near the surface boundary layer, mixing is controlled by a balance between buoyancy input (e.g., heat and freshwater fluxes) and mechanical forcing (e.g., wind) that establish the surface boundary layer and fluxes through it. Climate models of this era used boundary layer schemes such as Gaspar et al. (1990) and Large et al. (1994). In the stably stratified ocean interior, both shear-driven mixing (Pacanowski and Philander 1981; Large et al. 1994) and double-diffusive processes (Large et al. 1994) were parameterized. Gravitational instabilities giving rise to vertical convection were accounted for through a large vertical diffusivity (Large et al. 1994; Klinger et al. 1996) or a convective adjustment scheme (Rahmstorf 1993).

In the deep ocean, a prognostic parameterization for internal tide-driven mixing was introduced by St. Laurent et al. (2002), who combined an estimate of internal tide generation over rough topography with an empirical vertical decay scale for the enhanced turbulence (see the section on "Nearfield tidal mixing"). State-of-the-art ocean climate simulations prior to the CPT, as represented by the Geophysical Fluid Dynamics Laboratory (GFDL) and National Center for Atmospheric Research (NCAR) phase 5 of the Coupled Model Intercomparison Project (CMIP5) simulations (Dunne et al. 2012; Danabasoglu et al. 2012), included a version of Eq. (3) (see the section on "Near-field tidal mixing"), along with parameterizations of mixing in the surface (Large et al. 1994) and bottom boundary layers and/ or overflows (Legg et al. 2006; Danabasoglu et al. 2010) and mixing from resolved shear (Large et al. 1994; Jackson et al. 2008). These parameterizations produced spatially and temporally varying diapycnal diffusivities, with bottom enhancement and stratification dependence. However, these simulations did not include an energetically consistent representation of internal tide breaking away from the generation site, explicit representation of mixing from internal waves generated by winds and subinertial flows, nor spatial and temporal variability in the dissipation vertical profile. The work described here has revolved around developing and testing energetically consistent, spatially and temporally variable mixing parameterizations. The resulting parameterizations are based upon internal gravity wave dynamics and the patterns of wave generation, propagation, and dissipation.

Overall strategy and philosophy of the CPT approach. As with previous CPTs, we have found that parameterizations are most productively developed when there is a broad base of knowledge that is in a state of readiness to be consolidated, implemented, and tested. Much of the basic research described here was published or nearing completion at the time this project started, allowing for a focused effort on parameterization development, model implementation, and global model testing. A key CPT component was the inclusion of four dedicated postdoctoral scholars, who formed "the glue" to bridge the expertise of different principal investigators, promoting projects at the intersection of theory and models, observations, and simulations, while gaining valuable broad training and networking.

One of the important tenets of the CPT is the consistent use of energy, power, and the turbulent kinetic energy (KE) dissipation rate [epsilon] (dimensions of [L.sup.2][T.sup.-3]), rather than diapycnal diffusivity, as the currency of turbulent mixing; [epsilon] describes the rate at which turbulence dissipates mechanical energy at the smallest scales. It is typically related to a diapycnal diffusivity through a dimensionless mixing efficiency [GAMMA], following Osborn (1980):

[kappa] = [GAMMA][epsilon]/[N.sup.2], (2)

where [N.sup.2] is the squared buoyancy frequency. Equation (2) shows that keeping the diffusivity fixed in a world with changing stratification implies changes in energy dissipation in ways that are not always consistent with the physical processes supplying energy for dissipation. We can overcome this problem by formulating parameterizations directly in terms of [epsilon]. This approach also has the advantage of providing a transparent connection to dynamical processes driving mixing, since the downscale energy cascade can be directly linked to constraints of total power available for turbulence and other facets of ocean energetics (e.g., St. Laurent and Simmons 2006; Ferrari and Wunsch 2009). The topic of an appropriate value for mixing efficiency has had a resurgence of interest in recent years. Some theoretical and numerical studies suggest that a mixing efficiency that is systematically lower in areas of low ocean stratification might bias the type of global mixing estimates presented here and require modifications to model parameterizations (Mashayek et al. 2013; Venayagamoorthy and Koseff 2016; Salehipour et al. 2016). A careful evaluation of mixing efficiency was not part of the CPT work, and a thorough discussion is beyond the scope of this paper. Interested readers are instead referred to recent reviews such as Peltier and Caulfield (2003) and Gregg et al. (2018).

GLOBAL PATTERNS AND CONSTRAINTS.

Many of the early parameterizations described in the section titled "A brief history of vertical mixing parameterizations used by ocean models" were motivated by individual process experiments or observational studies. At the same time, the novel observations, theories, and model results that fundamentally drive the field forward frequently arise unexpectedly from programs funded by many agencies. For example, the long-range propagation of coherent internal tides was discovered in both the Acoustic Thermometry of Ocean Climate (ATOC; Dushaw et al. 1995) and satellite altimeter (Ray and Mitchum 1996) datasets fortuitously; neither mission was set up with a focus on internal tides.

Another factor contributing to the readiness of this CPT was the increased use of new techniques to infer mixing rates indirectly from a wide variety of data sources, allowing the rich patterns like those in Fig. 2 to emerge. There are now enough direct microstructure and indirect estimates of turbulent dissipation rates and diapycnal diffusivities to examine depth and geographical patterns, temporal variability, and global budgets (Waterhouse et al. 2014). These patterns in turn have inspired new insights on the underlying dynamics driving and energetically supplying small-scale turbulence and provided valuable constraints on modeled turbulent mixing rates. Compilation of both direct microstructure measurements and indirect estimates of turbulence is discussed in the section titled "Tools and techniques." Here, we briefly describe recent results related to global patterns and statistics.

The average strength of turbulent diapycnal mixing appears to be roughly consistent, within error bars, with that required to raise the deep waters of the global meridional overturning circulation (MOC). Using the most comprehensive-to-date collection of full-depth microstructure data, Waterhouse et al. (2014) report a globally averaged diapycnal diffusivity below 1000-m depth of O([10.sup.-4]) [m.sup.2] [s.sup.-1] and above 1000-m depth of O([10.sup.-5]) [m.sup.2] [s.sup.-1]. These values are consistent with the global inverse estimate of Lumpkin and Speer (2007). Using an indirect finescale approach (see the section on "Observational data analysis: The finescale parameterizations"), but with a much larger dataset, Kunze (2017) finds a global depth-averaged value of 0.3-0.4 x [bar.x]-4 [m.sup.2] [s.sup.-1]. It is unclear whether any remaining differences between these estimates are due to sampling biases of the more limited microstructure data, to method biases of the finescale technique, or to assumptions of a fixed mixing efficiency.

The associated globally averaged turbulent dissipation rates inferred from these observations cluster around 2 [+ or -] 0.6 TW (Waterhouse et al. 2014; Kunze 2017). Given an assumed mixing efficiency, these rates are roughly consistent with estimates of power going through the three primary mechanisms of internal wave generation: barotropic tidal flow over topography leading to internal tides (~1 TW; see the sections on "Near-field tidal mixing" and "Far-field internal tides"), low-frequency flows over topography producing internal lee waves (0.2-0.7 TW; see section on "Internal lee waves"), and variable wind forcing producing near-inertial internal waves (0.3-1 TW; see section on "Wind-driven near-inertial motions").

Much more striking than average values is the enormous range and richness of the patterns visible in Fig. 2. Both the turbulent dissipation rate and diapycnal diffusivity vary by several orders of magnitude across ocean basins. Understanding how such patterns convolve with pathways of water mass movement, air-sea heat gain/loss, greenhouse gas input, and nutrient availability is the next frontier in interpreting diapycnal mixing in the ocean.

Many of these patterns (in space and time) can be interpreted in terms of the geography of internal wave generation, propagation, and dissipation (Fig. 1). Patterns immediately visible in Fig. 2 include elevated values associated with more complex topography such as that associated with the western Indian Ocean, western and central Pacific Ocean, and slowmidocean spreading ridges (Wijesekera et al. 1993; Polzin et al. 1997; Kunze et al. 2006; Decloedt and Luther 2010; Wu et al. 2011; Whalen et al. 2012; Waterhouse et al. 2014). Over rough or steep topography, turbulence is frequently bottom enhanced (Polzin et al. 1997; Waterhouse et al. 2014) but sometimes extends all the way up through the pycnocline (Kunze 2017). The temporal variability of diapycnal mixing shows seasonal (Whalen et al. 2012) and tidal cycles related to the two major internal wave energy sources, the winds and tides, as well as isolated events.

What follows in the sections below concerns first the main science efforts to consolidate our understanding of turbulence from 1) mixing elevated over rough topography related to internal wave generation by tides, 2) low-frequency flows that generate internal lee waves, and 3) near-inertial internal wave generation by winds. In each section, we describe the consequences of parameterizing these processes in ocean climate models. For tides, we subdivide our efforts into turbulence in the near field of internal tide generation sites (loosely within one mode-1 bounce) and the far field (waves that have propagated considerably farther before breaking). Following that we describe tools developed through the CPT now made available to the wider community, namely, 1) a uniquely comprehensive database of microstructure data, 2) techniques for analyzing observational data, and 3) new parameterizations of turbulence available for a variety of model implementations. We also briefly discuss the state of the art for high-resolution ocean models, which are beginning to partially resolve the internal gravity wave continuum on a global scale. We conclude this paper with thoughts for the future.

NEAR-FIELD TIDAL MIXING. Physical motivation. Tidal frequency internal waves, generated by barotropic tidal flow over topographic obstacles in a stably stratified fluid, lead to local mixing near the generation site, both due to direct wave breaking (close to topography) and enhanced rates of interaction with other internal waves (well above topography). The formulation of St. Laurent et al. (2002) represented the enhanced turbulent dissipation rate as the product of the rate of conversion of barotropic tidal energy into internal waves C, the fraction of that energy that is locally dissipated q (note that consequently 1 - q propagates away as low-mode internal tides), and a vertical distribution function of that local dissipation F(z). Through the Osborn relation in Eq. (2) (Osborn 1980), the enhanced turbulence is then related to a diffusivity as

[kappa] = [[kappa].sub.b] + q[GAMMA]C(x, y)F(z)/[rho][N.sup.2], (3)

where [[kappa].sub.b] is a placeholder background diffusivity. The conversion rate C is dependent on topographic roughness, tidal velocity, and bottom stratification (Bell 1975; Jayne and St. Laurent 2001; Garrett and Kunze 2007; Fig. 3c). St. Laurent et al. (2002) proposed a value of q = 1/3 and a function F(z) that decayed exponentially with height above topography, with a 500-m e-folding scale. They based these choices on analysis from several deep-ocean microstructure datasets. These values were used in climate model implementations, such as Simmons et al. (2004b), Jayne (2009), Dunne et al. (2012), and Danabasoglu et al. (2012). The background diffusivity [[kappa].sub.b] accounts for the mixing associated with energy that radiates from internal tide generation sites as well as other internal wave processes. Treatments of [[kappa].sub.b] have varied, including 1) a constant value of 1 x [10.sup.-5] [m.sup.2] [s.sup.-1] (Simmons et al. 2004b; Jayne 2009), 2) a latitudinal function capturing the equatorward decrease in wave-wave interactions (Henyey et al. 1986; Harrison and Hallberg 2008; Jochum 2009; Danabasoglu et al. 2012), and 3) a stratification-dependent function after Gargett (1984) [used in Dunne et al. (2012)]. Because of the sensitivity of the simulations to the different parameterizations, a major goal of the CPT has been to better understand and represent the physical processes that determine spatial and temporal variations in the parameters in Eq. (3).

A few estimates of q have been obtained, involving synthesis of observations and models. The radiated portion 1 - q maybe computed as the energy radiated out of a control volume [integral]J x [??]dA, where J is the internal wave energy flux, divided by an estimate of the conversion rate C. Alternately, a direct estimate is from the integrated dissipation rate over that same volume [integral][rho][GAMMA][epsilon]dV/C. The observational sampling requirements for both estimates, particularly the second, are considerable. At the Hawaiian Ridge, Klymak et al. (2006) obtained q = 0.15 using the second method, as compared to an estimate of q < 0.5 obtained with the first (Rudnick et al. 2003).

Existing theoretical predictions for C, summarized in Garrett and Kunze (2007) and Green and Nycander (2013), show dependence on topographic steepness relative to the internal tide characteristic steepness [gamma] = (dh/dx)/s [where

s = [square root of ([f.sup.2] - [[omega].sup.2])/([N.sup.2] - [[omega].sup.2]))],

dh/dx is the topographic gradient, [omega] is the wave frequency, and/ is the Coriolis parameter] as well as the ratio of tidal excursion distance to topographic width. At supercritical rough topography ([gamma] > 1) the conversion rate saturates (Balmforth and Peacock 2009; Zhang and Swinney 2014) compared to linear theory applicable at subcritical topography ([gamma] < 1) (Bell 1975). Estimates of C need to include the contribution of abyssal hill topography on scales O(10) km not resolved by current topography products. Small-scale topography may increase C by 10% globally and 100% regionally (Melet et al. 2013b; see Fig. 3c).

A global constraint on the near-field internal tide dissipation can be obtained from comparisons of satellite observations of internal tides with global simulations at 0(10)-km resolution that include realistic surface tidal forcing (Simmons et al. 2004a; Arbic et al. 2004, 2010; Niwa and Hibiya 2011; Muller et al. 2012; Shriver et al. 2012; Niwa and Hibiya 2014; Shriver et al. 2014; Waterhouse et al. 2014; Ansong et al. 2015; Buijsman et al. 2016; Rocha et al. 2016). All of these model runs explicitly simulate generation of low-mode tides, with horizontal scales > O(10) km. Some studies conducted since 2010 have also included concurrent atmospheric forcing, allowing for a more realistic, geographically varying background stratification field. In some of the models above, conversion to unresolved high modes, assumed to dissipate locally, is performed by a linear wave drag based on linear theory (Bell 1975). Buijsman et al. (2016) find that modeled and observed internal tides show the most agreement when about 60% of the energy converted to both low and high modes is dissipated close to the generation sites.

The vertical structure of associated turbulence appears to vary between deep rough topography and tall steep topography, reflecting differences in the underlying physics driving turbulence. At tall steep ridges, much of the baroclinic energy is contained in larger length scales that propagate away horizontally without breaking (St. Laurent and Nash 2004). Local mixing occurs through tidally generated, transient, arrested lee waves (Legg and Klymak 2008; Klymak et al. 2010; Alford et al. 2014; Fig. 3b), which might imply a q scaling with the barotropic flow speed U and an exponentially decaying vertical dissipation profile with length scale U/N. At the Kaena Ridge, Hawaii, this theory suggests q ~ 7%, which is less than the q ~ 15% values estimated from observations (Klymak et al. 2006). Interference with remotely generated internal tides modifies the local dissipation (Buijsman et al. 2012, 2014; Klymak et al. 2013); resonance between internal tides generated at adjacent ridges (e.g., Luzon Straits) can increase local dissipation up to 40% (Alford et al. 2015). The percentage of local dissipation may be systematically higher in marginal seas or areas where lower modes are not free to escape (St. Laurent 2008; Nagai and Hibiya 2015). Similarly, near-field tidal dissipation can be increased by topographically trapped internal waves generated by subinertial tidal constituents (Tanaka et al. 2013), that is, the diurnal constituents at latitudes >30[degrees] and the semidiurnal constituents at latitudes >74.5[degrees]. The energy density in such trapped motions increases with latitude and is all dissipated locally (Musgrave et al. 2016).

At deep rough topography a variety of processes facilitate local wave breaking (Fig. 3a). Wave-wave interactions can transfer energy to smaller-scale waves that are more likely to break (McComas 1977; Muller et al. 1986; Henyey et al. 1986). This process is modeled in Polzin (2004b) with a one-dimensional radiation balance equation, resulting in an algebraically decaying dissipation profile with a spatially varying decay scale that matches Brazil Basin observations (Polzin et al. 1997; Fig. 3d). For small-scale waves generated over subcritical abyssal hill topography, overturning of the upward-propagating waves (Muller and Buhler 2009) predicts a bottom-intensified dissipation, with a steeper than exponential decay with height and a local dissipation fraction as large as 60%. At and just below a critical latitude where the Coriolis frequency is half the tidal frequency, particularly efficient wave-wave interactions of a parametric subharmonic instability type lead to a dissipation profile with high values extending several hundred meters above the bottom, before decaying rapidly to background levels, and q > 0.4 (MacKinnon and Winters 2003; Ivey et al. 2008; Nikurashin and Legg 2011). Internal tide energy can also be transferred to smaller scales in the pycnocline and by scattering from rough topography following reflection from the upper surface (Buhler and Holmes-Cerfon 2011). The value of q = 0.3 used in existing parameterizations is therefore likely to be an underestimate in many places, while an overestimate in some.

New parameterizations. A major effort in the CPT and elsewhere has been to build upon the work of Jayne and St. Laurent (2001) and St. Laurent et al. (2002) by deriving more dynamically variable and accurate representations of both the decay profile F{z) and the fraction of locally dissipated wave energy q. For deep rough topography, Polzin (2009) formulates a parameterization of internal tide dissipation based on ID radiation balance equations with nonlinear closure. His formulation yields a dissipation that scales like [epsilon] = [[epsilon].sub.0]/[(1 + z/[z.sub.p]).sup.2], where z is the height above bottom (Fig. 3d). In Melet et al. (2013a), the scale height [z.sub.p] is written in the form

[z.sub.p] = [mu][U[([N.sub.b.sup.ref]).sup.2]/[h.sup.2][k.sup.2][N.sub.b.sup.3], (4)

where [mu] is a nondimensional constant, [N.sub.b.sup.ref] is a reference bottom buoyancy frequency, and U, h, k, and [N.sub.b] are, respectively, the barotropic velocity, topographic roughness, topographic wavenumber, and bottom buoyancy frequency for the particular location. Wentzel-Kramers-Brillouin (WKB) scaling contributes to the role of stratification in Eq. (4). Another global map of q and vertical profile of dissipation for small-scale rough topography has been generated by Lefauve et al. (2015) using the overturn mechanism of Muller and Biihler (2009).

For turbulence at tall, steep slopes, a new parameterization of the near-field mixing due to transient arrested lee waves (Klymak et al. 2010) uses linear theory for knife-edge ridge topography to estimate baroclinic energy conversion into each mode (Llewellyn Smith and Young 2003). Those modes with phase speeds less than the barotropic velocity at the top of the ridge are assumed to be arrested, leading to local dissipation. Combining the total energy loss with a vertical length scale of U/N produces a dissipation rate that decays exponentially away from the ridge top.

Consequences for large-scale circulation. Melet et al. (2013a) compare two simulations with the same formulation for internal tide energy input but using different vertical profiles of dissipation [the St. Laurent et al. (2002) and Polzin (2009) formulations, also included in the Community Earth System Model (CESM)]. They used the GFDL Climate Model, version 2G (CM2G), coupled climate model with an isopycnal vertical coordinate in the ocean (Dunne et al. 2012). With the Polzin formulation, diffusivities are higher around 1000-1500 m and lower in the deep ocean, resulting in modifications to the ocean stratification and changes of O(10%) in the meridional overturning circulation (Fig. 3e).

Additional enhancements in the CESM ocean component, meant to improve the representation of tidally driven mixing, include separate treatment of diurnal and semidiurnal tidal constituents and implementation of a subgrid-scale bathymetry parameterization that better resolves the vertical distribution of the barotropic energy flux, following Schmittner and Egbert (2014); alternative tidal dissipation energy datasets from Egbert and Ray (2003) and Green and Nycander (2013); and introduction of the 18.6-yr lunar nodal cycle on the tidal energy fields. The global climate impacts of these new enhancements are found to be rather small. However, there are local improvements such as a reduction in the warm bias in the upper ocean in the Kuril Strait region.

Future work. Ongoing work is synthesizing existing ideas for the dependence of q on topographic and flow parameters into a single global model for a spatially and temporally varying q and incorporating these ideas into simulations. Comparison with additional observations of the strength and vertical decay scale of turbulence over rough topography is also desirable. For example, Kunze (2017) finds that inferred dissipation rates over some topographic features extend upward well into the thermocline without appreciable decay. Parameterization of mixing by trapped tidally forced waves (perhaps especially important in the Arctic and Antarctic) also deserves dedicated attention.

FAR-FIELD INTERNAL TIDES. About 20%-80% of the internal tide energy is not dissipated near topographic sources (see the section on "Near-field tidal mixing") and instead radiates away as low-mode internal waves. Satellite altimetry shows that these low-mode internal tides may propagate for thousands of kilometers from sources such as the Hawaiian Ridge (Fig. 4a; Zhao et al. (2016)). This section examines where and how these low modes dissipate and parameterizations of this dissipation. Several mechanisms have been hypothesized as potential dissipators of far-field internal tides, including interactions with rough topography (Johnston and Merrifield 2003; Mathur et al. 2014); interactions with mean flows and eddies (St. Laurent and Garrett 2002; Rainville and Pinkel 2006; Dunphy and Lamb 2014; Kerry et al. 2014); cascade to smaller scales via wave-wave interactions (McComas 1977; Muller et al. 1986; Henyey et al. 1986; Lvov et al. 2004; Polzin 2004a), including the particular subset of wave interactions known as parametric subharmonic instability (PSI; Staquet and Sommeria 2002; MacKinnon and Winters 2005; Alford et al. 2007; Alford 2008; Hazewinkel and Winters 2011; MacKinnon et al. 2013a,b; Simmons 2008; Sun and Pinkel 2012, 2013); or evolution on continental slopes and shelves (Nash et al. 2004,2007; Martini et al. 2011; Kelly et al. 2013; Waterhouse et al. 2014). Here, we summarize current understanding from theoretical and process studies and observational campaigns, recent parameterization developments, and consequences of far-field dissipation for global ocean models.

Observations. The reflection, scattering, and dissipation of long-range, low-mode internal tides have been observed at a few large topographic features. Satellite altimetry indicates scattering of mode-1 tides to higher modes along the Line Islands Ridge, 1000 km south of Hawaii (Johnston and Merrifield 2003). Moored observations show significant reflection for mode-1 diurnal internal tides (but weak reflection for semidiurnal) at the South China Sea continental shelf (Klymak et al. 2011). Scattering of internal tides from low to high modes and associated mixing have been observed on the Virginia and Oregon continental slopes (Nash et al. 2004; Kelly et al. 2012; Martini et al. 2013). In contrast, at the steeper Tasmanian continental slope, mode-1 internal tides appear to reflect without significant energy loss (Johnston et al. 2015).

Theory and numerical simulations. The interaction between low-mode internal waves and largeamplitude topography, such as continental slopes or tall isolated ridges, is strongly dependent on the steepness of the topography (Cacchione and Wunsch 1974; Johnston and Merrifield 2003; Legg and Adcroft 2003; Venayagamoorthy and Fringer 2006; Helfrich and Grimshaw 2008; Hall et al. 2013; Legg 2014; Mathur et al. 2014). Shoaling subcritical topography can increase wave amplitude, increasing the Froude number (defined in the section on "Internal lee waves") and causing wave breaking. Supercritical topography reflects low-mode waves back toward deeper water, with only small energy loss to dissipation (Klymak et al. 2013). Near-critical topography scatters incident low-mode energy to much smaller wavelengths, leading to wave breaking and turbulence (Wunsch 1969; Ivey and Nokes 1989; Slinn and Riley 1996; Ivey et al. 2000) concentrated near the sloping topography. Kelly et al. (2013) estimated the fraction of incoming mode-1 energy flux transmitted, reflected and scattered into higher modes for two-dimensional sections across the continental slope for the entire global coastline. Three-dimensional topographic variations such as canyons, cross-slope ridges and troughs, and bumps may enhance the local dissipation of the low-mode tide.

Parameterizing far-field tides: A wave drag approach. In global simulations of the Hybrid Coordinate Ocean Model (HYCOM) with realistic atmospheric and tidal forcing (Arbic et al. 2010), the resolved internal waves lose energy to a wave drag applied to flow in the bottom 500 m (see the section on "Near-field tidal mixing"). This drag can be regarded as a parameterization of low- to high-mode scattering, and these high modes are assumed to dissipate at the generation site, within 500 m above the bottom topography. Comparison of the simulated M2 internal tide SSH amplitudes in 1/12.5[degrees] HYCOM with satellite altimetry (Shriver et al. 2012; Ansong et al. 2015; Buijsman et al. 2016) shows that the open-ocean wave drag is necessary to achieve agreement between modeled and observed barotropic and baroclinic tides, confirming the need for deep-ocean dissipation of the low-mode internal tides. Figures 4b and 4c, taken from Ansong et al. (2017), display the internal tide conversion rates and fluxes in HYCOM, and the comparison of HYCOM fluxes to fluxes in high-vertical-resolution moorings in the North Pacific (Zhao et al. 2010). Consistent with earlier studies, such as Simmons et al. (2004a), the conversion map shows that internal tides are generated in areas of rough topography such as the Hawaiian Ridge. The H YCOM-mooring comparison map in Fig. 4c indicates that the HYCOM simulations are able to predict tidal fluxes with some reasonable degree of accuracy. Buijsman et al. (2016) found that about 12% of these low modes reach the continental slopes, compared to 31% found by Waterhouse et al. (2014). The HYCOM results cited above suggest the necessity of parameterized energy loss, but the current wave drag formulation used in HYCOM is based only upon topographic scattering, motivating additional studies to understand a greater number of relevant physical mechanisms implicated in the damping of far-field internal tides.

Parameterizing far-field internal tides: A ray-tracing approach. To represent the geography of far-field internal tide dissipation in a physically based manner, the propagation, reflection, and dissipation of low-mode energy must be parameterized in a GCM. A new numerical framework employs a vertically integrated radiation balance equation to predict the horizontal propagation of low-mode energy, simplifying earlier surface and internal wave modeling (e.g., Hasselmann et al. 1988; Muller and Natarov 2003). In this approach, only the lowest modes are considered. Energy in each mode of each relevant tidal frequency is considered independently (or adiabatically), assuming minimal mode-mode energy transfer. Waves propagate horizontally with refraction due to variations in Coriolis, depth, and stratification, invoking classic ray-tracing equations for long internal gravity waves (Lighthill 1978). Effects of background flow (Rainville and Pinkel 2006) are currently neglected but will be included in future versions. The 1 - q fraction of the outgoing internal tide energy that does not dissipate locally (see the section on "Near-field tidal mixing") forms the source term in the radiation balance equation, and various parameterizations for dissipation can be plugged into the framework as sink terms. Dissipation mechanisms currently considered include scattering at small-scale roughness (Jayne and St. Laurent 2001), quadratic bottom drag [similar to some of the simulations in Ansong et al. (2015)], and Froude number-based breaking (Legg 2014). A scheme for partial reflection at continental slopes uses the reflection coefficients of Kelly et al. (2013). This framework, currently implemented in GFDL's Modular Ocean Model (MOM6) can be adapted or extended to incorporate new parameterizations of sink and source phenomena. Eden and Olbers (2014) have developed a similar approach for propagating low-mode energy, with scattering to a high-mode continuum due to wave-wave interaction and topographic roughness (not including reflection at continental slopes).

Consequences of far-field dissipation in GCMs. To examine the sensitivity of large-scale ocean circulation to the location of far-field internal tide dissipation, a series of simulations were performed with the Geophysical Fluid Dynamics Laboratory Earth System Model with GOLD component (GFDL-ESM2G) coupled climate model (Dunne et al. 2012). These simulations (Melet et al. 2016) all have the same total energy input into the internal tide field and the same magnitude and location of near-field dissipation, with q = 0.2 and the bottom-intensified vertical profile described in St. Laurent and Garrett (2002). The remaining 80% of energy dissipation is distributed at one of three horizontal locations--deep basins, continental slope, coastal shelves--with one of three vertical dissipation profiles; dissipation that decays exponentially with height above bottom scales like the buoyancy frequency N or like [N.sup.2] [see Melet et al. (2016) for more detail]. The resulting ocean circulation shows a significant dependence on the vertical profile of dissipation (Figs. 4e,f). In particular, more dissipation in the upper ocean leads to stronger subtropical overturning cells, a broader thermocline, and higher thermosteric sea level; more dissipation in the deep ocean leads to stronger deep meridional overturning circulation [more evidence of these impacts is shown in Melet et al. (2016)]. In addition, the geographic location of the far-field dissipation influences the large-scale circulation notably when it impacts dense-water formation regions: more dissipation on the slopes and shelves near the descending overflows tends to weaken the meridional overturning cell for which the lower branch is supplied by the overflows.

Future work. Future work on the ray-tracing approach should include refinement of the directional spectrum of radiated low-mode waves, including refraction by background flow and evaluation of its impact in GCMs. Further work is also needed to understand and incorporate some of the detailed mechanisms of internal tide dissipation. One of these mechanisms is PSI, which may be especially important near and equatorward of the diurnal turning latitudes ~29[degrees]N/S. Note that the tide energy pathways via the tide constituents [S.sub.2], [O.sub.1], and [K.sub.1], which collectively account for an amount of energy comparable to that of [M.sub.2] (even greater, in some regions), need to be better understood. In particular, internal tides of various frequencies may have different responses to the same bottom topography and time-varying background flow. Progress here will involve a combination of relevant theory and observations with both idealized simulations and realistic tidally forced global simulations. Another dissipation pathway worthy of close attention is wave breaking and turbulence on continental slopes and shelves, where the vertical structure may be heavily influenced by details of wave dynamics in the presence of small-scale coastal topography in ways that are not yet fully understood (e.g., Nash et al. 2007; Kunze et al. 2012; Wain et al. 2013; Pinkel et al. 2015; Waterhouse et al. 2017).

INTERNAL LEE WAVES. Theory and observations. As with tides, mean flows over rough topography can generate internal waves that can remove energy and momentum from the large-scale circulation and, when they break, produce turbulent mixing (Fig. 5a). Quasi-steady flow over small-amplitude bathymetry ([gamma] [??] 1/2; Nikurashin et al. 2014) gives rise to vertically propagating internal lee waves of frequency Uk, where k is the topographic horizontal wavenumber and U is the mean flow speed. For large-amplitude topography ([gamma] [??] 1/2), the Froude number of the flow F = U/NH is O(1), such that topographic flow blocking and splitting becomes prominent; the flow transits the bump generating a nonpropagating disturbance that converts parts of the flow kinetic energy to dissipation. Most of the real ocean lies between these two end cases (Bretherton 1969; Bell 1975; Pierrehumbert and Bacmeister 1987; St. Laurent and Garrett 2002). The drag due to the combination of internal lee-wave generation and topographic flow blocking and splitting is commonly denoted as wave drag in the atmospheric literature. Parameterizations of wave drag have been used for a long time in the atmospheric community (e.g., Palmer et al. 1986) but are less common in the ocean community.

Available global estimates for the energy conversion rate from geostrophic flows into internal lee waves range from 0.2 to 0.75 TW and highlight a prominent role of the Southern Ocean (Bell 1975; Nikurashin and Ferrari 2011; Scott et al. 2011; Wright et al. 2014). Several lines of evidence have suggested the existence of propagating lee waves (e.g., Naveira Garabato et al. 2004; St. Laurent et al. 2012; Waterman et al. 2013; Sheen et al. 2013, 2014; Clement et al. 2016; Fig. 5a). Yet, lee waves have not been definitively identified in ocean observations until recently, with Cusack et al. (2017) reporting unambiguous evidence of a lee wave in the Drake Passage (the search is complicated in part by the difficulty of observing motions with zero Eulerian frequency). Sparse observations also make it difficult to determine the fate of propagating lee waves. Nonpropagating lee waves have been observed in a variety of fracture zones and deep passages (Ferron et al. 1998; Thurnherr et al. 2005; MacKinnon 2013; Alford et al. 2013), but their integrated importance to abyssal mixing is unknown.

Parameterizations and consequences of lee wave-driven mixing on the ocean state. The sensitivity of large-scale ocean circulation to lee wave-driven mixing has been investigated in simulations with the GFDL ESM2G coupled climate model (Melet et al. 2014) using the estimated global map of energy conversion into lee waves of Nikurashin and Ferrari (2011; Fig. 5b). The St. Laurent et al. (2002) exponential vertical structure was used as an initial placeholder for the structure of dissipation associated with breaking lee waves. Although most estimates put the global energy input into lee waves smaller than that into internal tides, Melet et al. (2014) showed that lee wave-driven mixing significantly impacts the ocean state, yielding a reduction of the ocean stratification associated with a warming of the abyssal ocean. The lower cell of the MOC is also slightly lightened and increased in strength (Fig. 5c). The different spatial distribution of the internal tide and lee-wave energy input is largely responsible for the sensitivity described in Melet et al. (2014), highlighting the previously reported importance of the patchiness of internal wave-driven mixing in the ocean (e.g., Simmons et al. 2004a; Jayne 2009; Friedrich et al. 2011). Using a hydrographie climatology and a similar parameterization for lee wave-driven mixing, Nikurashin and Ferrari (2013) and De Lavergne et al. (2016) also show substantial water mass transformation in the Southern Ocean due to internal lee wave-driven mixing.

Trossman et al. (2013, 2016) implemented an inline wave drag parameterization (for both propagating and nonpropagating lee waves) from the atmospheric community (Garner 2005) into a high-resolution ocean general circulation model (Fig. 5d). The inline implementation allows for feedbacks between wave drag and the low-frequency flows that produce the lee waves. They found that the wave drag dissipated a substantial fraction of the wind energy input, significantly reduced both kinetic energy and stratification near the bottom, and reduced the model sea surface height variance and geostrophic surface kinetic energy by measurable amounts of ~20%, while the performance of the model relative to in situ and altimetric measurements of eddy kinetic energy was not negatively impacted. Trossman et al. (2015) showed that dissipations predicted by the Garner (2005) scheme are not inconsistent with microstructure observations within the bottom 500 m in two Southern Ocean regions.

Future work. More observations are needed, especially in the Southern Ocean, to provide definitive evidence of the extent of propagating lee waves in the ocean and further to explore 1) the fraction of local dissipation and the vertical profile of dissipation of the propagating drag, 2) the relative importance of the propagating and nonpropagating lee-wave drag, and 3) the observed mismatch between estimates of lee-wave energy generation and near-bottom dissipation of lee waves.

Enhancing our knowledge of the near-bottom stratification and velocity fields and using a more accurate representation of topographic blocking are crucial for reducing our uncertainty about the global conversion rate into lee waves. Indeed, Wright et al. (2014) found that the use of different stratification products yields a difference of up to 0.25 TW in the global conversion rate into lee waves. Conversion rates are even more sensitive to the near-bottom velocity field (Trossman et al. 2013; Melet et al. 2015), which can vary drastically with model resolution (Thoppil et al. 2011) and should take into account mesoscale eddy velocities. Topographic blocking accounts for most of the predicted dissipation by the Garner (2005) scheme in the bottom 1000 m of two Southern Ocean domains (Trossman et al. 2015). Recent laboratory experiments by Dossmann et al. (2016) have shown that, for most forcing parameters they considered, nonlinear mixing mechanisms close to abyssal topography, such as topographic blocking, dominate the remote mixing mechanism by lee waves. Yet, theoretical conversion rates are highly sensitive to the choice of uncertain parameters related to the representation of topographic blocking and splitting (Nikurashin et al. 2014).

As parameterized lee-wave drag makes a significant impact on the ocean state (Trossman et al. 2013, 2016), it should be included inline within climate models in a dynamically accurate manner to ensure credible ocean representation in a changing climate. Using linear theory and modeled resolved and parameterized bottom velocities and stratification, Melet et al. (2015) showed that the energy flux into lee waves exhibits a clear annual cycle in the Southern Ocean and that the global energy flux is projected to decrease by ~20% from preindustrial to future climate conditions under the representative concentration pathway 8.5 (RCP8.5) scenario. This time variability is primarily due to changes in bottom velocities (Melet et al. 2015). Ultimately, models should aspire to a full coupling between wind power, eddies and geostrophic circulations, stratification, and lee-wave drag and induced mixing. Such a coupling requires a state-dependent, time-evolving parameterization for the effects of lee waves.

WIND-DRIVEN NEAR-INERTIAL MOTIONS. Theory and observations. Much of what is known about wind-generated near-inertial waves (NIWs) builds on the observations and model studies of the Ocean Storms Experiment (D'Asaro et al. 1995; Dohan and Davis 2011); for a summary of the outcomes, other generation mechanisms, and additional studies, see a review by Alford et al. (2016). Inertial oscillations of the boundary layer are a free mode of the ocean and are its first response to changes in the wind stress (e.g., DAsaro 1985). Part of the inertial oscillation energy is dissipated in the boundary layer through shear instability, thus converting kinetic energy to heat and potential energy (Large and Crawford 1995), with the remainder radiated away downward (Fig. 6a) and equatorward (Fig. 6b) in the form of propagating near-inertial internal waves (Alford 2003a; Plueddemann and Farrar 2006; Alford et al. 2012; Simmons and Alford 2012). The partition between high and low modes and the energy lost to dissipation at the mixed-layer base is unknown. In the Ocean Storms Experiment, approximately one-third of the energy input by the wind was carried away equatorward in modes 1 and 2. Another study (Alford et al. 2012) found a similar fraction was carried downward in higher modes, while a modeling study by Furuichi et al. (2008) found that only 10% reached past 150 m. Inferred global upper-ocean dissipation rates show a clear seasonal cycle (Whalen et al. 2012), particularly in storm-track latitudes (Whalen et al. 2015). Near-inertial KE at all depths also shows a clear seasonal cycle, indicating that some of the energy makes it deep into the ocean (Alford and Whitmont 2007; Silverthorne and Toole 2009).

Parameterizations and consequences. The CPT tackled the upper-ocean portion of the NIW-related mixing with a three-step process, described in Jochum et al. (2013), suitable for general use in coupled atmosphere-ocean models. First, atmosphere and ocean models are coupled more frequently (e.g., 2 h instead of daily) to allow resonant generation of near-inertial motions in the oceanic surface boundary layer. Even with high-frequency coupling, the near-inertial speeds can be too weak by 50% if the frontal structure of storms is not properly resolved by the atmospheric component of climate models. In such cases, the missing amplitude of the NIWs must be computed during the integration and added to the shear calculation of the boundary layer parameterization. The online computation of the near-inertial part of the velocity is not trivial because during the integration the ocean model only has information about adjacent time steps. Fortunately, however, outside the deep tropics, velocity fluctuations from one model time step (e.g., 1 h) to the next are mostly due to NIWs, which allow the accurate determination of near-inertial velocity during the integration [see Jochum et al. (2013) for details and method verification], Last, the air-sea flux of inertial wave energy into the boundary layer is determined, and 30% of it (Rimac et al. 2016) is used to increase the background diffusivity below the boundary layer. The energy in the last step is distributed with an exponential decay scale of 2000 m (Alford and Whitmont 2007). The resultant turbulent mixing from near-inertial motions changes the heat distribution in the upper ocean significantly enough to modify tropical SST patterns and leads to a 20% reduction in tropical precipitation biases (Jochum et al. 2013; for the sensitivity of precipitation to the strength of near-inertial waves, see Figs. 6c and 6d).

Ongoing and future work. Much hinges on the appropriate representation of NIWs. The largest uncertainties are associated with the poorly known high-frequency and high-wavenumber part of the wind spectrum and the partitioning between locally dissipated energy and the amount radiated away. Thus, the energy available for NIW-induced mixing in the surface boundary layer ranges from 0.3 to 1.0 TW (Alford 2001, 2003b; Simmons and Alford 2012; Rimac et al. 2013). The Jochum et al. (2013) study was based on 0.34 TW; allowing for 0.68 TW in the Community Climate System Model would remove the spurious southern intertropical convergence zone (ITCZ) and would result in a realistically shaped South Pacific convergence zone (Fig. 6c). Thus, ongoing work focuses on the detailed analysis of moorings with collocated wind and ocean velocity measurements (e.g., Plueddemann and Farrar 2006; Alford et al. 2012).

TOOLS AND TECHNIQUES. Microstructure database. The CPT worked in conjunction with the CLIVAR and Carbon Hydrographie Data Office (CCHDO) at Scripps Institution of Oceanography to develop a standardized format for archiving microstructure data. Data have been archived as climate and forecast (CF)-compliant network Common Data Form (netCDF) files with 1-m binned data (where possible). The database contains the following variables: time, depth, pressure, temperature, salinity, latitude, longitude, and bottom depth. The database also contains the newly designated variables: epsilon (W kg-1; ocean turbulent kinetic energy dissipation rate) and, when available, chi-t ([degrees][C.sup.2] [s.sup.-1]; ocean dissipation rate of thermal variance from microtemperature) and chi-c ([degrees][C.sup.-2] [s.sup.-1]; ocean dissipation rate of thermal variance from microconductivity). Database entries include names of the project, project principal investigators (Pis), and cruise information (research ship, ports of entry and exit, cruise dates, and chief scientist). Database entries have project-specific DOIs to cite the data in publications. Relevant cruise reports, project-related papers, and other documents are also contained in the data archive. At present, the database consists of 25 separate projects and can be accessed online (at http://microstructure.ucsd.edu). Newly obtained microstructure data can be uploaded to the microstructure database by sending 1-m binned data to the CCHDO office (at http://cchdo.ucsd.edu/submit).

A repository for ocean mixing analysis tools, methods, and code. The availability of commercially manufactured turbulence profilers, along with an increased use of mixing proxies, have expanded the size of the mixing community and the number of publications that use mixing observations. Many variants of processing code have thus been developed in parallel by different groups. Some variants have subtle differences in methodology that can potentially lead to significant quantitative differences in the results. We thus sought to establish a community-based online repository for best-practice data analysis tools used for ocean mixing and internal wave calculations. Analysis code from many independent groups is available for download from the repository, thus facilitating comparison of techniques in an open, objective way. To accomplish this goal, a Github mixing repository was created (https://github.com /OceanMixingCommunity/) and populated with standard algorithms and process methods.

The goals of the public repository are to 1) enable reproducibility of analyses, 2) allow for comparison of different datasets using the same code, 3) provide a means for easy reanalysis if a bug is identified or a best-practice change is suggested, 4) allow testing of one code against another version, and 5) provide a well-documented and version-controlled repository suitable for citation of techniques employed in publications. The code is primarily (but not exclusively) MATLAB based and includes routines for calculation of Thorpe scales, [N.sup.2], finescale parameterizations, generic and instrument-specific turbulence processing code, and sample data files.

Observational data analysis: The finescale parameterizations. Many of the insights described in this paper were inspired in part by the vast expansion of mixing data (e.g., Fig. 2) that has come from widespread use of the finescale parameterization for ocean mixing rates. Its increasing popularity warrants a few comments here. Finescale parameterizations produce the average dissipation rate expected over several wave periods and therefore are helpful in assessing the spatial- and temporal-mean dissipation rate or diffusivity. Inferences of mixing from finescale parameterizations are more extensive than instantaneous observations of turbulence from microstructure measurements (e.g., Polzin et al. 1996; Kunze et al. 2006; Whalen et al. 2012).

Finescale parameterizations rely on the fact that the observed shear and strain variance in the thermocline and below is mainly caused by internal waves. The parameterizations also assume that the energy dissipation rate is primarily due to nonlinear interactions between internal waves that transfer energy from the finescale toward smaller-scale waves that subsequently break into turbulence. As discussed in Polzin et al. (2014), an expression of the downspectrum energy cascade in the open ocean has been developed (Henyey et al. 1986; Muller et al. 1986; Henyey and Pomphrey 1983) in terms of the shear and strain spectra. This expression allows for estimates of the dissipation rate as a function of the spectra.

Parameterizations using finescale shear and strain profiles have been tested in a variety of contexts, consistently demonstrating a factor-of-2-3 agreement with microstructure inferences in open-ocean conditions (Gregg 1989; Polzin et al. 1995; Winkel et al. 2002; Polzin et al. 2014) and with strain-only inferences in a variety of locations (Wijesekera et al. 1993; Frants et al. 2013; Waterman et al. 2014; Whalen et al. 2015). The shear- and strain-based parameterization is known to be less effective in regions where the underlying assumptions behind the parameterization do not apply (Polzin et al. 2014). These regions include continental shelves (Mackinnon and Gregg 2003), strong geostrophic flow regimes over rough topography (Waterman et al. 2014), and regions with very large overturning internal waves (Klymak et al. 2008). Implementation of the parameterizations in the open-ocean have revealed reasonable patterns and insight into the geography of diapycnal mixing using shear (Polzin et al. 1997; Kunze et al. 2006; Huussen et al. 2012) and strain (Polzin et al. 1997; Kunze et al. 2006; Huussen et al. 2012) and strain (Kunze et al. 2006; Wu et al. 2011; Whalen et al. 2012). A global dissipation rate product that is based on both finestructure estimates and microstructure measurements is currently in preparation that will be made publicly available (C. Whalen 2017, unpublished manuscript).

Global internal wave models. It has only been in the last decade that global models of internal waves have been developed (Arbic et al. 2004; Simmons et al. 2004a). As described above, several global internal wave models used in the community now include atmospheric and tidal forcing, enabling examination of many issues of interest such as the global three-dimensional internal wave geography, internal wave-mesoscale interactions, and an internal gravity wave continuum spectrum that approaches the observed continuum more closely as model resolution is refined (Muller et al. 2015).

The CVMix package. The Community Ocean Vertical Mixing (CVMix) package is a software package that provides transparent, robust, flexible, well-documented, and shared FORTRAN source codes for use in parameterizing vertical mixing processes in numerical ocean models. The project is focused on developing software for a consensus of first-order closures that return a vertical diffusivity, viscosity, and possibly a nonlocal transport [e.g., as in the K-profile parameterization (KPP) scheme of Large et al. (1994)], with each quantity dependent on the tracer or velocity being mixed. CVMix provides a software framework for the physical parameterizations arising from the internal wave-driven mixing CPT. For example, the Simmons et al. (2004b) tidal mixing scheme, available in CVMix, serves as a useful example for other tidal mixing schemes such as Melet et al. (2013a). Code development occurs within a community of scientists and engineers who make use of CVMix modules for a variety of ocean climate models [e.g., Model for Prediction across Scales Ocean (MPAS-O) used at Los Alamos National Laboratory, Parallel Ocean Program (POP) used at NCAR, and MOM6 used at GFDL], CVMix modules are freely available to the community under the GNU General Public License, version 2 (GPLv2), using an open development approach on Github (https://github .com/CVMix). We solicit further contributions of parameterizations, thus enabling a very broad group of climate modelers to make use of the schemes.

SUMMARY AND FUTURE SCIENCE DIRECTIONS. A frequently asked question related to this work is "Which mixing processes matter most for climate?" As with many alluringly comprehensive sounding questions, the answer is "it depends." Deep-ocean mixing matters for the decadal to centennial time scales on which the deep global circulation evolves. The mixing process most important for the deep circulation is the one with the most power, namely, the tides. The distribution of mixing above deep rough topography from near-field tidal dissipation is the most fully developed aspect of our work, both in terms of dynamical understanding and parameterization implementation (see the section on "Near-field tidal mixing" Fig. 3). As detailed in the section titled "Far-field internal tides," our understanding of far-field tidal dissipation is less complete. Lee waves may also contain significant power and play an important role in places like the Southern Ocean; preliminary results hint at a substantial role in water mass modification in this globally important region, but more observations and data-model-theory comparison is needed before we are confident of how best to represent them (see the section on "Internal lee waves"; Fig. 5). Nonpropagating form drag is known to be important for momentum budgets in the atmosphere but has just begun to receive significant oceanographic attention (Trossman et al. 2016); it may be not only locally important for mixing tracers and momentum wherever strong flow encounters sharp or rough topography but a globally important drain of mesoscale energy.

Mixing in the main pycnocline can impact heat distribution and steric sea level rise on decadal time scales, which makes it a compelling societal problem. Turbulent mixing in this depth range is controlled by a combination of downward-propagating, near-inertial waves (see the section on "Wind-driven near-inertial motions"; Fig. 6); low-mode, long-range-propagating internal tides breaking on continental slopes (see the section on "Far-field internal tides"; Fig. 4); and by near-field breaking of upward-propagating internal tides or lee waves through nonlinear interactions. Double-diffusion processes may also be significant in the main pycnocline (e.g., Schmitt et al. 2005) but are not covered here. For forward progress, a better understanding of low-mode wave breaking on slopes, with particular focus on the vertical structure of resultant dissipation (Carter and Gregg 2002; Nash et al. 2004, 2007; Martini et al. 2011; Kunze et al. 2012; Pinkel et al. 2015; Waterhouse et al. 2017), will help to constrain mixing rates.

It is increasingly clear that near-inertial-wave-driven mixing both below the surface boundary layer and down into the main thermocline is significantly mediated by the presence of mesoscale eddies. Areas of enhanced diffusivities have been linked to regions of elevated eddy kinetic energy, though the mechanisms are not always clear (e.g., Kunze et al. 1995; Whalen et al. 2012). In turn, interactions with internal waves may be a significant energy loss term for eddies (Buhler and Mclntyre 2005; Polzin 2010; Whalen 2015; Barkan et al. 2017).

Mixing in the upper ocean matters to climate phenomena of seasonal to interannual, and perhaps even longer, time scales. Turbulence beneath the surface boundary layer has a strong effect on upper-ocean freshwater content and heat and, through SST changes, on a variety of coupled air-sea interactions ranging from the MJO to ENSO (e.g., Moum et al. 2016). In this depth range (of order 100 m below the boundary layer), turbulence from breaking NIW plays a dominant role (see the section on "Wind-driven near-inertial motions"; Fig. 6). Again, the interaction with mesoscale eddies, and in particular mesoscale vorticity, may play a large role in setting the patterns and rates of wave propagation and dissipation in ways that are poorly constrained. We hope that continued work in this field will be closely coupled with the many active research programs focused on mixing parameterizations within the surface boundary layer, which may also be ripe for a CPT-style renovation.

Upper-ocean mixing takes on a unique relevance at high latitudes. The presence of ice (either ice shelves or sea ice) significantly changes both the dynamics and thermodynamics of turbulence near the poles, particularly in the near-surface ocean. Yet accurate representation of mixing in these environments is crucial if we are to accurately forecast everything from ice melt rates to high-latitude C[O.sub.2] absorption/ outgassing, to deep-water formation, to ecosystem responses to climate change. Multiple U.S. funding agencies are increasingly putting substantial resources into process studies, long-term observations, and modeling. A formalized CPT-like framework might help bring these components together.

BEST PRACTICES FOR CONTINUING SUCCESS. Once a field is in a state of readiness, where substantial observations, theory, and dynamical understanding exist, the Climate Process Team structure or similar programs provide a productive template for progress. The CPT framework allows for 1) motivation to bring some parts of that research to a state of closure and 2) the opportunity to bring together observationalists, theorists, and modelers to work through details of synthesizing observational reality, theoretical insights, and modeling efforts. The formal charge of CPT funding was essential to initiate this process and sustain it for the years necessary to bring such collaboration to productive fruition. A crucial component of this successful interaction has been the presence of dedicated personnel who pull together the state of observational science and/or are embedded within modeling centers; postdocs or early career scientists fit well into this role. Similar facilitated cross-field collaborations are increasingly built into the structure of other multi-PI projects, best practices for which are well described by Cronin et al. (2009). At the same time, the epiphanies, new ideas, and novel observations that fundamentally drive the field forward frequently come not from big science but from a cornucopia of much smaller exploratory efforts and the continued small-scale development of innovative observing technology and numerical techniques. We must not lose the ability to be surprised.

AFFILIATIONS: MacKinnon, Waterhouse, Pinkel, Merchant, Diggs, Barna, and Alford--Scripps Institution of Oceanography, La Jolla, California; Zhao and Whalen--Applied Physics Laboratory, University of Washington, Seattle, Washington; Trossman--Goddard Earth Sciences Technology and Research, Greenbelt, and Department of Earth and Planetary Sciences, Johns Hopkins University, Baltimore, Maryland; Sun, St. Laurent, Polzin, and Jayne--Woods Hole Oceanographie Institution, Woods Hole, Massachusetts; Simmons--University of Alaska Fairbanks, Fairbanks, Alaska; Pickering and Nash--Oregon State University, Corvallis, Oregon; Norton, Large, Danabasoglu, Bryan, and Briegleb--National Center for Atmospheric Research, * Boulder, Colorado; Musgrave--Massachusetts Institute of Technology, Cambridge, Massachusetts; Melet--Program in Atmospheric and Oceanic Sciences, Princeton University, Princeton, New Jersey, and Mercator Ocean, Ramonville St. Agne, France; Mater and Legg--Program in Atmospheric and Oceanic Sciences, Princeton University, Princeton, New Jersey; Kunze--Northwest Research Associates, Seattle, Washington; Klymak--University of Victoria, Victoria, British Columbia, Canada; Jochum--Niels Bohr Institute, Copenhagen, Denmark; Hallberg and Griffies--NOAA/Geophysical Fluid Dynamics Laboratory, Princeton, New Jersey; Chassignet--Center for Ocean-Atmospheric Prediction Studies, Florida State University, Tallahassee, Florida; Buijsman--The University of Southern Mississippi, Hattiesburg, Mississippi, and Division of Marine Science, John C. Stennis Space Center, Hancock County, Mississippi; Arbic and Ansong--Department of Earth and Environmental Sciences, University of Michigan, Ann Arbor, Michigan

* The National Center for Atmospheric Research is sponsored by the National Science Foundation.

CURRENT AFFILIATIONS: Musgrave--Woods Hole Oceanographic Intitution, Woods Hole, Massachusetts; Melet--Mercator Ocean, Ramonville Saint-Agne, France

CORRESPONDING AUTHOR: Jennifer A. MacKinnon, jmackinnon@ucsd.edu

The abstract for this article can be found in this issue, following the table of contents.

DOI: 10.1175/BAMS-D-16-0030.1

In final form 14 March 2017

ACKNOWLEDGMENTS. We are grateful to U.S. CLIVAR for their leadership in instigating and facilitating the Climate Process Team program. We are indebted to NSF and NOAA for sponsoring the CPT series, and Eric Itsweire in particular for his steadfast support and enthusiasm. We thank Alistair J. Adcroft, Mike Levy, Brandon Reichl, Todd Ringler, and Luke Van Roekel for their contributions to the CVMix project; Peter Gent for broad insight; and Andreas Schmittner and David Ullman for their contributions to the advances in tidal mixing parameterizations in the CESM ocean component.

REFERENCES

Alford, M. H., 2001: Internal swell generation: The spatial distribution of energy flux from the wind to mixed layer near-inertial motions. J. Phys. Oceanogr., 31, 2359-2368, doi:10.1175/1520-0485(2001)031<2359: ISGTSD>2,O.CO;2.

--, 2003a: Energy available for ocean mixing redistributed through long-range propagation of internal waves. Nature, 423, 159-163, doi:10.1038 /nature01628.

--, 2003b: Improved global maps and 54-year history of wind-work on ocean inertial motions. Geophys. Res. Lett., 30,1424-1427, doi:10.1029/2002GL016614.

--, 2008: Observations of parametric subharmonic instability of the diurnal internal tide in the South China Sea. Geophys. Res. Lett., 35, L15602, doi:10.1029 /2008GL034720.

--, and M. Whitmont, 2007: Seasonal and spatial variability of near-inertial kinetic energy from historical moored velocity records. J. Phys. Oceanogr., 37, 2022-2037, doi:10.1175/JP03106.1.

--, J. A. MacKinnon, Z. Zhao, R. Pinkel, J. Klymak, and T. Peacock, 2007: Internal waves across the Pacific. Geophys. Res. Lett., 34, L24601, doi:10.1029 /2007GL031566.

--, M. F. Cronin, and J. M. Klymak, 2012: Annual cycle and depth penetration of wind-generated near-inertial internal waves at Ocean Station Papa in the northeast Pacific. J. Phys. Oceanogr., 42, 889-909, doi:10 .1175/JPO-D-11-092.1.

--, J. B. Girton, G. Voet, G. S. Carter, J. B. Mickett, and J. M. Klymak, 2013: Turbulent mixing and hydraulic control of abyssal water in the Samoan Passage. Geophys. Res. Lett., 40, 4668-4674, doi:10.1002/grl .50684.

--, J. M. Klymak, and G. S. Carter, 2014: Breaking internal lee waves at Kaena Ridge, Hawaii. Geophys. Res. Lett., 41, 906-912, doi:10.1002/2013GL059070.

--, and Coauthors, 2015: The formation and fate of internal waves in the South China Sea. Nature, 521, 65-69, doi:10.1038/nature14399.

--, J. A. MacKinnon, H. L. Simmons, and J. D. Nash, 2016: Near-inertial internal gravity waves in the ocean. Annu. Rev. Mar. Sci., 8, 95-123, doi:10.1146 /annurev-marine-010814-015746.

Ansong, J. K, B. K. Arbic, M. C. Buijsman, J. G. Richman, J. F. Shriver, and A. J. Wallcraft, 2015: Indirect evidence for substantial damping of low-mode internal tides in the open ocean. J. Geophys. Res. Oceans, 120, 6057-6071, doi:10.1002/2015JC010998.

--, and Coauthors, 2017: Semidiurnal internal tide energy fluxes and their variability in a global ocean model and moored observations. J. Geophys. Res. Oceans, 122, 1882-1900, doi:10.1002/2016JC012184.

Arbic, B. K., S. T. Garner, R. Hallberg, and H. L. Simmons, 2004: The accuracy of surface elevations in forward global barotropic and baroclinic tide models. Deep-Sea Res. II, 51, 3069-3101, doi:10.1016/j.dsr2 .2004.09.014.

--, A. J. Wallcraft, and E. J. Metzger, 2010: Concurrent simulation of the eddying general circulation and tides in a global ocean model. Ocean Modell., 32, 175-187, doi:10.1016/j.ocemod.2010.01.007.

Balmforth, N. J., and T. Peacock, 2009: Tidal conversion by supercritical topography. J. Phys. Oceanogr., 39, 1965-1974, doi:10.1175/2009JPO4057.1.

Barkan, R., K. B. Winters, and J. C. McWilliams, 2017: Stimulated imbalance and the enhancement of eddy kinetic energy dissipation by internal waves. J. Phys. Oceanogr., 47, 181-198, doi:10.1175/JPO-D -16-0117.1.

Bell, T., 1975: Topographically generated internal waves in the open ocean. J. Geophys. Res., 80, 320-327, doi:10.1029/JC080i003p00320.

Bretherton, F. P., 1969: Momentum transport by gravity waves. Quari. J. Roy. Meteor. Soc., 95, 213-243, doi:10.1002/qj.49709540402.

Bryan, K., and L. J. Lewis, 1979: A water mass model of the World Ocean. J. Geophys. Res., 84, 2503-2517, doi:10.1029/JC084iC05p02503.

Buhler, O., and M. Mclntyre, 2005: Wave capture and wave-vortex duality. J. Fluid Mech., 534, 67-96, doi:10.1017/S0022112005004374.

--, and M. Holmes-Cerfon, 2011: Decay of an internal tide due to random topography in the ocean. J. Fluid Mech., 678, 271-293, doi:10.1017/jfm.2011.115.

Buijsman, M. C., S. Legg, and J. Klymak, 2012: Double ridge internal tide interference and its effect on dissipation in Luzon Strait. J. Phys. Oceanogr., 42, 1337-1356, doi: 10.1175/JPO-D-11-0210.1.

--, and Coauthors, 2014: Three-dimensional double ridge internal tide resonance in Luzon Strait. J. Phys. Oceanogr., 44,850-869, doi:10.1175/JPO-D-13 -024.1.

--, and Coauthors, 2016: Impact of internal wave drag on the semidiurnal energy balance in a global ocean circulation model. J. Phys. Oceanogr., 46, 1399-1419, doi:10.1175/JPO-D-15-0074.1.

Cacchione, D., and C. Wunsch, 1974: Experimental study of internal waves over a slope. J. Fluid Mech., 66, 223-239, doi:10.1017/S0022112074000164.

Carter, G. S., and M. C. Gregg, 2002: Intense, variable mixing near the head of Monterey Submarine Canyon. J. Phys. Oceanogr., 32, 3145-3165, doi:10.1175 /1520-0485(2002)032<3145:IVMNTH>2,O.CO;2.

Clement, L., E. Frajka-Williams, K. L. Sheen, J. A. Brearley, and A. C. Naveira Garabato, 2016: Generation of internal waves by eddies impinging on the western boundary of the North Atlantic. J. Phys. Oceanogr., 46, 1067-1079, doi:10.1175/JPO-D-14 -0241.1.

Cronin, M. F., S. Legg, and P. Zuidema, 2009: Climate research: Best practices for process studies. Bull. Amer. Meteor. Soc., 90, 917-918, doi:10.1175 /2009BAMS2622.1.

Cusack, J. M., A. C. Naveira Garabato, D. A. Smeed, and J. B. Girton, 2017: Observation of a large lee wave in the Drake Passage. J. Phys. Oceanogr., 47, 793-810, doi:10.1175/JPO-D-16-0153.1.

Danabasoglu, G., W. Large, and B. Briegleb, 2010: Climate impacts of parameterized Nordic Sea overflows. J. Geophys. Res., 115, C11005, doi:10.1029 /2010JC006243.

--, S. Bates, B. Briegleb, S. Jayne, M. Jochum, W. Large, S. Peacock, and S. Yeager, 2012: The CCSM4 ocean component. J. Climate, 25, 1361-1389, doi:10.1175 /JCLI-D-11-00091.1.

D'Asaro, E. A., 1985: The energy flux from the wind to near-inertial motions in the mixed layer. J. Phys. Oceanogr., 15, 943-959, doi:10.1175/1520 -0485(1985)015<0943:UOTSIC>2.0.CO;2.

--, C. E. Eriksen, M. D. Levine, P. Niiler, C. A. Paulson, and P. V. Meurs, 1995: Upper-ocean inertial currents forced by a strong storm. Part I: Data and comparisons with linear theory. J. Phys. Oceanogr., 25, 2909-2936, doi:10.1175/1520-0485(1995)025<2909: UOICFB>2-O.CO;2.

Decloedt, T., and D. Luther, 2010: On a simple empirical parameterization of topography-catalyzed diapycnal mixing in the abyssal ocean. J. Phys. Oceanogr., 40, 487-508, doi:10.1175/2009JP04275.1.

De Lavergne, C., G. Madec, J. Le Sommer, A. G. Nurser, and A. C. Naveira Garabato, 2016: The impact of a variable mixing efficiency on the abyssal overturning. J. Phys. Oceanogr., 46, 663-681, doi:10.1175 /JPO-D-14-0259.1.

Dohan, K., and R. E. Davis, 2011: Mixing in the transition layer during two storm events. J. Phys. Oceanogr., 41, 42-66, doi:10.1175/2010JP04253.1.

Dossmann, Y., M. G. Rosevear, R. W. Griffiths, A. McC. Hogg, G. O. Hughes, and M. Copeland, 2016: Experiments with mixing in stratified flow over a topographic ridge. J. Geophys. Res. Oceans, 121, 6961-6977, doi:10.1002/2016JC011990.

Dunne, J. P., and Coauthors, 2012: GFDL's ESM2 global coupled climate-carbon Earth System Models. Part I: Physical formulation and baseline simulation characteristics. J. Climate, 25, 6646-6665, doi: 10.1175 /JCLI-D-11-00560.1.

Dunphy, M., and K. G. Lamb, 2014: Focusing and vertical mode scattering of the first mode internal tide by mesoscale eddy interaction. J. Geophys. Res. Oceans, 119, 523-536, doi:10.1002/2013JC009293.

Dushaw, B., B. Howe, B. Cornuelle, P. Worcester, and D. Luther, 1995: Barotropic and baroclinic tides in the central North Pacific Ocean determined from long-range reciprocal acoustic transmissions. J. Phys. Oceanogr., 25, 631-647, doi:10.1175/1520 -0485(1995)025<0631:BABTIT>2.0.CO;2.

Eden, C., and D. Olbers, 2014: An energy compartment model for propagation, nonlinear interaction, and dissipation of internal gravity waves. J. Phys. Oceanogr., 44, 2093-2106, doi:10.1175/JPO -D-13-0224.1.

Egbert, G. D., and R. D. Ray, 2003: Semi-diurnal and diurnal tidal dissipation from TOPEX/ Poseidon altimetry. Geophys. Res. Lett., 30, 1907, doi:10.1029/2003GL017676.

Ferrari, R., and C. Wunsch, 2009: Ocean circulation kinetic energy: Reservoirs, sources, and sinks. Annu. Rev. Fluid Mech., 41, 253-282, doi:10.1146/annurev .fluid.40.111406.102139.

Ferron, B. H., H. Mercier, K. Speer, A. Gargett, and K. Polzin, 1998: Mixing in the Romanche Fracture Zone. J. Phys. Oceanogr., 28, 1929-1945, doi:10.1175/1520 -0485(1998)028<1929:MITRFZ>2.0.CO;2.

Frants, M., G. M. Damerell, S. T. Gille, K. J. Heywood, J. A. Mackinnon, and J. Sprintall, 2013: An assessment of density-based finescale methods for estimating diapycnal diffusivity in the Southern Ocean. J. Atmos. Oceanic Technol., 30, 2647-2661, doi:10.1175 /JTECH-D-12-00241.1.

Friedrich, T., A. Timmermann, T. Decloedt, D. S. Luther, and A. Mouchet, 2011: The effect of topography-enhanced diapycnal mixing on ocean and atmospheric circulation and marine biogeochemistry. Ocean Modell., 39, 262-274, doi:10.1016/j .ocemod.2011.04.012.

Furuichi, N., T. Hibiya, and Y. Niwa, 2008: Model-predicted distribution of wind-induced internal wave energy in the world's oceans. J. Geophys. Res., 113, C09034, doi: 10.1029/2008JC004768.

Gargett, A. E., 1984: Vertical eddy diffusivity in the ocean interior. J. Mar. Res., 42, 359-393, doi:10.1357 /002224084788502756.

Garner, S. T., 2005: A topographic drag closure built on an analytical base flux. J. Atmos. Sci., 62, 2302-2315, doi:10.1175/JAS3496.1.

Garrett, C., and E. Kunze, 2007: Internal tide generation in the deep ocean. Annu. Rev. Fluid Mech., 39, 57-87, doi:10.1146/annurev.fluid.39.050905.110227.

Gaspar, P., Y. Gregoris, and J. Lefevre, 1990: A simple eddy kinetic energy model for simulations of the oceanic vertical mixing: Tests at station Papa and long-term upper ocean study site. J. Geophys. Res., 95, 16 179-16 193, doi:10.1029/JC095iC09p16179.

Green, J. A. M., and J. Nycander, 2013: A comparison of tidal conversion parameterizations for tidal models. J. Phys. Oceanogr., 43,104-119, doi:10.1175/JPO-D-12 -023.1.

Gregg, M. C., 1989: Scaling turbulent dissipation in the thermocline. J. Geophys. Res., 94, 9686-9698, doi:10.1029/JC094iC07p09686.

--, E. A. D'Asaro, and J. J. Riley, 2018: Mixing coefficients and mixing efficiency in the ocean. Annu. Rev. Mar. Sci., doi:10.1146/annurev-marine -121916-063643, in press.

Griffies, S. M., R. C. Pacanowski, and R. W. Hallberg, 2000: Spurious diapycnal mixing associated with advection in a z-coordinate ocean model. Mon. Wea. Rev., 128, 538-564, doi:10.1175/1520 -0493(2000)128<0538:SDMAWA>2.0.CO;2.

Hall, R., J. Huthnance, and R. Williams, 2013: Internal wave reflection on shelf slopes with depth-varying stratification. J. Phys. Oceanogr., 43,248-258, doi:10.1175 /JPO-D-11-0192.1.

Harrison, M., and R. Hallberg, 2008: Pacific subtropical cell response to reduced equatorial dissipation. J. Phys. Oceanogr., 38, 1894-1912, doi:10.1175/2008JP03708.1.

Hasselmann, S., and Coauthors, 1988: The WAM model--A third generation ocean wave prediction model. J. Phys. Oceanogr., 18, 1775-1810, doi:10.1175/15200485(1988)018<1775:TWMTGO>2.0.CO;2.

Hazewinkel, J., and K. B. Winters, 2011: PSI of the internal tide on a plane: Flux divergence and near inertial wave propagation. J. Phys. Oceanogr., 41, 1673-1682, doi:10.1175/2011JP04605.1.

Helfrich, K. R., and R. H. J. Grimshaw, 2008: Nonlinear disintegration of the internal tide. J. Phys. Oceanogr., 38, 686-701, doi:10.1175/2007JP03826.1.

Henyey, F. S., and N. Pomphrey, 1983: Eikonal description of internal wave interactions: A non-diffusive picture of "induced diffusion." Dyn. Atmos. Oceans, 7, 189-219, doi:10.1016/0377-0265(83)90005-2.

--, J. Wright, and S. M. Flatte, 1986: Energy and action flow through the internal wave field. J. Geophys. Res., 91, 8487-8495, doi:10.1029/JC091iC07p08487.

Huussen, T. N., A. C. Naveira-Garabato, H. L. Bryden, and E. L. McDonagh, 2012: Is the deep Indian Ocean MOC sustained by breaking internal waves? J. Geophys. Res., 117, C08024, doi:10.1029 /2012JC008236.

Ilicak, M., A. Adcroft, S. Griffies, and R. Hallberg, 2012: Spurious dianeutral mixing and the role of momentum closure. Ocean Modell., 45-46, 37-58, doi:10.1016/j.ocemod.2011.10.003.

Ivey, G., and R. Nokes, 1989: Vertical mixing due to the breaking of critical internal waves on sloping boundaries. J. Fluid Mech., 204,479-500, doi:10.1017 /S0022112089001849.

--, K. Winters, and I. de Silva, 2000: Turbulent mixing in a sloping benthic boundary layer energized by internal waves. J. Fluid Mech., 418,59-76, doi:10.1017 /S0022112000008788.

--, --, and J. Koseff, 2008: Density stratification, turbulence, but how much mixing? Annu. Rev. Fluid Mech., 40, 169-184, doi:10.1146/annurev.fluid .39.050905.110314.

Jackson, L., R. Hallberg, and S. Legg, 2008: A parameterization of shear-driven turbulence for ocean climate models. J. Phys. Oceanogr., 38, 1033-1053, doi:10.1175/2007JP03779.1.

Jayne, S. R., 2009: The impact of abyssal mixing parameterizations in an ocean general circulation model. J. Phys. Oceanogr., 39,1756-1775, doi:10.1175 /2009JP04085.1.

--, and L. C. St. Laurent, 2001: Parameterizing tidal dissipation over rough topography. Geophys. Res. Lett., 28, 811-814, doi:10.1029/2000GL012044.

Jochum, M., 2009: Impact of latitudinal variations in vertical diffusivity on climate simulations. J. Geophys. Res., 114, C01010, doi:10.1029 /2008JC005030.

--, B. P. Briegleb, G. Danabasoglu, W. G. Large, N. J. Norton, S. R. Jayne, M. H. Alford, and F. O. Bryan, 2013: The impact of oceanic near-inertial waves on climate. J. Climate, 26, 2833-2844, doi:10.1175/JCLI -D-12-00181.1.

Johnston, T. M. S., and M. A. Merrifield, 2003: Internal tide scattering at seamounts, ridges, and islands. J. Geophys. Res., 108,3180, doi:10.1029/2002JC001528.

--, D. L. Rudnick, and S. M. Kelly, 2015: Standing internal tides in the Tasman Sea observed by gliders. J. Phys. Oceanogr., 45, 2715-2737, doi:10.1175/JPO -D-15-0038.1.

Kelly, S. M., J. D. Nash, K. I. Martini, M. H. Alford, and E. Kunze, 2012: The cascade of tidal energy from low to high modes on a continental slope. J. Phys. Oceanogr., 42, 1217-1232, doi:10.1175/JPO -D-ll-0231.1.

--, N. L. Jones, J. D. Nash, and A. F. Waterhouse, 2013: The geography of semidiurnal mode-1 internal-tide energy loss. Geophys. Res. Lett., 40, 4689-4693, doi:10.1002/grl.50872.

Kerry, C., B. Powell, and G. Carter, 2014: The impact of subtidal circulation on internal-tide-induced mixing in the Philippine Sea. J. Phys. Oceanogr., 44, 3209-3224, doi:10.1175/JPO-D-13-0249.1.

Klinger, B. A., J. Marshall, and U. Send, 1996: Representation of convective plumes by vertical adjustment. /. Geophys. Res., 101, 18 175-18 182, doi:10.1029 /96JC00861.

Klymak, J. M., and Coauthors, 2006: An estimate of tidal energy lost to turbulence at the Hawaiian Ridge. J. Phys. Oceanogr., 36, 1148-1164, doi:10.1175 /JP02885.1.

--, R. Pinkel, and L. Rainville, 2008: Direct breaking of the internal tide near topography: Kaena Ridge, Hawaii. J. Phys. Oceanogr., 38, 380-399, doi:10.1175 /2007JP03728.1.

--, S. Legg, and R. Pinkel, 2010: A simple parameterization of turbulent tidal mixing near supercritical topography. J. Phys. Oceanogr., 40, 2059-2074, doi:10 .1175/2010JP04396.1.

--, M. H. Alford, R. Pinkel, R. C. Lien, and Y. J. Yang, 2011: The breaking and scattering of the internal tide on a continental slope. J. Phys. Oceanogr., 41, 926-945, doi:10.1175/2010JP04500.1.

--, M. Buijsman, S. Legg, and R. Pinkel, 2013: Parameterizing surface and internal tide scattering and breaking on supercritical topography: The one- and two-ridge cases. J. Phys. Oceanogr., 43, 1380-1397, doi:10.1175/JPO-D-12-061.1.

Kunze, E., 2017: Internal-wave-driven mixing: Geography and budgets. J. Phys. Oceanogr., 47, 1325-1345, doi:10.1175/JPO-D-16-0141.1.

--, R. W. Schmitt, and J. M. Toole, 1995: The energy balance in a warm-core ring's near-inertial critical layer. J. Phys. Oceanogr., 25, 942-957, doi: 10.1175/1520-0485(1995)025<0942:TEBI AW >2,O.CO;2.

--, E. Firing, J. Hummon, T. K. Chereskin, and A. Thurnherr, 2006: Global abyssal mixing inferred from lowered ADCP shear and CTD strain profiles. J. Phys. Oceanogr., 36, 1553-1576, doi:10.1175/JP02926.1.

--, C. MacKay, E. E. McPhee-Shaw, K. Morrice, J. B. Girton, and S. R. Terker, 2012: Turbulent mixing and exchange with interior waters on sloping boundaries. J. Phys. Oceanogr., 42, 910-927, doi:10.1175 /JPO-D-11-075.1.

Large, W., and G. Crawford, 1995: Observations and simulations of upper-ocean response to wind events during the ocean storms experiment. J. Phys. Oceanogr., 25, 2831-2852, doi:10.1175/1520-0485 (1995)025<2831:OASOUO>2,O.CO;2.

--, J. McWilliams, and S. Doney, 1994: Oceanic vertical mixing: A review and a model with a nonlocal boundary layer parameterization. Rev. Geophys., 32, 363-403, doi:10.1029/94RG01872.

Lefauve, A., C. Muller, and A. Melet, 2015: A threedimensional map of tidal dissipation over abyssal hills. J. Geophys. Res. Oceans, 120, 4760-4777, doi:10 .1002/2014JC010598.

Legg, S., 2014: Scattering of low-mode internal waves at finite isolated topography. J. Phys. Oceanogr., 44, 359-383, doi: 10.1175/JPO-D-12-0241.1.

--, and A. Adcroft, 2003: Internal wave breaking at concave and convex continental slopes. J. Phys. Oceanogr., 33, 2224-2246, doi:10.1175/1520-0485 (2003)033<2224:IWBACA>2,O.CO;2.

--, and J. M. Klymak, 2008: Internal hydraulic jumps and overturning generated by tidal flow over a steep ridge. J. Phys. Oceanogr., 38, 1949-1964, doi:10.1175 /2008JP03777.1.

--, R. Hallberg, and J. Girton, 2006: Comparison of entrainment in overflows simulated by z-coordinate, isopycnal and non-hydrostatic models. Ocean Modell, 11, 69-97, doi:10.1016/j.ocemod.2004.11.006.

Lighthill, J., 1978: Waves in Fluids. Cambridge University Press, 504 pp.

Llewellyn Smith, S. G., and W. R. Young, 2003: Tidal conversion at a very steep ridge. J. Fluid Mech., 495, 175-191, doi:10.1017/S0022112003006098.

Lumpkin, R., and K. Speer, 2007: Global ocean meridional overturning. J. Phys. Oceanogr., 37, 2550-2562, doi:10.1175/JP03130.1.

Lvov, Y., K. L. Polzin, and E. G. Tabak, 2004: Energy spectra of the ocean's internal wave field: Theory and observations. Phys. Rev. Lett., 92, 128501, doi:10.1103 /PhysRevLett.92.128501.

MacKinnon, J., 2013: Mountain waves in the deep ocean. Nature, 501, 321-322, doi:10.1038/501321a.

--, and M. C. Gregg, 2003: Shear and baroclinic energy flux on the summer New England shelf. J. Phys. Oceanogr., 33, 1462-1475, doi:10.1175/1520 -0485(2003)033<1462:SA BEFO>2.0.CO;2.

--, and K. Winters, 2003: Spectral evolution of bottom-forced internal waves. Near-Boundary Processes and Their Parameterization: Proc. 'Aha Huliko'a Hawaiian Winter Workshop, Honolulu, HI, University of Hawai'i at Manoa, 73-83.

--, and --, 2005: Subtropical catastrophe: Significant loss of low-mode tidal energy at 28.9[degrees]. Geophys. Res. Lett., 32, L15605, doi:10.1029/2005GL023376.

--, M. H. Alford, R. Pinkel, J. Klymak, and Z. Zhao, 2013a: The latitudinal dependence of shear and mixing in the Pacific transiting the critical latitude for PSI. J. Phys. Oceanogr.,43, 3-16, doi:10.1175/JPO-D-ll-0107.1.

--, --, O. Sun, R. Pinkel, Z. Zhao, and J. Klymak, 2013b: Parametric subharmonic instability of the internal tide at 29[degrees]N. J. Phys. Oceanogr., 43, 17-28, doi:10.1175/JPO-D-l 1-0108.1.

--, L. St. Laurent, and A. C. Naveira Garabato, 2013c: Diapycnal mixing processes in the ocean interior. Ocean Circulation and Climate: A 21st Century Perspective, 2nd ed. G. Siedler et al., Eds., International Geophysics Series, Vol. 103, Academic Press, 159-183.

Martini, K. I., M. H. Alford, E. Kunze, S. M. Kelly, and J. D. Nash, 2011: Observations of internal tides on the Oregon continental slope. J. Phys. Oceanogr., 41, 1772-1794, doi:10.1175/2011JP04581.1.

--, --, --, --, and --, 2013: Internal bores and breaking internal tides on the Oregon continental slope. J. Phys. Oceanogr., 43, 120-139, doi:10.1175 /JPO-D-12-030.1.

Mashayek, A., C. Caulfield, and W. Peltier, 2013: Time-dependent, non monotonic mixing in stratified turbulent shear flows: Implications for oceanographic estimates of buoyancy flux. J. Fluid Mech., 736, 570-593, doi:10.1017/jfm.2013.551.

Mathur, M., G. S. Carter, and T. Peacock, 2014: Topographic scattering of the low-mode internal tide in the deep ocean. J. Geophys. Res. Oceans, 119, 2165-2182, doi:10.1002/2013JC009152.

McComas, C. H., 1977: Equilibrium mechanisms within the oceanic internal wave field. J. Phys. Oceanogr., 7, 836-845, doi:10.1175/1520-0485(1977)007<0836: EMWTOI>2.0.CO;2.

Melet, A., R. Hallberg, S. Legg, and K. L. Polzin, 2013a: Sensitivity of the ocean state to the vertical distribution of internal-tide-driven mixing. J. Phys. Oceanogr., 43,602-615, doi:10.1175/JPO-D-12-055.1.

--, M. Nikurashin, C. J. Muller, S. Falahat, J. Nycander, P. G. Timko, B. K. Arbic, and J. A. Goff, 2013b: Internal tide generation by abyssal hills using analytical theory. J. Geophys. Res. Oceans, 118, 6303-6318, doi:10.1002/2013JC009212.

--, R. Hallberg, S. Legg, and M. Nikurashin, 2014: Sensitivity of the ocean state to lee wave-driven mixing. J. Phys. Oceanogr., 44, 900-921, doi:10.1175 /JPO-D-13-072.1.

--, --, A. Adcroft, M. Nikurashin, and S. Legg,

2015: Energy flux into internal lee waves: Sensitivity to future climate changes using linear theory and a climate model./. Climate, 28, 2365-2384, doi:10.1175 /JCLI-D-14-00432.1.

--, S. Legg, and R. Hallberg, 2016: Climatic impacts of parameterized local and remote tidal mixing. J. Climate, 29, 3473-3500, doi:10.1175/JCI.I-D-15-0153.1.

Mourn, J. N., K. Pujiana, R.-C. Lien, and W. D. Smyth, 2016: Ocean feedback to pulses of the Madden-Julian oscillation in the equatorial Indian Ocean. Nat. Commun., 7, 13 203, doi:10.1038/ncomms13203.

Muller, C. J., and O. Biihler, 2009: Saturation of the internal tides and induced mixing in the abyssal ocean. J. Phys. Oceanogr., 39, 2077-2096, doi:10.1175/2009JP04141.1.

Muller, M., J. Cherniawsky, M. Foreman, and J.-S. von Storch, 2012: Global map of M2 internal tide and its seasonal variability from high resolution ocean circulation and tide modelling. Geophys. Res. Lett., 39, L19607, doi:10.1029/2012GL053320.

--, B. K. Arbic, J. G. Richman, J. F. Shriver, E. L. Kunze, R. B. Scott, A. J. Wallcraft, and L. Zamudio, 2015: Toward an internal gravity wave spectrum in global ocean models. Geophys. Res. Lett., 42, 3474-3481, doi:10.1002/2015GL063365.

Muller, P., and A. Natarov, 2003: The Internal Wave Action Model (IWAM). Near-Boundary Processes and Their Parameterization: Proc.'Aha Huliko'a Winter Workshop, Honolulu, HI, University of Hawai'i at Manoa, 95-105.

--, G. Holloway, F. Henyey, and N. Pomphrey, 1986: Nonlinear interactions among internal gravity waves. Rev. Geophys., 24, 493-536, doi:10.1029 /RG024i003p00493.

Musgrave, R., R. Pinkel, J. MacKinnon, M. R. Mazloff, and W. Young, 2016: Stratified tidal flow over a tall ridge above and below the turning latitude. J. Fluid Mech., 793, 933-957, doi:10.1017/jfm.2016.150.

Nagai, T., and T. Hibiya, 2015: Internal tides and associated vertical mixing in the Indonesian Archipelago. J. Geophys. Res. Oceans, 120, 3373-3390, doi:10 .1002/2014JC010592.

Nash, J. D., E. Kunze, J. M. Toole, and R. W. Schmitt, 2004: Internal tide reflection and turbulent mixing on the continental slope. J. Phys. Oceanogr., 34, 1117-1134, doi:10.1175/1520-0485(2004)034<1117: ITRATM>2,O.CO;2.

--, M. H. Alford, E. Kunze, K. I. Martini, and S. Kelly, 2007: Hotspots of deep ocean mixing on the Oregon continental slope. Geophys. Res. Lett., 34, L01605, doi:10.1029/2006GL028170.

Naveira Garabato, A. C., K. L. Polzin, B. A. King, K. J. Heywood, and M. Visbeck, 2004: Widespread intense turbulent mixing in the Southern Ocean. Science, 303, 210-213, doi:10.1126/science. 1090929.

Nikurashin, M., and R. Ferrari, 2011: Global energy conversion rate from geostrophic flows into internal lee waves in the deep ocean. Geophys. Res. Lett., 38, L08610, doi:10.1029/2011GL046576.

--, and S. Legg, 2011: A mechanism for local dissipation of internal tides generated at rough topography. J. Phys. Oceanogr., 41, 378-395, doi:10.1175/2010JP04522.1.

--, and R. Ferrari, 2013: Overturning circulation driven by breaking internal waves in the deep ocean. Geophys. Res. Lett., 40, 3133-3137, doi:10.1002 /grl.50542.

--,--, N. Grisouard, and K. Polzin, 2014: The impact of finite-amplitude bottom topography on internal wave generation in the Southern Ocean. J. Phys. Oceanogr., 44, 2938-2950, doi:10.1175/JPO -D-13-0201.1.

Niwa, Y., and T. Hibiya, 2011: Estimation of baroclinic tide energy available for deep ocean mixing based on three-dimensional global numerical simulations. J. Oceanogr., 67, 493-502, doi:10.1007/s10872-011 -0052-1.

--, and--, 2014: Generation of baroclinic tide energy in a global three-dimensional numerical model with different spatial grid resolutions. Ocean Modell., 80,59-73, doi:10.1016/j.ocemod.2014.05.003.

Osborn, T. R., 1980: Estimates of the local rate of vertical diffusion from dissipation measurements. J. Phys. Oceanogr., 10, 83-89, doi:10.1175/1520 -0485(1980)010<0083:EOTLRO>2.0.CO;2.

Pacanowski, R. C., and G. Philander, 1981: Parameterization of vertical mixing in numerical models of the tropical ocean. J. Phys. Oceanogr., 11, 1443-1451, doi:10.1175/1520-0485(1981)011<1443: POVMIN>2,0.CO;2.

Palmer, W. R., G. J. Shutts, and R. Swinbank, 1986: Alleviation of systematic westerly bias in general circulation and numerical weather prediction models through an orographic gravity wave drag parameterization. Quart. J. Roy. Meteor. Soc., 112, 1001-1039, doi:10.1002/qj.49711247406.

Peltier, W. R., and C. P. Caulfield, 2003: Mixing efficiency in stratified shear flows. Annu. Rev. Fluid Mech., 35, 135-167, doi:10.1146/annurev.fluid .35.101101.161144.

Pierrehumbert, R., and J. Bacmeister, 1987: On the realizability of Long's Model solutions for nonlinear stratified flow over an obstacle. Stratified Flows, ASCE, 99-112.

Pinkel, R., and Coauthors, 2015: Breaking internal tides keep the ocean in balance. Eos, 96, doi:10.1029/2015E0039555.

Plueddemann, A. J., and J. T. Farrar, 2006: Observations and models of the energy flux from the wind to mixed-layer inertial currents. Deep-Sea Res., 53, 5-30, doi:10.1016/j.dsr2.2005.10.017.

Polzin, K. L., 2004a: A flux representation of internal wave spectral transports. J. Phys. Oceanogr., 34, 214-230, doi: 10.1175/1520-0485(2004)034<0214: AHDOIW>2.0.CO;2.

--, 2004b: Idealized solutions for the energy balance of the finescale internal wave field. J. Phys. Oceanogr., 34, 231-246, doi:10.1175/1520-0485(2004)034<0231: ISFTEB>2.0.CO;2.

--, 2009: An abyssal recipe. Ocean Modell., 30, 298-309, doi:10.1016/j.ocemod.2009.07.006.

--, 2010: Mesoscale eddy-internal wave coupling. Part II: Energetics and results from PolyMode. J. Phys. Oceanogr., 40, 789-801, doi:10.1175/2009JP04039.1.

--, J. M. Toole, and R. W. Schmitt, 1995: Finescale parameterizations of turbulent dissipation. J. Phys. Oceanogr., 25, 306-328, doi:10.1175/1520-0485(1995) 025<0306:FPOTD>2.0.CO;2.

--, N. S. Oakey, J. M. Toole, and R. W. Schmitt, 1996: Fine structure and microstructure characteristics across the northwest Atlantic subtropical front. J. Geophys. Res., 101, 14111-14121, doi:10.1029/96JC01020.

--, J. M. Toole, J. R. Ledwell, and R. W. Schmitt, 1997: Spatial variability of turbulent mixing in the abyssal ocean. Science, 276, 93-96, doi:10.1126 /science.276.5309.93.

--, A. C. Naveira Garabato, T. N. Huussen, B. M. Sloyan, and S. Waterman, 2014: Finescale parameterizations of turbulent dissipation. J. Geophys. Res. Oceans, 119, 1383-1419, doi:10.1002/2013JC008979.

Rahmstorf, S., 1993: A fast and complete convection scheme for ocean models. Ocean Modell., 101, 9-11.

Rainville, L., and R. Pinkel, 2006: Baroclinic energy flux at the Hawaiian Ridge: Observations from the R/P FLIP. J. Phys. Oceanogr., 36, 1104-1122, doi:10.1175 /JP02882.1.

Ray, R. D., and G. T. Mitchum, 1996: Surface manifestation of internal tides generated near Hawaii. Geophys. Res. Lett., 23, 2101-2104, doi:10.1029/96GL02050.

Rimac, A., J.-S. von Storch, C. Eden, and H. Haak, 2013: The influence of high-resolution wind stress field on the power input to near-inertial motions in the ocean. Geophys. Res. Lett., 40, 4882-4886, doi:10.1002/grl.50929.

--, --, and --, 2016: The total energy flux leaving the ocean's mixed layer. J. Phys. Oceanogr., 46, 1885-1900, doi:10.1175/JPO-D-15-OH5.1.

Rocha, C. B., T. K. Chereskin, S. T. Gille, and D. Menemenlis, 2016: Mesoscale to submesoscale wavenumber spectra in Drake Passage. J. Phys. Oceanogr., 46,601-620, doi:10.1175/JPO-D-15-0087.1.

Rudnick, D., and Coauthors, 2003: From tides to mixing along the Hawaiian Ridge. Science, 301, 355-357, doi:10.1126/science.1085837.

Salehipour, H., W. R. Peltier, C. B. Whalen, and J. A. Mackinnon, 2016: A new characterization of the turbulent diapycnal diffusivities of mass and momentum in the ocean. Geophys. Res. Lett., 43, 3370-3379, doi:10.1002/2016GL068184.

Schmitt, R. W., J. R. Ledwell, E. T. Montgomery, K. L. Polzin, and J. M. Toole, 2005: Enhanced diapycnal mixing by salt fingers in the thermocline of the tropical Atlantic. Science, 308, 685-688, doi:10.1126 /science.1108678.

Schmittner, A., and G. Egbert, 2014: An improved parameterization of tidal mixing for ocean models. Geosci. Model Dev., 7, 211-224, doi:10.5194/gmd-7 -211-2014.

Scott, R., J. Goff, A. Garabato, and A. Nurser, 2011: Global rate and spectral characteristics of internal gravity wave generation by geostrophic flow over topography. J. Geophys. Res., 116, C09029, doi:10.1029/2011JC007005.

Sheen, K. L., and Coauthors, 2013: Rates and mechanisms of turbulent dissipation and mixing in the Southern Ocean: Results from the Diapycnal and Isopycnal Mixing Experiment in the Southern Ocean (DIMES). J. Geophys. Res. Oceans, 118, 2774-2792, doi: 10.1002/jgrc.20217.

--, and Coauthors, 2014: Eddy-induced variability in southern ocean abyssal mixing on climatic timescales. Nat. Geosci., 7, 577-582, doi:10.1038 /ngeo2200.

Shriver, J., B. K. Arbic, J. Richman, R. Ray, E. Metzger, A. Wallcraft, and P. Timko, 2012: An evaluation of the barotropic and internal tides in a high-resolution global ocean circulation model. J. Geophys. Res., 117, C10024, doi:10.1029/2012JC008170.

--, J. Richman, and B. Arbic, 2014: How stationary are the internal tides in a high-resolution global ocean circulation model? J. Geophys. Res. Oceans, 119, 2769-2787, doi:10.1002/2013JC009423.

Silverthorne, K. E., and J. M. Toole, 2009: Seasonal kinetic energy variability of near-inertial motions. J. Phys. Oceanogr., 39, 1035-1049, doi:10.1175/2008JP03920.1.

Simmons, H. L., 2008: Spectral modification and geographic redistribution of the semi-diurnal internal tide. Ocean Modell., 21, 126-138, doi:10.1016 /j.ocemod.2008.01.002.

--, and M. H. Alford, 2012: Simulating the long range swell of internal waves generated by ocean storms. Oceanography, 25, 30-41, doi:10.5670/oceanog .2012.39.

--, R. W. Hallberg, and B. K. Arbic, 2004a: Internal wave generation in a global baroclinic tide model. Deep-Sea Res. II, 51, 3043-3068, doi:10.1016 /j.dsr2.2004.09.015.

--, S. R. Jayne, L. C. St. Laurent, and A. J. Weaver, 2004b: Tidally driven mixing in a numerical model of the ocean general circulation. Ocean Modell, 6, 245-263, doi:10.1016/S1463-5003(03)00011-8.

Slinn, D. N., and J. J. Riley, 1996: Turbulent mixing in the oceanic boundary layer caused by internal wave reflection from sloping terrain. Dyn. Atmos. Oceans, 24, 51-62, doi:10.1016/0377-0265(95)00425-4.

Staquet, C., and J. Sommeria, 2002: Internal gravity waves: From instability to turbulence. Annu. Rev. Fluid Mech., 34, 559-593, doi:10.1146/annurev .fluid.34.090601.130953.

St. Laurent, L., 2008: Turbulent dissipation on the margins of the South China Sea. Geophys. Res. Lett., 35, L23615, doi:10.1029/2008GL035520.

--, and C. Garrett, 2002: The role of internal tides in mixing the deep ocean. J. Phys. Oceanogr., 32, 28822899, doi:10.1175/1520-0485(2002)032<2882:TROITI >2.0.CO;2.

--, and J. D. Nash, 2004: An examination of the radiative and dissipative properties of deep ocean internal tides. Deep-Sea Res. II, 51, 3029-3042, doi:10.1016 /j.dsr2.2004.09.008.

--, and H. L. Simmons, 2006: Estimates of power consumed by mixing in the ocean interior. J. Climate, 19, 4877-4890, doi:10.1175/JCLI3887.1.

--, H. Simmons, and S. Jayne, 2002: Estimating tidally driven mixing in the deep ocean. Geophys. Res. Lett., 29, 2106, doi:10.1029/2002GL015633.

--, A. C. Naveira Garabato, J. R. Ledwell, A. M. Thurnherr, J. M. Toole, and A. J. Watson, 2012: Turbulence and diapycnal mixing in Drake Passage. J. Phys. Oceanogr., 42, 2143-2152, doi:10.1175/JPO -D-12-027.1.

Sun, O. M., and R. Pinkel, 2012: Energy transfer from high-shear, low-frequency internal waves to high-frequency waves near Kaena Ridge, Hawai'i. J. Phys. Oceanogr., 42, 1524-1547, doi:10.1175/JPO -D-ll-0117.1.

--, and --, 2013: Subharmonic energy transfer from the semidiurnal internal tide to near-diurnal motions over Kaena Ridge, Hawai'i. J. Phys. Oceanogr., 43, 766-789, doi:10.1175/JPO-D-12-0141.1.

Tanaka, T., I. Yasuda, Y. Tanaka, and G. S. Carter, 2013: Numerical study on tidal mixing along the shelf break in the Green Belt in the southeastern Bering Sea. J. Geophys. Res. Oceans, 118, 6525-6542, doi:10.1002/2013JC009113.

Thoppil, P., J. Richman, and P. Hogan, 2011: Energetics of a global ocean circulation model compared to observations. Geophys. Res. Lett., 38, L15607, doi:10.1029/2011GL048347.

Thurnherr, A. M., L. C. St. Laurent, K. G. Speer, J. M. Toole, and J. R. Ledwell, 2005: Mixing associated with sills in a canyon on the midocean ridge flank. J. Phys. Oceanogr., 35, 1370-1381, doi:10.1175/JP02773.1.

Trossman, D. S., B. K. Arbic, S. T. Garner, J. A. Goff, S. R. Jayne, E. J. Metzger, and A. J. Wallcraft, 2013: Impact of parameterized lee wave drag on the energy budget of an eddying global ocean model. Ocean Modell., 72, 119-142, doi:10.1016/j.ocemod.2013.08.006.

--, S. Waterman, K. L. Polzin, B. K. Arbic, S. T. Garner, A. C. Naveira-Garabato, and K. L. Sheen, 2015: Internal lee wave closures: Parameter sensitivity and comparison to observations. J. Geophys. Res. Oceans, 120, 7997-8019, doi:10.1002/2015JC010892.

--, B. K. Arbic, J. G. Richman, S. T. Garner, S. R. Jayne, and A. J. Wallcraft, 2016: Impact of topographic internal lee wave drag on an eddying global ocean model. Ocean Modell., 97, 109-128, doi:10.1016 /j.ocemod.2015.10.013.

Venayagamoorthy, S. K., and O. Fringer, 2006: Numerical simulations of the interaction of internal waves with a shelf break. Phys. Fluids, 18, 076603, doi:10.1063/1.2221863.

--, and J. R. Koseff, 2016: On the flux Richardson number in stably stratified turbulence. J. Fluid Mech., 798, R1, doi:10.1017/jfm.2016.340.

Wain, D. J., M. C. Gregg, M. H. Alford, R. C. Lien, G. S. Carter, and R. A. Hall, 2013: Propagation and dissipation of the internal tide in upper Monterey Canyon. J. Geophys. Res. Oceans, 118, 4855-4877, doi:10.1002/jgrc.20368.

Waterhouse, A. F., and Coauthors, 2014: Global patterns of diapycnal mixing from measurements of the turbulent dissipation rate. J. Phys. Oceanogr., 44, 1854-1872, doi:10.1175/JPO-D-13-0104.1.

--, J. A. MacKinnon, R. C. Musgrave, S. M. Kelly, A. I. Pickering, and J. Nash, 2017: Internal tide convergence and mixing in a submarine canyon. J. Phys. Oceanogr., 47, 303-322, doi:10.1175/JPO -D-16-0073.1.

Waterman, S., A. C. Naveira Garabato, and K. L. Polzin, 2013: Internal waves and turbulence in the Antarctic Circumpolar Current. J. Phys. Oceanogr., 43, 259-282, doi:10.1175/JPO-D-l 1-0194.1.

--, K. L. Polzin, A. C. Naveira Garabato, K. L. Sheen, and A. Forryan, 2014: Suppression of internal wave breaking in the Antarctic Circumpolar Current near topography. J. Phys. Oceanogr., 44, 1466-1492, doi: 10.1175/JPO-D-12-0154.1.

Whalen, C. B., 2015: Illuminating spatial and temporal patterns of ocean mixing as inferred from Argo profiling floats. Ph.D. thesis, University of California, San Diego, 97 pp.

--, L. D. Talley, and J. A. MacKinnon, 2012: Spatial and temporal variability of global ocean mixing inferred from Argo profiles. Geophys. Res. Lett., 39, L18612, doi:10.1029/2012GL053196.

--, J. A. MacKinnon, L. D. Talley, and A. F. Waterhouse, 2015: Estimating the mean diapycnal mixing using a finescale strain parameterization. J. Phys. Oceanogr., 45, 1174-1188, doi:10.1175/JPO-D-14-0167.1.

Wijesekera, H. W., L. Padman, T. Dillon, M. Levine, C. Paulson, and R. Pinkel, 1993: The application of internal-wave dissipation models to a region of strong mixing. J. Phys. Oceanogr., 23, 269-286, doi:10.1175 /1520-0485(1993)023<0269:TAOIWD>2.0.CO;2.

Winkel, D. P., M. C. Gregg, and T. B. Sanford, 2002: Patterns of shear and turbulence across the Florida Current. J. Phys. Oceanogr., 32,3269-3285, doi:10.1175 /1520-0485(2002)032<3269:POSATA>2.0.CO;2.

Wright, C. J., R. B. Scott, P. Ailliot, and D. Furnival, 2014: Lee wave generation rates in the deep ocean. Geophys. Res. Lett., 41, 2434-2440, doi:10.1002/2013GL059087.

Wu, L., Z. Jing, S. Riser, and M. Visbeck, 2011 : Seasonal and spatial variations of Southern Ocean diapycnal mixing from Argo profiling floats. Nat. Geosci., 4, 363-366, doi:10.1038/ngeo1156.

Wunsch, C., 1969: Progressive internal waves on slopes. J. Fluid Mech., 35, 131-145, doi:10.1017 /S0022112069001005.

Zhang, L., and H. L. Swinney, 2014: Virtual seafloor reduces internal wave generation by tidal flow. Phys. Rev. Lett., 112, 104502, doi:10.1103/PhysRevLett .112.104502.

Zhao, Z., M. H. Alford, J. A. MacKinnon, and R. Pinkel, 2010: Long-range propagation of the semidiurnal internal tide from the Hawaiian Ridge. J. Phys. Oceanogr., 40,713-736, doi:10.1175/2009JP04207.1.

--, --, J. B. Girton, L. Rainville, and H. L. Simmons, 2016: Global observations of open-ocean mode-1 [M.sub.2] internal tides. J. Phys. Oceanogr., 46, 1657-1684, doi:10.1175/JPO-D-15-0105.1.

Caption: Fig. 1. Schematic of internal wave mixing processes in the open ocean that are considered as part of this CPT. Tides interact with topographic features to generate high-mode internal waves (e.g., at midocean ridges) and low-mode internal waves (e.g., at tall steep ridges such as the Hawaiian Ridge). Deep currents flowing over topography can generate lee waves (e.g., in the Southern Ocean). Storms cause inertial oscillations in the mixed layer, which can generate both low- and high-mode internal waves (e.g., beneath storm tracks). In the open ocean, these internal waves can scatter off of rough topography and potentially interact with mesoscale fronts and eddies until they ultimately dissipate through wave-wave interactions. Internal waves that reach the shelf and slope can scatter or amplify as they propagate toward shallower water.

Caption: Fig. 2. Depth-averaged diffusivity [kappa] from (a) the upper ocean (from MLD to 1000-m depth) and (b) the full water column, updated from Waterhouse et al. (2014). The background diffusivity map in (a) comes from the strain-based inferences of diffusivity from Argo floats, updated from Whalen et al. (2015) with observations included from 2006 to 2015. (c) Compiled observations of mixing measurements with blue and green squares and diamonds denoting microstructure measurements. Green represents full-depth profiles, while blue denotes microstructure profiles. Purple circles represent inferred diffusivity from a finescale parameterization using lowered acoustic Doppler current profiler (LADCP)/conductivity-temperature-depth (CTD) profiles [dark purple, Kunze et al. (2006); medium purple, Huussen et al. (2012)] and High Density Sounding System (HDSS) shipboard shear (light orange). Dark orange circles are diffusivities from density overturns in moored profiles.

Caption: Fig. 3. (a) A snapshot of baroclinic velocity (m [s.sup.-1]) from a two-dimensional numerical simulation of internal tides forced by [M.sub.2] (semidiurnal) tidal velocities over rough topography for parameters corresponding to the Brazil Basin (Nikurashin and Legg 2011), (b) Observational time series of internal wave breaking over tall steep topography; here, we see (top) northward velocity and (bottom) turbulent dissipation rate oscillate twice a day as the tide flows over Kaena Ridge, Hawaii (Klymak et al. 2008). (c) Global energy flux from the [M.sub.2] tide into internal tides ([log.sub.10] W [m.sup.-2]) estimated using (top) the topography resolved in the Shuttle Radar Topography Mission (SRTM) global bathymetry and elevation data at 30 arc s resolution with data voids filled (SRTM30_PLUS) bathymetry database and (bottom) a statistical representation of unresolved abyssal hill topography estimates (Melet et al. 2013b). (d) The vertical structure of dissipation from Brazil Basin observations (thick solid curve) and the Polzin (2009) [Eq. (4)] parameterization of near-field internal tide dissipation (thin solid curve), as well as associated observed values of stratification ([N.sup.2]) and diapycnal diffusivity (K[rho]). (e) The impact of the Polzin parameterization in the GFDL CM2G coupled climate model: (top) the Indo-Pacific meridional overturning streamfunction (Sv; 1 Sv = [10.sup.6] [m.sup.3] [s.sup.-1]; averaged over the final 100 years of a 1000-yr simulation) using the Polzin (2009) parameterization and (bottom) the differences in Indo- Pacific meridional overturning streamfunction (Sv) between the simulations with the Polzin (2009) parameterization and the St. Laurent et al. (2002) parameterization as implemented by Simmons et al. (2004b) (from Melet et al. 2013a).

Caption: Fig. 4. Far-field internal tide: (a) SSH amplitude (mm) of global mode- 1 [M.sub.2] internal tides from multisatellite altimetry (Zhao et al. 2016). The light blue color indicates regions of high mesoscale activity, which make extraction of internal tides from altimetry difficult. Modeled semidiurnal tidal fluxes and comparison to observations: (b) HYCOM-modeled semidiurnal internal tide barotropic-to-baroclinic conversion rates (background color) and vertically integrated energy flux vectors (black arrows, plotted every 768th grid point for clarity) and (c) depth-integrated semidiurnal mode-1 energy fluxes in HYCOM (red arrows) and high-resolution mooring observations to the north of Hawaii (blue arrows) (Ansong et al. 2017). Impact on thermosteric sea level of using different spatial distributions of remote internal tide energy dissipation in GFDL ESM2G climate model: (d) thermosteric sea level (m) in a reference simulation using a constant background diapycnal diffusivity for remote internal tide dissipation. Anomalies (m) of thermosteric sea level from the reference case in (d) for simulations where (e) all internal tide energy is dissipated locally, over the generation site and (f) 20% of the internal tide energy is dissipated locally and 80% is dissipated uniformly over the ocean basins with a vertical profile proportional to buoyancy squared [N.sup.2] (Melet et al. 2016).

Caption: Fig. 5. Internal lee waves: (a) observations from DIMES showing (left) turbulent dissipation rates (logarithmic scales from [10.sup.-10] to [10.sup.-7] W [kg.sup.-1]) for the Phoenix Ridge (circles in right inset) and (right) average height above bottom profiles of turbulent kinetic energy dissipation [see details in St. Laurent et al. (2012)]. (b) Power conversion into lee waves [Nikurashin and Ferrari (2011) used in Melet et al. (2014)],(c) consequences of parameterized lee-wave mixing on the global ocean meridional overturning circulation [Sv; averaged over the final 100 years of 1000-yr simulations, from Melet et al. (2014)], and (d) global map of depth- integrated dissipation due to parameterized topographic wave drag inserted inline into global 1/25[degrees] HYCOM simulation, from Trossman et al. (2016).

Caption: Fig. 6. Near-inertial internal waves: (a) observational example from Alford et al. (2012) showing a (top) 2-yr record of wind work and (bottom) near-inertial kinetic energy in the northeastern Pacific, (b) One estimate of global power input (shading) and low-mode NIW energy fluxes (arrows; Simmons and Alford 2012). (c) Impact of near-inertial waves on annual-mean precipitation in ocean climate models: (top) the mean precipitation (mm [day.sup.-1]) from an experiment where the NI flux is set to 0.34 TW and (bottom) the same experiment, but with a doubling of the NI flux to 0.68 TW. The total tropical precipitation in the two experiments differs by less than 1%. An increase in near-inertial energy flux within observational uncertainties ameliorates the double ITCZs in the Atlantic and Pacific Oceans and creates the South Pacific convergence zone, three significant improvements for climate simulations of tropical precipitation.
COPYRIGHT 2017 American Meteorological Society
No portion of this article can be reproduced without the express written permission from the copyright holder.
Copyright 2017 Gale, Cengage Learning. All rights reserved.

Article Details
Printer friendly Cite/link Email Feedback
Author:MacKinnon, Jennifer A.; Zhao, Zhongxiang; Whalen, Caitlin B.; Waterhouse, Amy F.; Trossman, David S.
Publication:Bulletin of the American Meteorological Society
Article Type:Report
Date:Nov 1, 2017
Words:17587
Previous Article:CHRONOLOGICAL DEVELOPMENT OF TERRESTRIAL MEAN PRECIPITATION: Beginning with regional observations in preindustrial times, we focus on methods used...
Next Article:Comments on "Reanalyses and Observations: What's the Difference?".
Topics:

Terms of use | Privacy policy | Copyright © 2019 Farlex, Inc. | Feedback | For webmasters