Operational Dst index prediction model based on combination of artificial neural network and empirical model

In this paper, an operational Dst index prediction model is developed by combining empirical and Artificial Neural Network (ANN) models. ANN algorithms are widely used to predict space weather conditions. While they require a large amount of data for machine learning, large-scale geomagnetic storms have not occurred sufficiently for the last 20 years, Advanced Composition Explorer (ACE) and Deep Space Climate Observatory (DSCOVR) mission operation period. Conversely, the empirical models are based on numerical equations derived from human intuition and are therefore applicable to extrapolate for large storms. In this study, we distinguish between Coronal Mass Ejection (CME) driven and Corotating Interaction Region (CIR) driven storms, estimate the minimum Dst values, and derive an equation for describing the recovery phase. The combined Korea Astronomy and Space Science Institute (KASI) Dst Prediction (KDP) model achieved better performance contrasted to ANN model only. This model could be used practically for space weather operation by extending prediction time to 24 h and updating the model output every hour.


Introduction
Large-scale interplanetary disturbances from the Sun interact with Earth's magnetic field, resulting in severe space weather events, such as geomagnetic storms (Gonzalez et al., 1994;Ohtani et al., 2000;Bhaskar & Vichare, 2019). The Dst index is a representative index of geomagnetic activity in the space weather community (Sugiura, 1964;Rangarajan, 1989;Daglis et al., 1999;Wanliss & Showalter, 2006), and there have been numerous attempts to predict the Dst index (Lyons, 1998;Birn et al., 2001;Raeder & Maynard, 2001;Rastätter et al., 2013;Eastwood et al., 2017). Empirical and artificial neural network (ANN) models have been handled as important methods for predicting the Dst index (Murayama, 1982;Feldstein, 1992;Boyle et al., 1997;Kugblenu et al., 1999;Weigel, 2010). Many previous studies established that the southward component of the interplanetary magnetic field (IMF) presents the main cause of geomagnetic storms (e.g., Burton et al., 1975;Gonzalez et al., 1994;Echer et al., 2005). Burton et al. (1975) described the Dst index using the following empirical formula: where Dst* is the Dst index corrected by solar wind dynamic pressure. The term Q(t) is a function that expresses the rate at which the ring current is intensified by the dusk-ward solar wind electric field (represented in the geocentric solar magnetospheric coordinate system); P dyn is the solar wind dynamic pressure, s is the decay time when the ring current weakens by Topical Issue -Space Weather research in the Digital Age and across the full data lifecycle 1/e, and b and c are coefficients determined by observational data. As Q(t) is determined solely by the solar wind, we can accurately predict the Dst index using these formulas if we can observe or predict solar wind states. Many Dst models have been derived from this model (Lundstedt et al., 2002), and there have been numerous studies based on the above empirical formula (O'Brien & McPherron, 2000;Wang et al., 2003) since Burton et al. (1975). The model for predicting the Dst index 1 h ahead has been significantly improved by Temerin & Li (2002. In the 1990s, advances in computer technology led to ANN for predicting the Dst index (Freeman & Nagai, 1993;Lundstedt & Wintoft, 1994;Watanabe et al., 2002). These studies, which used solar wind data as input, initially utilized a simple feedforward ANN. Wu & Lundstedt (1996) applied a more complex Elman recurrent neural network. More data from solar wind observatories, Wind, and Advanced Composition Explorer (ACE), have been archived and used to improve the results obtained from neural networks. Currently, we can predict the Dst index 1 h ahead with relatively high accuracy. However, it is not clear whether this is sufficient for space weather operation, as it is not easy to determine the forecast time required for space weather applications. Obviously, more than 1 h of forecast time is required for space weather applications, changing scheduled satellite operations and the air routes in a polar region, and sending a warning signal to a complex communication system.
As Earth's magnetic field is disturbed almost simultaneously with solar wind variations, it is not easy to predict long-term changes in Earth's magnetic field using only the solar wind data observed at Lagrange point, L1, without forecasting solar winds. Lazzús et al. (2017) created a model that combined ANN with particle swarm optimization and predicted the Dst index 1-6 h ahead. Gruet et al. (2018) combined a long short-term memory recurrent neural network with a Gaussian process to predict the Dst index up to 6 h ahead. Naturally, the prediction accuracy tends to decrease sharply as prediction time increases. In this study, we also tried to extend the leading time of prediction up to 24 h differently.
The ACE mission, launched in 1997 and replaced by the Deep Space Climate (DSCOVR) satellite in 2015 for real-time operations, collected considerable solar wind data that can be used by ANN algorithms to predict the Dst index. Nevertheless, it is doubtful whether sufficient data is available for artificial intelligence and machine learning. Since geomagnetic storms with the Dst index of less than À100 nT occur only a few times a year during the solar maximum and significantly fewer times during the solar minimum, more than 20 years of solar wind data appear to be insufficient for predicting geomagnetic storms (Watari, 2017;Makarov, 2018). This paper presents three Dst prediction models, simple ANN, empirical model, and combination model (Korea Astronomy and Space Science Institute (KASI) Dst Prediction, KDP). The simple ANN is the baseline model for comparison. Here, we just tried to follow the existing neural network algorithm, so there are no significant advances and show a comparable prediction performance to the previous models, such as shown in Section 3. The empirical model (equations) indicates the intensity of a magnetic storm according to the degree of solar wind changes. Two kinds of geomagnetic storms are considered separately; coronal mass ejection (CME) driven and corotating interaction region (CIR) driven storms in Section 4. This preemptive information can help improve Dst prediction accuracy because a space weather operator can determine in advance whether CIR or CME causes future geomagnetic storms. Finally, the present study proposes a combination of empirical and ANN models in Section 5. The predicted values by the empirical model are used as an input parameter to the ANN model. Thus we could improve the prediction performance better than other models.
The method combining empirical and ANN models is similar to that of Revallo et al. (2014Revallo et al. ( , 2015; the difference is in the empirical model. They used time histories of solar wind and magnetosphere interaction proposed by Romashets et al. (2008) a neural network input data. They constructed an analytical representation of magnetic fields in the region where the solar wind interacted with the earth's magnetosphere and used the results as an input parameter for predicting the Dst index. On the other hand, we directly derived the Dst index's prediction values from the empirical equations and applied them to the ANN model. Thus, the KDP model could extend the prediction time up to 24 h while Revallo et al. (2014Revallo et al. ( , 2015 predict the Dst index just 1 h ahead. See and compare the model performance in Figure 9.

Data
For the development of a model, we use the solar wind data collected by the ACE and DSCOVR satellite and the Dst index as input values for the prediction model from 1999 to 2017. The Dst index is obtained from NASA (https://omniweb.gsfc.nasa. gov/ow.html), and the hourly average solar wind data are acquired from CDAWeb (https://cdaweb.gsfc.nasa.gov/index. html). For the model's operation, current Dst and solar wind data are downloaded from World Data Center for Geomagnetism, Kyoto (http://wdc.kugi.kyoto-u.ac.jp/dst_realtime/presentmonth/ index.html) and NOAA (https://services.swpc.noaa.gov/products/ solar-wind/). Tables A.1 and A.2 in Appendix list the CME-and CIRdriven geomagnetic storms from February 1999 to September 2017. Here, the geomagnetic storm is defined with a minimum Dst index of less than À50 nT. The event start time is defined as the maximum Dst time at the main phase, and the event end time is the recovery time exceeds À30 nT. These storm lists are mainly obtained from published papers (referenced therein), and some CIR-driven storms are added in this study. The effects of CME and CIR are known to be different (Liemohn et al., 2010). CME can potentially generate large geomagnetic storms, while CIR generally triggers minor storms. The empirical model developed in this study distinguishes the CME and CIR-driven storms, and the empirical model results are used as the input parameter of the combination model.

Artificial neural network (ANN) model
This section introduces a simple ANN model and shows its performance. The model is combined with an empirical model in Section 5. The ANN used in this study is the feed-forward network (Gardner & Dorling, 1998) known as the most simple ANN-type algorithm (Haykin, 1998). The training technique of weights and biases uses an error backpropagation learning algorithm with a gradient descent learning method. The basic structure of this ANN is composed of an input layer (I), a hidden layer (H), and an output layer (O) as shown in Figure 1. The hidden layer uses a nonlinear hyperbolic tangent function, H = tanh(I), and the output layer uses the O = H type linear function. The following equations express these layers: where H j and O k are the nodes of the hidden and output layers, respectively, wh i,j , wo j,k are the weights of the hidden and output layers, and b i and b k are biases, respectively. The suffixes i, j, and k denote the node numbers in the input, hidden, and output layers, and N and M denote the total number of nodes in the input and hidden layers, respectively. We use the input parameters listed in Table 1 for the input layer. Note that the simple ANN model uses the data collected for storm and quiet time and does not distinguish CME and CIR-driven storms. We use the current observations, and the differences between the current observations and the observations obtained 1 or 2 h earlier. The use of these differences as input data helps predict the Dst index trend and, therefore, plays an important role in predicting the Dst value more than 1 h ahead. The differences are estimated based on hourly averaged solar wind data. We also predict 1-24 h later Dst index as the target values of the output layer. The output layer has one node, as shown in Figure 1, which means the simple ANN model consists of 24 models for each hour predictions in advance. We train the model by adjusting hyperparameters such as the node numbers of the input and hidden layers, learning rate, and learning cycle to optimize the neural network. The optimization selects the best condition by calculating the correlation coefficient and root mean square error (RMSE) between the predicted and actual values for all hyperparameters. The equations are given as follows: where Y pre is a predicted value, Y real is a measured value, and N is the number of values.  Parameter Description N(tc) Current solar wind density N(tc) À N(tc À 1) Difference between current solar wind density and value obtained 1 h earlier V(tc) Current solar wind speed V(tc) À V(tc À 1) Difference between current solar wind speed and value obtained 1 h earlier Bt Current total IMF Bt(tc) À Bt(tc À 1) Difference between current total IMF and value obtained 1 h earlier Bz(tc) Current IMF Bz Bz(tc À 1) IMF Bz obtained 1 h earlier Bz(tc À 2) IMF Bz obtained 2 h earlier Bz(tc) À Bz(tc À 1) Difference between current IMF Bz and value obtained 1 h earlier Bz(tc À 1) À Bz(tc À 2) Difference between IMF Bz values obtained one and 2 h earlier Dst(tc) Current measured Dst index Dst(tc À 1) Measured Dst index values obtained 1 h earlier Dst(tc À 2) Measured Dst index values obtained 2 h earlier W. Park et al.: J. Space Weather Space Clim. 2021, 11, 38 For performance optimization, all hidden layers with node numbers of 7-50, a learning rate of 0.0001-0.01, and learning cycles of 1200-1500 were trained and validated (see Table 2). For the ANN training, we used the data observed in , 2007, 2011, 2012, 2016. For the model validation, we use 2000, and for the test, we use 2003,2005,2009,2015. A total of 53,839, 28,526, and 25,421 data points are used for training, validation, and testing, respectively. The data covered storm and quiet times.
The red dots (ANN model) in Figure 2 show the RMSE and the correlation coefficients obtained by comparing the prediction results with the observed Dst index during the test period. The statistical analysis is conducted for the entire period without dividing it into the storm and quiet period regardless of distinguishing between CME-and CIR-driven geomagnetic storms. Compared with other models, this figure also shows the results from the various ANN models reported in the literature (Lazzús et al., 2017). As seen in Figure 2, many models predict the Dst index just 1 h ahead, and more advanced models (Wu & Lundstedt, 1997;Stepanova et al., 2005) can predict up to 6 h, while our model attempts to predict the Dst index up to 24 h. The initial version of the ANN shows similar or slightly lower performance than other models because no improved techniques are used.

Empirical model
The relationship between the solar wind and the Dst index could be derived using an arbitrary empirical formula (McPherron & O'Brien, 2001). In deriving the equations, we regarded a geomagnetic storm as having only one main phase and one recovery phase to simplify the relationship between the measured Dst index and input parameters. However, this model could also be used to predict the multi-peak storms when combined with the ANN model. The empirical model consists of three equations: predicting the minimum Dst index, finding the time reaching Dst Minimum, and expressing the recovery phases.
Equations (7) and (8) are used to determine the minimum Dst index (Dst min ) from the current solar wind condition caused by CME and CIR, respectively. While Burton et al. (1975) expressed the ring current injection rate in terms of the solar wind speed (V) and the southward IMF (Bs), in this work, we add the solar wind density (N), total IMF (Bt), and current Dst index to complete the equation. The solar wind parameters on the right-hand side of equations (7) and (8) are measurements when the model starts to run each hour. Recall that this model is designed to produce new prediction values every hour for the space weather operation. These equations were derived by a polynomial fitting algorithm in commercial software for the storm period and added exponential terms for adjusting large storm trends. In deriving these equations, we fixed some arbitrary coefficients and fit a function to the measurements to find the unfixed coefficients. We do not think we found the optimum values describing minimum Dst because the equations have six input parameters. Note that there is still a large inconsistency between prediction and measurements in Figure 3. Nevertheless, these empirical equations effectively build the combination model, as shown in Section 5. At the initial storm phase, the predicted Dst min is mainly controlled by solar wind measurements because Dst current is generally a small value, and the term, e jDst current j 450 does not significantly contribute to the Dst min prediction. However, as a big storm progresses, only solar wind data are not enough to describe geomagnetic storm developments. The exponential term adjusts the estimation of Dst min . Figure 3 shows the minimum Dst predicted by using equations (7) and (8) (blue symbols), measured minimum Dst (black symbols), and trend (red line) of equations (7) and (8) under specific conditions, Bt = Bs, V = 500 km/s, and N = 30/cm 3 . In these figures, the predicted Dst min is calculated at the initial phase of the storm when V Â Bs is just over 0.49 mV/m, and several hours later, minimum Dst measurement is obtained. As shown in Figure 3, the frequency of smaller geomagnetic storms is higher. For the prediction algorithms based on artificial intelligence, a larger weighting factor is applied to more frequent storms. Therefore, the ANN algorithm shows the worse prediction accuracy for the greater storms. In contrast, the empirical models are derived by human intuition and attempt to find a correlation between minimum Dst index and input parameters, regardless of the amount of data.
Next, we calculate the time (T min ) required to reach the minimum Dst obtained from equations (7) and (8). If we could estimate the Dst variation (dDst) for 1 h, we can estimate T min as the following equation: As the Dst index represents the variation in the geomagnetic field due to ring current development, the amount of change in the Dst index is proportional to the injection rate of the ring current particles. In this study, we adopt Wang et al. (2003) empirical model to derive the rate of change in Dst per hour (dDst) during the storm main phase.
Dst Ã ¼ Dst À 7:26 ffiffiffi P p þ 11 nT; where P (in nPa) is the solar wind dynamic pressure and P 0 and c are parameters provided by Wang et al. (2003). These parameters depend on the Dst index, as shown in Table 3. Here, the Dst index refers to a measurement at the current time.
For more information about the model, refer to Wang et al. (2003). Now, we obtain the dDst, and the prediction equation for the main phase can be expressed as follows:   (7) and (8) with Bt = Bs, V = 500 km/s, and N = 30/cm 3 . W. Park et al.: J. Space Weather Space Clim. 2021, 11, 38 Here, we assumed the Dst index decrease linearly to Dst min from Dst current . While the Dst index changes during the main phase can be predicted using equation (11), the recovery phase's equation is required to be derived. Figure 4 is a superposed epoch plot of the Dst index for the storms' recovery phases listed in Tables A.1 and A.2 (Appendix). Here the minimum Dst is normalized to À1. We examined several fitting functions to derive an equation describing the Dst index trends for the recovery phase. We found that an equation represented by two exponential functions, such as equation (12), effectively follows the recovery phase trend. Here, equation (12) presented by the red line is not the averaged Figure 4 but just a plausible expression derived by intuition. Other researchers can get a different equation.
From equation (11), we can get the profile of Dst variation for the storm main phase and from equation (12) for the recovery phase from 1 h to 24 h ahead. During a storm is progressing, Dst min is updated every time, and if the predicted Dst min is less than the current Dst, the empirical model recognizes the storm is in the recovery phase. Otherwise, if the Dst min is larger than the current Dst, then the algorithm predicts another storm main phase. Thus, this model can generate multi-peak storm prediction. These values are used as input parameters of the ANN model.

Combination model
The KASI Dst prediction (KDP) model effectively predicts the Dst index by combining the empirical and ANN models. First, we estimate the Dst index from 1 h to 24 h ahead using empirical equations (11) and (12). The prediction results from the empirical model are combined with the measurements listed in Table 1. For example, the 6th ANN model's input parameters are the 6 h ahead prediction value from equations (11) and (12) and the measurements in Table 1. Hence, it is as if there are 24 independent models for predicting Dst for each hour. Even if the predicted Dst from the empirical model is just one among 15 inputs, note ANN gives greater weight to the empirical predictions. Therefore, just adding the empirical model outputs on the ANN's input data improves the final prediction performance. The method of training and validation is the same as the original ANN model. Table 4 shows the combination model's hyperparameters. Figure 5 shows an example of how the KDP model predicts the Dst index during geomagnetic storms. A halo CME was observed at 08:06 UT on 18 November 2003 by Solar and Heliospheric Observatory (SOHO) Large Angle and Spectrometric Coronagraph (LASCO) C2. At this time, the CME's speed was approximately 1150 km/s, and it would reach the Earth after approximately 32 h by simply dividing the Sun-Earth distance (150,000,000 km) by the solar wind speed (1150 km/s). However, there were big uncertainties at this moment whether the geomagnetic storm would occur by this CME, how much the Dst index would decrease, how fast it would reach the minimum Dst, and how long it would take recovery. At this moment, space weather operators set the KDP as CME model and wait for the CME shock arrival. The ACE satellite observed the CME shock at 07:26 UT on 20 November, indicating that the arrival time was approximately 47 h, approximately 15 h later than expected. However, the model did not predict geomagnetic storms because the IMF B z component was northward at the shock arrival. In the bottom panel of Figure 5, each color dots indicate the Dst values predicted 1 h at each hour, and the lines show the following predicted values for 24 h with one-hour interval. Note the first red dot indicates around zero, and no significant Dst changes are expected for 24 h. As time goes on, the colored dots and lines show the magnetic storm proceeding.
As the IMF B z component rapidly turned southward at 10:45 UT on 20 November, the KDP model started to predict a minor storm event, where a Dst minimum of À140 nT would be reached around 13:00 UT. As the IMF enhanced southward Table 3. Parameters which used for calculating the ring current injection rate, Q.

Dst index
P 0 c > À50 nT 3.3 0.2 À100 nT < Dst < À50 nT 3.1 0.19 À150 nT < Dst < À100 nT 3.5 0.25 < À150 nT 3.5 0.18 and the solar wind speed increased, the KDP model produced new predictions over time and provided results similar to actual Dst values. At 12:00 UT, the KDP model predicted that a Dst minimum of À430 nT would be reached around 20:00 UT. This value was similar to the measured value of À422 nT at 21:00 UT. At 21:00 UT on 20 November, the KDP model predicted that the storm would recover to À100 nT after one day. In this manner, space weather operators can deliver information about space weather conditions to users several hours before the minimum Dst is reached. Figure 6 shows another example of KDP model prediction where the magnetic storm occurred by CIR on 5 August 2019. Generally, CIR could be identified by the corona holes on the solar UV images and produces minor magnetic storms. From 03:00 UT on 5 August, the solar wind speed increased slowly to about 700 km/s, similar speed to the CME storm shown in Figure 5. However, the IMF B z was larger than À10 nT during the storm was progressing and caused different scale storms from the CME-driven storm described in Figure 5. At the initial stage of the storm, the KDP model predicted the Dst index well. However, the solar wind speed increased extraordinary, and the model overestimates the minimum Dst with À80 nT while the measurement was larger than À50 nT. Nevertheless, the model successfully predicted the minor storm caused by CIR.
Regarding the KDP model performance, Figure 7 shows the correlation coefficient and RMSE derived from equations (5) and (6)     model for space weather operation. In other literature, the space weather model's performance has been shown for both storm and quiet periods. These algorithm performances would be governed by the quiet period because it is considerably longer than the storm period. However, from a space weather operator's perspective, the important thing is how storm events are accurately predicted. We believe a better operational model should have good performance in predicting storm events. Thus, these results are obtained for the storm events in the test period. Figure 8 shows the KDP model's correlation coefficients and RMSEs for training, validation, and test group to illustrate the model's prediction capability. The higher correlation coefficient and smaller RMSE show that the Dst index is predicted more accurately using the combination model than the empirical and ANN models.

Discussion and conclusion
This paper introduces the KDP model that combined an ANN and an empirical model. The combination model showed better performance than each model in predicting the Dst index from 1 h to 24 h ahead, as shown in Figure 7. Therefore, we conclude that better space weather prediction can be achieved by using the combined models.
While this study aims to show an ANN model's improvement by combining an empirical model, as shown in Figure 7, there have been demands for comparing the KDP model's performance with others. While the ANN in the KDP model is trained only during the storm periods, we can find the input parameters even in quiet time. Thus, the KDP model can produce non-storm predictions. Figure 9 shows the storm and quiet time results, implying the combination model has superior prediction performance than other developed models. However, such a comparison does not tell us which model is better or worse because the test periods and conditions are different. The persistence model shown in Figure 9 assumes the prediction values are the same as the current Dst. Figure 9 clearly shows the KDP model's performance is better than the persistence model, and this model might be useful in space weather operation.
In addition to the persistence model, prediction efficiency is also used to validate the space weather model. Figure 10 shows the prediction efficiency of the KDP model calculated as below: W. Park et al.: J. Space Weather Space Clim. 2021, 11, 38 where Dst s , Dst nn s , and Dst stand for the measured, predicted, and mean Dst values respectively, and M is the length of the record. The prediction efficiency tests the ability of the model to predict the variation of Dst around the mean. A prediction efficiency of 1 is perfect agreement at all times. Prediction efficiencies less than or equal to zero do not provide useful predictions of the observations' time variation. Figure 10 shows the prediction efficiency decrease as prediction time increases, while the model is still useful in predicting the 24-hour Dst index.
It should be noted that the ANN and empirical models used in this study are developed with a conventional method. The ANN is a commonly used feed-forward network, and roughly estimated equations implement the empirical models. As shown in Figure 2, the ANN model is not superior to other existing models. As shown in Figure 3, the minimum Dst values are spread widely, and the empirical equations are not fitted well to the data. Nevertheless, the combined model yielded better results, which means the combination model could be significantly improved by adopting modernized ANN and sophisticate  empirical models. Further research will improve the empirical equations for predicting the Dst index. If two improved models are combined, the prediction accuracy would be more increase. Thus, this paper encourages researchers to improve their neural network models by adopting empirical models.
Besides, by categorizing geomagnetic storms into CMEdriven and CIR-driven storms, we improved the model performance. Through such categorization, we can increase the correlation between the Dst index and input parameters. In a future study, the solar cycle and seasonal effects could also be considered to make an improved empirical model. Space weather operators can easily recognize whether CME or CIR will cause magnetic storms from solar images. However, sometimes CMEs could be embedded in a CIR and resulting in a compound storm event (Al-Shakarchi & Morgan, 2018). In this case, the decision for selecting CME or CIR mode remains to the operators. Even if they fail in identifying CME or CIRdriven storms, the KDP model predicts the Dst index anyhow by solar wind and current Dst.
The KDP model has been implemented to predict the Dst index 24 h ahead, and the prediction output is updated every hour for space weather operations. In general, space weather information is delivered to a professional group of users, such as radio communication companies, satellite operators, and flight dispatchers. Space weather users state that they can cope with severe space weather conditions in reasonable prediction time. For example, flight dispatchers (Privet communication with Korea air dispatchers) say they can change routes 3 h before a flight. It is well known that satellite launches can be delayed just a few minutes before launch. GEO satellite operators monitor their satellites for 24 h, and they can act quickly for emergencies. As shown in Figure 7, the KDP model can predict the Dst storms 6 h ahead with a correlation coefficient of 0.8 and an RMSE of 24 nT or less. Figure 11 shows the correlations between measurements and prediction values for 3-hour, 6-hour, 12-hour, and 24-hour ahead. The KDP model shows good correlations up to 6-hour ahead. We emphasize these are acceptable characteristics for space weather operation. Thus, this model can be used as an operational model for delivering warning signals to space weather users.
We have noted that space weather users do not want to know only the minimum Dst index but the recovery time of storms. Users like to know when they can return to their daily lives. Even though the 24-hour prediction correlation coefficient of the KDP model is not satisfactory at 0.38, we attempted to develop a model to track the storm's progress using the 24-hour prediction model and predicted the recovery phase. While the Dst index's recovery does not mean that a space storm event has ended completely, it would be helpful information for space weather alerts.
W. Park et al.: J. Space Weather Space Clim. 2021, 11, 38 downloaded from World Data Center for Geomagnetism, Kyoto (http://wdc.kugi.kyoto-u.ac.jp/dst_realtime/presentmonth/index. html) and NOAA (https://services.swpc.noaa.gov/products/ solar-wind/). This work was supported by the National Meteorological Satellite Centre (NMSC) of the Korea Meteorological Administration (KMA) through the "Geostationary Meteorological Satellite Ground Segment Development" research project. The editor thanks two anonymous reviewers for their assistance in evaluating this paper.