Browse
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
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.
Abstract
Large uncertainties remain with respect to the representation of atmospheric gravity waves (GWs) in general circulation models (GCMs) with coarse grids. Insufficient parameterizations result from a lack of observational constraints on the parameters used in GW parameterizations as well as from physical inconsistencies between parameterizations and reality. For instance, parameterizations make oversimplifying assumptions about the generation and propagation of GWs. Increasing computational capabilities now allow GCMs to run at grid spacings that are sufficiently fine to resolve a major fraction of the GW spectrum. This study presents the first intercomparison of resolved GW pseudomomentum fluxes (GWMFs) in global convection-permitting simulations and those derived from satellite observations. Six simulations of three different GCMs are analyzed over the period of one month of August to assess the sensitivity of GWMF to model formulation and horizontal grid spacing. The simulations reproduce detailed observed features of the global GWMF distribution, which can be attributed to realistic GWs from convection, orography, and storm tracks. Yet the GWMF magnitudes differ substantially between simulations. Differences in the strength of convection may help explain differences in the GWMF between simulations of the same model in the summer low latitudes where convection is the primary source. Across models, there is no evidence for a systematic change with resolution. Instead, GWMF is strongly affected by model formulation. The results imply that validating the realism of simulated GWs across the entire resolved spectrum will remain a difficult challenge not least because of a lack of appropriate observational data.
Abstract
Large uncertainties remain with respect to the representation of atmospheric gravity waves (GWs) in general circulation models (GCMs) with coarse grids. Insufficient parameterizations result from a lack of observational constraints on the parameters used in GW parameterizations as well as from physical inconsistencies between parameterizations and reality. For instance, parameterizations make oversimplifying assumptions about the generation and propagation of GWs. Increasing computational capabilities now allow GCMs to run at grid spacings that are sufficiently fine to resolve a major fraction of the GW spectrum. This study presents the first intercomparison of resolved GW pseudomomentum fluxes (GWMFs) in global convection-permitting simulations and those derived from satellite observations. Six simulations of three different GCMs are analyzed over the period of one month of August to assess the sensitivity of GWMF to model formulation and horizontal grid spacing. The simulations reproduce detailed observed features of the global GWMF distribution, which can be attributed to realistic GWs from convection, orography, and storm tracks. Yet the GWMF magnitudes differ substantially between simulations. Differences in the strength of convection may help explain differences in the GWMF between simulations of the same model in the summer low latitudes where convection is the primary source. Across models, there is no evidence for a systematic change with resolution. Instead, GWMF is strongly affected by model formulation. The results imply that validating the realism of simulated GWs across the entire resolved spectrum will remain a difficult challenge not least because of a lack of appropriate observational data.
Abstract
As present weather forecast codes and increasingly many atmospheric climate models resolve at least part of the mesoscale flow, and hence also internal gravity waves (GWs), it is natural to ask whether even in such configurations subgrid-scale GWs might impact the resolved flow and how their effect could be taken into account. This motivates a theoretical and numerical investigation of the interactions between unresolved submesoscale and resolved mesoscale GWs, using Boussinesq dynamics for simplicity. By scaling arguments, first a subset of submesoscale GWs that can indeed influence the dynamics of mesoscale GWs is identified. Therein, hydrostatic GWs with wavelengths corresponding to the largest unresolved scales of present-day limited-area weather forecast models are an interesting example. A large-amplitude WKB theory, allowing for a mesoscale unbalanced flow, is then formulated, based on multiscale asymptotic analysis utilizing a proper scale-separation parameter. Purely vertical propagation of submesoscale GWs is found to be most important, implying inter alia that the resolved flow is only affected by the vertical flux convergence of submesoscale horizontal momentum at leading order. In turn, submesoscale GWs are refracted by mesoscale vertical wind shear while conserving their wave-action density. An efficient numerical implementation of the theory uses a phase-space ray tracer, thus handling the frequent appearance of caustics. The WKB approach and its numerical implementation are validated successfully against submesoscale-resolving simulations of the resonant radiation of mesoscale inertia GWs by a horizontally as well as vertically confined submesoscale GW packet.
Abstract
As present weather forecast codes and increasingly many atmospheric climate models resolve at least part of the mesoscale flow, and hence also internal gravity waves (GWs), it is natural to ask whether even in such configurations subgrid-scale GWs might impact the resolved flow and how their effect could be taken into account. This motivates a theoretical and numerical investigation of the interactions between unresolved submesoscale and resolved mesoscale GWs, using Boussinesq dynamics for simplicity. By scaling arguments, first a subset of submesoscale GWs that can indeed influence the dynamics of mesoscale GWs is identified. Therein, hydrostatic GWs with wavelengths corresponding to the largest unresolved scales of present-day limited-area weather forecast models are an interesting example. A large-amplitude WKB theory, allowing for a mesoscale unbalanced flow, is then formulated, based on multiscale asymptotic analysis utilizing a proper scale-separation parameter. Purely vertical propagation of submesoscale GWs is found to be most important, implying inter alia that the resolved flow is only affected by the vertical flux convergence of submesoscale horizontal momentum at leading order. In turn, submesoscale GWs are refracted by mesoscale vertical wind shear while conserving their wave-action density. An efficient numerical implementation of the theory uses a phase-space ray tracer, thus handling the frequent appearance of caustics. The WKB approach and its numerical implementation are validated successfully against submesoscale-resolving simulations of the resonant radiation of mesoscale inertia GWs by a horizontally as well as vertically confined submesoscale GW packet.
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.
