Next Article in Journal
Mapping Floristic Composition Using Sentinel-2A and a Case Study Evaluation of Its Application in Elephant Movement Ecology in Sagalla, Kenya
Previous Article in Journal
Sequential Processing of Inter-Satellite Doppler Tracking for a Dual-Spacecraft Configuration
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Technical Note

A Novel Method to Derive Exospheric Temperatures from Swarm Thermospheric Densities during Quiet Times

1
State Key Laboratory of Space Weather, National Space Science Center, Chinese Academy of Sciences, Beijing 100190, China
2
University of Chinese Academy of Sciences, Beijing 100049, China
3
Hainan National Field Science Observation and Research Observatory for Space Weather, National Space Science Center, Chinese Academy of Sciences, Beijing 100190, China
4
Faculty of Aerospace Engineering, Delft University of Technology, 2629 HS Delft, The Netherlands
*
Author to whom correspondence should be addressed.
Remote Sens. 2022, 14(21), 5382; https://0-doi-org.brum.beds.ac.uk/10.3390/rs14215382
Submission received: 13 August 2022 / Revised: 20 October 2022 / Accepted: 24 October 2022 / Published: 27 October 2022
(This article belongs to the Section Atmospheric Remote Sensing)

Abstract

:
One of the most important parameters in the atmosphere, the neutral temperature, becomes difficult to measure at high altitudes such as the exosphere. Therefore, based on the assumption of static equilibrium and isothermal atmosphere, a new method was developed to derive quiet-time exospheric temperatures using neutral atmospheric densities from 470 km to 550 km, which were obtained from the Swarm satellites. The derived neutral temperatures were obtained at an altitude of approximately 500 km in the low and middle latitudes from mid-April 2014 to early August 2014. The results were evaluated with nighttime temperatures from ground-based Fabry Perot Interferometers at 250 km. The mean deviation between the derived temperatures and FPI was 30.80 K and the standard deviation of the mean was 106.20 K. The diurnal variations of the exospheric temperatures, which tended to reach their maximum in the late afternoon, were in good consistency with the NRL-MSISE00 model simulations. This novel method performs well at low and middle latitudes. The greatest source of uncertainty is the mean molecular mass, which is also not well determined at these high altitudes. Hence, a measurement called “Satellite-tethered Mass Spectrometer Detection” was proposed to address this.

Graphical Abstract

1. Introduction

Exospheric temperature is a fundamental physical quantity of the neutral atmosphere, and as such representing temperature and temperature variations in a realistic way is a key requirement for any model of the atmospheric state. To achieve this, observations are needed. However, the exospheric temperature is difficult to measure. The atmospheric density decreases with height, and the dominant neutral components become O, H, and He at higher altitudes, which increases the difficulty of detecting neutral atmospheric temperatures in the exosphere. Some methods used at lower altitudes, such as infrared detection using carbon dioxide and optical detection using airglow, are unsuitable for exospheric measurements. The ground-based Fabry Perot Interferometer (FPI) has become the predominant atmospheric detection equipment in the upper thermosphere, and there are few space-based approaches that exist for obtaining exospheric temperatures. However, FPI observations have limitations in local times, heights, and global distribution, because the number of FPI stations is limited and they only observe at a fixed altitude during nighttime.
Due to the lack of exospheric temperature observations, previous studies have presented a few methods to derive exospheric temperatures from other neutral parameters. For example, using the 300–500 km neutral density detected by the GRACE satellites, Burke determined a linear function of exospheric temperature with height and neutral density and used this relationship in the Jacchia77 model [1]. Furthermore, based on the diffusion function of the Jacchia77 model and Burke’s method, Wise used the neutral density from CHAMP/GRACE and performed temperature iteration to derive exospheric temperatures by forcing the model density to approximate the measured density [2]. Subsequently, he compared the derived temperatures with NRL-MSISE00 (MSIS) model [3] temperatures and exospheric temperatures from the HASDM model, which were all in good agreement. Additionally, Forbes used temperature iteration as well; he changed the input parameter F10.7 flux to ensure that the model density approximated the measured value from CHAMP/GRACE, which resulted in improved model exospheric temperatures [4]. Weimer also used the density from the CHAMP/GRACE observations, together with the MSIS model and a mapping method. They obtained exospheric temperatures and corrected them according to theoretical disturbances based on atmospheric formulas [5]. Weng et al. constructed an exospheric temperature model using the CHAMP density and the MSIS model [6]. The input parameters of Weng’s model were latitude, longitude, local time, geomagnetic index, and F10.7 flux. They obtained temperatures from 2002 to 2020, which were well correlated with those from MSIS and GRACE. In addition, Evans used the N2-LBH radiation band of the GOLD satellite, diffusion equilibria, and isothermal atmospheric theory to calculate exospheric temperatures by obtaining the integral of the volume emission rate along the line of sight using the Chapman function and the scale height of N2 [7]. These attempts provide us with several applicable methods and new thinking about how to derive exospheric temperatures. Nevertheless, they have limits. For example, the methods of Burke, Wise, Forbes, and others depend strongly on model information and are therefore biased to the model results; while the method of Evans has fitting errors.
In view of the importance of a good knowledge of the neutral temperature at and above the thermosphere, and the difficulty to obtain good measurements, alternative methods to derive exospheric temperature are worth exploring. The previously mentioned methods obtain the exospheric temperature through other parameters. However, those methods are greatly depending on the model information. Hence, a new method was developed to derive exospheric temperatures from neutral densities at two different altitudes, to reduce the reliance on models. This method relies on the assumption that the atmosphere is at static equilibrium with isothermal conditions between the two altitudes. The exospheric temperatures at approximately 500 km were calculated using the neutral densities from the Swarm A/B/C satellites, which fly at two different altitudes separated by about 50 km. The derived temperatures are then compared with temperatures from ground-based FPI and model simulations. However, the role of mean molecular mass cannot be ignored in the inversion. There is a large dispersion of mean molecular mass both in temporal and spatial distribution. Assuming a mean molecular mass of 16 and a temperature of 1000 K, a small change of only 0.1 in mean molecular mass results in a temperature change of 7 K, which shows that the mean molecular mass has to be carefully selected. In this work, MSIS model simulations are used to obtain the mean molecular mass, but this might introduce errors, due to the model uncertainties. Therefore, to solve this problem, a detection method called “Satellite tethered Mass Spectrometer Detection” was proposed to obtain densities of specific components in order to determine the mean molecular mass.
Section 2 describes the data, model and method, and also provides an error estimate of the method. Section 3 details the comparisons between the derived temperatures, the ground-based FPI temperatures, and MSIS model simulations to evaluate this method. Section 4 proposes the detection idea, its estimation of the design, and the error evaluation. Finally, Section 5 summarizes this work.

2. Data and Methods

2.1. Data and Model

2.1.1. Swarm Satellites

The Swarm [8] satellites were launched on 22 November 2013, and are still operating today. These three satellites, named Swarm A, Swarm B, and Swarm C, have similar original weights, volumes, and 510 km orbiting heights. Due to the limited availability of accelerometer-derived densities due to accelerometer instrument issues, we use the derived neutral mass densities [9,10] from GPS precision orbit determination data instead in this study. The temporal and spatial resolution of the data points of the GPS-derived densities are 30 s in time, and 2.2° in latitude. However, the 30 s sampling does not reflect the actual signal information of the relatively smooth densities. High-frequency signals are not well recovered and therefore the temporal resolution of the GPS-derived density signal is about 20 min. Since the Swarm satellites are in near-polar orbits, at middle and low latitudes, the satellite tracks have small variations in longitude; at high latitudes, this varies much more rapidly. The daily average altitudes and neutral densities of the Swarm satellites during 2013–2020 are provided in Figure 1. This shows that the satellites were flying closely together initially. Two months later, they separated and were maneuvered to their selected orbits until 15 April 2014. After this, Swarm A and C flew together at an altitude of approximately 462 km with a 1.4° longitude interval; Swarm B was flying at a higher altitude of approximately 511 km. The satellites have gradually declined in orbit over time. In this study, the data range we used is from 15 April 2014 to 11 August 2014. During this period, the relative errors of the GPS-derived densities for the high-flying and lower-flying satellites are estimated to be, respectively 7% and 4% [10].

2.1.2. The Ground-Based Fabry-Perot Interferometer

The Fabry-Perot Interferometer is a multi-beam Interferometer [11,12] that uses the principle of Doppler interference to derive temperature through observed interference fringes. As the light of these airglows is weaker than the sunlight, it is difficult to observe during the daytime. Nighttime ground-based FPI temperatures are, however, reliable, and widely used in thermospheric research. There are many emission bands for FPI, which are located at different altitudes. In this work, the airglow data in the O630.0 nm band was used, because only the detection altitude of this band can reach the thermospheric altitude. The emissions of O630.0 nm are integrated from 200 km to 300 km to derive temperatures, and because the emission of O630.0 nm airglow reaches its maximum at the altitude of 250 km, it is generally considered that the altitude of the detection data is 250 km. We use the measured FPI temperatures at 250 km to evaluate the derived results at 500 km. According to the atmospheric profile provided by MSIS in Figure 2, temperature increases strongly with altitude until a certain altitude, after which this levels off until the variations in temperature become very small. According to MSIS, this altitude was just below 250 km in 2014; hence, the deviations in temperatures between 250 km and 500 km are considered relatively small when compared to the background temperatures.
Table 1 shows the geographical locations of FPIs, and the number of data points used in the comparison. In the comparison between FPI and derived temperature results, the nighttime temperatures from seven ground-based FPI stations in 2014 were used, consisting of the Xinglong station of the Meridian Program in Hebei, China, and six stations of the Madrigal Database, such as Millstone Hill station. FPI observations are strongly affected by the weather and other environmental conditions and become unreliable when it’s rainy or foggy. In such weather, it is difficult to observe the airglow emission clearly. The Madrigal data have a validity indicator, when validity 0 means the quality of the data is good, validity 1 means the data needs to be considered, and validity 2 means the quality of the data is bad. Unreliable data were removed using the data validity indicator from the Madrigal Database and using the FPI cloud map for the Xinglong station. Only the data with validity 0 or 1 from Madrigal and sunny days from Xinglong were used in this work.

2.1.3. NRL-MSISE00 Model

NRL-MSISE00(MSIS) [3] is an empirical model of the neutral atmosphere. The model is based on fundamental equations and observations from non-coherent scatter radars and satellites. It can provide global temperature, density, and number density of different atmosphere components from ground to 2000 km altitude. The input parameters of MSIS are the day of the year, UT, latitude, longitude, f10.7 flux, and geomagnetic Ap index, and the output are neutral atmospheric temperature, atmospheric mass density, and number density of O, O2, N2, H, H2, N, Ar.

2.2. New Method for Temperature Calculation

Thermospheric density ρ varies with height, as follows:
ρ ( z ) = ρ 0 e z H ( z ) ,
H ( z ) = R T v ( z ) m ¯ ( z ) g ( z )
where z is the altitude, H is the scale height, T v is the atmospheric temperature, m ¯ is the mean molecular mass, R is the gas constant, and g is the gravitational acceleration.
The gravitational acceleration g at high altitude was calculated as follows:
g ( r ) = G M r 2
where G is the gravitational constant, M is the mass of the Earth, and r is the distance from the core of the Earth to the target point.
We assume the atmosphere is in static equilibrium, so the scale heights and temperatures are constant in the area between the two Swarm satellites. When substituting ρ 1 and ρ 0 for the densities of the higher and lower flying satellites, Equation (1) becomes:
ρ 1 = ρ 0 e h 1 h 0 H
We used the average of the mean molecular mass and the gravity acceleration of the upper-flying satellite and lower-flying satellites as regional approximations in the calculations, and combined with Equation (2), the temperature in this region was determined by:
T v = H m ¯ g R = ( h 1 h 0 ) · m ¯ g l n ( ρ 1 ρ 0 ) · R
Because in-situ observations of the mean molecular mass are unavailable for the Swarm satellites, we used two different approaches to derive mean molecular mass:(1) Based on Equation (5), the mean molecular mass can be obtained from the known temperatures and densities at two different altitudes, so we simulated the temperatures and densities of Swarm A/C and Swarm B using MSIS, and thus derived the corresponding mean molecular mass m ¯ s along the Swarm orbits; (2) We calculated the global mean molecular masses in 2014 using the number densities of different atmospheric components from MSIS. The selected time period was from April 2014 to June 2014, which is comparable with the time period of the FPI observations, and the latitude and longitude ranges were from −60° to 60° geomagnetic latitude, 0°−360° geomagnetic longitude, 5°, and 10° grids, respectively. We averaged all the values of mean molecular mass to obtain a constant m ¯ 0 = 16.0318 with a standard deviation of 1.2%. It should be noted that the constant m ¯ 0 approach was used in order to further reduce the dependence on the MSIS model input.
In both cases, the gravitational acceleration was computed along the Swarm orbit using the average altitude of the upper-flying and lower-flying satellites.
To find conjunctions between Swarm B and the lower satellite pair, a single-point sliding window was used. The window was centered on Swarm B, with a latitudinal width of ±2.5°, a longitudinal width of ±7.5°, and a time range of ±5.5 min. When Swarm A or Swarm C entered the window, they were regarded as having the same observation conditions as Swarm B.
Using the sliding window mentioned above, we obtained 81011 points of temperature from 2014 to 2020 at an altitude of approximately 500 km. Due to complex disturbances and a lack of FPI observations in the high latitudes and polar regions, only the temperatures at the low and middle latitudes were considered, which are shown in Figure 3. Low latitude conjunctions between the higher-flying satellite and the lower pair only occurred from mid-April to early August in 2014; they otherwise occurred in the polar regions. Due to the orbit plane configuration of Swarm during this time span, a period of up to 19 h of low latitude conjunctions occurred every 6 days. Therefore, the inversion results were obtained every sixth day from mid-April to early August in 2014, and the frequency of occurrence of the derived results is half an hour during the 19 h of the conjunction. In the following section, the data for these three months are further analyzed.

2.3. Sources of Uncertainty

If ( ρ 1 ρ 0 ) is set to x, the Taylor expansion of Equation (5) is given as
T = T 0 + Δ z · m · g l n 2 ( x ) · R · 1 x · Δ x + Δ z · g l n ( x ) · R · Δ m + Δ z · m l n ( x ) · R · Δ g ,
Thus, the relative error of this method can be expressed as
Δ T T = 1 l n ( x ) Δ x x + Δ m m + Δ g g
Equation (7) shows that the relative error consists of three parts, the measurement error of the Swarm densities, the error of the mean molecular mass, and the error of gravitational acceleration. For the Swarm densities in 2014, the relative error is estimated to be about 7% and 4% for, respectively, Swarm B and the lower pair. The mean value of the relative error of ( ρ 1 ρ 0 ) is then about 3.01% ± 11.015%, and 1 l n ( x ) becomes 1.30.
The mean molecular mass is obtained using MSIS model information. MSIS is a climatological empirical atmospheric model, and it shows a good long-term variation of atmospheric temperatures and neutral compositions. However, there are deviations between the MSIS model and the actual atmosphere and this introduces uncertainties. It is difficult to determine the uncertainty of the modeled mean molecular mass, because there are no actual data for comparison. The error of the mean molecular mass due to the regional approximation, which is a result of the averaging between altitudes, can be estimated. To access this, the global mean molecular masses at the average altitudes of upper-flying and lower-flying satellites (about 477 km and 525 km) are obtained from MSIS similar to (1) in Section 2.2. The deviations between the nominal mean molecular masses and the global mean values at these two altitudes can be used to determine the relative error of the mean molecular mass, which is 3.25% ± 9.89%. In the same way, the gravitational acceleration at the average altitudes of the upper-flying and lower-flying satellites are determined, which leads to a relative error of gravitational acceleration of 0.7%. Substituting these three errors into Equation (7), the total relative error of this method becomes 7.86% ± 17.4%.

3. Temperature Comparisons between Observations, Model Simulations, and Inversion Temperatures

The derived temperatures from 15 April 2014 to 11 August 2014 were used for analysis, and their variations in latitude and local time are shown in Figure 3. It should be noted that temperature data are not available every day, but at a frequency of every 6 days, so there are only 30 days of data in this period. The results cover 0°−360° in longitudes and −60° to 60° in latitude and do not cover 24 h completely. Figure 4 shows there is a lack of temperatures in LT 2:00−6:00 and LT 14:00−16:00. In this section, we evaluate the derived temperatures using statistical deviations and diurnal variation characteristics. For brevity of the text, we refer to the derived temperatures calculated with m ¯ s by T ( m ¯ s ) , the derived temperatures calculated with m ¯ 0 by T ( m ¯ 0 ) , the temperatures of FPI by T ( F P I ) and the temperatures of MSIS by T ( M S I S ) .

3.1. Quantitative Comparison between Derived Temperatures, Ground-Based FPI Measured Temperatures and MSIS Simulations

Since there were no derived temperatures at the exact locations of the FPI stations, a window was again used to obtain conjunctions. The window was centered on the locations of the FPIs, with latitudinal and longitudinal widths of ±2.5° and ±7.5°, respectively. If the derived temperatures were located inside the window, they were considered to be under the same observation conditions. Additionally, cubic-spline interpolation was used to obtain the FPI temperatures at the same time as the derived temperatures. To minimize interpolation errors, a relatively small two-hour range was selected to fit a smooth trend. In addition, since the atmospheric conditions during the disturbance period do not support our assumption, those data points with Kp > 2 or Dst ≤ −30 are discarded. It should be noted that the effect of geomagnetic disturbances on neutral density can last for a few hours [13,14], therefore, ideally the geomagnetic indexes of both the current day and the previous day need to be considered. However, in this case, the remaining data are not enough for subsequent analysis. Therefore, only the geomagnetic index of the current day was considered. In this way, 203 data points were obtained for comparison (see Appendix A), and they were all at nighttime.
For the comparison between the derived temperatures and the MSIS comparison, MSIS model temperatures were determined corresponding to each derived temperature. In this case, the whole database of 81011 points is used for the comparison, including daytime and nighttime data.
The comparison was evaluated by the mean deviation and the standard deviation of the mean. First, a point-to-point deviation of the derived temperatures, FPI temperatures, and MSIS temperatures was obtained; from these differences, the mean deviations and the standard deviations of the mean were calculated. Table 2 shows the results of the statistical analysis of the derived temperatures calculated with different mean molecular masses compared with FPI and MSIS temperatures, respectively.
The comparison results of the derived temperatures and the FPI temperatures are analyzed first. As shown in Appendix A, the range of derived temperatures and the FPI temperatures was in both cases between 850 K to 1100 K. Table 2 shows that the mean deviation between T ( m ¯ s ) and T ( F P I ) is 30.80 K, indicating that T ( m ¯ s ) was larger than T ( F P I ) in general. According to MSIS, the expected temperature deviations between 250 km and 500 km are small. The mean temperature deviation between 250 km and 500 km simulated by MSIS is 23.52 K, which is close to the mean deviation between T ( m ¯ s ) and T ( F P I ) . However, the standard deviation between T ( m ¯ s ) and T ( F P I ) is 106.20 K, while that of MSIS is 14.53 K. Such a large great difference is probably due to the fact that temperature variations in MSIS simulations are much smaller than the variations in the actual atmosphere. The differences between T ( m ¯ s ) and T ( F P I ) are considered to be caused both by the temperature deviation between 250 km and 500 km altitude and the uncertainties of our method, which were estimated in Section 2.3.
The expected error of the method can be determined, as explained in Section 2.3. Using the average temperature of the samples in Table 2, the error is 73.98 K ± 179.78 K. The deviation between T ( m ¯ s ) and T(FPI) is 30.80 K ± 106.20 K, which is smaller than the estimated method error.
The negative mean deviation between T ( m ¯ 0 ) and T(FPI) indicates that the derived temperatures are lower, which does not conform to the theoretical relation that the temperature at 500 km is larger than at 250 km. The standard deviation between T ( m ¯ s ) and T(FPI) is 108.03 K, which is close to that between T ( m ¯ s ) and T(FPI), but slightly larger. From these results, we conclude that the use of the constant m ¯ 0 increases the bias. It is shown that the mean molecular masses have a great influence on the derived temperature, and more reliable results are obtained with m ¯ s .
Moreover, the last two rows of Table 2 show that the mean deviation between T ( m ¯ s ) and T(MSIS) is 16.52 K ± 84.46 K, while the mean deviation between T ( m ¯ 0 ) and T(MSIS) is 45.90 K ± 140.73 K. The better agreement between T ( m ¯ s ) and T(MSIS), compared to T ( m ¯ 0 ) , also supports the use of m ¯ s .

3.2. Comparison between Inversion Temperatures Calculated with m ¯ s , Ground-Based FPI Measured Temperatures and MSIS Simulations

The derived temperatures can also be evaluated by analyzing the diurnal variations of FPI temperatures, MSIS simulations, and derived temperatures at different latitudes. It should be noted that because the neutral temperature is sensitive to solar radiation, there are different variation characteristics in different months. Unfortunately, there are not enough data to show a complete diurnal variation per month, so all the data from April to August were put together to obtain an almost 24-h LT coverage. Due to this, our results are expected to be slightly different when compared to MSIS, and therefore only the general trends are analyzed in this evaluation.
Based on the distribution of temperatures in latitudes, the results were divided into 8 latitudinal bands, with about 10,000 points per latitude band. In order to get clear variation trends, the temperatures were averaged every 0.5 h from LT0−LT24 in each latitudinal band and plotted with an error bar. Figure 5 and Figure 6 show the diurnal variations in temperatures of MSIS, FPI, and derived temperatures calculated with m ¯ s and m ¯ 0 , where the horizontal axis is local time, and the vertical axis is temperature. The black dots indicate the derived temperatures, the red dots are the MSIS simulated temperatures, the red smoothed curves are the variation trends fitted with red dots, and the black triangles are the FPI observations.
Figure 5 shows that during the daytime, the diurnal trends of the derived temperature are greatly consistent with the point-by-point MSIS simulations at latitudes except for LT8−LT9. At night, in the northern hemisphere, the derived temperatures in the first half of the night at low latitudes are generally lower than those obtained with MSIS simulations, and in the second half of the night, they are generally higher. At middle latitudes, they are higher than the MSIS temperatures. On other hand, in the southern hemisphere, the derived temperatures are smaller than the MSIS simulations in both the first and second half of the night. In addition, though the FPI observations are mainly distributed in the 30°N−45°N latitude zone, the derived temperatures in this zone also show a similar nighttime variation trend.
The enhancement of LT8−LT9 shown in Figure 5 is stronger at low latitudes than at middle latitudes, and stronger in the southern hemisphere than in the northern hemisphere. Based on earlier findings, the same enhancement of LT8−LT9 was observed by the incoherent scattering radar [15] in Millstone Hill which is located at 42°N. These observations show the relation between enhancement and solar activities. At present, however, there are not enough data to further analyze this.
In Figure 6, T ( m ¯ 0 ) also showed diurnal variation, as well as similar nighttime characteristics as T ( m ¯ s ) in Figure 4. However, in the daytime, the derived temperatures are very different from the MSIS simulations, especially in the southern hemisphere. This might be caused by the use of mean molecular mass m ¯ 0 . m ¯ 0 was calculated with the mean molecular masses in the spring and summer of the northern hemisphere and in the autumn and winter of the southern hemisphere. The resulting average m ¯ 0 is then much larger than the value for the southern hemisphere.
In general, the derived temperatures all have diurnal trends and reach their maximum in the afternoon (around LT14). As there was no data in LT2−LT6, and the theoretical minimum usually occurs during this period, the minimum of the diurnal variation was not analyzed. The trends of T ( m ¯ s ) and T ( m ¯ 0 ) are very different during the day. In general, T ( m ¯ s ) is greatly consistent with MSIS simulations but differs in some details. This indicates that T ( m ¯ s ) performs qualitatively well, while T ( m ¯ 0 ) shows more obvious problems.
The two derived temperatures vary considerably due to the difference between the mean molecular mass. Figure 7 shows the diurnal variation of m ¯ s . It can be found that m ¯ s reaches its maximum around LT18 and its minimum around LT8. Moreover, the difference of m ¯ s between the northern and southern hemispheres is considerable, with lower m ¯ s in the southern hemisphere. Figure 7 shows that the center of the maximum is in the northern hemisphere and the center of the minimum is in the southern hemisphere. Comparing Figure 5 and Figure 6, the difference between the two derived temperatures is mainly in the local time from morning to before noon, during which the value of T ( m ¯ 0 ) is larger than that of T ( m ¯ s ) . According to Figure 7, the difference is obviously caused by the variation of mean molecular mass. Thus, it should be noted that the mean molecular mass plays an important role in the inversion, and its variation with time and location is not negligible.

4. A New Concept “Satellite-Tethered Mass Spectrometer Detection”

Information on the mean molecular mass is required for our method to derive temperatures. Since in situ observations are not available, we had to use MSIS simulations. In order to reduce the model dependency, a constant mean molecular mass was also used to derive temperatures. However, this value causes a large deviation in the obtained results. To determine the mean molecular mass, a detection idea, “Satellite-tethered Mass Spectrometer Detection”, is proposed. Its schematic is shown in Figure 8.
The concept consists of one satellite that carries two independent mass spectrometers that are separated by a certain distance and are used to detect the density at two different altitudes for multiple atmospheric components, such as O, H, etc. From these observations, the mean molecular mass can be determined, and combined with Equation (5) in Section 2.2, the neutral temperature of these different components can be derived.
Figure 9 shows the relationship between the estimation error, the wire length, and the densities measurement errors. The error estimation is based on Equation (7), and the influence of both the wire length, varying between 1−50 km, and a relative density measurement error in the 1−20% range on the total error were investigated. If the total error should be limited to 10%, and if the measurement error of the mass spectrometers is set to be 3%, under the worsts estimation (the relative error of ( ρ 1 ρ 0 ) reaches its maximum), the relative error of densities could reach 6%. In this case, the wire length should be at least 15 km. This is a preliminary estimation, and the wire length should be reconsidered according to the actual situation and the requirements for the total error. With this detection method, there is no need to use a model to simulate the mean molecular mass. Possibly, some extra observations can also be obtained, such as the neutral temperatures of different components.

5. Conclusions

The exospheric temperature is an important parameter, which is difficult to obtain. Therefore, a new method was developed to derive exospheric temperatures using densities at two different altitudes from the Swarm satellites, which reduces the reliance on the model information. Based on the assumption of atmospheric static equilibrium and an isothermal atmosphere, 81,011 points of neutral temperatures from mid-April 2014 to early August 2014 were obtained. These temperatures were mainly distributed at low and middle latitudes, at altitudes of around 500 km. By comparing with ground-based FPIs’ nighttime temperatures at 250 km and corresponding NRL-MSISE00 model temperatures, the derived results were evaluated, and the following conclusions were drawn:
(1) The derived temperatures and FPI temperatures both lie between 850 K to 1100 K at nighttime. The deviation between the derived temperatures calculated with m ¯ s and FPI is 30.80 K ± 106.20 K. The deviation between the derived temperatures calculated with m ¯ 0 and FPI is −30.67 K ± 108.03 K. We also simulated the temperatures deviation between 250 km and 500 km by MSIS, the deviation is 23.52 K ± 14.53 K. The mean deviation between T ( m ¯ s ) and T ( F P I ) is close to the expected deviation between 250 km and 500 km obtained by the MSIS model, while the deviation between T ( m ¯ 0 ) and T ( F P I ) is larger. This indicates that T ( m ¯ s ) is more reliable. Moreover, the deviation between T ( m ¯ s ) and T ( M S I S ) is 16.52 K ± 84.46 K; while the deviation between T ( m ¯ 0 ) and T ( M S I S ) is 45.90 K ± 140.73 K. This shows that T ( m ¯ s ) agrees better with MSIS.
(2) The diurnal temperature trends from April to August 2014 were used to evaluate the derived results. The diurnal variation trends show that both T ( m ¯ s ) and T ( m ¯ 0 ) reach a maximum in the late afternoon. The diurnal variation trends of T ( m ¯ s ) are greatly consistent with MSIS, with few differences in details. A noticeable difference is the temperature enhancement present in T ( m ¯ s ) at LT8−LT9; moreover, the T ( m ¯ s ) is slightly smaller than MSIS in the nighttime before midnight and slightly larger than MSIS in the nighttime after midnight at low latitudes in the Northern hemisphere. At middle latitudes, in the Northern hemisphere, T ( m ¯ s ) is higher than MSIS at nighttime. However, in the Southern Hemisphere, it is lower than MSIS at nighttime. Moreover, in the band of 30°N–45°N, the derived temperatures and FPI temperatures also have similar variation trends at nighttime. In addition, T ( m ¯ 0 ) differs significantly from MSIS during the daytime in the Northern Hemisphere and throughout the Southern Hemisphere.
(3) The mean molecular mass has a great influence on the derived temperatures. As this information is difficult to obtain from observations, it is the largest source of uncertainty in the method. In order to improve this and obtain more accurate exospheric temperatures, a measurement concept called “Satellite-tethered Mass Spectrometer Detection” was proposed. It detects the number densities of different atmospheric components like O, H, He et al. through a satellite carrying two mass spectrometers at a certain separational distance, and uses this information to derive exospheric temperatures. This detection does not need to use the simulations of the mean molecular mass in the calculation, also can obtain the exospheric temperatures of different atmospheric components, and therefore can be used to improve atmospheric empirical models.
The exospheric temperature retrieval method performs well in low and middle latitudes during quiet times. However, for polar regions and disturbance periods, there are currently not sufficient data for analysis. The orbits of the Swarm constellation change over time, and a period of regular low latitude conjunctions between Swarm B and the lower pair was again obtained during the summer of 2021. Further studies can be performed when more conjunction data become available.

Author Contributions

J.X.: Conceptualization, methodology; Y.L.: calculation, formal analysis, writing-original draft preparation; W.Y.: inversion FPI temperature data provision; J.v.d.I.: inversion Swarm neutral density provision, writing—review and editing. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (41831073, and 42174196), the Chinese Meridian Project, and the Specialized Research Fund for State Key Laboratories, the Project of Stable Support for Youth Team in Basic Research Field, CAS (YSBR-018) and the National Key R&D program of China (2021YFE0110200), the International Partnership Program of Chinese Academy of Sciences. Grant No. 183311KYSB20200003.

Data Availability Statement

Publicly available datasets were analyzed in this study. They can be found here: the thermospheric neutral density data of Swarm satellites are available online (https://swarm-diss.eo.esa.int/#swarm%2FLevel2daily%2FLatest_baselines%2FDNS%2FPOD); the global atmospheric model NRL-MSISE00 can be obtained online (through the website https://ccmc.gsfc.nasa.gov/pub/modelweb/atmospheric/msis/nrlmsise00/). [Dataset] National Space Science Center, Chinese Academy of Sciences. 2010-06-25. Inversion data from FPI at Xinglong station. National Space Science Data Center. V1.0. DOI:10.12176/01.05.009.

Acknowledgments

We thank the Chinese Meridian Project for the support fund. We thank Madrigal database for inversion FPI temperature data and Google Earth for the earth picture. We are also grateful to the anonymous reviewers for their constructive comments and suggestions to improve this manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Comparative data of inversion temperature and FPI temperature calculated with m ¯ 0 in 2014.
Table A1. Comparative data of inversion temperature and FPI temperature calculated with m ¯ 0 in 2014.
Name and Location of FPI StationsDOYUT, hLT, hLocation of Swarm B,°Average Location,
°
ALT-B, kmAverageALT, kmTemperature of FPI, kTemperature of Swarm, k
Location of Swarm A/C, °ALT-A/C, km
Arecibo (18.34°,
−66.75°)
1125.45001.383716.30, −6215.85,
−61
515.664921088.68887.57
15.39 −59.99468.33
1125.45421.339217.66, −63.4616.53,
−61.73
515.66492.011088.73837.66
15.39 −59.99468.35
1125.45421.338515.73, −63.4316.51,
−61.74
515.674921088.73894.91
17.29, −60.04468.33
1125.45421.386216.3, −6216.8,
−61.02
515.674921088.73922.71
17.29, −60.04468.33
1125.45831.341717.66, −63.4617.48,
−61.75
515.67492.011088.77868.91
17.29, −60.04468.35
1125.45831.389318.22, −62.0317.76,
−61.04
515.67492.021088.77893.26
17.29, −60.04468.36
1125.46251.344819.58, −63.4918.44,
−61.77
515.67492.041088.82843.63
17.29, −60.04468.4
1125.46251.344517.66, −63.4618.43,
−61.77
515.71492.031088.82903.36
19.2, −60.08468.35
1125.46251.392218.22, −62.0318.71,
−61.06
515.71492.041088.82929.69
19.2, −60.08468.36
1125.46671.347719.58, −63.4919.39,
−61.79
515.71492.061088.86876.11
19.2, −60.08468.4
1125.46671.395320.15, −62.0619.68,
−61.07
515.71492.071088.86898.58
19.2, −60.08468.43
1125.47081.351221.51, −63.5120.36,
−61.8
515.71492.11088.9849.12
19.2, −60.08468.49
1125.47081.350519.58, −63.4920.34
−61.81
515.79492.11088.9912.78
21.1, −60.12468.4
1125.47081.398220.15, −62.0620.63,
−61.09
515.71492.111088.9937.18
21.1, −60.12468.43
Easter Kentucky (37.75°,
−84.29°)
1433.858322.697334.74, −80.5535.34,
−77.42
517.39493.021038.311063.6
35.93, −74.28468.65
1433.862522.701836.66, −80.5436.3,
−77.41
517.39493.181038.621037.24
35.93, −74.28468.97
1433.866722.705736.66, −80.5437.25,
−77.42
517.73493.351038.481062.43
37.83, −74.29468.97
1433.870822.710538.58, −80.5238.21,
−77.41
517.73493.521038.341035.98
37.83, −74.29469.31
1433.87522.714738.58, −80.5239.16,
−77.41
518.09493.71038.21060.56
39.73, −74.29469.31
1433.879222.719540.5, −80.540.12,
−77.4
518.09493.881038.061034.47
39.73, −74.29469.66
1682.054220.562534.51, −85.2935.43,
−82.38
517.96493.271081.03944.58
36.34, −79.46468.57
1682.054220.563236.43, −85.2835.44,
−82.37
517.65493.271081.03903.84
34.44, −79.45468.89
1682.058320.518335.67, −86.7436.01,
−83.1
517.96493.351081.15910.72
36.34, −79.46468.74
1682.058320.56736.43, −85.2836.39,
−82.37
517.96493.431081.15926.49
36.34, −79.46468.89
1682.062520.523237.59, −86.7236.97,
−83.09
517.96493.521081.2893.85
36.34, −79.46469.07
1682.062520.570836.43, −85.2837.34,
−82.38
518.28493.591081.2949.9
38.24, −79.47468.89
1682.062520.571538.35, −85.2737.35,
−82.37
517.96493.591081.2909.78
36.34, −79.46469.22
1682.066720.52737.59, −86.7237.92,
−83.1
518.28493.681081.25915.91
38.24, −79.47469.07
1682.066720.575338.35, −85.2738.3,
−82.37
518.28493.751081.25932.6
38.24, −79.47469.22
1682.070820.531839.51, −86.738.88,
−83.09
518.28493.851081.3898.69
38.24, −79.47469.41
1682.070820.579538.35, −85.2739.25,
−82.37
518.62493.921081.3955.56
40.14, −79.47469.22
1682.070820.580240.27, −85.2539.26,
−82.36
518.28493.931081.3915.65
38.24, −79.47469.57
1682.07520.53639.51, −86.739.83,
−83.09
518.62494.021081.34920.31
40.14, −79.47469.41
1682.07520.584440.27, −85.2540.21,
−82.36
518.62494.11081.34938.06
40.14, −79.47469.57
1742.087520.037935.28, −94.1835.57,
−90.75
517.65493.16990.56920.98
35.86, −87.31468.67
1742.091720.041735.28, −94.1836.52,
−90.75
517.97493.32990.68930.27
37.76, −87.32468.67
1742.091720.042437.2, −94.1736.53,
−90.74
517.65493.32990.68910.9
35.86, −87.31468.99
1742.095820.046237.2, −94.1737.48,
−90.75
517.97493.48990.88920.11
37.76, −87.32468.99
1742.120.0537.2, −94.1738.43,
−90.75
518.3493.65991.08928.83
39.66, −87.33468.99
1742.120.050739.12, −94.1638.44,
−90.74
517.97493.65991.08910.08
37.76, −87.32469.33
1742.104220.054539.12, −94.1639.39,
−90.75
518.3493.82991.27918.73
39.66, −87.33469.33
Peach mountain (42.27°,
−83.75°)
1682.095820.560243.35, −86.6444.59,
−83.04
519.7494.921036.22949.87
45.83, −79.43470.13
1682.095820.560945.26, −86.5944.6,
−83.03
519.33494.921036.22912.71
43.93, −79.46470.51
1742.12520.078742.95, −94.144.15,
−90.7
519.35494.7914.79920.76
45.35, −87.29470.04
1742.12520.079744.87, −94.0544.16,
−90.68
518.99494.7914.79905.27
43.45, −87.31470.41
Pisgah (35.2°,
−82.85°)
1127.11.389733.11, −87.1732.88,
−85.66
516.91493.28898.091001.83
32.64, −84.14469.65
1127.11.437733.67, −85.7333.16,
−84.94
516.91493.33898.09991.8
32.64, −84.14469.74
1127.10421.393835.03, −87.1733.84,
−85.66
516.91493.43898.36965.39
32.64, −84.14469.94
1127.10421.393233.11, −87.1733.83,
−85.67
517.18493.42898.361044.64
34.54, −84.16469.65
1127.10421.441233.67, −85.7334.11,
−84.95
517.18493.46898.361033.88
34.54, −84.16469.74
1127.10831.397335.03, −87.1734.79,
−85.67
517.18493.56898.291005.55
34.54, −84.16469.94
1127.10831.445335.59, −85.7335.07,
−84.95
517.18493.61898.29994.54
34.54, −84.16470.03
1127.11251.401836.95, −87.1635.75,
−85.66
517.18493.71898.22968.96
34.54, −84.16470.24
1127.11251.401235.03, −87.1735.74,
−85.67
517.48493.71898.221048.81
36.44, −84.17469.94
1127.11251.449235.59, −85.7336.02,
−84.95
517.48493.76898.221037
36.44, −84.17470.03
1127.11671.405736.95, −87.1636.7,
−85.67
517.48493.86898.151009.58
36.44, −84.17470.24
1127.11671.453737.51, −85.7236.98,
−84.95
517.48493.91898.15997.47
36.44, −84.17470.34
1127.12081.410538.87, −87.1437.66,
−85.66
517.48494.03898.08972.57
36.44, −84.17470.57
1127.12081.409536.95, −87.1637.65,
−85.67
517.79494.02898.081053.3
38.34, −84.18470.24
1433.845822.733531.54, −79.1132.79,
−76.69
517.06492.62893.311076.3
34.03, −74.26468.17
1433.845822.734233.46, −79.1132.8,
−76.68
516.76492.61893.311025.52
32.13, −74.24468.45
1433.8522.689732.82, −80.5533.43,
−77.41
517.06492.71893.051064.71
34.03, −74.26468.35
1433.8522.737733.46, −79.1133.75,
−76.69
517.06492.76893.051049.78
34.03, −74.26468.45
1433.854222.693834.74, −80.5534.39,
−77.41
517.06492.86892.781038.25
34.03, −74.26468.65
1433.854222.741133.46, −79.1134.7,
−76.7
517.39492.92892.781075.51
35.93, −74.28468.45
1433.854222.742235.38, −79.134.71,
−76.68
517.06492.91892.781023.8
34.03, −74.26468.76
1433.858322.697334.74, −80.5535.34,
−77.42
517.39493.02891.821063.6
35.93, −74.28468.65
1433.858322.745635.38, −79.135.66,
−76.69
517.39493.08891.821048.6
35.93, −74.28468.76
1433.862522.701836.66, −80.5436.3,
−77.41
517.39493.18891.571037.24
35.93, −74.28468.97
1433.862522.749535.38, −79.136.61,
−76.7
517.73493.25891.571074.14
37.83, −74.29468.76
1433.862522.750137.3, −79.0936.62,
−76.69
517.39493.24891.571022.68
35.93, −74.28469.08
1433.866722.705736.66, −80.5437.25,
−77.42
517.73493.35891.321062.43
37.83, −74.29468.97
1433.866722.75437.3, −79.0937.57,
−76.69
517.73493.41891.321047.32
37.83, −74.29469.08
1682.045820.506933.75, −86.7433.15,
−83.09
517.36492.91000.05883.9
32.54, −79.43468.43
1682.045820.554532.59, −85.2933.52,
−82.37
517.65492.971000.05939.43
34.44, −79.45468.28
1682.045820.555134.51, −85.2933.53,
−82.36
517.36492.971000.05898.59
32.54, −79.43468.57
1682.0520.510433.75, −86.7434.1,
−83.1
517.65493.04999.8905.9
34.44, −79.45468.43
1682.0520.558734.51, −85.2934.48,
−82.37
517.65493.11999.8921.31
34.44, −79.45468.57
1682.054220.514535.67, −86.7435.06,
−83.1
517.65493.2999.22888.8
34.44, −79.45468.74
1682.054220.562534.51, −85.2935.43,
−82.38
517.96493.27999.22944.58
36.34, −79.46468.57
1682.054220.563236.43, −85.2835.44,
−82.37
517.65493.27999.22903.84
34.44, −79.45468.89
1682.058320.518335.67, −86.7436.01,
−83.1
517.96493.35998.98910.72
36.34, −79.46468.74
1682.058320.56736.43, −85.2836.39,
−82.37
517.96493.43998.98926.49
36.34, −79.46468.89
1682.062520.523237.59, −86.7236.97,
−83.09
517.96493.52998.73893.85
36.34, −79.46469.07
1682.062520.570836.43, −85.2837.34,
−82.38
518.28493.59998.73949.9
38.24, −79.47468.89
1682.062520.571538.35, −85.2737.35,
−82.37
517.96493.59998.73909.78
36.34, −79.46469.22
Urbana (40.13°,
−88.2°)
1127.12081.410538.87, −87.1437.66,
−85.66
517.48494.03994.39972.57
36.44, −84.17470.57
1127.12081.409536.95, −87.1637.65,
−85.67
517.79494.02994.391053.3
38.34, −84.18470.24
1127.12081.457537.51, −85.7237.93,
−84.95
517.79494.07994.391040.3
38.34, −84.18470.34
1127.1251.414338.87, −87.1438.61,
−85.66
517.79494.18993.781013.64
38.34, −84.18470.57
1127.1251.462339.43, −85.738.89,
−84.94
517.79494.23993.781000.54
38.34, −84.18470.67
1127.12921.419240.78, −87.1239.56,
−85.65
517.79494.35993.74977.73
38.34, −84.18470.9
1127.12921.418538.87, −87.1439.56,
−85.66
518.12494.35993.741056.14
40.24, −84.18470.57
1127.12921.466539.43, −85.739.84,
−84.94
518.12494.4993.741042.11
40.24, −84.18470.67
1127.13331.423340.78, −87.1240.51
−85.65
518.12494.51993.71017.76
40.24, −84.18470.9
1127.13331.471741.34, −85.6740.79,
−84.93
518.12494.57993.71002.4
40.24, −84.18471.01
1127.13751.428842.7, −87.0841.47,
−85.63
518.12494.69993.66982.72
40.24, −84.18471.25
1127.13751.427840.78, −87.1241.46,
−85.65
518.47494.69993.661059.24
42.14, −84.17470.9
1127.13751.476241.34, −85.6741.74,
−84.92
518.47494.74993.661042.83
42.14, −84.17471.01
1127.14171.433342.7, −87.0842.42
−85.63
518.47494.86993.461021.89
42.14, −84.17471.25
1375.395823.275538.97, −93.4937.8,
−91.81
517.17493.29996.41018.72
36.63, −90.12469.4
1375.395823.274937.05, −93.537.79,
−91.82
517.52493.29996.41095.5
38.53, −90.13469.06
1375.423.231438.36, −94.9338.45,
−92.53
517.52493.42996.241055.53
38.53, −90.13469.31
1375.423.279438.97, −93.4938.75,
−91.81
517.52493.46996.241056.11
38.53, −90.13469.4
1375.404223.236240.28, −94.9139.41,
−92.52
517.52493.59996.011021.46
38.53, −90.13469.66
1375.404223.235538.36, −94.9339.4,
−92.53
517.89493.6996.011094.6
40.43, −90.13469.31
1375.404223.284540.88, −93.4639.71,
−91.8
517.52493.64996.011019.53
38.53, −90.13469.75
1375.404223.283538.97, −93.4939.7,
−91.81
517.89493.65996.011095.3
40.43, −90.13469.4
1375.408323.240440.28, −94.9140.36,
−92.52
517.89493.78995.851058.55
40.43, −90.13469.66
1375.408323.288740.88, −93.4640.66,
−91.8
517.89493.82995.851056.57
40.43, −90.13469.75
1375.412523.245542.2, −94.8841.32,
−92.51
517.89493.96995.71025.08
40.43, −90.13470.02
1375.412523.244540.28, −94.9141.31,
−92.52
518.26493.96995.71096.58
42.33, −90.13469.66
1375.412523.293942.8, −93.4341.62,
−91.78
517.89494.01995.71019.87
40.43, −90.13470.12
1375.412523.292840.88, −93.4641.61,
−91.8
518.26494.01995.71094.54
42.33, −90.13469.75
1375.416723.249742.2, −94.8842.27,
−92.51
518.26494.14995.541061.25
42.33, −90.13470.02
1375.416723.29842.8, −93.4342.57,
−91.78
518.26494.19995.541055.78
42.33, −90.13470.12
1563.566721.571638.77, −92.3437.78,
−89.93
517.8493.57934.33874.53
36.78, −87.51469.33
1563.570821.575938.77, −92.3438.72,
−89.93
518.15493.74933.98894.65
38.67, −87.51469.33
1563.570821.624539.5, −90.8839.09,
−89.2
518.15493.8933.98897.76
38.67, −87.51469.45
1563.57521.58140.69, −92.3139.68,
−89.91
518.15493.92933.57877.68
38.67, −87.51469.69
1563.579221.585240.69, −92.3140.63,
−89.91
518.51494.1933.16897.64
40.57, −87.51469.69
1563.579221.633941.42, −90.8541,
−89.18
518.51494.17933.16900.69
40.57, −87.51469.82
1563.583321.590442.61, −92.2841.59,
−89.9
518.51494.29932.75881.13
40.57, −87.51470.06
1563.587521.594542.61, −92.2842.54,
−89.9
518.88494.47932.34900.72
42.47, −87.51470.06
Virginia (37.2°,
−80.42°)
1127.10831.397335.03, −87.1734.79,
−85.67
469.94493.56997.971005.55
34.54, −84.16493.56
1127.10831.445335.59, −85.7335.07,
−84.95
470.03493.61997.97994.54
34.54, −84.16493.61
1127.11251.401836.95, −87.1635.75,
−85.66
470.24493.71997.55968.96
34.54, −84.16493.71
1127.11251.401235.03, −87.1735.74,
−85.67
469.94493.71997.551048.81
36.44, −84.17493.71
1127.11251.449235.59, −85.7336.02,
−84.95
470.03493.76997.551037
36.44, −84.17493.76
1127.11671.405736.95, −87.1636.7,
−85.67
470.24493.86997.571009.58
36.44, −84.17493.86
1127.11671.453737.51, −85.7236.98,
−84.95
470.34493.91997.57997.47
36.44, −84.17493.91
1127.12081.410538.87, −87.1437.66,
−85.66
470.57494.03997.59972.57
36.44, −84.17494.03
1127.12081.409536.95, −87.1637.65,
−85.67
470.24494.02997.591053.3
38.34, −84.18494.02
1127.12081.457537.51, −85.7237.93,
−84.95
470.34494.07997.591040.3
38.34, −84.18494.07
1127.1251.414338.87, −87.1438.61,
−85.66
470.57494.18998.091013.64
38.34, −84.18494.18
1127.1251.462339.43, −85.738.89,
−84.94
470.67494.23998.091000.54
38.34, −84.18494.23
1127.12921.419240.78, −87.1239.56,
−85.65
470.9494.35998.11977.73
38.34, −84.18494.35
1127.12921.418538.87, −87.1439.56,
−85.66
470.57494.35998.111056.14
40.24, −84.18494.35
1433.858322.697334.74, −80.5535.34,
−77.42
468.65493.02984.751063.6
35.93, −74.28493.02
1433.858322.745635.38, −79.135.66,
−76.69
468.76493.08984.751048.6
35.93, −74.28493.08
1433.862522.701836.66, −80.5436.3,
−77.41
468.97493.18984.661037.24
35.93, −74.28493.18
1433.862522.749535.38, −79.136.61,
−76.7
468.76493.25984.661074.14
37.83, −74.29493.25
1433.862522.750137.3, −79.0936.62,
−76.69
469.08493.24984.661022.68
35.93, −74.28493.24
1433.866722.705736.66, −80.5437.25,
−77.42
468.97493.35984.571062.43
37.83, −74.29493.35
1433.866722.75437.3, −79.0937.57,
−76.69
469.08493.41984.571047.32
37.83, −74.29493.41
1433.870822.710538.58, −80.5238.21,
−77.41
469.31493.52984.481035.98
37.83, −74.29493.52
1433.870822.758237.3, −79.0938.52,
−76.69
469.08493.59984.481072.23
39.73, −74.29493.59
1433.870822.758539.22, −79.0838.53,
−76.69
469.43493.58984.481021.35
37.83, −74.29493.58
1433.87522.714738.58, −80.5239.16,
−77.41
469.31493.7984.391060.56
39.73, −74.29493.7
1433.87522.762739.22, −79.0839.48,
−76.69
469.43493.76984.391045.4
39.73, −74.29493.76
1622.016721.085734.53, −76.3534.92,
−73.97
468.71493.241048.98980.04
35.3, −71.58493.24
1622.020821.041835.67, −77.7935.49,
−74.69
468.87493.321048.52957.52
35.3, −71.58493.32
1622.020821.090236.45, −76.3435.88,
−73.96
469.03493.41048.52958.83
35.3, −71.58493.4
1622.02521.046337.58, −77.7836.44,
−74.68
469.2493.481048.2937.99
35.3, −71.58493.48
1622.02521.045735.67, −77.7936.44,
−74.69
468.87493.481048.2982.09
37.2, −71.59493.48
1622.02521.09436.45, −76.3436.83,
−73.97
469.03493.561048.2983.54
37.2, −71.59493.56
1622.029221.050237.58, −77.7837.39,
−74.69
469.2493.641047.88961.85
37.2, −71.59493.64
1622.029221.098938.36, −76.3237.78,
−73.96
469.36493.721047.88962.74
37.2, −71.59493.72
1622.033321.05539.5, −77.7638.35,
−74.68
469.55493.821047.56942.48
37.2, −71.59493.82
1622.033321.05437.58, −77.7838.34,
−74.69
469.2493.811047.56986.25
39.1, −71.6493.81
1622.033321.102738.36, −76.3238.73,
−73.96
469.36493.891047.56987.27
39.1, −71.6493.89
1622.037521.058939.5, −77.7639.3,
−74.68
469.55493.991047.25966.22
39.1, −71.6493.99
1622.037521.107540.28, −76.339.69,
−73.95
469.71494.071047.25966.66
39.1, −71.6494.07
1682.054220.514535.67, −86.7435.06,
−83.1
468.74493.21023.73888.8
34.44, −79.45493.2
1682.054220.562534.51, −85.2935.43,
−82.38
468.57493.271023.73944.58
36.34, −79.46493.27
1682.054220.563236.43, −85.2835.44,
−82.37
468.89493.271023.73903.84
34.44, −79.45493.27
1682.058320.518335.67, −86.7436.01,
−83.1
468.74493.351023.57910.72
36.34, −79.46493.35
1682.058320.56736.43, −85.2836.39,
−82.37
468.89493.431023.57926.49
36.34, −79.46493.43
1682.062520.523237.59, −86.7236.97,
−83.09
469.07493.521023.42893.85
36.34, −79.46493.52
1682.062520.570836.43, −85.2837.34,
−82.38
468.89493.591023.42949.9
38.24, −79.47493.59
1682.062520.571538.35, −85.2737.35,
−82.37
469.22493.591023.42909.78
36.34, −79.46493.59
1682.066720.52737.59, −86.7237.92,
−83.1
469.07493.681023.26915.91
38.24, −79.47493.68
1682.066720.575338.35, −85.2738.3,
−82.37
469.22493.751023.26932.6
38.24, −79.47493.75
1682.070820.531839.51, −86.738.88,
−83.09
469.41493.851022.9898.69
38.24, −79.47493.85
1682.070820.579538.35, −85.2739.25,
−82.37
469.22493.921022.9955.56
40.14, −79.47493.92
1682.070820.580240.27, −85.2539.26,
−82.36
469.57493.931022.9915.65
38.24, −79.47493.93
Xing Long (40.2°,
117.4°)
11816.63330.897338.19, 123.2337.88,
123.96
517.47493.81959.93922.82
37.57, 124.69470.15
11816.63750.854239.55, 121.8138.56,
123.25
517.47493.93959.84895.31
37.57, 124.69470.39
11816.63750.853537.63, 121.7938.55,
123.24
517.8493.93959.84959.91
39.47, 124.69470.06
11816.63750.901538.19, 123.2338.83,
123.96
517.8493.98959.84958.39
39.47, 124.69470.15
11816.64170.858339.55, 121.8139.51,
123.25
517.8494.1959.74929.13
39.47, 124.69470.39
11816.64170.906740.11, 123.2639.79,
123.98
517.8494.14959.74927.44
39.47, 124.69470.48
11816.64580.863541.47, 121.8440.47,
123.27
517.8494.27959.65900.2
39.47, 124.69470.73
11816.64580.862539.55, 121.8140.46,
123.25
518.15494.27959.65964.76
41.37, 124.69470.39
11816.64580.910840.11, 123.2640.74,
123.98
518.15494.32959.65963.02
41.37, 124.69470.48
11816.650.867741.47, 121.8441.42,
123.27
518.15494.44959.56934.08
41.37, 124.69470.73
11816.650.915742.03, 123.2841.7,
123.99
518.15494.49959.56932.15
41.37, 124.69470.82
11816.65420.873243.38, 121.8842.38,
123.29
518.15494.62959.47905.01
41.37, 124.69471.09
11816.65420.872241.47, 121.8442.37,
123.27
518.51494.62959.47969.58
43.26, 124.7470.73
11816.65420.920242.03, 123.2842.65,
123.99
518.51494.67959.47967.59
43.26, 124.7470.82
14314.887522.714938.39, 115.8737.74,
117.41
517.61493.471118.16900.95
37.08, 118.95469.32
14314.891722.670737.76, 114.4338.37,
116.69
517.96493.591117.94927.64
38.98, 118.94469.22
14314.891722.718638.39, 115.8738.69,
117.41
517.96493.641117.94930.77
38.98, 118.94469.32
14314.895822.675239.68, 114.4439.33,
116.69
517.96493.771117.73899.96
38.98, 118.94469.57
14314.895822.723540.31, 115.8939.65,
117.42
517.96493.821117.73902.19
38.98, 118.94469.68
14314.895822.722938.39, 115.8739.64,
117.41
518.33493.831117.73961.69
40.88, 118.94469.32
14314.922.679439.68, 114.4440.28,
116.69
518.33493.951117.51929.23
40.88, 118.94469.57
14314.922.727740.31, 115.8940.6,
117.42
518.33494.011117.51931.65
40.88, 118.94469.68
14314.904222.684541.59, 114.4741.24,
116.71
518.33494.131117.3902.39
40.88, 118.94469.93
14314.904222.732942.23, 115.9241.56,
117.43
518.33494.191117.3903.65
40.88, 118.94470.04
14314.904222.732240.31, 115.8941.54,
117.42
518.71494.21117.3961.81
42.77, 118.95469.68
14314.908322.68941.59, 114.4742.18,
116.71
518.71494.321117.09931.03
42.77, 118.95469.93
14314.908322.737442.23, 115.9242.5,
117.44
518.71494.381117.09932.43
42.77, 118.95470.04

References

  1. Burke, W.J.; Lin, C.S.; Hagan, M.P.; Huang, C.Y.; Weimer, D.R.; Wise, J.O.; Gentile, L.C.; Marcos, F.A. Storm time global thermosphere: A driven-dissipative thermodynamic system. J. Geophys. Res.-Space Phys. 2009, 114, A06306. [Google Scholar] [CrossRef]
  2. Wise, J.O.; Burke, W.J.; Sutton, E.K. Globally averaged exospheric temperatures derived from CHAMP and GRACE accelerometer measurements. J. Geophys. Res.-Space Phys. 2012, 117. [Google Scholar] [CrossRef]
  3. Picone, J.M.; Hedin, A.E.; Drob, D.P.; Aikin, A.C. Nrlmsise-00 Empirical Model of the Atmosphere: Statistical Comparisons and Scientific Issues. J. Geophys. Res.-Space Phys. 2002, 107, SIA 15-1–SIA 15-16. [Google Scholar] [CrossRef]
  4. Forbes, J.M.; Zhang, X.; Bruinsma, S.L. New perspectives on thermosphere tides: 2. Penetration to the upper thermosphere. Earth Planets Space 2014, 66, 1–11. [Google Scholar] [CrossRef] [Green Version]
  5. Weimer, D.R.; Sutton, E.K.; Mlynczak, M.G.; Hunt, L.A. Intercalibration of neutral density measurements for mapping the thermosphere. J. Geophys. Res. Space Phys. 2016, 121, 5975–5990. [Google Scholar] [CrossRef] [Green Version]
  6. Weng, L.; Lei, J.; Sutton, E.; Dou, X.; Fang, H. An exospheric temperature model from CHAMP thermospheric density. Space Weather 2017, 15, 343–351. [Google Scholar] [CrossRef]
  7. Evans, J.S.; Lumpe, J.D.; Correira, J.; Veibell, V.; Kyrwonos, A.; McClintock, W.E.; Solomon, S.C.; Eastes, R.W. Neutral Exospheric Temperatures From the GOLD Mission. J. Geophys. Res.-Space Phys. 2020, 125, e2020JA027814. [Google Scholar] [CrossRef]
  8. Friis-Christensen, E.; Luhr, H.; Knudsen, D.; Haagmans, R. Swarm—An Earth Observation Mission investigating Geospace. Adv. Space Res. 2008, 41, 210–216. [Google Scholar] [CrossRef]
  9. Doornbos, E.; van den Ijssel, J.; Luhr, H.; Forster, M.; Koppenwallner, G. Neutral Density and Crosswind Determination from Arbitrarily Oriented Multiaxis Accelerometers on Satellites. J. Spacecr. Rockets 2010, 47, 580–589. [Google Scholar] [CrossRef] [Green Version]
  10. Van den Ijssel, J.; Doornbos, E.; Iorfida, E.; March, G.; Siemes, C.; Montenbruck, O. Thermosphere densities derived from Swarm GPS observations. Adv. Space Res. 2020, 65, 1758–1771. [Google Scholar] [CrossRef]
  11. Yuan, W.; Xu, J.; Ma, R.; Wu, Q.; Jiang, G.; Gao, H.; Liu, X.; Chen, S. First observation of mesospheric and thermospheric winds by a Fabry-Perot interferometer in China. Chin. Sci. Bull. 2010, 55, 4046–4051. [Google Scholar] [CrossRef]
  12. Yuan, W.; Liu, X.; Xu, J.; Zhou, Q.; Jiang, G.; Ma, R. FPI observations of nighttime mesospheric and thermospheric winds in China and their comparisons with HWM07. Ann. Geophys. 2013, 31, 1365–1378. [Google Scholar] [CrossRef] [Green Version]
  13. Chen, G.-m.; Xu, J.; Wang, W.; Burns, A.G. A comparison of the effects of CIR- and CME-induced geomagnetic activity on thermospheric densities and spacecraft orbits: Statistical studies. J. Geophys. Res. Space Phys. 2014, 119, 7928–7939. [Google Scholar] [CrossRef]
  14. Cai, X.; Burns, A.G.; Wang, W.; Qian, L.; Solomon, S.C.; Eastes, R.W.; Pedatella, N.; Daniell, R.E.; McClintock, W.E. The Two-Dimensional Evolution of Thermospheric ΣO/N2 Response to Weak Geomagnetic Activity During Solar-Minimum Observed by GOLD. Geophys. Res. Lett. 2020, 47, e2020GL088838. [Google Scholar] [CrossRef]
  15. Buonsanto, M.J.; Pohlman, L.M. Climatology of neutral exospheric temperature above Millstone Hill. J. Geophys. Res.-Space Phys. 1998, 103, 23381–23392. [Google Scholar] [CrossRef]
Figure 1. Daily mean altitudes and neutral densities of the Swarm satellites (The colors of Swarm A, Swarm B, Swarm C are pink, green and red. The densities of Swarm A and Swarm C overlap because the two satellites are flying closely together).
Figure 1. Daily mean altitudes and neutral densities of the Swarm satellites (The colors of Swarm A, Swarm B, Swarm C are pink, green and red. The densities of Swarm A and Swarm C overlap because the two satellites are flying closely together).
Remotesensing 14 05382 g001
Figure 2. Variation of daytime and nighttime neutral temperatures with altitudes in 2008 and 2014. The green dots are the temperature of Swarm B, the red dots are the temperature of Swarm A and Swarm C, the black triangles are the temperatures at FPIs’ altitude.
Figure 2. Variation of daytime and nighttime neutral temperatures with altitudes in 2008 and 2014. The green dots are the temperature of Swarm B, the red dots are the temperature of Swarm A and Swarm C, the black triangles are the temperatures at FPIs’ altitude.
Remotesensing 14 05382 g002
Figure 3. Time distribution of inversion temperatures from January 2014 to October 2014, the result number is the number of points in the low and middle latitudes.
Figure 3. Time distribution of inversion temperatures from January 2014 to October 2014, the result number is the number of points in the low and middle latitudes.
Remotesensing 14 05382 g003
Figure 4. The diurnal variation of derived temperatures (K) from 15 April 2014 to 11 August 2014.
Figure 4. The diurnal variation of derived temperatures (K) from 15 April 2014 to 11 August 2014.
Remotesensing 14 05382 g004
Figure 5. Diurnal variation of derived temperatures calculated with m ¯ s , model temperatures and ground-based FPI temperatures. The duration is from mid-April to early August in 2014. Black dots represent the derived temperatures, red dots represent the MSIS temperatures, red line represents the fitted curve of red dots; and black solid triangles represent the ground-based FPI temperatures during this time. The left side represents the four latitude bands of the Southern Hemisphere, and the right represents the Northern Hemisphere.
Figure 5. Diurnal variation of derived temperatures calculated with m ¯ s , model temperatures and ground-based FPI temperatures. The duration is from mid-April to early August in 2014. Black dots represent the derived temperatures, red dots represent the MSIS temperatures, red line represents the fitted curve of red dots; and black solid triangles represent the ground-based FPI temperatures during this time. The left side represents the four latitude bands of the Southern Hemisphere, and the right represents the Northern Hemisphere.
Remotesensing 14 05382 g005
Figure 6. Diurnal variation of derived temperatures calculated with m ¯ 0 , model temperatures and ground-based FPI temperatures. The duration is from mid-April to early August in 2014. Black dots represent the derived temperatures, red dots represent the MSIS temperatures, red line represents the fitted curve of red dots; and black solid triangles represent the ground-based FPI temperatures during this time. The left side represents the four latitude bands of the Southern Hemisphere, and the right represents the Northern Hemisphere.
Figure 6. Diurnal variation of derived temperatures calculated with m ¯ 0 , model temperatures and ground-based FPI temperatures. The duration is from mid-April to early August in 2014. Black dots represent the derived temperatures, red dots represent the MSIS temperatures, red line represents the fitted curve of red dots; and black solid triangles represent the ground-based FPI temperatures during this time. The left side represents the four latitude bands of the Southern Hemisphere, and the right represents the Northern Hemisphere.
Remotesensing 14 05382 g006
Figure 7. The diurnal variation of mean molecular mass m ¯ s (g/mol) from 15 April 2014 to 11 August 2014 (The same calculation as Figure 4).
Figure 7. The diurnal variation of mean molecular mass m ¯ s (g/mol) from 15 April 2014 to 11 August 2014 (The same calculation as Figure 4).
Remotesensing 14 05382 g007
Figure 8. The schematic diagram of Satellite tethered Mass Spectrometer Detection.
Figure 8. The schematic diagram of Satellite tethered Mass Spectrometer Detection.
Remotesensing 14 05382 g008
Figure 9. The relationship between the relative errors, wire length and densities’ measurement error. The horizontal axis indicates the wire length, the vertical axis indicates the relative error of ( ρ 1 ρ 0 ), and the color indicates the error of the derived temperatures. The dashed line shows the different error of derived temperatures.
Figure 9. The relationship between the relative errors, wire length and densities’ measurement error. The horizontal axis indicates the wire length, the vertical axis indicates the relative error of ( ρ 1 ρ 0 ), and the color indicates the error of the derived temperatures. The dashed line shows the different error of derived temperatures.
Remotesensing 14 05382 g009
Table 1. Geographical locations of FPI stations and their number of selected observations. Positive latitude and longitude represent north and east.
Table 1. Geographical locations of FPI stations and their number of selected observations. Positive latitude and longitude represent north and east.
FPI StationLatitudeLongitudeUsed Observations in the Comparison
Arecibo18.34°−66.75°14
Pisgah35.2°−82.85°41
Virginia37.2°−80.42°52
Eastern Kentucky37.75°−84.29°27
Urbana40.13°−88.2°38
Xinglong40.2°117.4°27
Peach mountain42.27°−83.75°4
Table 2. Comparative statistical analysis of inversion temperatures calculated with m ¯ 0 and m ¯ s , FPI temperatures and MSIS simulations.
Table 2. Comparative statistical analysis of inversion temperatures calculated with m ¯ 0 and m ¯ s , FPI temperatures and MSIS simulations.
Temperatures for
Comparison
Sample Points Mean Deviation, KStandard Deviation, K
T ( m ¯ s ) & T(FPI)20330.80106.20
T ( m ¯ 0 ) & T(FPI)203−30.67108.03
T (MSIS at 500 km)
& T (MSIS at 250 km)
20323.5214.53
T ( m ¯ s ) & T(MSIS)8101116.5284.46
T ( m ¯ 0 ) & T(MSIS)8101145.90140.73
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lin, Y.; Xu, J.; van den IJssel, J.; Yuan, W. A Novel Method to Derive Exospheric Temperatures from Swarm Thermospheric Densities during Quiet Times. Remote Sens. 2022, 14, 5382. https://0-doi-org.brum.beds.ac.uk/10.3390/rs14215382

AMA Style

Lin Y, Xu J, van den IJssel J, Yuan W. A Novel Method to Derive Exospheric Temperatures from Swarm Thermospheric Densities during Quiet Times. Remote Sensing. 2022; 14(21):5382. https://0-doi-org.brum.beds.ac.uk/10.3390/rs14215382

Chicago/Turabian Style

Lin, Yingyun, Jiyao Xu, Jose van den IJssel, and Wei Yuan. 2022. "A Novel Method to Derive Exospheric Temperatures from Swarm Thermospheric Densities during Quiet Times" Remote Sensing 14, no. 21: 5382. https://0-doi-org.brum.beds.ac.uk/10.3390/rs14215382

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop