Space Weather

On 3 February 2022, 38 satellites launched by SpaceX re-entered the atmosphere and were subsequently destroyed. An investigation found that a minor geomagnetic storm occurred on 3–4 February 2022 led to a neutral density enhancement and large atmospheric drag. To better understand the responses of the thermosphere to geomagnetic storms, the method proposed by Li et al. (2023, https://doi.org/10.1029/2022ja030988) was employed to extract exospheric temperature (Tex) from ionosonde electron density profiles (∼150–200 km) in Beijing (geolocation: 39.56°N; 116.2°E; geomagnetic location: 30.16°N; 172.08°W) station. The retrieved Tex was plugged into the NRLMSISE-00 model to calculate the corresponding neutral density. Derived results showed a ∼2%–7% enhancement in Tex and a ∼15%–38% enhancement in neutral density at 430 km. The relative deviation in neutral density on the satellites’ orbital trajectory ranges from ∼10% (210 km) to ∼35% (500 km) on 3 February, and from ∼13% (210 km) to ∼60% (500 km) on 4 February. Furthermore, the neutral density reproduced the variations observed by the SWARM-C satellite fairly well both on quiet and disturbed days. These results suggest that even a minor geomagnetic storm can cause significant changes in neutral temperature and neutral density at middle latitudes. Additionally, the application of our inversion method, combined with the global, long-term and real-time ionospheric observations from ionosondes, provides an opportunity to improve the capability of thermosphere forecasting and nowcasting.

Surface charging is one of the most common causes of spacecraft anomalies. When and to what potential the spacecraft is charged are two important questions in space weather. Here, for a Chinese geosynchronous navigation satellite, we infer the extreme negative surface charging potentials from the ion differential fluxes measured by a low-energy ion spectrometer. Without the solar eclipse effect away from the midnight, the charging potentials are found to have a negative limit which is determined by the maximum SuperMAG electrojet index in the preceding 2 hr. Such an empirical relation can be reasonably explained by the dependence of 1–50 keV electron fluxes on substorm strength. Similar relations may also exist for other inner magnetospheric spacecraft in the non-eclipse region, which would be useful for spacecraft engineering and space weather alerts.

Space weather can have major impacts on electrical infrastructure. Multiple instances of transformer damage have been attributed to geomagnetic storms in recent decades, for example, the Hydro Quebec incident of 1989 and the November 2001 storm in New Zealand. While many studies exist on the impacts of geomagnetic storms on power transformers in New Zealand, no studies exist that employ Dissolved Gas Analysis (DGA) techniques to relate geomagnetic storms to transformer gassing. A relationship has been reported between geomagnetic activity and DGA for South Africa, while none was found in a recent study in Great Britain. This paper attempts to examine this research question by examining dissolved gas data across eight power transformers in different substations in New Zealand from 2016 to 2019. Case studies were conducted which analyzed the DGA readings of each transformer alongside horizontal magnetic field component rate of change measurements at Eyrewell across six geomagnetic storms. These case studies were then augmented with an analysis of the entire data set where magnetic field measurements were compared with individual gas rates to establish a correlation between gas production and geomagnetic activity. Analysis of the results of this study concluded that no link had been found between the production of combustible gasses in a transformer and geomagnetic activity during the observation period. However, we note our dissolved gas analysis was largely in a geomagnetically quieter period, which may limit our analysis. The production of combustible gasses is not correlated to geomagnetic storms for the time period and transformers analyzed.

No abstract is available for this article.

Global ionospheric total electron content (TEC) maps exhibit TEC intensifications and depletions of various sizes and shapes. Characterizing key features on TEC maps and understanding their dynamic coupling with external drivers can significantly benefit space weather forecasting. However, comprehensive analysis of ionospheric structuring over decades of TEC maps is currently lacking due to large data volume. We develop feature extraction software based on image processing techniques to extract TEC intensification regions, that is, contiguous regions with sufficiently elevated TEC values than surrounding areas, from global TEC maps. Applying the software to the Jet Propulsion Laboratory Global Ionospheric Map data, we generate a TEC intensification data set for years 2003–2022 and carry out a statistical study on the number and strength of TEC intensifications. We find that the majority of the TEC maps (about 86%) are characterized with one or two intensification(s), while the rest of the TEC maps have three or more intensifications. Both the number and strength of TEC intensifications exhibit semi-annual variation that peaks near equinoxes and dips near solstices, as well as an annual asymmetry with larger values around December solstice compared to June solstice. The number and strength of intensifications increase with enhanced solar extreme-violet irradiance. The strength of intensifications also increases with elevated geomagnetic activity, but the number of intensifications does not. In addition, the number of intensifications is not correlated with the strength of intensifications.

An intelligent high-definition and short-term prediction of ionograms with/without Spread-F for the observation at Hainan (19.5°N, 109.1°E, magnetic 11°N) is presented in this paper, which comprises a spatio-temporal ConvGRU network and a super-resolution EDSR network. Our prediction is based on spatio-temporal features in the ionogram graph only. There are 469,227 ionograms classified into 5 categories, that is, frequency/range/mix/strong range/no Spread F, over a solar cycle (14 years) labeled manually by the research group, and we process these ionograms into two data sets for training the two networks mentioned above. A series of comprehensive experiments have been designed and conducted to determine the optimal super-parameters. Our method inputs 8 consecutive authentic ionograms (lasting 2 hr) and generates the next 2 figures (next 30 min). Remarkably, all predicted figures achieve a high accuracy rate of over 94% in predicting the occurrence of Spread-F.

Low Earth Orbit satellites offer extensive data of the radiation belt region, but utilizing these observations is challenging due to potential contamination and difficulty of intercalibration with spacecraft measurements at Highly Elliptic Orbit that can observe all equatorial pitch-angles. This study introduces a new intercalibration method for satellite measurements of energetic electrons in the radiation belts using a Data assimilation (DA) approach. We demonstrate our technique by intercalibrating the electron flux measurements of the National Oceanic and Atmospheric Administration (NOAA) Polar-orbiting Operational Environmental Satellites (POES) NOAA-15,-16,-17,-18,-19, and MetOp-02 against Van Allen Probes observations from October 2012 to September 2013. We use a reanalysis of the radiation belts obtained by assimilating Van Allen Probes and Geostationary Operational Environmental Satellites observations into 3-D Versatile Electron Radiation Belt (VERB-3D) code simulations via a standard Kalman filter. We compare the reanalysis to the POES data set and estimate the flux ratios at each time, location, and energy. From these ratios, we derive energy and L* dependent recalibration coefficients. To validate our results, we analyze on-orbit conjunctions between POES and Van Allen Probes. The conjunction recalibration coefficients and the data-assimilative estimated coefficients show strong agreement, indicating that the differences between POES and Van Allen Probes observations remain within a factor of two. Additionally, the use of DA allows for improved statistics, as the possible comparisons are increased 10-fold. Data-assimilative intercalibration of satellite observations is an efficient approach that enables intercalibration of large data sets using short periods of data.

Predicting the impacts of coronal mass ejections (CMEs) is a major focus of current space weather forecasting efforts. Typically, CME properties are reconstructed from stereoscopic coronal images and then used to forward model a CME's interplanetary evolution. Knowing the uncertainty in the coronal reconstructions is then a critical factor in determining the uncertainty of any predictions. A growing number of catalogs of coronal CME reconstructions exist, but no extensive comparison between these catalogs has yet been performed. Here we develop a Living List of Attributes Measured in Any Coronal Reconstruction (LLAMACoRe), an online collection of individual catalogs, which we intend to continually update. In this first version, we use results from 24 different catalogs with 3D reconstructions using Solar Terrestrial Relations Observatory observations between 2007 and 2014. We have collated the individual catalogs, determining which reconstructions correspond to the same events. LLAMACoRe contains 2,954 reconstructions for 1,862 CMEs. Of these, 511 CMEs contain multiple reconstructions from different catalogs. Using the best-constrained values for each CME, we find that the combined catalog reproduces the generally known solar cycle trends. We determine the typical difference we would expect between two independent reconstructions of the same event and find values of 4.0° in the latitude, 8.0° in the longitude, 24.0° in the tilt, 9.3° in the angular width, 0.1 in the shape parameter κ, 115 km/s in the velocity, and 2.5 × 1015 g in the mass. These remain the most probable values over the solar cycle, though we find more extreme outliers in the deviation toward solar maximum.

In this work, a model for calculating the galactic cosmic rays (GCRs) radiation dose rate is developed. The model is based on a GCR modulation model, which is established by Shen and Qin, and the fluence-dose conversion coefficients (FDCCs) published by the International Commission on Radiological Protection (ICRP). With the model, the radiation absorbed dose rate of GCRs near the lunar surface over long time periods is calculated and compared with the observation data from the Cosmic Ray Telescope for the Effects of Radiation and the Lunar Lander Neutron and Dosimetry. First, the energy spectrum of GCRs at 1 AU in the ecliptic, where the lunar orbit is located, is computed using the GCR modulation model. Then, using the FDCCs of ICRP 123, the absorbed dose rates of 15 human organs/tissues at the lunar orbit position are calculated to represent the general absorbed dose rate of the body (in water). Furthermore, considering the albedo radiation (excluding neutrons) and using the water-silicon conversion coefficients, the total absorbed dose rates of GCRs near the lunar surface (in silicon) are calculated, it is shown that our modeling results generally agree with the observations from spacecraft. This work is useful for future manned space exploration to the Moon or other celestial bodies in the solar system.

Current magnetohydrodynamics (MHD) models largely rely on synoptic magnetograms, such as the ones produced by the Global Oscillation Network Group (GONG). Magnetograms are currently available mostly from the front side of the Sun, which significantly reduces the accuracy of MHD modeling. Extreme Ultraviolet (EUV) images can instead be obtained from other vantage points. To investigate the potential, we explore the possibility of using EUV information from the Atmospheric Imaging Assembly (AIA) to directly generate the input for the state-of-the-art 3D MHD model European Heliospheric FORecasting Information Asset (EUHFORIA). Toward this goal, we develop a method called Transfer-Solar-GAN which combines a conditional generative adversarial network with a transfer learning approach to overcome training data set limitations. The source domain data set is constructed from multiple pairs of the central portion of co-registered AIA and Helioseismic and Magnetic Imager (HMI) line of sight (LOS) full-disk images, while the target domain is constructed from pairs of portions of AIA and GONG sine-latitude synoptic maps that we call segments. We evaluate Transfer-Solar-GAN by comparing modeled and measured solar wind velocity and magnetic field density parameters at the L 1 Lagrange point and along the Parker Solar Probe (PSP) trajectory which were determined with EUHFORIA using both empirical GONG and artificial-intelligence (AI)-synthetic synoptic magnetograms as inputs. Our results demonstrate that the Transfer-Solar-GAN model can provide the necessary information to run solar physics models by EUV information. Our proposed model is trained with only 528 paired image segments and enforces a reliable data division strategy.

Signals recorded by two stations in the Brazilian region: [Fortaleza (3.74°S, 38.57°W) and Inconfidentes (22.31°S, 46.32°W)], receiving L1 transmissions from different geostationary satellites, were used to evaluate the amplitude scintillation index S 4 and several characteristics of scintillation events (continuous record with S 4 > 0.2) during nighttime hours (18:00 LT–02:00 LT) in the years 2014–2016. The effects from solar activity, season, and local time on the number of scintillation events per night, maximum scintillation, scintillation event duration, and spacing between consecutive events will be discussed. The results indicate that: (a) scintillation occurs from September to March in both links; (b) the most likely numbers of observed scintillation events per night were two or three, particularly during the first 2 years; (c) on average, the first scintillation event usually had larger maximum S 4 values when compared to those of the later ones along the night; (d) the first scintillation event had a longer mean duration than the succeeding ones in a given night; (e) the durations of scintillation events, regardless of their numbers per night and the location, decreased with local time; (f) the opposite dependence of spacings between consecutive events on local time was observed; (g) the cumulative distribution functions of the scintillation onset time indicated a strong dependence on the dip latitude of the station; and (h) early occurrences of onset times are directly related to the increased probability of the occurrence of multiple scintillation events.