Global MGS TES data and Mars-GRAM validation

Mars Global Reference Atmospheric Model (Mars-GRAM 2001) is an engineering-level Mars atmosphere model widely used for many Mars mission applications. From 0 to 80 km, it is based on NASA Ames Mars General Circulation Model (MGCM), while above 80 km it is based on University of Michigan Mars Thermospheric General Circulation Model. Mars-GRAM 2001 and MGCM use surface topography from Mars Global Surveyor Mars Orbiter Laser Altimeter (MOLA). Validation studies are described comparing Mars-GRAM with a global summary data set of Mars Global Surveyor Thermal Emission Spectrometer (TES) data. TES averages and standard deviations were assembled from binned TES data which covered surface to ∼40 km, over more than a full Mars year (February 1999–June 2001, just before start of a Mars global dust storm). TES data were binned in 10° × 10° latitude–longitude bins (36 longitude bins, centered at 5°–355°, by 18 latitude bins, centered at −85° to +85°), and 12 seasonal bins (based on 30° increments of L_s angle). Bin averages and standard deviations were assembled at 23 data levels (temperature at 21 pressure levels, plus surface temperature and surface pressure). Two time-of-day bins were used: local time near 2 or 14 h. Two dust optical depth bins were used: infrared optical depth, either less than or greater than 0.25 (which corresponds to visible optical depth less than or greater than about 0.5). For interests in aerocapture and precision entry and landing, comparisons focused on atmospheric density. TES densities versus height were computed from TES temperature versus pressure, using assumptions of perfect gas law and hydrostatics. Mars-GRAM validation studies used density ratio (TES/Mars-GRAM) evaluated at data bin center points in space and time. Observed average TES/Mars-GRAM density ratios were generally 1 ± 0.05, except at high altitudes (15–30 km, depending on season) and high latitudes (>45°N), or at most altitudes in the southern hemisphere at L_s ∼ 90° and 180°. Compared to TES averages for a given latitude and season, TES data had average density standard deviation about the mean of ∼2.5% for all data, or ∼1–4%, depending on time of day and dust optical depth. Average standard deviation of TES/Mars-GRAM density ratio was 8.9% for local time 2 h and 7.1% for local time 14 h. Thus standard deviation of observed TES/Mars-GRAM density ratio, evaluated at matching positions and times, is about three times the standard deviation of TES data about the TES mean value at a given position and season.     

Latitudinal features of Total Electron Content over the African and European longitude sector following the St. Patrick’s day storm of 2015

GNSS TEC values have been obtained from 18 stations distributed from the magnetic equator to nearly 80°N magnetic dip in the African and west-European longitude sector corresponding to the March 17–18, 2015 geomagnetic storm. Significantly depleted ionosphere have been observed at stations north of 50°N geographic on March 18, 2015 following the above storm over a longitude swath 11.9°–21°E covering the Eastern Africa and Western European longitude sector. High ROTI values were noted on March 17^th at locations around 80°N magnetic dip. Two prominent peaks in PCN were noted around 09:00 UT and 14:00 UT on March 17, 2015 and around 15:00 UT on March 18, 2015. Daytime thermospheric (O/N₂) ratio was markedly less on March 18th at latitudes above 60°N geographic which is suggested to be the major driver behind depleted high latitude ionosphere during the recovery phase of the storm on March 18, 2015.     

The Venus ionosphere

Physical properties of the Venus ionosphere obtained by experiments on the US Pioneer Venus and the Soviet Venera missions are presented in the form of models suitable for inclusion in the Venus International Reference Atmosphere. The models comprise electron density (from 120 km), electron and ion temperatures, and relative ion abundance in the altitude range from 150 km to 1000 km for solar zenith angles from 0° to 180°. In addition, information on ion transport velocities, ionopause altitudes, and magnetic field characteristics of the Venus ionosphere, are presented in tabular or graphical form. Also discussed is the solar control of the physical properties of the Venus ionosphere.     

Observation of aphelion cloud belt over Martian tropics,its evolution,and associated dust distribution from MCS data

The present study uses five Martian years of observations from Mars Climate Sounder onboard Mars Reconnaissance Orbiter for investigating the Aphelion Cloud Belt (ACB) over the tropics. Analysis of zonal mean water ice column opacity suggests that the spatial extension of the ACB is mainly confined over the tropics and mid-latitudes (-20 – 40°N) during L_S ~ 45 – 135° (L_S = 0° signifies northern spring equinox). The ACB is seen primarily in the nighttime only due to the truncation of the daytime profile observations at significantly higher altitudes (at ~30 km). Zonal mean ice extinction profiles show ACB’s altitudinal range within ~10 – 40 km, and the existence of a thin cloud band in the absence of a thick ACB during aphelion season. Three phases of the ACB could be identified as the formation phase during L_S = 45 – 75° (phase 1), the peak phase during L_S = 76 – 105° (phase 2), and the decaying phase during L_S = 106 – 135° (phase 3). Observation of the cloud latitude belt shows a northward movement starting from phase 2, prominent over regions nearby Lunae Planum and Xanthe Terra. During this phase, the top level of thick clouds within the ACB decreases to ~20 km in the southern hemisphere, while it increases a little over the northern hemisphere (NH). The decreasing tendency continues in phase 3 over the entire region ?10 – 10°N, and the thick cloud base moves higher over the NH, though the vertical depth of it becomes narrower than phase 2. Temperature profiles do not show any noticeable influence on the northward evolution of the ACB. However, the study at a regional level indicates a possible association of upper tropospheric dustiness with the ACB’s evolution. The mechanism is evident in the correlation analysis mostly at an altitude range of ~18 – 35 km. The migrating semidiurnal tide (SMD) as a proxy of dust or water ice forcing, and the calculated upper tropospheric dust radiative heating, shows an apparent northward movement of their peak amplitude within the three phases of the ACB. This match between the spatiotemporal variations of the SMD and the water ice was not observed previously. However, the correlating behavior seems to be prominent in the areas nearby Lunae Planum and Xanthe Terra and the upper-tropospheric region of the atmosphere.     

火星表面风沙物理过程及沙尘暴的起动风速

通过建立风吹动地表沙粒运动的模型并根据大气湍流边界层风速廓线规律,计算了火星上沙尘暴的起动风速和沙尘暴发生时空中悬浮沙尘粒径的大小,并对沙粒从地面跳起进入气流的方式进行了验证.发现当大气为中性层结时,火星沙尘暴的起动需要离地面2m高处的风速达到28.7m·s-1.在起动临界风速下,地表沙粒需要滚动一个粒径的距离才能跳起,沙尘暴发生后,火星大气中悬浮沙尘的粒径小于30μm.      

由电离层电位所建立的大气电状态

本文在已知电离层电位分布下, 解析地计算了大气电位, 电场和电流强度的全球分布.结果表明, 在大气导电率随高度呈指数增加的情况下, 100km高度上的电离层电位, 几乎无衰减地扫到25km以下.大气电场较强的区域主要在20km以下的低层大气区, 其垂直分量比水平分量大4个数量级.而中高层大气电场较弱, 且两分量量级相当.本文还提出了一种考虑地面形状对大气电场影响的解析方法.      

The lower ionosphere response to its disturbances by powerful radio waves

The paper presents data from some campaigns at Sura heating facility in 2011–1016. The experiments on probing of the artificial disturbed region of the lower ionosphere were carried out at two observation sites. One of them was located near Vasil’sursk 1 km from Sura facility (56.1°N; 46.1°E) and the other site was located at the Observatory (55.85°N; 48.8°E) of Kazan State University, 170 km to the East. Investigation of the features of the disturbed region of the lower ionosphere based on its diagnostics by the methods of the vertical sounding and oblique backscattering is the main goal of this paper. Ionosphere disturbance was fulfilled by the effect of the powerful radio wave of the ordinary or extraordinary polarization emitted by transmitters of the Sura facility with effective radiated power ERP = 50–120 MW at the frequency of 4.3, 4.7 and 5.6 MHz. Pumping waves were emitted with period from 30 s to 15 min. The disturbed region of the ionosphere in Vasil’sursk was probed by the vertical sounding technique using the partial reflexion radar at the frequency of 2.95 and 4.7 MHz. For the oblique sounding of the disturbed region the modified ionosonde Cyclon-M, operating at ten frequencies from 2.01 to 6.51 MHz was used at the Observatory site. On many heating sessions simultaneous variations of the probing partial reflection signals in Vasil’sursk and backscattered signals in Observatory were observed at the height at 40–100 km below the reflection height of the pumping wave. These observations were correlated with the pumping periods of the Sura facility. Possible mechanisms of the appearance of the disturbance in the lower ionosphere and its effect on the probing radio waves are discussed.     

Effects of CME and CIR induced geomagnetic storms on low-latitude ionization over Indian longitudes in terms of neutral dynamics

This paper presents the response of the ionosphere during the intense geomagnetic storms of October 12–20, 2016 and May 26–31, 2017 which occurred during the declining phase of the solar cycle 24. Total Electron Content (TEC) from GPS measured at Indore, Calcutta and Siliguri having geomagnetic dips varying from 32.23°N, 32°N and 39.49°N respectively and at the International GNSS Service (IGS) stations at Lucknow (beyond anomaly crest), Hyderabad (between geomagnetic equator and northern crest of EIA) and Bangalore (near magnetic equator) in the Indian longitude zone have been used for the storms. Prominent peaks in diurnal maximum in excess of 20–45 TECU over the quiet time values were observed during the October 2016 storm at Lucknow, Indore, Hyderabad, Bangalore and 10–20 TECU for the May 2017 storm at Siliguri, Indore, Calcutta and Hyderabad. The GUVI images onboard TIMED spacecraft that measures the thermospheric O/N₂ ratio, showed high values (O/N₂ ratio of about 0.7) on October 16 when positive storm effects were observed compared to the other days during the storm period. The observed features have been explained in terms of the O/N₂ ratio increase in the equatorial thermosphere, CIR-induced High Speed Solar Wind (HSSW) event for the October 2016 storm. The TEC enhancement has also been explained in terms of the Auroral Electrojet (AE), neutral wind values obtained from the Horizontal Wind Model (HWM14) and equatorial electrojet strength from magnetometer data for both October 2016 and May 2017 storms. These results are one of the first to be reported from the Indian longitude sector on influence of CME- and CIR-driven geomagnetic storms on TEC during the declining phase of solar cycle 24.     

A three-dimensional modelling study of the processes leading to mid latitude nitric oxide increases in the lower thermosphere following periods of high geomagnetic activity

The processes leading to enhancements in mid latitude nitric oxide (NO) densities following geomagnetic storms have been investigated using the University College London (UCL) Coupled Middle Atmosphere and Thermosphere (CMAT) general circulation model. A comparison of calculated storm time and quiet time NO densities at 110 km altitude reveals the presence of aurorally produced NO at both high and mid latitudes for several days after subsidence of activity. At 150 km, the NO enhancements are shorter lived and remain for up to approximately 2 days after the storm. By separating the contribution of chemical production and loss, horizontal and vertical advection, and molecular and eddy diffusion in the calculation of NO densities, we show that at 150 km altitude, horizontal transport must be taken into consideration if post-storm mid latitude enhancements are to be reproduced. Chemical production of NO at high latitudes continues for up to 2 days after subsidence of a storm at altitudes of around 150 km. We show that equatorward winds at this altitude are sufficiently strong to transport the aurorally produced NO to mid latitudes. Vertical diffusion transports NO from altitudes of 150 km and above, to lower altitudes where it is longer lived. At 110 km altitude, chemical, diffusive and advective terms must all be included in the calculation of NO density in order to simulate realistic mid latitude enhancements. We propose that it is the combined effects of increased chemical production, downward diffusion from altitudes of 150 km and above, and transport by winds that lead to increases in mid latitude NO density at altitudes of around 110 km. This is the first detailed study of the causes of post-storm mid latitude NO enhancements to use a three-dimensional general circulation model.     

Impact of geomagnetic storm of September 7–8, 2017 on ionosphere and HF propagation: A multi-instrument study

The present study reveals the features of ionospheric parameters variations during the geomagnetic storm of September 7–8, 2017. In particular, parameters of vertical (foF2, foEs) and oblique ionospheric sounding (MOF, modes), absorption level, Total Electron Content (TEC) and particle fluxes at high altitudes were under analysis. The storm was characterized by two Dst-index mimima and can be considered as a sequence of two storms: first - with Dstmin?=??142?nT at 02 UT on September 8th and second - with Dstmin?=??122?nT and at 15 UT on September 8th. It was found that these two storms had different impacts on the ionosphere and HF propagation at mid- and high-latitudes of Northern Hemisphere. The signals of vertical and oblique ionospheric sounding were present in all ionograms before the first storm. Further, at the maximum of the first storm these signals were totally absorbed. Then, before the second storm and during its maximum the signals were detected again in the ionograms due to the low absorption. GOES satellite data showed the significant burst of electrons and protons only during the first storm and small particle fluxes - during the second storm. This feature was also confirmed with GPS data: TEC increased during the first storm and decreased during the second storm.     

Effects of geomagnetic storms on the bottomside ionospheric F region

The geomagnetic storm is a complex process of solar wind/magnetospheric origin. The variability of the ionospheric parameters increases substantially during geomagnetic storms initiated by solar disturbances. Various features of geomagnetic storm act at various altitudes in the ionosphere and neutral atmosphere. The paper deals with variability of the electron density of the ionospheric bottomside F region at every 10 km of altitude during intense geomagnetic storms with attention paid mainly to the distribution of the F1 region daytime ionisation. We have analysed all available electron density profiles from some European middle latitude stations (Chilton, Pruhonice, Ebro, Arenosillo, Athens) for 36 events that occurred in different seasons and under different levels of solar activity (1995–2003). Selected events consist of both depletion and increase of the F2 region electron density. For European higher middle and middle latitude the F1 region response to geomagnetic storm was found to be negative (decrease of electron density) independent on the storm effect on the F2 region. For lower middle latitude the F1 response is weaker and less regular. Results of the analysis also show that the maximum of the storm effect may sometimes occur below the height of the maximum of electron density (NmF2).     

中低纬电离层不同季节地磁暴响应的事例分析

利用中国中低纬台站漠河（53.5°N,122.3°E）、北京（40.3°N,116.2°E）、武汉（30.5°N,114.2°E）和三亚（18.3°N,109.6°E）的电离层观测数据,对比分析了4个台站电离层参数在2015年不同季节4个地磁扰动事件期间的变化特征.结果表明,4个磁暴事件期间电离层的响应特征并不完全一致,有着明显的季节特征,春季、夏季和秋季电离层以负相扰动为主,冬季以正相扰动为主.分析发现,中性成分O/N₂的降低与电离层负相扰动有关,但三亚地区的负相扰动还与扰动发电机电场相关.正相扰动的机制在不同事件中并不相同,穿透电场可能是引起春季磁暴事件期间电离层短时正暴效应的原因,而冬季长时间的正暴效应则是扰动电场和中性风共同作用的结果.      

Recent results on the Venus atmosphere from pioneer Venus radio occultations

Radio occultation measurements of the temperature structure of the Venus atmosphere have been obtained during seven occultation “seasons” extending from December 1978 to December 1983. Approximately 123 vertical profiles of temperature from about 40 km to about 85 km altitudes have been derived. Since these measurements cover latitudes from both poles to the equator, they have shown the latitudinal dependence of thermal structure. There is a smooth transition from the troposphere to the mesosphere at latitudes below about 45°, with the tropopause at about 56 km. The troposphere then rises to about 62 km in the “collar cloud” region between about 60° and 80° latitude, where a strong temperature inversion (up to 30 K) is present. In the polar areas, 80°–90°, the mesosphere becomes isothermal and there is no inversion. This latitudinal behavior is related to the persistent circulation pattern, in which a predominantly zonal retrograde motion at latitudes below 45° gradually changes to a circumpolar vortex at the “collar cloud” latitudes. Indeed, the radio occultation data have been used in a cyclostrophic balance model to derive zonal winds in the Venus atmosphere, which showed a mid-latitude (50°–55°) jet with a speed of about 120–140 ms^?1 at about 70 km altitude /1,2/. The observations obtained in 1983 and 1984 have shown that above the tropopause there is considerable temporal variability in the detailed thermal structure, suggesting that the persistent circulation pattern is subject to weather-like variability.     

Statistical ionospheric E layer properties measured with the Cyprus Digisonde and comparisons with IRI predictions

This paper reports the diurnal, seasonal, and long term variability of the E layer critical frequency (foE) and peak height (hmE) derived from Digisonde measurements from 2009 to 2016 at the low-middle latitude European station of Nicosia, Cyprus (geographical coordinates: 35°N, 33°E, geomagnetic lat. 29.38°N, I = 51.7°). Manually scaled monthly median values of foE and hmE are compared with IRI-2012 predictions with a view to assess the predictability of IRI. Results show that in general, IRI slightly overestimates foE values both at low and high solar activity. At low solar activity, overestimations are mostly limited to 0.25?MHz (equivalent electron density, 0.775?×?10³?el/m^?3) but can go as high as 0.5?MHz (equivalent electron density, 3.1?×?10³?el/m^?3, during noon) around equinox. In some months, underestimations, though sporadic in nature, up to 0.25?MHz are noted (mostly during sunrise and sunset). At high solar activity, a similar pattern of over-/underestimation is evident. During the entire period of study, over-/under estimations are mostly limited to 0.25?MHz. In very few cases, these exceed 0.25?MHz but are limited to 0.5?MHz. Analysis of hmE reveals that: (1) hmE remains almost constant during ±2 to ±4?h around local noon, (2) hmE values are higher in winter than in spring, summer and autumn, (3) there are two maxima near sunrise and sunset with a noontime minimum in between. During the entire period of study, significant differences between observed hmE and the IRI predictions have been noted. IRI fails to predict hmE and outputs a constant value of 110?km, which is higher than most of the observed values. Over- and under estimations range from 3 to 13?km and from 0 to 3?km respectively.     

Ionospheric total electron content response to annular solar eclipse on June 21, 2020

The effects of physical events on the ionosphere structure is an important field of study, especially for navigation and radio communication. The paper presents the spatio-temporal ionospheric TEC response to the recent annular solar eclipse on June 21, 2020, which spans across two continents, Africa and Asia, and 14 countries. This eclipse took place on the same day as the June Solstice. The Global Navigation Satellite System (GNSS) based TEC data of the Global Ionosphere Maps (GIMs), 9 International GNSS Service (IGS) stations and FORMOSAT-7/COSMIC-2 (F7/C2) were utilized to analyze TEC response during the eclipse. The phases of the TEC time series were determined by taking the difference of the observed TEC values on eclipse day from the previous 5-day median TEC values. The results showed clear depletions in the TEC time series on June 21. These decreases were between 1 and 9 TECU (15–60%) depending on the location of IGS stations. The depletions are relatively higher at the stations close to the path of annular eclipse than those farther away. Furthermore, a reduction of about ?10 TECU in the form of an equatorial plasma bubble (EPB) was observed in GIMs at ~20° away from the equator towards northpole, between 08:00–11:00 UT where its maximum phase is located in southeast Japan. Additionally, an overall depletion of ~10% was observed in F7/C2 derived TEC at an altitude of 240 km (hmF2) in all regions affected by the solar eclipse, whereas, significant TEC fluctuations between the altitudes of 100 km ? 140 km were analyzed using the Savitzky-Golay smoothing filter. To prove TEC depletions are not caused by space weather, the variation of the sunspot number (SSN), solar wind (V_SW), disturbance storm-time (Dst), and Kp indices were investigated from 16th to 22nd June. The quiet space weather before and during the solar eclipse proved that the observed depletions in the TEC time series and profiles were caused by the annular solar eclipse.     

Response of the lower atmosphere to intense geomagnetic storms

In this short paper we examine the possible connection between atmospheric parameters measured at low and middle altitudes and geomagnetic storms occurred in 2000 and 2003. For that, from a chain of stations located near the meridian 60°W we compare the storm time values of temperature and wind speed with their standard deviation 2σ obtained from quiet time values. We observed statistically significant variations at several altitudes during the storm recovery phase and after it, both in neutral wind speed and temperature. The results obtained suggest that atmospheric parameters could be affected by geomagnetic storms.     

Deformation of the ionosphere structure during the space weather events on October–November 2003

We present the spatial maps of the ionosphere–plasmasphere slab thickness τ (ratio of the vertical total electron content, TEC, to the F-region peak electron density, NmF2) during the intense ionospheric storms of October–November 2003. The model-assisted technology for estimate of the upper boundary of the ionosphere, hup, from the slab thickness components in the bottomside and topside ionosphere – eliminating the plasmasphere contribution of τ – is applied at latitudes 35° to 70°N and longitudes −10° to 40°E, from the data of 20 observatories of GPS-TEC and ionosonde networks, for selected days and hours of October and November 2003. The daily–hourly values of NmF2, hmF2 and TECgps are used as the constrained parameters for the International Reference Ionosphere extended to the plasmasphere, IRI-Plas, during the ionospheric quiet days, positive and negative storm phases for estimate of τ and hup. Good correlation has been found between the slab thickness and the upper boundary of the ionosphere for the intense ionospheric storms at October–November 2003. During the negative phase of the ionospheric storm, when the ionospheric plasma density is exhausted, the nighttime upper boundary of the ionosphere is greatly uplifted towards the magnetosphere tail, while the daytime upper boundary of the ionosphere is reduced below 500 km over the Earth.     

Seasonal variations in vertical distribution of meteor decay time as observed from meteor radars at 8.5°N and 80°N

The decay times of meteor radar echoes have been used for decades to investigate characteristics of the mesosphere and lower thermosphere (MLT) region. As the meteor echo decay time depends on background atmospheric parameters, in the present communication, we examine the seasonal variation of the vertical distributions of underdense meteor echo decay times with respect to echo strength. Observations from two similar radars located at two distinct geographical locations, Thumba (8.5°N, 77°E) and Eureka (80°N, 85.8°W) were used for the present study. Here, the radar received signal power is categorized into strong and weak echoes and vertical profiles of their decay times are constructed. It has been noticed that the monthly mean decay time vertical profile turning altitude (i.e., inflection point) varies in the range of 80–87?km of altitude depending on latitude. The turning altitude is observed at relatively lower heights in the winter than in summer at both the latitudes. The present analysis shows that the meteor decay time below the mean turning altitude follows a decreasing trend with decreasing altitude, which is quite distinct to the behaviour of ambipolar diffusion. It is also observed that there is a difference in mean decay time of strong and weak echoes below 90?km of altitude, which is very prominently seen at lower altitudes. This difference shows a seasonal pattern at high latitude, but does not show any seasonal variation at low latitude. The present results are discussed in light of current understanding of the meteor decay time.     

Modelled ionospheric Te profiles at mid-latitudes for possible IRI application

The electron temperature (T_e) variation in the mid-latitude ionosphere at altitudes between 120 – 800 km has been modelled for various seasonal and solar-cycle conditions. The calculated electron temperatures are consistent with plasma densities and ion temperatures computed from a time-dependent ionospheric model. The T_e distribution can be represented by a subset of standard T_e profiles. T_e above 200 km is controlled by the magnetospheric heat flux into the ionosphere. For realistic values of the magnetospheric heat flux, the maximum electron temperature ranges from 3000 to 10,000 K at 800 km. The effect of increasing the heat flux is to increase the topside temperature but retain the profile shape. Hence, given a topside T_e observation and selection of an appropriate profile shape, the entire T_e distribution can be computed.