Browse
Abstract
The Greenland Sea produces a significant portion of the dense water from the Nordic seas that supplies the lower limb of the Atlantic meridional overturning circulation. Here, we use a continuous 10-yr hydrographic record from moored profilers to examine dense-water formation in the central Greenland Sea between 1999 and 2009. Of primary importance for dense-water formation is air–sea heat exchange, and 60%–80% of the heat lost to the atmosphere during winter occurs during intense, short-lived events called cold-air outbreaks (CAOs). The long duration and high temporal resolution of the moored record has for the first time facilitated a statistical quantification of the direct impact of CAOs on the wintertime mixed layer in the Greenland Sea. The mixed layer development can be divided into two phases: a cooling phase and a deepening phase. During the cooling phase (typically between November and January), CAOs cooled the mixed layer by up to 0.08 K day−1, depending on the intensity of the events, while the mixed layer depth remained nearly constant. Later in winter (February–April), heat fluxes during CAOs primarily led to mixed layer deepening of up to 38 m day−1. Considerable variability was observed in the mixed layer response, indicating that lateral fluxes of heat and salt were also important. The magnitude and vertical distributions of these fluxes were quantified, and idealized mixed layer simulations suggest that their combined effect is a reduction in the mixed layer depth at the end of winter of up to several hundred meters.
Abstract
The Greenland Sea produces a significant portion of the dense water from the Nordic seas that supplies the lower limb of the Atlantic meridional overturning circulation. Here, we use a continuous 10-yr hydrographic record from moored profilers to examine dense-water formation in the central Greenland Sea between 1999 and 2009. Of primary importance for dense-water formation is air–sea heat exchange, and 60%–80% of the heat lost to the atmosphere during winter occurs during intense, short-lived events called cold-air outbreaks (CAOs). The long duration and high temporal resolution of the moored record has for the first time facilitated a statistical quantification of the direct impact of CAOs on the wintertime mixed layer in the Greenland Sea. The mixed layer development can be divided into two phases: a cooling phase and a deepening phase. During the cooling phase (typically between November and January), CAOs cooled the mixed layer by up to 0.08 K day−1, depending on the intensity of the events, while the mixed layer depth remained nearly constant. Later in winter (February–April), heat fluxes during CAOs primarily led to mixed layer deepening of up to 38 m day−1. Considerable variability was observed in the mixed layer response, indicating that lateral fluxes of heat and salt were also important. The magnitude and vertical distributions of these fluxes were quantified, and idealized mixed layer simulations suggest that their combined effect is a reduction in the mixed layer depth at the end of winter of up to several hundred meters.
Abstract
Mesoscale eddies are ubiquitous features of the global ocean circulation. Traditionally, anticyclonic eddies are thought to be associated with positive temperature anomalies while cyclonic eddies are associated with negative temperature anomalies. However, our recent study found that about one-fifth of the eddies identified from global satellite observations are cold-core anticyclonic eddies (CAEs) and warm-core cyclonic eddies (WCEs). Here we show that in the tropical oceans where the probabilities of CAEs and WCEs are high, there are significantly more CAEs and WCEs in summer than in winter. We conduct a suite of idealized numerical model experiments initialized with composite eddy structures obtained from Argo profiles as well as a heat budget analysis. The results highlight the key role of relative wind-stress-induced Ekman pumping, surface mixed layer depth, and vertical entrainment in the formation and seasonal cycle of these unconventional eddies. The relative wind stress is found to be particularly effective in converting conventional eddies into CAEs or WCEs when the surface mixed layer is shallow. The abundance of CAEs and WCEs in the global ocean calls for further research on this topic.
Abstract
Mesoscale eddies are ubiquitous features of the global ocean circulation. Traditionally, anticyclonic eddies are thought to be associated with positive temperature anomalies while cyclonic eddies are associated with negative temperature anomalies. However, our recent study found that about one-fifth of the eddies identified from global satellite observations are cold-core anticyclonic eddies (CAEs) and warm-core cyclonic eddies (WCEs). Here we show that in the tropical oceans where the probabilities of CAEs and WCEs are high, there are significantly more CAEs and WCEs in summer than in winter. We conduct a suite of idealized numerical model experiments initialized with composite eddy structures obtained from Argo profiles as well as a heat budget analysis. The results highlight the key role of relative wind-stress-induced Ekman pumping, surface mixed layer depth, and vertical entrainment in the formation and seasonal cycle of these unconventional eddies. The relative wind stress is found to be particularly effective in converting conventional eddies into CAEs or WCEs when the surface mixed layer is shallow. The abundance of CAEs and WCEs in the global ocean calls for further research on this topic.
Abstract
Recent studies suggest that the eddy kinetic energy is localized in the lee of significant topographic features in the Antarctic Circumpolar Current (ACC). Here we explore the importance of the local dynamics quantitatively using the outputs from the realistic ocean general circulation model hindcast with the aid of the modified Lorentz energy cycle. Results confirm the importance of energy transfer among reservoirs in the downstream region of standing meanders, showing that the major five standing meanders are responsible for more than 70% of the kinetic energy transfer to eddies and dissipation over the Antarctic Circumpolar Current region. The eddy kinetic energy is generated in the upper 3000-m depth downstream of the standing meanders and transported due to the vertical energy redistribution governed by the vertical pressure flux toward the deeper layer where the eddy energy is dissipated. Moreover, we also calculate the work done by the Ekman transport to confirm that the wind energy input works as the dominant energy source for the baroclinic energy pathway. The advantage of this quantity against the vertical mean density flux is that it is independent of the reference states defined arbitrarily. It is shown that the westerlies can supply sufficient energy locally to initiate baroclinic instability in the Indian and Pacific sectors of the ACC, whereas the nonlocal process is important in the Atlantic sector. Our results suggest that the five narrow regions associated with significant topography play key roles in the energy balance of the ACC region.
Significance Statement
The purpose of this study is to understand the eddy–mean flow interactions in the Antarctic Circumpolar Current from the energetic viewpoint. Our results show that the five narrow regions called “hotspots” in our study are responsible for the energy transfer from the mean flow to eddies. It is also found that the hotspots are important for the energy sink in the Southern Ocean. These findings suggest that the five hotspots are likely to play key roles in the responses of the Antarctic Circumpolar Current to the changes in westerlies in these decades.
Abstract
Recent studies suggest that the eddy kinetic energy is localized in the lee of significant topographic features in the Antarctic Circumpolar Current (ACC). Here we explore the importance of the local dynamics quantitatively using the outputs from the realistic ocean general circulation model hindcast with the aid of the modified Lorentz energy cycle. Results confirm the importance of energy transfer among reservoirs in the downstream region of standing meanders, showing that the major five standing meanders are responsible for more than 70% of the kinetic energy transfer to eddies and dissipation over the Antarctic Circumpolar Current region. The eddy kinetic energy is generated in the upper 3000-m depth downstream of the standing meanders and transported due to the vertical energy redistribution governed by the vertical pressure flux toward the deeper layer where the eddy energy is dissipated. Moreover, we also calculate the work done by the Ekman transport to confirm that the wind energy input works as the dominant energy source for the baroclinic energy pathway. The advantage of this quantity against the vertical mean density flux is that it is independent of the reference states defined arbitrarily. It is shown that the westerlies can supply sufficient energy locally to initiate baroclinic instability in the Indian and Pacific sectors of the ACC, whereas the nonlocal process is important in the Atlantic sector. Our results suggest that the five narrow regions associated with significant topography play key roles in the energy balance of the ACC region.
Significance Statement
The purpose of this study is to understand the eddy–mean flow interactions in the Antarctic Circumpolar Current from the energetic viewpoint. Our results show that the five narrow regions called “hotspots” in our study are responsible for the energy transfer from the mean flow to eddies. It is also found that the hotspots are important for the energy sink in the Southern Ocean. These findings suggest that the five hotspots are likely to play key roles in the responses of the Antarctic Circumpolar Current to the changes in westerlies in these decades.
Abstract
The energy and momentum balance of an abyssal overflow across a major sill in the Samoan Passage is estimated from two highly resolved towed sections, set 16 months apart, and results from a two-dimensional numerical simulation. Driven by the density anomaly across the sill, the flow is relatively steady. The system gains energy from divergence of horizontal pressure work
Abstract
The energy and momentum balance of an abyssal overflow across a major sill in the Samoan Passage is estimated from two highly resolved towed sections, set 16 months apart, and results from a two-dimensional numerical simulation. Driven by the density anomaly across the sill, the flow is relatively steady. The system gains energy from divergence of horizontal pressure work
Abstract
The high-resolution mooring observations reported here reveal a cascade process from internal solitary waves (ISWs) to turbulent mixing via high-frequency internal waves near the maximum local buoyancy frequency (near-N waves) in the deep water of the northern South China Sea (SCS). Riding on the parent ISW, near-N waves with a peak frequency of 20 cph emerged at the trough of the ISW and extended to the rear face of the ISW. Most of the near-N waves occurred around the thermocline, where the isothermal displacements induced by the near-N waves were largest with an amplitude of 12 m. The energy of near-N waves was 5% of that of the parent ISW, and instability investigations showed that due to the strong shear, Ri in the region of strong near-N waves was less than 1/4, suggesting that the near-N waves were unstable and might dissipate rapidly. Simulations based on the Korteweg–de Vries (KdV)–Burgers equation reproduced the formation of observed near-N waves due to the energy cascade from ISWs. Our observational results demonstrate a new energy cascade route from ISWs to turbulence in the deep water, deepening the understanding of the energy dissipation process of ISWs and their roles in the enhanced mixing in the northern SCS.
Abstract
The high-resolution mooring observations reported here reveal a cascade process from internal solitary waves (ISWs) to turbulent mixing via high-frequency internal waves near the maximum local buoyancy frequency (near-N waves) in the deep water of the northern South China Sea (SCS). Riding on the parent ISW, near-N waves with a peak frequency of 20 cph emerged at the trough of the ISW and extended to the rear face of the ISW. Most of the near-N waves occurred around the thermocline, where the isothermal displacements induced by the near-N waves were largest with an amplitude of 12 m. The energy of near-N waves was 5% of that of the parent ISW, and instability investigations showed that due to the strong shear, Ri in the region of strong near-N waves was less than 1/4, suggesting that the near-N waves were unstable and might dissipate rapidly. Simulations based on the Korteweg–de Vries (KdV)–Burgers equation reproduced the formation of observed near-N waves due to the energy cascade from ISWs. Our observational results demonstrate a new energy cascade route from ISWs to turbulence in the deep water, deepening the understanding of the energy dissipation process of ISWs and their roles in the enhanced mixing in the northern SCS.
Abstract
An existing approximately neutral surface, the ω surface, minimizes the neutrality error and hence also exhibits very small fictitious dianeutral diffusivity Df
that arises when lateral diffusion is applied along the surface, in nonneutral directions. However, there is also a spurious dianeutral advection that arises when lateral advection is applied nonneutrally along the surface; equivalently, lateral advection applied along the neutral tangent planes creates a vertical velocity e
sp through the ω surface. Mathematically, e
sp = u ⋅ s, where u is the lateral velocity and s is the slope error of the surface. We find that e
sp produces a leading-order term in the evolution equations of temperature and salinity, being similar in magnitude to the influence of cabbeling and thermobaricity. We introduce a new method to form an approximately neutral surface, called an ω
u
·
s
surface, that minimizes e
sp by adjusting its depth so that the slope error is nearly perpendicular to the lateral velocity. The e
sp on a surface cannot be reduced to zero when closed streamlines contain nonzero neutral helicity. While e
sp on the ω
u
·
s
surface is over 100 times smaller than that on the ω surface, the fictitious dianeutral diffusivity on the ω
u
·
s
surface is larger, nearly equal to the canonical 10−5 m2 s−1 background diffusivity. Thus, we also develop a method to minimize a combination of e
sp and Df
, yielding the
Abstract
An existing approximately neutral surface, the ω surface, minimizes the neutrality error and hence also exhibits very small fictitious dianeutral diffusivity Df
that arises when lateral diffusion is applied along the surface, in nonneutral directions. However, there is also a spurious dianeutral advection that arises when lateral advection is applied nonneutrally along the surface; equivalently, lateral advection applied along the neutral tangent planes creates a vertical velocity e
sp through the ω surface. Mathematically, e
sp = u ⋅ s, where u is the lateral velocity and s is the slope error of the surface. We find that e
sp produces a leading-order term in the evolution equations of temperature and salinity, being similar in magnitude to the influence of cabbeling and thermobaricity. We introduce a new method to form an approximately neutral surface, called an ω
u
·
s
surface, that minimizes e
sp by adjusting its depth so that the slope error is nearly perpendicular to the lateral velocity. The e
sp on a surface cannot be reduced to zero when closed streamlines contain nonzero neutral helicity. While e
sp on the ω
u
·
s
surface is over 100 times smaller than that on the ω surface, the fictitious dianeutral diffusivity on the ω
u
·
s
surface is larger, nearly equal to the canonical 10−5 m2 s−1 background diffusivity. Thus, we also develop a method to minimize a combination of e
sp and Df
, yielding the
Abstract
Simplified descriptions of the ocean are useful both for formulating explanatory theories and for conveying meaningful global attributes. Here, using a 26-yr average of a global state estimate from ECCO, the basis for Munk’s “abyssal recipes” is evaluated on a global scale between 1000- and 3000-m depth. The two specific hydrographic stations he used prove untypical, with potential temperature and salinity more generally displaying different vertical scale heights, and thus differing in one-dimensional (in the vertical) values of mixing coefficients and/or vertical velocities. The simplest explanation is that the circulation is fully three-dimensional with temperature and salinity fields not describable with a one-dimensional steady balance. In contrast, the potential density and buoyancy are quantitatively describable through a one-dimensional exponential balance, and which calls for an explanation in terms of turbulent mixing processes.
Abstract
Simplified descriptions of the ocean are useful both for formulating explanatory theories and for conveying meaningful global attributes. Here, using a 26-yr average of a global state estimate from ECCO, the basis for Munk’s “abyssal recipes” is evaluated on a global scale between 1000- and 3000-m depth. The two specific hydrographic stations he used prove untypical, with potential temperature and salinity more generally displaying different vertical scale heights, and thus differing in one-dimensional (in the vertical) values of mixing coefficients and/or vertical velocities. The simplest explanation is that the circulation is fully three-dimensional with temperature and salinity fields not describable with a one-dimensional steady balance. In contrast, the potential density and buoyancy are quantitatively describable through a one-dimensional exponential balance, and which calls for an explanation in terms of turbulent mixing processes.
Abstract
Subinertial, topographically trapped diurnal internal tides are an important energy source for turbulent mixing in the subarctic oceans. However, their generation may not be estimated by the conventional barotropic-to-baroclinic conversion because their vertical structure is sometimes barotropic, unlike superinertial internal tides that are always baroclinic. Here, a new energy diagram is presented, in which the barotropic mode is decomposed into the surface and topographic modes, with the latter being classified as part of the internal modes together with the baroclinic mode. The energy equation for the newly defined topographic mode is then derived, providing an appropriate formulation of the energy conversion rate from the subinertial surface tides to the topographically trapped internal tides. A series of numerical experiments confirm that the formulation successfully predicts the energy conversion rate for various cases, with the relative contribution of the baroclinic and topographic modes varying significantly depending on the bottom topography and stratification. Furthermore, this surface-to-internal conversion is demonstrated to give a significantly larger estimate than the barotropic-to-baroclinic conversion for subinertial tides. Applying the formulation to the results of a realistic numerical simulation in the Kuril Straits, an area with the strongest mixing due to subinertial diurnal tides, shows that the surface mode is converted into the baroclinic and topographic modes with comparable magnitudes, responsible for most of the energy dissipation in this area. These results indicate the need to reestimate the global distribution of the generation rate of the subinertial internal tides using our new formulation and to clarify their dissipation mechanisms.
Significance Statement
Diurnal internal tides in mid- to high-latitude oceans are very different from semidiurnal internal tides, in that they are trapped by topographic features and characterized by uniform rotation throughout the water column rather than by vertical oscillations within the water column. Focusing on this character, we formulate for the first time the generation rate of diurnal internal tides trapped over variable bottom topography. A series of idealized numerical experiments and application to the Kuril Straits show the validity and usefulness of this formulation, which provides a significantly larger generation estimate than previous studies. The results of this study are important for accurate global mapping of turbulent mixing induced by the breaking of internal tides.
Abstract
Subinertial, topographically trapped diurnal internal tides are an important energy source for turbulent mixing in the subarctic oceans. However, their generation may not be estimated by the conventional barotropic-to-baroclinic conversion because their vertical structure is sometimes barotropic, unlike superinertial internal tides that are always baroclinic. Here, a new energy diagram is presented, in which the barotropic mode is decomposed into the surface and topographic modes, with the latter being classified as part of the internal modes together with the baroclinic mode. The energy equation for the newly defined topographic mode is then derived, providing an appropriate formulation of the energy conversion rate from the subinertial surface tides to the topographically trapped internal tides. A series of numerical experiments confirm that the formulation successfully predicts the energy conversion rate for various cases, with the relative contribution of the baroclinic and topographic modes varying significantly depending on the bottom topography and stratification. Furthermore, this surface-to-internal conversion is demonstrated to give a significantly larger estimate than the barotropic-to-baroclinic conversion for subinertial tides. Applying the formulation to the results of a realistic numerical simulation in the Kuril Straits, an area with the strongest mixing due to subinertial diurnal tides, shows that the surface mode is converted into the baroclinic and topographic modes with comparable magnitudes, responsible for most of the energy dissipation in this area. These results indicate the need to reestimate the global distribution of the generation rate of the subinertial internal tides using our new formulation and to clarify their dissipation mechanisms.
Significance Statement
Diurnal internal tides in mid- to high-latitude oceans are very different from semidiurnal internal tides, in that they are trapped by topographic features and characterized by uniform rotation throughout the water column rather than by vertical oscillations within the water column. Focusing on this character, we formulate for the first time the generation rate of diurnal internal tides trapped over variable bottom topography. A series of idealized numerical experiments and application to the Kuril Straits show the validity and usefulness of this formulation, which provides a significantly larger generation estimate than previous studies. The results of this study are important for accurate global mapping of turbulent mixing induced by the breaking of internal tides.
Abstract
The deep channel north of New Guinea (NG) is the choke site for the upper deep branches of the Pacific meridional overturning circulation (U-PMOC). The U-PMOC is a crucial element of the ocean’s climate and biogeochemical systems. It carries the mixed water of the Upper Circumpolar Water and North Pacific Deep Water with a potential temperature over 1.2°–2.2°C. The pathway and volume transport of U-PMOC through the deep channel north of NG are revealed by mooring measurements from 2014 to 2019. Mean U-PMOC is located at ∼2000–3500 m with a velocity core at 2550 m and is directed eastward. The U-PMOC shows a strong seasonal variability with a direction reversal from June to September. The oceanic reanalysis product GLORYS12V1 well reproduces the observed U-PMOC and is thus used to estimate the mean and standard deviation of U-PMOC’s volume transport as 2.19 ± 11.4 Sv (1 Sv ≡ 106 m3 s−1) and to explore the underlying dynamics of the U-PMOC. The seasonality of U-PMOC is induced by the vertical propagation of the Rossby energy through the upper ocean in the eastern Pacific to the deep ocean in the western Pacific. The mean eastward U-PMOC transport is forced by the zonal deep pressure gradient, which is mainly determined by the local upper-ocean processes above 500 m.
Abstract
The deep channel north of New Guinea (NG) is the choke site for the upper deep branches of the Pacific meridional overturning circulation (U-PMOC). The U-PMOC is a crucial element of the ocean’s climate and biogeochemical systems. It carries the mixed water of the Upper Circumpolar Water and North Pacific Deep Water with a potential temperature over 1.2°–2.2°C. The pathway and volume transport of U-PMOC through the deep channel north of NG are revealed by mooring measurements from 2014 to 2019. Mean U-PMOC is located at ∼2000–3500 m with a velocity core at 2550 m and is directed eastward. The U-PMOC shows a strong seasonal variability with a direction reversal from June to September. The oceanic reanalysis product GLORYS12V1 well reproduces the observed U-PMOC and is thus used to estimate the mean and standard deviation of U-PMOC’s volume transport as 2.19 ± 11.4 Sv (1 Sv ≡ 106 m3 s−1) and to explore the underlying dynamics of the U-PMOC. The seasonality of U-PMOC is induced by the vertical propagation of the Rossby energy through the upper ocean in the eastern Pacific to the deep ocean in the western Pacific. The mean eastward U-PMOC transport is forced by the zonal deep pressure gradient, which is mainly determined by the local upper-ocean processes above 500 m.
Abstract
The spectral description of the energy of oceanic internal gravity waves is generally represented by the Garrett–Munk (GM) model, a function with a power-law decrease of spectral energy in wavenumber–frequency space. Besides the slopes of these power laws, the spectrum is expressed as a function of energy and a bandwidth parameter that fixes the range of vertical modes excited in the respective state. Whereas concepts have been developed and agreed upon of what processes feed the wave spectrum and what dissipates energy, there is no explanation of what shapes the spectral distribution, i.e., how the power laws come about and what sets the bandwidth. The present study develops a parametric spectral model of energy and bandwidth from the basic underlying energy balance in terms of forcing, propagation, refraction, spectral transfer, and dissipation. The model is an extension of the IDEMIX (Internal Wave Dissipation, Energy and Mixing) models where bandwidth was taken as a constant parameter. The current version of the model is restricted to single-column mode and the slopes of the spectral power laws are fixed. A coupled system of predictive equations for energy and bandwidth (for up- and downward propagating waves) results. The equations imply that bandwidth relates to energy by a power law with an exponent given by the dynamical parameters. It agrees favorably with energy, bandwidth, and slope data from previously published fits of the GM model to Argo float observations. Numerical solutions of the coupled energy–bandwidth model in stand-alone modus are presented.
Abstract
The spectral description of the energy of oceanic internal gravity waves is generally represented by the Garrett–Munk (GM) model, a function with a power-law decrease of spectral energy in wavenumber–frequency space. Besides the slopes of these power laws, the spectrum is expressed as a function of energy and a bandwidth parameter that fixes the range of vertical modes excited in the respective state. Whereas concepts have been developed and agreed upon of what processes feed the wave spectrum and what dissipates energy, there is no explanation of what shapes the spectral distribution, i.e., how the power laws come about and what sets the bandwidth. The present study develops a parametric spectral model of energy and bandwidth from the basic underlying energy balance in terms of forcing, propagation, refraction, spectral transfer, and dissipation. The model is an extension of the IDEMIX (Internal Wave Dissipation, Energy and Mixing) models where bandwidth was taken as a constant parameter. The current version of the model is restricted to single-column mode and the slopes of the spectral power laws are fixed. A coupled system of predictive equations for energy and bandwidth (for up- and downward propagating waves) results. The equations imply that bandwidth relates to energy by a power law with an exponent given by the dynamical parameters. It agrees favorably with energy, bandwidth, and slope data from previously published fits of the GM model to Argo float observations. Numerical solutions of the coupled energy–bandwidth model in stand-alone modus are presented.