Issue 
J. Space Weather Space Clim.
Volume 6, 2016



Article Number  A7  
Number of page(s)  11  
DOI  https://doi.org/10.1051/swsc/2015045  
Published online  05 February 2016 
Research Article
Validation of Earth atmosphere models using solar EUV observations from the CORONAS and PROBA2 satellites in occultation mode
^{1}
Lebedev Physical Institute of Russian Academy of Sciences, Moscow, Russia
^{2}
Institute for Physics of Microstructures of Russian Academy of Sciences, N. Novgorod, Russia
^{3}
Royal Observatory of Belgium, Brussels, Belgium
^{*} Corresponding author: slem@sci.lebedev.ru
Received:
12
February
2015
Accepted:
15
December
2015
Aims: Knowledge of properties of the Earth’s upper atmosphere is important for predicting the lifetime of loworbit spacecraft as well as for planning operation of space instruments whose data may be distorted by atmospheric effects. The accuracy of the models commonly used for simulating the structure of the atmosphere is limited by the scarcity of the observations they are based on, so improvement of these models requires validation under different atmospheric conditions. Measurements of the absorption of the solar extreme ultraviolet (EUV) radiation in the upper atmosphere below 500 km by instruments operating on lowEarth orbits (LEO) satellites provide efficient means for such validation as well as for continuous monitoring of the upper atmosphere and for studying its response to the solar and geomagnetic activity.
Method: This paper presents results of measurements of the solar EUV radiation in the 17 nm wavelength band made with the SPIRIT and TESIS telescopes on board the CORONAS satellites and the SWAP telescope on board the PROBA2 satellite in the occulted parts of the satellite orbits. The transmittance profiles of the atmosphere at altitudes between 150 and 500 km were derived from different phases of solar activity during solar cycles 23 and 24 in the quiet state of the magnetosphere and during the development of a geomagnetic storm. We developed a mathematical procedure based on the Tikhonov regularization method for solution of illposed problems in order to retrieve extinction coefficients from the transmittance profiles. The transmittance profiles derived from the data and the retrieved extinction coefficients are compared with simulations carried out with the NRLMSISE00 atmosphere model maintained by Naval Research Laboratory (USA) and the DTM2013 model developed at CNES in the framework of the FP7 project ATMOP.
Results: Under quiet and slightly disturbed magnetospheric conditions during high and low solar activity the extinction coefficients calculated by both models agreed with the measurements within the data errors. The NRLMSISE00 model was not able to predict the enhancement of extinction above 300 km observed after 14 h from the beginning of a geomagnetic storm whereas the DTM2013 model described this variation with good accuracy.
Key words: Thermosphere / XUV/EUV flux / Storm / Modelling / Solar activity
© V. Slemzin et al., Published by EDP Sciences 2016
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1. Introduction
The influence of the Earth’s atmosphere is very important for spacecraft and scientific instruments operating in nearEarth space. The lifetime of a satellite at altitudes of several hundred km strongly depends on the density of the atmosphere and its variation with solar activity, so the prediction accuracy of atmospheric conditions has a significant economic aspect. The solar shortwavelength (Xray and EUV) radiation is appreciably attenuated in the upper atmosphere at altitudes below 500 km and fully absorbed below 100 km, which should be taken into account in planning solar observations from satellites in the lowEarth orbits (LEO). On the other hand, due to the high sensitivity of solar observations to atmospheric absorption, the data obtained in the presence of absorption can be used for validation and improvement of atmosphere models.
Spectrometric investigations of the transparency of the Earth atmosphere in short wavelengths as a function of altitude began with rocket launches and first satellite flights during the 1950s to 1970s with the aim of studying absorption of the atmosphere under different observational conditions (see the review by Smith & Hunten 1990 and references therein). For example, Hinteregger (1962) describes such measurements at low zenith angles in the spectral band 28.3–121.6 nm for altitudes of 100–200 km. Later, measurements were carried out from satellites during occultations at the sunset or sunrise parts of the orbit. Below the 100 km level, the measurements of radiation flux variations during solar or stellar occultations are used for monitoring of ozone, nitric, and halogen molecules as well as other minor atmosphere constituents in many spectral bands from UV to infrared (experiments Nimbus7/TOMS, Bowman & Krueger 1985; SAGE I–III, Chu et al. 2003; Damadeo et al. 2013; UARS/HALOE, Russell et al. 1993; ENVISAT/GOMOS, Bertaux et al. 2010). The solar occultation technique was widely used for remote sensing of the atmospheres of planets and satellites such as Mars, Saturn, and Titan (Smith et al. 1983; Smith & Hunten 1990; Krasnopolsky et al. 1991; Korablev 2002; Maltagliati et al. 2013). The foundation of the occultation method including applications to radio waves and the GPSbased localization were treated in detail in the reviews of Elliott (1979), Smith & Hunten (1990) and in the book “Occultations for Probing Atmosphere and Climate” (Kirchengast et al. 2004).
Systematic measurements of the absorption of solar EUV radiation in the Earth’s atmosphere are not numerous. Since 2002, the Thermosphere Ionosphere Mesosphere Energetics and Dynamics (TIMED) satellite is monitoring the EUV solar radiation. The Solar EUV Experiment (SEE) on board TIMED (Woods et al. 2005) provides occultation measurements of solar irradiance through the upper atmosphere at different EUV wavelengths, but due to telemetry constraints measurements have been performed at only one single altitude on each orbit. Thus, the whole profile in the range 100–400 km is summarized from measurements on many orbits.
Due to the integration of the measured absorption along the lineofsight (LOS), occultation measurements have the evident drawback, as a method for validation of atmosphere models, of not permitting to obtain local densities directly. At the same time, this method provides selfcalibration of the data, high height resolution, and short time of scanning through a wide range of altitudes which allows measuring the atmosphere absorption on a wide spatial scale.
The ordinary method of solar occultation uses measurements of the total flux integrated over the full solar disk as a function of the socalled tangential altitude h (see Fig. 1), i.e. the nearest distance between the LOS and the Earth surface (subLOS point). Therefore, the transmittance profile is averaged over the altitude interval corresponding to the angular size of the Sun (~0.53°), which for a typical satellite orbit with an altitude of about 600 km varies from 1–2 km at h ~ 400 km to ~20 km at h ~ 200 km. The use of highresolution solar images, instead of integrated fluxes, permits to diminish this uncertainty to the size of one pixel, which corresponds to only several arc seconds, thus increasing the spatial resolution to scales below 1 km. Additionally, bright active regions, where the flux is the most variable, may easily be isolated by intensity and excluded from the integration (Slemzin et al. 2003).
Figure 1. Scheme of the occultation measurements. R_{0} – radius of the Earth, H_{orb} – altitude of the satellite orbit. Other notations are given in the text. 
Since the launch of the first spaceborne EUV solar telescope Phobos/TEREK in 1988 (Zhitnik et al. 1998) and the next CORONASI/TEREKC telescope in 1994 (Sobelman et al. 1996), it was realized that atmospheric absorption effects are very important for regular solar imaging in the EUV spectral range. One can completely avoid the absorption effects if the spacecraft is in the Lagrangian L1 point like the SOHO observatory, or sufficiently diminish them in time if the satellite operates in a sunsynchronous orbit. Due to precession, all lowEarth orbits (LEO) have seasonal periods, when the lineofsight crosses the Earth atmosphere. In most cases, observational programs of solar instruments do not include such periods in order to avoid distortion of the results by attenuated data. However, the scientific programs of the CORONASF/SPIRIT (Slemzin et al. 2005), CORONASPhoton/TESIS (Kuzin et al. 2009), and PROBA2/SWAP (Seaton et al. 2013) telescopes foresaw special periods of observations of the Sun in occultations to investigate the response of the upper atmosphere to longterm or shortterm variations of solar and geomagnetic activity. The first results of determination of atmospheric transmittance carried out by SPIRIT in the 0.842, 17.5, and 30.4 nm spectral bands as a function of tangential altitude in solar maximum (October–November 2001) were published by Zhitnik et al. (2003).
So far, several models intended for the calculation of characteristics of the upper atmosphere such as densities of the main constituents and temperature have been developed. These models allow calculating densities of the main atmosphere constituents at different altitudes, for given geographic coordinates, date, time, and solar and geomagnetic conditions. The first models (e.g. Jacchia 1971) were based on satellite drag data only. The series of the MassSpectrometerIncoherentScatter (MSIS) models progressively improved accuracy through the use of the mass spectrometer and incoherent scatter radar data. These models are named MSISE83, MSISE86, MSISE90 (Hedin 1983, 1987, 1991) and the latest NRLMSISE00 (Picone et al. 2002), which allows calculating densities of the main constituents and kinetic temperatures at altitudes from 0 to 1000 km. The NRLMSISE00 model can be run online at http://ccmc.gsfc.nasa.gov/modelweb/models/nrlmsise00.php.
Comparisons of the NRLMSISE00 model calculations with the drag measurements on the satellites ARGOS in the altitude range 150–200 km (2000–2001, Aksnes et al. 2007) and CHAMP at about 400 km (2002, Zhou et al. 2009) have shown that under quiet geomagnetic conditions the temperature and density height distributions given by the model were in satisfactory agreement with the measurements. However, it was found, that the model was unable to describe well the rapid enhancements in the thermospheric mass density and temperature during geomagnetic storms in October–November 2003. Density enhancements during storms are very important for the LEO spacecraft such as the International Space Station, because they can seriously alter the orbit and reduce altitude and useful lifetime of the mission.
A series of semiempirical Drag Temperature Models (DTMs) were developed since 1978 (Barlier et al. 1978). The latest advanced model DTM2013 has been developed under the ATMOP FP7 project with the goal of meeting the requirements for operational orbit computations (Dudok de Wit & Bruinsma 2011; Bruinsma et al. 2012; Bruinsma 2015). This model was constructed using drag observations of 10 satellites, gathered during the period of 1969–2011. The data include drag measurements obtained with accelerometers on board the CHAMP (Bruinsma et al. 2004), GRACE (Bruinsma & Forbes 2008), and GOCE (Bruinsma et al. 2014) satellites and use the F30 (30 cm radio flux) solar activity proxy which proved to be more representative than the F10.7 index (Dudok de Wit et al. 2014). A detailed description of the DTM2013 model and its comparison with the previous models DTM2009 (Bruinsma et al. 2012) and JB2008 are given in Bruinsma (2015). The online version of the model can be run at http://www.atmop.eu/index.php/models). In our work, we employed the model DTM2013 for comparison with absorption measurements.
Measurements of atmospheric attenuation of solar EUV radiation constitute an independent way of validating these atmosphere models at different levels of solar and geomagnetic activity. The local extinction coefficients can be retrieved mathematically from the measured absorption and compared with the model predictions. Here we describe the results of occultation measurements of the atmospheric transmittance in the 17 nm range at different altitudes in the quiet state of the atmosphere and during the development of a geomagnetic storm. The data have been obtained with the SPIRIT and TESIS telescopes on board the CORONAS satellites as well as with the SWAP telescope on board PROBA2. We compare the transmittance profiles and the retrieved extinction coefficients as a function of altitude with calculations by the NRLMSISE00 and DTM2013 atmosphere models.
2. Data and analysis
Figure 1 shows the scheme of solar observations from satellites on LEO orbits in occultation mode. The lineofsight (LOS) crosses the atmosphere at altitude h (tangential height) and forms an angle z with the direction to the Earth’s center. The rays from the satellite to different points on the solar disk deviate from the central ray by angles θ within half of the angular size of the Sun (~0.53°). This corresponds to variations dh(θ) of the tangential height h.
Due to the finite size of the Sun, the highest and the lowest rays differ in height by Δh ≈ α(R_{0} + H_{orb})cos z, where α is the angular size of the Sun, H_{orb} is the altitude of the orbit, z is the angle between LOS and direction to the Earth’s center, and R_{0} is the Earth’s radius. For instance, if H_{orb} = 600 km and h = 300 km, ∆h ≈ 18 km.
Transmittance T of the atmosphere at a specific wavelength λ along the ray corresponding to a LOS height h and angle θ is given by the integral over the LOS(1)where k_{λ}(x) is the wavelengthdependent extinction coefficient. Generally, extinction (or attenuation) of the light is caused by three processes – absorption, refraction, and scattering of light in the medium. Due to the very small density of the atmosphere at altitudes above 100 km (more than seven orders less than that on the ground), refraction and scattering produce negligible effects, so photoabsorption caused by photoionization of the atmospheric constituents through solar EUV radiation dominates. Thus, the extinction coefficient can be expressed in the form:(2)where d_{i}(h) and σ_{i}(λ) are the density of the ith constituent (provided by the atmosphere model) and its photoabsorption crosssection at wavelength λ.
In general, using solar imaging telescopes one can obtain transmittance profiles in two ways: by integrating the solar flux over the whole disk (the integrated mode), or by dividing the disk into many horizontal strips corresponding to discrete heights and by constructing the total transmittance height profile from partial profiles for these strips (the imaging mode).
In the integrated mode, which was used in this work, the signaltonoise ratio is the highest, but the resulting profile is averaged over the height divergence interval. However, it can be shown, that for a quasiuniformly radiant solar disk (which is the case in the EUV 17 nm band) the resulting profile coincides with that for the central ray within ∆h ~ 0.5 km. In the imaging mode, the height resolution is defined by the size of the resolved image cell which can be set as small as possible only limited by the signal level and by the growing influence of local temporal variations of the solar emission. The tests have shown that by using the imaging mode during quiet solar condition one can easily achieve a geometrical height resolution better than 1 km (Slemzin et al. 2003).
It is worth noting that even in the integrated mode the height divergence due to the solar size is less than the atmosphere scale height H_{s}. Calculations with DTM2013 model for a 600 km orbit in the tangential height range h = 200–500 km give values of H_{s} that vary from 42 to 104 km at solar maximum (the SPIRIT 1st data set) and from 24 to 120 km at minimum (the TESIS data set). Taking into account smooth dependence of density on height in the given range, the finite solar size cannot produce significant errors in the transmittance profiles.
We constructed transmittance profiles using solar images obtained in the sunset or sunrise sections of the orbit. To calculate the positions of the satellites and the subLOS points (at the Earth surface) at the moment of image acquisition we used twoline element (TLE) ballistic data provided by the Joint Space Operations Center of the US Strategic Command (https://www.spacetrack.org), processed with the TRAKSTAR code developed by T.S. Kelso (http://celestrak.com/software/tskelsosw.asp). As a rule, the TLE catalog contains one or two ballistic data sets per day for a given satellite from which we selected the nearest data set in time to the measured image. For quiet Sun periods (the TESIS and SWAP data sets), it resulted in an accuracy in the determination of the tangential height of ~25–50 m (this value was derived from the difference by two subsequent TLEs for the same instant of time). Even at solar maximum (the SPIRIT case), the error in h did not exceed 100 m.
The initial images were processed by standard procedures for each instrument which included corrections for the background, the flatfield distortion, and the nonlinearity. Then the intensity was integrated over the disk area excluding the corona above the limb where the EUV flux can significantly vary in time due to solar activity. To obtain transmittance profiles, the integrated light curves were determined as a function of tangential height and then normalized by the mean value of intensity above 500 km. At this altitude the absorption is comparable to or less than the relative variation of the solar flux due to solar activity and can therefore be neglected. In the studied cases, the dispersion of the solar flux at h > 500 km in the considered spectral range did not exceed 0.3% at maximum activity (SPIRIT) and 0.1% at minimum activity (TESIS and SWAP).
We validated the models by confronting them with the atmosphere absorption observed at different levels of solar activity (near maximum, at minimum, and in the growing phase of the solar cycle) under quiet and disturbed geomagnetic conditions. In the quiet or slightly disturbed geomagnetic state, we analyzed the following data sets: the SPIRIT images in the 17.5 nm band obtained on March 5, 2002, between 02:45:35 UT and 03:11:48 UT (data set 1 with cadence 32 s near maximum of solar cycle 23, geomagnetic index ap = 32 nT); the TESIS images in the 17.1 nm band (near solar minimum between cycles 23 and 24, July 17, 2009, 13:45:30 UT–July 19, 2009, 15:20:35 UT, cadence 5 min, ap < 3 nT); the SWAP images in the 17.4 nm band (the growing phase of cycle 24, December 22, 2010, 08:27:19 UT–08:34:39 UT, cadence 20 s, ap = 0). The response of the models to disturbed geomagnetic conditions was studied during the development of a geomagnetic storm using the SPIRIT data set 2 obtained before the storm on 23 May 2002 01:54:20–02:20:33 UT (ap = 12 nT, Dst = −1 nT), and data set 3 during the storm, 15:59:59–16:48:56 UT (ap = 236 nT, Dst = −75 nT). The SPIRIT and SWAP data were collected in the sunrise parts of single orbits, the TESIS data were combined from the sunrise and sunset parts of 30 orbits. Spectral functions for all instruments have halfwidths of 0.6–1 nm which is wider than the differences in the peak wavelengths, so they operated in similar wavelength bands.
Figure 2A presents the temporal positions of these data superimposed on the evolution curve of solar activity as testified by the sunspot cycle taken from the WDCSILSO at the Royal Observatory of Belgium in Brussels (http://sidc.oma.be/silso). Figure 2B shows the transmittance profiles for SPIRIT, TESIS, and SWAP corresponding to solar maximum of the cycle 23, and the minimum and growing phase of the solar cycle 24. Points in Figure 2B correspond to single measurements, solid curves to the fitted profiles. With growing solar activity, the increase of the transmittance from ~0 to ~1 occurs at higher altitudes and extends over a wider altitude range. Difference in transmittance is maximal at the heights 220–350 km, where the atmosphere absorbs the major part of solar EUV flux. The map in Figure 3 shows the geographic locations of the subsatellite and subLOS points for the input data. During the measurements of transmittance profiles, the subLOS points shift in latitude and longitude by tens of degrees, so the obtained profiles characterize the distribution of atmospheric densities over a wide spatial range.
Figure 2. (A) Sunspot numbers in solar cycles 23 and 24 and observation times of the SPIRIT, TESIS, and SWAP data (explanations in the text); (B) transmittance profiles for SPIRIT (solar maximum, data set 1, red line), TESIS (solar minimum, black line), and SWAP (rising phase of cycle 24, blue line). 
Figure 3. Map of occultation locations: red symbols – locations of subsatellite points for SPIRIT (data set 1), TESIS and SWAP, blue symbols – locations of subLOS points. 
3. Comparison of experimental transmittance profiles with model calculations
Figure 4 displays the height profiles of extinction coefficients at 17 nm wavelength of the main atmospheric constituents in the height range 200–500 km – atomic oxygen, molecular nitrogen, molecular oxygen, and helium – calculated with the atmosphere models NRLMSISE00 and DTM2013 at solar maximum, representative for the SPIRIT data set 1, and at solar minimum, representative for the TESIS data set. The photoabsorption crosssections for the dominant atomic components (O, N, He) were taken from Verner et al. (1996) and Fennelly & Torr (1992), for O_{2} and N_{2} from Verner et al. (1996) and Cole & Dexter (1978).
Figure 4. Height profiles of extinction coefficients of the main atmospheric constituents calculated with the NRLMSISE00 and DTM2013 models for the EUV wavelength 17 nm in the quiet geomagnetic state. (A) At solar maximum (the SPIRIT data set 1), (B) at solar minimum (the TESIS data set). 
A comparison has shown that in the spectral interval 5–35 nm the difference between the data from different sources does not exceed 2% for O and N and 10% for O_{2} and N_{2} which is close to the theoretical data accuracy (5–8%). Extinction coefficients for different components calculated with densities provided by both models depend differently on height and level of solar activity. The total extinction calculated with the NRLMSISE00 model at solar maximum is higher than that calculated with the DTM2013 model (the ratio averaged over the most important height range 200–500 km is 1.49), at solar minimum the former model gives the smaller value (the ratio is 0.74).
The transmittance profiles corresponding to the abovementioned SPIRIT, TESIS, and SWAP data were calculated according to Eqs. (1) and (2). The transmittance along LOS as a function of height was determined by numerical integration of the extinction under the assumption of spherical symmetry of the atmosphere. In Figure 5, the modeled profiles for the SPIRIT, TESIS, and SWAP data sets are compared with those derived from measurements.
Figure 5. Comparison of the transmittance height profiles, presented in Figure 2B, with those calculated by the NRLMSISE00 and DTM2013 models. The black lines display the averaged measurement points (for TESIS individual points are also shown), the red/blue lines show simulations by the NRLMSISE00 and DTM2013 models. In the SPIRIT and SWAP cases, the error bars correspond to the transmittance determination errors, in the TESIS case to the dispersion of data averaged over multiple profiles. 
The difference between the measured and modeled profiles can be expressed in mean displacement of the transmittance height profile (in km) averaged over the transmittance range 0.1–0.9 (Table 1). The profiles calculated with the NRLMSISE00 model are closer to the data fitting functions for periods of high solar activity (SPIRIT set 1) and moderate activity (SWAP), the DTM2013 model profile is closer for solar minimum (TESIS), although the mean displacements in the SPIRIT and TESIS cases are below the data errors, which included the intensity measuring error and variation of the unabsorbed solar flux during measurements.
Displacement of the modeled transmittance height profiles from those derived from the data presented in Figure 5: <∆h_{MSISE}>, <∆h_{DTM}> – mean displacements in height for the NRLMSISE00 and DTM2013 models correspondingly; <∆h_{dat}> – mean error of the derived height profile (the values were averaged over the transmittance interval 0.1–0.9).
4. Comparison of the retrieved extinction coefficients with the modeled values
As it was shown in Section 2, a comparison of the modeled transmittance height profiles with those derived from measurements gives an indication of how the models describe the largescale distribution of atmospheric densities, but does not permit a detailed estimation of the accuracy of their local modeled values. For this purpose, in the cases of SPIRIT and SWAP we calculated the local extinction coefficients as a function of height at subLOS points with both models and compared them with the values retrieved from the measured profiles.
We have developed an approach for solving the inverse problem of retrieval of the extinction coefficients from measured profiles based on Tikhonov’s theory of illposed problems and regularization method (Tikhonov & Arsenin 1977), which allows overcoming the inevitable limitations inherent in the foregoing modelbased consideration; in particular, this approach makes it possible to retrieve variations of the absorption coefficient with fine altitude resolution. The basic principles of this method and its application to the analysis of the attenuation profiles are presented in Annex A. While the attenuation profiles characterize the density distribution in the atmosphere in the wide scale, the extinction coefficients are directly defined by the local density near the subLOS point at a given height, so they are very useful for validation of the atmosphere models. Figures 6 and 7 illustrate a comparison of the extinction coefficients (at the subLOS points) retrieved by the method described in Annex A with those calculated using the densities provided by the NRLMSISE00 and DTM2013 models.
Figure 6. Comparison of the modeled and retrieved extinction coefficients as a function of height for the SPIRIT data set 1 (March 5, 2002, left) and the SWAP data set (December 22, 2010, right). Black line – the data retrieved from measurements, red/blue lines – calculations by the NRLMSISE00 and DTM2013 models. 
Figure 7. Variation of the retrieved and modeled extinction coefficients during the development of the geomagnetic storm of 23 May 2002: (A) before the storm (the SPIRIT data set 2, Dst = −1 nT), (B) 15 h later (the SPIRIT data set 3, Dst = −75 nT). The line colors are the same as in Figure 6. 
4.1. Quiet atmosphere at different levels of solar activity
Figure 6 shows a comparison of the extinction coefficients as a function of height retrieved from the SPIRIT data (set 1, March 5, 2002) and the SWAP data (December 22, 2010) with those modeled with NRLMSISE00 and DTM2013. The SPIRIT data correspond to a slightly disturbed geomagnetic state (the geomagnetic index ap = 32 nT), the SWAP data to the quiet state (ap = 0 nT).
As seen in Figure 6, at altitudes from 220 to 460 km, in the SPIRIT case the NRLMSISE00 model underestimates the extinction coefficients of ~30%, whereas the DTM2013 model agrees with the data within 5%. In the SWAP case, the models show mean deviations from the data that reach 48% (NRLMSISE00) and 39% (DTM2013) at heights above 300 km. It can probably be explained by the fact that the PROBA 2 satellite operates in a dawndusk orbit, so the subLOS points for SWAP are located near the terminator where density gradients and the uncertainty of the simulation are the largest. The accuracy of the retrieved extinction coefficients in all cases was about 5% due to the selfconsistent mathematical procedure (see Annex A2), which is better than the accuracy of the individual data values.
4.2. Disturbed atmosphere (geomagnetic storm)
Figure 7 shows the variation of the retrieved extinction coefficients during the development of the moderate geomagnetic storm on May 23, 2002, 01–15 h UT (the SPIRIT data sets 2 and 3). During this period, the ap index rose from 12 nT to 236 nT, and the Dst index fell from −1 nT to −75 nT. The storm began at ~9 h UT (http://omniweb.gsfc.nasa.gov). The largest variation of the extinction occurred at the heights 300–450 km (increase at 400 km by a factor of 1.7) which was probably produced by density enhancements in the ionosphere associated with magnetic disturbances during the storm (Hedin et al. 1977).
The NRLMSISE00 model describes well (within the data accuracy) the variation of the extinction coefficient in the quiet state before the storm, but fails to describe its growth during the storm (difference about 70%). The DTM2013 model describes well the behavior of the extinction coefficients above 300 km height in the quiet state (accuracy ~12%) and under storm conditions (~1%). Below 300 km the extinction coefficient given by this model exceeds the one retrieved from measurements, at 200 km by a factor of 2. At the moment, the reasons for such a discrepancy are not clear.
5. Summary and conclusions

The study of the atmospheric attenuation of solar EUV radiation using telescopes and photometers on board LEO satellites during occultations provides valuable information about the properties and condition of the Earth’s upper atmosphere and can be used for validating atmosphere models. The EUV wavelength band is well suited for the investigation of the atmosphere at heights 200–500 km where the influence of solar and geomagnetic activity is the largest.

The transmittance height profiles during solar occultations obtained from the total integrated flux over the whole disk have a limited height resolution of 10–20 km due to averaging over the finite angular size of the Sun. However, under the assumption of a quasiuniform distribution of brightness over the solar disk (which is valid in the EUV 17 nm band), the resulting profile coincides with that of the ray from the solar center within ±0.5 km. Earlier (Slemzin et al. 2003) it was shown that the height resolution of the profiles can be increased to less than 1 km if using solar images segmented into narrow strips.

For the analysis of the atmospheric absorption we used solar observation data obtained during orbital occultations by the telescopes CORONASF/SPIRIT (17.5 nm), CORONASPhoton/TESIS (17.1 nm), and PROBA2/SWAP (17.4 nm). The selected data were obtained at different levels of solar activity: at the maximum of cycle 23 (SPIRIT, March–May 2002), in the deep minimum between cycles 23 and 24 (TESIS, July 2009), and at moderate solar activity during the growing phase of the cycle 24 (SWAP, December 2010).

We compared the obtained transmittance profiles at the heights 200–500 km with simulations based on the atmosphere models NRLMSISE00 and DTM2013. The profiles calculated with the NRLMSISE00 model are closer to the data fitting functions for the periods of high solar activity (SPIRIT set 1, the mean displacement about 1 km) and during a moderately active Sun (SWAP data set, displacement 4.5 km), the DTM2013 model for solar minimum (TESIS data set, displacement 1.7 km).

Extinction coefficients were retrieved with an accuracy of about 5% from the measured transmittance profiles by solving the inverse problem using Tikhonov’s theory of illposed problems and regularization method based on the generalized discrepancy principle. It was found that at altitudes between 220 and 460 km, under slightly disturbed magnetospheric conditions (in the SPIRIT case) the NRLMSISE00 model underestimates the extinction coefficients by ~30%, whereas the DTM2013 model agrees with them within 5%. In the SWAP case (quiet magnetosphere), the mean deviations from the data at heights above 300 km are higher (48% for NRLMSISE00 and 39% for DTM2013), which can probably be explained by the uncertainty of the simulations in regions of high density gradients.

The response of the upper atmosphere to geomagnetic disturbances was demonstrated by the comparison of the retrieved extinction coefficients for the SPIRIT orbits 4508 and 4524 (23 May 2002, with starting times 01:54 and 15:59 UT) corresponding to transition from the quiet magnetospheric conditions (ap = 12 nT, Dst = −1 nT) to a moderate storm (ap = 236 nT, Dst = −75 nT). Under geomagnetic storm conditions we found that the retrieved extinction coefficients at 17 nm at heights between 280 and 450 km were larger by up to 30–60% in comparison with those in the quiet geomagnetic state. Below 280 km, the extinction remained the same during the storm. The DTM2013 model describes this variation of extinction during the storm very well, whereas the NRLMSISE00 model failed to reproduce it.
In conclusion, our investigation has shown that the method of solar occultations in the EUV range is a very sensitive tool to study the properties of the upper atmosphere at heights 200–500 km and its response to solar and geomagnetic activity.
Acknowledgments
The research leading to these results has received funding from the European Commission’s Seventh Framework Programme (FP7/20072013) under the Grant Agreement eHeroes (project No 284461, http://www.eheroes.eu). SWAP is a project of the Centre Spatial de Liège and the Royal Observatory of Belgium funded by the Belgian Federal Science Policy Office (BELSPO). The SPIRIT and TESIS experiments have been developed and carried out by the Laboratory of Xray Astronomy of the Sun, P.N. Lebedev Physical Institute of the Russian Academy of Sciences, under the CORONAS International solar space program. The editor thanks two anonymous referees for their assistance in evaluating this paper.
References
 Aksnes, A., R. Eastes, S. Budzien, and K. Dymond. Dependence of neutral temperatures in the lower thermosphere on geomagnetic activity. J. Geophys. Res., 112, A06302, 2007, DOI: 10.1029/2006JA012214. [CrossRef] [Google Scholar]
 Barlier, F., C. Berger, J.L. Falin, G. Kockarts, and G. Thuillier. A thermospheric model based on satellite drag data. Ann. Geophys., 34, 9–24, 1978. [Google Scholar]
 Bertaux, J.L., E. Kyrölä, D. Fussen, A. Hauchecorne, F. Dalaudier, et al. Global ozone monitoring by occultation of stars: an overview of GOMOS measurements on ENVISAT. Atmos. Chem. Phys., 10, 12091–12148, 2010, DOI: 10.5194/acp10120912010. [NASA ADS] [CrossRef] [Google Scholar]
 Bowman, K.P., and A.J. Krueger. A global climatology of total ozone from the Nimbus7 Total Ozone Mapping Spectrometer. J. Geophys. Res., 90, 7967–7976, 1985, DOI: 10.1029/JD090iD05p07967. [CrossRef] [Google Scholar]
 Bruinsma, S. The DTM2013 thermosphere model, J. Space Weather Space Clim., 5, A1, 2015, DOI: 10.1051/swsc/2015001. [CrossRef] [EDP Sciences] [Google Scholar]
 Bruinsma, S., D. Tamagnan, and R. Biancale. Atmospheric densities derived from CHAMP/STAR accelerometer observations. Planet. Space Sci., 52, 297–312, 2004, DOI: 10.1016/j.pss.2003.11.004. [CrossRef] [Google Scholar]
 Bruinsma, S.L., and J.M. Forbes. Medium to largescale density variability as observed by CHAMP. Space Weather, 6, S08002, 2008, DOI: 10.1029/2008SW000411. [CrossRef] [Google Scholar]
 Bruinsma, S.L., N. SánchezOrtiz, E. Olmedo, and N. Guijarro. Evaluation of the DTM2009 thermosphere model for benchmarking purposes. J. Space Weather Space Clim., 2, A04, 2012, DOI: 10.1051/swsc/2012005. [CrossRef] [EDP Sciences] [Google Scholar]
 Bruinsma, S.L., E. Doornbos, and B.R. Bowman. Validation of GOCE densities and thermosphere model evaluation. Adv. Space Res., 54, 576–585, 2014, DOI: 10.1016/j.asr.2014.04.008. [CrossRef] [Google Scholar]
 Chu, W.P., C.R. Trepte, and G. Taha. Initial comparison of SAGE III Data with GOMOS and SCIAMACHY, in: Proc. of Envisat Validation Workshop, Frascati, Italy, ESA SP531, 2003. [Google Scholar]
 Cole, B.E., and R.N. Dexter. Photoabsorption and photoionisation measurements on some atmospheric gases in the wavelength region 50–340 A. J. Phys. B: At. Mol. Phys., 11, 1011–1023, 1978, DOI: 10.1088/00223700/11/6/013. [CrossRef] [Google Scholar]
 Damadeo, R.P., J.M. Zawodny, L.W. Thomason, and N. Lyer. SAGE version 7.0 algorithm: application to SAGE II. Atmos. Meas. Tech., 6, 3539–3561, 2013, DOI: 10.5194/amt635392013. [CrossRef] [Google Scholar]
 Dudok de Wit, T., and S. Bruinsma. Determination of the most pertinent EUV proxy for use in thermosphere modelling. Geophys. Res. Lett., 38, L19102, 2011, DOI: 10.1029/2011GL049028. [CrossRef] [Google Scholar]
 Dudok de Wit, T., S. Bruinsma, and K. Shibasaki. Synoptic radio observations as proxies for upper atmosphere modelling. J. Space Weather Space Clim., 4, A06, 2014, DOI: 10.1051/swsc/2014003. [CrossRef] [EDP Sciences] [Google Scholar]
 Elliott, J.L. Stellar occultation studies of the solar system. Ann. Rev. Astron. Astrophys., 17, 445–475, 1979, DOI: 10.1146/annurev.aa.17.090179.002305. [CrossRef] [Google Scholar]
 Fennelly, J.A., and D.G. Torr. Photoionization and photoabsorption cross sections of O, N2 O2, and N for aeronomic calculations. At. Data Nucl. Data Tables, 51, 321–363, 1992, DOI: 10.1016/0092640X(92)900042. [NASA ADS] [CrossRef] [Google Scholar]
 Gaikovich, K.P. Inverse Problems in Physical Diagnostics. Nova Science Publishers Inc., New York, 2004. [Google Scholar]
 Gaikovich, K.P., A.S. Gurevich, and A.P. Naumov. On a reconstruction of meteorological parameters from intraatmospheric measurements of optical refraction of cosmic sources. Izvestiya, Atmospheric and Oceanic Physics, 19 (7), 507–512, 1983. [Google Scholar]
 Hedin, A.E. A revised thermospheric model based on mass spectrometer and incoherent scatter data  MSIS83. J. Geophys. Res., 88 (A12), 10170–10188, 1983, DOI: 10.1029/JA088iA12p10170. [CrossRef] [Google Scholar]
 Hedin, A.E. MSIS86 thermospheric model. J. Geophys. Res., 92, 4649, 1987, DOI: 10.1029/JA088iA12p10170. [NASA ADS] [CrossRef] [Google Scholar]
 Hedin, A.E. Extension of the MSIS thermosphere model into the middle and lower atmosphere. J. Geophys. Res., 96, 1159, 1991, DOI: 10.1029/90JA02125. [NASA ADS] [CrossRef] [Google Scholar]
 Hedin, A.E., H.G. Mayr, L.H. Brace, H.C. Brinton, D.T. Pelz, P. Bauer, G.R. Carignan, and A.D. Parks. Observations of neutral composition and related ionospheric variations during a magnetic storm in February 1974. J. Geophys. Res., 82, 3183–3189, 1977, DOI: 10.1029/JA082i022p03183. [CrossRef] [Google Scholar]
 Hinteregger, H.E. Absorption spectrometric analysis of the upper atmosphere in the EUV region. J. Atmos. Sci., 19, 351–368, 1962, DOI: 10.1175/15200469(1962)019<0351:ASAOTU>2.0.CO;2. [CrossRef] [Google Scholar]
 Jacchia, L.G. Revised Static Models of the Thermosphere and Exosphere with Empirical Temperature Profiles. SAO Special Report #332, 1971. [Google Scholar]
 Kirchengast G., U. Foelsche, and A. Steiner, Editors. Occultations for Probing Atmosphere and Climate, Springer Science & Business Media, SpringerVerlag Berlin Heidelberg, ISBN: 3540341218, 9783540341215, 2004. [CrossRef] [Google Scholar]
 Korablev, O.I. Solar occultation measurements of the Martian atmosphere on the Phobos spacecraft: water vapor profile, aerosol parameters, and other results. Sol. Syst. Res., 36, 12–34, 2002. [CrossRef] [Google Scholar]
 Krasnopolsky, V.A., O.B. Likin, F. Farnik, and B. Valnicek. Solar occultation observations of the Martian atmosphere in the ranges of 2–4 and 4–8 keV measured by PHOBOS 2. Icarus, 89, 147–151, 1991, DOI: 10.1016/00191035(91)90094A. [CrossRef] [Google Scholar]
 Kuzin, S.V., S.A. Bogachev, I.A. Zhitnik, et al. TESIS experiment on EUV imaging spectroscopy of the Sun. Adv. Space Res., 43, 1001, 2009, DOI: 10.1016/j.asr.2008.10.021. [NASA ADS] [CrossRef] [Google Scholar]
 Maltagliati, L., F. Montmessin, O. Korablev, A. Fedorova, F. Forget, A. Määttänen, F. Lefèvre, and J.L. Bertaux. Annual survey of water vapor vertical distribution and wateraerosol coupling in the martian atmosphere observed by SPICAM/MEx solar occultations. Icarus, 223, 942–962, 2013, DOI: 10.1016/j.icarus.2012.12.012. [CrossRef] [Google Scholar]
 Picone, J.M., A.E. Hedin, and D.P. Drob. NRLMSISE00 empirical model of the atmosphere: Statistical comparisons and scientific issues. J. Geophys. Res., 107 (A12), 1468, 2002, DOI: 10.1029/2002JA009430. [CrossRef] [Google Scholar]
 Russell, J.M. III, L.L. Gordley,, J.H. Park, S.R. Drayson, W.D. Hesket, et al. The Halogen occultation experiment. J. Geophys. Res., 98, 10777–10797, 1993, DOI: 10.1029/93JD00799. [CrossRef] [Google Scholar]
 Seaton, D.B., D. Berghmans, B. Nicula, J.P. Halain, A. De Groof, et al. The SWAP EUV imaging telescope part I: instrument overview and preflight testing. Solar Phys., 286, 43, 2013, DOI: 10.1007/s1120701201146. [NASA ADS] [CrossRef] [Google Scholar]
 Slemzin, V., O. Bugaenko, A. Ignatiev, V. Krutov, S. Kuzin, et al. Investigation of absorption of solar EUVradiation in the Earth’s atmosphere at altitudes of 100–500 km using solar images in the experiments TEREKC (CoronasI) and SPIRIT (CoronasF). In: A. Wilson, Editor, Solar Variability as an Input to the Earth’s Environment. International Solar Cycle Studies (ISCS) Symposium, Tanranská Lomnica, Slovak Republic, ESASP535, ESA Publications Division, Noordwijk, 389–392, ISBN 929092845X, 2003. [Google Scholar]
 Slemzin, V.A., S.V. Kuzin, I.A. Zhitnik, J.P. Delaboudinière, F. Auchère, et al. Observations of solar EUV radiation with the CORONASF/SPIRIT and SOHO/EIT instruments. Sol. Syst. Res., 39, 489–500, 2005. [NASA ADS] [CrossRef] [Google Scholar]
 Smith, G.R., and D.M. Hunten. Study of planetary atmospheres by absorptive occultations. Rev. Geophys., 28, 117–143, 1990, DOI: 10.1029/RG028i002p00117. [CrossRef] [Google Scholar]
 Smith, G.R., D.E. Shemansky, J.B. Holberg, A.L. Broadfoot, B.R. Sandel, and J.C. McConnell. Saturn’s upper atmosphere from the Voyager 2 EUV solar and stellar occultations. J. Geophys. Res., 88, 8667–8678, 1983, DOI: 10.1029/JA088iA11p08667. [NASA ADS] [CrossRef] [Google Scholar]
 Sobelman, I.I., I.A. Zhitnik, A.P. Ignat’ev, V.V. Korneev, V.Y. Klepikov, and V.V. Krutov. Xray spectroscopy of the Sun in the 0.84–30.4 nm band in the TEREKK and RESK experiments on the KORONAS satellite. Astron. Lett., 22, 539–554, 1996. [Google Scholar]
 Tikhonov, A.N., and V.Y. Arsenin. Solutions of illposed problems, Winston, New York, ISBN 0470991240, 1977. [Google Scholar]
 Verner, D.A., G.J. Ferland, K.T. Korista, and D.G. Yakovlev. Atomic data for astrophysics. II. New analytic FITS for photoionization cross sections of atoms and ions. Astrophys. J., 465, 487–498, 1996, DOI: 10.1086/177435. [NASA ADS] [CrossRef] [Google Scholar]
 Woods, T.N., F.G. Epavier, S.M. Bailey, P.C. Chamberlin, J. Lean, G.J. Rottman, S.C. Solomon, W.K. Tobiska, and D.L. Woodraska. Solar EUV Experiment (SEE): mission overview and first results. J. Geophys. Res., 110, A01312, 2005, DOI: 10.1029/2004JA010765. [NASA ADS] [CrossRef] [Google Scholar]
 Zhitnik, I.A., A.P. Ignatiev, V.V. Korneev, V.V. Krutov, S.V. Kuzin, et al. Instruments for imaging XUV spectroscopy of the Sun on board the CORONASI satellite. Proc. SPIE, 3406, 1–19, 1998. [CrossRef] [Google Scholar]
 Zhitnik, I.A., K.A. Boyarchuk, O.I. Bugaenko, G.S. IvanovKholodnyi, A.P. Ignat’ev, et al. Effects of absorption of solar XUV radiation by the Earth’s atmosphere at altitudes of 100–500 km in the Xray solar images obtained onboard the CORONASI (TEREK telescope) and CORONASF (SPIRIT Xray complex) satellites, Sol. Syst. Res., 37, 296–301, 2003. [CrossRef] [Google Scholar]
 Zhou, Y.L., S.Y. Ma, H. Lühr, C. Xiong, and C. Reigber. An empirical relation to correct stormtime thermospheric mass density modeled by NRLMSISE00 with CHAMP satellite air drag data. Adv. Space Res., 43, 819, 2009, DOI: 10.1016/j.asr.2008.06.016. [CrossRef] [Google Scholar]
Annex A: Retrieval of extinction coefficients from measured profiles
A1. Reduction of Sun measurements to a point source
Let us consider the occultation scheme (Fig. 1) and rewrite Eq. (1) with notations more convenient for the analysis. The received solar radiant intensity I^{m} (r_{0}, y_{0}) depends on the ray perigee position of the center of the solar disk r_{0} and it does not depend on the transversal position (to the xr plane) y_{0} of the Sun. This intensity is an integral over the disk surface D of the surface density of the radiant intensity J(r, y) (the radiant intensity of the received emission per unit disk surface). The value of J is also independent of the transversal position y of the corresponding element of the solar surface. Taking this into account, one has(A1)where r_{s} is the radius of the Sun projection around the tangential point (i.e., r_{s} = dh in Fig. 1).
This expression is valid for a uniformly radiating solar disk. Multiplying J(r) by the disk area, one obtains the intensity dependence I(r), such as it would be for the whole Sun located at the ray perigee position r. The relationship between the measured dependence I^{m}(r_{0}) and I(r) is expressed by the integral equation:(A2)
The kernel of (A2) is normalized to unity, hence in the case of slow variation of I(r) over the solar disk one has I^{m}(r_{0}) ≈ I(r). The solution of Eq. (A2) gives the dependence of radiance intensity I(r) of an equivalent point source located at a fixed ray perigee r. It is a Volterrakind equation with variable integration limits, and it is known as a weakly illposed problem. This equation can also be expressed as a Fredholm integral equation of the 1st kind with the kernel F (Eq. (A2)) that can easily be solved using Tikhonov’s method of generalized discrepancy. In the case of an inhomogeneous distribution of solar radiation over the disk, this distribution should be taken into account in the formation of the kernel F. The retrieved dependence I(r) satisfies the evident relationship following from (1):(A3)where τ is the optical depth along the ray of perigee r and I_{0} is the solar radiant intensity outside the atmosphere. From (A3) one can express the difference of the optical depths corresponding to the rays with perigee r and r_{m} as(A4)
Under the assumption that the Earth atmosphere is spherically symmetric, this difference is determined by the absorption altitude profile γ(r) as:(A5) where r_{H} is the radial distance of the receiver.
In our case, we can restrict the analysis to the situation where r < r_{H}. In this case, the corresponding inverse problems can be solved using the Abel inverse transform (in the refraction inverse problem) or the inversion of the Volterra integral equation of the 1st kind in its classical statement (Gaikovich 2004). If there would be some atmospheric contribution above the level r = r_{H}, Eq. (A5) transforms into a Volterra integral equation of the 1st kind in the illposed statement (Gaikovich 2004, p. 21).
A highaccuracy solution (like that demonstrated in this paper) can be achieved for the layers below the height of the receiver since the contribution of the perigee vicinity is strongly dominant because of the weak peculiarity in the kernel of (A5).
All our results were obtained by data at negative elevation angles (as illustrated in Fig. 1). It is nevertheless worth noting that by performing measurements also at the respective symmetric positive angles and looking at the difference between the two, one can guarantee that the solution does not depend on the atmosphere above the height of the receiver. In this case, instead of (A5), one has an ordinary Volterra equation of the 1st kind with a weak peculiarity that can easily be solved (Gaikovich 2004):(A6)
In general, we have to deal with the much more difficult problem of finding the solution of (A5) both below and above the altitude of the receiver. A previous application of the described method to the reconstruction of the refraction parameters of the atmosphere shows that at altitudes much lower than the receiver altitude one has a sufficiently accurate solution (Gaikovich et al. 1983; Gaikovich 2004); but at heights comparable to and higher than the receiver altitude, the kernel of (A5) changes – it has no weak peculiarity above this level. Hence, the solution quality in this range is inherent to illposed problems related to Fredholm integral equations of the 1st kind: the solution has a weak convergence when reducing the measurement errors, and the retrieval accuracy is sharply reduced with the increase of the altitude. An estimation of the retrieval accuracy can only be obtained, as in any illposed problem, from numerical simulation taking into consideration the convergence properties of the regularized solution.
A2. Method of absorption profiling and results of the numerical simulation
The considered problem consists in solving the integral Eqs. (A2) and (A5). Both equations can also be expressed as an integral Fredholm equation of the 1st kind(A7)where K_{h} is the approximation of the equation kernel K. The upper limit should be chosen so that the contribution of the range above this limit is smaller than the measurement errors. The lower limit can be chosen arbitrarily in the range r_{0min} – r_{s} ≤ r_{1} < r_{2}, where r_{0min} is the lower value of the perigee of the center of the solar disk in measurements. This makes it possible to solve the problem with various values of r1. This property can be useful when measurement errors are strongly altitudedependent.
To solve Eq. (A7), it is necessary to use regularization methods which involve additional a priori information about the function to be retrieved. The Tikhonov’s method of generalized discrepancy considered here (Tikhonov & Arsenin 1977) uses the information that the solution belongs to the set of squareintegrable functions with squareintegrable derivatives. In this method, an approximate solution y^{α} of (A3) is found using the parameter of measurement errors δ and the parameter of the kernel error h_{k} that satisfy the conditions:(A8)where z is the lefthand side of (A7) that corresponds to the exact solution y; z_{m} represents the measurements (affected by measurement errors). According to Tikhonov’s method, the approximate solution of (A7) is found from minimizing the generalized discrepancy:(A9)with the additional condition (principle of generalized discrepancy)(A10)and makes it possible to find the value of the regularization parameter at given values of δ and h.
An important advantage of this method is the convergence of the approximate solution to the exact one in the metric W_{1}^{2}, and, hence, the uniform convergence in the metric C (with the maximum modulus as a norm) at (δ, h_{k}) → 0 estimated in the integral metric L_{2}. So there is no need to use large statistical ensembles to obtain representative estimations of the retrieval accuracy as it is necessary in other methods with integral or mean square convergence. It is sufficient to run the numerical simulation for typical and, maybe, for extreme exact solutions to estimate the maximum error of the retrieved distribution. There is some uncertainty in the measurement error estimation related to the random origin of errors that should be taken into account. It is of advantage to search for a solution as the deviation from a reasonable model function. This deviation is typically a smoother and simpler function, and such functions are much better retrieved in solving illposed problems.
The accuracy of the proposed method has been studied by numerical simulations. These simulations involved the following steps: (a) for an initial simulated distribution γ(r) the exact lefthand side of (A3) I(r) is calculated; (b) random Gaussdistributed errors with rms that correspond to the real measurement accuracy are added to I(r) to obtain “measurement data”; (c) the inverse problem (A3) is solved to obtain Δτ(r_{0}); (d) these data are used to retrieve the distribution γ(r) from the solution of the inverse problem (A6); (e) the retrieved distribution γ(r) is compared to the initial one.
In numerical simulations, we used the initial distributions γ(h) obtained from deviations from the atmosphere model:(A11)
Where the atmospheric model(A12)is a function of h the altitude above the earth surface, as well as of several coefficients here obtained from the NRLMSISE00 model: k_{1} = −5.94, k_{2} = −0.00797, k_{3} = 2.36 × 10^{−7} for wavelength 17.5 nm.
The measure of retrieval error in solving (A6) (A13)was estimated from measurements as δ = 0.05.
From results of the abovedescribed numerical simulation with the level of random errors δ = 0.05 (in the integral metric L_{2}), true estimations of the accuracy of the retrieved profiles γ(h) were obtained in the maximum modulus metric. They demonstrate that the accuracy of retrieval depends on the depth scale ∆h of the profile perturbation: it changes from 10% for largescale variations at ∆h ~ 50 km or more up to 30% for smallscale variations at ∆h ~ 10 km.
Cite this article as: Slemzin V, Ulyanov A, Gaikovich K, Kuzin S, Pertsov A, et al. Validation of Earth atmosphere models using solar EUV observations from the CORONAS and PROBA2 satellites in occultation mode. J. Space Weather Space Clim., 6, A7, 2016, DOI: 10.1051/swsc/2015045.
All Tables
Displacement of the modeled transmittance height profiles from those derived from the data presented in Figure 5: <∆h_{MSISE}>, <∆h_{DTM}> – mean displacements in height for the NRLMSISE00 and DTM2013 models correspondingly; <∆h_{dat}> – mean error of the derived height profile (the values were averaged over the transmittance interval 0.1–0.9).
All Figures
Figure 1. Scheme of the occultation measurements. R_{0} – radius of the Earth, H_{orb} – altitude of the satellite orbit. Other notations are given in the text. 

In the text 
Figure 2. (A) Sunspot numbers in solar cycles 23 and 24 and observation times of the SPIRIT, TESIS, and SWAP data (explanations in the text); (B) transmittance profiles for SPIRIT (solar maximum, data set 1, red line), TESIS (solar minimum, black line), and SWAP (rising phase of cycle 24, blue line). 

In the text 
Figure 3. Map of occultation locations: red symbols – locations of subsatellite points for SPIRIT (data set 1), TESIS and SWAP, blue symbols – locations of subLOS points. 

In the text 
Figure 4. Height profiles of extinction coefficients of the main atmospheric constituents calculated with the NRLMSISE00 and DTM2013 models for the EUV wavelength 17 nm in the quiet geomagnetic state. (A) At solar maximum (the SPIRIT data set 1), (B) at solar minimum (the TESIS data set). 

In the text 
Figure 5. Comparison of the transmittance height profiles, presented in Figure 2B, with those calculated by the NRLMSISE00 and DTM2013 models. The black lines display the averaged measurement points (for TESIS individual points are also shown), the red/blue lines show simulations by the NRLMSISE00 and DTM2013 models. In the SPIRIT and SWAP cases, the error bars correspond to the transmittance determination errors, in the TESIS case to the dispersion of data averaged over multiple profiles. 

In the text 
Figure 6. Comparison of the modeled and retrieved extinction coefficients as a function of height for the SPIRIT data set 1 (March 5, 2002, left) and the SWAP data set (December 22, 2010, right). Black line – the data retrieved from measurements, red/blue lines – calculations by the NRLMSISE00 and DTM2013 models. 

In the text 
Figure 7. Variation of the retrieved and modeled extinction coefficients during the development of the geomagnetic storm of 23 May 2002: (A) before the storm (the SPIRIT data set 2, Dst = −1 nT), (B) 15 h later (the SPIRIT data set 3, Dst = −75 nT). The line colors are the same as in Figure 6. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.