Numerical modeling of inhomogeneous orographic wave influence on planetary waves in the middle atmosphere

Parameterization of dynamical and thermal effects of stationary orographic gravity waves (OGWs) generated by the Earth’s surface topography is incorporated into a numerical model of general circulation of the middle and upper atmosphere. Responses of atmospheric general circulation and characteristics of planetary waves at altitudes from the troposphere up to the thermosphere to the effects of OGWs propagating from the earth surface are studied. Changes in atmospheric circulation and amplitudes of planetary waves due to variations of OGW generation and propagation in different seasons are considered. It is shown that during solstices the main OGW dynamical and heat effects occur in the middle atmosphere of winter hemispheres, where changes in planetary wave amplitudes due to OGWs may reach up to 50%. During equinoxes OGW effects are distributed more homogeneously between northern and southern hemispheres.     

Simulating planetary wave propagation to the upper atmosphere during stratospheric warming events at different mountain wave scenarios

Parameterization schemes of atmospheric normal modes (NMs) and orographic gravity waves (OGWs) have been implemented into the mechanistic Middle and Upper Atmosphere Model (MUAM) simulating atmospheric general circulation. Based on the 12-members ensemble of runs with the MUAM, a composite of the stratospheric warming (SW) has been constructed using the UK Met Office data as the lower boundary conditions. The simulation results show that OGW amplitudes increase at altitudes above 30 km in the Northern Hemisphere after the SW event. At altitudes of about 50 km, OGWs have largest amplitudes over North American and European mountain systems before and during the composite SW, and over Himalayas after the SW. Simulations demonstrate substantial (up to 50–70%) variations of amplitudes of stationary planetary waves (PWs) during and after the SW in the mesosphere-lower thermosphere of the Northern Hemisphere. Westward travelling NMs have amplitude maxima not only in the Northern, but also in the Southern Hemisphere, where these modes have waveguides in the middle and upper atmosphere. Simulated variations of PW and NM amplitudes correspond to changes in the mean zonal wind, EP-fluxes and wave refractive index at different phases of the composite SW events. Inclusion of the parameterization of OGW effects leads to decreases in amplitudes (up to 15%) of almost all SPWs before and after the SW event and their increase (up to 40–60%) after the SW in the stratosphere and mesosphere at middle and high northern latitudes. It is suggested that observed changes in NM amplitudes in the Southern Hemisphere during SW could be caused by divergence of increased southward EP-flux. This EP-flux increases due to OGW drag before SW and extends into the Southern Hemisphere.     

A numerical model characterizing internal gravity wave propagation into the upper atmosphere

A new two-dimensional, time-dependent and fully nonlinear model is developed to numerically simulate plane wave motions for internal gravity waves in a non-isothermal and windy atmosphere that accounts for the dissipation due to eddy and molecular processes. The atmosphere has been treated as a well mixed total gas with a constant mean molecular weight. The effects of Rayleigh friction and Newtonian cooling are applied near the upper boundary of the model to simulate the radiation conditions as well as to act as a sponge layer; lateral boundaries are periodic over a horizontal wavelength to simulate a horizontally infinite domain. The thermal excitation to initiate upward propagating waves is spatially localized in the troposphere and is a Gaussian function of time. A time-splitting technique is applied to the finite difference equations that are derived from the Navier–Stokes equations. The time integration for these highly coupled equations is performed using an explicit second order Lax–Wendroff scheme and an implicit Newton–Raphson scheme. The wave solutions derived from the model are found to be broadly agreeable with those derived from a Wentzel–Kramers–Brillouin theory.     

Semi-empirical model of middle atmosphere wind from the ground to the lower thermosphere

During recent years, special attention has been paid to understanding the background circulation of the middle atmosphere. Particularly in the mesosphere/lower thermosphere (MLT) region, this has involved including data from a range of new radar measurements. It has also involved the comparison of existing empirical middle atmosphere wind models, such as CIRA-86 and HWM-93 to the new data. This has led to the construction of empirical models of MLT winds such as the Global Empirical Wind Model (GEWM). Further investigations are aimed at the construction of new empirical and semi-empirical wind models of the entire middle atmosphere including these new experimental results. The results of a new wind climatology (0–100 km) are presented here, based upon the GEWM, a reanalysis of stratospheric data, and a numerical model which is used to fill the gap between data from the stratospheric and MLT regions.     

Response of the extratropical middle atmosphere to the September 2002 major stratospheric sudden warming

The effects of a major stratospheric sudden warming (SSW) at extratropical latitudes have been investigated with wind and temperature observations over a Brazilian station, Cachoeira Paulista (22.7°S, 45°W) during September–October 2002. In response to the warming at polar latitudes a corresponding cooling at tropical and extratropical latitudes is prominent in the stratosphere. A conspicuous signature of latitudinal propagation of a planetary wave of zonal wavenumbers 1 and 2 from polar to low latitude has been observed during the warming period. The polar vortex which split into two parts of different size is found to travel considerably low latitude. Significant air mass mixing between low and high latitudes is caused by planetary wave breaking. The meridional wind exhibits oscillations of period 2–4 days during the warming period in the stratosphere. No wave feature is evident in the mesosphere during the warming period, although a 12–14 day periodicity is observed after 2 weeks of the warming event, indicating close resemblance to the results of other simultaneous investigations carried out from high latitude Antarctic stations. Convective activity over the present extratropical station diminishes remarkably during the warming period. This behavior is possibly due to destabilization and shift of equatorial convective active regions towards the opposite hemisphere in response to changes in the mean meridional circulation in concert with the SSW.     

Regional 4-D modeling of the ionospheric electron density from satellite data and IRI

Accurate knowledge of the electron density is the key point in correcting ionospheric delays of electromagnetic measurements and in studying ionosphere physics. During the last decade Global Navigation Satellite Systems (GNSS) have become a promising tool for monitoring ionospheric parameters such as the total electron content (TEC). In this contribution we present a four-dimensional (4-D) model of the electron density consisting of a given reference part, i.e., the International Reference Ionosphere (IRI), and an unknown correction term expanded in terms of multi-dimensional base functions. The corresponding series coefficients are calculable from the satellite measurements by applying parameter estimation procedures. Since satellite data are usually sampled between GPS satellites and ground stations, finer structures of the electron density are modelable just in regions with a sufficient number of ground stations. The proposed method is applied to simulated geometry-free GPS phase measurements. The procedure can be used, for example, to study the equatorial anomaly.     

Response of low latitude D-region ionosphere to the total solar eclipse of 22 July 2009 deduced from ELF/VLF analysis

Response of the D-region of the ionosphere to the total solar eclipse of 22 July 2009 at low latitude, Varanasi (Geog. lat., 25.27° N; Geog. long., 82.98° E; Geomag. lat. = 14° 55’ N) was investigated using ELF/VLF radio signal. Tweeks, a naturally occurring VLF signal and radio signals from various VLF navigational transmitters are first time used simultaneously to study the effect of total solar eclipse (TSE). Tweeks occurrence is a nighttime phenomena but the obscuration of solar disc during TSE in early morning leads to tweek occurrence. The changes in D-region ionospheric VLF reflection heights (h) and electron density (n_e: 22.6–24.6 cm⁻³) during eclipse have been estimated from tweek analysis. The reflection height increased from ∼89 km from the first occurrence of tweek to about ∼93 km at the totality and then decreased to ∼88 km at the end of the eclipse, suggesting significant increase in tweek reflection height of about 5.5 km during the eclipse. The reflection heights at the time of totality during TSE are found to be less by 2–3 km as compared to the usual nighttime tweek reflection heights. This is due to partial nighttime condition created by TSE. A significant increase of 3 dB in the strength of the amplitude of VLF signal of 22.2 kHz transmitted from JJI-Japan is observed around the time of the total solar eclipse (TSE) as compared to a normal day. The modeled electron density height profile of the lower ionosphere depicts linear variation in the electron density with respect to solar radiation as observed by tweek analysis also. These low latitude ionospheric perturbations on the eclipse day are discussed and compared with other normal days.