2026-03-02 | | Total: 3
Intense solar eruptions occasionally accelerate solar energetic particles (SEPs) and can trigger ground-level enhancements (GLEs). Among the 77 known GLEs, the first four GLEs, #1 -- 4 in the 1940s took place before the advent of the standard neutron monitors and were missing from the International GLE Database. This data gap challenged their quantification. To overcome this difficulty, we systematically gathered, digitised, and quantified contemporaneous cosmic-ray records pertaining to these GLEs. These data allow us to reconstruct the temporal evolution, with the 1 -- 15 min resolutions, of these GLEs, and broaden their geographical coverage to a global scale. GLEs #1 and #3 exhibited gradual increases in their rise times, measured at 45 +/- 15 and 105 +/- 15 min, respectively. In contrast, GLEs #2 and #4 both exhibited abrupt increases of 15 +/- 15 min. We also compared integral ionisation increase on the standard ionisation chambers and their local geomagnetic cutoff rigidities Pc to qualitatively compare these GLE's spectral hardness: Our result indicates that their spectra are extremely hard for GLEs #2 and #4 and mildly hard for GLEs #1 and #3. GLE #3 showed the greatest integral ionisations for polar detectors among them.
We investigate the impact of a Coronal Mass Ejection (CME) on the transport and acceleration of relativistic protons in the solar wind using a coupled 3D Magnetohydrodynamics (MHD) simulation and a test-particle approach. The CME is driven by a spheromak injected into a Parker solar wind at a heliocentric distance of 0.139 AU. The trajectories of 5 GeV protons, injected towards the CME from 3 AU, are integrated in the guiding-centre approximation limit and scattered in velocity space with a mean free path $λ_{\|}$. Our results show that the CME can increment the protons energy by several GeV. The acceleration occurs during the time particles stream along the portion of a magnetic field line subject to compression downstream of the quasi-perpendicular portion of the CME-driven shock. In our configuration, the maximum energy gain, which is of the order of a few percent per shock crossing, occurs when the shock approaches 0.3 AU. Large energy gains require multiple passes through the acceleration region, which is made possible by the combined action of the mirror force and pitch angle scattering. The efficiency of the acceleration on time scales of the order of hours scales as $λ_{\|}^{-3/2}$. Energy spectra harden for decreasing parallel mean free path $λ_{\|}$.
Magnetic reconnection in partially ionized plasmas plays a crucial role in a wide range of solar, astrophysical, and laboratory environments. While reconnection in such plasmas is commonly characterized by the ion-neutral coupling strength and the ionization fraction $χ=n_{i}/(n_{i}+n_{n})$, most previous studies have focused primarily on the former. A systematic exploration of the ionization fraction, particularly in combination with ion-neutral coupling, is still lacking. This study presents the first systematic scan of the two-dimensional parameter space defined by ion-neutral collisionality and ionization fraction, enabling investigation of the transition from strongly coupled reconnection to faster, decoupled reconnection. To achieve this, we employ a new three-fluid, five-moment numerical model that treats electrons, ions, and neutrals as separate species on an equal footing. We find that in the strongly coupled regime, the reconnection rate is consistent with a $χ^{1/4}$ scaling. As collisionality decreases, the system transitions to a fast, ionization-independent regime. On the other hand, in all simulations, the current sheet thins down to the ion inertial length di, rather than the expanded hybrid scale $d_{i}χ^{-1/2}$ predicted by analytic fluid theories. The identified critical thickness and the resulting onset of fast reconnection agree reasonably well with recent fully kinetic simulations and laboratory experiments. In addition, we show that, over a wide range of coupling strengths, the ion outflow velocities remain Alfvénic, scaling with the appropriate ion or hybrid Alfvén speed, while the hybrid outflow velocity scales as $χ^{1/2}$ when normalized by ion Alfvén speed.