# Search Results

## Abstract

The spontaneous generation of inertia–gravity waves in idealized life cycles of baroclinic instability is investigated using the Weather Research and Forecasting Model. Two substantially different life cycles of baroclinic instability are obtained by varying the initial zonal jet. The wave generation depends strongly on the details of the baroclinic wave’s development. In the life cycle dominated by cyclonic behavior, the most conspicuous gravity waves are excited by the upper-level jet and are broadly consistent with previous simulations of O’Sullivan and Dunkerton. In the life cycle that is dominated by anticyclonic behavior, the most conspicuous gravity waves even in the stratosphere are excited by the surface fronts, although the fronts are no stronger than in the cyclonic life cycle. The anticyclonic life cycle also reveals waves in the lower stratosphere above the upper-level trough of the baroclinic wave; these waves have not been previously identified in idealized simulations. The sensitivities of the different waves to both resolution and dissipation are discussed.

## Abstract

The spontaneous generation of inertia–gravity waves in idealized life cycles of baroclinic instability is investigated using the Weather Research and Forecasting Model. Two substantially different life cycles of baroclinic instability are obtained by varying the initial zonal jet. The wave generation depends strongly on the details of the baroclinic wave’s development. In the life cycle dominated by cyclonic behavior, the most conspicuous gravity waves are excited by the upper-level jet and are broadly consistent with previous simulations of O’Sullivan and Dunkerton. In the life cycle that is dominated by anticyclonic behavior, the most conspicuous gravity waves even in the stratosphere are excited by the surface fronts, although the fronts are no stronger than in the cyclonic life cycle. The anticyclonic life cycle also reveals waves in the lower stratosphere above the upper-level trough of the baroclinic wave; these waves have not been previously identified in idealized simulations. The sensitivities of the different waves to both resolution and dissipation are discussed.

## Abstract

The gravity waves (GWs) produced by three-dimensional potential vorticity (PV) anomalies are examined under the assumption of constant vertical shear, constant stratification, and unbounded domain. As in the two-dimensional case analyzed in an earlier paper, the disturbance near the PV anomaly is well modeled by quasigeostrophic theory. At larger distances the nature of the disturbance changes across the two inertial layers that are located above and below the anomaly, and it takes the form of a vertically propagating GW beyond these.

For a horizontally monochromatic PV anomaly of infinitesimal depth, the disturbance is described analytically using both an exact solution and a WKB approximation; the latter includes an exponentially small term that captures the change of the solution near the PV anomaly induced by the radiation boundary condition in the far field. The analytical results reveal a strong sensitivity of the emission to the Richardson number and to the orientation of the horizontal wavenumber: the absorptive properties of the inertial layers are such that the emission is maximized in the Northern Hemisphere for wavenumbers at negative angles to the shear.

For localized PV anomalies, numerical computations give the temporal evolution of the GW field. Analytical and numerical results are also used to establish an explicit form for the Eliassen–Palm flux that could be used to parameterize GW sources in GCMs. The properties of the Eliassen–Palm flux vector imply that in a westerly shear, the GWs exert a drag in a southwest direction in the upper inertial layer, and in a northwest direction at the altitudes where the GWs dissipate aloft.

## Abstract

The gravity waves (GWs) produced by three-dimensional potential vorticity (PV) anomalies are examined under the assumption of constant vertical shear, constant stratification, and unbounded domain. As in the two-dimensional case analyzed in an earlier paper, the disturbance near the PV anomaly is well modeled by quasigeostrophic theory. At larger distances the nature of the disturbance changes across the two inertial layers that are located above and below the anomaly, and it takes the form of a vertically propagating GW beyond these.

For a horizontally monochromatic PV anomaly of infinitesimal depth, the disturbance is described analytically using both an exact solution and a WKB approximation; the latter includes an exponentially small term that captures the change of the solution near the PV anomaly induced by the radiation boundary condition in the far field. The analytical results reveal a strong sensitivity of the emission to the Richardson number and to the orientation of the horizontal wavenumber: the absorptive properties of the inertial layers are such that the emission is maximized in the Northern Hemisphere for wavenumbers at negative angles to the shear.

For localized PV anomalies, numerical computations give the temporal evolution of the GW field. Analytical and numerical results are also used to establish an explicit form for the Eliassen–Palm flux that could be used to parameterize GW sources in GCMs. The properties of the Eliassen–Palm flux vector imply that in a westerly shear, the GWs exert a drag in a southwest direction in the upper inertial layer, and in a northwest direction at the altitudes where the GWs dissipate aloft.

## Abstract

The gravity waves (GWs) generated by potential vorticity (PV) anomalies in a rotating stratified shear flow are examined under the assumptions of constant vertical shear, two-dimensionality, and unbounded domain. Near a PV anomaly, the associated perturbation is well modeled by quasigeostrophic theory. This is not the case at large vertical distances, however, and in particular beyond the two inertial layers that appear above and below the anomaly; there, the perturbation consists of vertically propagating gravity waves. This structure is described analytically, using an expansion in the continuous spectrum of the singular modes that results from the presence of critical levels.

Several explicit results are obtained. These include the form of the Eliassen–Palm (EP) flux as a function of the Richardson number *N* ^{2}/Λ^{2}, where *N* is the Brunt–Väisälä frequency and Λ the vertical shear. Its nondimensional value is shown to be approximately exp(−*πN*/Λ)/8 in the far-field GW region, approximately twice that between the two inertial layers. These results, which imply substantial wave–flow interactions in the inertial layers, are valid for Richardson numbers larger than 1 and for a large range of PV distributions. In dimensional form they provide simple relationships between the EP fluxes and the large-scale flow characteristics.

As an illustration, the authors consider a PV disturbance with an amplitude of 1 PVU and a depth of 1 km, and estimate that the associated EP flux ranges between 0.1 and 100 mPa for a Richardson number between 1 and 10. These values of the flux are comparable with those observed in the lower stratosphere, which suggests that the mechanism identified in this paper provides a substantial gravity wave source, one that could be parameterized in GCMs.

## Abstract

The gravity waves (GWs) generated by potential vorticity (PV) anomalies in a rotating stratified shear flow are examined under the assumptions of constant vertical shear, two-dimensionality, and unbounded domain. Near a PV anomaly, the associated perturbation is well modeled by quasigeostrophic theory. This is not the case at large vertical distances, however, and in particular beyond the two inertial layers that appear above and below the anomaly; there, the perturbation consists of vertically propagating gravity waves. This structure is described analytically, using an expansion in the continuous spectrum of the singular modes that results from the presence of critical levels.

Several explicit results are obtained. These include the form of the Eliassen–Palm (EP) flux as a function of the Richardson number *N* ^{2}/Λ^{2}, where *N* is the Brunt–Väisälä frequency and Λ the vertical shear. Its nondimensional value is shown to be approximately exp(−*πN*/Λ)/8 in the far-field GW region, approximately twice that between the two inertial layers. These results, which imply substantial wave–flow interactions in the inertial layers, are valid for Richardson numbers larger than 1 and for a large range of PV distributions. In dimensional form they provide simple relationships between the EP fluxes and the large-scale flow characteristics.

As an illustration, the authors consider a PV disturbance with an amplitude of 1 PVU and a depth of 1 km, and estimate that the associated EP flux ranges between 0.1 and 100 mPa for a Richardson number between 1 and 10. These values of the flux are comparable with those observed in the lower stratosphere, which suggests that the mechanism identified in this paper provides a substantial gravity wave source, one that could be parameterized in GCMs.

## Abstract

Normal modes of a linear vertical shear (Eady shear) are studied within the linearized primitive equations for a rotating stratified fluid above a rigid lower boundary. The authors' interest is in modes having an inertial critical layer present at some height within the flow. Below this layer, the solutions can be closely approximated by balanced edge waves obtained through an asymptotic expansion in Rossby number. Above, the solutions behave as gravity waves. Hence these modes are an example of a spatial coupling of balanced motions to gravity waves.

The amplitude of the gravity waves relative to the balanced part of the solutions is obtained analytically and numerically as a function of parameters. It is shown that the waves are exponentially small in Rossby number. Moreover, their amplitude depends in a nontrivial way on the meridional wavenumber. For modes having a radiating upper boundary condition, the meridional wavenumber for which the gravity wave amplitude is maximal occurs when the tilts of the balanced edge wave and gravity waves agree.

## Abstract

Normal modes of a linear vertical shear (Eady shear) are studied within the linearized primitive equations for a rotating stratified fluid above a rigid lower boundary. The authors' interest is in modes having an inertial critical layer present at some height within the flow. Below this layer, the solutions can be closely approximated by balanced edge waves obtained through an asymptotic expansion in Rossby number. Above, the solutions behave as gravity waves. Hence these modes are an example of a spatial coupling of balanced motions to gravity waves.

The amplitude of the gravity waves relative to the balanced part of the solutions is obtained analytically and numerically as a function of parameters. It is shown that the waves are exponentially small in Rossby number. Moreover, their amplitude depends in a nontrivial way on the meridional wavenumber. For modes having a radiating upper boundary condition, the meridional wavenumber for which the gravity wave amplitude is maximal occurs when the tilts of the balanced edge wave and gravity waves agree.

## Abstract

Previous simulations of dipole vortices propagating through rotating, stratified fluid have revealed small-scale inertia–gravity waves that are embedded within the dipole near its leading edge and are approximately stationary relative to the dipole. The mechanism by which these waves are generated is investigated, beginning from the observation that the dipole can be reasonably approximated by a balanced quasigeostrophic (QG) solution. The deviations from the QG solution (including the waves) then satisfy linear equations that come from linearization of the governing equations about the QG dipole and are forced by the residual tendency of the QG dipole (i.e., the difference between the time tendency of the QG solution and that of the full primitive equations initialized with the QG fields). The waves do not appear to be generated by an instability of the balanced dipole, as homogeneous solutions of the linear equations amplify little over the time scale for which the linear equations are valid. Linear solutions forced by the residual tendency capture the scale, location, and pattern of the inertia–gravity waves, although they overpredict the wave amplitude by a factor of 2. There is thus strong evidence that the waves are generated as a forced linear response to the balanced flow. The relation to and differences from other theories for wave generation by balanced flows, including those of Lighthill and Ford et al., are discussed.

## Abstract

Previous simulations of dipole vortices propagating through rotating, stratified fluid have revealed small-scale inertia–gravity waves that are embedded within the dipole near its leading edge and are approximately stationary relative to the dipole. The mechanism by which these waves are generated is investigated, beginning from the observation that the dipole can be reasonably approximated by a balanced quasigeostrophic (QG) solution. The deviations from the QG solution (including the waves) then satisfy linear equations that come from linearization of the governing equations about the QG dipole and are forced by the residual tendency of the QG dipole (i.e., the difference between the time tendency of the QG solution and that of the full primitive equations initialized with the QG fields). The waves do not appear to be generated by an instability of the balanced dipole, as homogeneous solutions of the linear equations amplify little over the time scale for which the linear equations are valid. Linear solutions forced by the residual tendency capture the scale, location, and pattern of the inertia–gravity waves, although they overpredict the wave amplitude by a factor of 2. There is thus strong evidence that the waves are generated as a forced linear response to the balanced flow. The relation to and differences from other theories for wave generation by balanced flows, including those of Lighthill and Ford et al., are discussed.

## Abstract

In this article, long-duration balloon and spaceborne observations, and mesoscale numerical simulations are used to study the intermittency of gravity waves in the lower stratosphere above Antarctica and the Southern Ocean; namely, the characteristics of the gravity wave momentum-flux probability density functions (pdfs) obtained with these three datasets are described. The pdfs consistently exhibit long tails associated with the occurrence of rare and large-amplitude events. The pdf tails are even longer above mountains than above oceanic areas, which is in agreement with previous studies of gravity wave intermittency in this region. It is moreover found that these rare, large-amplitude events represent the main contribution to the total momentum flux during the winter regime of the stratospheric circulation. In contrast, the wave intermittency significantly decreases when stratospheric easterlies develop in late spring and summer. It is also shown that, except above mountainous areas in winter, the momentum-flux pdfs tend to behave like lognormal distributions. Monte Carlo simulations are undertaken to examine the role played by critical levels in influencing the shape of momentum-flux pdfs. In particular, the study finds that the lognormal shape may result from the propagation of a wave spectrum into a varying background wind field that generates the occurrence of frequent critical levels.

## Abstract

In this article, long-duration balloon and spaceborne observations, and mesoscale numerical simulations are used to study the intermittency of gravity waves in the lower stratosphere above Antarctica and the Southern Ocean; namely, the characteristics of the gravity wave momentum-flux probability density functions (pdfs) obtained with these three datasets are described. The pdfs consistently exhibit long tails associated with the occurrence of rare and large-amplitude events. The pdf tails are even longer above mountains than above oceanic areas, which is in agreement with previous studies of gravity wave intermittency in this region. It is moreover found that these rare, large-amplitude events represent the main contribution to the total momentum flux during the winter regime of the stratospheric circulation. In contrast, the wave intermittency significantly decreases when stratospheric easterlies develop in late spring and summer. It is also shown that, except above mountainous areas in winter, the momentum-flux pdfs tend to behave like lognormal distributions. Monte Carlo simulations are undertaken to examine the role played by critical levels in influencing the shape of momentum-flux pdfs. In particular, the study finds that the lognormal shape may result from the propagation of a wave spectrum into a varying background wind field that generates the occurrence of frequent critical levels.

## 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.

## 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.

## Abstract

Quantification of inertia–gravity waves (IGWs) generated by upper-level jet–surface front systems and their parameterization in global models of the atmosphere relies on suitable methods to estimate the strength of IGWs. A harmonic divergence analysis (HDA) that has been previously employed for quantification of IGWs combines wave properties from linear dynamics with a sophisticated statistical analysis to provide such estimates. A question of fundamental importance that arises is how the measures of IGW activity provided by the HDA are related to the measures coming from the wave–vortex decomposition (WVD) methods. The question is addressed by employing the nonlinear balance relations of the first-order *δ*–*γ*, the Bolin–Charney, and the first- to third-order Rossby number expansion to carry out WVD. The global kinetic energy of IGWs given by the HDA and WVD are compared in numerical simulations of moist baroclinic waves by the Weather Research and Forecasting (WRF) Model in a channel on the *f* plane. The estimates of the HDA are found to be 2–3 times smaller than those of the optimal WVD. This is in part due to the absence of a well-defined scale separation between the waves and vortical flows, the IGW estimates by the HDA capturing only the dominant wave packets and with limited scales. It is also shown that the difference between the HDA and WVD estimates is related to the width of the IGW spectrum.

## Abstract

Quantification of inertia–gravity waves (IGWs) generated by upper-level jet–surface front systems and their parameterization in global models of the atmosphere relies on suitable methods to estimate the strength of IGWs. A harmonic divergence analysis (HDA) that has been previously employed for quantification of IGWs combines wave properties from linear dynamics with a sophisticated statistical analysis to provide such estimates. A question of fundamental importance that arises is how the measures of IGW activity provided by the HDA are related to the measures coming from the wave–vortex decomposition (WVD) methods. The question is addressed by employing the nonlinear balance relations of the first-order *δ*–*γ*, the Bolin–Charney, and the first- to third-order Rossby number expansion to carry out WVD. The global kinetic energy of IGWs given by the HDA and WVD are compared in numerical simulations of moist baroclinic waves by the Weather Research and Forecasting (WRF) Model in a channel on the *f* plane. The estimates of the HDA are found to be 2–3 times smaller than those of the optimal WVD. This is in part due to the absence of a well-defined scale separation between the waves and vortical flows, the IGW estimates by the HDA capturing only the dominant wave packets and with limited scales. It is also shown that the difference between the HDA and WVD estimates is related to the width of the IGW spectrum.

## Abstract

Machine learning (ML) provides a powerful tool for investigating the relationship between the large-scale flow and unresolved processes, which need to be parameterized in climate models. The current work explores the performance of the random forest regressor (RF) as a nonparametric model in the reconstruction of nonorographic gravity waves (GWs) over midlatitude oceanic areas. The ERA5 dataset from the European Centre for Medium-Range Weather Forecasts (ECMWF) model outputs is employed in its full resolution to derive GW variations in the lower stratosphere. Coarse-grained variables in a column-based configuration of the atmosphere are used to reconstruct the GWs variability at the target level. The first important outcome is the relative success in reconstructing the GW signal (coefficient of determination *R*
^{2} ≈ 0.85 for “E3” combination). The second outcome is that the most informative explanatory variable is the local background wind speed. This questions the traditional framework of gravity wave parameterizations, for which, at these heights, one would expect more sensitivity to sources below than to local flow. Finally, to test the efficiency of a relatively simple, parametric statistical model, the efficiency of linear regression was compared to that of random forests with a restricted set of only five explanatory variables. Results were poor. Increasing the number of input variables to 15 hardly changes the performance of the linear regression (*R*
^{2} changes slightly from 0.18 to 0.21), while it leads to better results with the random forests (*R*
^{2} increases from 0.29 to 0.37).

## Abstract

Machine learning (ML) provides a powerful tool for investigating the relationship between the large-scale flow and unresolved processes, which need to be parameterized in climate models. The current work explores the performance of the random forest regressor (RF) as a nonparametric model in the reconstruction of nonorographic gravity waves (GWs) over midlatitude oceanic areas. The ERA5 dataset from the European Centre for Medium-Range Weather Forecasts (ECMWF) model outputs is employed in its full resolution to derive GW variations in the lower stratosphere. Coarse-grained variables in a column-based configuration of the atmosphere are used to reconstruct the GWs variability at the target level. The first important outcome is the relative success in reconstructing the GW signal (coefficient of determination *R*
^{2} ≈ 0.85 for “E3” combination). The second outcome is that the most informative explanatory variable is the local background wind speed. This questions the traditional framework of gravity wave parameterizations, for which, at these heights, one would expect more sensitivity to sources below than to local flow. Finally, to test the efficiency of a relatively simple, parametric statistical model, the efficiency of linear regression was compared to that of random forests with a restricted set of only five explanatory variables. Results were poor. Increasing the number of input variables to 15 hardly changes the performance of the linear regression (*R*
^{2} changes slightly from 0.18 to 0.21), while it leads to better results with the random forests (*R*
^{2} increases from 0.29 to 0.37).

## Abstract

The parameterization of inertia–gravity waves (IGWs) is of considerable importance in general circulation models. Among the challenging issues faced in studies concerned with parameterization of IGWs is the estimation of diabatic forcing in a way independent of the physics parameterization schemes, in particular, convection. The requirement is to estimate the diabatic heating associated with balanced motion. This can be done by comparing estimates of balanced vertical motion with and without diabatic effects. The omega equation provides the natural method of estimating balanced vertical motion without diabatic effects, and several methods for including diabatic effects are compared. To this end, the assumption of spatial-scale separation between IGWs and balanced flows is combined with a suitable form of the balanced omega equation. To test the methods constructed for estimating diabatic heating, an idealized numerical simulation of the moist baroclinic waves is performed using the Weather Research and Forecasting (WRF) Model in a channel on the *f* plane. In overall agreement with the diabatic heating of the WRF Model, in the omega-equation-based estimates, the maxima of heating appear in the warm sector of the baroclinic wave and in the exit region of the upper-level jet. The omega-equation-based method with spatial smoothing for estimating balanced vertical motion is thus presented as the proper way to evaluate diabatic forcing for parameterization of IGWs.

## Abstract

The parameterization of inertia–gravity waves (IGWs) is of considerable importance in general circulation models. Among the challenging issues faced in studies concerned with parameterization of IGWs is the estimation of diabatic forcing in a way independent of the physics parameterization schemes, in particular, convection. The requirement is to estimate the diabatic heating associated with balanced motion. This can be done by comparing estimates of balanced vertical motion with and without diabatic effects. The omega equation provides the natural method of estimating balanced vertical motion without diabatic effects, and several methods for including diabatic effects are compared. To this end, the assumption of spatial-scale separation between IGWs and balanced flows is combined with a suitable form of the balanced omega equation. To test the methods constructed for estimating diabatic heating, an idealized numerical simulation of the moist baroclinic waves is performed using the Weather Research and Forecasting (WRF) Model in a channel on the *f* plane. In overall agreement with the diabatic heating of the WRF Model, in the omega-equation-based estimates, the maxima of heating appear in the warm sector of the baroclinic wave and in the exit region of the upper-level jet. The omega-equation-based method with spatial smoothing for estimating balanced vertical motion is thus presented as the proper way to evaluate diabatic forcing for parameterization of IGWs.