Browse
Abstract
Idealized integral studies of the dynamics of atmospheric inertia–gravity waves (IGWs) from their sources in the troposphere (e.g., by spontaneous emission from jets and fronts) to dissipation and mean flow effects at higher altitudes could contribute to a better treatment of these processes in IGW parameterizations in numerical weather prediction and climate simulation. It seems important that numerical codes applied for this purpose are efficient and focus on the essentials. Therefore, a previously published staggered-grid solver for f-plane soundproof pseudoincompressible dynamics is extended here by two main components. These are 1) a semi-implicit time stepping scheme for the integration of buoyancy and Coriolis effects, and 2) the incorporation of Newtonian heating consistent with pseudoincompressible dynamics. This heating function is used to enforce a temperature profile that is baroclinically unstable in the troposphere and it allows the background state to vary in time. Numerical experiments for several benchmarks are compared against a buoyancy/Coriolis-explicit third-order Runge–Kutta scheme, verifying the accuracy and efficiency of the scheme. Preliminary mesoscale simulations with baroclinic wave activity in the troposphere show intensive small-scale wave activity at high altitudes, and they also indicate there the expected reversal of the zonal-mean zonal winds.
Abstract
Idealized integral studies of the dynamics of atmospheric inertia–gravity waves (IGWs) from their sources in the troposphere (e.g., by spontaneous emission from jets and fronts) to dissipation and mean flow effects at higher altitudes could contribute to a better treatment of these processes in IGW parameterizations in numerical weather prediction and climate simulation. It seems important that numerical codes applied for this purpose are efficient and focus on the essentials. Therefore, a previously published staggered-grid solver for f-plane soundproof pseudoincompressible dynamics is extended here by two main components. These are 1) a semi-implicit time stepping scheme for the integration of buoyancy and Coriolis effects, and 2) the incorporation of Newtonian heating consistent with pseudoincompressible dynamics. This heating function is used to enforce a temperature profile that is baroclinically unstable in the troposphere and it allows the background state to vary in time. Numerical experiments for several benchmarks are compared against a buoyancy/Coriolis-explicit third-order Runge–Kutta scheme, verifying the accuracy and efficiency of the scheme. Preliminary mesoscale simulations with baroclinic wave activity in the troposphere show intensive small-scale wave activity at high altitudes, and they also indicate there the expected reversal of the zonal-mean zonal winds.
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 southern part of South America and the Antarctic peninsula are known as the world’s strongest hotspot region of stratospheric gravity wave (GW) activity. Large tropospheric winds are deflected by the Andes and the Antarctic Peninsula and excite GWs that might propagate into the upper mesosphere. Satellite observations show large stratospheric GW activity above the mountains, the Drake Passage, and in a belt centered along 60°S. This scientifically highly interesting region for studying GW dynamics was the focus of the Southern Hemisphere Transport, Dynamics, and Chemistry–Gravity Waves (SOUTHTRAC-GW) mission. The German High Altitude and Long Range Research Aircraft (HALO) was deployed to Rio Grande at the southern tip of Argentina in September 2019. Seven dedicated research flights with a typical length of 7,000 km were conducted to collect GW observations with the novel Airborne Lidar for Middle Atmosphere research (ALIMA) instrument and the Gimballed Limb Observer for Radiance Imaging of the Atmosphere (GLORIA) limb sounder. While ALIMA measures temperatures in the altitude range from 20 to 90 km, GLORIA observations allow characterization of temperatures and trace gas mixing ratios from 5 to 15 km. Wave perturbations are derived by subtracting suitable mean profiles. This paper summarizes the motivations and objectives of the SOUTHTRAC-GW mission. The evolution of the atmospheric conditions is documented including the effect of the extraordinary Southern Hemisphere sudden stratospheric warming (SSW) that occurred in early September 2019. Moreover, outstanding initial results of the GW observation and plans for future work are presented.
Abstract
The southern part of South America and the Antarctic peninsula are known as the world’s strongest hotspot region of stratospheric gravity wave (GW) activity. Large tropospheric winds are deflected by the Andes and the Antarctic Peninsula and excite GWs that might propagate into the upper mesosphere. Satellite observations show large stratospheric GW activity above the mountains, the Drake Passage, and in a belt centered along 60°S. This scientifically highly interesting region for studying GW dynamics was the focus of the Southern Hemisphere Transport, Dynamics, and Chemistry–Gravity Waves (SOUTHTRAC-GW) mission. The German High Altitude and Long Range Research Aircraft (HALO) was deployed to Rio Grande at the southern tip of Argentina in September 2019. Seven dedicated research flights with a typical length of 7,000 km were conducted to collect GW observations with the novel Airborne Lidar for Middle Atmosphere research (ALIMA) instrument and the Gimballed Limb Observer for Radiance Imaging of the Atmosphere (GLORIA) limb sounder. While ALIMA measures temperatures in the altitude range from 20 to 90 km, GLORIA observations allow characterization of temperatures and trace gas mixing ratios from 5 to 15 km. Wave perturbations are derived by subtracting suitable mean profiles. This paper summarizes the motivations and objectives of the SOUTHTRAC-GW mission. The evolution of the atmospheric conditions is documented including the effect of the extraordinary Southern Hemisphere sudden stratospheric warming (SSW) that occurred in early September 2019. Moreover, outstanding initial results of the GW observation and plans for future work are presented.
Abstract
We present gravity wave climatologies based on 7 years (2012–18) of lidar and Sounding of the Atmosphere using Broadband Emission Radiometry (SABER) temperatures and reanalysis data at 54° and 69°N in the altitude range 30–70 km. We use 9452 (5044) h of lidar observations at Kühlungsborn [Arctic Lidar Observatory for Middle Atmosphere Research (ALOMAR)]. Filtering according to vertical wavelength (λ z < 15 km) or period (τ < 8 h) is applied. Gravity wave potential energy densities (GWPED) per unit volume (E pV) and per unit mass (E pm) are derived. GWPED from reanalysis are smaller compared to lidar. The difference increases with altitude in winter and reaches almost two orders of magnitude around 70 km. A seasonal cycle of E pV with maximum values in winter is present at both stations in nearly all lidar and SABER measurements and in reanalysis data. For SABER and for lidar (with λ < 15 km) the winter/summer ratios are a factor of ~2–4, but are significantly smaller for lidar with τ < 8 h. The winter/summer ratios are nearly identical at both stations and are significantly larger for E pm compared to E pV. Lidar and SABER observations show that E pV is larger by a factor of ~2 at Kühlungsborn compared to ALOMAR, independent of season and altitude. Comparison with mean background winds shows that simple scenarios regarding GW filtering, etc., cannot explain the Kühlungsborn–ALOMAR differences. The value of E pV decreases with altitude in nearly all cases. Corresponding E pV-scale heights from lidar are generally larger in winter compared to summer. Above ~55 km, E pV in summer is almost constant with altitude at both stations. The winter–summer difference of E pV scale heights is much smaller or absent in SABER and in reanalysis data.
Abstract
We present gravity wave climatologies based on 7 years (2012–18) of lidar and Sounding of the Atmosphere using Broadband Emission Radiometry (SABER) temperatures and reanalysis data at 54° and 69°N in the altitude range 30–70 km. We use 9452 (5044) h of lidar observations at Kühlungsborn [Arctic Lidar Observatory for Middle Atmosphere Research (ALOMAR)]. Filtering according to vertical wavelength (λ z < 15 km) or period (τ < 8 h) is applied. Gravity wave potential energy densities (GWPED) per unit volume (E pV) and per unit mass (E pm) are derived. GWPED from reanalysis are smaller compared to lidar. The difference increases with altitude in winter and reaches almost two orders of magnitude around 70 km. A seasonal cycle of E pV with maximum values in winter is present at both stations in nearly all lidar and SABER measurements and in reanalysis data. For SABER and for lidar (with λ < 15 km) the winter/summer ratios are a factor of ~2–4, but are significantly smaller for lidar with τ < 8 h. The winter/summer ratios are nearly identical at both stations and are significantly larger for E pm compared to E pV. Lidar and SABER observations show that E pV is larger by a factor of ~2 at Kühlungsborn compared to ALOMAR, independent of season and altitude. Comparison with mean background winds shows that simple scenarios regarding GW filtering, etc., cannot explain the Kühlungsborn–ALOMAR differences. The value of E pV decreases with altitude in nearly all cases. Corresponding E pV-scale heights from lidar are generally larger in winter compared to summer. Above ~55 km, E pV in summer is almost constant with altitude at both stations. The winter–summer difference of E pV scale heights is much smaller or absent in SABER and in reanalysis data.
Abstract
Current gravity wave (GW) parameterization (GWP) schemes are using the steady-state assumption, in which an instantaneous balance between GWs and mean flow is postulated, thereby neglecting transient, nondissipative interactions between the GW field and the resolved flow. These schemes rely exclusively on wave dissipation, by GW breaking or near critical layers, as a mechanism leading to forcing of the mean flow. In a transient GWP, without the steady-state assumption, nondissipative wave–mean-flow interactions are enabled as an additional mechanism. Idealized studies have shown that this is potentially important, and therefore the transient GWP Multiscale Gravity Wave Model (MS-GWaM) has been implemented into a state-of-the-art weather and climate model. In this implementation, MS-GWaM leads to a zonal-mean circulation that agrees well with observations and increases GW momentum-flux intermittency as compared with steady-state GWPs, bringing it into better agreement with superpressure balloon observations. Transient effects taken into account by MS-GWaM are shown to make a difference even on monthly time scales: in comparison with steady-state GWPs momentum fluxes in the lower stratosphere are increased and the amount of missing drag at Southern Hemispheric high latitudes is decreased to a modest but nonnegligible extent. An analysis of the contribution of different wavelengths to the GW signal in MS-GWaM suggests that small-scale GWs play an important role down to horizontal and vertical wavelengths of 50 km (or even smaller) and 200 m, respectively.
Abstract
Current gravity wave (GW) parameterization (GWP) schemes are using the steady-state assumption, in which an instantaneous balance between GWs and mean flow is postulated, thereby neglecting transient, nondissipative interactions between the GW field and the resolved flow. These schemes rely exclusively on wave dissipation, by GW breaking or near critical layers, as a mechanism leading to forcing of the mean flow. In a transient GWP, without the steady-state assumption, nondissipative wave–mean-flow interactions are enabled as an additional mechanism. Idealized studies have shown that this is potentially important, and therefore the transient GWP Multiscale Gravity Wave Model (MS-GWaM) has been implemented into a state-of-the-art weather and climate model. In this implementation, MS-GWaM leads to a zonal-mean circulation that agrees well with observations and increases GW momentum-flux intermittency as compared with steady-state GWPs, bringing it into better agreement with superpressure balloon observations. Transient effects taken into account by MS-GWaM are shown to make a difference even on monthly time scales: in comparison with steady-state GWPs momentum fluxes in the lower stratosphere are increased and the amount of missing drag at Southern Hemispheric high latitudes is decreased to a modest but nonnegligible extent. An analysis of the contribution of different wavelengths to the GW signal in MS-GWaM suggests that small-scale GWs play an important role down to horizontal and vertical wavelengths of 50 km (or even smaller) and 200 m, respectively.
Abstract
In a companion paper, the Multiscale Gravity Wave Model (MS-GWaM) has been introduced and its application to a global model as a transient subgrid-scale parameterization has been described. This paper focuses on the examination of intermittency of gravity waves (GWs) modeled by MS-GWaM. To introduce the variability and intermittency in wave sources, convective GW sources are formulated, using diabatic heating diagnosed by the convection parameterization, and they are coupled to MS-GWaM in addition to a flow-independent source in the extratropics accounting for GWs due neither to convection nor to orography. The probability density function (PDF) and Gini index for GW pseudomomentum fluxes are assessed to investigate the intermittency. Both are similar to those from observations in the lower stratosphere. The intermittency of GWs over tropical convection is quite high and is found not to change much in the vertical direction. In the extratropics, where nonconvective GWs dominate, the intermittency is lower than that in the tropics in the stratosphere and comparable to that in the mesosphere, exhibiting a gradual increase with altitude. The PDFs in these latitudes seem to be close to the lognormal distributions. Effects of transient GW–mean-flow interactions on the simulated GW intermittency are assessed by performing additional simulations using the steady-state assumption in the GW parameterization. The intermittency of parameterized GWs over tropical convection is found to be overestimated by the assumption, whereas in the extratropics it is largely underrepresented. Explanation and discussion of these effects are included.
Abstract
In a companion paper, the Multiscale Gravity Wave Model (MS-GWaM) has been introduced and its application to a global model as a transient subgrid-scale parameterization has been described. This paper focuses on the examination of intermittency of gravity waves (GWs) modeled by MS-GWaM. To introduce the variability and intermittency in wave sources, convective GW sources are formulated, using diabatic heating diagnosed by the convection parameterization, and they are coupled to MS-GWaM in addition to a flow-independent source in the extratropics accounting for GWs due neither to convection nor to orography. The probability density function (PDF) and Gini index for GW pseudomomentum fluxes are assessed to investigate the intermittency. Both are similar to those from observations in the lower stratosphere. The intermittency of GWs over tropical convection is quite high and is found not to change much in the vertical direction. In the extratropics, where nonconvective GWs dominate, the intermittency is lower than that in the tropics in the stratosphere and comparable to that in the mesosphere, exhibiting a gradual increase with altitude. The PDFs in these latitudes seem to be close to the lognormal distributions. Effects of transient GW–mean-flow interactions on the simulated GW intermittency are assessed by performing additional simulations using the steady-state assumption in the GW parameterization. The intermittency of parameterized GWs over tropical convection is found to be overestimated by the assumption, whereas in the extratropics it is largely underrepresented. Explanation and discussion of these effects are included.
Abstract
The way the large-scale flow determines the energy of the nonorographic mesoscale inertia–gravity waves (IGWs) is theoretically significant and practically useful for source parameterization of IGWs. The relations previously developed on the f plane for tropospheric sources of IGWs including jets, fronts, and convection in terms of associated secondary circulations strength are generalized for application over the globe. A low-pass spatial filter with a cutoff zonal wavenumber of 22 is applied to separate the large-scale flow from the IGWs using the ERA5 data of ECMWF for the period 2016–19. A comparison with GRACILE data based on satellite observations of the middle stratosphere shows reasonable representation of IGWs in the ERA5 data despite underestimates by a factor of smaller than 3. The sum of the energies, which are mass-weighted integrals in the troposphere from the surface to 100 hPa, as given by the generalized relations is termed initial parameterized energy. The corresponding energy integral for the IGWs is termed the diagnosed energy. The connection between the parameterized and diagnosed IGW energies is explored with regression analysis for each season and six oceanic domains distributed over the globe covering the Northern and Southern Hemispheres and the tropics. While capturing the seasonal cycle, the domain area-average seasonal mean initial parameterized energy is weaker than the diagnosed energy by a factor of 3. The best performance in regression analysis is obtained by using a combination of power and exponential functions, which suggests evidence of exponential weakness.
Abstract
The way the large-scale flow determines the energy of the nonorographic mesoscale inertia–gravity waves (IGWs) is theoretically significant and practically useful for source parameterization of IGWs. The relations previously developed on the f plane for tropospheric sources of IGWs including jets, fronts, and convection in terms of associated secondary circulations strength are generalized for application over the globe. A low-pass spatial filter with a cutoff zonal wavenumber of 22 is applied to separate the large-scale flow from the IGWs using the ERA5 data of ECMWF for the period 2016–19. A comparison with GRACILE data based on satellite observations of the middle stratosphere shows reasonable representation of IGWs in the ERA5 data despite underestimates by a factor of smaller than 3. The sum of the energies, which are mass-weighted integrals in the troposphere from the surface to 100 hPa, as given by the generalized relations is termed initial parameterized energy. The corresponding energy integral for the IGWs is termed the diagnosed energy. The connection between the parameterized and diagnosed IGW energies is explored with regression analysis for each season and six oceanic domains distributed over the globe covering the Northern and Southern Hemispheres and the tropics. While capturing the seasonal cycle, the domain area-average seasonal mean initial parameterized energy is weaker than the diagnosed energy by a factor of 3. The best performance in regression analysis is obtained by using a combination of power and exponential functions, which suggests evidence of exponential weakness.
Abstract
Inertia–gravity waves (IGWs) play an essential role in the terrestrial atmospheric dynamics as they can lead to energy and momentum flux when propagating upward. An open question is to what extent IGWs contribute to the total energy and to the flattening of the energy spectrum observed at the mesoscale. In this work, we present an experimental investigation of the energy distribution between the large-scale balanced flow and the small-scale imbalanced flow. Weakly nonlinear IGWs emitted from baroclinic jets are observed in the differentially heated rotating annulus experiment. Similar to the atmospheric spectra, the experimental kinetic energy spectra reveal the typical subdivision into two distinct regimes with slopes k −3 for the large scales and k −5/3 for the small scales. By separating the spectra into the vortex and wave components, it emerges that at the large-scale end of the mesoscale the gravity waves observed in the experiment cause a flattening of the spectra and provide most of the energy. At smaller scales, our data analysis suggests a transition toward a turbulent regime with a forward energy cascade up to where dissipation by diffusive processes occurs.
Abstract
Inertia–gravity waves (IGWs) play an essential role in the terrestrial atmospheric dynamics as they can lead to energy and momentum flux when propagating upward. An open question is to what extent IGWs contribute to the total energy and to the flattening of the energy spectrum observed at the mesoscale. In this work, we present an experimental investigation of the energy distribution between the large-scale balanced flow and the small-scale imbalanced flow. Weakly nonlinear IGWs emitted from baroclinic jets are observed in the differentially heated rotating annulus experiment. Similar to the atmospheric spectra, the experimental kinetic energy spectra reveal the typical subdivision into two distinct regimes with slopes k −3 for the large scales and k −5/3 for the small scales. By separating the spectra into the vortex and wave components, it emerges that at the large-scale end of the mesoscale the gravity waves observed in the experiment cause a flattening of the spectra and provide most of the energy. At smaller scales, our data analysis suggests a transition toward a turbulent regime with a forward energy cascade up to where dissipation by diffusive processes occurs.
