Next Article in Journal
Characteristics of a Heat Exchanger in a Liquid Rocket Engine Using Conjugate Heat Transfer Coupling with Open-Source Tools
Previous Article in Journal
Adaptive Feed-Forward Control for Gust Load Alleviation on a Flying-Wing Model Using Multiple Control Surfaces
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Proton Flux Prediction Method Based on an Attention Mechanism and Long Short-Term Memory Network

1
College of Physical Science and Technology, Yangzhou University, Yangzhou 225002, China
2
China Institute of Atomic Energy, Beijing 102413, China
3
Beijing Institute of Tracking and Telecommunications Technology, Beijing 100094, China
4
National Space Science Center, Chinese Academy of Sciences, Beijing 100190, China
5
Shanghai SastSpace Technology Co., Ltd., Shanghai 201109, China
6
Yangzhou Polytechnic Institute, Yangzhou 225127, China
7
College of Electrical, Energy and Power Engineering, Yangzhou University, Yangzhou 225127, China
*
Authors to whom correspondence should be addressed.
Aerospace 2023, 10(12), 982; https://doi.org/10.3390/aerospace10120982
Submission received: 9 October 2023 / Revised: 13 November 2023 / Accepted: 15 November 2023 / Published: 22 November 2023

Abstract

:
Accurately predicting proton flux in the space radiation environment is crucial for satellite in-orbit management and space science research. This paper proposes a proton flux prediction method based on a hybrid neural network. This method is a predictive approach for measuring proton flux profiles via a satellite during its operation, including crossings through the SAA region. In the data preprocessing stage, a moving average wavelet transform was employed to retain the trend information of the original data and perform noise reduction. For the model design, the TPA-LSTM model was introduced, which combines the Temporal Pattern Attention mechanism with a Long Short-Term Memory network (LSTM). The model was trained and validated using 4,174,202 proton flux data points over a span of 12 months. The experimental results indicate that the prediction accuracy of the TPA-LSTM model is higher than that of the AP-8 model, with a logarithmic root mean square error (logRMSE) of 3.71 between predicted and actual values. In particular, an improved accuracy was observed when predicting values within the South Atlantic Anomaly (SAA) region, with a logRMSE of 3.09.

1. Introduction

When it comes to the safety and operational stability of satellite spacecraft, the significance of the space radiation environment cannot be underestimated. The radiation present in space, such as high-energy particles, solar storms, and cosmic rays, profoundly impacts various aspects of satellite spacecraft.
Numerous studies focused on assessing the impact of space radiation on satellite electronic components. For instance, as demonstrated by N Ya’acob et al.’s [1] research, high-energy particle impacts may result in Single-Event Effects (SEE) and Total Ionizing Dose (TID) in electronic elements. These events can lead to temporary or permanent malfunctions of electronic devices, subsequently affecting satellite performance and functionality. DJ Cochran et al. [2] extensively examined the sensitivity of various candidate spacecraft electronic components to Total Ionizing Dose and displacement damage.
In the context of investigating the impact of space radiation on satellite spacecraft, proton flux stands as a crucial component of space radiation that holds undeniable significance. Both its energy and quantity can exert notable effects on satellite electronic components, thereby influencing satellite performance and operational stability. Numerous scholars have conducted pertinent research in this realm. For instance, R Uzel et al. [3] explored the efficiency of localized shielding structures in capturing protons to meet reliability requirements. S Katz et al. [4] observed an elevated single-event upset rate in the low Earth orbit polar satellite Eros B as high-energy proton flux increased.
However, merely understanding variations in the radiation environment falls short; accurate predictions and warnings are equally crucial. The proton flux data during satellite operation represent a stable time series. Predicting this time series not only provides crucial data references for space science research but also offers vital insights for satellite operation and management, thereby preventing malfunctions and instability in high-radiation environments.
In this regard, a substantial amount of research has been devoted to exploring how to construct more accurate radiation environment prediction and warning methods. For instance, J Chen et al. [5] introduced an embedded approach for forecasting Solar Particle Events (SPE) and Single-Event Upset (SEU) rates. This approach enables prediction through the increase in SEU count rates one hour prior to the forthcoming flux changes and corresponding SPE events.
In the realm of space radiation environments, particularly in the prediction of proton flux, although machine learning methods offer a promising avenue, they also encounter a series of challenges. Firstly, the complexity of the data amplifies the difficulty of prediction tasks [6]. Secondly, the presence of various noises could disrupt predictive models, impacting the accuracy of forecasts [7]. Most significantly, the accuracy of predictive models stands as a key challenge [8].
In recent years, deep learning methods have made significant progress in handling time series data for this issue. These methods include Convolutional Neural Networks (CNN) [9,10,11], ARMA Neural Networks [12], and Long Short-Term Memory Networks (LSTM) [13,14,15], among others [16]. For instance, L Wei et al. [17] used a Long Short-Term Memory Network (LSTM) to develop a model for predicting daily > 2 MeV electron integral flux in geostationary orbit. Their model takes inputs such as geomagnetic and solar wind parameters, as well as the values of the >1 MeV electron integral flux itself over the past consecutive five days. Their experimental results indicated a substantial improvement in predictive efficiency using the LSTM model compared to some earlier models. The predictive efficiencies for the years 2008, 2009, and 2010 were 0.833, 0.896, and 0.911, respectively.
In the realm of deep learning, attention mechanisms stand as powerful tools [18,19]. These mechanisms allow models to focus on crucial time segments, enhancing accuracy in predicting future trends. In the field of space radiation environment prediction, the application of attention mechanisms is gaining traction. For instance, X Kong et al. [20] introduced a Time Convolutional Network model based on an attention mechanism (TCN-Attention) for predicting solar radiation. Their experimental results underscore the strong predictive performance of their proposed model, with a Nash-–Sutcliffe Efficiency coefficient (NSE) of 0.922.
Building upon the foundation of attention mechanisms, SY Shih et al. [21] introduced Temporal Pattern Attention (TPA), a technique with distinct advantages in time series prediction. TPA can be employed to identify pivotal time segments within proton flux time series which might correspond to specific intervals when satellites traverse certain spatial positions. By accentuating these key time segments, the model becomes more adept at capturing the changing patterns of radiation events, consequently enhancing predictive performance.
In the prediction of space radiation environments, the combination of wavelet transform and neural network models can yield excellent predictive results [22]. Data processed through wavelet transform can exhibit strong predictive performance in the realm of space radiation, enhancing the model’s accuracy in forecasting radiation events [23]. For instance, SS Sharifi et al. [24] employed Continuous Wavelet Transform (CWT) in conjunction with a Multi-Layer Perceptron Neural Network (MLPNN) to predict Solar Radiation (SR). The results showcased the outstanding predictive performance of the CWT-MLPNN approach, with the R2 coefficient increasing from 0.86 to 0.89 compared to using a standalone MLPNN.
This paper presents a hybrid neural network model for predicting proton flux profiles measured via a satellite during its operation, including crossings through the SAA region. The model integrates the data preprocessing, hybrid neural network prediction, and accuracy verification steps. The data are processed using moving average wavelet transform, and the model incorporates latitude and longitude features, along with a Temporal Pattern Attention (TPA) mechanism. For this paper, building upon the LSTM neural network, a hybrid neural network model was constructed, and our experimental results demonstrate that the TPA-LSTM model exhibits good predictive performance in forecasting proton flux profiles.

2. Methods

2.1. Overall Method and Evaluation Criteria

The objective of this study was to predict space radiation proton flux using a hybrid neural network model. The model integrates multiple techniques, as depicted in Figure 1, encompassing data preparation and data preprocessing, hybrid neural network prediction, and accuracy verification and discussion.
We initially divided the dataset into 26 cases, randomly selecting 15 for the training set while leaving the rest for validation. During the data preprocessing phase, we employed moving average wavelet transform, marking our first innovation. This involved calculating the trend component of the data using the moving average method and denoising and decomposing the data using wavelet analysis. Subsequently, we analyzed proton flux data and identified the South Atlantic Anomaly (SAA) region as an area of significant fluctuation. This became our entry point for introducing latitude and longitude features and attention mechanisms, our second innovation. Following this, we constructed a hybrid neural network model, encompassing two LSTM layers, a TPA layer, and a dropout layer. In this step, our third innovation was the adoption of the TPA-LSTM model. Finally, we designed a series of experiments to verify the model’s accuracy. We primarily focused on the model’s accuracy in predicting proton flux for all time steps, particularly emphasizing its accuracy in predicting within the SAA region. Through comparisons with other models, we demonstrated the TPA-LSTM model’s accuracy in forecasting space radiation proton flux and showcased its performance across various time series prediction tasks.
For this paper, we employed two evaluation metrics to assess the predictive performance of the model: logRMSE (Logarithm of Root Mean Squared Error) and logMAE (Logarithm of Mean Absolute Error). Given the substantial differences in magnitude between the proton flux data in this study, we used logRMSE and logMAE as evaluation metrics. logRMSE is calculated by computing the mean of the logarithmic differences between predicted values and observed values. This mitigates the impact of data magnitude differences, ensuring that errors at larger scales do not unduly influence the evaluation metrics. logMAE, on the other hand, first takes the logarithm of each prediction error and then calculates their averages. Logarithmic error can shift the focus of the evaluation process towards relative error rather than absolute error. The formulas for both metrics are as follows:
l o g   R M S E = ln ( 1 n × i = 1 n ( y i y ^ i ) 2 )
l o g   M A E = ln ( 1 n × i = 1 n | y i y ^ i | )
In Equations (1) and (2), n is the number of samples, y ^ i is the predicted value, y i is the actual value.

2.2. Data Preparation and Dataset Partition

The in-orbit satellite-measured proton flux data in this study were provided by Shanghai SastSpace Technology Co., Ltd. The energy range of high-energy protons was 100.00 MeV to 275.00 MeV. The satellite was identified as 52085, and its orbital elements are presented in Table 1. The satellite employs Si detectors, with a shielding effect equivalent to 1 mm of Al. This study considers the measured proton flux data for this satellite from April 2022 to March 2023.
We partitioned the dataset into 26 distinct time periods, providing the time period number, start time, end time, and duration for each time period, and all time periods involve crossings through the SAA region, as outlined in Appendix A Table A1. Among the 26 cases, a subset of 15 cases were randomly selected for training purposes, while the remaining 11 were exclusively reserved for validation in each predictive iteration. To mitigate the influence of inherent randomness stemming from the selection of training and testing sets, along with the variability in parameter configuration, we conducted a minimum of 50 independent runs for each specific method to accurately portray their predictive accuracy. All evaluation metrics were derived from the predictive data generated by the validation set and compared against actual observations.

2.3. Data Preprocessing

During the data preprocessing stage, we employed the moving average wavelet denoising method. This process involves calculating the data’s trend component using the moving average technique, applying wavelet transformation for noise reduction, and selecting appropriate wavelet bases and decomposition levels. The result is a processed time series dataset. The processing procedure is illustrated in Figure 2.
In the data preprocessing process, we initially employed the moving average method to compute the trend component of the data. The moving average method is a commonly used smoothing technique that reduces noise and retains trend information by calculating the mean within a sliding window. By calculating the trend component, we can better grasp the overall trend of the data and reduce interference in subsequent analysis and modeling. The formula for the moving average method is as follows:
T t = m / j = k k y t + j
In the above equation, y t represents the time series data, m = 2 k + 1 , meaning that the trend at time t is derived by calculating the average of y t over the preceding and succeeding k cycles. This method is referred to as “m-MA”, denoting an m-order moving average. To achieve an even-order moving average, we performed two consecutive odd-order moving averages.
We proceeded with wavelet transformation, harnessing the potent capabilities of this signal processing tool. During the satellite’s operational phases, it may traverse specific zones, such as the South Atlantic Anomaly (SAA), or encounter space weather events like solar proton events, resulting in sudden magnitudes of proton flux data increase. Wavelet transformation facilitates the retention of these crucial data points while enabling a beneficial decomposition, thus aiding the learning process of neural networks. We adopted the minimax thresholding approach of wavelet transformation to establish the thresholds. Data points with trend values computed through the moving average method that surpass the threshold were designated as noise. Rooted in the distribution traits of wavelet coefficients, this method gauges noise by contrasting the amplitude of wavelet coefficients with the threshold. It constitutes a simple yet effective signal processing technique, proficiently eliminating noise while retaining vital information. The threshold is determined by an empirical formula:
λ = { 0 , N 32 0.3936 + 0.1829 ln ( N 2 ) , N > 32 }
where N represents the length of the input sequence.

2.4. Long Short-Term Memory Network (LSTM)

Long Short-Term Memory (LSTM) is a variant of the traditional Recurrent Neural Network (RNN). In handling longer time sequences, conventional RNNs often encounter issues with long-term dependencies, leading to the omission of information between distant time steps and subsequently causing performance degradation. LSTM, however, exhibits enhanced memory capabilities, allowing it to learn and capture long-term dependencies effectively, thereby facilitating more efficient information transmission within the time sequence. When dealing with time sequences, LSTM outperforms, enabling the model to better capture intricate relationships within time sequence data. This improvement in capturing complex associations heightens the accuracy and efficacy of time series modeling.
Figure 3 illustrates the structure of the LSTM model comprising three LSTM units, each denoted as Unit A. Within each LSTM layer, the neural units receive inputs from the current input X t , the output h t 1 of the preceding LSTM unit, and the cell state vector C t 1 of the previous LSTM unit. Meanwhile, the output of the neural units encompasses the cell state vector C t of the current LSTM layer’s memory unit at time t, as well as the output h t of the current LSTM unit at time t.
The internal neural structure of the LSTM network is depicted in Figure 3. Each LSTM layer acts as a robust memory cell, regulating the flow and retention of information through a forget gate f t , an input gate i t , and an output gate o t .
The forget gate plays a crucial role in LSTM, determining which information from the current time step input X t and the output h t 1 of the previous LSTM unit should be retained and which should be forgotten. By passing through a sigmoid activation function, the forget gate outputs a value between 0 and 1. A value close to 1 indicates that the information will be retained, while a value close to 0 signifies that the information will be forgotten. This allows LSTM to selectively keep important information and discard less relevant details, effectively addressing long-term dependencies in time series and enhancing sequence modeling and prediction performance. The formula is as follows:
f t = σ ( W f [ X t , h t 1 + b f ] )
In Equation (5), W f is the weight of the forget gate, and b f is the offset of the forget gate.
The input gate in LSTM is utilized to update the cell state. Its calculation process is as follows: Firstly, the current input X t is passed through a sigmoid activation function, along with the output h t 1 of the previous LSTM unit, to obtain the control value i t . This control value’s output ranges from 0 to 1, where 0 signifies unimportant and 1 represents important. Subsequently, the current input i t and the output h t 1 of the previous LSTM unit are put through a tanh function to generate a new state information C ˜ . The value of C ˜ also lies between 0 and 1, representing the effective information of the current LSTM unit. The control value i t determines the degree to which the information in C ˜ is integrated into the cell state C ˜ . Then, the previous cell state C t 1 is updated to the new cell state C t . Initially, the forget gate f t is multiplied by the previous cell state C t 1 , allowing for selective forgetting of certain information. Next, the input gate it is multiplied by C ˜ to obtain the updated cell state vector C t . This process achieves the selective updating of the cell state, efficiently transmitting vital information. The formulas are as follows:
i t = σ ( W i [ X t , h t 1 + b i ] )
C ˜ = tan h ( W c [ X t , h t 1 + b c ] )
C t = f t × C t 1 + i t × C ˜
In Equations (6) and (7), W i and b i are the weight and bias of the input gate; W c and b c are the weight and bias of the control unit.
The output gate, denoted as o t , is employed to regulate the output h t of the current LSTM unit at the given time step. Initially, the extent of the output gate’s opening is computed using a sigmoid activation function, resulting in an output control value ranging between 0 and 1. Then, the new cell state C t obtained from the input gate is passed through the tanh function, producing an output between −1 and 1. Lastly, the output of the tanh function is multiplied by the output control value of the output gate to yield the final output h t . This way, the output gate can flexibly determine which components of the cell state C t will influence the output h t at the current time step. The formulas are as follows:
o t = σ ( W o [ X t , h t 1 + b o ] )
h t = o t × tan ( C t )
In Equation (9), W o and b o are the weights and offsets of the output gate.

2.5. Temporal Pattern Attention (TPA)

Temporal Pattern Attention (TPA) is an advanced attention mechanism designed for capturing temporal dependencies and patterns in time series data. The model structure is illustrated in Figure 4. In this study, TPA demonstrates the capability to selectively focus on variations in latitude and longitude features during the South Atlantic Anomaly (SAA) region, which represents time periods of greater significance for changes in radiation proton flux. This adaptability enhances the model’s precision in predicting variations in radiation proton flux.
TPA is typically connected after an LSTM layer. Initially, the hidden state matrix H = { h t w , h t w 1 , h t 1 } , where w represents the sliding window length, is obtained from the LSTM layer’s output. Sequence features are extracted by performing convolutional operations using CNN filters C = { C 1 , C 2 , C T } , with a total of T filters. The convolution operation between H and C results in the temporal pattern matrix H C . The formula is as follows:
H i , j C = l = 1 w h i , ( t w 1 + l ) × C j , T w + l
where H i , j C represents the result of the action of the i th row vector and the j th filter in the hidden state matrix.
Next, the similarity between time steps is calculated through a scoring function. In TPA, the attention mechanism employs a “Query-Key-Value” mechanism to compute attention weights. In this mechanism, Query is used to inquire about the relevance of specific time steps, Key is utilized to compute the similarity between different time steps, and Value encompasses the features of time steps. By employing a scoring function to compute the similarity between Query and Key, attention weights for each time step can be derived, thus achieving weighted aggregation across different time steps. In this context, the LSTM layer’s output h t serves as the Query, the temporal pattern matrix H C functions as the Key, the calculation of relevance between h t and H C is conducted, and normalization using the sigmoid function yields the attention weights α i . The formulas are as follows:
f ( H i C , h t ) = ( H i C ) T h t
α i = sigmoid ( f ( H i C , h t ) )
The third step involves weighted aggregation. Weighted aggregation refers to the process of computing the weighted sum of Values based on attention weights to obtain the final context vector V t . Here, Value signifies the hidden state information of each time step, which corresponds to the temporal pattern matrix H C . This weighted representation highlights significant time steps, thereby extracting crucial temporal patterns. The formula is as follows:
V t = i = 1 N α i H C
where N is the length of the sequence.
Lastly, the weighted representation obtained from TPA, V t , is integrated into further computations of the model. The influence of the context vector V t on the current output h t is incorporated into h t , resulting in the ultimate output h t . The formula is as follows:
h t = W h h t + W v V t
where W h is the weight of the current output state, and W v is the weight of the context vector.

3. Results and Discussion

3.1. Moving Average Wavelet Transform

Firstly, we applied the moving average method to obtain the trend component of the original proton flux data. To select an appropriate moving average order, denoted as m, we experimented with different orders, including 2, 3, and 5. As observed from Figure 5, the trend component curve derived using the moving average method effectively captures the primary directional trend of the time series. As m increases, the moving average curve gradually approaches a smoother state due to larger m values, leading to more data smoothing. To better preserve the fluctuation pattern of the original curve, we opted for a smaller m value, specifically m = 2. This approach assists in accurately capturing the features of the time series data and ensuring that the trend component reflects the changing trends of the data.
The second step involves denoising the trend component. We experimented with five different wavelet bases: Haar, db2, db3, db4, and db5. Since excessive decomposition levels can result in the loss of some original signal information, a one-level decomposition was applied to all five wavelet bases.
Table 2 displays the logRMSE and logMAE values after performing wavelet transformations using different wavelet bases. Upon comparison, it is evident that the db2 and db5 wavelet bases exhibit the best performance in terms of evaluation metrics, although the differences are not substantial. This indicates that these wavelet bases possess similar capabilities in denoising and capturing data trends. The strong correlation between the wavelet-transformed data and the original data highlights the ability of the transformed data to retain the fluctuations present in the original data.
By contrasting the processing outcomes of various wavelet bases, we obtained the results of the wavelet transformation, as depicted in Figure 6. Among these wavelet bases, we observed that excessively high levels of decomposition lead to information loss. For our dataset, employing the db2 wavelet base proved sufficient to achieve denoising effects while also retaining the essential fluctuation trends of the parent time series.

3.2. Data Analysis

Figure 7 illustrates the distribution of high proton flux values across all operating cases after data preprocessing, with latitude and longitude serving as identifiers. From the figure, it can be observed that the overall data predominantly occupies quadrants 1, 2, and 3, with the majority of data points located in quadrant 3, namely the longitude range of (−90, 0) and latitude range of (−60, 0). This distribution pattern can be attributed to the presence of the South Atlantic Anomaly (SAA) on Earth. The SAA’s geographic coordinates span from longitude (−100, 20) and latitude (−60, 10). Within this region, there is an abnormal increase in high-energy particle flux, leading to sudden spikes in the detected proton flux data. Consequently, it is evident that data points with high values within the SAA area significantly deviate from the normal data distribution, forming a relatively concentrated anomaly region.
The data for Figure 8 were derived from proton flux data from all datasets, with Figure 8a representing the World map of proton flux, and Figure 8b representing the 3D-View of proton flux. The plotting of Figure 8 is based on proton flux data detected by the satellite as it passes through certain points. Through Figure 8, we can observe the spatiotemporal distribution characteristics of proton flux throughout the satellite’s entire operational period. We observed a concentrated burst of proton flux within the SAA region, while proton flux in other regions remains relatively stable. This observation aligns closely with our previous analysis of the spatiotemporal distribution of proton flux.
Building upon these findings, we have incorporated latitude and longitude features into our future proton flux predictions. Additionally, we have introduced the Temporal Pattern Attention (TPA) mechanism, focusing specifically on the impact of the SAA region. This approach aims to enhance the accuracy and reliability of proton flux profile forecasting.

3.3. Construction and Optimization of TPA-LSTM Neural Network

The hybrid neural network architecture of this study, as illustrated in Figure 9, follows a sequential arrangement: data input layer, LSTM layer, TPA layer, dropout layer, LSTM layer, and output layer. The data input layer takes longitude, latitude, and proton flux feature values as inputs into the network. The first LSTM layer captures temporal dependencies in the time series data, handling the time series data of proton flux feature values. Following the LSTM layer, the TPA layer is introduced and applied to the proton flux feature values. A dropout layer is added after the TPA layer to prevent overfitting. The second LSTM layer enhances the prediction accuracy of proton flux. Finally, an output layer is set up to perform predictions on the proton flux feature values.
In this neural network model, longitude and latitude feature values serve as additional inputs, providing extra contextual information within the LSTM layer. This aids the model in better comprehending time series data and enhancing predictive performance. Specifically, when longitude and latitude fall within the SAA range, the model focuses more on the predictive outcome for that period. Proton flux, serving as the primary feature input, furnishes the core data information to the model. It directly participates in predicting proton flux, thereby achieving more precise predictive outcomes.
In the TPA-LSTM neural network model we employed, there are two crucial hyperparameters: hidden state dimensions and epochs. Properly configuring the number of hidden neurons can not only enhance predictive accuracy but also effectively prevent overfitting. To determine the appropriate hidden state dimensions, during the experimentation process, we initially utilized an empirical formula to obtain a rough range. The empirical formula is as follows:
N h = N s / ( α ( N i + N o ) )
where N s is the number of samples in the training set, N i is the hidden state dimension of the input layer, N o is the hidden state dimension of the output layer, and α is an arbitrary value variable that can be taken by itself, usually 2–10.
For our neural network model, the training dataset size for each case was approximately 80,000 samples, while the hidden state dimensions for both the input and output layers were set to 16. Through calculations, we determined that the approximate range for the hidden state dimensions should be between 250 and 1200. During the tuning process, we experimented within this range to obtain the optimal hidden state dimensions.
The number of epochs, as one of the hyperparameters, also has a significant impact on the predictive results of the neural network. Insufficient epochs can lead to the underfitting of the model, while too many epochs can result in overfitting.
Figure 10 illustrates the impact of hidden state dimensions and epochs on the logRMSE value. With 40 epochs and a hidden state dimension of 800, the predicted logRMSE value reaches 3.00. This set of hyperparameters strikes a suitable balance between prediction accuracy and computational cost.

3.4. Experiments and Analysis

Figure 11 depicts a comparison of proton flux predictions among the TPA-LSTM model and the AP-8 model (model version: solar maximum, threshold flux for exposure: 1 cm2/s). The AP-8 model is primarily used to describe electron and proton fluxes between low Earth orbit (LEO) and the Earth’s atmosphere. It provides important parameters for characterizing these particle flows, such as flux and energy spectrum distribution [25]. Figure 11a represents the prediction scenario for case#21, which includes 16 peaks, each of them caused by the satellite transiting the SAA region. Appendix A Table A2 provides detailed information on the satellite transits through the SAA region, including the entrance time, exit time, and duration of each crossing. Figure 11b depicts the prediction scenario for the third peak within Figure 11a.
To mitigate the influence of prediction randomness, we conducted a minimum of 50 runs for each model across various scenarios, and the average logRMSE and logMAE values are listed in Table 3. From the results listed in Table 3, the following conclusions can be drawn: All two models exhibit commendable performance, with the TPA-LSTM model demonstrating the highest proficiency, achieving a logRMSE of 3.71. Notably, the TPA-LSTM model outperforms the AP-8 model in terms of predictive accuracy, boasting an approximately 2.03 higher logRMSE compared to the AP-8 model.
When focusing on the South Atlantic Anomaly (SAA) region, the predictive performance of the AP-8 model experiences a decline, while the TPA-LSTM model’s predictive performance is enhanced. As shown in Table 4, the TPA-LSTM model’s logRMSE value decreases by 0.62, and the logMAE value decreases by 1.13 in the SAA region. The TPA-LSTM model demonstrates good predictive capability for proton flux in the SAA region, meeting the expected performance criteria. However, the AP-8 model exhibits poorer predictive performance in the SAA region. Its logRMSE value increases by 1.55, reaching 4.19, and its logMAE value increases by 2.04. This predictive performance falls short of our requirements for proton flux forecasting. The application of TPA distinguishes our model from the AP-8 model, and our model exhibits enhanced predictive performance in the South Atlantic Anomaly (SAA) region.
Figure 12 presents scatter plots of predicted values versus actual values for 26 cases for both the TPA-LSTM model and the AP-8 model. Figure 12a,c show scatter plots of predicted versus actual values for the entire time series for the TPA-LSTM model and AP-8 model, respectively, with the red line indicating a 5% error range. Figure 12b,d depict scatter plots of predicted versus actual values for the SAA region for the TPA-LSTM model and AP-8 model, respectively, with the red line indicating a 10% error range. By comparing the two models, we can clearly see that the data points from the TPA-LSTM model are more tightly distributed within the confidence region, both for the entire time series and within the SAA region. These two confidence regions encompass the majority of data points, showcasing how our predictive model accurately captures the variations in proton flux in the vast majority of cases. By examining Figure 12a, we can further observe that in regions with low proton flux intensity (typically, low-intensity proton flux is mostly found in non-SAA regions), the distribution of predicted points is more scattered. Conversely, in regions with high proton flux intensity, the distribution of predicted points is more clustered. This further validates our model’s grasp of spatial characteristics, leading to better predictive performance in the SAA region. Although a few data points lie outside the confidence region, overall, our predictive outcomes exhibit a high level of accuracy across the entire time series.
Additionally, in the error distribution plot within the SAA region, we can observe that the data points are more scattered compared to the entire time series. However, the majority of data points from the TPA-LSTM model still fall within the confidence region, indicating that our predictive model maintains a high level of accuracy even within this challenging zone. However, we acknowledge the greater challenge posed by predicting within the SAA region. Hence, in our future research, a heightened focus on optimizing predictions within the SAA region will be imperative to further enhance the model’s performance in this critical area.
In general, the results obtained through multiple runs and averaging indicate that the TPA-LSTM model demonstrates a high level of stability and accuracy in predicting proton flux. This model utilizes the TPA mechanism and LSTM architecture to capture the time- dependence of proton flux profiles measured during satellite operation more comprehensively and with heightened sensitivity, particularly within the SAA region. This sensitivity proves beneficial for enhancing the precision of proton flux predictions. However, it is important to note that forecasting within the SAA region remains a challenging task, necessitating further refinement and optimization to meet higher predictive demands.

4. Conclusions

This paper explores the use of a novel hybrid neural network model named TPA-LSTM for predicting proton flux profiles along a specific satellite orbit mentioned in the article. We applied the moving average wavelet transform method to preprocess the proton flux data. Subsequently, leveraging insights from the analysis of proton flux data, we introduced latitude and longitude features, along with the TPA attention mechanism, to enhance the focus on proton flux variations within the South Atlantic Anomaly (SAA) region. Finally, the TPA-LSTM hybrid neural network model was designed and trained and validated using proton flux data during the operational period of the specified satellite orbit. This model demonstrates higher accuracy in predicting proton flux compared to the AP-8 model. The key findings of this study are as follows:
  • The moving average wavelet transform method demonstrates effective data preprocessing performance, proficiently retaining trend information while also effectively denoising the original dataset.
  • The hybrid neural network model presented in this paper demonstrates a high level of predictive performance, with an average logRMSE of 3.4. Even in the case of the lowest predictive performance, the logRMSE still reaches 3.71.
  • The neural network model, enhanced with latitude and longitude features and the TPA mechanism, demonstrates improved accuracy in capturing crucial time intervals influenced by latitude and longitude features within the SAA region. This enhancement effectively boosts the predictive performance of proton flux, with an average logRMSE of 3.09 for the SAA region.
Hence, the TPA-LSTM neural network model demonstrates a high level of accuracy in the field of predicting radiation proton flux. The implementation of the TPA attention mechanism allows for the comprehensive utilization of spatiotemporal features. Our approach can provide valuable assistance in space science research, such as the prediction of solar proton events. However, it is important to note that predicting space weather events requires not only a focus on particle flux data but also on features such as solar activity, magnetic field status, and more, and this is one of our future research directions.

Author Contributions

Z.Z.: calculations and writing. Y.X.: conceptualization and methodology. X.Z.: conceptualization, methodology, writing—review and editing, and supervision. L.L.: Methodology. L.Q.: Formal analysis and Methodology. G.S.: Resources and Methodology. R.Z.: Resources. Y.J.: Methodology. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (61474096, 12004329), the Yangzhou Science and Technology Bureau (YZ2020263), the Open Project of State Key Laboratory of Intense Pulsed Radiation Simulation and Effect (SKLIPR2115), the and Foundation of National Key Laboratory of Materials Behavior and Evaluation Technology in Space Environment (WDZC-HGD-2022-11). The APC was funded by Yangzhou University.

Data Availability Statement

The datasets used and analyzed during the current study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare that they have no known competing financial interest or personal relationships that could have appeared to influence the work reported in this paper.

Appendix A

Table A1. Detailed information of time periods.
Table A1. Detailed information of time periods.
No.Time RangeDuration (min)
12022/04/17 22:37–2022/04/18 12:09671
22022/04/24 10:33–2022/04/25 11:071394
32022/05/07 09:30–2022/05/08 11:261405
42022/05/09 08:58–2022/05/10 10:581440
52022/08/11 08:28–2022/08/13 07:562867
62022/08/13 07:56–2022/08/15 07:242871
72022/08/15 07:24–2022/08/17 06:522877
82022/08/17 06:52–2022/08/19 06:202887
92022/08/19 06:20–2022/08/21 05:482891
102022/08/25 04:44–2022/08/27 04:112870
112022/08/27 04:11–2022/08/29 03:392870
122022/09/02 02:35–2022/09/04 02:032870
132022/09/04 02:03–2022/09/06 01:312877
142022/09/06 01:31–2022/09/08 00:592877
152022/09/17 01:33–2022/09/19 00:412877
162022/11/16 10:33–2022/11/18 10:002877
172022/11/18 10:00–2022/11/20 09:282668
182022/11/20 09:28–202211/22 08:562672
192022/11/22 08:56–2022/11/24 08:242672
202022/11/24 08:24–2022/11/26 07:522668
212022/12/02 23:51–2022/12/04 23:192668
222022/12/06 22:47–2022/12/08 22:152668
232022/12/18 19:34–2022/12/20 10:111757
242023/01/01 21:08–2023/01/03 20:362870
252023/02/16 12:46–2023/02/18 12:132877
262023/02/18 12:13–2023/02/20 06:222528
Table A2. SAA transit times of case#21.
Table A2. SAA transit times of case#21.
No.Entrance TimeExit TimeDuration (min)
12022/12/04 00:042022/12/04 00:2419.6
22022/12/04 01:492022/12/04 02:1121.6
32022/12/04 03:362022/12/04 03:5922.4
42022/12/04 05:242022/12/04 05:4622.1
52022/12/04 07:102022/12/04 07:3222.3
62022/12/04 07:552022/12/04 07:593.2
72022/12/04 08:542022/12/04 09:0813.1
82022/12/04 09:392022/12/04 09:5111.9
92022/12/04 10:392022/12/04 10:466.8
102022/12/04 11:202022/12/04 11:4019.4
112022/12/04 13:042022/12/04 13:2722.9
122022/12/04 14:502022/12/04 15:1322.1
132022/12/04 16:392022/12/04 16:5616.4
142022/12/04 18:312022/12/04 18:408.9
152022/12/04 20:212022/12/04 20:253.4
162022/12/04 21:132022/12/04 21:2613.0
172022/12/04 22:562022/12/04 23:1518.4

References

  1. Ya, N.; Zainudin, A.; Magdugal, R.; Naim, N.F. Mitigation of space radiation effects on satellites at Low Earth Orbit (LEO). In Proceedings of the 2016 6th IEEE International Conference on Control System, Computing and Engineering (ICCSCE), Penang, Malaysia, 25–27 November 2016; IEEE: Piscataway, NJ, USA, 2016; pp. 56–61. [Google Scholar]
  2. Cochran, D.J.; Chen, D.; Oldham, T.R.; Sanders, A.B.; Kim, H.S.; Campola, M.J.; Buchner, S.P.; LaBel, K.A.; Marshall, C.J.; Pellish, J.A.; et al. Total ionizing dose and displacement damage compendium of candidate spacecraft electronics for NASA. In Proceedings of the 2010 IEEE Radiation Effects Data Workshop, Quebec, QC, Canada, 20–24 July 2009; IEEE: Piscataway, NJ, USA, 2010; p. 8. [Google Scholar]
  3. Uzel, R.; Özyildirim, A. A study on the local shielding protection of electronic components in space radiation environment. In Proceedings of the 2017 8th International Conference on Recent Advances in Space Technologies (RAST), Istanbul, Turkey, 19–22 June 2017; IEEE: Piscataway, NJ, USA, 2017; pp. 295–299. [Google Scholar]
  4. Katz, S.; Goldvais, U.; Price, C. The connection between space weather and Single Event Upsets in polar low earth orbit satellites. Adv. Space Res. 2021, 67, 3237–3249. [Google Scholar] [CrossRef]
  5. Chen, J.; Lange, T.; Andjelkovic, M.; Simevski, A.; Lu, L.; Krstic, M. Solar Particle Event and Single Event Upset Prediction from SRAM-based Monitor and Supervised Machine Learning. IEEE Trans. Emerg. Top. Comput. 2022, 10, 564–580. [Google Scholar] [CrossRef]
  6. Nagatsuma, T.; Sakaguchi, K.; Kubo, Y.; Belgraver, P.; Chastellain, F.; Muff, R.; Otomo, T. Space environment data acquisition monitor onboard Himawari-8 for space environment monitoring on the Japanese meridian of geostationary orbit. Earth Planets Space 2017, 69, 75. [Google Scholar] [CrossRef]
  7. Turhan, B. On the dataset shift problem in software engineering prediction models. Empir. Softw. Eng. 2012, 17, 62–74. [Google Scholar] [CrossRef]
  8. Sajid, M.; Chechenin, N.G.; Torres, F.S.; Khan, E.U.; Agha, S. Space radiation environment prediction for VLSI microelectronics devices onboard a LEO satellite using OMERE-TRAD software. Adv. Space Res. 2015, 56, 314–324. [Google Scholar] [CrossRef]
  9. Bamisile, O.; Oluwasanmi, A.; Ejiyi, C.; Yimen, N.; Obiora, S.; Huang, Q. Comparison of machine learning and deep learning algorithms for hourly global/diffuse solar radiation predictions. Int. J. Energy Res. 2022, 46, 10052–10073. [Google Scholar] [CrossRef]
  10. Raju, H.; Das, S. CNN-based deep learning model for solar wind forecasting. Sol. Phys. 2021, 296, 134. [Google Scholar] [CrossRef]
  11. Li, H.; Gou, L.; Li, H.; Liu, Z. Physics-Guided Neural Network Model for Aeroengine Control System Sensor Fault Diagnosis under Dynamic Conditions. Aerospace 2023, 10, 644. [Google Scholar] [CrossRef]
  12. David, M.; Ramahatana, F.; Trombe, P.J.; Lauret, P. Probabilistic forecasting of the solar irradiance with recursive ARMA and GARCH models. Solar Energy 2016, 133, 55–72. [Google Scholar] [CrossRef]
  13. Mirzaei, M.; Yu, H.; Dehghani, A.; Galavi, H.; Shokri, V.; Mohsenzadeh Karimi, S.; Sookhak, M. A novel stacked long short-term memory approach of deep learning for streamflow simulation. Sustainability 2021, 13, 13384. [Google Scholar] [CrossRef]
  14. Dey, S.; Fuentes, O. Predicting solar X-ray flux using deep learning techniques. In Proceedings of the 2020 International Joint Conference on Neural Networks (IJCNN), Glasgow, UK, 19–24 July 2020; IEEE: Piscataway, NJ, USA, 2020; pp. 1–7. [Google Scholar]
  15. Wei, X.; Li, Y.; Shang, R.; Ruan, C.; Xing, J. Airport Cluster Delay Prediction Based on TS-BiLSTM-Attention. Aerospace 2023, 10, 580. [Google Scholar] [CrossRef]
  16. Yildirim, A.; Bilgili, M.; Ozbek, A. One-hour-ahead solar radiation forecasting by MLP, LSTM, and ANFIS approaches. Meteorol. Atmos. Phys. 2023, 135, 10. [Google Scholar] [CrossRef]
  17. Wei, L.; Zhong, Q.; Lin, R.; Wang, J.; Liu, S.; Cao, Y. Quantitative prediction of high-energy electron integral flux at geostationary orbit based on deep learning. Space Weather 2018, 16, 903–916. [Google Scholar] [CrossRef]
  18. Niu, Z.; Zhong, G.; Yu, H. A review on the attention mechanism of deep learning. Neurocomputing 2021, 452, 48–62. [Google Scholar] [CrossRef]
  19. Zhu, T.; Li, Y.; Li, Z.; Guo, Y.; Ni, C. Inter-hour forecast of solar radiation based on long short-term memory with attention mechanism and genetic algorithm. Energies 2022, 15, 1062. [Google Scholar] [CrossRef]
  20. Kong, X.; Du, X.; Xu, Z.; Xue, G. Predicting solar radiation for space heating with thermal storage system based on temporal convolutional network-attention model. Appl. Therm. Eng. 2023, 219, 119574. [Google Scholar] [CrossRef]
  21. Shih, S.Y.; Sun, F.K.; Lee, H. Temporal pattern attention for multivariate time series forecasting. Mach. Learn. 2019, 108, 1421–1441. [Google Scholar] [CrossRef]
  22. Belkadhi, K.; Manai, K. Dose calculation using a numerical method based on Haar wavelets integration. Nucl. Instrum. Methods Phys. Res. Sect. A Accel. Spectrometers Detect. Assoc. Equip. 2016, 812, 73–80. [Google Scholar] [CrossRef]
  23. Singla, P.; Duhan, M.; Saroha, S. A hybrid solar irradiance forecasting using full wavelet packet decomposition and bi-directional long short-term memory (BiLSTM). Arab. J. Sci. Eng. 2022, 47, 14185–14211. [Google Scholar] [CrossRef]
  24. Sharifi, S.S.; Rezaverdinejad, V.; Nourani, V.; Behmanesh, J. Multi-time-step ahead daily global solar radiation forecasting: Performance evaluation of wavelet-based artificial neural network model. Meteorol. Atmos. Phys. 2022, 134, 50. [Google Scholar] [CrossRef]
  25. Jordan, C.E. NASA Radiation Belt Models AP-8 and AE-8; Geophysics Laboratory, Air Force Systems Command, US Air Force: Washington, DC, USA, 1989; p. 0026. [Google Scholar]
Figure 1. Overall methodology.
Figure 1. Overall methodology.
Aerospace 10 00982 g001
Figure 2. Moving average wavelet transform processing flow.
Figure 2. Moving average wavelet transform processing flow.
Aerospace 10 00982 g002
Figure 3. LSTM neural network cell structure diagram.
Figure 3. LSTM neural network cell structure diagram.
Aerospace 10 00982 g003
Figure 4. The model structure of Temporal Pattern Attention (TPA).
Figure 4. The model structure of Temporal Pattern Attention (TPA).
Aerospace 10 00982 g004
Figure 5. Results of moving average method (m = 2, 3, 5).
Figure 5. Results of moving average method (m = 2, 3, 5).
Aerospace 10 00982 g005
Figure 6. The results of processing different wavelet bases. ((a) is the curve of the actual data, (b) is the curve of the data after processing with Haar wavelet basis, and (cf) are the curves of the data after processing with db2, db3, db4, db5 wavelet bases).
Figure 6. The results of processing different wavelet bases. ((a) is the curve of the actual data, (b) is the curve of the data after processing with Haar wavelet basis, and (cf) are the curves of the data after processing with db2, db3, db4, db5 wavelet bases).
Aerospace 10 00982 g006
Figure 7. Scatter distribution of high value in proton flux.
Figure 7. Scatter distribution of high value in proton flux.
Aerospace 10 00982 g007
Figure 8. (a) World map of the proton flux; (b) 3D-view of the proton flux.
Figure 8. (a) World map of the proton flux; (b) 3D-view of the proton flux.
Aerospace 10 00982 g008
Figure 9. The structure of TPA-LSTM.
Figure 9. The structure of TPA-LSTM.
Aerospace 10 00982 g009
Figure 10. Hyperparameter optimization of neural network model (the effects of hidden state dimensions and epochs on logRMSE).
Figure 10. Hyperparameter optimization of neural network model (the effects of hidden state dimensions and epochs on logRMSE).
Aerospace 10 00982 g010
Figure 11. (a) Comparison of prediction results from different models (including 16 peaks, all caused by satellite crossings through the SAA region). (b) Comparison of prediction results for the third peak.
Figure 11. (a) Comparison of prediction results from different models (including 16 peaks, all caused by satellite crossings through the SAA region). (b) Comparison of prediction results for the third peak.
Aerospace 10 00982 g011aAerospace 10 00982 g011b
Figure 12. The scatter plots between the predicted values and actual values of the (a,b) TPA-LSTM model and (c,d) AP-8 model (the entire time series and the SAA region).
Figure 12. The scatter plots between the predicted values and actual values of the (a,b) TPA-LSTM model and (c,d) AP-8 model (the entire time series and the SAA region).
Aerospace 10 00982 g012
Table 1. Six orbital elements of the satellite 52085.
Table 1. Six orbital elements of the satellite 52085.
ElementValue
Semi-major axis7395 km
Eccentricity0.01°
Inclination63.398°
Perigee argument41.145°
RAAN132.892°
Mean anomaly319.739°
Orbital period105.47 min
Table 2. Evaluation criteria of different wavelet bases after wavelet transform.
Table 2. Evaluation criteria of different wavelet bases after wavelet transform.
Metrics/Basisdb2db3db4db5Haar
logRMSE3.263.573.743.573.80
logMAE1.632.092.091.642.10
Table 3. Evaluation criteria of prediction results of different models.
Table 3. Evaluation criteria of prediction results of different models.
Metrics/ModelTPA-LSTMAP-8
logRMSE3.714.19
logMAE2.253.26
Table 4. Evaluation criteria of prediction results of different models in SAA.
Table 4. Evaluation criteria of prediction results of different models in SAA.
Metrics/ModelTPA-LSTMAP-8
logRMSE3.095.74
logMAE1.125.30
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhang, Z.; Liu, L.; Quan, L.; Shen, G.; Zhang, R.; Jiang, Y.; Xue, Y.; Zeng, X. A Proton Flux Prediction Method Based on an Attention Mechanism and Long Short-Term Memory Network. Aerospace 2023, 10, 982. https://doi.org/10.3390/aerospace10120982

AMA Style

Zhang Z, Liu L, Quan L, Shen G, Zhang R, Jiang Y, Xue Y, Zeng X. A Proton Flux Prediction Method Based on an Attention Mechanism and Long Short-Term Memory Network. Aerospace. 2023; 10(12):982. https://doi.org/10.3390/aerospace10120982

Chicago/Turabian Style

Zhang, Zhiqian, Lei Liu, Lin Quan, Guohong Shen, Rui Zhang, Yuqi Jiang, Yuxiong Xue, and Xianghua Zeng. 2023. "A Proton Flux Prediction Method Based on an Attention Mechanism and Long Short-Term Memory Network" Aerospace 10, no. 12: 982. https://doi.org/10.3390/aerospace10120982

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

Article Metrics

Back to TopTop