Vertical Wave Propagation Research Articles

Simulation of Propagation of Acoustic-Gravity Waves Generated by Tropospheric Front Instabilities into the Upper Atmosphere

A numerical model of the propagation of acoustic-gravity waves excited by pressure fluctuations on the Earth’s surface is developed. Propagation of acoustic-gravity waves generated by instabilities of tropospheric fronts into the upper atmosphere is simulated. The experimental data on atmospheric pressure variations during 2016 year registered on a net of four microbarographs located in the Moscow region are processed. A case of very significant pressure fluctuations (up to 30 times larger than the average level) is selected, which were caused by an atmospheric front arrival. Observed surface pressure field variations for this field were approximated and used as the lower boundary condition for simulating the vertical wave propagation. The numerical simulations showed that just after the boundary source activation, the infrasonic waves in the upper atmosphere may have amplitudes of perturbations of temperature up to 100 K, and horizontal velocity up to 60 m/s. Internal gravity waves come into the upper atmosphere later and far horizontally away from the wave source. The influence of the limited dimensions of the computational domain on the simulation results is investigated. The conditions at the horizontal boundaries of the computational domain, which allow the runaway of waves beyond the domain are proposed. The frequency spectrum of waves in the non-isothermal atmosphere is analyzed.

Recent weakening of the southern stratospheric polar vortex and its impact on the surface climate over Antarctica

The variability in the southern stratospheric polar vortex (SSPV) and its downward coupling with the troposphere are known to play a crucial role in driving climate variability over Antarctica. In this study, SSPV weakening events and their impacts on the surface climate of Antarctica are examined using in-situ observation and reanalysis data. Combining criteria from several previous studies, we introduce a new detection method for SSPV weakening events. Based on the new criteria, the occurrence frequency of SSPV weakening events has exhibited a systematic increasing trend since the 2000 s. However, the weakened anomalies of individual SSPV events are not statistically different (95% confidence level) between the earlier (1979–1999) and later (2000–2017) periods examined in this study. The recent increase in the occurrence of SSPV weakening events is largely controlled by tropospheric mechanisms, i.e. the poleward heat flux carried by southern hemispheric planetary waves and associated vertical wave propagation. Among the various scales of planetary waves, the wavenumber 1 contributes most of the poleward eddy heat flux. We show that SSPV weakening events induce statistically significant cooling over the Antarctic Peninsula (AP) region and warming over the rest of Antarctica. Typically, surface air temperature anomalies with large negative values smaller than − 0.6 °C and positive values larger than + 0.8 °C are observed over the east coast of the tip of the AP and King Edward VII Land, respectively. The influence of an SSPV weakening event on the surface lasts for approximately three months with higher height anomalies off western Antarctica, providing favorable conditions for the atmosphere to transport cold air from the interior of Antarctica to the AP via the Weddell Sea. Distinct positive surface air temperature anomalies over the rest of Antarctica are associated with the northerly circulation anomaly from the eastern Weddell Sea to east Antarctica.

A Global Nonhydrostatic Atmospheric Model with a Mass- and Energy-conserving Vertically Implicit Correction (VIC) Scheme

Abstract Global nonhydrostatic atmospheric models are becoming increasingly important for studying the climates of planets and exoplanets. However, such models suffer from computational difficulties due to the large aspect ratio between the horizontal and vertical directions. To overcome this problem, we developed a global model using a vertically implicit correction (VIC) scheme in which the integration time step is no longer limited by the vertical propagation of acoustic waves. We proved that our model, based on the Athena++ framework and its extension for planetary atmospheres—SNAP (Simulating Nonhydrostatic Atmospheres on Planets), rigorously conserves mass and energy in finite-volume simulations. We found that traditional numerical stabilizers such as hyperviscosity and divergence damping are not needed when using the VIC scheme, which greatly simplifies the numerical implementation and improves stability. We present simulation results ranging from 1D linear waves to 3D global circulations with and without the VIC scheme. These tests demonstrate that our formulation correctly tracks local turbulent motions, produces Kelvin–Helmholtz instability, and generates a super-rotating jet on hot Jupiters. Employing this VIC scheme improves the computational efficiency of global simulations by more than two orders of magnitude compared to an explicit model and facilitates the capability of simulating a wide range of planetary atmospheres both regionally and globally.

Influence of Concentration on Sedimentation of a Dense Suspension in a Viscous Fluid

Macroscopic properties of sedimenting suspensions have been studied extensively and can be characterized using the Galileo number (Ga), solid-to-fluid density ratio (pi _p) and mean solid volume concentration ({bar{phi }}). However, the particle–particle and particle–fluid interactions that dictate these macroscopic trends have been challenging to study. We examine the effect of concentration on the structure and dynamics of sedimenting suspensions by performing direct numerical simulation based on an Immersed Boundary Method of monodisperse sedimenting suspensions of spherical particles at fixed Ga=144, pi _p=1.5, and concentrations ranging from {bar{phi }}=0.5 to {bar{phi }}=30%. The corresponding particle terminal Reynolds number for a single settling particle is Re_T = 186. Our simulations reproduce the macroscopic trends observed in experiments and are in good agreement with semi-empirical correlations in literature. From our studies, we observe, first, a change in trend in the mean settling velocities, the dispersive time scales and the structural arrangement of particles in the sedimenting suspension at different concentrations, indicating a gradual transition from a dilute regime ({bar{phi }} lesssim 2%) to a dense regime ({bar{phi }} gtrsim 10%). Second, we observe the vertical propagation of kinematic waves as fluctuations in the local horizontally-averaged concentration of the sedimenting suspension in the dense regime.

Vertical Propagation of Acoustic Gravity Waves from the Lower Atmosphere during a Solar Eclipse

This paper presents the results of the analysis of a numerical experiment on the generation of acoustic gravity waves (AGWs) in the lower atmosphere and their propagation into the upper atmosphere. The observed variations in atmospheric pressure during the passage of a solar eclipse on March 20, 2015, over Kaliningrad (54° N, 20° E) are considered the be a wave source. The calculation results showed that the passage of a solar eclipse is accompanied in the lower atmosphere by an increase in wave activity in the AGW range with periods of ~4–20 min. In the upper atmosphere, at altitudes of ~200 km, disturbances are formed due to the dissipation of AGWs coming from the lower atmosphere; these disturbances propagate with characteristic periods of ~50 min. The propagation direction of these disturbances was opposite to the propagation direction of the solar eclipse region along the Earth’s surface. The results of the numerical experiment are in qualitative agreement with the results of observations of ionospheric disturbances during the solar eclipse on March 20, 2015.

Modulation of the Westerly and Easterly Quasi-Biennial Oscillation Phases on the Connection between the Madden–Julian Oscillation and the Arctic Oscillation

Previous studies have revealed the relationship between the Madden–Julian oscillation (MJO) and the Arctic Oscillation (AO). The MJO phase 2/3 is followed by the positive AO phase, and the MJO phase 6/7 is followed by the negative AO phase. This study reveals that the MJO phase 6/7–AO connection is modulated by the Quasi-Biennial Oscillation (QBO) through both tropospheric and stratospheric pathways during boreal winter. The MJO 2/3 phase and AO relationship is favored in both QBO easterly (QBOE) and westerly (QBOW) years because of the MJO-triggered tropospheric Rossby wave train from the tropics toward the polar region. The AO following the MJO 6/7 phase shifts to negative in QBOW years, but the MJO–AO connection diminishes in QBOE years. In QBOW years, the Asian-Pacific jet is enhanced, leading to more evident poleward propagation of tropospheric Rossby wave train, which contributes to the tropospheric pathway of the AO–MJO 6/7 connection. Besides, the enhanced Asian-Pacific jet in QBOW years is favorable for vertical propagation of planetary waves into the stratosphere in MJO phase 6/7, leading to negative AO, which indicates the stratospheric pathway of the AO–MJO 6/7 connection.

Influence of the cloud-level neutral layer on the vertical propagation of topographically generated gravity waves on Venus

The reason for stationary gravity waves at Venus’ cloud top to appear mostly at low latitudes in the afternoon is not understood. Since a neutral layer exists in the lower part of the cloud layer, the waves should be affected by the neutral layer before reaching the cloud top. To what extent gravity waves can propagate vertically through the neutral layer has been unclear. To examine the possibility that the variation of the neutral layer thickness is responsible for the dependence of the gravity wave activity on the latitude and the local time, we investigated the sensitivity of the vertical propagation of gravity waves on the neutral layer thickness using a numerical model. The results showed that stationary gravity waves with zonal wavelengths longer than 1000 km can propagate to the cloud-top level without notable attenuation in the neutral layer with realistic thicknesses of 5–15 km. This suggests that the observed latitudinal and local time variation of the gravity wave activity should be attributed to processes below the cloud. An analytical approach also showed that gravity waves with horizontal wavelengths shorter than tens of kilometers would be strongly attenuated in the neutral layer; such waves should originate in the altitude region above the neutral layer.

Intermodel Spread in the Northern Hemisphere Stratospheric Polar Vortex Response to Climate Change in the CMIP5 Models

AbstractThis study investigates the intermodel spread of the Northern Hemisphere winter stratospheric polar vortex change to anthropogenic greenhouse gas increase in Coupled Model Intercomparison Project Phase 5 (CMIP5) models. Previous proposed mechanisms for the polar vortex response to climate change, based on analysis of atmosphere‐only models, are found inadequate to explain the intermodel spread in the coupled models in CMIP5. It is further found that resolved stationary wave driving in the polar vortex region accounts for less than 30% of the intermodel spread, and intermodel differences in both the vertical and meridional wave propagation contribute to differences in the wave driving. The results call for a detailed budget analysis of the stratospheric circulation response by including both the resolved and parameterized processes through the Dynamics and Variability Model Intercomparison Project. The results also highlight a need for an improved theoretical understanding of future projected polar vortex change and intermodel spread.

Stratospheric Initial Condition for Skillful Surface Prediction in the ECMWF Model

AbstractThe subseasonal‐to‐seasonal prediction of boreal‐winter sea‐level pressure (SLP) is investigated by examining 19‐year reforecast data sets from the European Centre for Medium‐Range Weather Forecasts ensemble prediction system. Based on a pattern anomaly correlation coefficient of 0.6, the model reliably predicted the extratropical SLP anomalies with a maximum forecast lead time of 9.8 days. This prediction skill, however, varies significantly from one forecast to another. The forecasts with relatively high skills and those with low skills exhibit significant zonal Wavenumber 1 patterns in the stratospheric initial conditions with an opposite sign. The high‐skill forecasts present a relatively weak westward Wavenumber 1 phase tilt with height compared with the low‐skill forecasts, dampening an upward wave propagation in the stratosphere. This result suggests that not only the zonal mean but also the zonally asymmetric stratospheric conditions play an important role in extended surface prediction by modulating vertical wave propagation or wave reflection.

Observed Near Inertial Waves in the Wake of Typhoon Linfa (2015) in the Northern South China Sea

This study examined the characteristics and vertical propagation of near inertial waves (NIWs) induced by Typhoon Linfa (2015), based on in situ observations conducted southeast of Dongsha Islands in the South China Sea. The results demonstrate that the near inertial currents induced by Linfa had velocities up to 35cm s−1 in the mixed layer and 20cm s−1 in the ocean interior. The near inertial currents were polarized with predominantly clockwise-rotating components, the magnitudes of which were about 10 times larger than the counter-clockwise rotating components. The energy density spectrum showed that the emergence of NIWs resulted in energy redistribution from the diurnal band to the near inertial band. The wavenumber spectrum and the downward/ upward current decomposition demonstrated that the NIWs and energy flux propagated mainly downward. The estimated vertical phase velocity and group velocity are 1.44 and 0.48m h−1, respectively, corresponding to a vertical wavelength of 49.7m. The e-folding time scale was 7.5 d based on the near inertial kinetic energy in the ocean interior. We found no obvious wave-wave interaction during the decay process of the NIWs. The frequency was blue-shifted, being 0.03f0 higher than the local inertial frequency, which was caused by the background vorticity. The normal mode analysis suggests that the higher mode plays a dominant role in the propagation stage of the NIWs.

Middle atmospheric planetary waves in contrasting QBO phases over the Indian low latitude region

The present study primarily focused on the characteristics of the planetary scale waves in the middle atmosphere during different phases of Quasi-Biennial Oscillation (QBO), a dominant oscillation in the low-latitude stratospheric region. The temperature profiles retrieved from the Rayleigh lidar measurements have been utilized for the 11 winter periods between 1998 and 2009 over a low-latitude station, Gadanki (13.5oN 79.2oE). The spectral analysis of temperature anomalies indicates two dominant planetary-scale modes, namely, quasi 12-day and quasi 16-day waves. The existence of these waves in the middle atmosphere is strongly controlled by the westerly and easterly phases of QBO. For instance, the 12-day wave is mainly observed in the QBO westerly phase, while the 16-day wave peaks at two heights; 30–40 km and above 60 km with large spread in the mesosphere due to Doppler shifting in presence of westerly winds. The QBO easterly phase indicates low wave activity with 16-day wave indicating appreciable amplitudes in the upper stratosphere and mesosphere and is absent in the lower altitudes. This indicates that the 16-day wave might be generated through some in-situ mechanism due to gravity breaking or instability in the mesospheric region. Moreover, the refractive index of the two dominant planetary waves are strongly negative in the easterly phase relative to the westerly phase of QBO in the lower troposphere and stratosphere. This indicates the high probability of the planetary wave vertical propagation in the westerly phase of QBO. Therefore, the presented report re-emphasizes the importance of QBO controlling the middle atmospheric dynamics through vertical propagation of planetary scale waves in low-latitudes.

Vertical Propagation of Acoustic-Gravity Waves from Atmospheric Fronts into the Upper Atmosphere

Vertical Propagation of Acoustic-Gravity Waves from Atmospheric Fronts into the Upper Atmosphere

Site dependent and spatially varying response spectra

The goal of this study is to provide a stochastic method to investigate the effects of the randomness of soil properties due to their natural spatial variability on the response spectra spatial variation at sites with varying conditions. For this purpose, Monte Carlo Simulations are used to include the variability of both incident ground motion and soil parameters in the response spectra by mean of an appropriate coherency loss function and a site-dependent transfer function, respectively. The approach is built on the assumption of vertical propagation of SH type waves in soil strata with uncertain parameters. The response spectra are obtained by numerical integration of the governing equation of a single-degree-of-freedom (SDOF) system under non-stationary site-dependent and spatially varying ground motion accelerations simulated with non-uniform spectral densities and coherency loss functions. Numerical examples showed that randomness of soil properties significantly affects the amplitudes of the response spectra, indicating that as the heterogeneity induced by the randomness of the parameters of the medium increases, the spectral ordinates attenuate.

Propagation to the upper atmosphere of acoustic-gravity waves from atmospheric fronts in the Moscow region

Abstract. The paper uses experimental data of pressure variations on the Earth's surface during the passage of an atmospheric front recorded by a network of four microbarographs in the Moscow region. Applying these experimental data, empirical approximations of atmospheric pressure field oscillations are suggested. The obtained approximating surface pressure functions are used as the lower boundary condition for simulating the vertical propagation of acoustic-gravity waves from a source in the lower troposphere. Estimates of the amplitude of temperature disturbances in the upper atmosphere caused by acoustic-gravity waves from a propagating atmospheric front are obtained. For the amplitude of wave temperature disturbances, values of about 200 K are obtained. The amplitude of temperature disturbances in the upper atmosphere caused by background pressure fluctuations on the Earth's surface is estimated at 4–5 K.

100 Years of Progress on Mountain Meteorology Research

ABSTRACTMountains significantly influence weather and climate on Earth, including disturbed surface winds; altered distribution of precipitation; gravity waves reaching the upper atmosphere; and modified global patterns of storms, fronts, jet streams, and climate. All of these impacts arise because Earth’s mountains penetrate deeply into the atmosphere. This penetration can be quantified by comparing mountain heights to several atmospheric reference heights such as density scale height, water vapor scale height, airflow blocking height, and the height of natural atmospheric layers. The geometry of Earth’s terrain can be analyzed quantitatively using statistical, matrix, and spectral methods. In this review, we summarize how our understanding of orographic effects has progressed over 100 years using the equations for atmospheric dynamics and thermodynamics, numerical modeling, and many clever in situ and remote sensing methods. We explore how mountains disturb the surface winds on our planet, including mountaintop winds, severe downslope winds, barrier jets, gap jets, wakes, thermally generated winds, and cold pools. We consider the variety of physical mechanisms by which mountains modify precipitation patterns in different climate zones. We discuss the vertical propagation of mountain waves through the troposphere into the stratosphere, mesosphere, and thermosphere. Finally, we look at how mountains distort the global-scale westerly winds that circle the poles and how varying ice sheets and mountain uplift and erosion over geologic time may have contributed to climate change.

Robustness of the stratospheric pathway in linking the Barents-Kara Sea sea ice variability to the mid-latitude circulation in CMIP5 models

This study investigates the robustness of the stratospheric pathway in linking the sea ice variability over the Barents-Kara Sea in late autumn and early winter to the mid-latitude circulation in the subsequent winter. Two groups of models from the Coupled Model Intercomparison Project phase 5 (CMIP5) archive, one with a well-resolved stratosphere (high-top models) and the other with a poorly-resolved stratosphere (low-top models) are explored to distinguish the role of the stratospheric pathway. The results show that, collectively, high-top models are able to capture the persistent mid-latitude circulation response in the subsequent winter. The response in low-top models is, however, weaker and not as long-lasting most likely due to lack of stratospheric variability. Analysis of eddy heat flux reveals that stronger vertical wave propagation leads to a stronger response in stratospheric polar vortex in high-top models. In particular, it shows that zonal wave-2 eddy heat flux is crucial in leading to a stronger linear constructive interference with the climatological waves in high-top models. The results find that multi-model ensemble of CMIP5 high-top models is able to capture the prolonged impact of sea ice variability on the mid-latitude circulation and out performs the low-top models in this regard. Our study suggests that the representation of the stratosphere in climate models plays an important role in amplifying and extending the mid-latitude circulation response.

Sources, Sinks, and Propagation Characteristics of the Quasi 6‐Day Wave and Its Impact on the Residual Mean Circulation

AbstractThis study employs a troposphere to lower thermosphere assimilation model data set generated by the Whole Atmosphere Community Climate Model with data assimilation provided by the Data Assimilation Research Test Bed (WACCM + DART) to explore the sources, sinks, and propagation characteristics of the quasi 6‐day wave (Q6DW) in the year 2007. WACCM + DART reproduces the burst‐like Q6DW and compares well with Sounding of the Atmosphere using Broadband Emission Radiometry and Thermosphere, Ionosphere, Mesosphere Energetics and Dynamics Doppler Interferometer observations. The most prominent Q6DW took place in later February and mid‐October, while the Q6DW was absent during solstice conditions in 2007. The occurrence of a large Q6DW in the equinoctial mesosphere and lower thermosphere is highly dependent on wave amplification and overreflection processes associated with barotropic/baroclinic instabilities and wave critical layers defined by the zonal mean zonal winds. During solstices, the winter hemisphere waveguide is negative and prevents the vertical wave propagation from the source region into the mesosphere and lower thermosphere. Meanwhile, the critical layer for the Q6DW encloses the unstable region in the summer hemisphere and thus blocks the energy conversion from the mean flow to the wave. The resulting circulation pattern due to the Q6DW momentum deposition is upward and poleward in both hemispheres and thus weakens the residual mean circulation in the summer hemisphere but strengthens it in the winter hemisphere. Also, the Q6DW impact on the residual mean circulation points to broader implications for the mean state of the upper atmosphere, for example, the thermospheric O/N2 ratio due to upward constituent transport and related changes in the ionospheric plasma.

Mantle Seismic Anisotropy beneath the Amur Plate According to the Data of ScS Waves from Deep-Focus Earthquakes

Mantle seismic anisotropy beneath the Amur Plate is studied using the data of ScS waves reflected from the outer core of the Earth from local deep earthquakes in the area of five stations located in Primorye and Priamurye regions. The results of measuring the parameters of the split ScS waves near the stations show the dominance of the polarization azimuths of the fast wave along the eastern directions and are consistent with azimuthal anisotropy and the direction of motion of the Pacific Plate (300°) and Amur Plate (~120°) depending on the epicenter–station direction. In Primorye, in the vicinity of the TEY station, the polarization azimuths of the fast ScS wave dominate in the interval of NE–E directions orthogonally to the lines of the mantle flow along the complex 3D surface of the subsiding Pacific Plate. It is revealed that the delay time of ScS waves increases to 2 s and 3.4 s in the upper mantle and in the mantle transition zone, respectively, as the depth of events increases. The highest degree of anisotropy manifests itself in the upper mantle. In the case of the vertical propagation of waves under conditions of horizontal mantle flow, the difference in the arrival times of the ScS waves is the lowest. The anisotropy in the upper zone of the transition mantle can be related to the wadsleyite texture with the polarization of the fast ScS wave parallel or orthogonally to the motion of the stagnant plate. The anisotropy in the upper part of the lower mantle is associated with the texture of perovskite and periclase with the orientation of the symmetry axis and the polarization of the ScS waves in parallel to the plate subsidence and the subduction direction.

Vertically Propagating Mountain Waves—A Hazard for High-Flying Aircraft?

AbstractStall warnings at flight level 410 (12.5 km) occurred unexpectedly during a research flight of the High Altitude and Long Range Research Aircraft (HALO) over Italy on 12 January 2016. The dangerous flight situation was mitigated by pilot intervention. At the incident location, the stratosphere was characterized by large horizontal variations in the along-track wind speed and temperature. On this particular day, strong northwesterly winds in the lower troposphere in concert with an aligned polar front jet favored the excitation and vertical propagation of large-amplitude mountain waves at and above the Apennines in Italy. These mountain waves carried large vertical energy fluxes of 8 W m−2 and propagated without significant dissipation from the troposphere into the stratosphere. While turbulence is a well-acknowledged hazard to aviation, this case study reveals that nonbreaking, vertically propagating mountain waves also pose a potential hazard, especially to high-flying aircraft. It is the wave-induced modulation of the ambient along-track wind speed that may decrease the aircraft speed toward the minimum needed stall speed.

Vertical propagation of acoustic waves in the solar internetworkas observed by IRIS

We investigate the Interface Region Imaging Spectrograph (IRIS) observations of the quiet-Sun (QS) to understand the propagation of acoustic waves in transition region (TR) from photosphere. We selected a few IRIS spectral lines, which include the photospheric (Mn~{\sc i} 2801.25~{\AA}), chromospheric (Mg~{\sc ii} k 2796.35~{\AA}) and TR (C~{\sc ii} 1334.53~{\AA}), to investigate the acoustic wave propagation.The wavelet cross-spectrum reveals significant coherence (about 70\% locations) between photosphere and chromosphere. Few minutes oscillations (i.e., period range from 1.6 to 4.0 minutes) successfully propagate into chromosphere from photosphere, which is confirmed by dominance of positive phase lags. However, in higher period regime (i.e., greater than $\approx$ 4.5 minutes), the downward propagation dominates is evident by negative phase lags. The broad spectrum of waves (i.e., 2.5-6.0 minutes) propagates freely upwards from chromosphere to TR. We find that only about 45\% locations (out of 70\%) show correlation between chromosphere and TR. Our results indicate that roots of 3 minutes oscillations observed within chromosphere/TR are located in photosphere. Observations also demonstrate that 5 minute oscillations propagate downward from chromosphere. \textbf{However, some locations within QS also show successful propagation of 5 minute oscillations as revealed by positive phase lags, which might be the result of magnetic field}. In addition, our results clearly show that a significant power, within period ranging from 2.5 to 6.0 minutes, of solar chromosphere is freely transmitted into TR triggering atmospheric oscillations. Theoretical implications of our observational results are discussed.