The interaction between the solar wind and the partially ionized gas of the local interstellar medium (ISM) creates a bubble known as the heliosphere. Classically, the shape of the heliosphere has been regarded as comet-like, with a long tail pointed in the direction opposite the Sun’s motion through the ISM. In this view, the solar magnetic field was assumed to have a negligible effect on the global structure of the heliosphere. Recent advances in numerical modeling have revealed the importance of the solar magnetic field in its ability to confine and collimate the solar wind plasma, and the shape of the heliosphere has been called into question. Energetic neutral atoms (ENAs) are created throughout the heliosphere via charge exchange. The separate contributions of the solar magnetic field topology and the solar wind structure to ENA observations is largely unexplored. The Interstellar Boundary Explorer (IBEX) has been providing a global perspective of the heliosphere through ENA maps with energies ranging from 0.2 to 6 keV. In this dissertation, three-dimensional magnetohydrodynamic simulations of the heliosphere are used as input to an ENA model designed to produce synthetic ENA maps. I compare modeled ENA maps with IBEX observations to investigate how different heliospheric conditions and properties affect ENAs created in the heliosphere, and therefore how ENA observations can be used to understand the heliosphere. First, I investigate the effect of the solar wind collimation by the solar magnetic field on ENA maps in the case of a solar wind without latitudinal variation. I find that even in the absence of variations of the solar wind, two lobes of strong ENA flux form at high latitudes, similar to what is observed by IBEX at high energies. Second, I test the effect of a latitudinally-varying solar wind on ENAs both with and without the inclusion of the solar magnetic field. I show that the latitudinal variations of the solar wind during solar minimum creates a structured ENA profile with latitude, corresponding to the profile observed at 1 AU, but that the solar magnetic field significantly enhances ENA flux in the region where the solar wind is confined. Lastly, I investigate the effect of the solar cycle on ENAs and how changing solar wind conditions (e.g. density, temperature, velocity) affect the heliosphere over time. I demonstrate that, given changes in the solar cycle, there is a significant evolution in the modeled ENA flux due to the changes in the solar wind profile and the solar magnetic field, which is also seen by ENA observations.

The modulation of galactic cosmic ray (GCR) Carbon in a north-south asymmetrical heliosphere is studied, using a two-dimensional numerical model that contains a solar wind termination shock (TS), a heliosheath, as well as particle drifts and diffusive shock re-acceleration of GCRs. The asymmetry in the geometry of the heliosphere is incorporated in the model by assuming a significant dependence on heliolatitude of the thickness of the heliosheath. As a result, the model allows comparisons of modulation in the north and south hemispheres during both magnetic polarity cycles of the Sun, and from solar minimum to moderate maximum conditions. When comparing the computed spectra between polar angles of 55o (approximating the Voyager 1 direction) and 125o (approximating the Voyager 2 direction), it is found that at kinetic energies E < 1:0 GeV/nuc the effects of the assumed asymmetry in the geometry of the heliosphere on the modulated spectra are insignificant up to 60 AU from the Sun, but become increasingly more significant with larger radial distances to reach a maximum inside the heliosheath. In contrast, with E > 1:0 GeV/nuc, these effects remained insignificant throughout the heliosphere even very close to the heliopause (HP). However, when the enhancement of both polar and radial perpendicular diffusion coefficients off the equatorial plane is assumed to differ from heliographic pole to pole, reflecting different modulation conditions between the two hemispheres, major differences in the computed intensities between the two Voyager directions are obtained throughout the heliosphere. The model is further improved by incorporating new information about the HP location and the relevant heliopause spectrum for GCR Carbon at E < 200 MeV/nuc based on the recent Voyager 1 observations. When comparing the computed solutions at the Earth with ACE observations taken during different solar modulation conditions, it is found that it is possible for the level of modulation at the Earth, when solar activity changes from moderate maximum conditions to solar minimum conditions, to exceed the total modulation between the HP and the Earth during solar minimum periods. In the outer heliosphere, reasonable compatibility with the corresponding Voyager observations is established when drifts are scaled down to zero in the heliosheath in both polarity cycles. The effects of neglecting drifts in the heliosheath are found to be more significant than neglecting the enhancement of polar perpendicular diffusion. Theoretical expressions for the scattering function required for the reduction of the drift coefficient in modulation studies are illustrated and implemented in the numerical model. It is found that when this scattering function decreases rapidly over the poles, the computed A < 0 spectra are higher than the A > 0 spectra at all energies at Earth primarily because of drifts, which is unexpected from a classical drift modeling point of view. Scenarios of this function with strong decreases over the polar regions seem realistic at and beyond the TS, where the solar wind must have a larger latitudinal dependence. / PhD (Space Physics), North-West University, Potchefstroom Campus, 2015

The modulation of galactic cosmic ray (GCR) Carbon in a north-south asymmetrical heliosphere is studied, using a two-dimensional numerical model that contains a solar wind termination shock (TS), a heliosheath, as well as particle drifts and diffusive shock re-acceleration of GCRs. The asymmetry in the geometry of the heliosphere is incorporated in the model by assuming a significant dependence on heliolatitude of the thickness of the heliosheath. As a result, the model allows comparisons of modulation in the north and south hemispheres during both magnetic polarity cycles of the Sun, and from solar minimum to moderate maximum conditions. When comparing the computed spectra between polar angles of 55o (approximating the Voyager 1 direction) and 125o (approximating the Voyager 2 direction), it is found that at kinetic energies E < 1:0 GeV/nuc the effects of the assumed asymmetry in the geometry of the heliosphere on the modulated spectra are insignificant up to 60 AU from the Sun, but become increasingly more significant with larger radial distances to reach a maximum inside the heliosheath. In contrast, with E > 1:0 GeV/nuc, these effects remained insignificant throughout the heliosphere even very close to the heliopause (HP). However, when the enhancement of both polar and radial perpendicular diffusion coefficients off the equatorial plane is assumed to differ from heliographic pole to pole, reflecting different modulation conditions between the two hemispheres, major differences in the computed intensities between the two Voyager directions are obtained throughout the heliosphere. The model is further improved by incorporating new information about the HP location and the relevant heliopause spectrum for GCR Carbon at E < 200 MeV/nuc based on the recent Voyager 1 observations. When comparing the computed solutions at the Earth with ACE observations taken during different solar modulation conditions, it is found that it is possible for the level of modulation at the Earth, when solar activity changes from moderate maximum conditions to solar minimum conditions, to exceed the total modulation between the HP and the Earth during solar minimum periods. In the outer heliosphere, reasonable compatibility with the corresponding Voyager observations is established when drifts are scaled down to zero in the heliosheath in both polarity cycles. The effects of neglecting drifts in the heliosheath are found to be more significant than neglecting the enhancement of polar perpendicular diffusion. Theoretical expressions for the scattering function required for the reduction of the drift coefficient in modulation studies are illustrated and implemented in the numerical model. It is found that when this scattering function decreases rapidly over the poles, the computed A < 0 spectra are higher than the A > 0 spectra at all energies at Earth primarily because of drifts, which is unexpected from a classical drift modeling point of view. Scenarios of this function with strong decreases over the polar regions seem realistic at and beyond the TS, where the solar wind must have a larger latitudinal dependence. / PhD (Space Physics), North-West University, Potchefstroom Campus, 2015

The data harvest from the Voyagers' (V1 and V2) Ultraviolet Spectrometers (UVS) covers encounters with the outer planets, measurements of the heliosphere sky-background, and stellar spectrophotometry. Because their period of operation overlaps with many ultraviolet missions, the calibration of V1. and V2 UVS with other spectrometers is invaluable. Here we revisit the UVS calibration to assess the intriguing sensitivity enhancements of 243% (V1) and 156% (V2) proposed recently. Using the Ly alpha airglow from Saturn, observed in situ by both Voyagers, and remotely by International Ultraviolet Explorer (IUE), we match the Voyager values to IUE, taking into account the shape of the Saturn Ly alpha line observed with the Goddard High Resolution Spectrograph on board the Hubble Space Telescope. For all known ranges of the interplanetary hydrogen density, we show that the V1 and V2 UVS sensitivities cannot be enhanced by the amounts thus far proposed. The same diagnostic holds for distinct channels covering the diffuse He I 58.4 nm emission. Our prescription is to keep the original calibration of the Voyager UVS with a maximum uncertainty of 30%, making both instruments some of the most stable EUV/FUV spectrographs in the history of space exploration. In that frame, we reassess the excess Ly alpha emission detected by Voyager UVS deep in the heliosphere, to show its consistency with a heliospheric but not galactic origin. Our finding confirms results obtained nearly two decades ago-namely, the UVS discovery of the distortion of the heliosphere and the corresponding obliquity of the local interstellar magnetic field (similar to 40 degrees from upwind) in the solar system neighborhood-without requiring any revision of the Voyager UVS calibration.

Thesis (M.Sc. (Physics))--North-West University, Potchefstroom Campus, 2008.

A three-dimensional (3D) steady-state electron modulation model based on Parker (1965) transport equation is applied to study the modelling of – 7 MeV galactic and Jovian electrons in the inner heliosphere. The latter is produced within Jupiter's magnetosphere which is situated at - 5 AU in the ecliptic plane. The heliospheric propagation of these particles is mainly described by the heliospheric diffusion tensor. Some elements of the tensor, such as the diffusion coefficient in the azimuthal direction, which were neglected in the previous two-dimensional modulation studies are investigated to account for the three-dimensional transport of Jovian electrons. Different anisotropic solar wind speed profiles that could represent solar minimum conditions were modelled and their effects were illustrated by computing the distribution of 7 MeV Jovian electrons in the equatorial regions. In particular, the electron intensity time-profile along the Ulysses spacecraft trajectory was calculated for these speed profiles and compared to the 3-10 MeV electron flux observed by the Kiel Electron Telescope (KET) on board the Ulysses spacecraft from launch (1990) up to end of its first out-of-ecliptic orbit (2000). It was found that the model solution computed with the solar wind profile previously assumed for typical solar minimum conditions produced good compatibility with observations up to 1998. After 1998 all model solutions deviated completely from the observations. In this study, as a further attempt to model KET observations more realistically, a new relation is established between the latitudinal dependence of the solar wind speed and the perpendicular polar diffusion. Based on this relation, a transition of an average solar wind speed from solar minimum conditions to intermediate solar activity and to solar maximum conditions was modelled based on the assumption of the time-evolution of large polar coronal holes and were correlated to different scenarios of the enhancement of perpendicular polar diffusion. Effects of these scenarios were illustrated, as a series of steady-state solutions, on the computed 7 MeV Jovian and galactic electrons in comparison with the 3-10 MeV electron observed by the KET instrument from the period 1998 up to the end of 2003. Subsequent effects of these scenarios were also shown on electron modulation in general. It was found that this approach improved modelling of the post-1998 discrepancy between the model and KET observations but it also suggested the need for a time-dependent 3D electron modulation model to describe modulation during moderate to extreme solar maximum conditions. / Thesis (M.Sc.)--North-West University, Potchefstroom Campus, 2004.

The interest in the role of the solar wind termination shock (TS) and heliosheath in cosmic ray (CR) modulation studies has increased sigm6cantly as the Voyager 1 and 2 spacecraft approach the estimated position of the TS. For this work the modulation of galactic CR protons, anti-protons, electrons with a Jovian source, positrons, Helium, and anomalous protons and Helium, and the consequent charge-sign dependence, are studied with an improved and extended two-dimensional numerical CR modulation model including a TS with diffusive shock acceleration, a heliosheath and drifts. The modulation is computed using improved local interstellar spectra (LIS) for almost all the species of interest to this study and new fundamentally derived diffusion coefficients, applicable to a number of CR species during both magnetic polarity cycles of the Sun. The model also allows comparisons of modulation with and without a TS and between solar minimum and moderate maximum conditions. The modulation of protons and Helium with their respective anomalous components are also studied to establish the consequent charge-sign dependence at low energies and the influence on the computed p/p, e-/p, and e-/He. The level of modulation in the simulated heliosheath, and the importance of this modulation 'barrier' and the TS for the different species are illustrated. From the computations it is possible to estimate the ratio of modulation occurring in the heliosheath to the total modulation between the heliopause and Earth for the mentioned species. It has been found that the modulation in the heliosheath depends on the particle species, is strongly dependent on the energy of the CRs, on the polarity cycle and is enhanced by the inclusion of the TS. The computed modulation for the considered species is surprisingly different and the heliosheath is important for CR modulation, although 'barrier' modulation is more prominent for protons, anti-protons and Helium, while the heliosheath cannot really be considered a modulation 'barrier' for electrons and positrons above energies of ~150 MeV. The effects of the TS on modulation are more pronounced for polarity cycles when particles are drifting primarily in the equatorial regions of the heliosphere along the heliospheric current sheet to the Sun, e.g. the A < 0 polarity cycle for protons, positrons, and Helium, and the A > 0 polarity cycle for electrons and anti-protons. This study also shows that the proton and Helium LIS may not be known at energies <~ 200 MeV until a spacecraft actually approaches the heliopause because of the strong modulation that occurs in the heliosheath, the effect of the TS, and the presence of anomalous protons and Helium. For anti-protons, in contrast, these effects are less pronounced. For positrons, with a completely different shape LIS, the modulated spectra have very mild energy dependencies <~ 300 MeV, even at Earth, in contrast to the other species. These characteristic spectral features may be helpful to distinguish between electron and positron spectra when they are measured near and at Earth. These simulations can be of use for future missions to the outer heliosphere and beyond. / Thesis (Ph.D. (Physics))--North-West University, Potchefstroom Campus, 2004.

A two dimensional hydrodynamic model describing the solar wind interaction with the local interstellar medium, which surrounds the solar system, is used to study the heliosphere both as a steady-state- and dynamic structure. The finite volume method used to solve the associated system of hydrodynamic equations numerically is discussed in detail. Subsequently the steady state heliosphere is studied for both the case where the solar wind and the interstellar medium are assumed to consist of protons only, as well as the case where the neutral hydrogen population in the interstellar medium is taken into account. It is shown that the heliosphere forms as three waves, propagating away from the initial point of contact between the solar wind and interstellar matter, become stationary. Two of these waves become stationary at sonic points, forming the termination shock and bow shock respectively. The third wave becomes stationary as a contact discontinuity, called the heliopause. It is shown that the position and geometry of the termination shock, heliopause and bow shock as well as the plasma flow characteristics of the heliosphere largely depend on the dynamic pressure of either the solar wind or interstellar matter. The heliosphere is modelled as a dynamic structure, including both the effects of the solar cycle and short term variations in the solar wind observed by a range of spacecraft over the past ~ 30 years. The dynamic model allows the calculation of an accurate record of the heliosphere state over the past ~ 30 years. This record is used to predict the time at which the Voyager 2 spacecraft will cross the termination shock. Voyager 1 observations of 10 MeV cosmic ray electrons are then used in conjunction with a cosmic ray modulation model to constrain the record of the heliosphere further. It is shown that the dynamic hydrodynamic model describes the heliosphere accurately within a margin of error of ±0.7 years and ±3 AU. The model predicts that Voyager 2 crossed the termination shock in 2007, corresponding to preliminary results from observations indicating that the crossing occurred in August 2007. / Thesis (M.Sc. (Physics))--North-West University, Potchefstroom Campus, 2008.

The time-dependent modulation of galactic cosmic rays in the heliosphere is studied by computing intensities using a two-dimensional, time-dependent modulation model. The compound approach of Ferreira and Potgieter (2004), which describes changes in the cosmic ray transport coefficients over a solar cycle, is improved by introducing recent theoretical advances in the model. Computed intensities are compared with Voyager 1 and 2, IMP 8 and Ulysses proton observations in search of compatibility. It is shown that this approach gives realistic cosmic ray proton intensities on a global scale at Earth and along both Voyager spacecraft trajectories. The results show that cosmic ray modulation, in particular during the present polarity cycle, is not just determined by changes in the drift coefficient but is also dependent on changes in the diffusion coefficients. Furthermore, a comparison of computations to observations along the Voyager 1 and Voyager 2 trajectories illustrates that the heliosphere is asymmetrical. Assuming the latter, E > 70 MeV and 133-242 MeV cosmic ray proton intensities along Voyager 1 and 2 trajectories are predicted from 2012 onwards. It is shown that the computed intensities along Voyager 1 can increase with an almost constant rate since the spacecraft is close to the heliopause. However, the model shows that Voyager 2 is still under the influence of temporal solar activity changes because of the relatively large distance to the heliopause when compared to Voyager 1. Along the Voyager 2 trajectory the intensities should remain generally constant for the next few years and then should start to steadily increase. It is also found that without knowing the exact location of heliopause and transport parameters one cannot conclude anything about local interstellar spectra. The effect of a dynamic inner heliosheath width on cosmic ray modulation is also studied by implementing a time-dependent termination shock position in the model. This does not lead to improved compatibility with spacecraft observations so that a time-dependent termination shock along with a time-dependent heliopause position is required. The variation of the heliopause position over a solar cycle is found to be smaller compared to that of the termination shock. The model predicts the heliopause and termination shock positions along Voyager 1 in 2012 at 119 AU and 88 AU respectively and along Voyager 2 at 100 AU and 84 AU respectively. / Thesis (PhD (Space Physics))--North-West University, Potchefstroom Campus, 2013

The time-dependent modulation of galactic cosmic rays in the heliosphere is studied by computing intensities using a two-dimensional, time-dependent modulation model. The compound approach of Ferreira and Potgieter (2004), which describes changes in the cosmic ray transport coefficients over a solar cycle, is improved by introducing recent theoretical advances in the model. Computed intensities are compared with Voyager 1 and 2, IMP 8 and Ulysses proton observations in search of compatibility. It is shown that this approach gives realistic cosmic ray proton intensities on a global scale at Earth and along both Voyager spacecraft trajectories. The results show that cosmic ray modulation, in particular during the present polarity cycle, is not just determined by changes in the drift coefficient but is also dependent on changes in the diffusion coefficients. Furthermore, a comparison of computations to observations along the Voyager 1 and Voyager 2 trajectories illustrates that the heliosphere is asymmetrical. Assuming the latter, E > 70 MeV and 133-242 MeV cosmic ray proton intensities along Voyager 1 and 2 trajectories are predicted from 2012 onwards. It is shown that the computed intensities along Voyager 1 can increase with an almost constant rate since the spacecraft is close to the heliopause. However, the model shows that Voyager 2 is still under the influence of temporal solar activity changes because of the relatively large distance to the heliopause when compared to Voyager 1. Along the Voyager 2 trajectory the intensities should remain generally constant for the next few years and then should start to steadily increase. It is also found that without knowing the exact location of heliopause and transport parameters one cannot conclude anything about local interstellar spectra. The effect of a dynamic inner heliosheath width on cosmic ray modulation is also studied by implementing a time-dependent termination shock position in the model. This does not lead to improved compatibility with spacecraft observations so that a time-dependent termination shock along with a time-dependent heliopause position is required. The variation of the heliopause position over a solar cycle is found to be smaller compared to that of the termination shock. The model predicts the heliopause and termination shock positions along Voyager 1 in 2012 at 119 AU and 88 AU respectively and along Voyager 2 at 100 AU and 84 AU respectively. / Thesis (PhD (Space Physics))--North-West University, Potchefstroom Campus, 2013