Mars Global Surveyor aerobraking: Atmospheric trends and model interpretation

Mars Global Surveyor (MGS) recently obtained coordinated lower-atmosphere (thermal and dust) measurements and simultaneous upper atmosphere accelerometer data (densities, scale heights and temperatures) for the purpose of safely aerobraking the spacecraft toward its mapping orbit (Keating et al. 1998). Much useful scientific information was also gleaned that describes the coupling of these atmospheric regions during this Phase I aerobraking period (September 1997–March 1998; Ls = 184–300). The major features of this aerobraking data are presented, and its trends elucidated in order to: (1) illustrate the aerobraking environment experienced by the spacecraft, and (2) decompose the processes responsible for the atmospheric variations observed. Coupled general circulation models of the Mars lower and upper atmospheres are exercised to investigate the solar-orbital, seasonal, wave, and dust variations observed during MGS aerobraking. The precession of the MGS periapsis position during Phase I permits longitudinal, latitudinal, local time, and vertical variations of the thermosphere to be monitored. Future aerobraking activities at Mars will benefit greatly from this MGS aerobraking data and its model interpretation.     

Comparison of the observed and modeled low- to mid-latitude thermosphere response to magnetic activity: Effects of solar cycle and disturbance time delay

The performance of JB2008 and NRLMSISE-00 models, in describing the response of the thermosphere to magnetic activity are evaluated against total mass density retrieved from accelerometer measurements made onboard CHAMP satellite during 5 years. We show that the global low- to mid-latitude disturbance amplitude is correctly described by the JB2008 model for low solar activity conditions and by both the JB2008 and the NRLMSISE-00 models for high solar activity conditions. For low solar activity conditions, statistics based on almost 3 years of data confirm the large underestimation by the NRLMSISE-00 model quantified by Lathuillère et al. (2008) for the year 2004. We also found that the time delay between low- to mid-latitude global thermosphere disturbance and magnetic activity is statistically well estimated by the NRLMSISE-00 and JB2008 models for disturbed conditions. For moderately disturbed conditions however, the time delay estimated by the JB2008 model is too large by about 3 h. For very disturbed conditions, we found different time delays during day-time and night-time, using new geomagnetic proxies with a 30-min time resolution.     

Characteristics of solar particle events at Mars relative to Earth

While not specifically designed to detect solar energetic particle radiation, the Electron Reflectometer onboard Mars Global Surveyor (MGS/ER) collected such data from January 1999 through October 2006. Energetic protons (?25 MeV) and other ions penetrated the MGS/ER shielding and registered counts within the instrument’s electronics. During solar particle events (SPE’s), prolonged enhancements in the particle background were observed at Mars with time intensity profiles similar to Earth based SPE observations. Throughout the lifespan of MGS/ER, 85 distinct SPE’s were observed. Basic characteristics of Mars based SPE observations and the frequency of SPE occurrences at Mars are compared to corresponding Earth based observations. Approximately 22% of SPE’s that occurred during MGS/ER operation were observed at Earth but not Mars. Similarly, 19% of SPE’s were observed at Mars but not Earth. Time intensity profiles at Earth and Mars match predictions provided in the literature, based on the physical location of the detector with respect to the motion of the interplanetary shock wave. Note: The work described herein was largely conducted as part of a doctoral dissertation produced by the author.     

Broglio Drag Balance for neutral thermosphere density measurement on UNICubeSAT

The nanosatellite UNICubeSAT is described, carrying a Broglio Drag Balance Instrument for neutral thermosphere density in situ measurements. The aim of the mission is to contribute to the development of accurate thermosphere models, achieving in situ, real time measurements of atmosphere density, that could be exploited for global atmosphere model validation and accurate short term (1–3 days) real time space weather forecasts. The satellite is inexpensive and swarms could be easily launched operating as a distributed sensor network to get simultaneous in situ local (not orbit averaged) measurements in multiple positions and orbit heights. The nanosatellite is based on the Cubesat standard architecture, weighing about 1 kg for 1-L volume. Atmospheric drag force is measured by the displacement of light plates exposed to the incoming particle flux seen by the spacecraft, applying the original three dimensional Broglio Drag Balance concept to a single nanosatellite axis. The instrument concept and its relation to the satellite bus is depicted, showing that many long term potential measurement error sources and biases can be removed in data processing if the spacecraft is spin stabilized. The expected accuracy in density measurements is 20%. The instrument cost is a fraction of that of accurate accelerometers. The onboard systems are based on commercial off the shelf components, in accordance with the short lifetime typical of aeronomy satellites.     

Mid-latitude thermospheric zonal winds during the equinoxes

The diurnal variation of the mid-latitude upper thermosphere zonal winds during equinoxes has been studied using data recently generated from CHAMP measurements from 2002 to 2004 using an iterative algorithm. The wind data was separated into two geomagnetic activity levels, representing high geomagnetic activity level (Ap > 8) and low geomagnetic activity level (Ap ? 8). The data were further separated into two solar flux levels; with F10.7 > 140 for high and F10.7 ? 140 for low. Geomagnetic activity is a correlator just as significant as solar activity. The response of mid-latitude thermospheric zonal winds to increases in geomagnetic disturbances and solar flux is evident. With increase in geomagnetic activity, midday to midnight winds are generally less eastward and generally more westward after the about midnight transitions. The results show that east west transitions generally occurred about midnight hours for all the situations analyzed. The west to east transition occurs from 1400–1500 MLT. Enhanced westward averaged zonal wind speeds going above 150 ms⁻¹ are observed in the north hemisphere mid-latitude about sunrise hours (∼0700–1100 MLT). Nighttime winds in the north hemisphere are in good agreement with previous single station ground observations over Millstone Hill. Improved ground observations and multi satellite observations from space will greatly improve temporal coverage of the Earth’s thermosphere.     

Nitric oxide density enhancements due to solar flares

A differential emission measure technique is used to determine flare spectra using solar observations from the soft X-ray instruments aboard the Thermosphere Ionosphere Mesosphere Energetics Dynamics and Solar Radiation and Climate Experiment satellites. We examine the effect of the solar flare soft X-ray energy input on the nitric oxide (NO) density in the lower thermosphere. The retrieved spectrum of the 28 October 2003 X18 flare is input to a photochemical thermospheric NO model to calculate the predicted flare NO enhancements. Model results are compared to Student Nitric Oxide Explorer Ultraviolet Spectrometer observations of this flare. We present results of this comparison and show that the model and data are in agreement. In addition, the NO density enhancements due to several flares are studied. We present results that show large solar flares can deposit the same amount of 0.1–2 and 0.1–7 nm energy to the thermosphere during a relatively short time as the Sun normally deposits in one day. The NO column density nearly doubles when the daily integrated energy above 5 J m⁻² is doubled.     

Neutral atmosphere and crustal magnetization as factors determining differences in plasma convection and magnetic fields distribution in the ionospheres of Mars and Venus

The pattern of the magnetic field/plasma convection can be, to some extent, recovered from the magnetic field measurements by employing either theoretical or numerical models. We use the MAG/ER day-time measurements of the magnetic field at the altitudes from 90 to 180 km during the elliptical orbits of MGS. Analysis of the altitude variation of the characteristics of the large-scale magnetic fields, which were measured some distance away from strong crustal magnetic anomalies, is summarized. The low density of the Martian atmosphere together with the crustal magnetization result in critical differences in plasma convection which are followed by remarkable differences of the magnetic field features within the ionosphere of Venus and Mars (even in its northern hemisphere where the crustal magnetization is, on the average, low) and distribution of currents.     

Effects of sector structure of the interplanetary magnetic field on the upper mesosphere–lower thermosphere dynamics

In this paper we study the influence of the interplanetary magnetic field (IMF) polarity changes caused by the Earth passing through the IMF sector boundary on the dynamic processes taking place in neutral atmosphere within the altitude interval of the upper mesosphere–lower thermosphere (83–101 km). The analysis has revealed the influence of the IMF sector structure on dynamics of the upper mesosphere–lower thermosphere. There has been a significant seasonal variation of the wind reaction to the IMF polarity changes observed. The influence of the IMF polarity changes on neutral atmosphere dynamics within the altitude range of 83–101 km is most pronounced in the zonal component of neural wind when the IMF polarity changes from negative to positive in all the seasons except for spring and when IMF polarity changes from positive to negative – in spring only.     

Variability study of the crest-to-trough TEC ratio of the equatorial ionization anomaly around 120°E longitude

In this paper, latitudinal profiles of the vertical total electron content (TEC) deduced from the dual-frequency GPS measurements obtained at ground stations around 120°E longitude were used to study the variability of the equatorial ionization anomaly (EIA). The present study mainly focuses on the analysis of the crest-to-trough TEC ratio (TEC-CTR) which is an important parameter representing the strength of EIA. Data used for the present study covered the time period from 01 January, 1998 to 31 December, 2004. An empirical orthogonal function analysis method is used to obtain the main features of the TEC-CTR’s diurnal and seasonal variations as well as its solar activity level dependency. Our results showed that: (1) The diurnal variation pattern of the TEC-CTR at 120°E longitude is characterized by two remarkable peaks, one occurring in the post-noon hours around 13–14 LT, and the other occurring in the post-sunset hours around 20–21 LT, and the post-sunset peak has a much higher value than the post-noon one. (2) Both for the north and south crests, the TEC-CTR at 120°E longitude showed a semi-annual variation with maximum peak values occurring in the equinoctial months. (3) TEC-CTR for the north crest has lower values in summer than in winter, whereas TEC-CTR for the south crest does not show this ‘winter anomaly’ effect. In other words, TEC-CTR for both the north and south crests has higher values in the northern hemispheric winter than in the northern hemispheric summer. (4) TEC-CTR in the daytime post-noon hours (12–14 LT) does not vary much with the solar activity, however, TEC-CTR in the post-sunset hours (19–21 LT) shows a clear dependence on the solar activity, its values increasing with solar activity.     

Seasonal patterns of condensation and sublimation cycles in the cryptic and non-cryptic regions of the South Pole

The South Pole of Mars is characterized by an asymmetric residual ice cap composed of water ice and CO₂ ice. On the opposite side of the residual cap, there exists an area called cryptic region which is relatively free of ice during summer time. Many fan-shaped km-scale structures apparently caused by a wind-blown system of dust-laden gas jets occurred dozens degrees of Ls before the complete sublimation of the CO₂ frost layer. We have examined the seasonal cycles of condensation and sublimation in the cryptic and non-cryptic regions by using the topographic data from the MOLA/MGS measurements. Using the MOLA topography data collected over one Martian year (1999–2001), we have studied the temporal elevation change and the seasonal cycle of the carbon dioxide frost on the southern polar caps. We have produced mapping of the seasonal CO₂ frost thickness variation for seven Ls (30°, 60°, 90°, 120°, 150°, 180°, 210°, 240°, 270° and 330°). It is found that the time variations of the CO₂ frost thickness in these two regions are quite similar. The greatest thickness of the CO₂ frost layer is about 0.76–0.78 m in both places occurs at Ls = 150°.     

Reevaluation of thermosphere heating by auroral electrons

This paper presents a new calculation of neutral gas heating by precipitating auroral electrons. It is found that the heating rate of the neutral gas is significantly lower than previous determinations below 200 km altitude. The neutral gas heating arises from the many exothermic chemical reactions that take place from the ions and excited species created by the energetic electrons. The calculations show that less than half the energy initially deposited ends up heating the neutral gases. The rest is radiated or lost in the dissociation of O₂ because the O atoms do not recombine in the thermosphere. This paper also presents a new way of calculating the heating rate per ionization that can be used for efficient determination of the overall neutral gas heating for global thermosphere models. The heating rates are relatively insensitive to the neutral atmosphere when plotted against pressure rather than altitude coordinates. At high altitudes, the heating rates are sensitive to the thermal electron density and long-lived species. The calculations were performed with the Field Line Interhemispheric Plasma (FLIP) model using a 2-stream auroral electron precipitation model. The heating rate calculations in this paper differ from previous heating rate calculations in the treatment of backscattered electrons to produce better agreement with observed flux spectra. This paper shows that more realistic model auroral electron spectra can be obtained by reflecting the up going flux back to the ionosphere at the upper boundary of the model. In this case, the neutral gas heating rates are 20%–25% higher than when the backscattered flux escapes from the ionosphere.     

Coordinated in situ measurements of plasma irregularities and ground based scintillation observations at the crest of equatorial anomaly

First comparison of in situ density fluctuations measured by the DEMETER satellite with ground based GPS receiver measurements at the equatorial anomaly station Bhopal (geographic coordinates (23.2°N, 77.6°E); geomagnetic coordinates (14.29°N, 151.12°E)) for the low solar activity year 2005, are presented in this paper. Calculation of the diurnal maximum of the strength of the equatorial electrojet, which can serve as precursor to ionospheric scintillations in the anomaly region is also done. The Langmuir Probe experiment and Plasma Analyzer onboard DEMETER measure the electron and ion densities respectively. Irregularities in electron density distribution cause scintillations on transionospheric links and there exists a close relationship between an irregularity and scintillation. In 40% of the cases, DEMETER detects the irregularity structures (dNe/Ne ? 5% and dNi/Ni (O⁺) ? 5%) and GPS L band scintillations (S4 ? 0.2) are also observed around the same time, for the low solar activity period. It is found that maximum irregularity intensity is obtained in the geomagnetic latitude range of 10–20° for both electron density and ion density. As the GPS signals pass through this irregularity structure, scintillations are recorded by the GPS receiver installed at the equatorial anomaly station, Bhopal it is interesting to note that in situ density fluctuations observed on magnetic flux tubes that pass over Bhopal can be used as indicator of ionospheric scintillations at that site. Many cases of density fluctuations and associated scintillations have been observed during the descending low solar activity period. The percentage occurrence of density irregularities and scintillations shows good correspondence with diurnal maximum of the strength of electrojet, however this varies with different seasons with maximum correspondence in summer (up to 66%) followed by equinox (up to 50%) and winter (up to 46%). Also, there is a threshold value of EEJ strength to produce density irregularities ((dNe/Ne)_max ? 5%) and for moderate to strong scintillations (S4 ? 0.3) to occur. For winter this value is found to be ∼40 nT whereas for equinox and summer it is around 50 nT.     

Effect of severe geomagnetic storm conditions on atomic oxygen greenline dayglow emission in mesosphere

Severe geomagnetic storms and their effects on the 557.7 nm dayglow emission are studied in mesosphere. This study is primarily based on photochemical model with the necessary input obtained from a combination of experimental observations and empirical models. The model results are presented for a low latitude station Tirunelveli (8.7°N, 77.8°E). The volume emission rates are calculated using MSISE-90 and NRLMSISE-00 neutral atmospheric models. A comparison is made between the results obtained from these two models. A positive correlation amongst volume emission rate (VER), O, O₂ number densities and Dst index has been found. The present results indicate that the variation in emission rate is more for MSISE-90 than in NRLMSISE-00 model. The maximum depletion in the VER of greenline dayglow emission is found to be about 30% at 96 km during the main phase of the one of the geomagnetic storms investigated in the case of MSISE-90 (which is strongest with Dst index −216 nT). The O₂ density decreases about 22% at 96 km during the main phase of the same geomagnetic storm.The NRLSMSISE-00 model does not show any appreciable change in the number density of O during any of the two events. The present study also shows that the altitude of peak emission rate is unaffected by the geomagnetic storms. The effect of geomagnetic storm on the greenline nightglow emission has also been studied. It is found that almost no correlation can be established between the Dst index and variations in the volume emission rates using the NRLMSISE-00 neutral model atmosphere. However, a positive correlation is found in the case of MSISE-90 and the maximum depletion in the case of nightglow is about 40% for one of the storms. The present study shows that there are significant differences between the results obtained using MSISE-90 and NRLMSISE-00.     

Signatures of the Sudden Stratospheric Warming events of January–February 2008 in Seoul,S. Korea

The period January–February 2008 was characterized by four Sudden Stratospheric Warmings (SSWs) in the Northern Hemisphere, of which the last warming, at the end of February 2008, was a major warming. A significant decrease in mesospheric water vapour (H₂O) of more than 2 ppmv (∼40%) was observed by the ground-based microwave (GBMW) radiometer in Seoul, S. Korea [37.3°N, 126.3°E] during the major SSW. A comparison with ground-based mesospheric H₂O observations from the mid-latitude station in Bern [46.9°N, 7°E] revealed an anticorrelation in the mesospheric H₂O data during the major SSW. In addition, prior to the major warming, strong periodic fluctuations were recorded in the Aura MLS vertical temperature distribution between 15 and 0.05 hPa at Seoul. The mesospheric temperature oscillation was found to have a period of ∼10–14 days with a persistency of 3–4 cycles.     

Lower-mesospheric inversion layers over brazilian equatorial region using TIMED/SABER temperature profiles

Lower-mesospheric inversion layers (MILs) were studied using the temperature profiles observed by TIMED/SABER over Cariri (7.5°S, 36.5°W), Brazil, in 2005. A total 175 MILs were identified with the maximum occurrence in April and October and the minimum in January and July. The lower MIL is located in a height region from 70 to 90 km, with the peak at around 83 ± 4 km with the temperature of 205 ± 5 K, and the thickness of 4–10 km. The results show large amplitudes of MILs during equinoxes and minimum in solstices, with a clear semiannual variation. A general feature of lower MIL in monthly mean profile was observed twice a year, one from February to May, and the other from August to October with a downward shift of the top level. These results suggest that formation and long persistence of MIL is an important factor to investigate propagation of atmospheric gravity waves in the mesosphere-lower thermosphere (MLT) region.     

An investigation of the solar cycle impact on the lower thermosphere O(S) nightglow emission as observed by WINDII/UARS

The yearly variation of the integrated emission rate of the O(¹S) nightglow in the lower thermosphere is studied and the solar cycle impact is examined from the observations of the Wind Imaging Interferometer (WINDII) operated on the Upper Atmosphere Research Satellite (UARS). More than 300,000 volume emission rate profiles of the O(¹S) nightglow observed by WINDII for 40°S–40°N latitudes during November 1991–August 1997 over half of a solar cycle are utilized. These profiles are vertically integrated for the altitude range of 80–100 km and the equivalent column integrated emission rates are then zonally averaged for bins with 10° latitude and 3 month intervals. It is found that for each latitude the O(¹S) nightglow emission rate appears to increase with increasing solar F10.7 cm flux, following a linear relationship. This characterizes the solar cycle impact on the O(¹S) nightglow, while the solar influence is modulated by a seasonal variation. Based on these variations, an empirical formula is derived for predicting the three-month averages of the O(¹S) nightglow integrated emission rate. The standard error of the estimated values from the formula is smaller than 30 Rayleigh.     

The asymmetrical features in electron density during extreme solar minimum

The variations of plasma density in topside ionosphere during 23rd/24th solar cycle minimum attract more attentions in recently years. In this analysis, we use the data of electron density (Ne) from DEMETER (Detection of Electromagnetic Emissions Transmitted from Earthquake Regions) satellite at the altitude of 660–710 km to investigate the solstitial and equinoctial asymmetry under geomagnetic coordinate system at LT (local time) 1030 and 2230 during 2005–2010, especially in solar minimum years of 2008–2009. The results reveal that ΔNe (December–June) is always positive over Southern Hemisphere and negative over northern part whatever at LT 1030 or 2230, only at 0–10°N the winter anomaly occurs with ΔNe (December–June) > 0, and its amplitude becomes smaller with the declining of solar flux from 2005 to 2009. The ΔNe between September and March is completely negative during 2005–2008, but in 2009, it turns to be positive at latitudes of 20°S–40°N at LT 1030 and 10°S–20°N at LT 2230. Furthermore, the solstitial and equinoctial asymmetry index (AI) are calculated and studied respectively, which all depends on local time, latitude and longitude. The notable differences occur at higher latitudes in solar minimum year of 2009 with those in 2005–2008. The equinoctial AI at LT 2230 is quite consistent with the variational trend of solar flux with the lowest absolute AI occurring in 2009, the extreme solar minimum, but the solstitial AI exhibits abnormal enhancement during 2008 and 2009 with bigger AI than those in 2005–2007. Compared with the neutral compositions at 500 km altitude, it illustrates that [O/N₂] and [O] play some roles in daytime and nighttime asymmetry of Ne at topside ionosphere.     

Earthquake induced dynamics at the ionosphere in presence of magnetic storm

The modifications induced in the dynamics of the ionosphere by the major Japan earthquake (EQ) of March 11, 2011 (epicenter at 38.322°N, 142.369°E, M = 8.9) in presence of a magnetic storm are examined by mapping latitudinal variations of F-layer ionization density (NmF2) from 22 stations covering the epicenter zone. The changes forced into the Total Electron Content (TEC) by the major EQ in the magnetic storm ambiance are also examined from the GPS data collected at Guwahati (26° 10′ N, 91° 45’ E), situated in the major fault system of East Asia. The contributions of pre-seismic electric field as well as of magnetic storm time electric field in the observed density variations are brought into the ambit of discussion. The influence of lower atmosphere in shaping TEC features during the study case is highlighted. The effects of solar activity on density variations during such complex ambiances are also addressed.     

Stability of the seasonal variations in diurnal and semidiurnal components of mid-latitude F2 layer parameters

We report work utilizing 15-min resolution ionospheric data obtained with DPS-4 digisonde in 2003–2011 to study the seasonal variations in amplitudes and phases of the most powerful spectral components of the F2 layer critical frequency (foF2) and peak height (hmF2) fluctuations over Irkutsk (52.5°N, 104.0°E). We show that fluctuations of both parameters contain quasi-harmonic components with periods of T_n = 24/n h (n = 1–7). The number of distinct spectral peaks varies from 3 in summer to 7 in winter. Amplitude and phase characteristics of the diurnal (n = 1) and semidiurnal (n = 2) components is studied using the data sets extracted from the original data sets with band-pass filter. It has been found that the amplitudes of diurnal/semidiurnal foF2 and diurnal hmF2 components are maximum in winter and minimum in summer. Amplitudes of the diurnal components vary gradually; those of the foF2 semidiurnal one, abruptly, thus forming a narrow winter maximum in November–January. The phase (local time of maximum) of the diurnal foF2 component increases gradually by 4–6 h from winter to summer. The phase of the semidiurnal foF2 component is nearly stable in winter/summer and sharply decreases (increases) by 2–3 h near the spring (autumn) equinox. The phase of the diurnal component of hmF2 (local time of minimum) varies slightly between 1130 and 1300 LT; that of the semidiurnal one decreases (increases) by 4–6 h from January to March (from September to November). The results obtained show that the main features of seasonal variations in the diurnal and semidiurnal components of the mid-latitude F2 layer parameters recur consistently during the solar activity growth and decline phases.