scholarly journals Abyssal Mixing through Critical Reflection of Equatorially Trapped Waves off Smooth Topography

2019 ◽  
Vol 49 (2) ◽  
pp. 519-542 ◽  
Author(s):  
Bertrand L. Delorme ◽  
Leif N. Thomas

AbstractThe inferred diapycnal upwelling in the abyssal meridional overturning circulation (AMOC) is intensified near the equator, but little is known as to why this is so. In this study, it is shown that the reflection of equatorially trapped waves (ETWs) off the bottom leads to seafloor-intensified mixing and substantial diapycnal upwelling near the equator when the full Coriolis force and the so-called nontraditional effects are taken into account. Using idealized simulations run with the MITgcm of downward-propagating ETWs of various types (i.e., inertia–gravity, Yanai, Kelvin, and Rossby waves) accounting for nontraditional effects, it is demonstrated that the reflection of ETWs off a flat seafloor generates beams of short inertia–gravity waves with strong vertical shear and low Richardson numbers that result in bottom-intensified, persistent, zonally invariant mixing at the inertial latitude of the ETW through the mechanism of critical reflection. The beams are more intense with weaker stratification and, for a given wave type, are stronger for waves with shorter periods and longer vertical wavelengths. The intensity of the beams also differs between wave types because their distinct meridional structures modulate the amount of energy fluxed to the bottom at the inertial latitude. As a result, equatorial inertia–gravity, Rossby, and eastward-propagating Yanai waves yield stronger mixing than Kelvin and westward-propagating Yanai waves in the simulations. It is estimated that this process can result in order 10 Sv (1 Sv ≡ 106 m3 s−1) of diapycnal upwelling per wavelength of ETW in the abyss and thus could play an important role in closing the AMOC.

2021 ◽  
Author(s):  
Adam Blaker ◽  
Michael Bell ◽  
Joel Hirschi ◽  
Amy Bokota

<p>Numerical model studies have shown the meridional overturning circulation (MOC) to exhibit variability on near-inertial timescales, and also indicate a region of enhanced variability on the equator. We present an analysis of a set of integrations of a global configuration of a numerical ocean model, which show very large amplitude oscillations in the MOCs in the Atlantic, Indian and Pacific oceans confined to the equatorial region. The amplitude of these oscillations is proportional to the width of the ocean basin, typically about 100 (200) Sv in the Atlantic (Pacific). We show that these oscillations are driven by surface winds within 10°N/S of the equator, and their periods (typically 4-10 days) correspond to a small number of low mode equatorially trapped planetary waves. Furthermore, the oscillations can be well reproduced by idealised wind-driven simulations linearised about a state of rest. Zonally integrated linearised equations of motion are solved using vertical normal modes and equatorial meridional modes representing Yanai and inertia-gravity waves. Idealised simulations capture between 85% and 95% of the variance of matching time-series segments diagnosed from the NEMO integrations. Similar results are obtained for the corresponding modes in the Atlantic and Indian Oceans. Our results raise questions about the roles of inertia-gravity waves near the equator in the vertical transfer of heat and momentum and whether these transfers will be explicitly resolved by ocean models or need to be parametrised.</p>


Author(s):  
Bertrand L. Delorme ◽  
Leif N. Thomas ◽  
Patrick Marchesiello ◽  
Jonathan Gula ◽  
Guillaume Roullet ◽  
...  

AbstractRecent theoretical work has shown that, when the so-called non-traditional effects are taken into account, the reflection of Equatorially Trapped Waves (ETWs) off the seafloor generates strong vertical shear that results in bottom-intensified mixing at the inertial latitude of the ETW via a mechanism of critical reflection. It has been estimated that this process could play an important role in driving diapycnal upwelling in the Abyssal Meridional Overturning Circulation (AMOC). However, these results were derived under an idealized configuration with a monochromatic ETW propagating through a flat ocean at rest. To test the theory in a flow that is more representative of the ocean, we contrast a set of realistic numerical simulations of the Eastern Equatorial Pacific run using either the hydrostatic or quasi-hydrostatic approximation, the latter of which accounts for non-traditional effects. The simulations are nested into a Pacific-wide hydrostatic parent solution forced with climatological data and realistic bathymetry, resulting in an ETW field and a deep circulation consistent with observations. Using these simulations, we observe enhanced abyssal mixing in the quasi-hydrostatic run, even over smooth topography, that is absent in the hydrostatic run. The mixing is associated with inertial shear that has spatio-temporal properties consistent with the critical reflection mechanism. The enhanced mixing results in a weakening of the abyssal stratification and drives diapycnal upwelling in our simulation, in agreement with the predictions from the idealized simulations. The diapycnal upwelling is on the order of O(10) Sv and thus could play an important role in closing the AMOC.


2013 ◽  
Vol 43 (2) ◽  
pp. 382-401 ◽  
Author(s):  
Julien Jouanno ◽  
Frédéric Marin ◽  
Yves du Penhoat ◽  
Jean-Marc Molines

Abstract A regional numerical model of the tropical Atlantic Ocean and observations are analyzed to investigate the intraseasonal fluctuations of the sea surface temperature at the equator in the Gulf of Guinea. Results indicate that the seasonal cooling in this region is significantly shaped by short-duration cooling events caused by wind-forced equatorial waves: mixed Rossby–gravity waves within the 12–20-day period band, inertia–gravity waves with periods below 11 days, and equatorially trapped Kelvin waves with periods between 25 and 40 days. In these different ranges of frequencies, it is shown that the wave-induced horizontal oscillations of the northern front of the mean cold tongue dominate the variations of mixed layer temperature near the equator. But the model mixed layer heat budget also shows that the equatorial waves make a significant contribution to the mixed layer heat budget through modulation of the turbulent cooling, especially above the core of the Equatorial Undercurrent (EUC). The turbulent cooling variability is found to be mainly controlled by the intraseasonal modulation of the vertical shear in the upper ocean. This mechanism is maximum during periods of seasonal cooling, especially in boreal summer, when the surface South Equatorial Current is strongest and between 2°S and the equator, where the presence of the EUC provides a background vertical shear in the upper ocean. It applies for the three types of intraseasonal waves. Inertia–gravity waves also modulate the turbulent heat flux at the equator through vertical displacement of the core of the EUC in response to equatorial divergence and convergence.


2014 ◽  
Vol 755 ◽  
pp. 397-428 ◽  
Author(s):  
Peng Wang ◽  
James C. McWilliams ◽  
Claire Ménesguen

AbstractThe linear instability of several rotating, stably stratified, interior vertical shear flows $\def \xmlpi #1{}\def \mathsfbi #1{\boldsymbol {\mathsf {#1}}}\let \le =\leqslant \let \leq =\leqslant \let \ge =\geqslant \let \geq =\geqslant \def \Pr {\mathit {Pr}}\def \Fr {\mathit {Fr}}\def \Rey {\mathit {Re}}\overline{U}(z)$ is calculated in Boussinesq equations. Two types of baroclinic, ageostrophic instability, AI1 and AI2, are found in odd-symmetric $\overline{U}(z)$ for intermediate Rossby number ($\mathit{Ro}$). AI1 has zero frequency; it appears in a continuous transformation of the unstable mode properties between classic baroclinic instability (BCI) and centrifugal instability (CI). It begins to occur at intermediate $\mathit{Ro}$ values and horizontal wavenumbers ($k,l$) that are far from $l= 0$ or $k = 0$, where the growth rate of BCI or CI is the strongest. AI1 grows by drawing kinetic energy from the mean flow, and the perturbation converts kinetic energy to potential energy. The instability AI2 has inertia critical layers (ICL); hence it is associated with inertia-gravity waves. For an unstable AI2 mode, the coupling is either between an interior balanced shear wave and an inertia-gravity wave (BG), or between two inertia-gravity waves (GG). The main energy source for an unstable BG mode is the mean kinetic energy, while the main energy source for an unstable GG mode is the mean available potential energy. AI1 and BG type AI2 occur in the neighbourhood of $A-S= 0$ (a sign change in the difference between absolute vertical vorticity and horizontal strain rate in isentropic coordinates; see McWilliams et al., Phys. Fluids, vol. 10, 1998, pp. 3178–3184), while GG type AI2 arises beyond this condition. Both AI1 and AI2 are unbalanced instabilities; they serve as an initiation of a possible local route for the loss of balance in 3D interior flows, leading to an efficient energy transfer to small scales.


2006 ◽  
Vol 63 (12) ◽  
pp. 3253-3276 ◽  
Author(s):  
Christoph Zülicke ◽  
Dieter Peters

Poleward-breaking Rossby waves often induce an upper-level jet streak over northern Europe. Dominant inertia–gravity wave packets are observed downstream of this jet. The physical processes of their generation and propagation, in such a configuration, are investigated with a mesoscale model. The study is focused on an observational campaign from 17 to 19 December 1999 over northern Germany. Different simulations with the fifth-generation Pennsylvania State University–National Center for Atmospheric Research (PSU–NCAR) Mesoscale Model (MM5) have been performed. For a high-resolution process study, three domains were set up that encompass the evolution of Rossby waves and that of inertia–gravity waves. To minimize the impact of model damping, the horizontal and vertical resolution has been adjusted appropriately. With a novel statistical approach, the properties of inertia–gravity wave packets have been estimated. This method uses the horizontal divergence field and takes into account the spatial extension of a wave packet. It avoids the explicit treatment of the background field and works for arbitrary wavelength. Two classes of inertia–gravity waves were found: subsynoptic waves with a horizontal wavelength of about 500 km and mesoscale waves with a horizontal wavelength of about 200 km. The subsynoptic structures were also detected in radiosonde observations during this campaign. The similarity between simulated and observed wavelengths and amplitudes suggests that the simulations can be considered as near realistic. Spontaneous radiation from unbalanced flow is an important process of inertia–gravity wave generation. Synoptic-scale imbalances in the exit region of the upper-tropospheric jet streak were identified with the smoothed cross-stream Lagrangian Rossby number. In a number of simulations with different physics, it was found that the inertia–gravity wave activity was related to the tropospheric jet, orography, and moist convection. The upward propagation of inertia–gravity waves was favored during this event of a poleward-breaking Rossby wave. The presence of the polar vortex induced background winds exceeding the critical line. Consequently, the activity of inertia–gravity waves in the lower stratosphere increased by an order of magnitude during the case study. The successful simulation of the complex processes of generation and propagation showed the important role of poleward Rossby wave breaking for the appearance of inertia–gravity waves in the midlatitudes.


2010 ◽  
Vol 664 ◽  
pp. 478-509 ◽  
Author(s):  
VICTOR I. SHRIRA ◽  
WILLIAM A. TOWNSEND

The paper is concerned with analytical study of inertia-gravity waves in rotating density-stratified ideal fluid confined in a spherical shell. It primarily aims at clarifying the possible role of these motions in deep ocean mixing. Recently, it was found that on the ‘non-traditional’ β-plane inertia-gravity internal waves can propagate polewards beyond their inertial latitude, where the wave frequency equals the local Coriolis parameter, by turning into subinertial modes trapped in the narrowing waveguides around the local minima of buoyancy frequency N. The behaviour of characteristics was established: wave horizontal and vertical scales decrease as the wave advances polewards and tend to zero at a latitude corresponding to an attractor of characteristics. However, the basic questions about wave evolution, its quantitative description and the possibility of its reflection from the critical latitude remain open. The present work addresses these issues by studying the linear inviscid evolution of finite bandwidth wavepackets on the ‘non-traditional’ β-plane past the inertial latitude for generic oceanic stratification. Beyond the inertial latitude, the wave field is confined in narrowing waveguides of three distinct generic types around different local minima of the buoyancy frequency. In the oceanic context, the widest is adjacent to the flat bottom, the thinnest is the upper mixed layer, and the middle one is located between the seasonal and main thermocline. We find explicit asymptotic solutions describing the wave field in the WKB approximation. As a byproduct, the conservation of wave action principle is explicitly formulated for all types of internal waves on the ‘non-traditional’ β-plane. The wave velocities and vertical shear tend to infinity and become singular at the attractor latitude or its vicinity for both monochromatic and finite bandwidth packets. We call this phenomenon singular focusing. These WKB solutions are shown to remain valid up to singularity for the bottom and mid-ocean waveguides. The main conclusion is that even in the inviscid setting the wave evolution towards smaller and smaller horizontal and vertical scales is irreversible: there is no reflection. For situations typical of deep ocean, a simultaneous increase in wave amplitude and decrease of vertical scale causes a sharp increase of vertical shear, which may lead to wave breaking and increased mixing.


2013 ◽  
Vol 70 (6) ◽  
pp. 1603-1615 ◽  
Author(s):  
Takenari Kinoshita ◽  
Kaoru Sato

Abstract A companion paper formulates the three-dimensional wave activity flux (3D-flux-M) whose divergence corresponds to the wave forcing on the primitive equations. However, unlike the two-dimensional wave activity flux, 3D-flux-M does not accurately describe the magnitude and direction of wave propagation. In this study, the authors formulate a modification of 3D-flux-M (3D-flux-W) to describe this propagation using small-amplitude theory for a slowly varying time-mean flow. A unified dispersion relation for inertia–gravity waves and Rossby waves is also derived and used to relate 3D-flux-W to the group velocity. It is shown that 3D-flux-W and the modified wave activity density agree with those for inertia–gravity waves under the constant Coriolis parameter assumption and those for Rossby waves under the small Rossby number assumption. To compare 3D-flux-M with 3D-flux-W, an analysis of the European Centre for Medium-Range Weather Forecasts (ECMWF) Interim Re-Analysis (ERA-Interim) data is performed focusing on wave disturbances in the storm tracks during April. While the divergence of 3D-flux-M is in good agreement with the meridional component of the 3D residual mean flow associated with disturbances, the 3D-flux-W divergence shows slight differences in the upstream and downstream regions of the storm tracks. Further, the 3D-flux-W magnitude and direction are in good agreement with those derived by R. A. Plumb, who describes Rossby wave propagation. However, 3D-flux-M is different from Plumb’s flux in the vicinity of the storm tracks. These results suggest that different fluxes (both 3D-flux-W and 3D-flux-M) are needed to describe wave propagation and wave–mean flow interaction in the 3D formulation.


2007 ◽  
Vol 64 (12) ◽  
pp. 4417-4431 ◽  
Author(s):  
Chris Snyder ◽  
David J. Muraki ◽  
Riwal Plougonven ◽  
Fuqing Zhang

Abstract Vortex dipoles provide a simple representation of localized atmospheric jets. Numerical simulations of a synoptic-scale dipole in surface potential temperature are considered in a rotating, stratified fluid with approximately uniform potential vorticity. Following an initial period of adjustment, the dipole propagates along a slightly curved trajectory at a nearly steady rate and with a nearly fixed structure for more than 50 days. Downstream from the jet maximum, the flow also contains smaller-scale, upward-propagating inertia–gravity waves that are embedded within and stationary relative to the dipole. The waves form elongated bows along the leading edge of the dipole. Consistent with propagation in horizontal deformation and vertical shear, the waves’ horizontal scale shrinks and the vertical slope varies as they approach the leading stagnation point in the dipole’s flow. Because the waves persist for tens of days despite explicit dissipation in the numerical model that would otherwise damp the waves on a time scale of a few hours, they must be inherent features of the dipole itself, rather than remnants of imbalances in the initial conditions. The wave amplitude varies with the strength of the dipole, with waves becoming obvious once the maximum vertical vorticity in the dipole is roughly half the Coriolis parameter. Possible mechanisms for the wave generation are spontaneous wave emission and the instability of the underlying balanced dipole.


1999 ◽  
Vol 17 (1) ◽  
pp. 115-121 ◽  
Author(s):  
L. Thomas ◽  
R. M. Worthington ◽  
A. J. McDonald

Abstract. Radar measurements at Aberystwyth (52.4° N, 4.1° W) of winds at tropospheric and lower stratospheric heights are shown for 12-13 March 1994 in a region of highly curved flow, downstream of the jet maximum. The perturbations of horizontal velocity have comparable amplitudes in the troposphere and lower stratosphere with downward and upward phase propagation, respectively, in these two height regions. The sense of rotation with increasing height in hodographs of horizontal perturbation velocity derived for hourly intervals show downwards propagation of energy in the troposphere and upward propagation in the lower stratosphere with vertical wavelengths of 1.7 to 2.3 km. The results indicate inertia-gravity waves propagating in a direction similar to that of the jet stream but at smaller velocities. Some of the features observed contrast with those of previous observations of inertia-gravity waves propagating transverse to the jet stream. The interpretation of the hodographs to derive wave parameters has taken account of the vertical shear of the background wind transverse to the direction of wave propagation.Key words. Meteorology and atmospheric dynamics (mesoscale meteorology; middle atmosphere dynamics; waves and tides)


2013 ◽  
Vol 70 (6) ◽  
pp. 1577-1602 ◽  
Author(s):  
Takenari Kinoshita ◽  
Kaoru Sato

Abstract The three-dimensional (3D) residual mean flow is expressed as the sum of the Eulerian-mean flow and the Stokes drift. The present study derives formulas that are approximately equal to the 3D Stokes drift for the primitive equation (PRSD) and for the quasigeostrophic equation (QGSD) using small-amplitude theory for a slowly varying time-mean flow. The PRSD has a broad utility that is applicable to both Rossby waves and inertia–gravity waves. The 3D wave activity flux whose divergence corresponds to the wave forcing is also derived using PRSD. The PRSD agrees with QGSD under the small-Rossby-number assumption, and it agrees with the 3D Stokes drift derived by S. Miyahara and by T. Kinoshita et al. for inertia–gravity waves under the constant-Coriolis-parameter assumption. Moreover, a phase-independent 3D Stokes drift is derived under the QG approximation. The 3D residual mean flow in the upper troposphere in April is investigated by applying the new formulas to the European Centre for Medium-Range Weather Forecasts (ECMWF) Interim Re-Analysis (ERA-Interim) data. It is observed that the PRSD is strongly poleward (weakly equatorward) upstream (downstream) of the storm track. A case study was also made for dominant gravity waves around the southern Andes in the simulation by a gravity wave–resolving general circulation model. The 3D residual mean flow associated with the gravity waves is poleward (equatorward) in the western (eastern) region of the southern Andes. This flow is due to the horizontal structure of the variance in the zonal component of the mountain waves, which do not change much while they propagate upward.


Sign in / Sign up

Export Citation Format

Share Document