1. Introduction
Comparison of atmospheric parameters observed with ground-based and satellite instruments is important to study the Earth’s atmosphere, especially the upper atmosphere in particular (mesosphere, lower thermosphere), where direct observations of temperature, composition, wind speed, etc. are hampered and quite rarely conducted using geophysical rockets. On the one hand, such comparisons allow verifying the methods of remote indirect observation of atmospheric parameters, while on the other hand, they offer additional opportunities for a more detailed investigation into physical and chemical processes in the atmosphere. For example, article [
1] describes air density correction in the NRLMSISE-00 empirical model. The correction is based on data obtained using the SABER instrument aboard the TIMED satellite. The performed correction results in far better accuracy of the model at 80–100 km. In [
2], the authors compared spatial distributions of the hydroxyl airglow intensity and air temperature obtained using the special all-sky imaging system at Ranchi station in India (23.3°N, 85.3°E) and temperatures obtained with SABER above this location. The results of the experiment are consistent, and the dispersions of the values obtained are within admissible limits. The article [
3] compares the mesopause temperature at the height of hydroxyl emission detected with the ground-based infrared spectrograph at the mid-latitude Maimaga station in Yakutsk and the mesopause temperature obtained with the SABER satellite in the same region. In [
4], too, comparative analysis of temperature was performed using the high-latitude ground-based spectrometer SATI (Spectral Airglow Temperature Imager) in Sierra Nevada, Spain, and the SABER satellite instrument. In turn, ref. [
5] compared several instruments and calculated shifts in temperatures measured with the ground-based Davis spectrometer in the Antarctic and those obtained from the Aura/MLS and TIMED/SABER satellite data. The article [
6] discusses the optimal choice of Einstein coefficients to determine the rotational temperature using a spectrometer (Horiba model: IHR550) in Xinglong China (40°24′N, 117°35′E) and SABER data. Conclusions put forth in all the listed works denote fairly good correspondence between data obtained with satellite and ground-based instrumentation. The airglow of the 557.7 nm atomic oxygen line occurs near the mesopause region, a little higher than the hydroxyl airglow. There exists a potential for temperature and wind observations by analyzing the 557.7 nm line Doppler broadening and shifting using aeronomic Fabry–Pérot interferometers that operate in high latitudes [
7,
8]. However, no comparison has been conducted with satellite data for such devices. This is important both to check the performance of the method for ground-based observation of the Earth’s upper atmosphere parameters and to refine features of dynamics of the upper atmosphere parameters at the mesopause at the mid-latitudes of Eurasia. Ground-based instruments have a better time resolution in comparison with the satellite-borne devices. Therefore, the importance of the ground-based instruments for studying the tides internal gravity waves and other processes that require time resolution in minutes and tens of minutes is obvious. Thus, in our study, we assessed the extent of compliance between the data from ground-based (Fabry-Pérot interferometer) and satellite (limb-scanning radiometer) measurements of the upper atmosphere temperature. The intensity of the 557.7 nm atmospheric airglow line recorded with ground-based equipment resulted from the atomic oxygen
–
) transition at the height of 90–100 km. Height profiles of the upper atmosphere physical and chemical parameters (including temperature) were considered in the same range of heights using satellite limb-scanning radiometry.
3. Results and Discussion
With the atomic oxygen concentration known, it is possible to calculate its intensity 557.7 nm (hereinafter 558) using the formula from [
21]:
where
M is the atmosphere density,
is the coefficient of three bodies recombination rate,
is the coefficient of
extinction with oxygen,
,
are the Einstein coefficients. Values of these coefficients are presented in [
21].
Figure 6 shows the airglow profiles of the 557.7 nm atomic oxygen line averaged over the entire time interval considered. The profiles were synthesized from calculated oxygen profiles data with different sets of coefficients. Here, we can note that the intensity obtained with new coefficients [
24] has grown but the maximum of the airglow profile (red line) stays in the same place as the old set of coefficients [
20] (green line). Similarly, the oxygen concentration intensity increases in the 93–100 km range for the new coefficients. The variance in the red curve due to a lack of data apparently has increased. The sharp increase in airglow intensity above 100 km is conditioned by the incorrectly determined atomic oxygen concentration at these heights, as was described above. The 557.7 nm airglow profile was obtained using oxygen from [
25], plotted using a magenta line. One can see a very narrow peak at 97 km and then a sharp increase in the airglow intensity starting from 98 km. It is rather surprising to find a sharp resolved peak in the vertical distribution of airglow because the calculations of airglow in this case were based on the monotonically increasing atomic oxygen concentration. The peak position near the height where the methods become unstable point to some possible errors in calculations due to variable data. Nevertheless, we took this peak into account for further calculations.
It can be seen that the airglow peak of the atomic oxygen 557.7 nm line is at a height of about 94–95 km. This indicates that the calculations we performed are correct in the first approximation and that the result is in agreement with the previous study [
21]. In
Figure 7, the temporal behavior of the 557.7 nm airglow presented in the same way as the temporal behavior of the atomic oxygen concentration in
Figure 5 is shown. The top panel of
Figure 7 presents data calculated using oxygen by [
20], the middle panel presents data calculated using oxygen by [
24], and bottom panel presents data calculated using oxygen by [
25]. The main difference is the significant increase in the intensity at summer for the top and middle panels in comparison with absence of such significant seasonal intensification of 557.7 nm airglow in the bottom panel.
Worth noting is the fact that the weak variation in the airglow maximum is in antiphase with the height variation and temperature of the mesopause. The mechanism of triple collisions responsible for the formation of atomic oxygen in the state is inversely dependent on temperature. In this regard, the airglow maximum should follow the mesopause and be of increased intensity in summer.
We should integrate the intensity over 80–102 km to get the correct values intensity in order to compare the ground-based and satellite observations. The 557.7 nm emission-integrated intensity obtained from SABER data and the interferometer-observed emission show a good match, especially in winter and springtime (
Figure 8). The temporal behavior of integrated intensities in the bottom panel of
Figure 8 already exhibit the mentioned absence of a summer increase in the 557.7 nm intensity obtained using more clear data from [
25].
Here the 557.7 nm airglow intensity observed by FPI is in arbitrary units. We did not yet calibrate the intensity to bring it to the photons flux values, but we took into account the background. Thus, the arbitrary units intensity is the 557.7 nm line intensity, which should coincide rather well with the real 557.7 nm photons flux.
The height profile of the atomic oxygen airglow intensity reconstructed from SABER data allows obtaining the effective temperature using Equation (1) and comparing it with the results of ground-based observations.
Figure 9 shows the behavior of these parameters throughout observations above the GPO in Tory, where the top panel presents data calculated using [
20], the middle panel presents data calculated using [
24], and the bottom panel presents data calculated using [
25].
The change in the observed characteristics over time (
Figure 9) shows that the effective temperature variations from SABER data have decreased their variance compared with the behavior of the average temperature in
Figure 2, but they do not have a pronounced seasonal change, especially in the top and middle panels’ data. The FPI-observed temperature has seasonal variations, which in the first approximation corresponds to the dynamics of seasonal temperature in the mesopause region. Significant differences using the method of reconstruction of the atomic oxygen profiles using [
20,
24] are not seen, although there are some gaps in the data when using [
24]. The temperatures in the bottom panel are rather similar, but still there are some discrepancies in summer and autumn time, where the intensities also differ (see
Figure 8).
As it was mentioned, seasonal variation in the height of the 557.7 nm intensity maximum does not correlate with the mesopause position. In this regard, let us assume that the possible reason for different seasonal variations in temperatures obtained using different methods (
Figure 9) is due to our incorrect reconstruction of the
airglow height profile using the proposed methods and SABER data. The reason can be in the incorrect height of the airglow intensity maximum. To briefly check the assumption, we used a simple technical procedure similar to the one used in [
11]. By shifting vertically the airglow height profile obtained from SABER data, we calculated a new effective temperature using Equation (1) for every shift like this. To find the most appropriate new height of the airglow maximum position, we used the error weighting function that was calculated over the entire period of observation:
where
is the atomic oxygen intensity integrated in height,
—atmosphere temperature measured with Fabry–Pérot interferometer,
—atmosphere effective temperature determined from SABER/TIMED data for the profile whose maximum is the height
, index
is the point of observation time. The height corresponding to the minimum value (Equation (8)) means that for a given position of glow profile, temporary variations in FPI-measured temperature will demonstrate the most exact coincidence with temporary variations in SABER-detected effective temperature.
Figure 10 shows the error function (Equation (8)) calculated for different height shifts of the profile of atomic oxygen airglow. The error function was calculated for profiles obtained by three different methods.
The minimum of the error weighting function for the blue line falls on the height shift of +2 km (~97 km) and for the red line on the height of +3 km (98 km). The height shift obtained for the magenta line is −7 km (90 km). Such a difference apparently appears from the absence of pronounced seasonal variation in the 557.7 nm intensity obtained using [
25] because the essence of the described Equation (8) procedure is a minimization of the difference in seasonal variations of two raw data sets.
Seasonal behavior for the effective temperature obtained using shifted intensity height profiles matches better the behavior of the interferometer-measured temperature. In the top and middle panels of
Figure 11, we put the effective temperatures calculated from airglow profiles arising from [
20,
24] and shifted the values corresponding to their error functions’ minima (
Figure 10). One can see that the only deviations in the seasonal temperature variations for FPI and SABER appear in the middle of summer and autumn. The FPI temperature here is lower than the temperature obtained using SABER data. The integrated intensities for FPI and SABER also differ at their periods (
Figure 8). Intensity from FPI is higher than intensity from SABER data. Apparently this lack of additional intensity in SABER data results in the differences in the temperatures observed.
The variation in the effective temperature obtained with the shifted airglow profile by [
25] is (bottom panel of
Figure 11) increased, and the seasonal variations became less pronounced. It can be due to the above-mentioned difference in the shifting direction of the airglow profile. Apparently this is due to the above-mentioned absence of pronounced seasonal variations in the airglow intensity obtained by [
25]. One needs to mention the fact that without shifting the effective temperature obtained using the 557.7 airglow profile by [
25], the result better fits the temperature obtained by FPI (
Figure 9). The airglow profile in this case concentrates the intensity at higher altitudes in comparison with airglow profiles by [
20,
24]. Therefore, the reason for the similarity in temperatures in this case can again be the intensity profile placed at higher altitudes for effective temperature calculation.
Because it is complicated to see the details of the data in
Figure 8 and
Figure 9,
Figure 12 demonstrates the results using coefficients from [
20] because the amount of data is larger than using coefficients from [
24]. It can be mentioned that in winter–spring and in the first half of the summer period, we have better correspondences between intensity and temperature than in the last half of summer and all of the autumn period. However, should the winter period be considered in more detail, as an example, one can note that the interferometer-observed behavior of the 557.7 nm intensity coincides with the that of the intensity from SABER data, while the temperature behavior is still somewhat inconsistent (
Figure 12).
The obtained better agreement of seasonal temperature variations when the airglow profile synthesized from SABER data is shifted about 2–3 km up for both versions of calculated atomic oxygen airglow using [
20,
24] apparently means that the airglow profile determined is not quite correct, and that in fact, the 557.7 nm emission layer is higher. At least, it is fair for mid-latitudes where the FPI is located. In addition, it should be noted that local fast (non-seasonal) temperature variations obtained using the two instruments are in antiphase. Moreover, there is an inverse relationship between temperature and intensity for FPI, while for the SABER-synthesized 557.7 profile, these parameters are directly correlated [
12]. Note that one can observe significant differences in the behavior of MLT (mesosphere–lower thermosphere) temperature from satellite and ground-based data during sudden stratospheric warming events, when FPI data show a considerable temperature increase, and SABER data show its decrease [
12]. Thus, a possible reason for these variations can be the vertical dynamics of the 557.7 nm intensity profile, which occurs due to factors that cannot be taken into account based on SABER data. For example, a change in the atomic oxygen concentration, or in the effectiveness of the Barthes mechanism due to the arrival from the underlying atmosphere of chemical components not registered with SABER can lead both to a decrease in the number of airglow precursors and to suppression of emission itself. Additionally, one should not discard the possibility that along with the Barthes mechanism at the mesopause level, there is another mechanism generating 557.7 nm emission. The intensity of this assumed mechanism can be much lower than that of the Barthes mechanism, and the airglow profile maximum can be shifted higher. Thus, if the glow due to the Barthes mechanism will mainly be suppressed, one begins to observe parameters of the 557.7 nm line, whose glow is generated higher. This leads to an increase in the observed temperature when the glow intensity decreases. This happens because the temperature gradient above the mesopause is positive and considerably higher than underneath. As we can see from [
21], the dependence of coefficients governing the airglow is significant and inversely proportional to the temperature.
4. Conclusions
In our study we compared the intensities of the atmosphere night airglow 557.7 nm and the air temperatures at the mesopause level over local part of the Earth’s surface in Baikal natural territory. The data were obtained with the ground-based Fabry–Pérot interferometer and TIMED/SABER satellite instrument. Pre-comparison of the two data sets demonstrated good compliance of time variations in the mesopause temperature and in the temperature measured with the interferometer. However, taking into account the glow intensity height profile, which allows for correct comparison of data from ground- and satellite-based instruments, led to a discrepancy in the results: the effective integrated temperature calculated from the SABER parameter profiles did not reproduce the seasonal temperature variations detected with FPI. Additionally, variation in the position of the maximum of atomic oxygen airglow profile in time from SABER data did not follow the mesopause position. This was inconsistent with the recognized (Barthes) mechanism of generation of the 557.7 nm airglow, whose effectiveness becomes lower as the temperature increases.
We assumed that the emission layer height profile determined using simple aeronomical models and satellite data was incorrect—namely, that the position of the calculated layer airglow maximum was wrong. Therefore, shifting this synthesized profile of intensity, whose maximum was detected at the 94 km by using [
20] and 95 km by using [
24] coefficients, height relative to temperature, from initial SABER data, we re-calculated the effective temperatures. The best coincidence of temperature variations obtained with ground-based and satellite instruments was achieved when the intensity maximum was shifted to the height of 97 km by using [
20] and 98 km by using [
24] coefficients. The coefficients according to [
24] showed that the maximum atomic oxygen concentration and emission intensity were 1 km higher than the old version [
20]. But using [
24] we obtained less data on the concentration of atomic oxygen compared with [
20]. The same comparisons performed for 557.7 nm airglow profile obtained based on [
25] atomic oxygen data did not lead to the agreement in the temperatures obtained by FPI and SABER. The reason is the absence of distinct seasonal variations in atomic oxygen and in 557.7 nm airglow. However, the higher position of airglow obtained based on [
25] led to rather good agreement with the temperature from FPI and the effective temperature calculated without shifting of the airglow profile in comparison with the same values calculated based on [
20,
24].
The shift in the emission intensity maximum according to SABER data contributed to the almost correct reconstruction of seasonal temperature variations according to SABER. However, Fabry-Pérot and SABER data demonstrated discrepancies in short-term temperature variations, especially in summer, which suggests the existence of additional factors distorting the 557.7 nm line airglow profile and requires further study. In future research, one should take into account more complicated photochemical models for green airglow, where one should pay more attention to the height of the emission maximum or double-peak structure in the resulting calculations. We also do not discard other possibilities, such as a change in the width of the emission profile, for example.