Browse
Abstract
Internal gravity waves (GWs) are ubiquitous in the atmosphere, making significant contributions to the mesoscale motions. Since the majority of their spectrum is unresolved in global circulation models, their effects need to be parameterized. In recent decades GWs have been increasingly studied in high-resolution simulations, which, unlike direct observations, allow us to explore full spatiotemporal variations of the resolved wave field. In our study we analyze and refine a traditional method for GW analysis in a high-resolution simulation on a regional domain around the Drake Passage. We show that GW momentum drag estimates based on the Gaussian high-pass filter method applied to separate GW perturbations from the background are sensitive to the choice of a cutoff parameter. The impact of the cutoff parameter is higher for horizontal fluxes of horizontal momentum, which indicates higher sensitivity for horizontally propagating waves. Two modified methods, which choose the parameter value from spectral information, are proposed. The dynamically determined cutoff is mostly higher than the traditional cutoff values around 500 km, leading to larger GW fluxes and drag, and varies with time and altitude. The differences between the traditional and the modified methods are especially pronounced during events with significant drag contributions from horizontal momentum fluxes.
Significance Statement
In this study, we highlight that the analysis of gravity wave activity from high-resolution datasets is a complex task with a pronounced sensitivity to the methodology, and we propose modified versions of a classical statistical gravity wave detection method enhanced by the spectral information. Although no optimal methodology exists to date, we show that the modified methods improve the accuracy of the gravity wave activity estimates, especially when oblique propagation plays a role.
Abstract
Internal gravity waves (GWs) are ubiquitous in the atmosphere, making significant contributions to the mesoscale motions. Since the majority of their spectrum is unresolved in global circulation models, their effects need to be parameterized. In recent decades GWs have been increasingly studied in high-resolution simulations, which, unlike direct observations, allow us to explore full spatiotemporal variations of the resolved wave field. In our study we analyze and refine a traditional method for GW analysis in a high-resolution simulation on a regional domain around the Drake Passage. We show that GW momentum drag estimates based on the Gaussian high-pass filter method applied to separate GW perturbations from the background are sensitive to the choice of a cutoff parameter. The impact of the cutoff parameter is higher for horizontal fluxes of horizontal momentum, which indicates higher sensitivity for horizontally propagating waves. Two modified methods, which choose the parameter value from spectral information, are proposed. The dynamically determined cutoff is mostly higher than the traditional cutoff values around 500 km, leading to larger GW fluxes and drag, and varies with time and altitude. The differences between the traditional and the modified methods are especially pronounced during events with significant drag contributions from horizontal momentum fluxes.
Significance Statement
In this study, we highlight that the analysis of gravity wave activity from high-resolution datasets is a complex task with a pronounced sensitivity to the methodology, and we propose modified versions of a classical statistical gravity wave detection method enhanced by the spectral information. Although no optimal methodology exists to date, we show that the modified methods improve the accuracy of the gravity wave activity estimates, especially when oblique propagation plays a role.
Abstract
Several studies have reported vertical kinetic energy spectra almost white in horizontal wavenumber space with evidence of two maxima at synoptic scales and mesoscales, leaving the explanation of these maxima open. Processes known to influence the shape of the horizontal kinetic energy spectra include the superposition of quasi-linear inertia–gravity waves (IGWs), quasigeostrophic turbulence, and moist convection. In contrast, vertical kinetic energy has been discussed much less, as measuring vertical velocity remains challenging. This study compares the horizontal and vertical kinetic energy spectra and their relationships in global storm-resolving simulations from the DYAMOND experiment. The consistency of these relationships with linear IGW theory is tested by diagnosing horizontal wind fluctuations associated with IGW modes. Furthermore, it is shown that hydrostatic IGW polarization relations provide a quantitative prediction of the spectral slopes of vertical kinetic energy at large scales and mesoscales, where the intrinsic frequencies are inferred from the linearized vorticity equation. Our results suggest that IGW modes dominate the vertical kinetic energy spectra at most horizontal scales, whereas an incompressible, isotropic scaling of the continuity equation captures the relationship between horizontal and vertical kinetic energy spectra at small scales.
Abstract
Several studies have reported vertical kinetic energy spectra almost white in horizontal wavenumber space with evidence of two maxima at synoptic scales and mesoscales, leaving the explanation of these maxima open. Processes known to influence the shape of the horizontal kinetic energy spectra include the superposition of quasi-linear inertia–gravity waves (IGWs), quasigeostrophic turbulence, and moist convection. In contrast, vertical kinetic energy has been discussed much less, as measuring vertical velocity remains challenging. This study compares the horizontal and vertical kinetic energy spectra and their relationships in global storm-resolving simulations from the DYAMOND experiment. The consistency of these relationships with linear IGW theory is tested by diagnosing horizontal wind fluctuations associated with IGW modes. Furthermore, it is shown that hydrostatic IGW polarization relations provide a quantitative prediction of the spectral slopes of vertical kinetic energy at large scales and mesoscales, where the intrinsic frequencies are inferred from the linearized vorticity equation. Our results suggest that IGW modes dominate the vertical kinetic energy spectra at most horizontal scales, whereas an incompressible, isotropic scaling of the continuity equation captures the relationship between horizontal and vertical kinetic energy spectra at small scales.
Abstract
Four state-of-the-science numerical weather prediction (NWP) models were used to perform mountain wave (MW)-resolving hindcasts over the Drake Passage of a 10-day period in 2010 with numerous observed MW cases. The Integrated Forecast System (IFS) and the Icosahedral Nonhydrostatic (ICON) model were run at Δx ≈ 9 and 13 km globally. The Weather Research and Forecasting (WRF) Model and the Met Office Unified Model (UM) were both configured with a Δx = 3-km regional domain. All domains had tops near 1 Pa (z ≈ 80 km). These deep domains allowed quantitative validation against Atmospheric Infrared Sounder (AIRS) observations, accounting for observation time, viewing geometry, and radiative transfer. All models reproduced observed middle-atmosphere MWs with remarkable skill. Increased horizontal resolution improved validations. Still, all models underrepresented observed MW amplitudes, even after accounting for model effective resolution and instrument noise, suggesting even at Δx ≈ 3-km resolution, small-scale MWs are underresolved and/or overdiffused. MW drag parameterizations are still necessary in NWP models at current operational resolutions of Δx ≈ 10 km. Upper GW sponge layers in the operationally configured models significantly, artificially reduced MW amplitudes in the upper stratosphere and mesosphere. In the IFS, parameterized GW drags partly compensated this deficiency, but still, total drags were ≈6 times smaller than that resolved at Δx ≈ 3 km. Meridionally propagating MWs significantly enhance zonal drag over the Drake Passage. Interestingly, drag associated with meridional fluxes of zonal momentum (i.e.,
Significance Statement
This study had three purposes: to quantitatively evaluate how well four state-of-the-science weather models could reproduce observed mountain waves (MWs) in the middle atmosphere, to compare the simulated MWs within the models, and to quantitatively evaluate two MW parameterizations in a widely used climate model. These models reproduced observed MWs with remarkable skill. Still, MW parameterizations are necessary in current Δx ≈ 10-km resolution global weather models. Even Δx ≈ 3-km resolution does not appear to be high enough to represent all momentum-fluxing MW scales. Meridionally propagating MWs can significantly influence zonal winds over the Drake Passage. Parameterizations that handle horizontal propagation may need to consider horizontal fluxes of horizontal momentum in order to get the direction of their forcing correct.
Abstract
Four state-of-the-science numerical weather prediction (NWP) models were used to perform mountain wave (MW)-resolving hindcasts over the Drake Passage of a 10-day period in 2010 with numerous observed MW cases. The Integrated Forecast System (IFS) and the Icosahedral Nonhydrostatic (ICON) model were run at Δx ≈ 9 and 13 km globally. The Weather Research and Forecasting (WRF) Model and the Met Office Unified Model (UM) were both configured with a Δx = 3-km regional domain. All domains had tops near 1 Pa (z ≈ 80 km). These deep domains allowed quantitative validation against Atmospheric Infrared Sounder (AIRS) observations, accounting for observation time, viewing geometry, and radiative transfer. All models reproduced observed middle-atmosphere MWs with remarkable skill. Increased horizontal resolution improved validations. Still, all models underrepresented observed MW amplitudes, even after accounting for model effective resolution and instrument noise, suggesting even at Δx ≈ 3-km resolution, small-scale MWs are underresolved and/or overdiffused. MW drag parameterizations are still necessary in NWP models at current operational resolutions of Δx ≈ 10 km. Upper GW sponge layers in the operationally configured models significantly, artificially reduced MW amplitudes in the upper stratosphere and mesosphere. In the IFS, parameterized GW drags partly compensated this deficiency, but still, total drags were ≈6 times smaller than that resolved at Δx ≈ 3 km. Meridionally propagating MWs significantly enhance zonal drag over the Drake Passage. Interestingly, drag associated with meridional fluxes of zonal momentum (i.e.,
Significance Statement
This study had three purposes: to quantitatively evaluate how well four state-of-the-science weather models could reproduce observed mountain waves (MWs) in the middle atmosphere, to compare the simulated MWs within the models, and to quantitatively evaluate two MW parameterizations in a widely used climate model. These models reproduced observed MWs with remarkable skill. Still, MW parameterizations are necessary in current Δx ≈ 10-km resolution global weather models. Even Δx ≈ 3-km resolution does not appear to be high enough to represent all momentum-fluxing MW scales. Meridionally propagating MWs can significantly influence zonal winds over the Drake Passage. Parameterizations that handle horizontal propagation may need to consider horizontal fluxes of horizontal momentum in order to get the direction of their forcing correct.
Abstract
Stratospheric gravity waves observed during the DEEPWAVE research flight RF25 over the Southern Ocean are analyzed and compared with numerical weather prediction (NWP) model results. The quantitative agreement of the NWP model output and the tropospheric and lower-stratospheric observations is remarkable. The high-resolution NWP models are even able to reproduce qualitatively the observed upper-stratospheric gravity waves detected by an airborne Rayleigh lidar. The usage of high-resolution ERA5 data—partially capturing the long internal gravity waves—enabled a thorough interpretation of the particular event. Here, the observed and modeled gravity waves are excited by the stratospheric flow past a deep tropopause depression belonging to an eastward-propagating Rossby wave train. In the reference frame of the propagating Rossby wave, vertically propagating hydrostatic gravity waves appear stationary; in reality, of course, they are transient and propagate horizontally at the phase speed of the Rossby wave. The subsequent refraction of these transient gravity waves into the polar night jet explains their observed and modeled patchy stratospheric occurrence near 60°S. The combination of both unique airborne observations and high-resolution NWP output provides evidence for the one case investigated in this paper. As the excitation of such gravity waves persists during the quasi-linear propagation phase of the Rossby wave’s life cycle, a hypothesis is formulated that parts of the stratospheric gravity wave belt over the Southern Ocean might be generated by such Rossby wave trains propagating along the midlatitude waveguide.
Abstract
Stratospheric gravity waves observed during the DEEPWAVE research flight RF25 over the Southern Ocean are analyzed and compared with numerical weather prediction (NWP) model results. The quantitative agreement of the NWP model output and the tropospheric and lower-stratospheric observations is remarkable. The high-resolution NWP models are even able to reproduce qualitatively the observed upper-stratospheric gravity waves detected by an airborne Rayleigh lidar. The usage of high-resolution ERA5 data—partially capturing the long internal gravity waves—enabled a thorough interpretation of the particular event. Here, the observed and modeled gravity waves are excited by the stratospheric flow past a deep tropopause depression belonging to an eastward-propagating Rossby wave train. In the reference frame of the propagating Rossby wave, vertically propagating hydrostatic gravity waves appear stationary; in reality, of course, they are transient and propagate horizontally at the phase speed of the Rossby wave. The subsequent refraction of these transient gravity waves into the polar night jet explains their observed and modeled patchy stratospheric occurrence near 60°S. The combination of both unique airborne observations and high-resolution NWP output provides evidence for the one case investigated in this paper. As the excitation of such gravity waves persists during the quasi-linear propagation phase of the Rossby wave’s life cycle, a hypothesis is formulated that parts of the stratospheric gravity wave belt over the Southern Ocean might be generated by such Rossby wave trains propagating along the midlatitude waveguide.
Abstract
The dynamical and thermodynamical features of Amazonian 2-day westward-propagating inertia–gravity (WIG) waves are examined. On the basis of a linear regression analysis of satellite brightness temperature and data from the 2014–15 Observations and Modeling of the Green Ocean Amazon (GoAmazon) field campaign, it is shown that Amazonian WIG waves exhibit structure and propagation characteristics consistent with the n = 1 WIG waves from shallow water theory. These WIG waves exhibit a pronounced seasonality, with peak activity occurring from March to May and a minimum occurring from June to September. Evidence is shown that mesoscale convective systems over the Amazon are frequently organized in 2-day WIG waves. Results suggest that many of the Amazonian WIG waves come from preexisting 2-day waves over the Atlantic, which slow down when coupled with the deeper, more intense convection over tropical South America. In contrast to WIG waves that occur over the ocean, Amazonian 2-day WIG waves exhibit a pronounced signature in surface temperature, moisture, and heat fluxes.
Abstract
The dynamical and thermodynamical features of Amazonian 2-day westward-propagating inertia–gravity (WIG) waves are examined. On the basis of a linear regression analysis of satellite brightness temperature and data from the 2014–15 Observations and Modeling of the Green Ocean Amazon (GoAmazon) field campaign, it is shown that Amazonian WIG waves exhibit structure and propagation characteristics consistent with the n = 1 WIG waves from shallow water theory. These WIG waves exhibit a pronounced seasonality, with peak activity occurring from March to May and a minimum occurring from June to September. Evidence is shown that mesoscale convective systems over the Amazon are frequently organized in 2-day WIG waves. Results suggest that many of the Amazonian WIG waves come from preexisting 2-day waves over the Atlantic, which slow down when coupled with the deeper, more intense convection over tropical South America. In contrast to WIG waves that occur over the ocean, Amazonian 2-day WIG waves exhibit a pronounced signature in surface temperature, moisture, and heat fluxes.
Abstract
Planetary waves disturbed the hitherto stable Arctic stratospheric polar vortex in the middle of January 2016 in such a way that unique tropospheric and stratospheric flow conditions for vertically and horizontally propagating mountain waves developed. Coexisting strong low-level westerly winds across almost all European mountain ranges plus the almost zonally aligned polar-front jet created these favorable conditions for deeply propagating gravity waves. Furthermore, the northward displacement of the polar night jet resulted in a widespread coverage of stratospheric mountain waves trailing across Northern Europe. This paper describes the particular meteorological setting by analyzing the tropospheric and stratospheric flows based on the ERA5 data. The potential of the flow for exciting internal gravity waves from nonorographic sources is evaluated across all altitudes by considering various indices to indicate flow imbalances as δ, Ro, Ro ζ , Ro⊥, and ΔNBE. The analyzed gravity waves are described and characterized. The main finding of this case study is the exceptionally vast extension of the mountain waves trailing to high latitudes originating from the flow across the mountainous sources that are located at about 45°N. The magnitudes of the simulated stratospheric temperature perturbations attain values larger than 10 K and are comparable to values as documented by recent case studies of large-amplitude mountain waves over South America. The zonal means of the resolved and parameterized stratospheric wave drag during the mountain wave event peak at −4.5 and −32.2 m s−1 day−1, respectively.
Abstract
Planetary waves disturbed the hitherto stable Arctic stratospheric polar vortex in the middle of January 2016 in such a way that unique tropospheric and stratospheric flow conditions for vertically and horizontally propagating mountain waves developed. Coexisting strong low-level westerly winds across almost all European mountain ranges plus the almost zonally aligned polar-front jet created these favorable conditions for deeply propagating gravity waves. Furthermore, the northward displacement of the polar night jet resulted in a widespread coverage of stratospheric mountain waves trailing across Northern Europe. This paper describes the particular meteorological setting by analyzing the tropospheric and stratospheric flows based on the ERA5 data. The potential of the flow for exciting internal gravity waves from nonorographic sources is evaluated across all altitudes by considering various indices to indicate flow imbalances as δ, Ro, Ro ζ , Ro⊥, and ΔNBE. The analyzed gravity waves are described and characterized. The main finding of this case study is the exceptionally vast extension of the mountain waves trailing to high latitudes originating from the flow across the mountainous sources that are located at about 45°N. The magnitudes of the simulated stratospheric temperature perturbations attain values larger than 10 K and are comparable to values as documented by recent case studies of large-amplitude mountain waves over South America. The zonal means of the resolved and parameterized stratospheric wave drag during the mountain wave event peak at −4.5 and −32.2 m s−1 day−1, respectively.
Abstract
Horizontally dispersing gravity waves with horizontal wavelengths of 30–40 km were observed at mesospheric altitudes over Auckland Island by the airborne advanced mesospheric temperature mapper during a Deep Propagating Gravity Wave Experiment (DEEPWAVE) research flight on 14 July 2014. A 3D nonlinear compressible model is used to determine which propagation conditions enabled gravity wave penetration into the mesosphere and how the resulting instability characteristics led to widespread momentum deposition. Results indicate that linear tunneling through the polar night jet enabled quick gravity wave propagation from the surface up to the mesopause, while subsequent instability processes reveal large rolls that formed in the negative shear above the jet maximum and led to significant momentum deposition as they descended. This study suggests that gravity wave tunneling is a viable source for this case and other deep propagation events reaching the mesosphere and lower thermosphere.
Abstract
Horizontally dispersing gravity waves with horizontal wavelengths of 30–40 km were observed at mesospheric altitudes over Auckland Island by the airborne advanced mesospheric temperature mapper during a Deep Propagating Gravity Wave Experiment (DEEPWAVE) research flight on 14 July 2014. A 3D nonlinear compressible model is used to determine which propagation conditions enabled gravity wave penetration into the mesosphere and how the resulting instability characteristics led to widespread momentum deposition. Results indicate that linear tunneling through the polar night jet enabled quick gravity wave propagation from the surface up to the mesopause, while subsequent instability processes reveal large rolls that formed in the negative shear above the jet maximum and led to significant momentum deposition as they descended. This study suggests that gravity wave tunneling is a viable source for this case and other deep propagation events reaching the mesosphere and lower thermosphere.
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.
Abstract
Stationary gravity waves, such as mountain lee waves, are effectively described by Grimshaw’s dissipative modulation equations even in high altitudes where they become nonlinear due to their large amplitudes. In this theoretical study, a wave-Reynolds number is introduced to characterize general solutions to these modulation equations. This nondimensional number relates the vertical linear group velocity with wavenumber, pressure scale height, and kinematic molecular/eddy viscosity. It is demonstrated by analytic and numerical methods that Lindzen-type waves in the saturation region, that is, where the wave-Reynolds number is of order unity, destabilize by transient perturbations. It is proposed that this mechanism may be a generator for secondary waves due to direct wave–mean-flow interaction. By assumption, the primary waves are exactly such that altitudinal amplitude growth and viscous damping are balanced and by that the amplitude is maximized. Implications of these results on the relation between mean-flow acceleration and wave breaking heights are discussed.
Abstract
Stationary gravity waves, such as mountain lee waves, are effectively described by Grimshaw’s dissipative modulation equations even in high altitudes where they become nonlinear due to their large amplitudes. In this theoretical study, a wave-Reynolds number is introduced to characterize general solutions to these modulation equations. This nondimensional number relates the vertical linear group velocity with wavenumber, pressure scale height, and kinematic molecular/eddy viscosity. It is demonstrated by analytic and numerical methods that Lindzen-type waves in the saturation region, that is, where the wave-Reynolds number is of order unity, destabilize by transient perturbations. It is proposed that this mechanism may be a generator for secondary waves due to direct wave–mean-flow interaction. By assumption, the primary waves are exactly such that altitudinal amplitude growth and viscous damping are balanced and by that the amplitude is maximized. Implications of these results on the relation between mean-flow acceleration and wave breaking heights are discussed.
Abstract
This paper compares two different approaches for the efficient modeling of subgrid-scale inertia–gravity waves in a rotating compressible atmosphere. The first approach, denoted as the pseudomomentum scheme, exploits the fact that in a Lagrangian-mean reference frame the response of a large-scale flow can only be due to forcing momentum. Present-day gravity wave parameterizations follow this route. They do so, however, in an Eulerian-mean formulation. Transformation to that reference frame leads, under certain assumptions, to pseudomomentum-flux convergence by which the momentum is to be forced. It can be shown that this approach is justified if the large-scale flow is in geostrophic and hydrostatic balance. Otherwise, elastic and thermal effects might be lost. In the second approach, called the direct scheme and not relying on such assumptions, the large-scale flow is forced both in the momentum equation, by anelastic momentum-flux convergence and an additional elastic term, and in the entropy equation, via entropy-flux convergence. A budget analysis based on one-dimensional wave packets suggests that the comparison between the abovementioned two schemes should be sensitive to the following two parameters: 1) the intrinsic frequency and 2) the wave packet scale. The smaller the intrinsic frequency is, the greater their differences are. More importantly, with high-resolution wave-resolving simulations as a reference, this study shows conclusive evidence that the direct scheme is more reliable than the pseudomomentum scheme, regardless of whether one-dimensional or two-dimensional wave packets are considered. In addition, sensitivity experiments are performed to further investigate the relative importance of each term in the direct scheme, as well as the wave–mean flow interactions during the wave propagation.
Abstract
This paper compares two different approaches for the efficient modeling of subgrid-scale inertia–gravity waves in a rotating compressible atmosphere. The first approach, denoted as the pseudomomentum scheme, exploits the fact that in a Lagrangian-mean reference frame the response of a large-scale flow can only be due to forcing momentum. Present-day gravity wave parameterizations follow this route. They do so, however, in an Eulerian-mean formulation. Transformation to that reference frame leads, under certain assumptions, to pseudomomentum-flux convergence by which the momentum is to be forced. It can be shown that this approach is justified if the large-scale flow is in geostrophic and hydrostatic balance. Otherwise, elastic and thermal effects might be lost. In the second approach, called the direct scheme and not relying on such assumptions, the large-scale flow is forced both in the momentum equation, by anelastic momentum-flux convergence and an additional elastic term, and in the entropy equation, via entropy-flux convergence. A budget analysis based on one-dimensional wave packets suggests that the comparison between the abovementioned two schemes should be sensitive to the following two parameters: 1) the intrinsic frequency and 2) the wave packet scale. The smaller the intrinsic frequency is, the greater their differences are. More importantly, with high-resolution wave-resolving simulations as a reference, this study shows conclusive evidence that the direct scheme is more reliable than the pseudomomentum scheme, regardless of whether one-dimensional or two-dimensional wave packets are considered. In addition, sensitivity experiments are performed to further investigate the relative importance of each term in the direct scheme, as well as the wave–mean flow interactions during the wave propagation.
