Abstract This thesis studies the structure and evolution of the large scale heliospheric magnetic field. The work covers the space age, the period when satellite measurements revolutionized our knowledge about space. Now, this period is known to be the declining phase of the grand modern maximum of solar activity. The thesis addresses how the hemispherical asymmetry of solar activity is seen in the photospheric magnetic field and how it appears in the corona and in the heliosphere until the termination shock. According to geomagnetic and heliospheric observations, the heliospheric current sheet has been southward shifted around the solar minima since 1930s. Using Ulysses probe observations, we derive an accurate estimate of 2° for the southward shift of the heliospheric current sheet during two very different solar minimum in the mid 1990s and 2000s. The overall structure of the heliospheric magnetic field has changed significantly now when the grand modern maximum has come to an end. During the present low solar activity the polar fields are weaker and the heliospheric current sheet covered a wide latitudinal range during the previous minimum. When the heliospheric current sheet is wide the asymmetry is less visible at the Earth’s orbit. We extend our study to the outer heliosphere using measurements made by Voyager and Pioneer probes and show that the hemispherical asymmetry in the coronal hole evolution, and the related southward shift of the heliospheric current sheet, are seen until the termination shock. In order to understand the origin of the hemispherical asymmetry, we complete a multipole analysis of the solar magnetic field since 1976. We find that the minimum time southward shift of the heliospheric current sheet is due to the quadrupole component of the coronal magnetic field. The quadrupole term exists because the generation and transport of the magnetic flux in the Sun tends to proceed differently in the northern and southern hemispheres. During this and the following decade the Sun is most likely going to be less active than it has been since 1920s. Therefore it is probable that the hemispherical asymmetry of the heliospheric magnetic field will be less visible in the ecliptic plane in the near future. Now, when the Sun seems to be at the maximum of cycle 24, we are looking forward to see how the polar fields and the heliospheric magnetic field are formed when approaching the following solar minimum. It is possible that, as the activity rises again after the present and future low cycles, the hemispherical asymmetry will be opposite to that of the 20th century and the minimum time heliospheric current sheet would be northward shifted.

The extreme temperatures and nonthermal nature of the solar corona and solar wind arise from an unidentified physical mechanism that preferentially heats certain ion species relative to others. Spectroscopic indicators of unequal temperatures commence within a fraction of a solar radius above the surface of the Sun, but the outer reach of this mechanism has yet to be determined. Here we present an empirical procedure for combining interplanetary solar wind measurements and a modeled energy equation including Coulomb relaxation to solve for the typical outer boundary of this zone of preferential heating. Applied to two decades of observations by the Wind spacecraft, our results are consistent with preferential heating being active in a zone extending from the transition region in the lower corona to an outer boundary 20-40 solar radii from the Sun, producing a steady-state super-massproportional a-to-proton temperature ratio of 5.2-5.3. Preferential ion heating continues far beyond the transition region and is important for the evolution of both the outer corona and the solar wind. The outer boundary of this zone is well below the orbits of spacecraft at 1 au and even closer missions such as Helios and MESSENGER, meaning it is likely that no existing mission has directly observed intense preferential heating, just residual signatures. We predict that the Parker Solar Probe will be the first spacecraft with a perihelion sufficiently close to the Sun to pass through the outer boundary, enter the zone of preferential heating, and directly observe the physical mechanism in action.

Le but de cette thèse est l'étude du développement de la turbulence dans le vent solaire entre 0.2 et 1 unité astronomique (UA) du soleil (i.e. l'orbite terrestre). L'étude est faite en résolvant numériquement les équations de la MHD après soustraction de l'écoulement moyen radial. Les deux aspects de la turbulence qui nous intéressent sont la structure 3D des spectres d’énergie et le chauffage du plasma qui résulte de la dissipation turbulente des tourbillons et couches de courant emportés par le vent. On cherche à déterminer quelles sont les conditions du plasma près du soleil qui permettent d’aboutir à ce qu'on observe à 1 UA. Un but important de mon travail est aussi de déterminer si la physique qui est présente dans les équations que j'intègre (la MHD) suffit pour arriver à reproduire ce qu'on a déjà observé dans cet intervalle de distance. Nous introduisons le contexte de notre travail dans la première partie. On y trouve les équations de base, une introduction à la turbulence, un résumé sur la physique du vent solaire et de la couronne solaire. La partie 2 sera consacrée à l'étude de l'anisotropie de la cascade turbulente, et plus précisément du spectre 3D. Dans la zone inertielle, les mesures in-situ vers 1 UA montrent des figures complexes pour ces spectres qu'on peut interpréter de plusieurs façons : nos simulations numériques permettent de lever toute ambiguïté. Plus précisément, la question est de savoir quand intervient l'axe soleil-terre, et quand intervient l'axe du champ magnétique moyen. La partie trois est centrée sur le chauffage turbulent dans les vents rapides et lents. Entre 0.3 et 1 UA, la température des protons diminue anormalement lentement, ce qui indique une source de chauffage, qu'on suppose ici être la dissipation des tourbillons et couches de courant emportés par le vent. Pour démontrer que cette hypothèse est raisonnable, nous considérons d’abord le modèle de Burgers qui est un modèle pour l'évolution d’ondes sonores. Ensuite, nous passons à l'étude du cas plus complexe d'un volume de plasma 3D. Nous examinerons les conditions initiales correspondant aux vents lents et rapides. Dans les deux cas, on adoptera des anisotropies spectrales différentes. Dans la dernière partie, nous exposerons les conclusions de notre travail et proposerons d'introduire l'anisotropie de la température dans un travail futur. / The aim of this thesis is the study of the development of turbulence in the solar wind between 0.2 and 1 astronomical unit (AU) from the Sun (i.e. Earth’s orbit). The study is done by solving the magnetohydrodynamics equations (MHD) after subtracting the mean radial flow. The two aspects of turbulence that interest us are the 3D structure of the energy spectra and the heating of plasma that results from the turbulent dissipation of eddies and current layers transported by the wind. We want to determine which conditions of the plasma close to the Sun can result into what we observe at 1 AU. We have relatively detailed measurements of what happens between 0.3 and 1 AU. One important goal of this work is to determine if the physics present in the equations that are integrated (MHD) is sufficient to reproduce what is observed in this interval of distances. We introduce the context of our work in the first part. We give a summary of the physics concerning the solar wind and the solar corona, and the basic equations used to describe the solar wind plasma and an introduction to turbulence. Part 2 is dedicated to the study of anisotropy in the turbulent cascade, which characterizes 3D spectra. In the inertial range, in-situ measurements at 1 AU show complex figures for these spectra that we can interpret in several ways : numerical simulations allow to clear ambiguities. An important question is to know whether the Earth-Sun symmetry axis or the mean magnetic field axis is dominant.The third part focuses on turbulent heating in fast and slow winds. Between 0.3 and 1 AU, proton temperature decreases more slowly than expected, which requires a heating source. This source is supposed to be the continuous dissipation of eddies and current layers transported by the wind. To start with, we consider the simple case of Burgers equation, which is a one-dimensional model for shock formation. Thereupon, we switch to the 3-dimensional case, where we consider initial conditions appropriate for slow and fast winds. In the last part we expose our conclusions and propose the implementation of temperature anisotropy as future work.

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.

A numerical scheme has been developed to solve the quasilinear form of the transonic stream function equation. The method is applied to compute steady two-dimensional axisymmetric solar wind-type problems. A single, perfect, non-dissipative, homentropic and polytropic gas-dynamics is assumed. The four equations governing mass and momentum conservation are reduced to a single nonlinear second order partial differential equation for the stream function. Bernoulli's equation is used to obtain a nonlinear algebraic relation for the density in terms of stream function derivatives. The vorticity includes the effects of azimuthal rotation and Bernoulli's function and is determined from quantities specified on boundaries. The approach is efficient. The number of equations and independent variables has been reduced and a rapid relaxation technique developed for the transonic full potential equation is used. Second order accurate central differences are used in elliptic regions. In hyperbolic regions a dissipation term motivated by the rotated differencing scheme of Jameson is added for stability. A successive-line-overrelaxation technique also introduced by Jameson is used to solve the equations. The nonlinear equationfor the density is a double valued function of the stream function derivatives. The velocities are extrapolated from upwind points to determine the proper branch and Newton's method is used to iteratively compute the density. This allows accurate solutions with few grid points. The applications first illustrate solutins to solar wind models. The equations predict that the effects of vorticity must be confined near the surface and far away the streamlines must resemble the spherically symmetric solution. Irrotational and rotational flows show this behavior. The streamlines bend toward the rotation axis for rapidly rotating models because the coriolis force is much larger than the centrifugal force. Models of galactic winds are computed by considering the flow exterior to a surface which surrounds a uniform density oblate spheroid. Irrotational results with uniform outward mass flux show streamlines bent toward the equator and nearly spherical sonic surfaces. Rotating models for which Bernoulli's function is not constant show the sonic surface is deformed consistent with the one-dimensional theory.

Rapid changes in relative helium abundance in solar wind are often attributed to the flux tube crossings - thus between solar wind streams coming from different parts of the solar corona. Recent studies however show, that the changes in relative helium abundance are not always correlated with changes in other parameters and therefore the relative helium abundance could be changed even during the propagation of solar wind inside a flux tube. The topic of this thesis is an analysis of rapid changes of solar wind parameters and their mutual connection. Using a multipoint (we use data from the SPEKTR-R, WIND and the THEMIS spacecraft) case study of interplanetary shocks, we show that a possible explanation of those changes could be a magnetosonic instability excited by the differential flow between helium and proton populations.

The solar wind interaction with an icy comet is studied through a model problem. A hybrid simulation is done of a box with evenly distributed water ions and protons, where initially the water ions are stationary, and protons move with the speed of the solar wind. The purpose of the thesis is to investigate the interaction between the two species through the convective electric field, and focus is on early acceleration of pick-up ions, and deflection of the solar wind. It is relevant to the cometary case, because it enables study of the physics of this interaction, without involving other mechanisms, such as bow shock, magnetic field pile-up and draping. The species are found to exchange kinetic energy similar to a damped oscillator, where the dampening is caused by kinetic energy being transferred to the magnetic field. At early times, i.e. times smaller than the gyration time for the water ions, the solar wind does not lose much speed when it is deflected. For comparable number densities, the solar wind can be deflected more than 90° at early times, and loses more speed, and water ions are picked up faster. The total kinetic energy of the system decreases when energy builds up in the magnetic field. The nature of the energy exchange is strongly dependent on the number density ratio between water ions and protons. A density instability with behaviour similar to a plasma beam instability forms as energy in the magnetic field increases, and limits the amount of time the simulation preserves total energy, for the particular hybrid solver used. There is a discussion on the structure of the density instability, and it is compared to cometary simulations.

La thèse porte sur le rôle de la rotation planétaire dans la structure globale de l'interaction vent solaire/magnétosphère à partir de simulations magnétohydrodynamiques (MHD). Les magnétosphères planétaires du système solaire présentent une incroyable diversité, et notamment dans leurs configurations respectives de l'inclinaison de leur axe magnétique par rapport à leur axe de rotation. La durée des périodes de rotation par rapport au temps de relaxation de chaque magnétosphère diffère aussi d'une planète à l'autre. On distingue ainsi les rotateurs lents (Mercure et la Terre), pour lesquels le temps de relaxation est plus court que la période de rotation, des rotateurs rapides (Jupiter, Saturne, Uranus et Neptune). Dans le cas du rotateur lent Mercure, on s'intéresse à l'influence des paramètres du vent solaire sur la structure globale du champ magnétique et de l'écoulement. En appui à la mission spatiale BepiColombo, nous présentons des simulations effectuées pour deux modèles différents de champ magnétique herméen. Nous détaillons le rôle des fronts d'onde MHD stationnaires, en particulier les fronts stationnaires de mode lent dans la magnétogaine. Saturne présente la particularité d'avoir un axe magnétique parfaitement aligné avec son axe de rotation. C'est donc un cas de rotateur rapide stationnaire, qui nous permet d'étudier la structure globale du champ magnétique et de l'écoulement pour différentes orientations de l'IMF, mais aussi pour différentes vitesses de rotation de la planète. Enfin, le cas d'une configuration quelconque, avec un grand angle entre l'axe magnétique et l'axe de rotation planétaire, est étudié en présence d'un vent solaire magnétisé en s'inspirant de la configuration d'Uranus au solstice et à l'équinoxe. Dans la configuration « solstice », c'est à dire lorsque l'axe de rotation pointe vers le Soleil, on montre qu'une structure de nature alfvénique se développe en hélice dans la queue de la magnétosphère, et que les zones de reconnexion entre le champ magnétique planétaire et l'IMF, qui forment aussi une double hélice, ralentissent la progression de la structure alfvénique. A l'équinoxe, lorsque l'axe de rotation est toujours dans le plan de l’écliptique mais perpendiculaire à la direction Soleil-Uranus, la structure en hélice disparaît. / The topic of the thesis is the part of planetary rotation in the global structure of the solar wind interaction with planetary magnetospheres using MHD simulations. We discuss the distinction between slow and fast rotators from a MHD point of view. In the case of a non-rotating magnetosphere (as is the one of Mercury), the part of standing MHD modes is studied, along with a method to identify them in simulations. A fast-rotating but stationary magnetosphere (inspired by the case of Saturn) is presented in details and provides a good test to validate the new version of the AMRVAC code allowing for any configuration regarding the respective directions of the planetary spin axis, planetary magnetic axis, solar wind inflow direction, and IMF orientation. Finally, a random configuration, with a large angle between the planetary spin and magnetic axis, is analyzed for the first time in presence of a magnetized solar wind, using configurations inspired from the planet Uranus at solstice and equinox.

In space and astrophysical plasmas, turbulence is responsible for transferring energy from large scales driven by violent events or instabilities, to smaller scales where turbulent energy is ultimately converted into plasma heat by dissipative mechanisms. In the inertial range, the self-similar turbulent energy cascade to smaller spatial scales is driven by the nonlinear interaction between counterpropagating Alfvén waves, denoted Alfvén wave collisions. For the more realistic case of the collision between two initially separated Alfvén wavepackets (rather than previous idealized, periodic cases), we use a nonlinear gyrokinetic simulation code, AstroGK, to demonstrate three key properties of strong Alfvén wave collisions: they (i) facilitate the perpendicular cascade of energy and (ii) generate current sheets self-consistently, and (iii) the modes mediating the nonlinear interaction are simply Alfvén waves. Once the turbulent cascade reaches the ion gyroradius scale, the Alfvén waves become dispersive and the turbulent energy starts to dissipate, energizing the particles via wave-particle interactions with eventual dissipation into plasma heat. The novel Field-Particle Correlation technique determines how turbulent energy dissipates into plasma heat by identifying which particles in velocity-space experience a net gain of energy. By utilizing knowledge of discrete particle arrival times, we devise a new algorithm called PATCH (Particle Arrival Time Correlation for Heliophysics) for implementing a field-particle correlator onboard spacecraft. Using AstroGK, we create synthetic spacecraft data mapped to realistic phase-space resolutions of modern spacecraft instruments. We then utilize Poisson statistics to determine the threshold number of particle counts needed to resolve the velocity-space signature of ion Landau damping using the PATCH algorithm.

<p>The Neutral Particle Detector (NPD) is a new type of instrumentation for energetic neutral atom (ENA) diagnostics. This thesis deals with development of the NPD sensor designed as a part of the plasma and neutral particle packages ASPERA-3 and ASPERA-4 on board Mars Express and Venus Express, the European Space Agency (ESA) satellites to Mars and Venus, respectively. It describes how the NPD sensors were designed, developed, tested and calibrated. </p><p>It also presents the first scientific results obtained with NPD during its operation at Mars. </p><p>The NPD package consists of two identical detectors, NPD1 and NPD2. Each detector has a 9^o x 90^o intrinsic field-of-view divided into three sectors. The ENA detection principle is based on the surface interaction technique. NPD detects ENA differential fluxes within the energy range of 100 eV to 10 keV and is capable of resolving hydrogen and oxygen atoms by time-of-flight (TOF) measurements or pulse height analysis.</p><p>During the calibration process the detailed response of the sensor was defined, including properties such as an angular response function and energy dependent efficiency of each of the sensor sectors for different ENA species. </p><p>Based on the NPD measurements at Mars the main scientific results reported so far are:</p><p>- observation of the Martian H-ENA jet / cone and its dynamics, </p><p>- observations of ENA emissions from the Martian upper atmosphere, </p><p>- measurements of the hydrogen exosphere density profile at Mars, </p><p>- observations of the response of the Martian plasma environment to an interplanetary shock, </p><p>- observations of the H-ENA fluxes in the interplanetary medium.</p>