Improved Internal Wave Spectral Continuum in a Regional Ocean Model

Recent work demonstrates that high‐resolution global models forced simultaneously by atmospheric fields and the astronomical tidal potential contain a partial internal (gravity) wave (IW) spectral continuum. Regional simulations of the MITgcm forced at the horizontal boundaries by a global run that carries a partial IW continuum spectrum are performed at the same grid spacing as the global run and at finer grid spacings in an attempt to fill out more of the IW spectral continuum. Decreasing only the horizontal grid spacing from 2 to 0.25 km greatly improves the frequency spectra and slightly improves the vertical wavenumber spectra of the horizontal velocity. Decreasing only the vertical grid spacing by a factor of 3 does not yield any significant improvements. Decreasing both horizontal and vertical grid spacings yields the greatest degree of improvement, filling the frequency spectrum out to 72 cpd. Our results suggest that improved IW spectra in regional models are possible if they are run at finer grid spacings and are forced at their lateral boundaries by remotely generated IWs. Additionally, consistency relations demonstrate that improvements in the spectra are indeed due to the existence of IWs at higher frequencies and vertical wavenumbers when remote IW forcing is included and model grid spacings decrease. By being able to simulate an IW spectral continuum to 0.25 km scales, these simulations demonstrate that one may be able to track the energy pathways of IWs from generation to dissipation and improve the understanding of processes such as IW‐driven mixing.Plain Language SummaryModels of internal waves (IWs) may help us to better understand the spatial geography of mixing in the ocean and are playing an increasingly important role in the planning of satellite missions. Following recent work showing that high‐resolution global models contain a partial IW spectrum, this paper describes further improvements in the spectrum seen in a high‐resolution regional model forced at the boundaries by a previously performed global IW simulation. Decreasing only the horizontal grid spacing greatly improves the frequency spectra and slightly improves the vertical wavenumber spectra of velocity. Increasing only the number of vertical levels does not yield any significant improvements. Decreasing both horizontal and vertical grid spacings yields the greatest improvement in both spectra. Our results suggest that regional models can exhibit improved IW spectra over global models if two conditions are met—they must have higher horizontal and vertical resolutions, and they must have remotely generated IWs at their boundaries. Application of the so‐called consistency relations demonstrates that the model is indeed carrying a field of high‐frequency IWs. Being able to simulate a fuller IW spectrum demonstrates that one may be able to use these models to improve the understanding of IW‐driven processes and energy pathways.Key PointsInternal gravity wave spectra in regional models are more realistic as model grid spacing decreasesThe vertical wavenumber spectra improve less dramatically than the frequency spectraInternal gravity wave consistency relations are applied to modeled spectraPeer Reviewedhttps://deepblue.lib.umich.edu/bitstream/2027.42/154917/1/jgrc23947_am.pdfhttps://deepblue.lib.umich.edu/bitstream/2027.42/154917/2/jgrc23947.pd

Energy transport by nonlinear internal waves

Winter stratification on Oregon’s continental shelf often produces a near-bottom layer of dense fluid that acts as an internal waveguide upon which nonlinear internal waves propagate. Shipboard profiling and bottom lander observations capture disturbances that exhibit properties of internal solitary waves, bores, and gravity currents. Wavelike pulses are highly turbulent (instantaneous bed stresses are 1 N m−2), resuspending bottom sediments into the water column and raising them 30+ m above the seafloor. The wave cross-shelf transport of fluid often counters the time-averaged Ekman transport in the bottom boundary layer. In the nonlinear internal waves that were observed, the kinetic energy is roughly equal to the available potential energy and is O(0.1) megajoules per meter of coastline. The energy transported by these waves includes a nonlinear advection term 〈uE〉 that is negligible in linear internal waves. Unlike linear internal waves, the pressure–velocity energy flux 〈up〉 includes important contributions from nonhydrostatic effects and surface displacement. It is found that, statistically, 〈uE〉 ≃ 2〈up〉. Vertical profiles through these waves of elevation indicate that up(z) is more important in transporting energy near the seafloor while uE(z) dominates farther from the bottom. With the wave speed c estimated from weakly nonlinear wave theory, it is verified experimentally that the total energy transported by the waves is 〈up〉 + 〈uE〉 ≃ c〈E〉. The high but intermittent energy flux by the waves is, in an averaged sense, O(100) watts per meter of coastline. This is similar to independent estimates of the shoreward energy flux in the semidiurnal internal tide at the shelf break

Convectively Driven Mixing in the Bottom Boundary Layer

Closely spaced vertical profiles through the bottom boundary layer over a sloping continental shelf during relaxation from coastal upwelling reveal structure that is consistent with convectively driven mixing. Parcels of fluid were observed adjacent to the bottom that were warm (by several millikelvin) relative to fluid immediately above. On average, the vertical gradient of potential temperature in the superadiabatic (statically unstable) bottom layer was found to be −1.7 × 10⁻⁴ K m⁻¹, or 6.0 × 10⁻⁵ kg m⁻⁴ in potential density. Turbulent dissipation rates (ε) increased toward the bottom but were relatively constant over the dimensionless depth range 0.4–1.0z/D (where D is the mixed layer height). The Rayleigh number Ra associated with buoyancy anomalies in the bottom mixed layer is estimated to be approximately 10¹¹, much larger than the value of approximately 10³ required to initiate convection in simple laboratory or numerical experiments. An evaluation of the data in which the bottom boundary layer was unstably stratified indicates that the greater the buoyancy anomaly is, the greater the turbulent dissipation rate in the neutral layer away from the bottom will be. The vertical structures of averaged profiles of potential density, potential temperature, and turbulent dissipation rate versus nondimensional depth are similar to their distinctive structure in the upper ocean during convection. Nearby moored observations indicate that periods of static instability near the bottom follow events of northward flow and local fluid warming by lateral advection. The rate of local fluid warming is consistent with several estimates of offshore buoyancy transport near the bottom. It is suggested that the concentration of offshore Ekman transport near the bottom of the Ekman layer when the flow atop the layer is northward can provide the differential transport of buoyant bottom fluid when the density in the bottom boundary layer decreases up the slope

A modified law-of-the-wall applied to oceanic bottom boundary layers

Near the bottom, the velocity profile in the bottom boundary layer over the continental shelf exhibits a characteristic law-of-the-wall that is consistent with local estimates of friction velocity from near-bottom turbulence measurements. Farther from the bottom, the velocity profile exhibits a deviation from the law-of-the-wall. Here the velocity gradient continues to decrease with height but at a rate greater than that predicted by the law-of-the-wall with the local friction velocity. We argue that the shape of the velocity profile is made consistent with the local friction velocity by the introduction of a new length scale that, near the boundary, asymptotes to a value that varies linearly from the bottom. Farther from the boundary, this length scale is consistent with the suppression of velocity fluctuations either by stratification in the upper part of the boundary layer or by proximity to the free surface. The resultant modified law-of-the-wall provides a good representation of velocity profiles observed over the continental shelf when a local estimate of the friction velocity from coincident turbulence observations is used. The modified law-of-the- wall is then tested on two very different sets of observations, from a shallow tidal channel and from the bottom of the Mediterranean outflow plume. In both cases it is argued that the observed velocity profile is consistent with the modified law-of-the-wall. Implicit in the modified law-of-the-wall is a new scaling for turbulent kinetic energy dissipation rate. This new scaling diverges from the law-of-the-wall prediction above 0.2D (where D is the thickness of the bottom boundary layer) and agrees with observed profiles to 0.6D

Diapycnal mixing in the Antarctic Circumpolar Current

Author Posting. © American Meteorological Society, 2011. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 41 (2011): 241-246, doi:10.1175/2010JPO4557.1.The vertical dispersion of a tracer released on a density surface near 1500-m depth in the Antarctic Circumpolar Current west of Drake Passage indicates that the diapycnal diffusivity, averaged over 1 yr and over tens of thousands of square kilometers, is (1.3 ± 0.2) × 10−5 m2 s−1. Diapycnal diffusivity estimated from turbulent kinetic energy dissipation measurements about the area occupied by the tracer in austral summer 2010 was somewhat less, but still within a factor of 2, at (0.75 ± 0.07) × 10−5 m2 s−1. Turbulent diapycnal mixing of this intensity is characteristic of the midlatitude ocean interior, where the energy for mixing is believed to derive from internal wave breaking. Indeed, despite the frequent and intense atmospheric forcing experienced by the Southern Ocean, the amplitude of finescale velocity shear sampled about the tracer was similar to background amplitudes in the midlatitude ocean, with levels elevated to only 20%–50% above the Garrett–Munk reference spectrum. These results add to a long line of evidence that diapycnal mixing in the interior middepth ocean is weak and is likely too small to dictate the middepth meridional overturning circulation of the ocean.This material is based upon work supported by the National Science Foundation Grants OCE-0622825,OCE-0622670, OCE-0622630, and OCE-0623177

Ocean convergence and the dispersion of flotsam

Floating oil, plastics, and marine organisms are continually redistributed by ocean surface currents. Prediction of their resulting distribution on the surface is a fundamental, long-standing, and practically important problem. The dominant paradigm is dispersion within the dynamical context of a nondivergent flow: objects initially close together will on average spread apart but the area of surface patches of material does not change. Although this paradigm is likely valid at mesoscales, larger than 100 km in horizontal scale, recent theoretical studies of submesoscales (less than ∼10 km) predict strong surface convergences and downwelling associated with horizontal density fronts and cyclonic vortices. Here we show that such structures can dramatically concentrate floating material. More than half of an array of ∼200 surface drifters covering ∼20 × 20 km2 converged into a 60 × 60 m region within a week, a factor of more than 105 decrease in area, before slowly dispersing. As predicted, the convergence occurred at density fronts and with cyclonic vorticity. A zipperlike structure may play an important role. Cyclonic vorticity and vertical velocity reached 0.001 s−1 and 0.01 ms−1, respectively, which is much larger than usually inferred. This suggests a paradigm in which nearby objects form submesoscale clusters, and these clusters then spread apart. Together, these effects set both the overall extent and the finescale texture of a patch of floating material. Material concentrated at submesoscale convergences can create unique communities of organisms, amplify impacts of toxic material, and create opportunities to more efficiently recover such material

Ocean convergence and the dispersion of flotsam

Floating oil, plastics, and marine organisms are continually redistributed by ocean surface currents. Prediction of their resulting distribution on the surface is a fundamental, long-standing, and practically important problem. The dominant paradigm is dispersion within the dynamical context of a nondivergent flow: objects initially close together will on average spread apart but the area of surface patches of material does not change. Although this paradigm is likely valid at mesoscales, larger than 100 km in horizontal scale, recent theoretical studies of submesoscales (less than ∼10 km) predict strong surface convergences and downwelling associated with horizontal density fronts and cyclonic vortices. Here we show that such structures can dramatically concentrate floating material. More than half of an array of ∼200 surface drifters covering ∼20 × 20 km2 converged into a 60 × 60 m region within a week, a factor of more than 105 decrease in area, before slowly dispersing. As predicted, the convergence occurred at density fronts and with cyclonic vorticity. A zipperlike structure may play an important role. Cyclonic vorticity and vertical velocity reached 0.001 s−1 and 0.01 ms−1, respectively, which is much larger than usually inferred. This suggests a paradigm in which nearby objects form submesoscale clusters, and these clusters then spread apart. Together, these effects set both the overall extent and the finescale texture of a patch of floating material. Material concentrated at submesoscale convergences can create unique communities of organisms, amplify impacts of toxic material, and create opportunities to more efficiently recover such material

Parametric subharmonic instability of the internal tide at 29°N

Author Posting. © American Meteorological Society, 2013. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 43 (2013): 17–28, doi:10.1175/JPO-D-11-0108.1.Observational evidence is presented for transfer of energy from the internal tide to near-inertial motions near 29°N in the Pacific Ocean. The transfer is accomplished via parametric subharmonic instability (PSI), which involves interaction between a primary wave (the internal tide in this case) and two smaller-scale waves of nearly half the frequency. The internal tide at this location is a complex superposition of a low-mode waves propagating north from Hawaii and higher-mode waves generated at local seamounts, making application of PSI theory challenging. Nevertheless, a statistically significant phase locking is documented between the internal tide and upward- and downward-propagating near-inertial waves. The phase between those three waves is consistent with that expected from PSI theory. Calculated energy transfer rates from the tide to near-inertial motions are modest, consistent with local dissipation rate estimates. The conclusion is that while PSI does befall the tide near a critical latitude of 29°N, it does not do so catastrophically.This work was sponsored by NSF OCE 04-25283.2013-07-0

Climate Process Team on internal wave–driven ocean mixing

Author Posting. © American Meteorological Society, 2017. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Bulletin of the American Meteorological Society 98 (2017): 2429-2454, doi:10.1175/BAMS-D-16-0030.1.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.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.2018-06-0

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