Introduction

Cosmic ray modulation refers to the temporal, spatial, and energy-dependent variation in the intensity of galactic cosmic rays (GCRs) as they propagate through the heliosphere before reaching Earth. Unlike the local interstellar spectrum (LIS), which represents the undisturbed galactic cosmic ray intensity outside the heliosphere, the flux observed near Earth is significantly modified by interactions with the expanding solar wind and the heliospheric magnetic field (HMF) (Parker, 1965; Gleeson & Axford, 1968; Potgieter, 2013). These interactions involve four principal transport mechanisms spatial diffusion, solar wind convection, gradient and curvature drifts, and adiabatic energy changes which collectively determine the energy-dependent modulation of cosmic rays throughout the heliosphere (Jokipii, 1971; Potgieter, 2013).

The heliosphere, extending from the solar corona to the heliopause at approximately 120 AU acts as a dynamic magnetic shield that continuously modifies the propagation of energetic charged particles. The outward-flowing solar wind transports the frozen-in Parker spiral magnetic field, creating a turbulent plasma environment through which galactic cosmic rays must travel before reaching the inner solar system (Parker, 1958, 1965). The intensity of cosmic rays observed near Earth therefore depends not only on their galactic origin and energy spectrum but also on the evolving properties of the heliospheric magnetic field, solar wind speed, magnetic turbulence, and large-scale heliospheric structures such as the heliospheric current sheet (Burger et al., 2000; Potgieter, 2013).

One of the most striking characteristics of cosmic ray modulation is its strong inverse relationship with solar activity. Long-term observations from neutron monitor stations demonstrate that cosmic ray intensity reaches maximum values during solar minimum conditions, when heliospheric turbulence is relatively weak, and decreases during solar maxima due to enhanced solar wind disturbances, stronger interplanetary magnetic fields, frequent coronal mass ejections (CMEs), and increased magnetic turbulence (Forbush, 1954; Lockwood, 1971; Potgieter, 2013). Superimposed on the approximately 11-year solar cycle is the 22-year Hale magnetic cycle, during which periodic reversals of the Sun's global magnetic polarity produce charge-sign-dependent drift effects that significantly influence cosmic ray transport (Jokipii et al., 1977; Burger & Hattingh, 1998).

The modulation process occurs over multiple temporal scales ranging from hours to decades. Short-term variations include recurrent 27-day periodicities associated with solar rotation, sudden intensity decreases known as Forbush decreases following the passage of interplanetary shocks and magnetic clouds, and transient increases associated with solar energetic particle (SEP) events (Cane, 2000; Richardson, 2004). Long-term modulation reflects changes in heliospheric conditions over successive solar cycles and provides valuable insights into the coupling between solar magnetic activity, heliospheric plasma dynamics, and galactic cosmic ray transport (Potgieter, 2013).

The theoretical description of cosmic ray modulation is based primarily on the Parker Transport Equation (PTE), which combines the effects of anisotropic spatial diffusion, convection by the solar wind, gradient and curvature drifts, and adiabatic energy losses into a single transport framework (Parker, 1965). During the past six decades, this equation has been extended through stochastic differential equation (SDE) methods, finite-difference techniques, and three-dimensional magnetohydrodynamic (MHD) simulations, significantly improving our understanding of cosmic ray propagation throughout the heliosphere (Zhang, 1999; Pei et al., 2010; Strauss et al., 2011).

Ground-based neutron monitor observations and spacecraft missions have revolutionized the study of cosmic ray modulation. Long-term measurements from global neutron monitor networks, together with observations from Voyager 1 and 2, Ulysses, ACE, SOHO, STEREO, PAMELA, AMS-02, Parker Solar Probe, Solar Orbiter, and Aditya-L1, provide continuous measurements of energetic particles, solar wind plasma, and heliospheric magnetic fields over multiple solar cycles (Stone et al., 2013; Aguilar et al., 2015; Potgieter et al., 2014). These observations have enabled significant improvements in numerical transport models and have substantially refined estimates of diffusion coefficients, particle drift velocities, heliospheric boundary conditions, and modulation during Solar Cycles 20–25.

Recent investigations have increasingly incorporated artificial intelligence (AI) and machine learning (ML) techniques into heliophysics and space weather research. Deep learning models, including long short-term memory (LSTM) networks, convolutional neural networks (CNNs), transformer architectures, random forests, gradient boosting algorithms, and physics-informed neural networks (PINNs), have demonstrated considerable promise for forecasting neutron monitor counts, identifying Forbush decreases, reconstructing heliospheric parameters, and improving operational space weather prediction (Camporeale, 2019; Nishizuka et al., 2018). Integrating these data-driven approaches with classical transport theory offers new opportunities for developing accurate and computationally efficient forecasting systems.

The scientific significance of cosmic ray modulation extends well beyond heliophysics. Galactic cosmic rays contribute substantially to radiation exposure experienced by astronauts, high-altitude aviation crews, and spacecraft electronics; influence atmospheric ionization; produce cosmogenic isotopes such as ¹⁰Be and ¹⁴C used in paleoclimate reconstruction; and provide important diagnostics of heliospheric plasma turbulence and solar-terrestrial interactions (Usoskin, 2017). Consequently, understanding cosmic ray modulation is essential for astrophysics, planetary exploration, radiation protection, climate studies, and operational space weather forecasting.

Several important contributions to understanding cosmic ray modulation, particularly its relationship with solar activity and heliospheric parameters, have been made by Agarwal and Mishra. Their studies investigated the modulation of galactic cosmic rays over multiple solar cycles, examined rigidity-dependent modulation, and analyzed the influence of interplanetary parameters on recurrent 27-day cosmic ray variations, providing valuable observational evidence for understanding long-term heliospheric modulation processes (Agarwal & Mishra, 2013a, 2013b, 2023).

The present review synthesizes current knowledge of cosmic ray modulation by integrating classical transport theory, heliospheric plasma physics, observational datasets, numerical simulations, and recent developments in artificial intelligence. Particular emphasis is placed on the physical mechanisms governing modulation, variability during Solar Cycles 20–25, state-of-the-art computational models, and future opportunities for AI-enabled heliospheric forecasting and space weather applications.

1. Historical Evolution

The history of cosmic ray research spans more than a century and reflects remarkable progress in atmospheric physics, nuclear physics, astrophysics, plasma physics, and heliophysics. What began as an investigation into the unexplained ionization of Earth's atmosphere has evolved into a multidisciplinary field that investigates the origin, acceleration, propagation, and modulation of energetic charged particles throughout the Universe. Continuous improvements in experimental techniques, theoretical modeling, and space-based observations have significantly advanced our understanding of cosmic ray transport and its relationship with solar activity and heliospheric dynamics (Simpson, 2000; Potgieter, 2013).

The foundation of cosmic ray research was established in 1912, when Victor Franz Hess conducted a series of high-altitude balloon experiments to investigate the source of atmospheric ionization. Hess observed that the ionization rate initially decreased with altitude but subsequently increased rapidly above approximately 3 km, demonstrating that highly penetrating radiation was entering Earth's atmosphere from outer space rather than originating solely from terrestrial radioactivity (Hess, 1912). These pioneering observations earned Hess the 1936 Nobel Prize in Physics and established the extraterrestrial origin of cosmic rays. Subsequent balloon measurements by Werner Kolhörster confirmed Hess's findings and extended ionization measurements to even greater altitudes, providing compelling evidence that cosmic rays were a pervasive phenomenon throughout Earth's atmosphere (Kolhörster, 1914).

During the 1920s and 1930s, the invention of the Geiger–Müller counter revolutionized cosmic ray detection by allowing individual energetic particles to be counted with unprecedented accuracy. The coincidence techniques developed by Walther Bothe and Bruno Rossi demonstrated that cosmic rays consist primarily of charged particles rather than electromagnetic radiation (Rossi, 1930). At the same time, cloud chambers and nuclear emulsions became indispensable tools for particle detection, leading to the discovery of several fundamental particles, including the positron by Carl Anderson (1932), the muon, and later the pion, thereby establishing cosmic rays as one of the earliest laboratories for high-energy particle physics (Anderson, 1933; Powell et al., 1947).

The emergence of nuclear and particle physics during the mid-twentieth century shifted scientific interest from merely detecting cosmic rays to understanding their origin and acceleration mechanisms. Enrico Fermi proposed the stochastic acceleration mechanism in 1949, introducing the concept that charged particles could gain energy through repeated interactions with moving magnetic clouds, a process now known as Fermi acceleration (Fermi, 1949). This theory was later extended into the highly efficient diffusive shock acceleration (DSA) mechanism, which is now widely accepted as the principal acceleration process operating in supernova remnants and other astrophysical shock environments (Bell, 1978; Blandford & Ostriker, 1978). These developments established the theoretical framework for explaining the observed power-law energy spectrum of galactic cosmic rays.

A major breakthrough in heliophysics occurred with Eugene N. Parker's revolutionary theory of the solar wind in 1958. Parker demonstrated that the million-degree solar corona continuously expands into interplanetary space as a supersonic plasma flow carrying the solar magnetic field (Parker, 1958). This discovery fundamentally transformed the understanding of cosmic ray propagation because it established that galactic cosmic rays travel through a dynamic, magnetized heliosphere rather than empty space. Building upon this work, Parker formulated the Parker Transport Equation in 1965, incorporating spatial diffusion, solar wind convection, particle drifts, and adiabatic energy changes into a unified mathematical framework that remains the cornerstone of modern cosmic ray modulation theory (Parker, 1965).

The International Geophysical Year (1957–1958) marked another milestone with the establishment of global neutron monitor networks capable of continuously recording secondary particles produced by cosmic ray interactions in Earth's atmosphere. These observations revealed the pronounced anti-correlation between cosmic ray intensity and solar activity and led to the identification of recurrent 27-day variations, Forbush decreases, and long-term modulation associated with the approximately 11-year solar cycle (Simpson, 2000; Forbush, 1954). Ground-based neutron monitor observations have now accumulated over seven decades of continuous records, providing one of the longest datasets available for studying heliospheric variability and space weather.

The space age revolutionized cosmic ray research by enabling direct observations beyond Earth's atmosphere. Early spacecraft such as Explorer, IMP, and Pioneer provided the first in situ measurements of energetic particles in interplanetary space. Later missions, including Helios, ISEE, Voyager 1 and Voyager 2, and Ulysses, significantly expanded understanding of cosmic ray propagation throughout the heliosphere (Stone et al., 2013; McDonald et al., 2003). The Voyager missions provided particularly important discoveries by directly measuring galactic cosmic rays near the termination shock and, eventually, beyond the heliopause, thereby establishing the local interstellar spectrum (LIS) that serves as the outer boundary condition for modulation models (Stone et al., 2013).

During the twenty-first century, several advanced space missions transformed heliophysical research. ACE, SOHO, STEREO, PAMELA, AMS-02, Parker Solar Probe, Solar Orbiter, and India's Aditya-L1 mission have provided unprecedented measurements of solar wind plasma, energetic particles, magnetic fields, and heliospheric structures across multiple heliocentric distances (Aguilar et al., 2015; Fox et al., 2016; Müller et al., 2020). These missions have enabled detailed investigations of diffusion coefficients, magnetic turbulence, particle drifts, solar energetic particle events, and the influence of coronal mass ejections on cosmic ray modulation. Together, they have significantly refined numerical transport models and enhanced operational space weather forecasting capabilities.

Parallel advances in computational science have transformed theoretical investigations of cosmic ray transport. Early analytical solutions have been supplemented by finite-difference methods, Monte Carlo simulations, stochastic differential equation (SDE) approaches, and three-dimensional magnetohydrodynamic (MHD) heliospheric models capable of simulating realistic solar wind conditions and particle propagation (Zhang, 1999; Strauss et al., 2011). More recently, artificial intelligence (AI), machine learning (ML), deep learning, physics-informed neural networks (PINNs), and digital-twin technologies have emerged as powerful tools for analyzing large heliophysical datasets, forecasting neutron monitor counts, predicting Forbush decreases, and reconstructing heliospheric transport parameters (Camporeale, 2019). These developments are opening new opportunities for integrating first-principles physics with data-driven methodologies to improve the accuracy and computational efficiency of cosmic ray forecasting.

Significant contributions to the understanding of solar modulation have also been made by Agarwal and Mishra, whose studies investigated the long-term modulation of galactic cosmic rays over successive solar cycles, rigidity-dependent modulation, and the influence of interplanetary parameters on recurrent 27-day cosmic ray variations (Agarwal & Mishra, 2013a, 2013b, and 2023). Their work has strengthened the observational understanding of heliospheric modulation processes and their relationship with solar magnetic activity.

Today, cosmic ray research is a highly interdisciplinary field that integrates astrophysics, heliophysics, plasma physics, computational science, and artificial intelligence. With ongoing observations from Parker Solar Probe, Solar Orbiter, Aditya-L1, AMS-02, and future missions such as the Interstellar Mapping and Acceleration Probe (IMAP), together with advances in high-performance computing and AI-based forecasting, the field is entering a new era of precision heliophysics. These developments are expected to provide deeper insights into cosmic ray transport, heliospheric turbulence, space weather prediction, and the broader interactions between the Sun and the interstellar environment.

1. Origin and Classification of Cosmic Rays

Cosmic rays are highly energetic charged particles that originate from a wide range of astrophysical environments throughout the Universe. They constitute one of the most energetic components of the interstellar medium and encompass particles with kinetic energies extending from approximately 10⁶ eV (MeV) to beyond 10²⁰ eV (EeV). Since their discovery by Victor Hess in 1912, understanding the origin, acceleration mechanisms, chemical composition, and propagation of cosmic rays has remained a central challenge in astrophysics and heliophysics (Hess, 1912; Blasi, 2013). Unlike photons, cosmic rays are electrically charged and therefore undergo repeated deflections by galactic, interplanetary, and heliospheric magnetic fields, making the direct identification of their sources particularly difficult (Gaisser et al., 2016).

The observed cosmic ray spectrum represents a superposition of particles originating from different astrophysical environments. Based on their source regions and acceleration mechanisms, cosmic rays are generally classified into Galactic Cosmic Rays (GCRs), Solar Energetic Particles (SEPs), Anomalous Cosmic Rays (ACRs), and Extragalactic Cosmic Rays (EGCRs) (Potgieter, 2013). Each population exhibits distinct energy ranges, elemental compositions, acceleration processes, and modulation characteristics, thereby contributing differently to the radiation environment of the heliosphere.

• 1. Galactic Cosmic Rays (GCRs)

Galactic Cosmic Rays (GCRs) constitute the dominant component of cosmic rays observed near Earth under quiet solar conditions. They originate primarily within the Milky Way Galaxy and are believed to be accelerated predominantly in supernova remnants (SNRs) through the mechanism of diffusive shock acceleration (DSA), also known as first-order Fermi acceleration (Bell, 1978; Blandford & Ostriker, 1978). Additional potential sources include pulsars, magnetars, stellar winds, super bubbles, and binary star systems, all of which contribute to the acceleration of charged particles to relativistic energies (Blasi, 2013).

GCRs account for approximately 85–90% protons (hydrogen nuclei), 9–12% alpha particles (helium nuclei), and nearly 1% heavier nuclei, together with trace amounts of electrons, positrons, antiprotons, and ultra-heavy nuclei (Simpson, 1983; Gaisser et al., 2016). Their energies typically range from hundreds of MeV to about 10¹⁷ eV, although the highest-energy Galactic particles are thought to approach the "knee" of the cosmic ray spectrum near 3 × 10¹⁵ eV. Upon entering the heliosphere, GCRs undergo significant modulation caused by solar wind convection, diffusion, particle drifts, and adiabatic energy losses, making them the primary focus of heliospheric cosmic ray modulation studies (Potgieter, 2013).

Solar Energetic Particles (SEPs)

Solar Energetic Particles (SEPs) originate within the solar atmosphere and are accelerated during explosive solar events such as solar flares and coronal mass ejections (CMEs). Unlike Galactic Cosmic Rays, SEPs exhibit relatively lower energies, typically ranging from several keV to a few GeV, but can reach relativistic energies during exceptionally intense solar events (Reames, 2013).

SEPs are generally classified into impulsive and gradual events. Impulsive SEP events are associated primarily with magnetic reconnection during solar flares and are characterized by enhanced abundances of electrons and heavy ions. Gradual SEP events are produced by CME-driven shock waves propagating through the corona and interplanetary space, resulting in large proton-dominated particle populations capable of producing significant radiation hazards for spacecraft and astronauts (Reames, 1999; Desai & Giacalone, 2016). These particles are of particular importance in space weather because they can reach Earth within tens of minutes following major solar eruptions.

Anomalous Cosmic Rays (ACRs)

Anomalous Cosmic Rays represent a unique population that originates from neutral atoms in the local interstellar medium (LISM). Neutral atoms such as hydrogen, helium, oxygen, and nitrogen penetrate deep into the heliosphere because they are unaffected by magnetic fields. Once inside the heliosphere, these atoms become ionized through photoionization, charge exchange, or electron impact, forming pickup ions that are subsequently transported outward by the solar wind (Fisk et al., 1974).

These pickup ions are believed to undergo further acceleration at the termination shock, where the supersonic solar wind slows abruptly before entering the heliosheath, producing anomalous cosmic rays with energies typically between 10 and 100 MeV per nucleon (Stone et al., 2005). Measurements by the Voyager spacecraft have significantly improved our understanding of ACR acceleration and demonstrated that the acceleration process extends beyond the termination shock into the heliosheath (Stone et al., 2005; Decker et al., 2008).

Extragalactic Cosmic Rays (EGCRs)

Extragalactic Cosmic Rays (EGCRs) constitute the highest-energy component of the cosmic ray spectrum, generally exceeding 10¹⁸ eV. Their enormous energies exceed the acceleration capabilities of conventional Galactic sources, suggesting origins in extremely energetic astrophysical environments such as active galactic nuclei (AGN), gamma-ray bursts (GRBs), radio galaxies, and galaxy clusters (Kotera & Olinto, 2011).

Ultra-high-energy cosmic rays (UHECRs), with energies approaching 10²⁰ eV, remain among the greatest mysteries of modern astrophysics. Observatories such as the Pierre Auger Observatory and the Telescope Array continue investigating their origin, composition, and arrival directions. At these extreme energies, cosmic rays interact with photons of the cosmic microwave background (CMB), leading to the Greisen–Zatsepin–Kuzmin (GZK) cutoff, which limits their propagation distances to approximately 100 Mpc (Greisen, 1966; Zatsepin & Kuzmin, 1966). Consequently, observed ultra-high-energy cosmic rays must originate from relatively nearby extragalactic sources.

Chemical Composition of Cosmic Rays

The chemical composition of cosmic rays provides valuable information regarding their astrophysical origin and acceleration history. Direct measurements by balloon experiments and space missions such as AMS-02, PAMELA, and ACE indicate that approximately 89% of primary cosmic rays consist of protons, 10% are helium nuclei (alpha particles), and about 1% comprise heavier nuclei extending from lithium to uranium (Aguilar et al., 2015; Gaisser et al., 2016). Small fractions of electrons, positrons, antiprotons, and rare isotopes are also present.

The abundances of secondary nuclei such as lithium, beryllium, and boron, which are produced by spallation reactions during propagation through the interstellar medium, provide critical diagnostics of cosmic ray transport. In particular, the boron-to-carbon (B/C) ratio is widely used to estimate Galactic diffusion coefficients and propagation timescales (Strong et al., 2007). Measurements of isotopic abundances also constrain nucleosynthesis processes occurring within supernovae and massive stars.

Acceleration Mechanisms

The enormous energies of cosmic rays require highly efficient astrophysical acceleration mechanisms. The most widely accepted process is diffusive shock acceleration (DSA) operating at supernova shocks, where particles repeatedly cross expanding shock fronts and gain energy during each crossing (Bell, 1978; Blandford & Ostriker, 1978). Other important acceleration mechanisms include stochastic (second-order Fermi) acceleration, magnetic reconnection, turbulent plasma acceleration, and acceleration within relativistic jets associated with active galactic nuclei and gamma-ray bursts (Fermi, 1949; Lazarian et al., 2020). These mechanisms collectively explain the broad energy spectrum and diverse composition of cosmic rays observed throughout the Universe.

Importance in Cosmic Ray Modulation

The classification of cosmic rays is fundamental for understanding heliospheric modulation because each particle population interacts differently with the heliospheric environment. Galactic Cosmic Rays experience substantial long-term modulation due to solar activity and therefore serve as the primary tracers of heliospheric transport processes. Solar Energetic Particles dominate short-duration space weather events and represent major radiation hazards, whereas Anomalous Cosmic Rays provide unique insights into heliospheric boundary physics. Extragalactic Cosmic Rays, owing to their extremely high rigidities, are only minimally influenced by heliospheric modulation but remain essential for understanding high-energy astrophysical acceleration mechanisms (Potgieter, 2013). Together, these particle populations provide a comprehensive picture of energetic particle acceleration, propagation, and interaction across multiple astrophysical environments.

Acceleration Mechanisms of Cosmic Rays

The extraordinary energies attained by cosmic rays, extending from approximately 10⁶ eV (MeV) to beyond 10²⁰ eV (EeV), require highly efficient astrophysical acceleration mechanisms capable of transferring enormous amounts of energy to charged particles. Since thermal particles possess energies of only a few electron volts, cosmic rays must undergo acceleration by several orders of magnitude before reaching the relativistic energies observed in nature. Understanding these acceleration mechanisms is fundamental to explaining the origin, energy spectrum, elemental composition, and spatial distribution of cosmic rays throughout the Universe (Blasi, 2013; Gaisser et al., 2016). Observational evidence indicates that no single mechanism is capable of explaining the complete cosmic ray energy spectrum. Instead, multiple acceleration processes operate in diverse astrophysical environments, including supernova remnants, pulsars, active galactic nuclei, gamma-ray bursts, stellar winds, magnetars, and solar eruptions (Bell, 1978; Blandford & Eichler, 1987).

First-Order Fermi Acceleration (Diffusive Shock Acceleration)

The most widely accepted mechanism responsible for the acceleration of Galactic Cosmic Rays is first-order Fermi acceleration, commonly known as Diffusive Shock Acceleration (DSA). Originally proposed by Fermi (1949) and later refined by Bell (1978) and Blandford and Ostriker (1978), this mechanism operates at strong astrophysical shock waves, particularly those generated during supernova explosions. Charged particles repeatedly scatter across the shock front through interactions with magnetic irregularities on both sides of the shock. During each shock crossing, particles gain kinetic energy because the converging plasma flows compress the particle trajectories.

For a non-relativistic strong shock, the average fractional energy gain per shock crossing is approximately

ΔE/E ≈ (4/3)(V_s/c)

where (V_s) is the shock velocity and (c) is the speed of light (Bell, 1978). Since the energy gain is proportional to the first power of the shock velocity, this process is referred to as first-order Fermi acceleration. Repeated shock crossings naturally produce the observed power-law energy spectrum

N(E)∝E−γ,

or equivalently,

N(E)=kE−γ, N(E)=kE^{-\gamma}, N(E)=kE−γ,

where the spectral index γ is typically close to 2.0–2.3 for strong shocks (Blandford & Eichler, 1987). Supernova remnants remain the leading candidates for accelerating Galactic Cosmic Rays up to the "knee" region of the cosmic ray spectrum (~3 × 10¹⁵ eV).

Second-Order Fermi (Stochastic) Acceleration

Before proposing shock acceleration, Enrico Fermi (1949) introduced a stochastic acceleration mechanism in which charged particles interact randomly with moving magnetic clouds or plasma turbulence. Unlike DSA, particles may encounter both approaching and receding magnetic irregularities, resulting in a much smaller average energy gain.

The average fractional energy gain is given by

ΔE/E​ ∝ (V/c​)²

where (V) represents the characteristic velocity of turbulent magnetic clouds. Because the energy gain depends on the square of the velocity ratio, stochastic acceleration is considerably less efficient than shock acceleration. Nevertheless, second-order Fermi acceleration remains important in highly turbulent astrophysical environments such as the interstellar medium, galaxy clusters, pulsar wind nebulae, and the solar corona (Schlickeiser, 2002).

Magnetic Reconnection

Magnetic reconnection has emerged as another important mechanism for accelerating charged particles in magnetized plasmas. During reconnection, oppositely directed magnetic field lines break and reconnect, rapidly converting stored magnetic energy into plasma heating, bulk flows, and particle acceleration. This mechanism is particularly effective in environments characterized by intense magnetic fields and high plasma conductivity, including solar flares, coronal mass ejections (CMEs), magnetars, pulsar wind nebulae, and accretion disks surrounding black holes (Lazarian et al., 2020).

Observations from spacecraft such as Parker Solar Probe and Solar Orbiter have demonstrated that magnetic reconnection plays a major role in accelerating Solar Energetic Particles (SEPs) during solar eruptions. Unlike shock acceleration, reconnection can accelerate particles over relatively short timescales, producing rapid bursts of energetic electrons and ions observed during major solar flare events (Desai & Giacalone, 2016).

Acceleration by Compact Astrophysical Objects

Some of the highest-energy cosmic rays are believed to originate in extreme astrophysical environments associated with compact objects possessing enormous gravitational and magnetic fields. Pulsars and magnetars, with magnetic field strengths reaching 10¹²–10¹⁵ G, accelerate charged particles through rapidly rotating magnetospheres and relativistic winds. Similarly, active galactic nuclei (AGN) accelerate particles within relativistic jets powered by accretion onto super massive black holes, while gamma-ray bursts (GRBs) generate ultra-relativistic shocks capable of accelerating particles to ultra-high energies exceeding 10²⁰ eV (Kotera & Olinto, 2011).

These extreme environments provide sufficient magnetic confinement and acceleration times to explain the origin of Ultra-High-Energy Cosmic Rays (UHECRs), whose energies exceed the capabilities of conventional Galactic accelerators.

Acceleration of Solar Energetic Particles

Solar Energetic Particles are accelerated primarily through two distinct mechanisms. Impulsive SEP events result from magnetic reconnection during solar flares, where electrons and heavy ions are rapidly energized within the solar corona. Gradual SEP events, in contrast, are accelerated at CME-driven shock fronts propagating through interplanetary space. These shocks efficiently accelerate ambient solar wind ions and suprathermal particles to energies ranging from several MeV to a few GeV (Reames, 2013).

SEP acceleration is of considerable practical importance because these particles constitute one of the primary sources of radiation hazards for astronauts, satellites, and deep-space missions. Modern space weather forecasting therefore relies heavily on understanding shock formation and particle acceleration during solar eruptions.

Cosmic Ray Energy Spectrum and Spectral Features

The differential energy spectrum of cosmic rays follows an approximate power-law over more than twelve orders of magnitude in energy,

J(E) ∝ E^−2.7

Indicating that higher-energy particles become progressively rarer. However, several distinct spectral features provide important clues regarding acceleration mechanisms and source populations.

• The Knee (~3 × 10¹⁵ eV): Generally attributed to the maximum acceleration energy achievable in Galactic supernova remnants.
• Second Knee (~5 × 10¹⁷ eV): Marks the transition between Galactic and extragalactic cosmic ray populations.
• The Ankle (~3 × 10¹⁸ eV): Represents the increasing dominance of extragalactic cosmic rays.
• GZK Cutoff (~5 × 10¹⁹ eV): Produced by interactions between ultra-high-energy cosmic rays and photons of the cosmic microwave background (Greisen, 1966; Zatsepin & Kuzmin, 1966).

These spectral features provide strong observational constraints on both acceleration models and propagation theories.

Relationship between Acceleration and Cosmic Ray Modulation

Although cosmic ray acceleration occurs predominantly outside the heliosphere, the initial energy spectrum generated by astrophysical accelerators strongly influences subsequent heliospheric modulation. Lower-energy Galactic Cosmic Rays undergo substantial modification through diffusion, solar wind convection, gradient and curvature drifts, and adiabatic energy losses, whereas ultra-high-energy particles experience negligible heliospheric modulation because of their extremely large rigidities (Potgieter, 2013). Consequently, accurate models of cosmic ray modulation require realistic source spectra derived from acceleration theory together with detailed descriptions of heliospheric transport processes.

Recent studies by Agarwal and Mishra have demonstrated that the observed modulation of galactic cosmic rays during Solar Cycles 20–24 depends not only on heliospheric transport but also on the rigidity distribution of the incident Galactic Cosmic Ray spectrum, emphasizing the importance of integrating acceleration theory with modulation modeling (Agarwal & Mishra, 2013a, 2013b, 2023).

Heliosphere and Solar Wind

The heliosphere is a large magnetized plasma cavity generated by the continuous expansion of the solar wind from the Sun into interstellar space. It represents the dominant environment through which galactic cosmic rays (GCRs) propagate before reaching Earth and therefore plays a central role in determining the intensity, energy spectrum, and temporal variability of cosmic rays observed throughout the solar system. The heliosphere extends from the solar corona to the heliopause, located at approximately 120–130 AU, where the outward dynamic pressure of the solar wind balances the pressure of the local interstellar medium (LISM) (Parker, 1958; Stone et al., 2013). Within this vast region, charged particles interact continuously with the turbulent heliospheric magnetic field (HMF), solar wind plasma, and large-scale heliospheric structures, producing the phenomenon known as solar modulation of cosmic rays (Potgieter, 2013).

The concept of the heliosphere originated from Parker's (1958) solar wind theory, which demonstrated that the million-degree solar corona cannot remain in hydrostatic equilibrium but instead expands continuously as a supersonic plasma flow carrying the Sun's magnetic field into interplanetary space. As the Sun rotates, the frozen-in magnetic field forms the characteristic Parker spiral, producing a large-scale magnetic geometry that governs the propagation, scattering, and drift of energetic charged particles (Parker, 1958; Parker, 1965). Consequently, the heliosphere acts as a dynamic magnetic shield that substantially reduces the intensity of low- and intermediate-energy Galactic Cosmic Rays entering the inner solar system.

Structure of the Heliosphere

The heliosphere comprises several distinct plasma regions, each characterized by unique physical conditions that influence cosmic ray transport.

Inner Heliosphere

The inner heliosphere, extending from the solar corona to approximately 80–90 AU, contains the supersonic solar wind and the Parker spiral magnetic field. Within this region, the solar wind velocity ranges between 300 and 800 km s^-1, while the heliospheric magnetic field strength decreases approximately as the inverse square of heliocentric distance (Parker, 1958). Most modulation of galactic cosmic rays occurs within this region because particles undergo repeated diffusion, convection, gradient and curvature drifts, and adiabatic cooling as they propagate toward Earth (Jokipii, 1971; Potgieter, 2013).

Termination Shock

The termination shock represents the boundary where the supersonic solar wind abruptly slows to subsonic speeds owing to its interaction with the surrounding interstellar medium. Voyager observations indicate that the termination shock is located at approximately 84–94 AU, although its position varies with solar activity (Stone et al., 2005). This region is believed to accelerate Anomalous Cosmic Rays (ACRs) through diffusive shock acceleration and contributes significantly to energetic particle populations within the outer heliosphere (Decker et al., 2008).

Heliosheath

Beyond the termination shock lies the heliosheath, a turbulent region of compressed plasma situated between the termination shock and the heliopause. The heliosheath exhibits enhanced magnetic turbulence, increased plasma density, and reduced solar wind velocity, making it an important region for particle scattering and modulation (Stone et al., 2013). Measurements from Voyager 1 and Voyager 2 have demonstrated that cosmic ray intensities increase substantially while traversing the heliosheath because of reduced solar wind convection and changing heliospheric boundary conditions.

Heliopause

The heliopause constitutes the outer boundary of the heliosphere, separating the solar wind plasma from the local interstellar medium. Voyager 1 crossed the heliopause in 2012, providing the first direct measurements of the local interstellar spectrum (LIS) of galactic cosmic rays (Stone et al., 2013). Beyond this boundary, cosmic ray intensities approach their unmodulated interstellar values, making the heliopause the natural outer boundary condition for heliospheric transport models (Potgieter, 2013).

Local Interstellar Medium (LISM)

Outside the heliopause lies the local interstellar medium, composed of dilute plasma, neutral atoms, interstellar magnetic fields, and galactic cosmic rays. Neutral atoms from the LISM penetrate the heliosphere and become ionized through charge exchange and photoionization, forming pickup ions that later contribute to the production of Anomalous Cosmic Rays (Fisk et al., 1974). Consequently, the interaction between the heliosphere and the LISM plays a significant role in defining the boundary conditions for cosmic ray modulation.

Solar Wind

The solar wind is a continuous outward flow of ionized plasma consisting primarily of electrons, protons, alpha particles, and trace heavy ions emitted from the solar corona. Parker's solar wind theory demonstrated that thermal pressure gradients accelerate the coronal plasma to supersonic velocities, carrying the heliospheric magnetic field throughout the solar system (Parker, 1958). At Earth's orbit, typical solar wind properties include velocities of 300–800 km s^-1, proton densities of 3–10 cm^-3, temperatures near 10⁵ K, and magnetic field strengths of 3–10 nT (Schwenn, 2006).

Two principal solar wind regimes are recognized:

• Slow Solar Wind (300–450 km s^-1): Originates primarily near the heliospheric current sheet and active regions, exhibiting higher densities, stronger magnetic fluctuations, and greater variability.
• Fast Solar Wind (600–800 km s^-1): Emanates from coronal holes and is characterized by lower densities, reduced turbulence, and relatively stable plasma conditions (Schwenn, 2006).

Interactions between fast and slow solar wind streams generate Corotating Interaction Regions (CIRs), which produce recurrent disturbances responsible for the observed 27-day periodic modulation of galactic cosmic rays (Richardson, 2004).

Heliospheric Magnetic Field

The heliospheric magnetic field (HMF) is transported outward by the expanding solar wind while remaining magnetically connected to the rotating Sun. Solar rotation transforms the magnetic field into the characteristic Parker spiral, producing large-scale curvature and gradients that govern particle drift motions (Parker, 1958). The polarity of the HMF reverses approximately every 11 years, giving rise to the 22-year Hale magnetic cycle, which strongly influences charge-sign-dependent modulation of cosmic rays (Jokipii et al., 1977; Potgieter, 2013).

Large-scale heliospheric structures, including the heliospheric current sheet (HCS), play a critical role in cosmic ray transport. During solar minimum, the HCS remains relatively flat, allowing efficient particle drifts, whereas during solar maximum the sheet becomes highly warped, increasing particle scattering and reducing drift effects (Burger & Hattingh, 1998).

Solar Wind Turbulence and Cosmic Ray Transport

The solar wind is inherently turbulent across a broad range of spatial and temporal scales. Magnetic fluctuations generated by Alfvén waves, turbulence cascades, coronal mass ejections (CMEs), stream interaction regions, and plasma instabilities scatter energetic charged particles, thereby regulating the diffusion coefficients employed in cosmic ray transport models (Matthaeus et al., 2003).

Solar wind turbulence directly influences all four principal transport mechanisms described by the Parker Transport Equation:

• Spatial Diffusion: Magnetic irregularities scatter charged particles, reducing their effective mean free paths.
• Solar Wind Convection: The outward plasma flow convects cosmic rays away from the Sun.
• Gradient and Curvature Drifts: Large-scale magnetic field geometry produces systematic particle drifts whose magnitude depends upon particle charge and rigidity.
• Adiabatic Energy Changes: Expansion of the solar wind causes energetic particles to lose energy while propagating through the heliosphere (Parker, 1965; Potgieter, 2013).

During solar maximum, enhanced magnetic turbulence associated with frequent CMEs and interplanetary shocks significantly reduces diffusion coefficients, producing stronger modulation and lower galactic cosmic ray intensities at Earth. Conversely, quieter solar minimum conditions permit more efficient diffusion and correspondingly higher cosmic ray fluxes (Potgieter, 2013).

Space Missions Investigating the Heliosphere

Understanding of the heliosphere has advanced dramatically through observations from major space missions. Voyager 1 and Voyager 2 provided the first measurements of the termination shock, heliosheath, and heliopause (Stone et al., 2013). Ulysses revealed the three-dimensional structure of the heliosphere and solar wind at high heliographic latitudes. ACE, SOHO, and STEREO supplied continuous observations of solar wind plasma, energetic particles, and coronal mass ejections. More recently, Parker Solar Probe, Solar Orbiter, and Aditya-L1 have enabled detailed investigations of the near-Sun solar wind, magnetic reconnection, coronal heating, and heliospheric magnetic field evolution, substantially improving our understanding of cosmic ray transport and space weather processes (Fox et al., 2016; Müller et al., 2020).

Recent investigations by Agarwal and Mishra (2013a, 2013b, 2023) further demonstrated that long-term galactic cosmic ray modulation is closely linked to variations in solar wind velocity, interplanetary magnetic field strength, heliospheric magnetic polarity, and recurrent 27-day solar rotational effects, emphasizing the importance of heliospheric plasma dynamics in governing cosmic ray transport.

Significance for Cosmic Ray Modulation

The heliosphere is not merely a passive region through which cosmic rays travel but rather a dynamic plasma laboratory that continuously modifies energetic particle transport through diffusion, convection, magnetic drifts, and energy losses. Variations in solar wind properties, heliospheric magnetic field geometry, magnetic turbulence, and heliospheric boundaries determine the temporal evolution of galactic cosmic ray intensities observed near Earth. Consequently, accurate prediction of cosmic ray modulation requires comprehensive knowledge of heliospheric structure, solar wind dynamics, plasma turbulence, and magnetic field evolution over multiple solar cycles. Modern observational missions, numerical magnetohydrodynamic simulations, and artificial intelligence–based forecasting techniques are increasingly integrated to improve understanding of heliospheric transport and operational space weather prediction.

Parker Transport Equation

The Parker Transport Equation (PTE) is the fundamental mathematical framework describing the propagation and modulation of galactic cosmic rays (GCRs) within the heliosphere. First formulated by Eugene N. Parker (1965), the equation integrates the four principal physical processes governing cosmic ray transport: spatial diffusion, solar wind convection, gradient and curvature drifts, and adiabatic energy changes. These mechanisms collectively determine the temporal, spatial, and energy-dependent variations of cosmic ray intensity observed throughout the heliosphere and form the theoretical basis of virtually all modern heliospheric modulation models (Parker, 1965; Jokipii, 1971; Potgieter, 2013). Since its introduction, the Parker Transport Equation has become the cornerstone of space weather research, heliophysics, and cosmic ray modulation studies.

Unlike neutral radiation, cosmic rays are charged particles that interact strongly with the heliospheric magnetic field (HMF) embedded in the expanding solar wind. As galactic cosmic rays enter the heliosphere from the local interstellar medium (LISM), they experience repeated scattering by magnetic irregularities, outward convection by the solar wind, systematic drift motions caused by magnetic field gradients and curvature, and gradual energy losses due to the adiabatic expansion of the solar wind plasma. The combined influence of these transport processes causes the observed cosmic ray spectrum near Earth to differ significantly from the local interstellar spectrum (LIS), particularly for particles with rigidities below several tens of GV (Potgieter, 2013).

Mathematical Formulation of the Parker Transport Equation

The Parker Transport Equation is expressed as

∂f/∂t = ∇•(K•∇f) − (V_sw + v_d)•∇f + (1/3)(∇•V_sw) ∂f/∂ln p + Q

Where:

f: cosmic-ray distribution function

t: time

K: diffusion tensor

V_sw: solar wind velocity

v_d: particle drift velocity

p: particle momentum

Q: source term

This equation describes cosmic-ray transport through diffusion, solar-wind convection, particle drifts, adiabatic energy changes, and source terms in the heliosphere.where (J(E)) represents the differential cosmic ray intensity (Parker, 1965; Gleeson & Axford, 1968).

The Parker Transport Equation is a multidimensional convection–diffusion equation that simultaneously accounts for spatial transport, energy evolution, and temporal variability. Analytical solutions are possible only under highly simplified assumptions; therefore, realistic heliospheric modulation studies rely on advanced numerical techniques.

Spatial Diffusion

Spatial diffusion represents the random scattering of charged particles by magnetic turbulence embedded within the heliospheric magnetic field. Because the solar wind contains magnetic fluctuations over a broad range of scales, cosmic rays undergo repeated pitch-angle scattering that randomizes their trajectories. Diffusion allows particles to propagate inward against the outward-flowing solar wind and therefore constitutes the primary mechanism through which galactic cosmic rays penetrate the heliosphere (Jokipii, 1966; Potgieter, 2013).

The diffusion tensor is anisotropic,

| K∥ 0 0 |

K = | 0 K⊥r 0 |

| 0 0 K⊥θ |

Where:

• K∥: Diffusion coefficient parallel to the magnetic field

• K⊥r: Perpendicular diffusion coefficient in the radial direction

• K⊥θ: Perpendicular diffusion coefficient in the latitudinal direction

This tensor represents anisotropic diffusion of cosmic rays in the heliosphere, where diffusion is generally greater along magnetic field lines than across them.

Since particle transport occurs much more efficiently along magnetic field lines than across them,

K∥ ≫ K⊥

The diffusion coefficient depends upon rigidity, particle velocity, magnetic turbulence, and heliocentric distance. A commonly adopted parameterization is

K∥ = K₀ β (P/P₀)^a (B₀/B)

where (P) denotes particle rigidity, (B) the heliospheric magnetic field strength, and (\beta=v/c) (Burger et al., 2000).

Solar Wind Convection

The outward expansion of the solar wind continuously transports energetic particles away from the Sun. This convective transport is represented by

−Vsw​⋅∇f

The solar wind velocity generally ranges from 300–800 km s^-1, depending upon solar latitude and solar activity (Schwenn, 2006). Convection competes directly with inward diffusion, and the balance between these two processes largely determines the radial gradient of galactic cosmic rays throughout the heliosphere.

During solar maximum, stronger solar wind disturbances increase convective transport, thereby reducing cosmic ray intensities observed near Earth (Potgieter, 2013).

Gradient and Curvature Drifts

Charged particles moving through the large-scale Parker spiral magnetic field experience systematic drift motions arising from magnetic field gradients, field-line curvature, and the heliospheric current sheet (Jokipii et al., 1977).

The drift velocity is approximately

v_d = ∇ × (K_A B/B)

where (K_A) denotes the antisymmetric drift coefficient.

Drift effects produce:

• Gradient drifts
• Curvature drifts
• Current-sheet drifts

These motions depend upon particle charge, rigidity, and the polarity of the heliospheric magnetic field. Consequently, drift effects reverse direction every 11 years, giving rise to the 22-year Hale magnetic cycle observed in cosmic ray modulation (Burger & Hattingh, 1998).

Adiabatic Energy Changes

As cosmic rays propagate through the expanding solar wind, they lose energy because of adiabatic cooling. The corresponding term in the Parker equation is

(1/3)(∇·Vsw)(∂f/∂ln p)

This process decreases particle momentum continuously during outward propagation and has the greatest influence on low-energy cosmic rays below several GeV (Parker, 1965). Adiabatic cooling explains why the observed spectrum near Earth is softer than the local interstellar spectrum.

Boundary Conditions

Numerical solutions of the Parker equation require physically meaningful boundary conditions.

The outer boundary is generally placed near the heliopause (120–130 AU), where the Local Interstellar Spectrum (LIS) is prescribed from Voyager observations (Stone et al., 2013).

The inner boundary is usually located between 0.1 and 0.3 AU, depending on the numerical model, because transport processes close to the solar corona require additional physical descriptions beyond classical modulation theory (Potgieter, 2013).

Numerical Solution Techniques

Because the Parker Transport Equation is nonlinear and multidimensional, several numerical techniques have been developed.

These include

• Finite Difference Methods (FDM)
• Finite Volume Methods (FVM)
• Stochastic Differential Equation (SDE) methods
• Monte Carlo simulations
• Three-dimensional Magnetohydrodynamic (MHD) coupling

Modern heliospheric models such as HelMod, WSA–ENLIL, EUHFORIA, and CRONOS integrate realistic solar wind conditions with numerical solutions of the Parker equation, enabling accurate simulation of cosmic ray modulation over complete solar cycles (Bobik et al., 2012; Pei et al., 2010).

Artificial Intelligence and Hybrid Modeling

Recent advances in artificial intelligence have complemented classical Parker equation modeling. Machine learning algorithms including Long Short-Term Memory (LSTM) networks, Convolutional Neural Networks (CNNs), Transformer architectures, Random Forests, XGBoost, and Physics-Informed Neural Networks (PINNs) are increasingly employed to estimate diffusion coefficients, reconstruct heliospheric parameters, forecast neutron monitor counts, and predict Forbush decreases (Camporeale, 2019).

Unlike purely empirical methods, hybrid AI–physics frameworks preserve the governing transport physics while improving computational efficiency and predictive accuracy. These approaches are expected to play a major role in future operational space weather forecasting.

Several investigations by Agarwal and Mishra (2013a, 2013b, 2023) examined the modulation of galactic cosmic rays during Solar Cycles 20–24 and demonstrated that cosmic ray intensity variations are strongly correlated with solar wind velocity, heliospheric magnetic field strength, interplanetary parameters, and rigidity-dependent transport processes. Their work provides important observational validation for the physical mechanisms represented in the Parker Transport Equation and highlights the role of heliospheric conditions in long-term cosmic ray modulation.

The Parker Transport Equation remains the most comprehensive theoretical description of cosmic ray propagation within the heliosphere. By simultaneously incorporating spatial diffusion, solar wind convection, gradient and curvature drifts, adiabatic energy losses, and heliospheric boundary conditions, it successfully explains the observed temporal, spatial, and energy-dependent modulation of galactic cosmic rays. Continuous improvements in spacecraft observations, numerical transport models, high-performance computing, and artificial intelligence are transforming the Parker equation into a powerful predictive framework for heliophysics and space weather research.

Physical Modulation Processes

The modulation of galactic cosmic rays (GCRs) within the heliosphere is governed by a complex interplay of plasma and magnetic field processes that modify the intensity, energy spectrum, spatial distribution, and temporal variability of energetic charged particles before they reach Earth. Although the Parker Transport Equation provides the mathematical framework for describing cosmic ray transport, the underlying physical mechanisms responsible for modulation are determined by the properties of the heliospheric magnetic field (HMF), solar wind turbulence, particle rigidity, magnetic polarity, and large-scale heliospheric structures. These processes operate over multiple temporal scales, ranging from transient disturbances lasting hours to long-term variations associated with the 11-year solar cycle and the 22-year Hale magnetic cycle (Parker, 1965; Potgieter, 2013). Understanding these physical modulation mechanisms is essential for accurate prediction of cosmic ray intensities, radiation exposure in space, and operational space weather forecasting.

Anisotropic Diffusion

Spatial diffusion is the dominant transport mechanism allowing galactic cosmic rays to penetrate the heliosphere against the outward-flowing solar wind. Unlike classical isotropic diffusion, cosmic ray transport in the heliosphere is anisotropic because charged particles move preferentially along the heliospheric magnetic field rather than across it (Jokipii, 1966). The Parker spiral magnetic field provides an ordered large-scale geometry that causes the diffusion coefficient parallel to the magnetic field K∥ to exceed the perpendicular diffusion coefficients ((K_{\perp})) by approximately one order of magnitude or more.

The diffusion tensor is expressed as

[\mathbf{K}=\begin{bmatrix}K_{\parallel} & 0 & 0\0 & K_{\perp r} & 0\0 & 0 & K_{\perp\theta}

\end{bmatrix},]

where (K_{\parallel}) represents diffusion along magnetic field lines, while (K_{\perp r}) and (K_{\perp\theta}) describe radial and latitudinal perpendicular diffusion. Typical observations indicate

[K_{\parallel} \approx 5-20 \times K_{\perp}.]

Anisotropic diffusion depends strongly upon magnetic turbulence, heliocentric distance, particle rigidity, and solar activity. During solar minimum, reduced turbulence permits larger diffusion coefficients, allowing greater penetration of galactic cosmic rays into the inner heliosphere. Conversely, enhanced turbulence during solar maximum decreases diffusion coefficients and produces stronger modulation (Potgieter, 2013; Burger et al., 2000).

Magnetic Turbulence

The solar wind is inherently turbulent over an enormous range of spatial scales, from several solar radii to hundreds of astronomical units. Turbulence consists primarily of Alfvén waves, magnetohydrodynamic (MHD) fluctuations, stream interaction regions, coronal mass ejections (CMEs), and nonlinear plasma instabilities that continuously modify the heliospheric magnetic field (Matthaeus et al., 2003).

Magnetic turbulence influences cosmic ray transport by scattering charged particles through resonant interactions between particle gyro-motion and magnetic field fluctuations. The efficiency of scattering depends upon the relationship between particle gyroradius and the wavelength of magnetic turbulence. Consequently, diffusion coefficients become functions of both particle rigidity and the turbulence spectrum.

The solar wind turbulence spectrum approximately follows a Kolmogorov power law

[P(k)\propto k^{-5/3},]

where (k) denotes the wave number. Some heliospheric regions also exhibit an Iroshnikov–Kraichnan spectrum

[P(k)\propto k^{-3/2},]

particularly in magnetically dominated plasma environments (Matthaeus et al., 2003).

Increasing turbulence reduces particle mean free paths and enhances scattering, thereby decreasing cosmic ray intensities observed near Earth. During solar maximum, frequent CMEs and interplanetary shocks substantially increase turbulence levels, leading to stronger cosmic ray modulation and more frequent Forbush decreases.

Rigidity Dependence

One of the defining characteristics of heliospheric modulation is its dependence on particle rigidity, defined as

[R=\frac{pc}{Ze},]

where

• (p) = particle momentum,
• (Z) = atomic charge,
• (e) = elementary charge.

Rigidity determines the ability of charged particles to resist magnetic deflection. Low-rigidity particles are scattered efficiently by heliospheric magnetic irregularities and therefore undergo strong modulation, whereas high-rigidity particles propagate more freely through the heliosphere.

Observations indicate that diffusion coefficients approximately follow

[K\propto R^{a},]

where the spectral index (a) typically ranges between 0.3 and 1.2, depending upon turbulence conditions and solar activity (Burger et al., 2000).

This rigidity dependence explains why cosmic ray modulation primarily affects particles below approximately 20–30 GV, while ultra-high-energy cosmic rays experience negligible heliospheric influence.

Several observational investigations by Agarwal and Mishra (2013a, 2013b) demonstrated that cosmic ray modulation varies significantly with geomagnetic cutoff rigidity, confirming that rigidity-dependent transport is one of the principal characteristics of heliospheric modulation.

Gradient and Curvature Drifts

Charged particles propagating through the large-scale Parker spiral magnetic field experience systematic drift motions caused by magnetic field gradients and field-line curvature (Jokipii et al., 1977). These drift motions represent one of the four principal transport mechanisms incorporated into the Parker Transport Equation.

The total drift velocity is given approximately by

[\mathbf{v}_d\nabla\times\left(K_A\frac{\mathbf{B}}{B}\right),]

where (K_A) denotes the antisymmetric diffusion coefficient.

Three principal drift processes occur within the heliosphere:

• Gradient drift resulting from spatial variations in magnetic field strength.
• Curvature drift caused by the curvature of Parker spiral field lines.
• Heliospheric current sheet (HCS) drift occurring near the magnetic neutral sheet.

Drift motions strongly influence cosmic ray latitudinal gradients, radial intensity distributions, and solar-cycle-dependent modulation patterns.

Charge-Sign Effects

Because drift motions depend upon particle charge, positively and negatively charged particles follow different trajectories through the heliosphere. This phenomenon, known as the charge-sign effect, produces measurable differences between proton, electron, positron, and antiproton modulation (Potgieter, 2013).

During A > 0 polarity cycles, when the heliospheric magnetic field points outward in the northern hemisphere, positively charged particles drift inward primarily through the polar regions and outward along the heliospheric current sheet.

During A < 0 cycles, positively charged particles enter preferentially through the equatorial current sheet, while negatively charged particles drift inward through the polar regions.

Spacecraft observations from AMS-02, PAMELA, and Ulysses have confirmed these charge-dependent modulation effects through measurements of electron–positron and proton–antiproton intensity differences (Aguilar et al., 2015).

Charge-sign dependence represents one of the strongest observational validations of drift theory.

The 22-Year Hale Magnetic Cycle

The Sun undergoes an approximately 11-year sunspot cycle, during which solar activity increases and decreases periodically. However, the global solar magnetic field reverses polarity at every solar maximum, producing a complete magnetic cycle lasting approximately 22 years, known as the Hale Magnetic Cycle.

Because particle drift directions depend upon magnetic polarity, galactic cosmic ray modulation exhibits a corresponding 22-year periodicity (Jokipii et al., 1977).

The two magnetic polarity states are:

• A > 0 Cycle

o Magnetic field directed outward in the northern hemisphere.

o Positively charged particles enter mainly through polar regions.

• A < 0 Cycle

o Magnetic field directed inward in the northern hemisphere.

o Positively charged particles enter primarily along the heliospheric current sheet.

This polarity reversal produces observable differences in:

• Cosmic ray intensity profiles,

• Recovery following solar maxima,

• Latitudinal gradients,

• Drift velocities,

• Charge-sign modulation.

The Hale cycle is therefore one of the defining characteristics of heliospheric cosmic ray modulation and cannot be explained by diffusion and convection alone (Potgieter, 2013).

Solar Cycle Dependence

Cosmic ray modulation varies systematically throughout the approximately 11-year solar cycle.

During solar minimum:

• Heliospheric magnetic field is relatively weak.
• Magnetic turbulence is low.
• Diffusion coefficients are large.
• Drift effects dominate.
• Cosmic ray intensity reaches maximum values.

During solar maximum:

• Solar wind disturbances increase dramatically.
• Coronal mass ejections become frequent.
• Turbulence intensifies.
• Diffusion coefficients decrease.
• Drift effects are suppressed.
• Cosmic ray intensity reaches minimum levels.

This inverse relationship between sunspot number and galactic cosmic ray intensity has been consistently observed by neutron monitor networks for more than seven decades (Simpson, 2000).

Combined Physical Picture

Heliospheric modulation results from the simultaneous action of multiple transport mechanisms rather than any single physical process. Anisotropic diffusion governs large-scale particle transport, magnetic turbulence determines scattering efficiency, rigidity controls sensitivity to magnetic fluctuations, gradient and curvature drifts produce charge-dependent transport, and the Hale magnetic cycle introduces long-term polarity-dependent modulation. Together, these mechanisms generate the observed temporal, spatial, and spectral variations of galactic cosmic rays measured by ground-based neutron monitors and spacecraft including Voyager, ACE, Ulysses, PAMELA, AMS-02, Parker Solar Probe, Solar Orbiter, and Aditya-L1.

Recent observational analyses by Agarwal and Mishra (2013a, 2013b, 2023) demonstrated that rigidity dependence, solar wind velocity, interplanetary magnetic field strength, and recurrent 27-day heliospheric variations collectively regulate long-term galactic cosmic ray modulation over multiple solar cycles. Their findings support modern transport theory and highlight the importance of combining observational data with numerical and artificial intelligence–based modeling approaches.

Solar Cycles 20–25

The long-term modulation of galactic cosmic rays (GCRs) is closely linked to the approximately 11-year solar activity cycle and the associated 22-year Hale magnetic cycle. Variations in solar magnetic activity alter the heliospheric magnetic field (HMF), solar wind velocity, magnetic turbulence, coronal mass ejection (CME) occurrence, and heliospheric current sheet (HCS) geometry, thereby changing the transport conditions experienced by cosmic rays propagating through the heliosphere (Parker, 1965; Potgieter, 2013). Consequently, the intensity of galactic cosmic rays observed near Earth exhibits a well-established inverse correlation with solar activity, reaching maximum values during solar minima and minimum values during solar maxima (Simpson, 2000).

Continuous observations from global neutron monitor networks, together with spacecraft such as IMP, Helios, Voyager, Ulysses, ACE, SOHO, STEREO, PAMELA, AMS-02, Parker Solar Probe, Solar Orbiter, and Aditya-L1, have enabled detailed investigations of cosmic ray modulation across Solar Cycles 20–25. These observations provide valuable insight into the long-term evolution of heliospheric transport processes and validate modern numerical modulation models.

Solar Cycle 20 (1964–1976)

Solar Cycle 20 marked the beginning of the modern era of continuous cosmic ray observations. This cycle exhibited moderate solar activity compared with later cycles and coincided with significant advances in ground-based neutron monitor observations and early space missions.

During this period:

• Sunspot numbers increased from solar minimum to a maximum around 1968–1969.
• Solar wind measurements from early spacecraft confirmed Parker's solar wind theory.
• The anti-correlation between sunspot number and neutron monitor intensity became firmly established.
• Initial theoretical models describing cosmic ray diffusion and convection were developed.
• Parker's Transport Equation became the standard theoretical framework for heliospheric modulation (Parker, 1965).

Long-term neutron monitor observations during Cycle 20 demonstrated that increasing solar activity was accompanied by reduced galactic cosmic ray intensity due to enhanced heliospheric magnetic turbulence and stronger solar wind convection.

Solar Cycle 21 (1976–1986)

Solar Cycle 21 was characterized by relatively high solar activity, frequent solar flares, and numerous coronal mass ejections. Space missions including Helios, IMP, and ISEE provided improved measurements of solar wind plasma and energetic particles.

Important developments included:

• Observation of increased CME occurrence.
• Identification of recurrent 27-day cosmic ray variations associated with corotating interaction regions (CIRs).
• Improved understanding of heliospheric magnetic field structure.
• Development of drift theory explaining charge-sign dependence of cosmic ray modulation (Jokipii et al., 1977).

The high solar activity produced pronounced decreases in neutron monitor counts and numerous Forbush decreases, confirming the influence of transient heliospheric disturbances on cosmic ray transport.

Solar Cycle 22 (1986–1996)

Solar Cycle 22 exhibited stronger solar magnetic fields and frequent polarity reversals. This cycle coincided with the launch of the Ulysses mission, which provided the first observations of the heliosphere at high solar latitudes.

Major advances included:

• Three-dimensional mapping of the heliospheric magnetic field.
• Measurement of latitudinal cosmic ray gradients.
• Confirmation of drift-dominated modulation during solar minimum.
• Improved determination of diffusion coefficients.

Observations during Cycle 22 provided strong evidence supporting the role of gradient and curvature drifts, particularly during periods of low solar activity (Burger & Hattingh, 1998).

Solar Cycle 23 (1996–2008)

Solar Cycle 23 represents one of the most intensively studied solar cycles because of the availability of high-quality spacecraft observations from ACE, SOHO, Ulysses, and Voyager.

This cycle was characterized by:

• Strong solar maximum around 2000–2002.
• Frequent X-class solar flares.
• Numerous halo CMEs.
• Exceptional interplanetary disturbances.
• Extensive neutron monitor observations.

The cycle also produced some of the largest Forbush decreases ever recorded, demonstrating the strong influence of CME-driven shocks on galactic cosmic ray transport (Cane, 2000).

Near the end of Cycle 23, the Sun entered an unusually prolonged and quiet solar minimum (2008–2009). During this period:

• Solar wind pressure reached historically low values.
• Heliospheric magnetic field strength decreased significantly.
• Cosmic ray intensities measured by neutron monitors reached record highs.

These observations challenged existing modulation models and stimulated extensive improvements in heliospheric transport theory (Potgieter, 2013).

Solar Cycle 24 (2008–2019)

Solar Cycle 24 was one of the weakest solar cycles observed during the space age, exhibiting considerably lower sunspot numbers and reduced solar magnetic activity compared with Cycles 21–23.

Major characteristics included:

• Weak solar maximum around 2014.
• Reduced heliospheric magnetic field strength.
• Lower solar wind dynamic pressure.
• Reduced CME frequency compared with previous cycles.
• Record-breaking galactic cosmic ray intensities.

Several landmark observations occurred during this cycle.

Voyager and the Heliopause

In 2012, Voyager 1 crossed the heliopause and directly measured the Local Interstellar Spectrum (LIS) for the first time (Stone et al., 2013). These measurements provided the essential outer boundary condition for modern modulation models.

AMS-02 and PAMELA

The AMS-02 detector aboard the International Space Station and the PAMELA satellite measured proton, helium, electron, positron, and antiproton spectra with unprecedented precision, confirming:

• Rigidity-dependent modulation,
• Charge-sign effects,
• Drift-dominated transport during solar minimum.

Several investigations by Agarwal and Mishra (2013a, 2013b, 2023) analyzed long-term galactic cosmic ray modulation during Solar Cycles 20–24 using neutron monitor observations and interplanetary parameters. Their studies demonstrated significant correlations between cosmic ray intensity and:

• Solar wind velocity,
• Interplanetary magnetic field strength,
• Sunspot number,
• Solar magnetic polarity,
• Recurrent 27-day heliospheric variations.

These investigations reinforced the importance of heliospheric transport processes in explaining long-term cosmic ray modulation.

Solar Cycle 25 (2019–Present)

Solar Cycle 25 represents the current phase of solar activity and is being observed using an unprecedented combination of spacecraft and ground-based facilities.

Major observational platforms include:

• Parker Solar Probe
• Solar Orbiter
• Aditya-L1
• AMS-02
• Global neutron monitor networks
• Solar Dynamics Observatory (SDO)
• SOHO
• DSCOVR

Unlike previous cycles, Solar Cycle 25 benefits from continuous observations of:

• Near-Sun solar wind,
• Magnetic reconnection,
• Coronal mass ejections,
• Solar energetic particles,
• Heliospheric magnetic field evolution.

Recent observations indicate that Solar Cycle 25 is more active than initially predicted, with increasing sunspot numbers, frequent X-class flares, and numerous halo CMEs. Correspondingly, neutron monitor observations show decreasing galactic cosmic ray intensity, consistent with enhanced heliospheric modulation.

The availability of real-time observations from Parker Solar Probe, Solar Orbiter, and Aditya-L1 provides unprecedented opportunities to validate numerical transport models and artificial intelligence–based forecasting systems.

Sunspot Variability and Cosmic Ray Intensity

One of the most robust relationships in heliophysics is the inverse correlation between sunspot number and galactic cosmic ray intensity.

During solar minimum:

• Low sunspot number,
• Weak heliospheric magnetic field,
• Reduced turbulence,
• Larger diffusion coefficients,
• Maximum neutron monitor counts.

During solar maximum:

• High sunspot number,
• Frequent CMEs,
• Enhanced turbulence,
• Stronger convection,
• Reduced neutron monitor counts.

This anti-correlation has remained remarkably consistent throughout Solar Cycles 20–25 (Simpson, 2000; Potgieter, 2013).

Heliospheric Magnetic Polarity Reversals

Every solar maximum is accompanied by reversal of the Sun's global magnetic field.

These polarity reversals produce:

• Changes in gradient and curvature drifts,
• Altered particle entry paths,
• Different modulation for positive and negative particles,
• The 22-year Hale magnetic cycle.

The alternating A > 0 and A < 0 polarity epochs produce distinct cosmic ray intensity profiles, particularly during solar minima (Jokipii et al., 1977).

Coronal Mass Ejections and Forbush Decreases

Coronal Mass Ejections represent one of the most important transient modulation mechanisms.

Fast CMEs generate:

• Interplanetary shocks,
• Magnetic clouds,
• Enhanced turbulence,
• Temporary exclusion of galactic cosmic rays.

These disturbances produce Forbush decreases, characterized by sudden reductions in neutron monitor counts followed by gradual recovery over several days (Forbush, 1954).

The frequency and magnitude of Forbush decreases increase markedly during solar maximum because of enhanced CME activity.

Neutron Monitor Observations

Ground-based neutron monitor stations provide the longest continuous record of cosmic ray intensity available.

Networks including:

• Oulu,
• Newark,
• Moscow,
• McMurdo,
• Climax,
• South Pole,

have operated for several decades and have documented:

• Solar cycle modulation,
• 27-day recurrent variations,
• Forbush decreases,
• Ground Level Enhancements (GLEs),
• Long-term heliospheric variability.

These observations remain indispensable for validating Parker transport models and modern AI-based forecasting systems.

Overall Perspective

Analysis of Solar Cycles 20–25 demonstrates that galactic cosmic ray modulation is controlled by the combined influence of sunspot activity, heliospheric magnetic polarity, solar wind turbulence, coronal mass ejections, and long-term evolution of the heliosphere. Improvements in spacecraft observations—from early missions during Cycle 20 to Parker Solar Probe, Solar Orbiter, and Aditya-L1 during Cycle 25—have transformed understanding of heliospheric transport and substantially refined the Parker Transport Equation. Together with long-term neutron monitor observations and the contributions of Agarwal and Mishra on rigidity dependence and interplanetary modulation, these datasets provide a comprehensive observational foundation for contemporary cosmic ray modulation theory and future artificial intelligence–driven space weather forecasting.

Observational Data

Observational measurements form the foundation of modern cosmic ray modulation studies by providing continuous records of galactic cosmic ray (GCR) intensity, solar wind properties, heliospheric magnetic fields, and transient solar phenomena across multiple solar cycles. Since cosmic ray modulation is governed by dynamic interactions between energetic charged particles and the heliosphere, observations from both spacecraft missions and ground-based neutron monitor networks are essential for validating theoretical models, constraining numerical simulations, and improving space weather forecasting. Over the past five decades, a combination of interplanetary probes, near-Earth satellites, and global neutron monitor stations has enabled unprecedented investigations of cosmic ray transport from the solar corona to the local interstellar medium (Potgieter, 2013; Stone et al., 2013).

The integration of measurements from Voyager, ACE, Ulysses, SOHO, STEREO, PAMELA, AMS-02, Parker Solar Probe, Solar Orbiter, Aditya-L1, the Oulu Neutron Monitor, and the Neutron Monitor Database (NMDB) has transformed the understanding of heliospheric modulation by providing comprehensive observations spanning nearly the entire heliosphere. These datasets collectively measure solar wind plasma, magnetic fields, energetic particles, solar energetic particle (SEP) events, coronal mass ejections (CMEs), heliospheric current sheet dynamics, and long-term galactic cosmic ray intensity variations.

Voyager Missions

The Voyager 1 and Voyager 2 spacecraft, launched in 1977, have provided the longest continuous in situ observations of energetic particles within the heliosphere. Their primary contribution to cosmic ray research has been the direct exploration of the outer heliosphere, including the termination shock, heliosheath, and heliopause.

Major scientific contributions include:

• Discovery of the heliospheric termination shock.
• Measurement of anomalous cosmic rays (ACRs).
• Direct observations of heliosheath plasma.
• First measurement of the Local Interstellar Spectrum (LIS).
• Determination of heliopause boundary conditions.

In August 2012, Voyager 1 crossed the heliopause, becoming the first spacecraft to enter interstellar space. Measurements beyond the heliopause revealed galactic cosmic ray intensities largely free from heliospheric modulation, providing the benchmark Local Interstellar Spectrum used as the outer boundary condition in Parker transport models (Stone et al., 2013).

Advanced Composition Explorer (ACE)

Launched in 1997, the Advanced Composition Explorer (ACE) occupies the Sun–Earth L1 Lagrange point and continuously monitors the upstream solar wind before it reaches Earth.

ACE measures:

• Solar wind velocity.
• Proton density.
• Interplanetary magnetic field.
• Solar energetic particles.
• Galactic cosmic rays.
• Heavy-ion composition.

The spacecraft provides essential real-time solar wind parameters that are widely used in numerical cosmic ray modulation models and operational space weather forecasting (Stone et al., 1998).

Ulysses

The Ulysses mission (1990–2009) revolutionized heliophysics by becoming the first spacecraft to observe the heliosphere over the Sun's poles.

Major contributions include:

• Three-dimensional mapping of the heliosphere.
• Measurement of high-latitude solar wind.
• Determination of latitudinal cosmic ray gradients.
• Confirmation of particle drifts theory.
• Solar-cycle dependence of cosmic ray transport.

Ulysses demonstrated that galactic cosmic ray intensities vary significantly with heliographic latitude, confirming predictions of gradient and curvature drift models (Heber & Potgieter, 2006).

Solar and Heliospheric Observatory (SOHO)

The Solar and Heliospheric Observatory (SOHO), launched jointly by ESA and NASA in 1995, has continuously observed the Sun from the Earth–Sun L1 point.

SOHO contributes to cosmic ray modulation research by monitoring:

• Coronal mass ejections (CMEs).
• Solar flares.
• Coronal magnetic structures.
• Solar wind source regions.
• Helioseismology.

The LASCO coronagraph aboard SOHO has produced one of the most comprehensive CME catalogs available, enabling detailed studies of CME-driven Forbush decreases and solar energetic particle events (Domingo et al., 1995).

STEREO

NASA's Solar TErrestrial RElations Observatory (STEREO) mission consists of two spacecraft—STEREO-A and STEREO-B—launched in 2006 to provide stereoscopic observations of the Sun.

Its principal contributions include:

• Three-dimensional CME reconstruction.
• CME propagation tracking.
• Solar energetic particle observations.
• Heliospheric imaging.
• Magnetic cloud evolution.

STEREO observations have significantly improved understanding of the relationship between CMEs, interplanetary shocks, and cosmic ray modulation.

PAMELA

The Payload for Antimatter Matter Exploration and Light-nuclei Astrophysics (PAMELA) satellite operated between 2006 and 2016, providing high-precision measurements of charged cosmic rays.

PAMELA measured:

• Protons.
• Helium nuclei.
• Electrons.
• Positrons.
• Antiprotons.
• Light nuclei.

The mission confirmed rigidity-dependent modulation and provided valuable evidence for charge-sign effects predicted by heliospheric drift theory (Adriani et al., 2013).

AMS-02

The Alpha Magnetic Spectrometer (AMS-02) has operated aboard the International Space Station since 2011, providing the most precise measurements of charged cosmic rays ever obtained.

AMS-02 measures:

• Proton spectra.
• Helium nuclei.
• Electrons.
• Positrons.
• Antiprotons.
• Heavy nuclei.
• Isotopic abundances.

Its exceptional precision has enabled detailed investigations of:

• Rigidity dependence.
• Solar-cycle modulation.
• Charge-sign modulation.
• Cosmic ray composition.
• Time-dependent particle spectra.

AMS-02 data constitute one of the primary datasets for validating contemporary heliospheric transport models (Aguilar et al., 2015).

Parker Solar Probe

Launched in 2018, the Parker Solar Probe (PSP) is the first spacecraft designed to fly directly through the solar corona.

Its scientific objectives include:

• Solar wind acceleration.
• Magnetic reconnection.
• Coronal heating.
• Near-Sun turbulence.
• Solar energetic particle acceleration.

PSP has measured magnetic switchbacks, Alfvénic turbulence, and near-SSun plasma conditions that directly influence cosmic ray transport coefficients (Fox et al., 2016).

Solar Orbiter

The Solar Orbiter mission, launched in 2020 by ESA and NASA, combines remote sensing with in situ measurements.

The spacecraft investigates:

• Solar magnetic fields.
• Polar solar regions.
• Solar wind formation.
• Coronal dynamics.
• CME initiation.
• Solar energetic particles.

Solar Orbiter observations complement Parker Solar Probe measurements by linking near-Sun plasma processes with heliospheric evolution (Müller et al., 2020).

Aditya-L1

India's Aditya-L1 mission, launched by ISRO in 2023, is positioned near the Sun–Earth L1 point and represents India's first dedicated solar observatory.

Its scientific payloads monitor:

• Solar corona.
• Coronal magnetic fields.
• Solar wind plasma.
• Energetic particles.
• Solar flares.
• Coronal mass ejections.
• Interplanetary magnetic field.

Aditya-L1 provides continuous observations of the upstream solar wind and solar activity, contributing significantly to future studies of cosmic ray modulation and operational space weather forecasting.

Oulu Neutron Monitor

The Oulu Neutron Monitor, located in Finland, is one of the world's most widely used ground-based cosmic ray observatories.

The station continuously records secondary neutrons generated when primary galactic cosmic rays interact with Earth's atmosphere.

Its measurements provide:

• Long-term cosmic ray intensity.
• Solar cycle modulation.
• Forbush decreases.
• Ground Level Enhancements (GLEs).
• Daily cosmic ray variations.
• 27-day recurrent modulation.

Because of its high stability and nearly continuous operation, Oulu data are extensively employed in heliospheric modulation studies and numerical model validation (Usoskin et al., 2015).

Neutron Monitor Database (NMDB)

The Neutron Monitor Database (NMDB) integrates real-time observations from dozens of neutron monitor stations distributed worldwide.

NMDB provides:

• Minute-resolution cosmic ray intensity.
• Historical neutron monitors records.
• Global station comparisons.
• Ground Level Enhancement catalogs.
• Forbush decrease observations.
• Solar cycle analyses.

The global coverage of NMDB enables researchers to investigate cosmic ray anisotropies, geomagnetic cutoff rigidity effects, and worldwide modulation patterns. It has become the standard repository for neutron monitor observations used in space weather forecasting and heliospheric research.

Multi-Mission Synergy

Modern cosmic ray modulation research increasingly relies on combining measurements from multiple observational platforms.

A typical integrated analysis includes:

Long-term observational analyses by Agarwal and Mishra (2013a, 2013b, 2023) have made important contributions to understanding galactic cosmic ray modulation across Solar Cycles 20–24. Using neutron monitor observations together with interplanetary parameters, their studies examined rigidity-dependent modulation, 27-day recurrent variations, and the influence of solar wind velocity and interplanetary magnetic field strength on cosmic ray intensity. These investigations complement spacecraft observations by providing long-term statistical validation of heliospheric transport models and demonstrating the close relationship between solar activity and galactic cosmic ray variability.

Future Observational Directions

The integration of spacecraft observations, neutron monitor networks, heliospheric numerical simulations, and artificial intelligence is ushering in a new era of cosmic ray research. Real-time measurements from Parker Solar Probe, Solar Orbiter, Aditya-L1, AMS-02, ACE, and the NMDB, combined with advanced machine learning techniques, are enabling increasingly accurate predictions of cosmic ray modulation, solar energetic particle events, and radiation hazards for future human exploration of the Moon, Mars, and deep-space environments.

Numerical Models and Artificial Intelligence

Numerical modeling has become indispensable for understanding the transport and modulation of galactic cosmic rays (GCRs) within the heliosphere. Although the Parker Transport Equation (PTE) provides the theoretical framework describing cosmic ray propagation, its multidimensional, nonlinear, and time-dependent nature precludes analytical solutions under realistic heliospheric conditions. Consequently, a variety of numerical techniques have been developed to solve the Parker equation while incorporating realistic solar wind dynamics, heliospheric magnetic field (HMF) geometry, turbulence, particle drifts, and transient solar events (Parker, 1965; Potgieter, 2013). In recent years, artificial intelligence (AI) and machine learning (ML) have emerged as powerful complementary tools, enabling rapid prediction of cosmic ray intensity, optimization of transport parameters, and real-time space weather forecasting. Hybrid AI–physics models now represent one of the most promising directions in heliospheric research.

Evolution of Numerical Modeling

The earliest numerical studies of cosmic ray modulation employed simplified one-dimensional finite-difference solutions of Parker's Transport Equation during the late 1960s and early 1970s. As computational resources advanced, models evolved toward two-dimensional and fully three-dimensional simulations incorporating realistic heliospheric geometry, anisotropic diffusion, gradient and curvature drifts, heliospheric current sheet (HCS) dynamics, and solar-cycle-dependent boundary conditions (Jokipii, 1971; Potgieter, 2013).

Modern numerical models integrate observational data from spacecraft such as ACE, SOHO, STEREO, Voyager, Parker Solar Probe, Solar Orbiter, and Aditya-L1, together with neutron monitor observations, to simulate cosmic ray transport across the heliosphere under both steady-state and time-dependent conditions.

Finite-Difference Methods (FDM)

Finite-Difference Methods constitute one of the earliest and most widely used numerical approaches for solving the Parker Transport Equation. In this method, spatial and temporal derivatives are approximated using finite differences on a computational grid.

For a one-dimensional diffusion equation,

[\frac{\partial f}{\partial t}K\frac{\partial^2 f}{\partial r^2},]

the spatial derivative can be approximated as

[\frac{\partial^2 f}{\partial r^2}\approx\frac{f_{i+1}-2f_i+f_{i-1}}{\Delta r^2}.]

Finite-difference methods offer:

• High numerical accuracy.
• Straightforward implementation.
• Efficient treatment of steady-state modulation.
• Good convergence for smooth heliospheric conditions.

However, they become computationally demanding in fully three-dimensional, time-dependent simulations involving complex heliospheric geometries and transient disturbances.

Stochastic Differential Equation (SDE) Methods

Stochastic Differential Equation (SDE) techniques have become one of the most successful methods for solving the Parker Transport Equation. Rather than solving the transport equation directly, the method transforms the deterministic partial differential equation into an equivalent set of stochastic differential equations describing pseudo-particle trajectories (Zhang, 1999).

A simplified stochastic equation is

[dXA(X),dt+B(X),dW,]

where

• (A(X)) represents deterministic transport,
• (B(X)) describes diffusion,
• (dW) denotes Wiener random processes.

Advantages include:

• Naturally handles multidimensional transport.
• Efficient treatment of complex boundaries.
• Easily incorporates particle drifts.
• Well suited for realistic heliospheric geometries.
• Widely adopted in modern modulation studies.

SDE methods form the computational basis of several contemporary cosmic ray transport models.

Monte Carlo Simulations

Monte Carlo methods simulate the trajectories of large ensembles of particles undergoing random scattering within the heliosphere.

Each simulated particle experiences:

• Diffusion,
• Drift,
• Solar wind convection,
• Adiabatic cooling,

until reaching the observation point.

Advantages include:

• Excellent statistical representation.
• Flexible treatment of complex physics.
• Easy inclusion of turbulence.
• Suitable for nonlinear transport.

Limitations include:

• Large computational cost.
• Slow convergence for high-precision simulations.
• Requirement of millions of particle trajectories.

Despite these challenges, Monte Carlo techniques remain valuable for validating deterministic numerical models.

Magnetohydrodynamic (MHD) Models

Magnetohydrodynamic (MHD) models simulate the evolution of heliospheric plasma and magnetic fields by solving the coupled conservation equations of mass, momentum, magnetic flux, and energy.

The ideal MHD equations include:

• Continuity equation,
• Momentum equation,
• Induction equation,
• Energy conservation equation.

MHD models provide:

• Solar wind velocity,
• Plasma density,
• Magnetic field structure,
• Shock evolution,
• CME propagation,
• Heliospheric current sheet geometry.

These plasma parameters serve as inputs for cosmic ray transport simulations, thereby coupling heliospheric plasma dynamics with energetic particle propagation.

HelMod

HelMod is one of the most widely used numerical models specifically designed for galactic cosmic ray modulation.

Key characteristics include:

• Three-dimensional Parker equation solver.
• Drift-dependent transport.
• Time-dependent heliosphere.
• Realistic diffusion coefficients.
• Solar-cycle-dependent parameters.
• Validation using PAMELA and AMS-02 observations.

HelMod successfully reproduces long-term proton, helium, and electron spectra measured during Solar Cycles 23–25 (Bobik et al., 2012).

WSA–ENLIL

The WSA–ENLIL model combines the Wang–Sheeley–Arge (WSA) solar wind model with the ENLIL three-dimensional MHD heliospheric model.

Primary applications include:

• Solar wind prediction.
• CME propagation.
• Interplanetary shock evolution.
• Solar energetic particle forecasting.
• Operational space weather prediction.

Although originally developed for solar wind forecasting, WSA–ENLIL provides realistic heliospheric conditions that significantly improve cosmic ray modulation simulations.

EUHFORIA

The European Heliospheric Forecasting Information Asset (EUHFORIA) is a state-of-the-art three-dimensional MHD model developed for heliospheric forecasting.

Major capabilities include:

• Solar wind simulation.
• CME evolution.
• Interplanetary magnetic field prediction.
• Heliospheric current sheet modeling.
• Solar energetic particle propagation.

EUHFORIA has become an important framework for coupling heliospheric plasma simulations with energetic particle transport and operational space weather services (Pomoell & Poedts, 2018).

Artificial Intelligence Approaches

Artificial intelligence is increasingly transforming cosmic ray modulation research by learning nonlinear relationships between heliospheric conditions and energetic particle observations. Unlike traditional numerical models that explicitly solve the Parker Transport Equation, AI algorithms infer complex patterns directly from observational datasets.

These methods utilize measurements from:

• Voyager,
• ACE,
• SOHO,
• STEREO,
• AMS-02,
• PAMELA,
• Parker Solar Probe,
• Solar Orbiter,
• Aditya-L1,
• Oulu Neutron Monitor,
• NMDB.

AI provides:

• Rapid prediction,
• Parameter optimization,
• Missing-data reconstruction,
• Feature extraction,
• Real-time forecasting.

Long Short-Term Memory (LSTM)

LSTM networks are recurrent neural networks designed to model sequential time-series data.

Applications include:

• Cosmic ray intensity prediction.
• Neutron monitor forecasting.
• Solar wind time-series modeling.
• Forbush decrease prediction.
• Solar cycle forecasting.

Because heliospheric modulation exhibits long temporal dependencies, LSTM models significantly outperform conventional autoregressive approaches.

Convolutional Neural Networks (CNNs)

CNNs automatically extract spatial and temporal features from multidimensional datasets.

Applications include:

• CME classification.
• Solar image analysis.
• Coronal hole detection.
• Solar flare prediction.
• Heliospheric feature extraction.

CNN-derived heliospheric parameters can subsequently be incorporated into modulation models.

Transformer Networks

Transformer architectures employ self-attention mechanisms that capture long-range dependencies without recurrent connections.

Their advantages include:

• Parallel computation.
• Long-sequence learning.
• Multi-variable prediction.
• Improved forecasting accuracy.

Transformers are increasingly applied to:

• Solar activity prediction,
• Space weather forecasting,
• Cosmic ray time-series analysis,
• Multivariate heliospheric datasets.

Physics-Informed Neural Networks (PINNs)

Physics-Informed Neural Networks represent one of the most significant recent developments in scientific machine learning.

Rather than relying solely on observational data, PINNs incorporate governing physical equations directly into the neural network loss function.

For cosmic ray modulation, PINNs enforce:

• Parker Transport Equation,
• Boundary conditions,
• Initial conditions,
• Conservation laws.

Advantages include:

• Physically consistent predictions.
• Reduced data requirements.
• Improved extrapolation.
• Better interpretability than purely data-driven models.

PINNs offer a promising pathway toward hybrid numerical–AI transport models.

XGBoost

Extreme Gradient Boosting (XGBoost) is a decision-tree ensemble algorithm capable of modeling nonlinear relationships efficiently.

Applications include:

• Cosmic ray intensity regression.
• Space weather classification.
• Feature importance analysis.
• Solar flare prediction.
• Geomagnetic storm forecasting.

XGBoost performs particularly well on structured heliospheric datasets with heterogeneous predictors.

Random Forest

Random Forest algorithms combine multiple decision trees to improve predictive robustness.

Applications include:

• Forbush decrease classification.
• Solar wind prediction.
• Variable importance ranking.
• Missing-data estimation.
• Space weather event detection.

Random Forest models provide excellent interpretability and remain valuable baseline methods for heliospheric prediction.

Hybrid Physics–AI Frameworks

The next generation of cosmic ray modeling combines traditional numerical transport models with artificial intelligence.

Typical hybrid workflow:

1. MHD simulation generates heliospheric plasma.
2. Parker Transport Equation models particle propagation.
3. AI estimates diffusion coefficients.
4. Machine learning predicts cosmic ray intensity.
5. Real-time observations continuously update forecasts.

Such hybrid systems preserve the physical rigor of transport theory while exploiting the predictive capabilities of modern machine learning.

Comparative Assessment of Numerical Models

Future Directions

Future progress in cosmic ray modulation research will increasingly rely on the convergence of physics-based numerical models, high-performance computing, and artificial intelligence. Continuous observations from Parker Solar Probe, Solar Orbiter, Aditya-L1, AMS-02, ACE, Voyager, and global neutron monitor networks, combined with hybrid MHD–Parker–AI frameworks, are expected to enable near-real-time forecasting of galactic cosmic ray intensities, solar energetic particle events, and radiation hazards throughout the heliosphere. The integration of Physics-Informed Neural Networks, digital twins of the heliosphere, data assimilation techniques, and explainable AI will likely define the next generation of predictive heliophysics.

Space Weather Applications

Cosmic ray modulation is not merely a theoretical aspect of heliophysics but a phenomenon with profound practical implications for modern technological society, human space exploration, aviation safety, satellite operations, planetary science, and Earth's atmospheric system. Variations in galactic cosmic ray (GCR) intensity, driven by changes in solar activity and heliospheric conditions, influence radiation exposure, electronic system reliability, communication infrastructure, navigation systems, and the radiation environment encountered throughout the Solar System. Accurate prediction of cosmic ray modulation has therefore become a major objective of contemporary space weather research. Advances in observational capabilities, numerical transport models, and artificial intelligence (AI) have significantly improved forecasting of cosmic ray intensity, enabling more effective mitigation of radiation hazards and operational risks (Schrijver et al., 2015; Potgieter, 2013).

Cosmic Rays as a Space Weather Hazard

Space weather encompasses the dynamic conditions in the Sun–Earth system arising from solar activity, including solar flares, coronal mass ejections (CMEs), solar energetic particle (SEP) events, geomagnetic storms, and variations in galactic cosmic ray intensity. Unlike transient solar energetic particles, galactic cosmic rays constitute a continuous background radiation environment whose intensity varies systematically with the solar cycle and rapidly during interplanetary disturbances such as Forbush decreases.

The principal space weather hazards associated with cosmic rays include:

• Enhanced ionizing radiation.
• Single-event effects (SEEs) in electronics.
• Radiation damage to spacecraft materials.
• Increased radiation exposure during aviation.
• Health risks to astronauts.
• Communication system degradation.
• Navigation and positioning errors.

Consequently, accurate modeling of heliospheric modulation has become an essential component of operational space weather forecasting.

Astronaut Radiation Exposure

Radiation exposure remains one of the greatest challenges for long-duration human space exploration. Outside Earth's magnetosphere, astronauts are continuously exposed to:

• Galactic Cosmic Rays (GCRs),
• Solar Energetic Particles (SEPs),
• Secondary radiation generated within spacecraft shielding.

Galactic cosmic rays contain highly energetic heavy ions capable of penetrating conventional shielding materials and producing complex biological damage through dense ionization tracks.

Potential health effects include:

• DNA double-strand breaks.
• Increased cancer risk.
• Central nervous system impairment.
• Cataracts.
• Cardiovascular disease.
• Acute radiation syndrome during extreme SEP events.

Radiation dose equivalent during interplanetary missions can exceed several hundred millisieverts, making cosmic ray forecasting a critical requirement for future crewed missions to the Moon, Mars, and deep-space destinations (Cucinotta et al., 2013).

Modern transport models coupled with spacecraft observations from Parker Solar Probe, Solar Orbiter, ACE, and Aditya-L1 are increasingly employed to estimate astronaut radiation exposure under varying heliospheric conditions.

Aviation Radiation

Commercial aircraft operating on high-latitude and polar routes are exposed to elevated cosmic ray fluxes because Earth's magnetic shielding is weakest near the poles.

Radiation exposure depends upon:

• Flight altitude.
• Geomagnetic latitude.
• Solar activity.
• Flight duration.
• Aircraft trajectory.

During solar minimum, enhanced galactic cosmic ray intensity increases radiation doses received by:

• Airline passengers.
• Pilots.
• Cabin crew.

Conversely, solar energetic particle events can produce temporary radiation spikes during solar maximum.

Cosmic ray monitoring has therefore become increasingly important for:

• Polar aviation.
• Long-duration international flights.
• Supersonic and high-altitude aircraft.
• Future suborbital transportation.

Operational forecasting systems now integrate neutron monitor observations with heliospheric transport models to estimate aviation radiation exposure in real time.

Satellite Systems

Artificial satellites operate continuously within the space radiation environment and are therefore vulnerable to energetic charged particles.

Cosmic rays affect satellites through:

• Surface charging.
• Internal charging.
• Solar cell degradation.
• Detector background noise.
• Material damage.
• Electronic component failures.

High-energy heavy ions can penetrate spacecraft shielding and produce Single Event Effects (SEEs), including:

• Single Event Upsets (SEUs),
• Single Event Latchups (SELs),
• Single Event Burnout (SEB),
• Memory corruption,
• Processor resets.

Accurate prediction of galactic cosmic ray intensity assists spacecraft operators in:

• Radiation risk assessment.
• Mission planning.
• Payload protection.
• Spacecraft anomaly diagnosis.

Modern satellite missions increasingly incorporate onboard radiation monitors calibrated using heliospheric modulation models.

Electronics and Semiconductor Reliability

One of the most significant technological consequences of cosmic rays is their interaction with semiconductor devices.

Energetic particles generate electron–hole pairs within integrated circuits, producing transient or permanent electronic failures.

Affected systems include:

• Computer processors.
• Satellite electronics.
• Aircraft avionics.
• Autonomous vehicles.
• Medical devices.
• High-performance computing systems.
• Data centers.

As semiconductor feature sizes continue to decrease below several nanometers, device sensitivity to ionizing radiation increases substantially.

Consequently, radiation-aware chip design and cosmic ray forecasting have become increasingly important for next-generation electronic systems.

Communication Systems

Space weather significantly influences radio wave propagation throughout Earth's atmosphere.

Cosmic ray modulation contributes indirectly by affecting:

• Atmospheric ionization.
• D-region electron density.
• Radio-wave absorption.
• High-frequency (HF) communication.
• Satellite communication quality.

During intense solar disturbances, increased particle precipitation and geomagnetic activity can disrupt:

• HF radio communication,
• Satellite communication,
• Navigation systems,
• Emergency communication networks.

Combined monitoring of solar energetic particles and galactic cosmic rays improves operational forecasting for communication infrastructure.

Navigation and Positioning Systems

Global Navigation Satellite Systems (GNSS), including GPS, Galileo, GLONASS, and NavIC, rely on accurate propagation of radio signals through Earth's ionosphere.

Space weather disturbances associated with solar activity and energetic particles produce:

• Ionospheric scintillation.
• Signal delays.
• Positioning errors.
• Timing inaccuracies.

Cosmic ray observations complement ionospheric monitoring by providing information regarding heliospheric conditions associated with major space weather events.

Future integrated forecasting systems combine:

• Cosmic ray observations,
• Solar wind measurements,
• Magnetospheric models,
• AI prediction algorithms,

to improve navigation reliability.

Planetary Exploration

Future exploration of the Moon, Mars, and other planetary bodies requires accurate prediction of cosmic ray radiation environments.

Unlike Earth, most planetary bodies possess:

• Weak magnetic shielding,
• Thin or absent atmospheres.

Consequently, astronauts and robotic systems experience substantially higher radiation exposure.

Cosmic ray modulation influences:

• Habitat shielding requirements.
• Rover electronics.
• Scientific instruments.
• Human mission duration.
• Surface radiation forecasting.

Mars missions particularly depend upon accurate heliospheric transport models because radiation exposure represents one of the principal engineering and biomedical challenges for long-duration exploration.

Data from Curiosity's Radiation Assessment Detector (RAD), together with heliospheric observations from ACE, Voyager, Parker Solar Probe, and Solar Orbiter, have significantly improved understanding of interplanetary radiation environments.

Climate and Atmospheric Implications

Galactic cosmic rays contribute to atmospheric ionization through interactions with atmospheric nuclei, producing cascades of secondary particles including neutrons, muons, electrons, and gamma rays.

These interactions influence:

• Atmospheric chemistry.
• Ion production.
• Ozone chemistry.
• Electrical conductivity.
• Cloud microphysics.

One proposed mechanism, often referred to as the ion-aerosol nucleation hypothesis, suggests that enhanced atmospheric ionization may influence aerosol formation and cloud condensation nuclei under specific environmental conditions. This potential link has been investigated through laboratory studies such as the CLOUD experiment at CERN, although the magnitude of its contribution to observed climate variability remains an active area of research. Current assessments indicate that while cosmic rays measurably affect atmospheric ionization, evidence for a dominant influence on recent global climate change is limited compared with well-established anthropogenic greenhouse gas forcing (Kirkby et al., 2016; IPCC, 2021).

Cosmic ray observations nevertheless remain valuable for:

• Atmospheric radiation studies.
• Long-term atmospheric ionization monitoring.
• Paleoclimate reconstructions using cosmogenic isotopes (e.g., ¹⁰Be and ¹⁴C).
• Investigations of solar variability over centennial to millennial timescales.

Operational Space Weather Forecasting

Operational forecasting increasingly combines observations from:

• Voyager,
• ACE,
• SOHO,
• STEREO,
• Parker Solar Probe,
• Solar Orbiter,
• Aditya-L1,
• AMS-02,
• Oulu Neutron Monitor,
• NMDB,

with advanced numerical models including:

• HelMod,
• WSA–ENLIL,
• EUHFORIA,

and artificial intelligence techniques such as:

• Long Short-Term Memory (LSTM),
• Convolutional Neural Networks (CNNs),
• Transformers,
• Physics-Informed Neural Networks (PINNs),
• XGBoost,
• Random Forest.

These integrated systems enable:

• Real-time cosmic ray forecasting.
• Radiation dose estimation.
• Forbush decrease prediction.
• Solar energetic particle forecasting.
• Mission risk assessment.
• Autonomous spacecraft decision support.

The emergence of hybrid physics–AI frameworks is expected to substantially improve forecasting accuracy during Solar Cycle 25 and future solar cycles.

Long-term investigations by Agarwal and Mishra (2013a, 2013b, 2023) demonstrated that variations in galactic cosmic ray intensity are strongly correlated with solar wind velocity, interplanetary magnetic field strength, heliospheric magnetic polarity, and recurrent 27-day solar rotational effects. Their analyses provide important observational evidence supporting the application of heliospheric modulation studies to operational space weather forecasting, particularly for radiation environment assessment, long-term solar cycle prediction, and validation of numerical transport models.

Future space weather services will increasingly integrate global observational networks, high-performance numerical simulations, and explainable artificial intelligence to provide continuous predictions of the heliospheric radiation environment. Digital twins of the heliosphere, real-time data assimilation from Parker Solar Probe, Solar Orbiter, Aditya-L1, AMS-02, and global neutron monitor networks, together with Physics-Informed Neural Networks, are expected to enable highly accurate forecasts of cosmic ray modulation and radiation hazards. These advances will be critical for safeguarding astronauts, protecting satellites and electronic infrastructure, supporting aviation operations, and enabling sustainable human exploration of the Moon, Mars, and deep space.

Conclusions

Cosmic ray modulation is a complex, multiscale phenomenon arising from the continuous interaction between galactic cosmic rays and the dynamic heliosphere. As energetic charged particles propagate from the local interstellar medium toward the inner Solar System, their intensity, energy spectrum, and spatial distribution are modified by the combined effects of anisotropic diffusion, solar wind convection, gradient and curvature drifts, adiabatic energy changes, and interactions with heliospheric magnetic turbulence. These transport processes are governed by variations in solar activity over both the approximately 11-year sunspot cycle and the 22-year Hale magnetic cycle, making cosmic ray modulation one of the most important manifestations of Sun–heliosphere coupling. The Parker Transport Equation remains the fundamental theoretical framework for describing these processes, while decades of observational and computational research have substantially refined our understanding of the physical mechanisms responsible for long-term and transient modulation.

Long-term observations spanning Solar Cycles 20–25 have demonstrated the close relationship between solar magnetic activity and galactic cosmic ray intensity. Continuous measurements from Voyager, ACE, Ulysses, SOHO, STEREO, PAMELA, AMS-02, Parker Solar Probe, Solar Orbiter, Aditya-L1, and worldwide neutron monitor networks, including the Oulu Neutron Monitor and the Neutron Monitor Database (NMDB), have provided an unprecedented observational foundation for studying heliospheric transport. These datasets have confirmed the inverse correlation between sunspot number and cosmic ray intensity, quantified rigidity-dependent modulation, validated charge-sign effects associated with particle drifts, and characterized the influence of coronal mass ejections, magnetic clouds, and recurrent solar wind structures on transient phenomena such as Forbush decreases. Direct measurements of the Local Interstellar Spectrum by Voyager 1 beyond the heliopause have further established the outer boundary conditions required for modern transport models.

Advances in numerical simulation have transformed the study of cosmic ray modulation from simplified one-dimensional analytical approximations to sophisticated three-dimensional, time-dependent models capable of reproducing realistic heliospheric conditions. Finite-difference methods, stochastic differential equation (SDE) techniques, Monte Carlo simulations, and magnetohydrodynamic (MHD) models have enabled increasingly accurate solutions of the Parker Transport Equation. Dedicated heliospheric frameworks such as HelMod, WSA–ENLIL, and EUHFORIA now incorporate realistic solar wind structures, heliospheric magnetic fields, and transient disturbances to simulate cosmic ray transport across the heliosphere with high fidelity. These numerical approaches, supported by continuous spacecraft observations, have significantly improved the predictive capability of heliophysical models and enhanced our understanding of space weather processes.

Artificial intelligence has emerged as a powerful complement to traditional physics-based modeling. Machine learning techniques—including Long Short-Term Memory (LSTM) networks, Convolutional Neural Networks (CNNs), Transformer architectures, Random Forests, XGBoost, and Physics-Informed Neural Networks (PINNs)—are increasingly being applied to analyze large heliophysical datasets, forecast cosmic ray intensity, estimate transport coefficients, reconstruct missing observations, and predict space weather events. Among these methods, PINNs represent a particularly promising development because they combine observational data with the governing equations of cosmic ray transport, enabling physically consistent predictions while reducing dependence on extensive labeled datasets. The convergence of artificial intelligence with numerical simulation and observational data is rapidly redefining predictive heliophysics and opening new opportunities for real-time forecasting.

The practical significance of cosmic ray modulation extends well beyond fundamental astrophysics. Accurate prediction of the heliospheric radiation environment is essential for protecting astronauts during lunar and interplanetary missions, assessing radiation exposure in commercial aviation, safeguarding satellites and spacecraft electronics, ensuring the reliability of communication and navigation systems, and supporting the design of future planetary exploration missions. Although galactic cosmic rays also contribute to atmospheric ionization and influence several atmospheric processes, current evidence indicates that their role in recent global climate change is secondary compared with dominant anthropogenic forcing. Nevertheless, long-term records of cosmogenic isotopes and neutron monitor observations remain invaluable for reconstructing past solar activity and investigating Sun–Earth interactions over centennial and millennial timescales.

Observational investigations by Agarwal and Mishra have further demonstrated the importance of rigidity dependence, heliospheric magnetic polarity, solar wind velocity, interplanetary magnetic field strength, and recurrent 27-day solar rotational variations in governing long-term galactic cosmic ray modulation. Their analyses, based on multiple solar cycles and neutron monitor observations, complement international spacecraft missions and provide important validation of modern transport theories. Such long-duration observational studies underscore the necessity of combining ground-based measurements with spaceborne observations to achieve a comprehensive understanding of heliospheric modulation.

Looking ahead, the next generation of cosmic ray research will increasingly rely on the integration of heliophysics, high-resolution observations, advanced numerical simulation, and artificial intelligence. Continuous measurements from missions such as Parker Solar Probe, Solar Orbiter, Aditya-L1, AMS-02, and future deep-space observatories, combined with global neutron monitor networks and real-time data assimilation will provide increasingly complete descriptions of the evolving heliosphere. The development of digital twins of the heliosphere, hybrid physics–AI models, explainable machine learning algorithms, and high-performance computational frameworks is expected to substantially improve forecasting of cosmic ray modulation, solar energetic particle events, and radiation hazards. Such integrated approaches will be indispensable for operational space weather services and for the safe expansion of human activities beyond Earth.

In summary, cosmic ray modulation is inherently a multiscale, multidisciplinary process that cannot be fully understood through observations, theory, numerical simulation, or artificial intelligence alone. Continued progress requires the integration of heliophysical observations, plasma physics, magnetohydrodynamic modeling, cosmic ray transport theory, advanced computational methods, and data-driven artificial intelligence within a unified scientific framework. Such an interdisciplinary approach will not only deepen our understanding of energetic particle transport throughout the heliosphere but will also strengthen our ability to predict and mitigate space weather hazards in an era of expanding technological dependence and sustained human exploration of the Solar System.

Acknowledgements

The authors express their sincere gratitude to the scientific community whose pioneering contributions in heliophysics, cosmic ray physics, astrophysics, plasma physics, and space weather research have laid the foundation for the advancements reviewed in this article. The authors acknowledge the invaluable observational data provided by numerous international space missions and ground-based facilities, including Voyager 1 and 2, ACE, Ulysses, SOHO, STEREO, PAMELA, AMS-02, Parker Solar Probe, Solar Orbiter, Aditya-L1, the Oulu Neutron Monitor, and the Neutron Monitor Database (NMDB), which have significantly advanced the understanding of galactic cosmic ray modulation and heliospheric processes. The authors also appreciate the efforts of researchers and organizations that have developed numerical modeling frameworks and openly shared scientific datasets, thereby facilitating progress in heliophysics and space weather research. Finally, the authors thank their respective institutions for providing the academic environment and research support necessary for the preparation of this review article.

Funding

This research received no external funding.

Conflict of Interest

The authors declare no conflict of interest.

References

1. Adriani, O., Barbarino, G. C., Bazilevskaya, G. A., Bellotti, R., Boezio, M., Bogomolov, E. A., ... & PAMELA Collaboration. (2013). Time dependence of the proton flux measured by PAMELA during the 2006–2009 solar minimum. The Astrophysical Journal, 765(2), 91. https://doi.org/10.1088/0004-637X/765/2/91
2. Agarwal, R. and Mishra, R.K., Galactic cosmic ray modulation during last four solar cycles, 31st Int. Cosmic Ray Conf., Lodz (Poland), 2009, SH 3.4, icrc0015.
3. Agarwal R.; Mishra, Rajesh Kumar, Modulation of Cosmic Rays at Different Cut Off Rigidity, 33rd International Cosmic Ray Conference (ICRC 2013), Rio de Janeiro, Brazil 2-9 July 2013, Volume 1 of 5, ISBN: 978-1-5108-1008-2, 3684-3687.
4. Mishra, R.K. and Agarwal, R. (2023), Impact of interplanetary parameters on cosmic ray 27-day variation, International Journal of Scientific Research in Engineering and Management, ISSN: 2852-3939, 2023, 7 (11), 1- 7..
5. Aguilar, M., Ali Cavasonza, L., Alpat, B., Ambrosi, G., Arruda, L., Attig, N., ... AMS Collaboration. (2015). Precision measurement of the proton flux in primary cosmic rays from rigidity 1 GV to 1.8 TV with the Alpha Magnetic Spectrometer on the International Space Station. Physical Review Letters, 114(17), 171103. https://doi.org/10.1103/PhysRevLett.114.171103
6. Anderson, C. D. (1933). The positive electron. Physical Review, 43(6), 491–494.
7. Bell, A. R. (1978). The acceleration of cosmic rays in shock fronts—I. Monthly Notices of the Royal Astronomical Society, 182(2), 147–156.
8. Blandford, R. D., & Eichler, D. (1987). Particle acceleration at astrophysical shocks: A theory of cosmic ray origin. Physics Reports, 154(1), 1–75.
9. Blandford, R. D., & Ostriker, J. P. (1978). Particle acceleration by astrophysical shocks. The Astrophysical Journal Letters, 221, L29–L32.
10. Blasi, P. (2013). The origin of galactic cosmic rays. Astronomy and Astrophysics Review, 21, 70.
11. Bobik, P., Boschini, M. J., Della Torre, S., Gervasi, M., Grandi, D., Jóhannesson, G., ... Potgieter, M. S. (2012). Systematic investigation of solar modulation of galactic protons using the HelMod model. The Astrophysical Journal, 745, 132.
12. Burger, R. A., & Hattingh, M. (1998). Propagation of cosmic rays in the heliosphere. The Astrophysical Journal, 505, 244–252.
13. Burger, R. A., Potgieter, M. S., & Heber, B. (2000). Rigidity dependence of cosmic ray modulation. Journal of Geophysical Research, 105(A12), 27447–27456.
14. Camporeale, E. (2019). The challenge of machine learning in space weather. Space Weather, 17(8), 1166–1207.
15. Cane, H. V. (2000). Coronal mass ejections and Forbush decreases. Space Science Reviews, 93, 55–77.
16. Cucinotta, F. A., Kim, M. H. Y., Chappell, L. J., & Huff, J. L. (2013). How safe is safe enough? Radiation risk for a human mission to Mars. PLoS ONE, 8(10), e74988.
17. Decker, R. B., Krimigis, S. M., Roelof, E. C., Hill, M. E., Armstrong, T. P., Gloeckler, G., ... Stone, E. C. (2008). Mediation of the solar wind termination shock by non-thermal ions. Nature, 454, 67–70.
18. Desai, M., & Giacalone, J. (2016). Large gradual solar energetic particle events. Living Reviews in Solar Physics, 13, 3.
19. Domingo, V., Fleck, B., & Poland, A. I. (1995). The SOHO mission: An overview. Solar Physics, 162, 1–37.
20. Fermi, E. (1949). On the origin of the cosmic radiation. Physical Review, 75(8), 1169–1174.
21. Forbush, S. E. (1954). World-wide cosmic-ray variations, 1937–1952. Journal of Geophysical Research, 59(4), 525–542.
22. Fox, N. J., Velli, M. C., Bale, S. D., Decker, R., Driesman, A., Howard, R. A., ... Szabo, A. (2016). The Parker Solar Probe mission. Space Science Reviews, 204, 7–48.
23. Gaisser, T. K., Engel, R., & Resconi, E. (2016). Cosmic rays and particle physics (2nd ed.). Cambridge University Press.
24. Gleeson, L. J., & Axford, W. I. (1968). Solar modulation of galactic cosmic rays. The Astrophysical Journal, 154, 1011–1026.
25. Greisen, K. (1966). End to the cosmic-ray spectrum? Physical Review Letters, 16(17), 748–750.
26. Heber, B., & Potgieter, M. S. (2006). Cosmic rays at high heliolatitudes. Space Science Reviews, 127, 117–194.
27. Hess, V. F. (1912). Über Beobachtungen der durchdringenden Strahlung bei sieben Freiballonfahrten. Physikalische Zeitschrift, 13, 1084–1091.
28. Jokipii, J. R. (1966). Cosmic-ray propagation. I. Charged particles in a random magnetic field. The Astrophysical Journal, 146, 480–487.
29. Jokipii, J. R. (1971). Propagation of cosmic rays in the solar wind. Reviews of Geophysics, 9(1), 27–87.
30. Jokipii, J. R., Levy, E. H., & Hubbard, W. B. (1977). Effects of particle drift on cosmic-ray transport. The Astrophysical Journal, 213, 861–868.
31. Kirkby, J., et al. (2016). Ion-induced nucleation of pure biogenic particles. Nature, 533, 521–526.
32. Kotera, K., & Olinto, A. V. (2011). The astrophysics of ultrahigh-energy cosmic rays. Annual Review of Astronomy and Astrophysics, 49, 119–153.
33. Lazarian, A., Eyink, G. L., Jafari, A., Kowal, G., Li, H., Xu, S., & Vishniac, E. T. (2020). 3D magnetic reconnection and turbulent particle acceleration. Physics of Plasmas, 27, 012305.
34. Müller, D., St. Cyr, O. C., Zouganelis, I., Gilbert, H., Marsden, R., Nieves-Chinchilla, T., ... Walsh, A. (2020). The Solar Orbiter mission. Astronomy & Astrophysics, 642, A1.
35. Parker, E. N. (1958). Dynamics of the interplanetary gas and magnetic fields. The Astrophysical Journal, 128, 664–676.
36. Parker, E. N. (1965). The passage of energetic charged particles through interplanetary space. Planetary and Space Science, 13(1), 9–49.
37. Pei, C., Bieber, J. W., Burger, R. A., & Clem, J. (2010). A three-dimensional model of galactic cosmic ray modulation. Journal of Geophysical Research: Space Physics, 115, A12107.
38. Pomoell, J., & Poedts, S. (2018). EUHFORIA: European heliospheric forecasting information asset. Journal of Space Weather and Space Climate, 8, A35.
39. Potgieter, M. S. (2013). Solar modulation of cosmic rays. Living Reviews in Solar Physics, 10, 3.
40. Reames, D. V. (1999). Particle acceleration at the Sun and in the heliosphere. Space Science Reviews, 90, 413–491.
41. Reames, D. V. (2013). The two sources of solar energetic particles. Space Science Reviews, 175, 53–92.
42. Schlickeiser, R. (2002). Cosmic ray astrophysics. Springer.
43. Schrijver, C. J., Kauristie, K., Aylward, A. D., Denardini, C. M., Gibson, S. E., Glover, A., ... Vilmer, N. (2015). Understanding space weather to shield society. Advances in Space Research, 55(12), 2745–2807.
44. Simpson, J. A. (1983). Elemental and isotopic composition of the galactic cosmic rays. Annual Review of Nuclear and Particle Science, 33, 323–382.
45. Simpson, J. A. (2000). The cosmic ray nucleonic component: The invention and scientific uses of the neutron monitor. Space Science Reviews, 93, 11–32.
46. Stone, E. C., Cummings, A. C., McDonald, F. B., Heikkila, B., Lal, N., & Webber, W. R. (2013). Voyager 1 observes low-energy galactic cosmic rays beyond the heliopause. Science, 341(6142), 150–153.
47. Strauss, R. D., Potgieter, M. S., Büsching, I., & Kopp, A. (2011). Modeling galactic cosmic ray modulation using stochastic differential equations. The Astrophysical Journal, 735, 83.
48. Strong, A. W., Moskalenko, I. V., & Ptuskin, V. S. (2007). Cosmic-ray propagation and interactions in the Galaxy. Annual Review of Nuclear and Particle Science, 57, 285–327.
49. Usoskin, I. G. (2017). A history of solar activity over millennia. Living Reviews in Solar Physics, 14, 3.
50. Usoskin, I. G., Mursula, K., Kovaltsov, G. A., et al. (2015). The worldwide neutron monitor network: Current status and applications. Journal of Space Weather and Space Climate, 5, A28.
51. Zatsepin, G. T., & Kuzmin, V. A. (1966). Upper limit of the spectrum of cosmic rays. JETP Letters, 4, 78–80.
52. Zhang, M. (1999). A Markov stochastic process theory of cosmic-ray modulation. The Astrophysical Journal, 513(2), 409–420.