A standalone prediction model for atomic oxygen and coronal mass ejections

This paper presents a standalone predictive model for Atomic Oxygen (AO), Coronal Mass Ejections (CMEs) and other space-environment parameters. The prediction is based on the numerical method of Holt–Winter’s triple smooth exponential forecasting of atmospheric constituents. Solar cycle 25 is likely to show about the same activity as cycle 23. The corresponding AO-flux–solar-activity correlation coefficients for altitudes 100, 200, and 300 km are: 0.62, 0.53, and 0.48, respectively, while the correlation coefficients for higher altitudes are lower than 0.48, an advantage that makes them more favorable for LEOs due to the harmful corrosive effects.


Introduction
Satellite protection depends greatly on accurate predictions of the surrounding space-environment components and hazards. The space environment of the Low-Earth Orbits (LEOs) is highly affected by solar activity. Solar activity is the key factor affecting the space environment. Assessment of the space environment depends strictly on three subprocesses or consecutive steps. The first step is the study of all components of the environment of the space mission by ad-dressing all questions about that mission (orbit, start time, duration, solar activity, type of mission, etc.). Secondly, detailed investigations of all the mission components should be done. The final step it to use all the analysis done to protect the mission in order to fulfill its target. Satellite protection requires full prediction of all the constraints and worst-case scenarios related to the mission. The prediction may be developed using deep learning, machine learning or statistical analyses.
During solar-activity cycles particle, radiation, and magnetic fluxes in the heliosphere change, which cause different space-weather effects at Earth. Prediction of the solaractivity cycle is therefore the main step in protecting space missions and satellite technology (Bhowmik and Nandy 2018). Space weather influences the electromagnetic environment around Earth and human life. Many space-weather events that are caused by solar eruptions, are potential risks to the social infrastructure such as aviation, communications, artificial satellites, electric power, and positioning systems (Kusano et al. 2021). Solar cycles 23 and 24 were weak cycles at the bottom of the 80-120 years longterm Wolf-Gleissberg solar cycle (Yousef 2006), (Mawad 2017), (Yousef et al. 2018), and (Mawad and Abdel-Sattar 2019). 2015) and (Farid et al. 2020). Solar activity has the greatest impact on AO density in the ionospheric layers leading to the enhancement of the erosion depths for material surfaces directly exposed to AO. The density and flux values of AO are higher during maximum solar activity than at minimum. The average fluence of AO per year and the corresponding erosion depth varies in response to solar-activity variations within the solar cycle (Samwell 2014) and (Farid et al. 2015) as shown in Fig. 1. In general, AO, hence O+, are the dominant species in the LEO environment. It is anticipated that as solar-flare events reach Earth, they enhance O and O+. Oxygen atoms have high corrosive power during and after combining with the material. According to the satellite's orbital velocity of 7.8 km/s in LEO, the satellite is exposed to very strong streams of AO at energies in the range of 5 eV (Dooling and Finckenor 1999) and (Mahmoud et al. 2021).

Coronal Mass Ejections
Coronal Mass Ejections (CMEs) are powerful eruptions of magnetic flux and plasma from the Sun into interplanetary space (Liu et al. 2020). The most powerful CMEs and their associated flares have very strong impacts on the near-Earth environment, affecting the lifetime of space technology, (Baker et al. 2004) and (Mawad et al. 2014). Many efforts should be made to develop new models for forecasting CMEs, (Bobra and Ilonidis 2016) and (Inceoglu et al. 2018).
Predicting the lifetime of the LEO satellites is of great concern to the satellite industry and technology. All LEO satellites suffer orbital decay due to their interactions with the Earth's atmosphere. The single-exponential smoothing method and the single moving-average method were applied to the sample as used in Aknil (Hudaningsih et al. 2020) and (Khodairy et al. 2020). The two methods were used to compare the most accurate forecasting methods close to their actual values. Previous studies used empirical methods for prediction of the travel time of interplanetary Coronal Mass Ejection Shocks (ICME) such as (Youssef et al. 2011). Also, many studies used the empirical methods to investigate the CME occurrence in accordance with the solar flare. Other studies used artificial-intelligence networks to detect the arrival time of interplanetary coronal mass ejection shocks during solar cycles, . For this study, to present numerical predictions for CMEs characteristics empirical methods were used.
The main objective of this study is to develop a standalone predictive model for the space-environment parameters and compare the results with real data.

Methodology and data sources
Statistical analyses and forecasting techniques depend on gathering historical data and finding the pattern (trend or seasonal variation) and determining the forecasting method. Following forecasting calculations, a verification step is required to assure the accuracy and precision of the prediction model.
The Holt-Winters method is used by many companies to produce short-term forecasts as their sales data contain a trend or seasonal pattern. This method is simple, easily automated, and has low data-storage requirements. The Holt-Winters method was applied to the air-transportation industry. The time series was decomposed into three additive components: trend, seasonal, and remainder. Multiple regression may be recognized as a better forecasting method for daily and weekly short-term forecasting, whereas Holt-Winters methods presented better forecasting values for monthly and yearly long-term forecasting (Tratar and Strmčnik 2016). Due to strong seasonal variation detected in the series, the Holt-Winters method can be used (Dantas et al. 2017).
Our data source is obtained from the Naval Research Laboratory for Mass Spectrometer and Incoherent Scatter Radar (NRLMSISE-00) Atmosphere Model (5) during the period 1996 to 2020 for the atmosphere. The data includes sunspot number, density, temperature, and flux of H-atoms cm −3 , He-atoms cm −3 , N-atoms cm −3 , oxygen atoms cm −3 , AO atoms cm −3 , atomic N-atoms cm −3 , and Ar-atoms cm −3 . The altitude resolution is 50 km and the date step size is 7 days, i.e., three observations per month. The data under study covers the period of solar cycles 23 and 24.
Segregation of data by altitude levels is then done by skipping altitudes under 100 km out of the LEO region. The observations are also taken at latitude 0.0°from the Earth's surface and from 0 km to 1000 km above the Earth's surface. The specified location [equator] at the Earth's surface is due to the geoid shape of the Earth's atmosphere. In other words, the greater the atmospheric altitude, the more space hazards occur. The equator has the highest locations in the atmosphere. The data time series was plotted.
Data regression is required. Without equations governing the variation of atmospheric properties and interpolation with time, the best way is to apply numerical analyses methods to drive the regression for such properties. The building and developing of the model depend on about 60% of the data, while the remainder 40% are used in its verification.
Careful analyses of the data are required. The data are plotted in Fig. 2, showing seasonal variations. It shows the correlation relation between the solar activity and AO flux for altitudes 100 km, 500 km, and 1000 km. The relations between sunspot number (SSN) and the property of the atmosphere reveals that the variation is the solar cycle. The visual inspection of graphs indicates the strong relation between AO density and the solar activity for altitudes 500 km and 1000 km, while for 100 km, the relation became more stable during the solar cycle except for minimum solar activity. In other words, the AO flux is more favorable for lower altitudes, The solar cycle is mainly about 11 years [about 132 months], i.e., every 132 rows, there is a new cycle. Thus, the regression method must be a seasonal method. There is a variety of seasonal regression methods such as Holt-Winter's seasonal additive method, Holt-Winter's seasonal multiplication method, exponential regression, and the Triple-Exponential Smoothing method.
On trying most of these methods, the chosen one is the Triple-Exponential Smoothing method, depending on the lowest error value. The basic equations of this method are: We define: y is the observation -S is the smoothed observation b is the trend factor -I is the seasonal index -F is the forecast at m periods ahead t is an index denoting a time period It is evident from both Fig. 3 and Fig. 4 that cycle 25 is likely to be almost as active as cycle 23.
The prediction model for atmospheric density, temperature, and AO fluxes is based on differential equations as follows:

Atmospheric Pressure (P )
where, ρ is the density, g is the gravitational constant, h is the altitude in the atmosphere, F is the net force, A is the cross-sectional area, R is the universal gas constant, and T is the temperature where, v a is the average thermal speed, v rms is the rootmean-square speed, and v p is the most probable speed of an air molecule (m s −1 ). These speeds are related to absolute temperature.

Atomic
for density, atomic oxygen, and temperature, respectively, at 1000 km, where x is the SSN number and ρ is the AO density

Results and discussion
In order to initialize the Triple-Exponential Smoothing method, we need at least one complete set of a solar cycle's data to determine initial estimates of the solar cycle's indices I t−L . All the above-discussed methods have been used for forecasting atmospheric constituents (He, Ar, O molecules, AO, N molecules, N atoms, and/or H), solar irradiance, orbital lifetime (semimajor axis, eccentricity and/or density drag), coronal mass ejections (linear speed, central PA, mass and/or angular width), debris hazards, the orbital visualization as well as density and temperature, T . The results of the software (S/W) are verified by comparison of the predicted data from S/W with the real data from NRLMSISE-00. The last date in the database of the S/W is 01/09/2020 and the verification was done on 01/06/2021. RMS and the correlation coefficients are calculated for ρ, T , and H-atoms as shown in Table 1. Table 2 displays the crosscorrelation between real and forecasted data for CMEs' position central angle PA (deg), angular width (deg), linear speed (km/s), and measurement position angle MPA (deg). PA is the four categories of position angle 90°, 180°, 270°, and 360°.The optimum correlation and an optimum number of lags are shaded. The lag is equal 7 days. Detailed comparisons for CMEs' central PA (deg), angular width (deg), linear speed (km/s), and MPA (deg) are displayed in Figs. 5,6,7,and 8,respectively. For AO forecasting and comparison with real data, see Table 3 and Figs. 9, 10, and 11 for altitudes 100, 200, and 300 km, respectively. The rest of the altitudes (from 350 km to 1000 km) do not give a strong correlation.
The high correlation-coefficient values for lower altitudes indicates the presence of atomic oxygen in abundance at such altitudes, in contrast to its presence in limited quantities at higher altitudes.

Conclusion
Solar cycle number 25 is likely to have about the same activity as cycle 23. However, unexpected strengthening of solar cycle 25 is not excluded. If that is the case, then the applied predictive models should be modified accordingly to fit the observations. AO is the most abundant species for the near LEOs. The corresponding AO-flux-solar-activity correlation coef- The altitudes from 350 km to 1000 km do not give a strong correlation for the AO fluxes forecasting, an advantage that makes them more favorable for LEOs due to the harmful corrosive effects.
The model of AO forecasting needs more modifications and refinements by adding more parameters and other trending and seasonal variations.

Recommendation
This research recommends continuous development and improvement of predictive models for solar activity and space-environment components. In addition, many S/W verifications are required.