Simultaneous bi-station measurements of mesospheric waves from Indian low latitudes

Coordinated mesospheric night-airglow measurements have been carried out from two stations, Gadanki (13.5° N, 79.2° E) and Allahabad (25.5° N, 81.9° E), India during April 2009 to study the common gravity wave features. With two nights of coordinated measurements we find some of the wave periodicities to be similar at the two locations. Simultaneous OH and O₂ intensity measurements over Gadanki reveal these features to be upward propagating gravity waves while the coordinated OH intensity measurements of similar waves from Allahabad show the large spatial extent of these waves.     

Gravity wave propagation studies using the Indian MST radar observations

The Indian MST radar facility at Gadanki (13.5°N, 79.2°E) has been utilised to study the propagation of gravity waves from the troposphere/lower stratosphere to the mesosphere and their interaction with the radar backscattered signal variations. The main objective is to correlate vertically propagating gravity waves derived from the tropospheric velocity fields with the dynamics of mesospheric scattering centres. The tropospheric wind velocities and signal strengths over the entire height range have been subjected to power spectral and wavelet analysis to determine the predominant wave periods/amplitudes and the coupling between the lower atmosphere and mesosphere. Results show that (a) the gravity waves are clearly detectable near tropopause heights, (b) while relatively higher period gravity waves (20–50 min) interact with mesospheric scattering centres, the lower period waves (<20 min) are absorbed in the troposphere itself, (c) the mesospheric scattering layers are affected by gravity waves of complementary periods.     

Wave perturbations in the MLT at high northern latitudes in winter,observed by two SATI instruments

To investigate the Mesosphere and Lower Thermosphere (MLT) region, several ground-based instruments called SATI (Spectral Airglow Temperature Imager) were designed and built to measure airglow emission and temperature in the upper mesosphere. One SATI instrument was installed at Resolute Bay (74.7°N, 94.9°W) and has monitored the polar MLT region since November, 2001. In October 2007 another SATI instrument was installed at Eureka (80.0°N, 86.3°W) at the Polar Environment Atmospheric Research Laboratory (PEARL) as part of the Canadian Network for the Detection of Atmospheric Change (CANDAC) project. SATI is a spatial scanning Fabry–Perot spectrometer measuring column emission rates for several rotational lines of OH and O₂ airglow at 87 and 94 km height. The rotational temperatures are inferred from the ratios of these lines. The measurements are divided into 12 sectors with an annular field of view. The phase differences between the sectors yield information on the horizontal atmospheric wave direction and wavelength. Horizontal perturbations of 2–8 h period have correlatively been observed and investigated at both locations. Short-periodic oscillations identified as gravity waves with periods between 2 and 8 h propagate in southward and eastward directions, but in opposite directions in some cases. The wave propagation characteristics are often different at the two locations; the relationship with the lower mean wind is considered.     

Occurrence features of simultaneous H+- and He+-band EMIC emissions in the outer radiation belt

As an important loss mechanism of radiation belt electrons, electromagnetic ion cyclotron (EMIC) waves show up as three distinct frequency bands below the hydrogen (H⁺), helium (He⁺), and oxygen (O⁺) ion gyrofrequencies. Compared to O⁺-band EMIC waves, H⁺- and He⁺-band emissions generally occur more frequently and result in more efficient scattering removal of <～5?MeV relativistic electrons. Therefore, knowledge about the occurrence of these two bands is important for understanding the evolution of the relativistic electron population. To evaluate the occurrence pattern and wave properties of H⁺- and He⁺-band EMIC waves when they occur concurrently, we investigate 64 events of multi-band EMIC emissions identified from high quality Van Allen Probes wave data. Our quantitative results demonstrate a strong occurrence dependence of the multi-band EMIC emissions on magnetic local time (MLT) and L-shell to mainly concentrate on the dayside region of L?=?～4–6. We also find that the average magnetic field amplitude of H⁺-band waves is larger than that of He⁺-band waves only when L?<?4.5 and AE¹?<?300?nT, and He⁺-band emissions are more intense under all other conditions. In contrast to 5 events that have average H⁺-band amplitude over 2 nT, 19 events exhibit >2 nT He⁺-band amplitude, indicating that the He⁺-band waves can be more easily amplified than the H⁺-band waves under the same circumstances. For simultaneous occurrences of the two EMIC wave bands, their frequencies vary with L-shell and geomagnetic activity: the peak wave frequency of H⁺-band emissions varies between 0.25 and 0.8 f_cp with the average between 0.25 and 0.6 f_cp, while that of He⁺-band emissions varies between 0.03 and 0.23 f_cp with the average between 0.05 and 0.15 f_cp. These newly observed occurrence features of simultaneous H⁺- and He⁺-band EMIC emissions provide improved information to quantify the overall contribution of multi-band EMIC waves to the loss processes of radiation belt electrons.     

Vertical wavenumber spectra of atmospheric ozone measured from ozonesonde observations

Vertical profiles of ozone have been measured at balloon altitudes. Our purpose is to examine the character of vertical wavenumber spectra of ozone fluctuations, to assess the possible roles of gravity wave field in ozone fluctuations, and to determine dominant vertical wavelengths of ozone spectra. Vertical wavenumber spectra of 12 ozone fluctuations obtained during June–August 2003 are presented. Results indicate that mean spectral slopes in the wavenumber range from 4.69 × 10⁻⁴ to 2.50 × 10⁻³ cyc/m are about −2.91 in the troposphere and −2.87 in the lower stratosphere, which is close to the slope of −3 predicted by current gravity wave saturation models. The consistency of the observed spectral slopes with the value of −3 predicted by current gravity wave saturation models suggests that the observed ozone fluctuations are due primarily to atmospheric gravity waves. At m = 1/(1000 m) the mean spectral amplitude is over 30 times larger in the lower stratosphere than in the troposphere. Mean vertical wavenumber spectra in area-preserving form reveal dominant vertical wavelengths of ∼2.6 km in the troposphere and ∼2.7 km in the lower stratosphere, which is consistent with the values varying between 1.5 and 3.0 km estimated from the velocity field and temperature field at these heights.     

Estimates of mesospheric gravity wave activity over convection from a global model

A new convective gravity wave source spectrum parameterization has been implemented in the Whole Atmosphere Community Climate Model version 2 (WACCM2). This parameterization specifies the momentum flux phase speed spectrum of gravity waves in the Tropics based on the properties of underlying convection; Hence, this parameterization provides realistic global estimates of gravity wave activity. In this paper, we show the estimated gravity wave phase speed spectra in the Tropics from a WACCM2 simulation, at the source level and at 85 km. Spatial distribution of gravity wave activity at 85 km is also presented. Subsequently, we discuss the factors that are primarily responsible for the estimated differences in gravity wave distribution across phase speeds with latitude and asymmetries in direction of gravity wave propagation in the mesosphere. We also examine which of the model assumptions can lead to uncertainties in our estimates of mesospheric gravity wave activity and we discuss how these assumptions provide challenges for comparison with observations of gravity waves in the mesosphere.     

Forward and reverse CIR shocks at 4–5 AU: Ulysses

In this article, we study fast shocks at CIR boundaries during an extended interval of 15 consecutive major high speed solar wind streams in 1992–1993. Ulysses was 4–5 AU from the sun. The Abraham-Schrauner shock normal method and the Rankine-Hugoniot relations were used to determine fast shock directions and speeds. Out of 33 potential CIR shocks, 14 were determined to be fast forward shocks (FSs) and 14 were fast reverse shocks (RSs). Of the remaining 5 events, 2 were forward waves and 3 were reverse waves. CIR edges at latitudes below ∼30^o were, for the most part, bounded by fast magnetosonic shocks. The forward shocks were generally quasi-perpendicular (average θ_nBo = 67^o). The reverse shocks were more oblique (average θ_nBo = 52^o), but they extended to all angles. Both FSs and RSs had magnetosonic Mach numbers ranging from 1 to 5 or 6. The average Mach numbers were 2.4 and 2.6 for FSs and RSs, respectively. The shock Mach numbers were noted to generally decrease with increasing latitude. The non-shock events or waves were noted to occur preferentially at high (∼−30° to −35°) heliolatitudes where stream-stream interactions were presumably weaker. These results are consistent with expectations, indicating the general accuracy of the Abraham-Schrauner technique.     

Thermospheric F-region travelling disturbances detected at low latitude by an OI 630 nm digital imager system

Gravity wave effects in the nocturnal thermospheric F-region domain are seldom detected in the intertropical region by optical (airglow) techniques, especially during geomagnetically quiet times, in part because the low inclination of the magnetic field, as opposed to the case of the mid-latitude region, does not favor significant vertical excursions of ionospheric plasma in response to meridional winds. Such difficulty in detecting gravity wave signatures in the F-region by means of optical techniques tends to increase in the absence of geomagnetic storms because of the lack of strong forcing mechanisms necessary to generate high intensity gravity waves. The purpose of this work is to show that during the quiet day of 9 August 1999, the Terminator may have been a source region of wave-like disturbances in the nocturnal F-region at the low-latitude station Cachoeira Paulista (22°41'S; 45°00W, dip 30°). A digital all-sky OI 630nm imager system located at that station has shown propagating wave-like spatial structures in the airglow intensity near the Terminator. This observation supports a previous study on the evidence of the presence of gravity waves during the post-sunset period at the same location by means of a scanning photometer system (1997, Sobral, J. Atmos. Terr. Phys. 59, 1611–1623). The absence of range-type spread-F as monitored by a local digisonde and the absence of radio wave scintillation as monitored by a local GPS receiver, excludes the hypothesis that the wave-like airglow structures are associated with the occurrence of the ionospheric plasma bubbles. Downwards motion of the iso-density real height contours at 22.0 ms⁻¹ and 33.1 ms⁻¹ are observed. The wave detection by the imager system is reported and discussed here.     

The response of sporadic E-layer to the total solar eclipse of July 22, 2009 over the equatorial ionization anomaly region of the Indian zone

The digital ionosonde located in Bhopal (23.2°N, 77.2°E), India has been used to investigate the responses of the Es layer in the equatorial ionization anomaly (EIA) crest to the total solar eclipse (TSE) of July 22, 2009. Results show the presence of intense Es layer during and after the eclipse period. The gravity waves induced by the solar eclipse propagated upward in the Es layer and produced the periodic disturbance. The results of the wavelet analysis display the presence of dominant oscillation of about 24–32, 16–20 and 8 min. The appearance of intense sporadic-E concomitantly with the signatures of gravity wave suggests that the wind shear introduced by the solar eclipse induced gravity wave might be the plausible mechanism behind the intensification of Es-layer ionization.     

Inversion of satellite gradiometry data using statistically modified integral formulas for local gravity field recovery

The satellite gravity gradiometric data can be used directly to recover the gravity anomaly at sea level using inversion of integral formulas. This approach suffers by the spatial truncation errors of the integrals, but these errors can be reduced by modifying the formulas. It allows us to consider smaller coverage of the satellite data over the region of recovery. In this study, we consider the second-order radial derivative (SORD) of disturbing potential (T_rr) and determine the gravity anomaly with a resolution of 1° × 1° at sea level by inverting the statistically modified version of SORD of extended Stokes’ formula. Also we investigate the effect of the spatial truncation error on the quality of inversion considering noise of T_rr. The numerical investigations show satisfactory results when the area of T_rr coverage is the same with that of the gravity anomaly and the integral formula is modified by the biased least-squares modification. The error of recovery will be about 6 mGal after removing the regularization bias in the presence of 1 mE noise in T_rr measured on the orbit.     

Gravitational wave detection on the Moon and the moons of Mars

The Moon and the moons of Mars should be extremely quiet seismically and could therefore become sensitive gravitational wave detectors, if instrumented properly. Highly sensitive displacement sensors could be deployed on these planetary bodies to monitor the motion induced by gravitational waves. A superconducting displacement sensor with a 10-kg test mass cooled to 2 K will have an intrinsic instrument noise of 10⁻¹⁶ m Hz^−1/2. These sensors could be tuned to the lowest two quadrupole modes of the body or operated as a wideband detector below its fundamental mode. An interesting frequency range is 0.1–1 Hz, which will be missed by both the ground detectors on the Earth and LISA and would be the best window for searching for stochastic background gravitational waves. Phobos and Deimos have their lowest quadrupole modes at 0.2–0.3 Hz and could offer a sensitivity h_min ? 10⁻²² Hz^−1/2 within their resonance peaks, which is within two orders of magnitude from the goal of the Big Bang Observer (BBO). The lunar and Martian moon detectors would detect many interesting foreground sources in a new frequency window and could serve as a valuable precursor for BBO.     

Aura/MLS与TIMED/SABER观测全球重力波特性

大气重力波是临近空间环境主要大气波动之一,对全球环流具有重要影响。卫星上搭载的临边探测器能够探测临近空间大气温度,可用于临近空间大气重力波研究。利用2012－2014年Aura的微波临边探测器（MLS）和TIMED的红外临边探测器（SABER）的探测数据,对20~50 km高度的大气重力波扰动分布特征开展了分析研究,两种观测重力波活动基本一致,重力波随季节、纬度及高度的变化显著。冬季半球高纬度重力波扰动较强,赤道和夏季半球近赤道地区上空也存在明显重力波活动区域,夏季半球高纬度重力波扰动最弱。重力波扰动强度随高度增加。TIMED/SABER重力波扰动强度数值比 Aura/MLS略强。      

AIRS观测资料研究全球平流层重力波特性

利用2012—2014年1月和7月AIRS（Atmospheric Infrared Sounder）第79通道的观测数据,分析了平流层重力波活动强弱的全球分布以及重力波发生频率的全球分布;分析了重力波活动随纬度和经度的变化特征,给出了重力波活动在全球范围内的热点区域及其活动强度;对比了白天与夜间的重力波活动强度及发生频率.研究表明重力波活动强度呈现出随纬度变化的特征,在低纬度地区（0°—30°）,冬季半球重力波活动强度低,夏季半球重力波活动高;在中高纬度地区,冬季半球重力波活动强度高而夏季半球重力波活动强度低.在1月,全球重力波活动有4个突出的热点区域,分别为50°N附近欧洲大陆与大西洋交接地带、北美洲与大西洋交接地,20°S附近南美洲与大西洋交接地区、非洲与印度洋交接地区.在7月,重力波活动突出的地方为巴塔哥尼亚至南极半岛地区,以及50°S和75°E附近的印度洋区域.重力波活动强度在夜间大于白天,但是夜间的强重力波活动区域小于白天.      

Registration of wave disturbances over Yakutia

The results of investigations of wave processes with periods 2 hours on their influence and on the night sky airglow intensity are given. The observations were carried out by multichannel spectrometer for three seasons of 1985–1988 at the optical testing ground Maimaga (γ = 63°N; λ = 129, 5°E). The synchronous detection of two and sometimes of three emissions of night sky airglow yielded the oppotunity to track a vertical travel of waves and to estimate their parameters. In most cases the waves propagate upward, i.e. the sources of waves were below mesosphere. The estimated vertical velocity change within 0,9-3,3 m/s and vertical wave length - within 18–85 km. A horizontal velocity varies from 83 to 330 m/s. The wave activity (the occurence frequency) and their amplitude in winter is higher than in spring. The estimated energies transfered by waves to the upper atmosphere are in winter 3.8·10⁻³ W/m² and in spring 2.7·10⁻³ W/m².     

Ion dynamics associated with Alfven wave in the near-Earth magnetotail: Two-dimensional global hybrid simulation

Ion dynamics in the near-Earth magnetotail region is examined during periods of fast Earthward flow with a two-dimensional (2-D) global-scale hybrid simulation. The simulation shows that shear Alfven waves are generated at x ∼ −10R_E, where the strong earthward flow is arrested by the dipole field, and propagate along field lines from the equator to both southern and northern polar ionosphere. Non-gyrotropic ion velocity distributions occur where the large-amplitude Alfven waves are dominant. The simulation indicates that the Alfven waves are generated by interaction of the fast earthward flow with the stationary near-Earth plasma. Beam ions are found to be pitch-angle scattered and trapped in the wave field, leading to the non-gyrotropic ion distributions in the high-latitude plasma sheet boundary. In addition, significant particle heating and acceleration are found to occur behind the dipolarization front due to the effect of wave turbulence.     

Study of waves in the magnetotail region with cluster and DSP

The study of the neutral sheet is of fundamental importance in understanding the dynamics of the Earth’s magnetosphere. From the earliest observation of the magnetotail, it has been found that the neutral sheet frequently appears to be in motion due to changing solar wind conditions and geomagnetic activity. Multiple crossings of the neutral sheet by spacecraft have been attributed to a flapping motion of the neutral sheet in the north–south direction, a wavy profile either along the magnetotail or the dawn–dusk direction. Cluster observations have revealed that the flapping motions of the Earth’s magnetotail are of internal origin and that kink-like waves are emitted from the central part of the tail and propagate toward the tail flanks. This flapping motion is shown here to propagate at an angle of ∼45° with x_GSM. A possible assumption that the flapping could be created by a wake travelling away from a fast flow in the current sheet is rejected. Other waves in the magnetotail are found in the ULF range. One conjunction event between Cluster and DoubleStar TC1 is presented where all spacecraft show ULF wave activity at a period of approximately 5 min during fast Earthward flow. These waves are shown to be Kelvin–Helmholtz waves on the boundaries of the flow channel. Calculations show that the conversion of flow energy into magnetic energy through the Kelvin–Helmholtz instability can contribute to a significant part of flow breaking between Cluster and DoubleStar TC1.     

3D Imaging of the OH mesospheric emissive layer

A new and original stereo imaging method is introduced to measure the altitude of the OH nightglow layer and provide a 3D perspective map of the altitude of the layer centroid. Near-IR photographs of the OH layer are taken at two sites separated by a 645 km distance. Each photograph is processed in order to provide a satellite view of the layer. When superposed, the two views present a common diamond-shaped area. Pairs of matched points that correspond to a physical emissive point in the common area are identified in calculating a normalized cross-correlation coefficient (NCC). This method is suitable for obtaining 3D representations in the case of low-contrast objects. An observational campaign was conducted in July 2006 in Peru. The images were taken simultaneously at Cerro Cosmos (12°09′08.2″ S, 75°33′49.3″ W, altitude 4630 m) close to Huancayo and Cerro Verde Tellolo (16°33′17.6″ S, 71°39′59.4″ W, altitude 2272 m) close to Arequipa. 3D maps of the layer surface were retrieved and compared with pseudo-relief intensity maps of the same region. The mean altitude of the emission barycenter is located at 86.3 km on July 26. Comparable relief wavy features appear in the 3D and intensity maps. It is shown that the vertical amplitude of the wave system varies as exp (Δz/2H) within the altitude range Δz = 83.5–88.0 km, H being the scale height. The oscillatory kinetic energy at the altitude of the OH layer is comprised between 3 × 10⁻⁴ and 5.4 × 10⁻⁴ J/m³, which is 2–3 times smaller than the values derived from partial radio wave at 52°N latitude.     

Solar cycle effect on temporal and spatial variation of topside ion density measured by SROSS C2 and ROCSAT 1 over the Indian longitude sector

We investigated the diurnal, seasonal and latitudinal variations of ion density Ni over the Indian low and equatorial topside ionosphere within 17.5°S to 17.5°N magnetic latitudes by combining the data from SROSS C2 and ROCSAT 1 for the 9 year period from 1995 to 2003 during solar cycle 23. The diurnal maximum density is found in the local noon or in the afternoon hours and the minimum occurs in the pre sunrise hours. The density is higher during the equinoxes as compared to that in the June and December solstice. The local time spread of the daytime maximum ion density increases with increase in solar activity. A north south asymmetry with higher ion density over northern hemisphere in the June solstice and over southern hemisphere in December solstice has been observed in moderate and high solar activity years. The crest to crest distance increases with increase in solar flux. Ion density bears a nonlinear relationship with F_10.7 cm solar flux and EUV flux in general. The density increases linearly with solar flux up to ∼150 sfu (1 sfu = 10⁻²²Wm⁻²Hz⁻¹) and EUV flux up to ∼50 units (10⁹ photons cm⁻² s⁻¹). But beyond this the density saturates. Inverse saturation and linear relationship have been observed in some season or latitude also. Inter-comparison of the three solar activity indices F_10.7 cm flux, EUV flux and F_10.7P (= (F_10.7 + F_10.7A)/2, where F_10.7A is the 81 day running average value of F_10.7) shows that the ion density correlates better with F_10.7P and F_10.7 cm fluxes. The annual average daytime total ion density from 1995 to 2003 follows a hysteresis loop as the solar cycle reverses. The ion density at 500 km over the Indian longitude sector as obtained by the international reference ionosphere is in general lower than the measured densities during moderate and high solar activity years. In low solar activity years the model densities are equal or higher than measured densities. The IRI EIA peaks are symmetric (±10°) in equinox while densities are higher at 10°N in June solstice and at 10°S in the December solstice. The model density follows F_10.7 linearly up to about F_10.7 > ∼150 sfu and then saturates.     

Evidence for a 2200 km extended wave system at the mesopause level

A panoramic view of the nightglow atmospheric emission in the 780–1000 nm spectral range is constructed using CCD images taken at the Pic de Châteaurenard (Altitude 2989 m, Hautes-Alpes) on July 14–15, 1999. A set of 28 images each having a 36° × 36° field of view is assembled to form a panorama covering 360° in azimuth and extending from the horizon to the zenith. Each photograph is processed in order to invert the perpective effect assuming that the emission comes from a thin layer located at the altitude of 85 km. The effect of refraction is calculated and taken into account. The stars are removed using a numerical filter. The inverted panorama appears as a disk having a radius equal to 1100 km. It is comparable to a satellite view of the emissive layer. A wave system extends in the W-NW to E-SE direction over more than 2200 km. A second set of 30 successive images of the same field of view taken on May 18–19, 1998 is used to determine the wave parameters. The main horizontal wavelength is equal to 42 km and the horizontal phase velocity has a value of 40 ± 2 m.s⁻¹. The images show that the atmospheric OH emission is a tracer of the dynamics of the atmosphere at the level where the excited OH radicals are produced. The OH* radical population depends upon its quenching by O, O₂ and N₂. As a result, the emission intensity is a function of the air temperature and density which are subject to variations due to gravity and windshear waves and other dynamic processes such as tides and turbulence.