Interplanetary Space Research Articles

Three-Dimensional Guidance Laws for Spacecraft Propelled by a SWIFT Propulsion System

This paper discusses the optimal control law, in a three-dimensional (3D) heliocentric orbit transfer, of a spacecraft whose primary propulsion system is a Solar Wind Ion Focusing Thruster (SWIFT). A SWIFT is an interesting concept of a propellantless thruster, proposed ten years ago by Gemmer and Mazzoleni, which deflects, collects, and accelerates the charged particles of solar wind to generate thrust in the interplanetary space. To this end, the SWIFT uses a large conical structure made of thin metallic wires, which is positively charged with the aid of an electron gun. In this sense, a SWIFT can be considered as a sort of evolution of the Janhunen’s E-Sail, which also uses a (nominally flat) mesh of electrically charged tethers to deflect the solar wind stream. In the recent literature, the optimal performance of a SWIFT-based vehicle has been studied by assuming a coplanar orbit transfer and a two-dimensional scenario. The mathematical model proposed in this paper extends that result by discussing the optimal guidance laws in the general context of a 3D heliocentric transfer. In this regard, a number of different forms of the spacecraft state vectors are considered. The validity of the obtained optimal control law is tested in a simplified Earth–Venus and Earth–Mars transfer by comparing the simulation results with the literature data in terms of minimum flight time.

Occurrence of sudden storm commencement in interplanetary space

Sudden storm commencement (SSC) events detected in interplanetary space from 1869 to 2023, are subjectedto wavelet and Lomb-Scargle periodogram analyses to determine the potential frequency of their recurrence. This period includes SSC events from the preceding fourteen solar cycles. Using Lomb-Scargle periodogram analysis, we have identified SSC occurrences independently during the odd and even solar activity cycles. We also used the wavelet analysis approach throughout the positive (A > 0) and negative (A < 0) polarity states of the heliospheric magnetic fields, in addition to solar activity cycles. This study reveals some notable short-term periodic variations in the SSC time series. The findings indicate that SSC occurrences have a conspicuous and important ∼ 44.0-day period. We have also observed solar and extended solar rotation periods in the time series, along with some intermittent variations (e.g., ∼22.0-, ∼19.0-days). The observed periodicities in the SSC time series are discussed and compared with the previously detected periodicities of solar wind and plasma, geomagnetic activity, and cosmic rays’ intensity indicators.

Numerical Simulation of Equal Ratio Relations for the Peak Intensities of >10 MeV Energetic Protons

Previous studies have highlighted the significance of perpendicular diffusion in the decay phase of particle intensities for >10 MeV energetic protons. Recently, an observational study has indicated that the peak intensity ratios across different energy channels (13−64 MeV protons) remain almost constant as the spacecraft location varies in many solar proton events. This interesting phenomenon is referred to as equal ratio relations. These findings suggest that perpendicular diffusion not only affects particle intensity during the decay phase but also throughout the rising phase. In this study, we perform numerical simulations of >10 MeV energetic proton events observed by STEREO A, STEREO B, and the Solar and Heliospheric Observatory. Our findings demonstrate that perpendicular diffusion strongly affects the entire time profile of particle intensity for spacecraft not magnetically connected to the source region. The numerical simulation results indicate that in order to reproduce observations, we need to include perpendicular diffusion near the source region and in interplanetary space. Perpendicular diffusion leads to nearly uniform peak ratios at different locations and contributes to the formation of the reservoir phenomenon during the decay phase. Consequently, these numerical results support the significant role of perpendicular diffusion in the formation of the longitudinal distribution of >10 MeV proton events.

Large-amplitude transverse MHD waves prevailing in the Hα chromosphere of a solar quiet region revealed by MiHI integrated field spectral observations

The investigation of plasma motions in the solar chromosphere is crucial for understanding the transport of mechanical energy from the interior of the Sun to the outer atmosphere and into interplanetary space. We report the finding of large-amplitude oscillatory transverse motions prevailing in the non-spicular Hα chromosphere of a small quiet region near the solar disk center. The observation was carried out on 2018 August 25 with the Microlensed Hyperspectral Imager (MiHI) installed as an extension to the spectrograph at the Swedish Solar Telescope (SST). MiHI produced high-resolution Stokes spectra of the Hα line over a two-dimensional array of points (sampled every 0.066″ on the image plane) every 1.33 s for about 17 min. We extracted the Doppler-shift-insensitive intensity data of the line core by applying a bisector fit to Stoke I line profiles. From our time–distance analysis of the intensity data, we find a variety of transverse motions with velocity amplitudes of up to 40 km s−1 in fan fibrils and tiny filaments. In particular, in the fan fibrils, large-amplitude transverse MHD waves were seen to occur with a mean velocity amplitude of 25 km s−1 and a mean period of 5.8 min, propagating at a speed of 40 km s−1. These waves are nonlinear and display group behavior. We estimate the wave energy flux in the upper chromosphere at 3 × 106 erg cm−2 s−1. Our results contribute to the advancement of our understanding of the properties of transverse MHD waves in the solar chromosphere.

Мониторинг околоземного космического пространства, магнитосферы и атмосферы Земли в периоды Форбуш-эффектов в конце августа 2005 г.

We present the results of near-Earth interplanetary space, magnetosphere, and atmosphere monitoring during large-scale solar wind disturbances at the end of August 2005. The monitoring was carried out using ground-level cosmic ray (CR) observations made at the worldwide network of neutron monitors as well as muon telescopes in Yakutsk and Novosibirsk. As a result of the analysis performed by different methods, we have obtained variation properties of CRs of different rigidities in Earth’s orbit, their pitch angle anisotropy, orientation and configuration of the interplanetary magnetic field, changes in the planetary system of geomagnetic cutoff rigidities during geomagnetic disturbances, as well as mass average air temperature over CR stations equipped with muon telescopes. For the periods of geomagnetic disturbances, we have determined parameters of magnetospheric ring current and magnetopause currents.

Monitoring of near-Earth space, Earth’s magnetosphere and atmosphere during Forbush decreases in August 2005

We present the results of near-Earth interplanetary space, magnetosphere, and atmosphere monitoring during large-scale solar wind disturbances at the end of August 2005. The monitoring was carried out using ground-level cosmic ray (CR) observations made at the worldwide network of neutron monitors as well as muon telescopes in Yakutsk and Novosibirsk. As a result of the analysis performed by different methods, we have obtained variation properties of CRs of different rigidities in Earth’s orbit, their pitch angle anisotropy, orientation and configuration of the interplanetary magnetic field, changes in the planetary system of geomagnetic cutoff rigidities during geomagnetic disturbances, as well as mass average air temperature over CR stations equipped with muon telescopes. For the periods of geomagnetic disturbances, we have determined parameters of magnetospheric ring current and magnetopause currents.

How Switchbacks Can Maintain a Longer Time in the Interplanetary Space

Parker Solar Probe, the closest spacecraft to the Sun, has renewed our understanding of the solar corona and the interplanetary space. One of its important findings is the prevalence of switchbacks, which display localized magnetic reversals along the otherwise Parker spirals. While some switchbacks might disappear quickly, others can maintain for a long period of time, and there are indications that many switchbacks strengthen from the solar corona to the interplanetary space despite their magnetic tension force, which tends to straighten the magnetic field lines. Therefore, how these switchbacks could be maintained for a long period of time remains a mystery. In this paper, we employed a 3D data-driven global full magnetohydrodynamics numerical model to explore the evolution of switchbacks formed in the dynamic corona. Our simulations indicate that two factors can affect the lifetime of a switchback. One factor is the combination of angle and leg length ensures that the switchback with greater curvature after reconnection can last longer, and the greater the angle, the more magnetic field lines that can be reconnected, and thus the longer the duration. We call this influencing factor flux tube shape factor. The other factor is the velocity shear, i.e., when the solar wind at the convex-outward turning of a switchback is faster than that at the concave-outward turning, the switchback becomes enhanced during propagation, and it weakens when the velocity difference is opposite.

Multi-spacecraft study with the Icarus model. Modelling the propagation of CMEs to Mercury and Earth

Coronal mass ejections (CMEs) are the main drivers of the disturbances in interplanetary space. Earth-directed CMEs can be dangerous, and understanding the CME interior magnetic structure is crucial for advancing space weather studies. It is important to assess the capabilities of a numerical heliospheric model, as a firm understanding of the nature and extent of its limitations can be used to improve the model and the space weather predictions based on it. The aim of the present study is to test the capabilities of the recently developed heliospheric model Icarus and the linear force-free spheromak model that has been implemented in it. To validate the Icarus space weather modelling tool, two CME events were selected that were observed by two spacecraft located near Mercury and Earth, respectively. This enables us to test the heliospheric model computed with Icarus at two distant locations. The source regions for the CMEs were identified, and the CME parameters were determined and later optimised. Different adaptive mesh refinement levels were applied in the simulations to assess its performance by comparing the simulation results to in situ measurements. The first CME event erupted at 15:25 on July 9, 2013. The modelled time series were in good agreement with the observations both at MESSENGER and ACE. The second CME event started at 10:25 on February 16, 2014, and was more complicated, as three CME interactions occurred in this event. It was impossible to recover the observed profiles without modelling the other two CMEs that were observed, one before the main CME and one afterward. The parameters for the three CMEs were identified and the three CMEs were modelled in Icarus. For both CME studies, AMR level 3 was sufficient to reconstruct small-scale features near Mercury, while at Earth, AMR level 4 was necessary due to the radially stretched grid that was used. The profiles obtained at both spacecraft resemble the in situ measurements well. The current limitations of the space weather modelling tool result in an excessively small deceleration of the CME propagation during the CME--CME interaction as measured by MESSENGER and ACE.

Effects of Alpha Particle Fraction on Kinetic Plasma Turbulence

Solar activities have an extraordinary impact on interplanetary space, enriching the plasma dynamics including turbulent heating of various species. The small fraction of alpha particles is believed to play a significant role in the turbulent dynamics of the solar wind. Here we present fully kinetic particle-in-cell simulations to reveal the influences of the alpha particles in decaying plasma turbulence. Multiple run cases with different controlled variations of proton and alpha density are performed to compare and evaluate the energy conversion processes. It is found that the alpha particles can suppress the energy conversion rate with increasing density. Besides, the alpha particles show more heating intermittency than the proton species. Interestingly, on the other hand, the electrons do not show any change in their dynamics, including overall heating. These two positive charge species have more correlation in temperature anisotropy as their densities are comparable. Our results provide valuable insights on the turbulence with different species compositions to a certain extent, especially with abundant heavy particles.

Proton polytropic behavior of periodic density structures in the solar wind

Context. In recent years, mesoscales have gained scientific interest because they have been determined to be important in a broad range of phenomena throughout heliophysics. The solar wind mesoscale structures include periodic density structures (PDSs), which are quasi-periodic increases in the density of the solar wind that range from a few minutes to a few hours. These structures have been extensively observed in remote-sensing observations of the solar corona and in in situ observations out to 1 AU, where they manifest as radial length scales greater than or equal to the size of the Earth’s dayside magnetosphere, that is, from tens to hundreds of Earth radii (RE). While the precise mechanisms that form PDSs are still debated, recent studies confirmed that most PDSs are of solar origin and do not form through dynamics during their propagation in the interplanetary space. Aims. We further investigate the origin of PDSs by exploring the thermodynamic signature of these structures. To do this, we estimate the values of the effective polytropic index (Y) and the entropy of protons, which in turn are compared with the corresponding values found for the solar wind. Methods. We used an extensive list of PDS events spanning more than two solar cycles of Wind measurements (the entire Wind dataset from 1995 to 2022) to investigate the thermodynamic signatures of PDSs. With the use of wavelet methods, we classified these PDSs as coherent or incoherent, based on the shared periodic behavior between proton density and alpha-to-proton ratio, and we derive the proton polytropic index. Results. Our results indicate that the coherent PDSs exhibit lower Y values (Ῡ≈1.54) on average and a higher entropy than the values in the entire Wind dataset (Ῡ≈1.79), but also exhibit similarities with the magnetic cloud of an interplanetary coronal mass ejection. In contrast, incoherent PDSs exhibit the same Y values as those of the entire Wind dataset.

The winding number of coronal flux ropes

Context. Magnetic flux ropes are key structures in solar and solar-terrestrial research. Their magnetic twist is an important quantity for understanding their eruptivity, their evolution in interplanetary space, and their consequences for planetary space environments. The magnetic twist is expressed in terms of a winding number that describes how many times a field line winds about the axis of the flux rope (FR). Due to the complexity of calculating the winding number, current methods rely largely on its approximation. Aims. We use a data-driven simulated FR to investigate the winding number Tg in comparison to the commonly used twist proxy Tw, which describes a winding of two infinitesimally close field lines. We also estimate the magnetic flux enclosed in the resultant FR(s). Methods. We use the magnetic field analysis tools (MAFIAT) software to compute Tg and Tw for data-driven time-dependent magnetofrictional modelling of AR12473. Results. We find that the FR boundaries can significantly differ depending on whether they are defined using the twist approximation Tw or the winding number Tg. This also significantly affects the FR structure and the estimates of the enclosed magnetic flux. For the event analysed here, Tg also reveals that the twisted flux system consists of two separate intertwined FRs. Conclusions. The results of this study suggest that the computation of the winding number (Tg) is important for investigations of solar FRs.

Three-stage Acceleration of Solar Energetic Particles Detected by Parker Solar Probe

Coronal mass ejections (CMEs) drive powerful shocks and thereby accelerate solar energetic particles (SEPs) as they propagate from the corona into interplanetary space. Here we present the processes of three-stage particle acceleration by a CME-driven shock detected by the in situ spacecraft—Parker Solar Probe (PSP) on 2022 August 27. The onset of SEPs is produced by a fast CME with a speed of 1284 km s−1 when it propagates to ∼2.85 R ⊙. The second stage of particle acceleration occurs when the fast CME catches up and interacts with a preceding slow one in interplanetary space at ∼40 R ⊙ (∼0.19 au). The CME interaction is accompanied by an intense interplanetary type II radio enhancement. Such direct measurement of particle acceleration during interplanetary CME interaction/radio enhancement is rarely recorded in previous studies. The third stage of energetic storm particles is associated with the CME-driven shock passage of the PSP at ∼0.38 au. Obviously, harder particle spectra are found in the latter two stages than the first one, which can arise from a stronger shock produced by the CME interaction and the enriched seed particles inside the preceding CME.

Ulysses spacecraft in situ detections of cometary dust trails.

The Ulysses spacecraft was launched in 1990 and, after a Jupiter swing-by in 1992, became the first interplanetary spacecraft orbiting the Sun on a highly inclined trajectory with an inclination of [Formula: see text]. The spacecraft was equipped with an impact ionization dust detector which provided 17 years of in situ dust measurements in interplanetary space from 1990 to 2007. Cometary meteoroid streams (also referred to as trails) exist along the orbits of comets, forming fine structures of the interplanetary dust cloud. We use the Interplanetary Meteoroid Environment for eXploration (IMEX) dust streams in space model (Soja RH et al. 2015 Characteristics of the dust trail of 67P/Churyumov-Gerasimenko: an application of the IMEX model. Astron. Astrophys. 583, A18. (doi:10.1051/0004-6361/201526184)) to predict cometary stream traverses by Ulysses and re-analyse the Ulysses dust dataset in order to identify impacts of cometary stream particles detected during such trail traverses. We identify 19 particles compatible with three Ulysses trail traverses on 12 March 1995, 25-27 April 2001 and 16-19 May 2001. The particle origin is compatible with up to five comets, i.e. 10P/Tempel 2, 146P/Shoemaker-LINEAR, 267P/LONEOS and possibly 45P/Honda-Mrkos-Pajdušáková and P/1999 RO28 (LONEOS). We find a dust spatial density in these trails of approximately [Formula: see text]. The radii of the detected cometary stream particles derived from the dust instrument calibration are in the micrometre range. The in situ analysis of meteoroid trail particles in space, which can be traced back to their source bodies, opens a new opportunity for remote compositional analysis of comets and asteroids without the necessity to send a spacecraft to or even land on these celestial bodies, opening new opportunities for future space missions equipped with in situ dust analyzers. This article is part of the theme issue 'Dust in the Solar System and beyond'.

Statistical Reliability Analysis of Interplanetary Spacecraft Operating at Different Interplanetary Extremity

The success of spacecraft depends on their reliability which is influenced by the space environment complexity. This study examines the reliability of interplanetary spacecraft in different regions. We focused on spacecraft in interplanetary space, specifically within the interplanetary boundary. Our analysis excluded planetary landers, probes, and satellites in Earth orbit. 131 spacecraft were studied, with data on launch mass and lifespan recorded. Non-parametric analysis was initially conducted, followed by parametric analysis using the Weibull Distribution. Results showed higher reliability for spacecraft beyond Earth and Mars extremities. Analyses were also done based on spacecraft mass categories to validate reliability effects. Various factors affecting reliability were discussed, including design, integration processes, testing methodologies, and mass constraints.

How Does the Structure of Rarefaction Regions Develop?

A rarefaction region (RR) occurs at the trailing edge of the fast solar wind stream. It comes from an area of small longitudinal extent on the solar surface and exhibits a fine and complex structure. In our study, we did a superposed epoch analysis of the proton and α parameters across the RR and observed their gradual evolution. We did not find any clear boundary between the fast and slow solar winds inside the RR because a majority of our observations show that most of the RR plasma corresponds to the fast solar wind; only the α–proton drift velocity decreases from the beginning of the RR. We investigate different ways of its reduction in interplanetary space and show that this feature is likely associated with the mirroring of the multicomponent solar wind. Nevertheless, considering the observed solar wind characteristics and taking into account the mutual relations between the proton and α parameters, we define the composition boundary where the α relative abundance and α–proton temperature ratio change abruptly from the values typical for the fast wind toward slow wind values. This boundary is the most probable candidate for the stream interface. Based on these findings, we speculate that the RR formation starts already near the Sun and formulate two possible scenarios.

Atacama Large Aperture Submillimeter Telescope (AtLAST) science: Planetary and cometary atmospheres.

The study of planets and small bodies within our Solar System is fundamental for understanding the formation and evolution of the Earth and other planets. Compositional and meteorological studies of the giant planets provide a foundation for understanding the nature of the most commonly observed exoplanets, while spectroscopic observations of the atmospheres of terrestrial planets, moons, and comets provide insights into the past and present-day habitability of planetary environments, and the availability of the chemical ingredients for life. While prior and existing (sub)millimeter observations have led to major advances in these areas, progress is hindered by limitations in the dynamic range, spatial and temporal coverage, as well as sensitivity of existing telescopes and interferometers. Here, we summarize some of the key planetary science use cases that factor into the design of the Atacama Large Aperture Submillimeter Telescope (AtLAST), a proposed 50-m class single dish facility: (1) to more fully characterize planetary wind fields and atmospheric thermal structures, (2) to measure the compositions of icy moon atmospheres and plumes, (3) to obtain detections of new, astrobiologically relevant gases and perform isotopic surveys of comets, and (4) to perform synergistic, temporally-resolved measurements in support of dedicated interplanetary space missions. The improved spatial coverage (several arcminutes), resolution (~ 1.2'' - 12''), bandwidth (several tens of GHz), dynamic range (~ 10 5) and sensitivity (~ 1 mK km s -1) required by these science cases would enable new insights into the chemistry and physics of planetary environments, the origins of prebiotic molecules and the habitability of planetary systems in general.

Dynamics of Sputtered Neutral Sodium Atoms in the Near‐Mercury Space

AbstractThe solar wind sputtering in the magnetospheric polar cusp is an important source of heavy atoms in Mercury’s exosphere and magnetosphere. However, the majority of ejected atoms are neutral, undergoing an extended period before photoionization occurs. In this study, we employ an ab initio simulation to investigate the behavior of sodium (Na) atoms prior to their photoionization. Our results reveal that overall only approximately 2.7% of the sputtered atoms contribute to magnetospheric ions, while the vast majority of these ions (∼82.9%) escape into interplanetary space. The remaining fraction (14.4%) eventually returns to the planetary surface. For Na atoms ionized inside the magnetosphere, a larger proportion of Na+ (53.5%) is supplied to the magnetotail compared to the polar cusp (39.4%), which is due to the tailward acceleration caused by solar radiation. Additionally, the remaining Na+ (7.1%) contributes to the dayside ring current region, as demonstrated by the observation of the MErcury Surface, Space ENvironment, GEochemistry, and Ranging (MESSENGER) spacecraft. Our research introduces a perspective on Na+ transport in the magnetosphere that complements and coexists with traditional mechanisms.

Decentralization of Webhosting Using Blockchain

Abstract: This paper introduces a ground breaking approach to web hosting, leveraging the decentralized infrastructure of the Interplanetary File System (IPFS) and the secure framework of the Ethereum blockchain. In this novel model, Ethereum smart contracts govern the IPFS network and hosting services, ensuring robust management and reliability. IPFS serves as the backbone for data storage and website hosting, with storage miner nodes providing essential pinning services to maintain the longevity of stored content. Additionally, these nodes facilitate the Interplanetary Name Space (IPNS) service, enabling the creation and maintenance of mutable links to IPFS contents, thus enhancing accessibility and ease of updates. DNS TXT records are employed to seamlessly map domain names to IPNS addresses, seamlessly integrating decentralized hosting with traditional web infrastructure. To safeguard user privacy and data integrity, websites deploy advanced encryption algorithms, ensuring confidentiality in the decentralized environment. The proposed model represents a fusion of IPFS and blockchain networks, creating a robust platform for decentralized web hosting. Through experimental validation, the efficacy of this approach is demonstrated, showcasing its ability to surpass traditional web hosting models in terms of confidentiality, integrity, and availability. This innovative solution not only addresses existing challenges in web hosting but also heralds a new era of decentralized, secure, and accessible internet infrastructure

High-energy particle observations from the Moon.

The Moon is a unique natural laboratory for the study of the deep space plasma and energetic particles environment. During more than 3/4 of its orbit around the Earth it is exposed to the solar wind. Being an unmagnetized body and lacking a substantial atmosphere, solar wind and solar energetic particles bombardthe Moon's surface, interacting with the lunar regolith and the tenuous lunar exosphere. Energetic particles arriving at the Moon's surface can be absorbed, or scattered, or can remove another particle from the lunar regolith by sputtering or desorption. A similar phenomenon occurs also with the galactic cosmic rays, which have fluxes and energy spectra representative of interplanetary space. During the remaining part of its orbit the Moon crosses the tail of the terrestrial magnetosphere. It then provides the opportunity to study in-situ the terrestrial magnetotail plasma environment as well as atmospheric escape from the Earth's ionosphere, in the form of heavy ions accelerated and streaming downtail. The lunar environment is thus a unique natural laboratory for analysing the interaction of the solar wind, the cosmic rays and the Earth's magnetosphere with the surface, the immediate subsurface, and the surface-bounded exosphere of an unmagnetized planetary body. This article is part of a discussion meeting issue 'Astronomy from the Moon: the next decades (part 2)'.

Ground Electric Field, Atmospheric Weather and Electric Grid Variations in Northeast Greece Influenced by the March 2012 Solar Activity and the Moderate to Intense Geomagnetic Storms

In a recent paper, we extended a previous study on the solar solar influence to the generation of the March 2012 heatwave in the northeastern USA. In the present study we check the possible relationship of solar activity with the early March 2012 bad weather in northeast Thrace, Greece. To this end, we examined data from various remote sensing instrumentation monitoring the Sun (SDO satellite), Interplanetary space (ACE satellite), the Earth’s magnetosphere (Earth-based measurements, NOAA-19 satellite), the top of the clouds (Terra and Aqua satellites), and the near ground atmosphere. Our comparative data analysis suggests that: (i) the winter-like weather (rainfall, fast winds, decreased temperature) in Thrace started on 6 March 2012, the same day as the heatwave started in USA, (ii) during the March 2012 winter-like event in Thrace (6–15 March), the ACE satellite recorded enhanced fluxes of solar energetic particles (SEPs), while SOHO and PAMELA recorded solar protons at very high energies (&gt;500 MeV), (iii) Between 3–31 March, the temperature in Alexandoupoli and the ACE/EPAM solar high energy (1.88–4.70 MeV) proton flux were strongly anticorelated (r = −0.75, p = 0.5). (iv) Thrace experienced particularly intense cyclonic circulation, during periods of magnetic storms on 8–10 and 12–13 March, which occurred after the arrival at ACE of two interplanetary shock waves, on March 8 and March 11, respectively, (v) at the beginning of the two above mentioned periods large atmospheric electric fields were recorded, with values ranging between ~−2000 V/m and ~1800 V/m on 8 March, (vi) the winter-like weather on 8–10 March 2012 occurred after the detection of the main SEP event related with a coronal mass ejection released in interplanetary space as a result of intense solar flare activity observed by SDO on 7 March 2012, (vi) the 8–10 March weather was related with a deep drop of ~63 °C in the cloud top temperature measured by MODIS/Terra, which favors strong precipitation. Finally, we analyzed data from the electric power network in Thrace (~41°N) and we found, for the first time sudden voltage changes of ~3.5 kV in the electric grid in Greece, during the decay phase of the March 2012 storm series. We discuss the winter-like March 2012 event in Thrace regarding the influence of solar cosmic rays on the low troposphere mediated by positive North Atlantic Oscillation (NAO). Finally, we infer that the novel finding of the geomagnetic effects on the electric power grid in Thrace may open a new window into space weather applications research.