1. Introduction
With the development of space technology, satellites flying in space are becoming increasingly complex. In the absence of atmospheric protection, the complex space weather environment is very prone to affecting the satellite payload, endangering satellite safety and normal operation. There are many international and domestic satellites carrying space weather monitoring payloads [
1,
2,
3,
4]; among these, the NASA Global-scale Observations of the Limb and Disk (GOLD) mission has flown an ultraviolet-imaging spectrograph on SES-14, which uses two identical channels to observe Earth’s far-ultraviolet (FUV) airglow at a distance of ~134–162 nm [
5]. The Swarm mission of the European Space Agency’s Earth Exploration Program (ESA) uses three satellites in a lower altitude, two of which measure the gradient of the magnetic field. The measured magnetic field is the sum of many contributions including the magnetic fields and currents inside the Earth, as well as the currents in Earth’s space [
6]. The FY-4 series is China’s new-generation geosynchronous-orbit meteorological satellite. Currently in orbit, it includes Fengyun-4A (hereinafter referred to as FY-4A) and Fengyun-4B (hereinafter referred to as FY-4B). The space weather instruments carried by FY-4A include high-energy electron detectors, high-energy proton detectors, fluxgate magnetometers, and radiation dosimeters, respectively. The FY-4B space weather instrument is composed of a charging potential monitor, on the basis of inheriting the high-energy charged particles and space magnetic field detection capabilities of the FY-4A satellite space weather monitoring instrument package. It has the same ability to detect high-energy charged particles and space magnetic fields, and can detect the medium-energy charged particles, low-energy charged particles, multi-band ionospheric ultraviolet spectral image, solar polar ultraviolet image, and X-EUV flow. Double satellites conducting space exploration simultaneously in high orbit can greatly enhance the space weather operation capabilities, provide guarantees for space weather users such as the aerospace industry, adapt to the development of national space weather monitoring and warning products, enhance national space weather monitoring, warning, and service capabilities, provide space weather monitoring and warning products with independent monitoring data for the public and professional users, and improve China’s space weather monitoring, warning, and service capabilities [
7,
8]. The FY4-(02) batch ground application system serves as the full business process support system for the FY-4B satellite, which is divided into ten major systems according to function, including the Space Weather Application System (SWAS) and Mission Management & Control System (MCS) [
9,
10]. Among them, the SWAS uploads the status information of space weather products and subsystems at all levels to the MCS, and the MCS s is responsible for the integrated safety management and scheduling control of satellite and ground application systems [
11,
12,
13].
With the increasing number of Fengyun-series meteorological satellites in orbit, the probability of various types of payload devices being affected by space weather continues to increase. The high timeliness and reliability requirements for space weather forecasting and the early warning for satellite safety pose a challenge to the existing satellite safety and health management system. Many studies have been conducted on the impact of space weather factors on satellite safety and protective measures. Robin and John found that surface charging and deep dielectric charging are two common causes of satellite anomalies [
14]; Randall et al. found that high-energy particle precipitation during solar storms leads to the production of a large amount of nitrogen oxides in the upper atmosphere [
15]; Rozanov et al. found that the response amplitude of the atmosphere to high-energy electron precipitation events may exceed the influence of solar ultraviolet flux [
16]; Fennell et al. found that, when the differential potential value is large and negative, it usually occurs at night, and the number of occurrences during the day is less than 1% [
17]. The high-energy protons from progressive solar high-energy particle events may pose significant radiation hazards to astronauts, space equipment, and high-altitude aircraft passengers flying polar routes [
18]. W. DRGE et al. analyzed the 65–105 keV electrons observed in a solar electron event observed by STEREO-A, STEREO-B, and ACE, and proposed, based on a three-dimensional particle propagation model, that, in the energy range of approximately 60–300 keV, due to the non-radial divergence of magnetic field lines or particle diffusion, the lateral transmission of electrons partially occurs near the sun and partially occurs in the interstellar medium [
19]. Verkhogliadova et al. used the Monte Carlo method to model the transport of particles escaping from impact, obtained the time intensity distribution of protons and iron ions, and analyzed the radial correlation of the peak flux (J) of protons and iron ions from 0.5 to 2AU [
20]. Wang et al. studied the occurrence time of solar high-energy particle events through numerical simulation, using a three-dimensional focused transport model to calculate the proton flux observed in the ecliptic at 1 AU in the energy range between 10 MeV and 80 MeV, confirming that the VDA method is effective for the pulse source duration, low background, weak scattering in interplanetary space, or rapid diffusion in the solar atmosphere [
21].
China is still in the exploratory stage of predicting and warning us about the space environment index for satellite safety [
22,
23]. During the construction period of the MCS, in order to meet the risk prediction and real-time warning requirements for the impact of space weather factors on satellite safety [
24,
25], based on the SWAS subsystem space weather data, the space weather payload detection results of the FY-4A and FY-4B are used for the first time and the National Center for Space Weather (NCSW) prediction results, combined with some open-source data, are designed to predict the comprehensive risk index of multi-star high-energy protons, low-energy particles, and high-energy electrons, as well as real-time risk index warnings for multi-satellite SEU, surface charging, and deep charging. Furthermore, through the deployment of a full-process multi-satellite space environment risk prediction and real-time warning automation system in the MCS subsystem, a commercial 24 h prediction and real-time warning service for the safety of Fengyun-series high-orbit satellites is achieved for the first time, improving the timeliness, stability, and reliability of satellite safety space environment assessment automation.
2. Data and Methods
The multi-satellite space environmental risk prediction and real-time warning system mainly uses independent data sources from FY-4A and FY-4B satellites, provided by the National Space Weather Monitoring and Warning Center [
26]. Through the Fengyun-4 ground application system business operation platform, it can continuously and stably obtain all the data required for forecasting and warning level division functions. Part of the open-source data used is provided by the Space Weather Prediction Center (SWPC) and the World Data Center for Geomagnetism (WDC Kyoto) at Kyoto University in Japan, which can effectively improve the time accuracy of risk forecasting.
The Fengyun-4 meteorological satellite space weather business system constructed by NCSW is based on real-time and continuous monitoring of near-Earth space environment data by FY-4A and FY-4B satellites, achieving the full business process of space weather monitoring, data preprocessing, product generation, storage, and distribution, and, through a complete service system that includes service product generation, product application, and evaluation, it meets the different needs of decision-making for the public and the professional users [
27]. Space Weather Prediction Center (SWPC) is a laboratory and weather forecasting center of the National Oceanic and Atmospheric Administration (NOAA) of the United States [
28]; its website (
https://www.swpc.noaa.gov/, accessed on 7 May 2023) updates 1-day data every 5 min, including electron flux, geomagnetic index, proton flux X-ray flux, planetary K-index, real-time solar wind, and K-station and A-station indices. At the same time, the website can provide a 27-day outlook for 3-day geomagnetic forecasts, sunspot numbers, solar cycles, 10.7 cm radio flux, and geomagnetic indices, with proton event probability data and the high-energy proton Fp index, the solar wind, and Ap index data. World Data Center for Geomagnetism, Kyoto, Japan (
https://wdc.kugi.kyoto-u.ac.jp, accessed on 8 May 2023) provides geomagnetic field data on its official website, including AE index, ASY/SYM index, and Dst index [
29]. The introduction of automatic data acquisition by the multi-satellite space environmental risk prediction and real-time warning system is shown in
Table 1.
The design of a multi-satellite space environment risk prediction and real-time warning system focuses on various factors that affect satellite safety and space weather. The comprehensive risk prediction includes three types of risk prediction: high-energy protons, low-energy particles, and high-energy electrons. The comprehensive risk real-time warning includes three types of risk warnings: SEU, surface charging, and deep charging, all of which are expressed by risk indices.
2.1. High-Energy Proton Risk Prediction
The prediction probability of solar proton events is positively correlated with the flux of solar high-energy protons [
30,
31] (reference). This article uses the high-energy proton risk prediction at 02:00 UT every day based on the probability of solar proton events P1 published by NCSW on T + 0 day and the probability prediction of high-energy proton Fp value P2 based on SWPC, respectively, selecting the larger
p value for level classification to predict the hourly high-energy proton risk level for FY-4A and FY-4B satellites in the next 24 h [
32,
33]. The risk level is set using the Chinese national industry standard solar proton event intensity classification [
34] and combined with the author’s satellite business operation experience, as shown in
Table 2.
Based on engineering experience [
35], when SWPC shows a high-energy proton flux Fp with energy greater than 10 MeV at 22:00 on T − 1, which is higher than 10 cm
−2s
−1sr
−1, the high-energy proton risk of the T + 0 satellite is directly labeled as Level 3.
2.2. Low-Energy Particle Risk Prediction
High-flux low-energy particles are the source of surface charging for satellites. This article is based on the geomagnetic AE index for hourly calculation, and the Ap index prediction results for daily calculation. By predicting the relevant geomagnetic disturbance indices (AE and Ap), the higher value of the two risk levels is adopted, achieving the prediction of low-energy particle risk [
36,
37,
38,
39,
40,
41].
The increase in geomagnetic AE index is positively correlated with the corresponding increase in local plasma, and is related to the local time of the satellite. This article adopts the author’s early warning method for the surface charging effect of geosynchronous orbit satellites and the method patented by the early warning system [
42]. The local time of the satellite’s location is divided into four time domains, and the corresponding relationship between the AE index in each time domain and the surface charging potential is different. Among them, the safe period is from 6:00 to 21:00 local time at the Subsatellite Point (SSP), during which there will be no surface charging effect; i.e., the surface potential will always be 0.
SSP 21:00 LT–24:00 LT is transition zone I. Taking the AE time point as the satellite local time (TL), the calculated values of the surface charging potential model for the time period from TL − 1 to TL + 1 (time units in hours) are as follows:
In Equation (1), AE is the latest hourly average of the AE index, in nT; U1 calculates the surface charging potential in volts for the transition zone I mode.
SSP 0:00 LT–4:00 LT is the disturbance zone. Taking the AE time point as the satellite local time (TL), the calculated values of the surface charging potential model for the time period from TL − 1 to TL + 1 (time units in hours) are as follows:
In Equation (2), AE is the latest hourly average of the AE index, in nT; U2 calculates the surface charging potential in volts for the disturbance zone mode.
SSP 4:00 LT–6:00 LT is the transition zone II. Taking the AE time point as the satellite local time (TL), the calculated values of the surface charging potential model for the time period from TL − 1 to TL + 1 (time unit is hour) are as follows:
In Equation (3), AE is the latest 1 h average of the AE index, in nT; U3 calculates the surface charging potential in volts for the transition zone II mode.
Hourly risk prediction of low-energy particles is based on the surface charging potential calculation model provided by the surface charging detection results of Fengyun-4, where U = U
1 ∪ U
2 ∪ U
3, represents the time zone where different real-time satellite points are located. To adapt to the application of other satellites in geosynchronous orbit, the surface charging potential model calculation value U of AE index is used to classify the level of surface charging danger. The level setting is based on the author’s long-term operational experience in satellite management, as expressed in
Table 3.
The daily prediction adopts the FY-4A and FY-4B satellite low-energy particle risk prediction daily products, with daily 01:00 UT analysis. The NCSW product Ap index prediction data are read for 24 h prediction value for risk assessment. The setting of surface charging hazard level is based on the author’s engineering construction experience, as shown in
Table 4.
2.3. High-Energy Electron Risk Prediction
High-energy electron risk prediction obtains high-energy particle data, and performs prediction calculations and deep charge amplitude prediction calculations separately to determine the level of high-energy deep charge [
43,
44,
45,
46,
47]. This article adopts the author’s patented method for predicting high-energy electron burst events [
48] to achieve high-energy particle deep charge risk prediction.
The 2 MeV Daily high-energy particle risk prediction T + 0 daily E > 2 MeV high-energy electron daily flux is calculated using the following formula:
In Equation (4), V is the average speed of the solar wind on T − 1 day and is the high-energy electron flux on T = 1 day, and Ap is the SWPC on T − 1 day Ap index.
The daily high-energy particle risk prediction obtains the high-energy electron differential flux counting products of the FY-4A and FY-4B satellites on the T + 0 day, analyzes the daily integrated flux data of E > 2 MeV high-energy electron flux Ft at a sampling interval of 2 s on the T − 1 day, and accumulates it to double the calculation.
The high-energy electron risk forecast is launched once a day at 03:00 UT, and the deep charge level of high-energy particles can be determined by the high-energy electron flux of E > 2 MeV on that day. The risk level is determined using the Chinese meteorological industry standard electronic daily integral intensity grading with energy above 2 MeV at the geostationary orbit [
49], and is set based on the author’s engineering construction experience [
50], as shown in
Table 5.
2.4. Real-Time Warning of Single-Event Upset Risk
The probability of single-event upset significantly increases in the process of higher strong proton events, mainly caused by high-energy protons with medium to high energy. The main characteristic parameter is the high-energy proton flux of 80–165 MeV [
51,
52]. Therefore, after obtaining the FY-4A and FY-4B satellite high-energy proton flux FP prediction data products, proton energy spectrum calculations are carried out. The real-time warning subsystem for single-event upset risk automatically analyzes high-energy proton flux data of 80–165 MeV, averages FP data over a 5 min period, calculates the probability of single-event upset, and statistically analyzes the corresponding relationship between high-energy proton flux and single particle event occurrence rate. The risk level setting for single-event upset is based on the author’s engineering construction experience [
53], as shown in
Table 6.
2.5. Real-Time Warning of Surface Charging Risk
The subsatellite point of FY-4A and FY-4B satellites are different, and the risk of low-energy particles is different [
54]. The real-time warning subsystem for surface charging risk in multi-satellite space environment obtains data from FY-4A and FY-4B satellite surface differential potential upward detection, automatically analyzes the differential potential value U of satellite upward surface potential probe C with a sampling interval of 2 s within 1 min at the latest time [
55], judges the extra satellite and other ionospheric environments, and divides the levels of charging events for each satellite. Based on long-term operational experience in satellite management, the author of this article sets the surface charging warning level recognition level, as shown in
Table 7.
2.6. Real-Time Warning of Deep Charging Risk
The risk prediction of deep charging in multi-satellite space environment is initiated every 5 min to obtain the latest high-energy electron differential flux counts of FY-4A and FY-4B satellites, analyze the E > 2 MeV high-energy electron flux Ft data at a sampling interval of 2 s within 5 min [
56], average the flux during this period, and obtain the net input charge of multi-satellite deep charging. By analyzing the corresponding relationship between the input charge and deep charging voltage, the deep charging warning level is determined, realizing real-time warning of deep charging [
57,
58]. The setting of the multi-star deep charging warning level is based on the author’s engineering construction experience, as shown in
Table 8.
2.7. Comprehensive Risk Index Strategy
The multi-satellite space environment risk prediction and warning system needs to visually display the impact of the space environment on satellites. For various space weather factors, a comprehensive risk index is used to express the risk. Among them, high-energy protons, low-energy particles, and high-energy electrons are the three risk predictions, and single-event upset, surface charging, and deep charging are the three risk warnings. The highest one is selected as the comprehensive risk prediction level and the comprehensive risk real-time warning level. For the convenience of business personnel, the comprehensive risk index consists of risk levels and textual prompts. The first-level risk prompt is for safety, the second-level risk prompt is for attention, and the third-level risk prompt is for caution. The process of spatial environmental comprehensive risk index is shown in
Figure 1. The comprehensive risk index guideline is shown in
Table 9 and
Table 10.
4. Conclusions
The multi-satellite space environmental risk prediction and early warning system is the first to address the safety of Fengyun-series high-orbit satellites in the face of complex spatial environmental factors. It is based on the FY-4A and FY-4B space weather detection results, NCSW prediction results, and some open-source data for the first time. It adopts Chinese national industry standards, Chinese meteorological industry standards, and invention patents, and combines the author’s years of engineering practice and satellite safety management experience to focus on the high-energy protons of FY-4A and FY-4B satellites, low-energy particle and high-energy proton risk prediction, the real-time warning for a single-event upset, surface charging, and deep charging to design and deploy an automatic system. At present, the main construction of the system has been completed and business trial operation services have been provided to users. For the first time, an intelligent spatial risk assessment has been implemented for the safety of multiple stationary meteorological satellites, generating risk index prediction products hour by hour and risk index warning products minute by minute. The error between the predicted results of various risk indices and the measured data is less than one level, with a prediction accuracy that is better than 80% and a warning accuracy that is better than 90%. During the typical space weather event on 21 April 2023, the FY-4A and FY-4B low-energy particles were correctly predicted and warned about, as well as the high-risk level of surface charging. The verification results showed that the system design theory was correct, the system operation was stable, and the prediction and real-time warning results were effective and reliable. There are many factors that affect satellite safety in complex space weather environments, with a wide variety of satellite platforms and payload types. The successful deployment of this system can play an important data support role in the pre-deployment of satellite security management during the active period of the future space environment, in response to the requirements of satellite security management.
Space weather prediction and the early warning for the safety of spacecraft are a fundamental capability directly related to the survival and development of human beings in the space age. This paper only focuses on the prediction and early warning of various elements of the space environment for the safety of high-orbit meteorological satellites, and preliminarily establishes a set of operational systems. In the future, the system can rely on the integrated monitoring data of the “causal chain” in solar–terrestrial space and the panoramic view of the earth space to build an intelligent risk assessment system for regional spatial environmental disturbances and the effects of high- and low-orbit satellites, aircraft, navigation, communication, radar, power, and other facilities. With the continuous deepening of research on the basic physical processes that affect space weather changes, the construction of intelligent prediction and early warning systems with machine deep-learning capabilities, the development of a “full chain” system for space weather monitoring, research, modeling, and forecasting, and the establishment of globalized and standardized services for different spatial regions. It is an important trend for the future development of space weather prediction and early warning, which can provide useful references for the development of society in the space environment.