Assessment of the consequences of the Fengyun-1C breakup in low Earth orbit

On 11 January 2007, the People’s Republic of China conducted a successful anti-satellite test against one of their defunct polar-orbiting weather satellites. The target satellite, called Fengyun-1C, had a mass of 880 kg and was orbiting at an altitude of about 863 km when the collision occurred. Struck by a direct-ascent interceptor at a speed of 9.36 km/s, the satellite disintegrated, spreading the cataloged fragments between 200 and 4000 km, with the highest concentration near the breakup height. By the end of April 2008, 2377 pieces of debris, including the original payload remnant, had officially been cataloged by the US Space Surveillance Network. Of these, nearly 1% had reentered the Earth’s atmosphere. This deliberate act is the largest debris-generating event on record, and its consequences will adversely affect circumterrestrial space for many years.     

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.     

Mesospheric doppler wind measurements from Aura Microwave Limb Sounder (MLS)

This paper describes a microwave limb technique for measuring Doppler wind in the Earth’s mesosphere. The research algorithm has been applied to Aura Microwave Limb Sounder (MLS) 118.75 GHz measurements where the O₂ Zeeman lines are resolved by a digital autocorrelation spectrometer. A precision of ∼17 m/s for the line-of-sight (LOS) wind is achieved at 80–92 km, which corresponds to radiometric noise during 1/6 s integration time. The LOS winds from Aura MLS are mostly in the meridional direction at low- and mid-latitudes with vertical resolution of ∼8 km. This microwave Doppler technique has potential to obtain useful winds down to ∼40 km of the Earth’s atmosphere if measurements from other MLS frequencies (near H₂O, O₃, and CO lines) are used. Initial analyses show that the MLS winds from the 118.75 GHz measurements agree well with the TIDI (Thermosphere Ionosphere Mesosphere Energetics and Dynamics Doppler Interferometer) winds for the perturbations induced by a strong quasi 2-day wave (QTDW) in January 2005. Time series of MLS winds reveal many interesting climatological and planetary wave features, including the diurnal, semidiurnal tides, and the QTDW. Interactions between the tides and the QTDW are clearly evident, indicating possible large tidal structural changes after the QTDW events dissipate.     

Consequences of the application of the streamer fluid model to the study of the sprite inception mechanism

We present the results of a streamer-fluid model used to investigate the electrodynamical coupling between the troposphere and upper atmosphere due to the penetration of lightning electric fields into the mesosphere and the lower ionosphere, generating sprites. The model solves the continuity equation for electrons and ions coupled to Poisson equation. The dominant physical response of the atmosphere is the formation of a screening-ionization wave. The wave shields the atmosphere above it from the action of the lightning field and, together with the conductivity reduction below it due to attachment, the wave amplifies the total field below it, allowing for the penetration of intense electric fields in the mesosphere as it propagates downwards into regions of higher density that compress the wave. This is the key physical mechanism for sprite inception. We evaluated the effects of the thundercloud charge geometry, lightning current waveshape, atmospheric conductivity, via different electron density profiles, and the effect of ionization, attachment and electron mobility coefficients in the electrical breakdown process, related to halo production, and sprite streamer initiation. The results showed that electrons with higher mobility are more efficient in shielding the lightning electric field before breakdown, causing delay, and they contribute to the formation of the streamer seed after breakdown, anticipating the sprite streamer inception. Similarly, a higher effective ionization rate, produced by modifications in the attachment and ionization coefficients, anticipates sprite inception. The simulations with 6 different electron density profiles, and therefore conductivities, spanning 4 orders of magnitude, showed that the altitude of breakdown and sprite initiation, as well as their time delays from the lightning discharge are directly related to atmospheric conductivity: higher conductivities produce halo and sprite inception at lower altitudes with longer delays and may hinder sprite formation. We document that variations of 30 times in the lightning current leads to sprite initiation altitudes in the range 66.0–73.5 km, with delays between 1.550 and 34.500 ms, while variations of 4 orders of magnitude in the conductivity profile lead to initiation altitudes 61.0–70.6 km, with delays in the range 3.825–9.825 ms. Consequently, we suggest that lightning characteristics dominate over atmospheric parameters in determining sprites’ initiation altitude and delay. The simulation of a −CG, with a constant current of 30 kA, did not produce a sprite seed, confirming an asymmetry in the response of the atmosphere to positive and negative lightning. This is due to the free electron drift direction that is away from the screening ionization wave, preventing the formation of the streamer seed for the great majority of −CGs. The same does not apply to halos, which depend on the occurrence of breakdown and can be produced by discharges of both polarities.     

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

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

Thermal neutrons’ response to the GLEs

Thermal neutrons’ flux near the Earth’s crust is very sensitive regarding different processes and phenomena both in the near-Earth space and in the Earth’s crust by reason of the dual nature of the thermal neutron flux. Its first source is associated with high-energy particles of cosmic rays penetrating into the Earth’s atmosphere and interacting with its elements. The second source originates from the radioactive gases contained in the Earth’s crust. So the contributions of these two sources are specified by the Earth’s crust conditions and its movements, on one hand and variations of high-energy particles flux near the Earth.     

The ground GNSS tomography – unconstrained approach

The GNSS signal along its more than 20,000 km line of sight is bended, attenuated and delayed. These effects are results of the Earth’s atmosphere, and Sun radiation. Amongst all mentioned effects the one considered in this paper is a signal phase delay in troposphere.     

Gravity waves generated by thunderstorms: South-East Asia tropical region study

Gravity waves are recognized as an integral part of earth’s atmosphere which are mainly responsible for energy and momentum distribution among different layers and regions in the atmosphere. Various sources present in land, ocean, and atmosphere such as mountains, convection, jets and fronts etc. are responsible for gravity waves generation. Thunderstorms (deep convection) are one of the major sources of gravity waves in the tropical region, capable of generating waves with a wide range of frequencies and scales and significantly affecting the existing waves. Previous numerical studies have characterised the wave properties that are generated from thunderstorms, but there are no statistically quantified studies. In this paper, we have modelled the relationship between the latent heat generated inside a thunderstorm and the gravity wave properties at the geo-collocated points. Gravity waves are identified over Singapore radiosonde station (with data available until 30?km altitude with 12?h temporal resolution) in the stratosphere using wavelet studies. Based on the GROGRAT ray tracing methods to identify the thunderstorm locations, and RAMS cloud-resolving models simulations to obtain the latent heating of the thunderstorm, a regression analysis is performed using 200 cases of gravity waves. Furthermore, cloud-top momentum flux analysis is performed for various cases latent heat. This study is expected to provide more quantified and concrete information on the coupling between the thunderstorm and gravity wave which includes the variance in these relationships due to wave frequency spectrum and generation mechanisms.     

Lower ionosphere response to external forcing: A brief review

There are two ways of external forcing of the lower ionosphere, the region below an altitude of about 100 km: (1) From above, which is directly or indirectly of solar origin. (2) From below, which is directly or indirectly of atmospheric origin. The external forcing of solar origin consists of two general factors – solar ionizing radiation variability and space weather. The solar ionization variability consist mainly from the 11-year solar cycle, the 27-day solar rotation and solar flares, strong flares being very important phenomenon in the daytime lower ionosphere due to the enormous increase of the solar X-ray flux resulting in temporal terminating of MF and partly LF and HF radio wave propagation due to heavy absorption of radio waves. Monitoring of the sudden ionospheric disturbances (SIDs – effects of solar flares in the lower ionosphere) served in the past as an important tool of monitoring the solar activity and its impacts on the ionosphere. Space weather effects on the lower ionosphere consist of many different but often inter-related phenomena, which govern the lower ionosphere variability at high latitudes, particularly at night. The most important space weather phenomenon for the lower ionosphere is strong geomagnetic storms, which affect substantially both the high- and mid-latitude lower ionosphere. As for forcing from below, it is caused mainly by waves in the neutral atmosphere, i.e. planetary, tidal, gravity and infrasonic waves. The most important and most studied waves are planetary and gravity waves. Another channel of the troposphere coupling to the lower ionosphere is through lightning-related processes leading to sprites, blue jets etc. and their ionospheric counterparts. These phenomena occur on very short time scales. The external forcing of the lower ionosphere has observationally been studied using predominantly ground-based methods exploiting in various ways the radio wave propagation, and by sporadic rocket soundings. All the above phenomena are briefly mentioned and some of them are treated in more detail.     

Impact of the time span selected to calibrate the ballistic parameter on spacecraft re-entry predictions

The trajectory modeling of satellites that are re-entering the Earth’s atmosphere, as a result of natural orbital decay, has always been a challenging task. Residual lifetime estimations and re-entry predictions are affected by substantial uncertainties, associated with atmospheric density models, with the forecasts of the relevant solar and geomagnetic activity indices and with tracking data, which for uncontrolled re-entries are usually sparse and not particularly accurate. Furthermore, modeling the aerodynamic forces that act on low altitude satellites is a formidable task, especially for objects of a complex shape and unknown attitude evolution.     

Long-term (50 years) measurements of cosmic ray fluxes in the atmosphere

Since the middle of 1957 till present time the group of researchers of P.N. Lebedev Physical Institute of the Russian Academy of Sciences has carried out the regular balloon borne measurements of charged particle fluxes in the atmosphere. The measurements are performed at polar (northern and southern) and middle latitudes and cover the interval of heights from the ground level up to 30–35 km. Standard detectors of particles (gas-discharged counters) have been used. More than 80,000 measurements of cosmic ray fluxes in the atmosphere have been performed to the present time. In the data analysis the geomagnetic field and the Earth’s atmosphere are used as cosmic ray spectrometers.     

Atmosphere-ionosphere coupling from convectively generated gravity waves

Ionospheric variability impacts operational performances of a variety of technological systems, such as HF communication, Global Positioning System (GPS) navigation, and radar surveillance. The ionosphere is not only perturbed by geomagnetic inputs but is also influenced by atmospheric tides and other wave disturbances propagating from the troposphere to high altitudes. Atmospheric Gravity Waves (AGWs) excited by meteorological sources are one of the largest sources of mesoscale variability in the ionosphere. In this paper, Total Electron Content (TEC) data from networks of GPS receivers in the United States are analyzed to investigate AGWs in the ionosphere generated by convective thunderstorms. Two case studies of convectively generated gravity waves are presented. On April 4, 2014 two distinct large convective systems in Texas and Arkansas generated two sets of concentric AGWs that were observed in the ionosphere as Traveling Ionospheric Disturbances (TIDs). The period of the observed TIDs was 20.8 min, the horizontal wavelength was 182.4 km, and the horizontal phase speed was 146.4 m/s. The second case study shows TIDs generated from an extended squall line on December 23, 2015 stretching from the Gulf of Mexico to the Great Lakes in North America. Unlike the concentric wave features seen in the first case study, the extended squall line generated TIDs, which exhibited almost plane-parallel phase fronts. The TID period was 20.1 min, its horizontal wavelength was 209.6 km, and the horizontal phase speed was 180.1 m/s. The AGWs generated by both of these meteorological events have large vertical wavelength (>100 km), which are larger than the F2 layer thickness, thus allowing them to be discernible in the TEC dataset.     

Limb Viewing Hyper Spectral Imager (LiVHySI) for airglow measurements onboard YOUTHSAT-1

The Limb Viewing Hyper Spectral Imager (LiVHySI) is one of the Indian payloads onboard YOUTHSAT (inclination 98.73°, apogee 817 km) launched in April, 2011. The Hyper-spectral imager has been operated in Earth’s limb viewing mode to measure airglow emissions in the spectral range 550–900 nm, from terrestrial upper atmosphere (i.e. 80 km altitude and above) with a line-of-sight range of about 3200 km. The altitude coverage is about 500 km with command selectable lowest altitude. This imaging spectrometer employs a Linearly Variable Filter (LVF) to generate the spectrum and an Active Pixel Sensor (APS) area array of 256 × 512 pixels, placed in close proximity of the LVF as detector. The spectral sampling is done at 1.06 nm interval. The optics used is an eight element f/2 telecentric lens system with 80 mm effective focal length. The detector is aligned with respect to the LVF such that its 512 pixel dimension covers the spectral range. The radiometric sensitivity of the imager is about 20 Rayleigh at noise floor through the signal integration for 10 s at wavelength 630 nm. The imager is being operated during the eclipsed portion of satellite orbits. The integration in the time/spatial domain could be chosen depending upon the season, solar and geomagnetic activity and/or specific target area. This paper primarily aims at describing LiVHySI, its in-orbit operations, quality, potential of the data and its first observations. The images reveal the thermospheric airglow at 630 nm to be the most prominent. These first LiVHySI observations carried out on the night of 21st April, 2011 are presented here, while the variability exhibited by the thermospheric nightglow at O(¹D) 630 nm has been described in detail.     

Model for the VLF/LF radio signal anomalies formation associated with earthquakes

The influence of quasi-static electric field of seismic origin on the characteristics of the internal gravity waves (IGWs) in the Earth’s ionosphere is considered. The electric field in the ionosphere arises due to the injection of charged aerosols into the atmosphere, formation of an EMF in the near Earth atmosphere and perturbation of the conductive electric current in the global electric circuit. Amplification of the electric current in seismic zone is accompanied by the formation of perturbation of the lower ionosphere that affects the amplitude and phase of VLF/LF signals. The action of the electric field on the IGWs is connected with the appearance of the Ampere’s force in the ionosphere. In the spectral range of these waves the latter acts on the neutral component of the ionosphere plasma. As the result of this interaction the ionosphere starts to support the discrete spectrum of oscillations. Periods of their maximums increase as numbers of natural sequence. The existence of such peculiarities of the waves in the ionosphere is confirmed by observations.     

Influence of Earth’s shadowing effects on space debris stability

Solar radiation pressure affects the evolution of high area-to-mass geostationary space debris. In this work, we extend the stability study of Valk et al. (2009) by considering the influence of Earth’s shadows on the short- and long-term time evolutions of space debris. To assess the orbits stability, we use the Mean Exponential Growth factor of Nearby Orbits (MEGNO), which is an efficient numerical tool to distinguish between regular and chaotic behaviors. To reliably compute long-term space debris motion, we resort to the Global Symplectic Integrator (GSI) of Libert et al. (2011) which consists in the symplectic integration of both Hamiltonian equations of motion and variational equations. We show how to efficiently compute the MEGNO indicator in a complete symplectic framework, and we also discuss the choice of a symplectic integrator, since propagators adapted to the structure of the Hamiltonian equations of motion are not necessarily suited for the associated variational equations. The performances of our method are illustrated and validated through the study of the Arnold diffusion problem. We then analyze the effects of Earth’s shadows, using the adapted conical and cylindrical Earth’s shadowing models introduced by Hubaux et al. (2012) as the smooth shadow function deriving from these models can be easily included into the variational equations. Our stability study shows that Earth’s shadows greatly affect the global behaviour of space debris orbits by increasing the size of chaotic regions around the geostationary altitude. We also emphasize the differences in the results given by conical or cylindrical Earth’s shadowing models. Finally, such results are compared with a non-symplectic integration scheme.     

Study of the 2013 Lushan M7.0 earthquake coseismic ionospheric disturbances

On April 20, 2013, an earthquake of M7.0 occurred in Lushan, Sichuan province, China. This paper investigates the coseismic ionospheric anomalies using GPS (Global Positioning System) data from 23 reference stations in Sichuan province that are a part of the Crustal Movement Observation Network of China (CMONOC). The recorded results show that a clear ionospheric anomaly occurred within 15 min after the earthquake near the epicenter, and the occurrence time of the anomalies recorded by various stations is related to the distance from the epicenter. The maximum anomaly is 0.25 TECu, with a 2 min duration and the distance of the recording station to the epicenter is 83 km. Acoustic waves generated by the crustal vertical movement during the earthquake propagate up to the height of the ionosphere lead to the ionospheric anomaly, and the propagation speed of the acoustic wave is calculated as 0.72 ± 0.04 km/s based on the propagation time and propagation distance, consistent with the average speed of sound waves within a 0–450 km atmospheric height.     

Neutral upper atmospheres of Venus and Mars

Numerous measurements of the neutral upper atmosphere above 100 km have been made from spacecraft over Venus and over Mars. The Venus exospheric temperatures are unexpectedly low (less than 300°K near noon and less than 130°K near midnight). These very low temperatures may be partially caused by collisional excitation of CO₂ vibrational states by atomic oxygen and partially by eddy cooling. The Venus atmosphere is unexpectedly insensitive to solar EUV variability. On the other hand, the Martian dayside exospheric temperature varies from 150°K to 400°K over the 11-year solar cycle, where CO₂ 15-μm cooling may be less effective because of lower atomic oxygen mixing ratios. On Venus, temperature increases with altitude on the dayside (thermosphere), but decreases with altitude from 100 to 150 km on the nightside (cryosphere). However, dayside Martian temperatures near solar minimum for maximum planet-sun distance and low solar activity are essentially isothermal from 40 km to 200 km. During high solar activity, the thermospheric temperatures of Mars sharply increase. The Venus neutral upper atmosphere contains CO₂, O, CO, C, N₂, N, He, H, D and hot nonthermal H, O, C, and N, while the dayside Mars neutral upper atmosphere contains CO₂, O, O₂, CO, C, N₂, He, H, and Ar. There is evidence on Venus for inhibited day-to-night transport as well as superrotation of the upper atmosphere. Both atmospheres have substantial wave activity. Various theoretical models used to interpret the planetary atmospheric data are discussed.     

The effect of the longitudinal propagating waves on the electron acceleration in the reconnecting current sheet

Based upon the most efficient electron acceleration near the midplane of 3D non-neutral driven reconnecting current sheet (RCS) and the electrostatic wave excitation by the drift Maxwellian distribution of electrons in Vlasov simulation, we assume that the electrostatic waves mainly propagate opposite to the reconnecting electric field and investigated how these waves affect the electron acceleration. The main results are: (1) when the electron’s velocity equals to the phase speed of the waves, they will be trapped and have the different accelerating characteristics from the untrapped electrons through solving the momentum equations of electrons analytically; (2) the test particle simulations further prove that the number of the energetic electrons decreases with the increasing intensity of unstable waves, and the distribution of the energetic electrons takes on the double power-law.     

Plasma and electromagnetic wave simulations of meteors

Every day billions of meteoroids impact and disintegrate in the Earth’s atmosphere. Current estimates for this global meteor flux vary from 2000 to 200,000 tons per year, and estimates for the average velocity range between 10 km/s and 70 km/s. The basic properties of this global meteor flux, such as the average mass, velocity, and chemical composition remain poorly constrained. We believe much of the mystery surrounding the basic parameters of the interplanetary meteor flux exists for the following reason, the unknown sampling characteristics of different radar meteor observation techniques, which are used to derive or constrain most models. We believe this arises due to poorly understood radio scattering characteristics of the meteor plasma, especially in light of recent work showing that plasma turbulence and instability greatly influences meteor trail properties at every stage of evolution. We present our results on meteor plasmas simulations of head echoes using particle in cell (PIC) ions, which show that electric fields strongly influence early stage meteor plasma evolution, by accelerating ions away from the meteoroid body. We also present the results of finite difference time domain electromagnetic simulations (FDTD), which can calculate the radar cross section of the simulated meteor plasmas. These simulations have shown that the radar cross section depends in a complex manner on a number of parameters. These include the angle between radar and meteor entry, a large dependence on radar frequency, which shows that for a given meteor plasma size and density, the reflectivity as a function of probing radar frequency varies, but typically peaks below 100 MHz.