Gravity Wave Momentum Flux Research Articles

Oblique propagation of monsoon gravity waves during the northern hemisphere 2007 summer

AbstractWe present a combination of satellite observation and high‐resolution model output to understand monsoon convection as a source of high‐latitude mesospheric gravity waves (GWs). The GWs generated over the Northern Hemisphere (NH) monsoon region during the 2007 summer and the role of the winds in focusing these GWs toward the high‐latitude middle atmosphere are analyzed using the Sounding of the Atmosphere using Broadband Emission Radiometry/Thermosphere Ionosphere Mesosphere Energetics and Dynamics (SABER/TIMED) satellite temperature data and the high‐resolution Navy Operational Global Atmospheric Prediction System‐Advanced Level Physics High Altitude (NOGAPS/ALPHA) model results. In the NH, above the stratosphere, the monsoon GW Momentum Flux (GWMF) exhibits a poleward tilt that follows the slanted structure of the easterly jet. The correlation coefficients (>0.5) between the time series of NH tropical stratospheric GWMF and the global winds also have a slanted structure that coincide with the easterly jet, confirming the modeling theory that stratospheric monsoon GWs are refracted into the summer easterly jet and can reach the high‐latitude mesosphere. Since Polar Mesospheric Clouds (PMCs) are sensitive indicators of changes in the polar summer mesosphere, we compared the time series of tropical stratospheric GWMF to the PMC occurrence frequency (OF) obtained from the Cloud Imaging and Particle Size/Aeronomy of Ice in the Mesosphere satellite data to assess the influence of this wave focusing in the mesosphere. There is a significant positive correlation between the high‐latitude PMC OF and the tropical stratospheric GWMF suggesting a definite influence of monsoon GWs on the high‐latitude mesosphere. The disagreement in correlation at the end of the PMC season is attributed to the enhancement of the quasi 5 day planetary wave dominating over the influence of monsoon GWs on PMCs.

On the Relation between Gravity Waves and Wind Speed in the Lower Stratosphere over the Southern Ocean

Abstract The relationship between gravity wave momentum fluxes and local wind speed is investigated for oceanic regions at high southern latitudes during austral spring. The motivation is to better describe the gravity wave field by identifying a simple relationship between gravity waves and the large-scale flow. The tools used to describe the gravity waves are probability density functions of the gravity wave momentum fluxes. Three independent datasets covering high latitudes in the Southern Hemisphere springtime are analyzed: simulations with a mesoscale model, analyses from the European Centre for Medium-Range Weather Forecasts, and observations from superpressure balloons of the Concordiasi campaign in 2010. A remarkably robust relation is found, with stronger momentum fluxes much more likely in regions of strong winds. The tails of the probability density functions are well described as lognormal. The median momentum flux increases linearly with background wind speed: for winds stronger than 50 m s−1, the median gravity wave momentum fluxes are about 4 times larger than for winds weaker than 10 m s−1. From model output, this relation is found to be relevant from the tropopause to the midstratosphere at least. The flux dependence on wind speed shows a somewhat steeper slope at higher altitude. Several different processes contribute to this relation, involving both the distribution of sources and the effects of propagation and filtering. It is argued that the location of tropospheric sources is the main contributor in the upper troposphere and lowermost stratosphere and that lateral propagation into regions of strong winds becomes increasingly important above.

Unexpected climatological behavior of MLT gravity wave momentum flux in the lee of the Southern Andes hot spot

AbstractThe Southern Argentina Agile MEteor Radar (SAAMER), located at Tierra del Fuego (53.7°S, 67.7°W), has been providing near‐continuous high‐resolution measurements of winds and high‐frequency gravity wave (GW) momentum fluxes of the mesopause region since May 2008. As SAAMER is located in the lee of the largest seasonal GW hot spot on Earth, this is a key location to study GWs and their interaction with large‐scale motions. GW momentum flux climatologies are shown for the first time for this location and discussed in light of these unique dynamics. Particularly, the large eastward GW momentum fluxes during local winter are surprising, as these observations cannot be explained by the direct upward propagation of expected large‐amplitude mountain waves (MWs) through the eastward stratospheric jet. Instead, these results are interpreted as secondary GWs propagating away from stratospheric sources over the Andes accompanying MW breaking over the Southern Andes.

Directional gravity wave momentum fluxes in the stratosphere derived from high‐resolution AIRS temperature data

AbstractIn order to reduce uncertainties in modeling the stratospheric circulation, global observations of gravity wave momentum flux (GWMF) vectors are required for comparison with distributions of resolved and parametrized GWMF in global models. For the first time, we derive GWMF vectors globally from data of a nadir‐viewing satellite instrument: we apply a 3‐D method to an Atmospheric Infrared Sounder (AIRS) temperature data set that was optimized for gravity wave (GW) analysis. For January 2009, the resulting distributions of GW amplitudes and of net GWMF highlight the importance of GWs in the polar vortex and the summertime subtropics. Net GWMF is preferentially directed opposite to the background wind, and, interestingly, it is dominated by large‐amplitude GWs of relatively long horizontal wavelength. For convective GW sources, these large horizontal scales are in contradiction with traditional thoughts. However, the observational filter effect needs to be kept in mind when interpreting the results.

Gravity wave momentum flux in the mesosphere measured by VHF radar at Davis, Antarctica

AbstractPolar mesospheric summer echoes observed in VHF MST radar measurements at Davis, Antarctica 68.5°S, 78.0°E, are used to measure high‐frequency gravity wave momentum flux in the vicinity of the mesopause during the Southern Hemisphere summer of 2014–2015. The momentum flux spectrum is generally isotropic with a slight bias in the seasonal mean toward the southeast. Zonal and meridional components of the seasonal mean momentum flux are an order of magnitude smaller than the absolute momentum flux. Diurnal and semidiurnal modulation of the momentum flux correspond to perturbations of the seasonal mean winds by atmospheric tides. Quantitative measures of gravity wave intermittency yield high values and increase with altitude through the mesopause region. These intermittency characteristics are attributed to the increase with altitude of wind variability through the mesopause region as a result of the tides.

Satellite observations of middle atmosphere gravity wave absolute momentum flux and of its vertical gradient during recent stratospheric warmings

Abstract. Sudden stratospheric warmings (SSWs) are circulation anomalies in the polar region during winter. They mostly occur in the Northern Hemisphere and affect also surface weather and climate. Both planetary waves and gravity waves contribute to the onset and evolution of SSWs. While the role of planetary waves for SSW evolution has been recognized, the effect of gravity waves is still not fully understood, and has not been comprehensively analyzed based on global observations. In particular, information on the gravity wave driving of the background winds during SSWs is still missing.We investigate the boreal winters from 2001/2002 until 2013/2014. Absolute gravity wave momentum fluxes and gravity wave dissipation (potential drag) are estimated from temperature observations of the satellite instruments HIRDLS and SABER. In agreement with previous work, we find that sometimes gravity wave activity is enhanced before or around the central date of major SSWs, particularly during vortex-split events. Often, SSWs are associated with polar-night jet oscillation (PJO) events. For these events, we find that gravity wave activity is strongly suppressed when the wind has reversed from eastward to westward (usually after the central date of a major SSW). In addition, gravity wave potential drag at the bottom of the newly forming eastward-directed jet is remarkably weak, while considerable potential drag at the top of the jet likely contributes to the downward propagation of both the jet and the new elevated stratopause. During PJO events, we also find some indication for poleward propagation of gravity waves. Another striking finding is that obviously localized gravity wave sources, likely mountain waves and jet-generated gravity waves, play an important role during the evolution of SSWs and potentially contribute to the triggering of SSWs by preconditioning the shape of the polar vortex. The distribution of these hot spots is highly variable and strongly depends on the zonal and meridional shape of the background wind field, indicating that a pure zonal average view sometimes is a too strong simplification for the strongly perturbed conditions during the evolution of SSWs.

The QBO, gravity waves forced by tropical convection, and ENSO

AbstractBy means of theory, a simplified cartoon illustrating wave forcing of the stratospheric quasi‐biennial oscillation (QBO), and general circulation modeling of the QBO, it is argued that the period of the QBO is mainly controlled by the magnitude of the gravity wave (GW) vertical fluxes of horizontal momentum (GWMF) forcing the QBO, while the QBO amplitude is mainly determined by the phase speeds of the GWs that make up this momentum flux. It is furthermore argued that it is the zonally averaged GWMF that principally determines the QBO period irrespective of the longitudinal distribution of this GW momentum flux. These concepts are used to develop a hypothesis for the cause of a previously reported El Niño–Southern Oscillation (ENSO) modulation of QBO periods and amplitudes. Some observational evidence is reported for the ENSO modulation of QBO amplitudes to have been different before the 1980s than after about 1990. A hypothesis is also given to explain this in terms of the different ENSO modulation of tropical deep convection that took place before the 1980s from that which occurred after about 1990. The observational evidence, while consistent with our hypotheses, does not prove that our hypotheses are correct given the small number of El Niños and La Niñas that occurred in the early and later periods. Further research is needed to support or refute our hypotheses.

An Analysis of Gravity Wave Spectral Characteristics in Moist Baroclinic Jet–Front Systems

Abstract This study investigates gravity wave spectral characteristics based on high-resolution mesoscale simulations of idealized moist baroclinic jet–front systems with varying degrees of convective instability, with the intent of improving nonorographic gravity wave parameterizations. In all experiments, there is a clear dominance of negative vertical flux of zonal momentum. The westward momentum flux is distributed around the estimated ground-based baroclinic wave phase velocity along the zonal direction, while strong moist runs indicate a dipole structure pattern with stronger westward momentum flux centers at slower phase velocity and weaker eastward momentum flux centers at faster phase velocity. The spectral properties of short-scale wave components (50–200 km) generally differ from those of medium-scale ones (200–600 km). Compared to the medium-scale wave components, the momentum flux phase speed spectra for the short-scale ones appear to be more sensitive to the increasing initial moisture content. The spectral behavior in horizontal wavenumber space or phase velocity space is highly anisotropic, with a noticeable preference along the jet direction, except for the short-scale components in strong moist runs. It is confirmed that the dry gravity wave source (i.e., upper jet and/or surface front) generates a relatively narrow and less symmetrical power spectrum (dominated by negative momentum flux) centered around lower phase velocity and horizontal wavenumber, whereas the moist gravity wave source (i.e., moist convection) generates a broader and more symmetrical power spectrum, with a broader range of phase speeds and horizontal wavenumbers. This study also shows that the properties of gravity wave momentum flux depend on the location relative to the baroclinic jet.

Modeling the QBO-Improvements resulting from higher-model vertical resolution.

Using the NASA Goddard Institute for Space Studies (GISS) climate model, it is shown that with proper choice of the gravity wave momentum flux entering the stratosphere and relatively fine vertical layering of at least 500 m in the upper troposphere‐lower stratosphere (UTLS), a realistic stratospheric quasi‐biennial oscillation (QBO) is modeled with the proper period, amplitude, and structure down to tropopause levels. It is furthermore shown that the specified gravity wave momentum flux controls the QBO period whereas the width of the gravity wave momentum flux phase speed spectrum controls the QBO amplitude. Fine vertical layering is required for the proper downward extension to tropopause levels as this permits wave‐mean flow interactions in the UTLS region to be resolved in the model. When vertical resolution is increased from 1000 to 500 m, the modeled QBO modulation of the tropical tropopause temperatures increasingly approach that from observations, and the “tape recorder” of stratospheric water vapor also approaches the observed. The transport characteristics of our GISS models are assessed using age‐of‐air and N2O diagnostics, and it is shown that some of the deficiencies in model transport that have been noted in previous GISS models are greatly improved for all of our tested model vertical resolutions. More realistic tropical‐extratropical transport isolation, commonly referred to as the “tropical pipe,” results from the finer vertical model layering required to generate a realistic QBO.

Spectral distribution of gravity wave momentum fluxes over the Antarctic Peninsula from Concordiasi superpressure balloon data

AbstractGravity waves generated by flow over the steep topography of the Antarctic Peninsula transport significant amounts of zonal and meridional momentum into the stratosphere. Quantitative determination of this transport has been carried out for wave periods of 1 h or greater using data from a previous Antarctic superpressure balloon campaign in austral spring 2005 (VORCORE). The present study uses data from the later Concordiasi campaign (2010) to extend the momentum flux determination to shorter periods. Maps of the vertical fluxes of meridional and zonal momentum are presented for periods down to 12 min. We find that the momentum fluxes for periods below 1 h are comparable to those at longer periods, despite larger variances at longer periods. The momentum fluxes in the vicinity of the peninsula provide a significant zonal acceleration of the lower stratosphere, confirming a conclusion from the VORCORE data. The geographical distribution of fluxes around the peninsula has peaks both leeward and windward of the main terrain features. Numerical simulations suggest that the separate peaks may be related to wave transience caused by unsteady winds over the peninsula. Momentum fluxes comprise a main distribution maximizing at moderate flux values and a secondary distribution maximizing at high values exhibiting a high degree of intermittency. The high flux events account for the largest part of the average flux and suggest that drag parameterizations should take them into account. It is found that waves generated by the jet stream are also a significant source of momentum flux.

Characteristics of Gravity Waves from Convection and Implications for Their Parameterization in Global Circulation Models

Abstract Characteristic properties of gravity waves from convection over the continental United States are derived from idealized high-resolution numerical simulations. In a unique modeling approach, waves are forced by a realistic thermodynamic source based on observed precipitation data. The square of the precipitation rate and the gravity wave momentum fluxes both show lognormal occurrence distributions, with long tails of extreme events. Convectively generated waves can give forces in the lower stratosphere that at times rival orographic wave forcing. Throughout the stratosphere, zonal forces due to convective wave drag are much stronger than accounted for by current gravity wave drag parameterizations, so their contribution to the summer branch of the stratospheric Brewer–Dobson circulation is in fact much larger than models predict. A comparison of these forces to previous estimates of the total drag implies that convectively generated gravity waves are a primary source of summer-hemisphere stratospheric wave drag. Furthermore, intermittency and strength of the zonal forces due to convective gravity wave drag in the lower stratosphere resemble analysis increments, suggesting that a more realistic representation of these waves may help alleviate model biases on synoptic scales. The properties of radar precipitation and gravity waves seen in this study lead to a proposed change for future parameterization methods that would give more realistic drag forces in the stratosphere without compromising mesospheric gravity wave drag.

Impact of cyclone Nilam on tropical lower atmospheric dynamics

A deep depression formed over the Bay of Bengal on 28 October 2012, and developed into a cyclonic storm. After landfall near the south coast of Chennai, cyclone Nilam moved north-northwestwards. Coordinated experiments were conducted from the Indian stations of Gadanki (13.5°N, 79.2°E) and Hyderabad (17.4°N, 78.5°E) to study the modification of gravity-wave activity and turbulence by cyclone Nilam, using GPS radiosonde and mesosphere–stratosphere–troposphere radar data. The horizontal velocities underwent large changes during the closest approach of the storm to the experimental sites. Hodograph analysis revealed that inertia gravity waves (IGWs) associated with the cyclone changed their directions from northeast (control time) to northwest following the path of the cyclone. The momentum flux of IGWs and short-period gravity waves (1–8 h) enhanced prior to, and during, the passage of the storm (±0.05 m2 s−2 and ±0.3 m2 s−2, respectively), compared to the flux after its passage. The corresponding body forces underwent similar changes, with values ranging between ±2–4 m s−1 d−1 and ±12–15 m s−1 d−1. The turbulence refractivity structure constant (C 2 ) showed large values below 10 km before the passage of the cyclone when humidity in the region was very high. Turbulence and humidity reduced during the passage of the storm when a turbulent layer at ~17 km became more intense. Turbulence in the lower troposphere and near the tropopause became weak after the passage of the cyclone.

Tuning of a convective gravity wave source scheme based on HIRDLS observations

Abstract. Convection as one dominant source of atmospheric gravity waves (GWs) has been the focus of investigation over recent years. However, its spatial and temporal forcing scales are not well known. In this work we address this open issue by a systematic verification of free parameters of the Yonsei convective GW source scheme based on observations from the High Resolution Dynamics Limb Sounder (HIRDLS). The instrument can only see a limited portion of the gravity wave spectrum due to visibility effects and observation geometry. To allow for a meaningful comparison of simulated GWs to observations, a comprehensive filter, which mimics the instrument limitations, is applied to the simulated waves. By this approach, only long horizontal-scale convective GWs are addressed. Results show that spectrum, distribution of momentum flux, and zonal mean forcing of long horizontal-scale convective GWs can be successfully simulated by the superposition of three or four combinations of parameter sets reproducing the observed GW spectrum. These selected parameter sets are different for northern and southern summer. Although long horizontal-scale waves are only part of the full spectrum of convective GWs, the momentum flux of these waves is found to be significant and relevant for the driving of the QBO (quasi-biennial oscillation). The zonal momentum balance is considered in vertical cross sections of GW momentum flux (GWMF) and GW drag (GWD). Global maps of the horizontal distribution of GWMF are considered and consistency between simulated results and HIRDLS observations is found. The latitude dependence of the zonal phase speed spectrum of GWMF and its change with altitude is discussed.

A two-dimensional Stockwell transform for gravity wave analysis of AIRS measurements

Abstract. Gravity waves (GWs) play a crucial role in the dynamics of the earth's atmosphere. These waves couple lower, middle and upper atmospheric layers by transporting and depositing energy and momentum from their sources to great heights. The accurate parameterisation of GW momentum flux is of key importance to general circulation models but requires accurate measurement of GW properties, which has proved challenging. For more than a decade, the nadir-viewing Atmospheric Infrared Sounder (AIRS) aboard NASA's Aqua satellite has made global, two-dimensional (2-D) measurements of stratospheric radiances in which GWs can be detected. However, one problem with current one-dimensional methods for GW analysis of these data is that they can introduce significant unwanted biases. Here, we present a new analysis method that resolves this problem. Our method uses a 2-D Stockwell transform (2DST) to measure GW amplitudes, horizontal wavelengths and directions of propagation using both the along-track and cross-track dimensions simultaneously. We first test our new method and demonstrate that it can accurately measure GW properties in a specified wave field. We then show that by using a new elliptical spectral window in the 2DST, in place of the traditional Gaussian, we can dramatically improve the recovery of wave amplitude over the standard approach. We then use our improved method to measure GW properties and momentum fluxes in AIRS measurements over two regions known to be intense hotspots of GW activity: (i) the Drake Passage/Antarctic Peninsula and (ii) the isolated mountainous island of South Georgia. The significance of our new 2DST method is that it provides more accurate, unbiased and better localised measurements of key GW properties compared to most current methods. The added flexibility offered by the scaling parameter and our new spectral window presented here extend the usefulness of our 2DST method to other areas of geophysical data analysis and beyond.

QBO modulation of the mesopause gravity wave momentum flux over Tierra del Fuego

AbstractThe interannual variability of the mesosphere and lower thermosphere (MLT) gravity wave momentum flux over southern midlatitudes (53.7°S) has been studied using more than 7 years of meteor radar observations at Río Grande, Argentina. A modulation, with periods similar to that of the equatorial stratospheric quasi‐biennial oscillation (QBO), is observed in the vertical flux of zonal as well as meridional momentum. The QBO signal is largest in the zonal component during summer and is in phase with the stratospheric QBO at 50 hPa (∼21 km). The relation between the stratospheric QBO and the QBO modulation in the MLT gravity wave forcing (derived from the divergence of the momentum flux) was found to be consistent with that expected from the Holton‐Tan effect coupled to the interhemispheric coupling mechanism. These results provide the first observational support for the existence of the midlatitude gravity wave forcing anomalies as hypothesized in the interhemispheric coupling mechanism.

Gravity-wave momentum fluxes in the mesosphere over Ascension Island (8° S, 14° W) and the anomalous zonal winds of the semi-annual oscillation in 2002

Abstract. Anomalously strong westward winds during the first phase of the equatorial mesospheric semi-annual oscillation (MSAO) have been attributed to unusual filtering conditions producing exceptional gravity-wave fluxes. We test this hypothesis using meteor-radar measurements made over Ascension Island (8° S, 14° W). An anomalous wind event in 2002 of −85.5 ms−1 occurred simultaneously with the momentum fluxes of high-frequency gravity waves reaching the largest observed westward values of −29 m2 s−2 and strong westward wind accelerations of −510 ms−1 day−1. However, despite this strong wave forcing during the event, no unusual filtering conditions or significant increases in wave-excitation proxies were observed. Further, although strong westward wave-induced accelerations were also observed during the 2006 MSAO first phase, there was no corresponding simultaneous response in westward wind. We thus suggest that strong westward fluxes/accelerations of high-frequency gravity waves are not always sufficient to produce anomalous first-phase westward MSAO winds and other forcing may be significant.

Evaluation of momentum flux with radar

The statistics of gravity wave momentum flux estimation are investigated using data from the MU radar at Shigariki, Japan (136°E, 35°N). The radar has been operating during campaign periods since 1986. The first part of the paper focuses on a multi-day campaign during October 13–31, 1986. The second part of the paper investigates data after 2006 when the radar was operated in a meteor scatter mode. Momentum fluxes are derived from both the turbulent scatter and the meteor scatter measurements, but the techniques are quite different. Probability Distribution Functions are formed using turbulent scatter data. These show that wave packets sometimes have momentum flux magnitudes in excess of 100m2s−2. The technique for meteor radars, introduced by Hocking (2005), has been widely adopted by the radar community in recent years. The momentum flux estimated using this technique is found to be anti-correlated with the background tidal winds. A validation investigation is carried out for periods with a high meteor echo data rate. The conclusion was that the method can be used to calculate the sign of momentum flux, but does not accurately specify the magnitude.

Southern Hemisphere Extratropical Gravity Wave Sources and Intermittency Revealed by a Middle-Atmosphere General Circulation Model

Abstract Southern Hemisphere extratropical gravity wave activity is examined using simulations from a free-running middle-atmosphere general circulation model called Kanto that contains no gravity wave parameterizations. The total absolute gravity wave momentum flux (MF) and its intermittency, diagnosed by the Gini coefficient, are examined during January and July. The MF and intermittency results calculated from the Kanto model agree well with results from satellite limb and superpressure balloon observations. The analysis of the Kanto model simulations indicates the following results. Nonorographic gravity waves are generated in Kanto in the frontal regions of extratropical depressions and around tropopause-level jets. Regions with lower (higher) intermittency in the July midstratosphere become more (less) intermittent by the mesosphere as a result of lower-level wave removal. The gravity wave intermittency is low and nearly homogeneous throughout the SH middle atmosphere during January. This indicates that nonorographic waves dominate at this time of year, with sources including continental convection as well as oceanic depressions. Most of the zonal-mean MF at 40°–65°S in January and July is due to gravity waves located above the oceans. The zonal-mean MF at lower latitudes in both months has a larger contribution from the land regions but the fraction above the oceans remains larger.

Intermittency of gravity wave momentum flux in the mesopause region observed with an all‐sky airglow imager

AbstractThe intermittency of gravity wave momentum flux (MF) near the OH airglow layer (∼87 km) in the mesopause region is investigated for the first time using observation of all‐sky airglow imager over Maui, Hawaii (20.7°N, 156.3°W), and Cerro Pachón, Chile (30.3°S, 70.7°W). At both sites, the probability density function (pdf) of gravity wave MF shows two distinct distributions depending on the magnitude of the MF. For MF smaller (larger) than ∼16 m2 s−2 (0.091 mPa), the pdf follows a lognormal (power law) distribution. The intermittency represented by the Bernoulli proxy and the percentile ratio shows that gravity waves have higher intermittency at Maui than at Cerro Pachón, suggesting more intermittent background variation above Maui. It is found that most of the MF is contributed by waves that occur very infrequently. But waves that individually contribute little MF are also important because of their higher occurrence frequencies. The peak contribution is from waves with MF around ∼2.2 m2 s−2 at Cerro Pachón and ∼5.5 m2 s−2 at Maui. Seasonal variations of the pdf and intermittency imply that the background atmosphere has larger influence on the observed intermittency in the mesopause region.

A Grid Transformation Method for a Quasi-Uniform, Circular Fine Region Using the Spring Dynamics

Regionally enhanced meshes that have quasi-uniformly fine circular region is proposed using a new transformation method with icosahedral grids to obtain a cost-effective simulation for waves, transports, and mixing processes, the behaviors of which depend strongly on the horizontal resolution. The target region, which is designed to be composed of a finer mesh, is connected to a coarser mesh region, which is generated with the Schmidt transformation to maintain an isotropy of grid shapes. To realize these requirements, the spring dynamics method can be used and the characteristic length of the spring connecting grid nodes should be determined through three parameters: (i) the number of grid points placed in the target region, (ii) the area of the target region, and (iii) a parameter of the Schmidt transformation. By introducing a set of mathematical formulae, the minimum grid interval in the target region can be uniquely determined as a function of the area of the target region only. It is confirmed that fine and quasi-homogeneous meshes in the target region are generated using the grid transformation proposed in this study. Numerical simulations under realistic atmospheric conditions are performed using a non-hydrostatic model with the grid system proposed in this study and in a previous study. As the new grid system has a more homogenous resolution in the target region compared with that of the previous study, the estimation of the momentum fluxes of gravity waves are less affected by their dependence of the grid resolution.